Search

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.

A computational approach is proposed for efficient design study of a reducer stent to be percutaneously implanted in enlarged right ventricular outflow tracts (RVOT). The need for such a device is driven by the absence of bovine or artificial valves which could be implanted in these RVOT to replace the absent or incompetent native valve, as is often the case over time after Tetralogy of Fallot repair. Hemodynamics are simulated in the stented RVOT via a reduce order model based on proper orthogonal decomposition (POD), while the artificial valve is modeled as a thin resistive surface. The reduced order model is obtained from the numerical solution on a reference device configuration, then varying the geometrical parameters (diameter) for design purposes. To validate the approach, forces exerted on the valve and on the reducer are monitored, varying with geometrical parameters, and compared with the results of full CFD simulations. Such an approach could also be useful for uncertainty quantification. Device design, percutaneous pulmonary valve replacement, proper orthogonal decomposition (POD), finite element method, blood flow CFD, repaired Tetralogy of Fallot.

International audience; Tetralogy of Fallot is a congenital heart disease characterized over time, after the initial repair, by the absence of a functioning pulmonary valve, which causes regurgitation, and by progressive enlargement of the right ventricle outflow tract (RVOT). Due to this pathological anatomy, available transcatheter valves are usually too small to be deployed there. To avoid surgical valve replacement, an alternative consists in implanting a reducer prior to or in combination with the valve. It has been shown in animal experiments to be promising, but with some limitations. The effect of a percutaneous pulmonary valve reducer on hemodynamics in enlarged RVOT is thus studied by computational modeling. To this aim, blood flow in the RVOT is modeled with CFD coupled to a simplified valve model and 0D downstream models. Simulations are performed in an image-based geometry and boundary conditions tuned to reproduce the pathological flow without the device. Different device designs are built and compared with the initial device-free state, or with the reducer alone. Results suggest that pressure loss is higher for the reducer alone than for the full device, and that the latter successfully restores hemodynamics to a healthy state and induces a more symmetric flow in the pulmonary arteries. Moreover, pressure forces on the reducer and on the valve have the same magnitudes. Migration would occur towards the right ventricle rather than the pulmonary arteries. Results support the thesis that the reducer does not introduce clinically significant pressure gradients, as was found in animal experiments. Such study could help transfer to patients.

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; 3D computational fluid dynamics (CFD) in patient-specific geometries provides complementary insights to clinical imaging, to better understand how heart disease, and the side effects of treating heart disease, affect and are affected by hemodynamics. This information can be useful in treatment planning for designing artificial devices that are subject to stress and pressure from blood flow. Yet, these simulations remain relatively costly within a clinical context. The aim of this work is to reduce the complexity of patient-specific simulations by combining image analysis, computational fluid dynamics and model order reduction techniques. The proposed method makes use of a reference geometry estimated as an average of the population, within an efficient statistical framework based on the currents representation of shapes. Snapshots of blood flow simulations performed in the reference geometry are used to build a POD (Proper Orthogonal Decomposition) basis, which can then be mapped on new patients to perform reduced order blood flow simulations with patient specific boundary conditions. This approach is applied to a data-set of 17 tetralogy of Fallot patients to simulate blood flow through the pulmonary artery under normal (healthy or synthetic valves with almost no backflow) and pathological (leaky or absent valve with backflow) conditions to better understand the impact of regurgitated blood on pressure and velocity at the outflow tracts. The model reduction approach is further tested by performing patient simulations under exercise and varying degrees of pathophysiological conditions based on reduction of reference solutions (rest and medium backflow conditions respectively).