14-17 May 2018
New Orleans
US/Central timezone

Numerical Simulation Study on Pore Scale Seepage of Porous Media Based on Finite Volume Method

15 May 2018, 18:30
New Orleans

New Orleans

Poster MS 2.10: Advanced finite-volume methods for flow and transport in porous media Poster 2


Ms Shuo Wen (Beijing Institute of Petrochemical Engineering)


Abstract: The numerical simulation of pore scale seepage in porous media is of great significance for the development of new energy sources, such as shale gas and geothermal energy. In the past few years, many scholars have developed the numerical method including finite volume method (FVM), Lattice Boltzmann Method (LBM) and molecular dynamics to achieve this simulation in different levels. It is easy for us to understand the basic idea of FVM, and then we can give a direct physical explanation for the above phenomena. FVM has stronger advantages in CFD (computational fluid dynamics) since the conservation law is meet in the whole computational domain. Based on the FVM, a large number of models on multi-phase flow and non-Newtonian fluid have been developed and it has showed more advantages in the aspect of solving complex pore flow problems. However, there are two shortcomings in the process of using FVM: (1) Due to the complexity of the pore geometry in the porous media and irregularly curved pore walls, the quality of the generated grid is not high; (2) It is relatively difficult to deal with a large number of pore boundaries in the porous media. To cure the above problems, a simple processing method is proposed as follows in this paper: On the one hand, in porous medium, both the penetrating space and solid frame are regarded as the fluid region, while, to describe solid boundary approximately, the viscosity of the solid skeleton structure is set to infinite. On the other hand, the whole region is discretized through the Cartesian orthogonal grid and the Immersed Boundary Method(IBM) is adopted to describe the smooth and complexed boundaries between the penetrating space and solid frame. In this paper, the proposed method stated above is conducted on OpenFOAM by the finite volume method in which the convection term is discreted in the QUICK scheme. After that, the corresponding systematic evaluations of the proposed method are performed and results illustrate that the method is of great accuracy and robustness. It can be used to study the flow in complex porous media.


Khan I, Aidun C K. Modeling the macroscopic behavior of saturated deformable porous media using direct numerical simulations[J]. International Journal of Multiphase Flow, 2015, 71:74-82.
M. Mößner, R. Radespiel. Flow Simulations over Porous Media – Comparisons with Experiments[J]. Computers & Fluids, 2017.
Mansour M A, El-Aziz M M A, Mohamed R A, et al. Numerical Simulation of Natural Convection in Wavy Porous Cavities Under the Influence of Thermal Radiation Using a Thermal Non-equilibrium Model[J]. Transport in Porous Media, 2011, 86(2):585-600.
De S, Kuipers J A M, Peters E A J F, et al. Viscoelastic flow simulations in random porous media[J]. Journal of Non-Newtonian Fluid Mechanics, 2017, 248:50-61.

Acceptance of Terms and Conditions Click here to agree

Primary authors

Ms Shuo Wen (Beijing Institute of Petrochemical Engineering) Mr Dongxu Han (School of Mechanical Engineering, Beijing Institute of Petrochemical Technology, Beijing 102617, China) Mr Yajun Deng (China University of Petroleum,Beijing) Prof. Bo Yu (School of Mechanical Engineering, Beijing Institute of Petrochemical Technology, Beijing 102617, China) Prof. Dongliang Sun (School of Mechanical Engineering, Beijing Institute of Petrochemical Technology, Beijing 102617, China)

Presentation Materials

There are no materials yet.