Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userImplicit linearization method for non-standard two-phase flow in porous media
120
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this paper, we consider a non-local (in time) two-phase flow model. The non-locality is introduced through the wettability alteration induced dynamic capillary pressure function. We present a monotone fixed-point iterative linearization scheme for the resulting non-standard model. The scheme treats the dynamic capillary pressure functions semi-implicitly and introduces an $L$-scheme type cite{List2016, Radu2015} stabilization term in the pressure as well as the transport equations. We prove the convergence of the proposed scheme theoretically under physically acceptable assumptions and verify the theoretical analysis with numerical simulations. The scheme is implemented and tested for a variety of reservoir heterogeneity in addition to the dynamic change of the capillary pressure function. The proposed scheme satisfies the predefined stopping criterion within a few numbers of iterations. We also compared the performance of the proposed scheme against the iterative IMplicit Pressure Explicit Saturation scheme
rate research
Read More
A high-performance gas kinetic solver using multi-level parallelization is developed to enable pore-scale simulations of rarefied flows in porous media. The Boltzmann model equation is solved by the discrete velocity method with an iterative scheme. The multi-level MPI/OpenMP parallelization is implemented with the aim to efficiently utilise the computational resources to allow direct simulation of rarefied gas flows in porous media based on digital rock images for the first time. The multi-level parallel approach is analyzed in details confirming its better performance than the commonly-used MPI processing alone for an iterative scheme. With high communication efficiency and appropriate load balancing among CPU processes, parallel efficiency of 94% is achieved for 1536 cores in the 2D simulations, and 81% for 12288 cores in the 3D simulations. While decomposition in the spatial space does not affect the simulation results, one additional benefit of this approach is that the number of subdomains can be kept minimal to avoid deterioration of the convergence rate of the iteration process. This multi-level parallel approach can be readily extended to solve other Boltzmann model equations.
This paper addresses how two time integration schemes, the Heuns scheme for explicit time integration and the second-order Crank-Nicolson scheme for implicit time integration, can be coupled spatially. This coupling is the prerequisite to perform a coupled Large Eddy Simulation / Reynolds Averaged Navier-Stokes computation in an industrial context, using the implicit time procedure for the boundary layer (RANS) and the explicit time integration procedure in the LES region. The coupling procedure is designed in order to switch from explicit to implicit time integrations as fast as possible, while maintaining stability. After introducing the different schemes, the paper presents the initial coupling procedure adapted from a published reference and shows that it can amplify some numerical waves. An alternative procedure, studied in a coupled time/space framework, is shown to be stable and with spectral properties in agreement with the requirements of industrial applications. The coupling technique is validated with standard test cases, ranging from one-dimensional to three-dimensional flows.
Deformable elastic bodies in viscous and viscoelastic media constitute a large portion of synthetic and biological complex fluids. We present a parallelized 3D-simulation methodology which fully resolves the momentum balance in the solid and fluid domains. An immersed boundary algorithm is exploited known as the immersed finite element method (IFEM) which accurately determines the internal forces in the solid domain. The scheme utilized has the advantages of requiring no costly re-meshing, handling finite Reynolds number, as well as incorporating non-linear viscoelasticity in the fluid domain. Our algorithm is designed for computationally efficient simulation of multi-particle suspensions with mixed structure types. The internal force calculation in the solid domain in the IFEM is coupled with a finite volume based incompressible fluid solver, both of which are massively parallelized for distributed memory architectures. We performed extensive case studies to ensure the fidelity of our algorithm. Namely, a series of single particle simulations for capsules, red blood cells, and elastic solid deformable particles were conducted in viscous and viscoelastic media. All of our results are in excellent quantitative agreement with the corresponding reported data in the literature which are based on different simulation platforms. Furthermore, we assess the accuracy of multi-particle simulation of blood suspensions (red blood cells in plasma) with and without platelets. Finally, we present the results of a novel simulation of multiple solid deformable objects in a viscoelastic medium.
Slow flow of a single fluid through a porous medium is well understood on a macroscopic level through Darcys law, a linear relation between flow rate and a combination of pressure differences, viscosity, and gravitational forces. Two-phase flow is complicated by the interface separating the fluids, but understanding of two-dimensional, two-phase flow has been obtained from experiments using transparent cells. In most three-dimensional media, however, visual observation is difficult. Here, we present preliminary results of experiments on a model medium consisting of randomly packed glass spheres, in which one fluorescent liquid invades another. By refractive index matching and scanning with a sheet-shaped laser beam, we obtain slices of the flow patterns, which we combine into three-dimensional pictures. We observe a compact region of invading fluid, surrounded by finger-like protrusions. The compact region becomes more dominant with increasing invader flow rate. The patterns are theoretically analyzed in terms of the interplay between gravitational, viscous, and capillary forces.
We present a theoretical framework for immiscible incompressible two-phase flow in homogeneous porous media that connects the distribution of local fluid velocities to the average seepage velocities. By dividing the pore area along a cross-section transversal to the average flow direction up into differential areas associated with the local flow velocities, we construct a distribution function that allows us not only to re-establish existing relationships between the seepage velocities of the immiscible fluids, but also to find new relations between their higher moments. We support and demonstrate the formalism through numerical simulations using a dynamic pore-network model for immiscible two-phase flow with two- and three-dimensional pore networks. Our numerical results are in agreement with the theoretical considerations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
