Experimental and numerical determination of Darcy's law for yield stress fluids in porous media

In this work, we studied experimentally and numerically the pressure–flow rate relationship for yield stress fluids in porous media. We developed and validated 3D numerical simulations of the velocity field via a lattice Boltzmann method based on the TRT scheme, and a specific experimental setup allowing yield stress fluids to flow in a closed-loop system to obtain stable rheological properties over a wide range of flow rates (4 decades). The porous medium studied experimentally is a sandstone. The flow properties were also simulated on a regular sphere packing and a random sphere packing as well as on a 3D geometry of a sandstone obtained by x-ray tomography. These different geometries allow highlighting the role of the heterogeneity of the pore structure on the flow properties. All results are expressed as the flow rate ˜Q versus the difference of the pressure gradient to the critical pressure (Δ˜P−Δ˜Pc); Pc defines the pressure below which there is no flow. We observed both numerically and experimentally three specific scaling regimes, already identified by Talon and Bauer [Eur. Phys. J. E 36, 139 (2013)] for a Bingham fluid in 2D porous media. We also evidenced the existence of two critical pressures: the “true” critical pressure ˜Pc defined as the pressure below which there is no flow and the “pseudo” critical pressure threshold ˜P∞c determined by fitting the data in the high-flow-rate regime. We show that the “true” critical pressure is always lower than the “pseudo” one in heterogeneous porous media and can be equal only in the case of regular porous structures. We explain these observations using an energy minimization principle.