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

• Published in: Transport in porous media 2017-08, Vol.119 (1), p.77-94
• Main Authors: Sinha, Santanu, Bender, Andrew T., Danczyk, Matthew, Keepseagle, Kayla, Prather, Cody A., Bray, Joshua M., Thrane, Linn W., Seymour, Joseph D., Codd, Sarah L., Hansen, Alex
• Format: Article
• Language: English
• Subjects: Beads; Civil Engineering; Classical and Continuum Physics; Computational fluid dynamics; Computer simulation; Deionization; Earth and Environmental Science; Earth Sciences; Flow velocity; Fluid flow; Geotechnical Engineering & Applied Earth Sciences; Hydrogeology; Hydrology/Water Resources; Industrial Chemistry/Chemical Engineering; Mathematical models; Newtonian fluids; Porous materials; Porous media; Pressure drop; Rheological properties; Rheology; Steady state; Stress concentration; Three dimensional flow; Transportation networks; Two phase flow; Wetting
• ISSN: 0169-3913; 1573-1634
• DOI: 10.1007/s11242-017-0874-4

We present an experimental and numerical study of immiscible two-phase flow of Newtonian fluids in three-dimensional (3D) porous media to find the relationship between the volumetric flow rate ( Q ) and the total pressure difference ( Δ P ) in the steady state. We show that in the regime where capillary forces compete with the viscous forces, the distribution of capillary barriers at the interfaces effectively creates a yield threshold ( P t ), making the fluids reminiscent of a Bingham viscoplastic fluid in the porous medium. In this regime, Q depends quadratically on an excess pressure drop ( Δ P - P t ). While increasing the flow rate, there is a transition, beyond which the overall flow is Newtonian and the relationship is linear. In our experiments, we build a model porous medium using a column of glass beads transporting two fluids, deionized water and air. For the numerical study, reconstructed 3D pore networks from real core samples are considered and the transport of wetting and non-wetting fluids through the network is modeled by tracking the fluid interfaces with time. We find agreement between our numerical and experimental results. Our results match with the mean-field results reported earlier.