Stochastic 3D Navier‐Stokes Flow in Self‐Affine Fracture Geometries Controlled by Anisotropy and Channeling

This study presents a probabilistic analysis of 3D Navier‐Stokes (NS) fluid flow through 30 randomly generated sheared fractures with equal roughness properties (Hurst exponent = 0.8). The results of numerous 3D NS realizations are compared with the highly simplified local cubic law (LCL) solutions regarding flow orientations and regimes. The transition between linear and nonlinear flow conditions cannot be described with a generally valid critical Reynolds number urn:x-wiley:00948276:media:grl62319:grl62319-math-0001, but rather depends on the individual fracture's void geometry. Over 10% reduction in flow is observed for increased global Re (>100) due to the increasing impact of nonlinear conditions. Furthermore, the fracture geometry promotes flow anisotropy and the formation of channels. Flow perpendicular to the shearing leads to increased channeling and fluid flow (∼40% higher) compared to flow parallel to the shearing. In the latter case, dispersed flow and irregular flow paths cause a reduction of LCL validity.