International audienceMultiphase 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 specifc volume-averaged equations and therefore specifc computational fluid dynamics (CFD) tools. In this work, we develop a toolbox for modeling multiphase ow in porous media with OpenFOAM®, an opensource 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, 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.