Development of depth-averaged nonlinear flow model for fractures considering heterogeneity and significant inertial effects

Fluid flow through fractures is central to many engineering applications and geophysical processes. Although the intertwined effects of aperture heterogeneity and inertial forces on fluid flow have been interrogated in last few decades, it remains challenging to accurately predict fracture flow field at high pressure gradient conditions when inertial effects are significant. To address this challenge, we rigorously derived a depth-averaged nonlinear flow model by introducing a quadratic velocity term into the linear flow model (i.e., Local Cubic Law that describes a linear relationship between pressure gradient and fluid flux) at the mechanistic level. The supporting laboratory experiments and direct numerical simulations through six rough-walled fractures with varying roughness further confirmed the robustness of the proposed model. We found that the predictive model outperforms existing alternatives with acceptable relative errors in mean velocity (