Abstract:The dynamic two-phase flow in porous media was theoretically developed based on mass, momentum conservation, and fundamental constitutive relationships for simulating immiscible fluid-fluid retention behavior and seepage in the natural geomaterial. The simulation of transient two-phase flow seepage is, therefore, dependent on both the hydraulic boundaries applied and the immiscible fluid-fluid retention behavior experimentally measured. Many previous studies manifested the velocity-dependent capillary pressure–saturation relationship (Pc-S) and relative permeability (Kr-S). However, those works were experimentally conducted on a continuum scale. To discover the dynamic effects from the microscale, the Computational Fluid Dynamic (CFD) is usually adopted as a novel method. Compared to the conventional CFD methods solving Naiver–Stokes (NS) equations incorporated with the fluid phase separation schemes, the two-phase Lattice Boltzmann Method (LBM) can generate the immiscible fluid-fluid interface using the fluid-fluid/solid interactions at a microscale. Therefore, the Shan–Chen multiphase multicomponent LBM was conducted in this study to simulate the transient two-phase flow in porous media. The simulation outputs demonstrate a preferential flow path in porous media after the non-wetting phase fluid is injected until, finally, the void space is fully occupied by the non-wetting phase fluid. In addition, the inter-relationships for each pair of continuum state variables for a Representative Elementary Volume (REV) of porous media were analyzed for further exploring the dynamic nonequilibrium effects. On one hand, the simulating outcomes reconfirmed previous findings that the dynamic effects are dependent on both the transient seepage velocity and interfacial area dynamics. Nevertheless, in comparison to many previous experimental studies showing the various distances between the parallelly dynamic and static Pc-S relationships by applying various constant flux boundary conditions, this study is the first contribution showing the Pc-S striking into the nonequilibrium condition to yield dynamic nonequilibrium effects and finally returning to the equilibrium static Pc-S by applying various pressure boundary conditions. On the other hand, the flow regimes and relative permeability were discussed with this simulating results in regards to the appropriateness of neglecting inertial effects (both accelerating and convective) in multiphase hydrodynamics for a highly pervious porous media. Based on those research findings, the two-phase LBM can be demonstrated to be a powerful tool for investigating dynamic nonequilibrium effects for transient multiphase flow in porous media from the microscale to the REV scale. Finally, future investigations were proposed with discussions on the limitations of this numerical modeling method.

Dynamic effects on capillary pressure–Saturation relationships for two‐phase porous flow: Implications of temperature

Determination of dynamic capillary effect on two-phase flow in porous media: A perspective from various methods

Bundle-of-Tubes Model for Calculating Dynamic Effects in the Capillary-Pressure- Saturation Relationship

Dimensional analysis of two-phase flow including a rate-dependent capillary pressure–saturation relationship

Effects of dynamic capillarity on the shear strength of sandy soils during transient two‐phase flow: Insights from non‐equilibrium triaxial simulations

Uniqueness of Specific Interfacial Area–Capillary Pressure–Saturation Relationship Under Non-Equilibrium Conditions in Two-Phase Porous Media Flow

Semi-analytical Model to Predict Dynamic Capillary Pressure–Saturation Relationship for Flows in Heterogeneous Porous Media

Pore-Scale Simulations of Simultaneous Steady-State Two-Phase Flow Dynamics Using a Lattice Boltzmann Model: Interfacial Area, Capillary Pressure and Relative Permeability

Application of the Darcy and Richards equations for modelling of water capillary rise in building materials

A Comprehensive Review of Factors Affecting Dynamic Capillary Effect in Two-Phase Flow

On the relationship between capillary pressure, saturation, and interfacial area for three-phase flow in water-wet porous media

Dynamic simulation of unsteady Casson–Carreau fluid flow considering temperature-dependent variable-properties and thermo-solutal mixed convection over flat and slendering sheets

The average equilibrium capillary pressure-saturation relationship in two-phase flow in porous media

Pore network modeling of the effects of viscosity ratio and pressure gradient on steady-state incompressible two-phase flow in porous media

Laboratory Experiments and Simulations on the Significance of Non-Equilibrium Effect in the Capillary Pressure-Saturation Relationship

Discovery of Dynamic Two-Phase Flow in Porous Media Using Two-Dimensional Multiphase Lattice Boltzmann Simulation

Pore-Scale Modeling of the Effect of Wettability on Two-Phase Flow Properties for Newtonian and Non-Newtonian Fluids

Homogenization of Two-Phase Flow in Porous Media From Pore to Darcy Scale: A Phase-Field Approach

Thermodynamic basis of capillary pressure in porous media

Effective Rheology of Two-Phase Flow in Three-Dimensional Porous Media: Experiment and Simulation

Capillary fluctuations and energy dynamics for flow in porous media