Abstract:Mineral dissolution in fractures is largely affected by the flow field within the fracture opening. The flow field maintains transport of the dissolved solutes and influences the concentration gradient and thermodynamic driving force for the mineral dissolution. In this study, we focus on investigating two factors that influence mineral dissolution in fractures: fracture wall roughness and interporosity fluid exchange or the flux through the fracture walls. The classical cubic law for the flow in parallel-plate channels cannot describe the flow in a fracture with rough permeable walls. The wavy, irregular shape of the rough boundaries can result in local flow features and might cause a shift of the overall flow and dissolution regime. At low Reynolds numbers, roughness simply increases the surface area available for mineral dissolution. Moreover, at higher Reynolds numbers, the inertial effects become important and the non-linear flow, flow instabilities, and reverse flow will form. With the creation of eddies and vorticities in the area adjacent to the walls, a flow dead zone forms that reduces the reactive surface area. In the presence of interporosity fluid exchange between the fracture and its surrounding rock, the non-zero fluid velocity on the fracture walls influences the local concentration gradient. In addition, the newly introduced fluid may facilitate or suppress the fracture wall mineral dissolution depending on its composition. This work highlights the compound effect of roughness and fluid flux through the walls on fracture mineral dissolution. For this purpose, the asymptotic solutions of the steady-state Navier–Stokes equations with non-zero velocity on the borders are used to determine the velocity field within the fracture opening. The quadratic and cubic corrections to Darcy’s law are expected as a result of the wall roughness and the flow through the walls. The flow field is coupled with a transport module and a geochemical model (PHREEQC). As a test case, we investigate calcite dissolution in a single fracture at a different influx rate. For each flow rate, a flat fracture and two sets of rough profiles with and without the permeable walls were compared. The simulation results showed that compared to impermeable walls, the pervious walls result in a non-uniform non-periodic mineral dissolution along the fracture, which is more focused at the inlet. At low Reynolds numbers, the hotspots of dissolution are slightly shifted from the smaller cross sectional area to the larger cross sectional area. The effective reaction rate for mineral dissolution of fracture walls increases with an increase in fracture surface roughness. It is shown that at low Reynolds numbers, the permeable fracture walls can improve the effective reaction rate significantly. At relatively larger Reynolds numbers, the impact of flux through the walls on the effective reaction rate is less important, but for some roughness profiles, still not negligible.

Direct numerical simulations of immiscible two-phase flow in rough fractures: Impact of wetting film resolution

Direct Numerical Simulation of Channel Flow with Real Surface Roughness Using a Ghost Cell Immersed Boundary Method

Computer simulation of two-phase immiscible fluid motion in unsaturated complex fractures using a volume of fluid method

Finite volume simulations of particle-laden viscoelastic fluid flows: application to hydraulic fracture processes

Numerical Simulation of Unsaturated Flow in Single Fracture Considering the Water Exchange Between Fracture and Matrix

Reduced‐Order Models Unravel the Joint Impact of Aperture Heterogeneity and Shear‐Thinning Rheology on Fracture‐Scale Flow Metrics

Simulation of hydro-mechanically coupled processes in rough rock fractures using an immersed boundary method and variational transfer operators

Nonlinear fluid flow through three-dimensional rough fracture networks: Insights from 3D-printing, CT-scanning, and high-resolution numerical simulations

A model for saturated–unsaturated flow with fractures acting as capillary barriers

Numerical Simulation of Fracture Flow Interaction Based on Discrete Fracture Model

The role of wetting on the flow of two immiscible fluids in porous media

Modelling of hydro-mechanical processes in heterogeneous fracture intersections using a fictitious domain method with variational transfer operators

Experimental and numerical analysis of flow through a natural rough fracture subject to normal loading

Effect of wettability on two-phase quasi-static displacement: Validation of two pore scale modeling approaches

Modeling of surface tension and contact angles with smoothed particle hydrodynamics

Numerical Simulations of Viscous Fingering in Fractured Porous Media

Pore-scale investigation of immiscible displacement in rough fractures

Numerical simulation of water injection into layered fractured carbonate reservoir analogs

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

Dissipative Particle Dynamics Simulation of Fluid Motion Through an Unsaturated Fracture and Fracture Junction

Semi-analytical Models of Mineral Dissolution in Rough Fractures with Permeable Walls