Comprehensive Darcy-type law for viscoplastic fluids: Framework

A comprehensive Darcy-type law for viscoplastic fluids is proposed. Different regimes of yield-stress fluid flow in porous media can be categorized based on the Bingham number (i.e., the ratio of the yield stress to the characteristic viscous stress). In a recent study [Chaparian, ], we addressed the yield/plastic limit (infinitely large Bingham number), namely, the onset of flow when the applied pressure gradient is just sufficient to overcome the yield-stress resistance. A purely geometrical universal scale was derived for the nondimensional critical pressure gradient, which was thoroughly validated against computational data. In the present work, we investigate the Newtonian limit (infinitely large pressure difference compared to the yield stress of the fluid—ultralow Bingham number) both theoretically and computationally. We then propose a Darcy-type law applicable across the entire range of Bingham numbers by combining the mathematical models of the yield/plastic and Newtonian limits. The Exhaustive computational data generated in this study (using an augmented Lagrangian method coupled with anisotropic adaptive mesh at the pore scale) confirm the validity of the theoretical proposed law.