Effects of inertia and directionality on flow and transport in a rough asymmetric fracture

[1] Two-dimensional flow and transport simulations are conducted for a 15-cm-long fracture mapped via X-ray computed tomography. The simulations consider either Navier-Stokes equation (NSE) based flow or Stokes equation (SE) based flow and are run for opposing directions. NSE and SE solutions deviate at larger bulk velocity, and the errors are sensitive to fracture geometry. Transport in all cases is non-Fickian, owing mainly to the presence of a large eddy, and exhibits power law residence time distributions (RTDs). The tailing is more persistent at higher Reynolds numbers (Re), with the exponents of the power law RTDs related to Re via a power function, and becomes asymptotic at higher Re. The late-time portion of the RTD changes with flow direction and is also sensitive to whether the flow is represented by the NSE or the SE. The sensitivity of transport to flow direction and formulation is primarily driven by varying eddy geometry under the different conditions. Our study opens the path to developing robust and mechanistic continuum-scale models of solute transport in fractured geologic media.

[1]  C. Tsang,et al.  Flow channeling in heterogeneous fractured rocks , 1998 .

[2]  Joel P. L. Johnson,et al.  Mixing at fracture intersections: Influence of channel geometry and the Reynolds and Peclet Numbers , 2001 .

[3]  Jean Schmittbuhl,et al.  Flow enhancement of a rough fracture , 2000 .

[4]  S. Ge,et al.  The effect of surface geometry on fracture permeability: A case study using a sinusoidal fracture , 1998 .

[5]  B. Kueper,et al.  Evaluation of cubic law based models describing single‐phase flow through a rough‐walled fracture , 2004 .

[6]  Brian Berkowitz,et al.  Computing “Anomalous” Contaminant Transport in Porous Media: The CTRW MATLAB Toolbox , 2005, Ground water.

[7]  P. Reimus,et al.  Using multiple experimental methods to determine fracture/matrix interactions and dispersion of nonreactive solutes in saturated volcanic tuff , 2000 .

[8]  B. Berkowitz Characterizing flow and transport in fractured geological media: A review , 2002 .

[9]  David J. Brush,et al.  Fluid flow in synthetic rough‐walled fractures: Navier‐Stokes, Stokes, and local cubic law simulations , 2003 .

[10]  Christopher C. Pain,et al.  Non-linear regimes of fluid flow in rock fractures , 2004 .

[11]  Georg Kosakowski,et al.  Anomalous transport of colloids and solutes in a shear zone. , 2004, Journal of contaminant hydrology.

[12]  Joel Koplik,et al.  Flow of power-law fluids in self-affine fracture channels. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  M. Bayani Cardenas,et al.  Hydrodynamics of coupled flow above and below a sediment-water interface with triangular bedforms , 2007 .

[14]  B. Berkowitz,et al.  Three‐dimensional flow measurements in rock fractures , 1999 .

[15]  W. Durham,et al.  An experimental exploration of the transport and capture of abiotic colloids in a single fracture , 2000 .

[16]  M. Becker,et al.  Tracer transport in fractured crystalline rock: Evidence of nondiffusive breakthrough tailing , 2000 .

[17]  Andrea Cortis Péclet-dependent memory kernels for transport in heterogeneous media. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[18]  Young-Jin Park,et al.  Analytical solutions for solute transfer characteristics at continuous fracture junctions , 1999 .

[19]  G. Bodvarsson,et al.  Evidence of Multi-Process Matrix Diffusion in a Single Fracture from a Field Tracer Test , 2006 .

[20]  M. Bayani Cardenas,et al.  Dunes, turbulent eddies, and interfacial exchange with permeable sediments , 2007 .

[21]  Arturo A. Keller,et al.  Effect of fracture aperture variations on the dispersion of contaminants , 1999 .

[22]  W. Durham,et al.  Fracture flow simulation using a finite-difference lattice Boltzmann method. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[23]  B. Armaly,et al.  Experimental and theoretical investigation of backward-facing step flow , 1983, Journal of Fluid Mechanics.

[24]  C. Chrysikopoulos,et al.  Transport of Neutrally Buoyant and Dense Variably Sized Colloids in a Two-Dimensional Fracture with Anisotropic Aperture , 2003 .

[25]  Jens T. Birkholzer,et al.  Continuous time random walk analysis of solute transport in fractured porous media , 2008 .

[26]  Stephen R. Brown,et al.  The effect of anisotropic surface roughness on flow and transport in fractures , 1991 .

[27]  M. Bayani Cardenas,et al.  Three‐dimensional vortices in single pores and their effects on transport , 2008 .

[28]  R. Ketcham,et al.  Acquisition, optimization and interpretation of X-ray computed tomographic imagery: applications to the geosciences , 2001 .

[29]  P. Kitanidis,et al.  Prediction of single phase transport parameters in a variable aperture fracture , 1995 .

[30]  Giovanni Grasselli,et al.  Trapping zones: The effect of fracture roughness on the directional anisotropy of fluid flow and colloid transport in a single fracture , 2006 .

[31]  Richard A. Ketcham,et al.  Navier‐Stokes flow and transport simulations using real fractures shows heavy tailing due to eddies , 2007 .

[32]  M. Sahimi,et al.  Fluid flow and conduction in two-dimensional fractures with rough, self-affine surfaces: A comparative study , 2003 .

[33]  Li Yuan-hui,et al.  Diffusion of ions in sea water and in deep-sea sediments , 1974 .

[34]  H. K. Moffatt Viscous and resistive eddies near a sharp corner , 1964, Journal of Fluid Mechanics.

[35]  Harold Auradou,et al.  Permeability anisotropy induced by the shear displacement of rough fracture walls , 2005 .

[36]  Stephen R. Brown,et al.  Applicability of the Reynolds Equation for modeling fluid flow between rough surfaces , 1995 .

[37]  Matthew D. Jackson,et al.  Characterization of multiphase electrokinetic coupling using a bundle of capillary tubes model , 2008 .