Influence of phase connectivity on the relationship among capillary pressure, fluid saturation, and interfacial area in two-fluid-phase porous medium systems

James E. McClure, Mark A. Berrill, William G. Gray, Cass T. Miller

Research output: Contribution to journalArticlepeer-review

48 Scopus citations

Abstract

Multiphase flows in porous medium systems are typically modeled at the macroscale by applying the principles of continuum mechanics to develop models that describe the behavior of averaged quantities, such as fluid pressure and saturation. These models require closure relations to produce solvable forms. One of these required closure relations is an expression relating the capillary pressure to fluid saturation and, in some cases, other topological invariants such as interfacial area and the Euler characteristic (or average Gaussian curvature). The forms that are used in traditional models, which typically consider only the relationship between capillary pressure and saturation, are hysteretic. An unresolved question is whether the inclusion of additional morphological and topological measures can lead to a nonhysteretic closure relation. Relying on the lattice Boltzmann (LB) method, we develop an approach to investigate equilibrium states for a two-fluid-phase porous medium system, which includes disconnected nonwetting phase features. A set of simulations are performed within a random close pack of 1964 spheres to produce a total of 42 908 distinct equilibrium configurations. This information is evaluated using generalized additive models to quantitatively assess the degree to which functional relationships can explain the behavior of the equilibrium data. The variance of various model estimates is computed, and we conclude that, except for the limiting behavior close to a single fluid regime, capillary pressure can be expressed as a deterministic and nonhysteretic function of fluid saturation, interfacial area between the fluid phases, and the Euler characteristic. To our knowledge, this work is unique in the methods employed, the size of the data set, the resolution in space and time, the true equilibrium nature of the data, the parametrizations investigated, and the broad set of functions examined. The conclusion of essentially nonhysteretic behavior provides support for an evolving class of two-fluid-phase flow in porous medium systems models.

Original languageEnglish
Article number033102
JournalPhysical Review E - Statistical, Nonlinear, and Soft Matter Physics
Volume94
Issue number3
DOIs
StatePublished - Sep 2 2016

Funding

This work was supported by Army Research Office Grant No. W911NF-14-1-02877, Department of Energy Grant No. DE-SC0002163, and National Science Foundation Grant No. 1619767. An award of computer time was provided by the Department of Energy INCITE program. This research also used resources of the Oak Ridge Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract No. DE-AC05-00OR22725

FundersFunder number
National Science Foundation1619767
U.S. Department of EnergyDE-SC0002163
Army Research OfficeW911NF-14-1-02877
Office of ScienceDE-AC05-00OR22725

    Fingerprint

    Dive into the research topics of 'Influence of phase connectivity on the relationship among capillary pressure, fluid saturation, and interfacial area in two-fluid-phase porous medium systems'. Together they form a unique fingerprint.

    Cite this