Search

International audience; An image-based effective medium approximation (EMA) is developed so as to permit very fast transport properties evaluations of 3D porous media. From an image-based porous network (IBPN) built upon digital image processing of 3D binary images, we focus on throat’s local geometrical properties at the pore scale, for being the most sensible structural units which build up the local pressure. This approach is a 3D image–based extension of the critical point approach proposed in 2D fractures. We show, from analyzing various core rock samples available in the literature, that the asymptotic assumptions associated with the preeminence of critical points in throats are indeed geometrically relevant. We then describe how the image-based EMA evaluated from the conductances computed from the discrete IBPN can be reliably evaluated. The proposed method is evaluated upon the estimation of core sample permeability from binarized image obtained using X-ray tomography. Since it combines digital image treatments with statistical data post-processing without the need of computational fluid dynamics (CFD) computation, it is extremely cost efficient. The results are compared with a micro-scale Stokes flow computation in various rock samples. The sensitivity to the pore discretization also is discussed and illustrated.

The IMplicit Pressure Explicit Saturation (IMPES) method is a prevalent way to simulate multiphase flows in porous media. The numerical stability of this sequential method implies limitations on the time step, which may depend on the flow regime studied. In this note, three stability criteria related to the IMPES method, that differ in their construction on the observed variables, are compared on homogeneous and heterogeneous configurations for different two-phase flow regimes (viscous/capillary/gravitational). This highlights that there is no single optimal criterion always ensuring stability and efficiency. For capillary dominated flows, the Todd's condition is the most efficient one, while the standard Coat condition should be preferred for viscous flows. When gravity effects are present, Coat's condition must be restricted, but remains more efficient than the Todd's condition.

International audience; Multiphase flow in porous media provides a wide range of applications: from the environmental understanding (aquifer, site-pollution) to industrial process improvements (oil production, waste management). Modeling of such flows involve specific volume-averaged equations and therefore specific computational fluid dynamics (CFD) tools. In this work, we develop a toolbox for modeling multiphase flow in porous media with OpenFOAM®, an open-source platform for CFD. The underlying idea of this approach is to provide an easily adaptable tool that can be used in further studies to test new mathematical models or numerical methods. The package provides the most common effective properties models of the literature (relative permeability, capillary pressure) and specific boundary conditions related to porous media flows. To validate this package, a solvers based on the IMplicit Pressure Explicit Saturation (IMPES) methodare developed in the toolbox. The numerical validation is performed by comparison with analytical solutions on academic cases. Then, a satisfactory parallel efficiency of the solver is shown on a more complex configuration.

In this work, several boundary value problems used to numerically evaluate the absolute permeability tensors of porous media using core-scale images are compared and discussed. The various configurations differ by the type of boundary conditions used to compute the flow at the micro-scale. The issue is the ability of the method to capture anisotropy correctly and to avoid possible percolation artifacts. This study is carried on two-dimensional synthetic, isotropic or anisotropic, porous media, that are chosen to illustrate the various difficulties mentioned above. A new method is proposed which consists in embedding the porous medium in question in a homogenized one. Using an iterative optimization procedure on the surrounding permeability, the method determines the absolute permeability tensor of the original medium. The equivalent permeability tensor that minimizes the effect on the surrounding porous medium is, unlike that of classical methods, de facto symmetrical due to the use of periodic boundary conditions and exhibits significantly lower permeabilities. The way in which non-diagonal terms of the permeability tensor are obtainedwith the various methods are thoroughly discussed.

We present a theoretical and numerical study of three-dimensional pulsatile confined flow between two rigid horizontal surfaces separated by an average gap h, and having three-dimensional wavy shapes with arbitrary amplitude σh where σ ~ O(1), but long-wavelength variations λ, with h/λ ≪ 1. We are interested in pulsating flows with moderate inertial effect arising from the Reynolds stress due to the cavity non-parallelism. We analyse the inertial steady-streaming and the second harmonic flows in a lubrication approximation. The dependence of the three-dimensional velocity field in the transverse direction is analytically obtained for arbitrary Womersley numbers and possibly overlapping Stokes layers. The horizontal dependence of the flow is solved numerically by computing the first two pressure fields of an asymptotic expansion in the small inertial limit. We study the variations of the flow structure with the amplitude, the channel's wavelength and the Womersley number for various families of three-dimensional channels. The steady-streaming flow field in the horizontal plane exhibits a quadrupolar vortex, the size of which is adjusted to the cavity wavelength. When increasing the wall amplitude, the wavelengths characterizing the channel or the Womersley number, we find higher-order harmonic flow structures, the origin of which can either be inertially driven or geometrically induced. When some of the channel symmetries are broken, a steady-streaming current appears which has a quadratic dependence on the pressure drop, the amplitude of which is linked to the Womersley number.

