Numerical modelling of flow and transport in rough fractures

Abstract Simulation of flow and transport through rough walled rock fractures is investigated using the lattice Boltzmann method (LBM) and random walk (RW), respectively. The numerical implementation is developed and validated on general purpose graphic processing units (GPGPUs). Both the LBM and RW method are well suited to parallel implementation on GPGPUs because they require only next-neighbour communication and thus can reduce expenses. The LBM model is an order of magnitude faster on GPGPUs than published results for LBM simulations run on modern CPUs. The fluid model is verified for parallel plate flow, backward facing step and single fracture flow; and the RW model is verified for point-source diffusion, Taylor-Aris dispersion and breakthrough behaviour in a single fracture. Both algorithms place limitations on the discrete displacement of fluid or particle transport per time step to minimise the numerical error that must be considered during implementation.

[1]  F. Durst,et al.  Lattice BGK direct numerical simulation of fully developed turbulence in incompressible plane channel flow , 2006 .

[2]  Gerhard Wellein,et al.  Performance engineering for the Lattice Boltzmann method on GPGPUs: Architectural requirements and performance results , 2011, ArXiv.

[3]  Hongbin Zhan,et al.  Experimental study of solute transport under non-Darcian flow in a single fracture , 2011 .

[4]  Jonas Tölke,et al.  Implementation of a Lattice Boltzmann kernel using the Compute Unified Device Architecture developed by nVIDIA , 2009, Comput. Vis. Sci..

[5]  Michael C. Sukop,et al.  Lattice Boltzmann Modeling: An Introduction for Geoscientists and Engineers , 2005 .

[6]  Carl E. Renshaw,et al.  On the relationship between mechanical and hydraulic apertures in rough-walled fractures , 1995 .

[7]  Andrew F. B. Tompson,et al.  Numerical simulation of solute transport in three-dimensional, randomly heterogeneous porous media , 1990 .

[8]  Brian Berkowitz,et al.  Non-Fickian Transport in Transparent Replicas of Rough-Walled Rock Fractures , 2013, Transport in Porous Media.

[9]  Haibo Huang,et al.  Evaluation of permeability and non‐Darcy flow in vuggy macroporous limestone aquifer samples with lattice Boltzmann methods , 2013 .

[10]  Richard A. Ketcham,et al.  Effects of inertia and directionality on flow and transport in a rough asymmetric fracture , 2009 .

[11]  Zhi Dou,et al.  Influence of wettability on interfacial area during immiscible liquid invasion into a 3D self-affine rough fracture: Lattice Boltzmann simulations , 2013 .

[12]  Jonas Latt,et al.  Hydrodynamic limit of lattice Boltzmann equations , 2007 .

[13]  Stacy E. Howington,et al.  Pore-scale simulation of dispersion , 2000 .

[14]  Scott C. James,et al.  Monodisperse and polydisperse colloid transport in water saturated fractures with various orientations: Gravity effects. , 2007 .

[15]  Y. Tsang,et al.  The Effect of Tortuosity on Fluid Flow Through a Single Fracture , 1984 .

[16]  Paul W. J. Glover,et al.  Fluid flow through rough fractures in rocks. II: A new matching model for rough rock fractures , 2006 .

[17]  Enhua Wu,et al.  Simulation of miscible binary mixtures based on lattice Boltzmann method , 2006, Comput. Animat. Virtual Worlds.

[18]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[19]  Stanley Reitsma,et al.  Laboratory measurement of capillary pressure‐saturation relationships in a rock fracture , 1994 .

[20]  Paul W. J. Glover,et al.  Synthetic rough fractures in rocks , 1998 .

[21]  Peter Bailey,et al.  Accelerating Lattice Boltzmann Fluid Flow Simulations Using Graphics Processors , 2009, 2009 International Conference on Parallel Processing.

[22]  Roger Beckie,et al.  The influence of surface sorption on dispersion in parallel plate fractures , 1997 .

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

[24]  Bruce Hunt,et al.  DISPERSIVE SOURCES IN UNIFORM GROUND-WATER FLOW , 1978 .

[25]  Gudmundur S. Bodvarsson,et al.  Hydraulic conductivity of rock fractures , 1996 .

[26]  H. Stockman,et al.  A lattice gas study of retardation and dispersion in fractures: Assessment of errors from desorption kinetics and buoyancy , 1997 .

[27]  Stephen R. Brown,et al.  Fluid flow through rock joints: The effect of surface roughness , 1987 .

[28]  Ahmed E. Hassan,et al.  On using particle tracking methods to simulate transport in single-continuum and dual continua porous media , 2003 .

[29]  H. Cheng,et al.  Application of Average Flow Model to Lubrication Between Rough Sliding Surfaces , 1979 .

[30]  C. R. Cole,et al.  Multicomponent mass transport model: theory and numerical implementation (discrete-parcel-random-walk version) , 1977 .

[31]  Jeffrey D Brewster.  Lattice-Boltzmann simulations of three-dimensional fluid flow on a desktop computer. , 2007, Analytical chemistry.

[32]  Harihar Rajaram,et al.  Accuracy and Computational Efficiency in 3D Dispersion via Lattice-Boltzmann: Models for Dispersion in Rough Fractures and Double-Diffusive Fingering , 1998 .

[33]  Shiyi Chen,et al.  Lattice-Boltzmann a versatile tool for multiphase Fluid Dynamics and other complicated flows , 1994 .

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

[35]  Paul W. J. Glover,et al.  Fluid flow in synthetic rough fractures and application to the Hachimantai geothermal hot dry rock test site , 1998 .

[36]  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 .

[37]  J. Boon.  The Lattice Boltzmann Equation for Fluid Dynamics and Beyond , 2003 .

[38]  Erdinc Eker,et al.  Lattice Boltzmann Simulation of Fluid Flow in Synthetic Fractures , 2006 .

[39]  Ian G. Main,et al.  Chapter 7 ▸ – Hydromechanical Behavior of Fractured Rocks , 2003 .

[40]  Scott C. James,et al.  Transport of polydisperse colloid suspensions in a single fracture , 1999 .

[41]  John Crank,et al.  The Mathematics Of Diffusion , 1956 .

[42]  Gudmundur S. Bodvarsson,et al.  Lubrication theory analysis of the permeability of rough-walled fractures , 1991 .

[43]  Giovanni Grasselli,et al.  A method to evaluate the three-dimensional roughness of fracture surfaces in brittle geomaterials. , 2009, The Review of scientific instruments.

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

[45]  Frederick Delay,et al.  Simulating Solute Transport in Porous or Fractured Formations Using Random Walk Particle Tracking: A Review , 2005 .

[46]  B. Sleep,et al.  Colloid transport in dolomite rock fractures: effects of fracture characteristics, specific discharge, and ionic strength. , 2012, Environmental science & technology.

[47]  Mutsuo Saito,et al.  Variants of Mersenne Twister Suitable for Graphic Processors , 2010, TOMS.

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

[49]  Paul W. J. Glover,et al.  Modelling fluid flow in rough fractures: Application to the Hachimantai geothermal HDR test site , 1997 .