A CSF-SPH method for simulating drainage and imbibition at pore-scale resolution while tracking interfacial areas

Research paper by Rakulan Sivanesapillai, Nadine Falkner, Alexander Hartmaier, Holger Steeb

Indexed on: 19 Mar '16Published on: 28 Aug '15Published in: Advances in Water Resources


We present a conservative smoothed particle hydrodynamics (SPH) model to study the flow of multiple, immiscible fluid phases in porous media using direct pore-scale simulations. Particular focus is put on continuously tracking the evolution of interfacial areas, which are considered to be important morphological quantities affecting multiphase transport in porous media. In addition to solving the Navier–Stokes equations, the model accounts for the effects of capillarity at interfaces and contact lines. This is done by means of incorporating the governing interfacial mass and momentum balances using the continuum surface force (CSF) method, thus rendering model calibration routines unnecessary and minimizing the set of constitutive and kinematic assumptions. We address the application of boundary conditions at rigid solid surfaces and study the predictive capability of the model as well as optimal choices for numerical parameters using an extensive model validation procedure. We demonstrate the applicability of the model to simulate multiphase flows involving partial wettability, dynamic effects, large density ratios (up to 1000), large viscosity ratios (up to 100), as well as fragmentation and coalescence of fluid phases. The model is used to study the evolution of fluid–fluid interfacial areas during saturation-controlled primary drainage and main imbibition of heterogeneous pore spaces at low capillary numbers. A variety of pore-scale effects, such as wetting phase entrapment and fragmentation due to snap-off, are observed. Specific fluid–fluid interfacial area is observed to monotonically increase during primary drainage and hysteretic effects are apparent during main imbibition.

Figure 10.1016/j.advwatres.2015.08.012.0.jpg
Figure 10.1016/j.advwatres.2015.08.012.1.jpg
Figure 10.1016/j.advwatres.2015.08.012.2.jpg
Figure 10.1016/j.advwatres.2015.08.012.3.jpg
Figure 10.1016/j.advwatres.2015.08.012.4.jpg
Figure 10.1016/j.advwatres.2015.08.012.5.jpg
Figure 10.1016/j.advwatres.2015.08.012.6.jpg
Figure 10.1016/j.advwatres.2015.08.012.7.jpg