The modeling of blood flows confined in micro-channels or micro-capillary beds depends on the interactions between the cell-phase, plasma and the complex geometry of the network. In the case of capillaries or channels having a high aspect ratio (their longitudinal size is much larger than their transverse one), this modeling is much simplified from the use of a continuous description of fluid viscosity as previously proposed in the literature. Phase separation or plasma skimming effect is a supplementary mechanism responsible for the relative distribution of the red blood cell’s volume density in each branch of a given bifurcation. Different models have already been proposed to connect this effect to the various hydrodynamics and geometrical parameters at each bifurcation. We discuss the advantages and drawbacks of these models and compare them to an alternative approach for modeling phase distribution in complex channels networks. The main novelty of this new formulation is to show that albeit all the previous approaches seek for a local origin of the phase segregation phenomenon, it can arise from a global non-local and nonlinear structuration of the flow inside the network. This new approach describes how elementary conservation laws are sufficient principles (rather than the complex parametric models previously proposed) to provide non local phase separation. Spatial variations of the hematocrit field thus result from the topological complexity of the network as well as nonlinearities arising from solving a new free boundary problem associated with the flux and mass conservation. This network model approach could apply to model blood flow distribution either on artificial micro-models, micro-fluidic networks, or realistic reconstruction of biological micro-vascular networks.

International audience; In this work, a complete work flow from pore-scale imaging to absolute permeability determination is described and discussed. Two specific points are tackled, concerning (1) the mesh refinement for a fixed image resolution and (2) the impact of the determination method used. A key point for this kind of approach is to work on enough large samples to check the representativity of the obtained evaluations, which requires efficient parallel capabilities. Image acquisition and processing are realized using a commercial micro-tomograph. The pore-scale flows are then evaluated using the finite volume method implemented in the open-source platform OpenFOAM®. For this numerical method, the influence of the different aspects mentioned above are studied. Moreover, the parallel efficiency is also tested and discussed. We observe that the level of mesh refinement has a non-negligible impact on permeability tensor. Moreover, increasing the refinement level tends to reduce the gap between the methods of computational measurements. The increase in computation time with the mesh is balanced with the good parallel efficiency of the platform.

International audience; Cerebral blood flow (CBF) is the most common parameter for the quantification of brain's function. Literature data indicate a widespread dispersion of values that might be related to some differences in the measurement conditions that are not properly taken into account in CBF evaluation. Using recent high-resolution imaging of the complete cortical microvasculature of primate brain, we perform extensive numerical evaluation of the cerebral perfusion. We show that blood perfusion associated with intravascular tracers should be normalized by the surface of the voxel rather than by its volume and we consistently test this result on the available literature data.

International audience; Cerebral micro-vascular networks control the blood pressure distribution when considering in vitro blood rheology models. Blood rheology is complex and non-linear. In small vessels, the effective viscosity variations are important due to red blood cells packing in capillaries, the so-called Fåhræus-Lindquist effect, whilst concomitantly phase segregation appears in bifurcations. Direct numerical simulations of different non-linear rheological models of the blood are performed on realistic three-dimensional micro-vascular networks. These simulations exhibit two significant results. First, various rheological models lead to very similar pressure distribution over the whole range of physiologically relevant hematocrits. Secondly, different models for phase segregation lead to very distinct hematocrit distributions in the micro-vacular network. Nevertheless, the hematocrit distribution very weakly affects the pressure distribution. Hence, our results suggest that the micro-vacular network structure mainly controls the pressure distribution in micro-circulation, whilst the effect of hematocrit distribution is weak.