In multiphase flow in porous media the consistent pore to Darcy scale description of two-fluid flow processes has been a long-standing challenge. Immiscible displacement processes occur at the scale of individual pores. However, the larger scale behavior is described by phenomenological relationships such as relative permeability, which typically uses only fluid saturation as a state variable. As a consequence pore scale properties such as contact angle cannot be directly related to Darcy scale flow parameters. Advanced imaging and computational technologies are closing the gap between the pore and Darcy scale, supporting the development of new theory. We utilize fast x-ray microtomography to observe pore-scale two-fluid configurations during immiscible flow and initialize lattice Boltzmann simulations that demonstrate that the mobilization of disconnected nonwetting phase clusters can account for a significant fraction of the total flux. We show that fluid topology can undergo substantial changes during flow at constant saturation, which is one of the underlying causes of hysteretic behavior. Traditional assumptions about fluid configurations are therefore an oversimplification. Our results suggest that the role of fluid connectivity cannot be ignored for multiphase flow. On the Darcy scale, fluid topology and phase connectivity are accounted for by interfacial area and Euler characteristic as parameters that are missing from our current models.

Multiphase flows in porous media are important in many natural and industrial processes. Pore-scale models for multiphase flows have seen rapid development in recent years and are becoming increasingly useful as predictive tools in both academic and industrial applications. However, quantitative comparisons between different pore-scale models, and between these models and experimental data, are lacking. Here, we perform an objective comparison of a variety of state-of-the-art pore-scale models, including lattice Boltzmann, stochastic rotation dynamics, volume-of-fluid, level-set, phase-field, and pore-network models. As the basis for this comparison, we use a dataset from recent microfluidic experiments with precisely controlled pore geometry and wettability conditions, which offers an unprecedented benchmarking opportunity. We compare the results of the 14 participating teams both qualitatively and quantitatively using several standard metrics, such as fractal dimension, finger width, and displacement efficiency. We find that no single method excels across all conditions and that thin films and corner flow present substantial modeling and computational challenges.

No abstract

A model formulated in terms of both conservation and kinematic equations for phases and interfaces in two-fluid-phase flow in a porous medium system is summarized. Macroscale kinematic equations are derived as extensions of averaging theorems and do not rely on conservation principles. Models based on both conservation and kinematic equations can describe multiphase flow with varying fidelity. When only phase-based equations are considered, a model similar in form to the traditional model for twofluid-phase flow results. When interface conservation and kinematic equations are also included, a novel formulation results that naturally includes evolution equations that express dynamic changes in fluid saturations, pressures, the capillary pressure, and the fluid-fluid interfacial area density in a two-fluid-system. This dynamic equation set is unique to this work, and the importance of the modeled physics is shown through both microfluidic experiments and high-resolution lattice Boltzmann simulations. The validation work shows that the relaxation of interface distribution and shape toward an equilibrium state is a slow process relative to the time scale typically allowed for a system to approach an apparent equilibrium state based upon observations of fluid saturations and external pressure measurements. Consequently, most pressuresaturation data intended to denote an equilibrium state are likely a sampling from a dynamic system undergoing changes of interfacial curvatures that are not typically monitored. The results confirm the importance of kinematic analysis in combination with conservation equations for faithful modeling of system physics.

Models that describe two-fluid flow in porous media suffer from a widely-recognized problem that the constitutive relationships used to predict capillary pressure as a function of the fluid saturation are non-unique, thus requiring a hysteretic description. As an alternative to the traditional perspective, we consider a geometrical description of the capillary pressure, which relates the average mean curvature, the fluid saturation, the interfacial area between fluids, and the Euler characteristic. The state equation is formulated using notions from algebraic topology and cast in terms of measures of the macroscale state. Synchrotron-based X-ray micro-computed tomography (µCT) and highresolution pore-scale simulation is applied to examine the uniqueness of the proposed relationship for six different porous media. We show that the geometric state function is able to characterize the microscopic fluid configurations that result from a wide range of simulated flow conditions in an averaged sense. The geometric state function can serve as a closure relationship within macroscale models to effectively remove hysteretic behavior attributed to the arrangement of fluids within a porous medium. This provides a critical missing component needed to enable a new generation of higher fidelity models to describe two-fluid flow in porous media. PACS numbers: 92.40Cy, 92.40.Kf, 91.60.Tn, 91.60.Fe