Macroscale lattice‐Boltzmann methods for low Peclet number solute and heat transport in heterogeneous porous media

[1] This paper introduces new methods for simulating subsurface solute and heat transport in heterogeneous media using large-scale lattice-Boltzmann models capable of representing both macroscopically averaged porous media and open channel flows. Previous examples of macroscopically averaged lattice-Boltzmann models for solute and heat transport are only applicable to homogeneous media. Here, we extend these models to properly account for heterogeneous pore-space distributions. For simplicity, in the majority of this paper we assume low Peclet number flows with an isotropic dispersion tensor. Nevertheless, this approach may also be extended to include anisotropic-dispersion by using multiple relaxation time lattice-Boltzmann methods. We describe two methods for introducing heterogeneity into macroscopically averaged lattice-Boltzmann models. The first model delivers the desired behavior by introducing an additional time-derivative term to the collision rule; the second model by separately weighting symmetric and anti-symmetric components of the fluid packet densities. Chapman-Enskog expansions are conducted on the governing equations of the two models, demonstrating that the correct constitutive behavior is obtained in both cases. In addition, methods for improving model stability at low porosities are also discussed: (1) an implicit formulation of the model; and (2) a local transformation that normalizes the lattice-Boltzmann model by the local porosity. The model performances are evaluated through comparisons of simulated results with analytical solutions for one- and two-dimensional flows, and by comparing model predictions to finite element simulations of advection isotropic-dispersion in heterogeneous porous media. We conclude by presenting an example application, demonstrating the ability of the new models to couple with simulations of reactive flow and changing flow geometry: a simulation of groundwater flow through a carbonate system.

[1]  Dieter Wolf-Gladrow,et al.  3. Lattice-gas cellular automata , 2000 .

[2]  A. Lachenbruch,et al.  Thermal conductivity of rocks from measurements on fragments and its application to heat‐flow determinations , 1971 .

[3]  I. Ginzburg.  Equilibrium-type and link-type lattice Boltzmann models for generic advection and anisotropic-dispersion equation , 2005 .

[4]  David Freed,et al.  Lattice-Boltzmann Method for Macroscopic Porous Media Modeling , 1998 .

[5]  Vedat Batu.  Two‐Dimensional Dispersion from Strip Sources , 1983 .

[6]  Vedat Batu,et al.  Applied Flow and Solute Transport Modeling in Aquifers: Fundamental Principles and Analytical and Numerical Methods , 2005 .

[7]  Robert S. Bernard,et al.  Simulation of flow through bead packs using the lattice Boltzmann method , 1998 .

[8]  F. Gabrovšek,et al.  Processes of Speleogenessis , 2005 .

[9]  M. Sukop,et al.  Lattice Boltzmann model for the elder problem , 2004 .

[10]  Marcel G. Schaap,et al.  Comparison of pressure‐saturation characteristics derived from computed tomography and lattice Boltzmann simulations , 2007 .

[11]  J. Bear.  Hydraulics of Groundwater , 1979 .

[12]  J. Bahr.  Kinetically influenced terms for solute transport affected by heterogeneous and homogeneous classical reactions , 1990 .

[13]  Karsten Pruess,et al.  Mineral Sequestration of Carbon Dixoide in a Sandstone-Shale System , 2004 .

[14]  M. Saar,et al.  Depth dependence of permeability in the Oregon cascades inferred from hydrogeologic, thermal, seismic, and magmatic modeling constraints , 2004 .

[15]  H. Rajaram,et al.  Dissolution of limestone fractures by cooling waters: Early development of hypogene karst systems , 2005 .

[16]  C. Pan,et al.  Lattice‐Boltzmann simulation of two‐phase flow in porous media , 2004 .

[17]  John McCloskey,et al.  Lattice Boltzmann scheme with real numbered solid density for the simulation of flow in porous media , 1998 .

[18]  C. Steefel,et al.  A new kinetic approach to modeling water-rock interaction: The role of nucleation, precursors, and Ostwald ripening , 1990 .

[19]  Qinjun Kang,et al.  Lattice Boltzmann simulation of chemical dissolution in porous media. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  J. Stefan,et al.  Ueber die Theorie der Eisbildung, insbesondere über die Eisbildung im Polarmeere , 1891 .

[21]  Stuart D. C. Walsh,et al.  A new partial-bounceback lattice-Boltzmann method for fluid flow through heterogeneous media , 2009, Comput. Geosci..

[22]  C. Neuzil,et al.  Groundwater in geologic processes, 2nd edition , 2006 .

[23]  F Verhaeghe,et al.  Lattice Boltzmann model for diffusion-controlled dissolution of solid structures in multicomponent liquids. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[24]  Y. Qian,et al.  Lattice BGK Models for Navier-Stokes Equation , 1992 .

[25]  M. Banda,et al.  Lattice Boltzmann simulation of dispersion in two‐dimensional tidal flows , 2009 .

[26]  Christoph Clauser,et al.  Numerical simulation of reactive flow in hot aquifers : SHEMAT and processing SHEMAT , 2003 .

[27]  Daniel H. Rothman,et al.  Lattice-Boltzmann simulations of flow through Fontainebleau sandstone , 1995 .

[28]  S. Chapman.  On the Kinetic Theory of a Gas. Part II: A Composite Monatomic Gas: Diffusion, Viscosity, and Thermal Conduction , 1918 .

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

[30]  A. Howard,et al.  Early development of karst systems: 1. Preferential flow path enlargement under laminar flow , 1994 .

[31]  D. Wolf-Gladrow.  Lattice-Gas Cellular Automata and Lattice Boltzmann Models: An Introduction , 2000 .

[32]  A. Ladd,et al.  Simulation of chemical erosion in rough fractures. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[33]  H. S. Fogler,et al.  The kinetics of calcite dissolution in acetic acid solutions , 1998 .

[34]  Sauro Succi,et al.  Exponential Tails in Two-Dimensional Rayleigh-Bénard Convection , 1993 .

[35]  John W. Crawford,et al.  A lattice BGK model for advection and anisotropic dispersion equation , 2002 .

[36]  Elaine S. Oran,et al.  Numerical Simulation of Reactive Flow , 1987 .

[37]  Dongxiao Zhang,et al.  Unified lattice Boltzmann method for flow in multiscale porous media. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  Taehun Lee,et al.  An Eulerian description of the streaming process in the lattice Boltzmann equation , 2003 .

[39]  S. Chapman,et al.  On the Law of Distribution of Molecular Velocities, and on the Theory of Viscosity and Thermal Conduction, in a Non-Uniform Simple Monatomic Gas , 1916 .

[40]  John McCloskey,et al.  Permeability porosity relationships from numerical simulations of fluid flow , 1998 .

[41]  S. Succi,et al.  Three-Dimensional Flows in Complex Geometries with the Lattice Boltzmann Method , 1989 .

[42]  E. LeBoeuf,et al.  Least-squares finite-element scheme for the lattice Boltzmann method on an unstructured mesh. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[43]  M. Person.  Groundwater in Geologic Processes , 2000 .

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

[45]  Remo Guidieri.  Res , 1995, RES: Anthropology and Aesthetics.

[46]  Shiyi Chen,et al.  Lattice Boltzmann computations for reaction‐diffusion equations , 1993 .

[47]  A. Palmer.  Origin and morphology of limestone caves , 1991 .

[48]  Shi Jin,et al.  Physical symmetry and lattice symmetry in the lattice Boltzmann method , 1997 .

[49]  G. L. Guymon.  Hydraulics of groundwater: Jacob Bear McGraw-Hill, New York, £18.55 , 1980 .

[50]  Peter Bailey,et al.  Accelerating geoscience and engineering system simulations on graphics hardware , 2009, Comput. Geosci..

[51]  Daniel H. Rothman,et al.  Lattice-Gas Cellular Automata: Simple Models of Complex Hydrodynamics , 1997 .

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

[53]  T. K. Perkins,et al.  A Review of Diffusion and Dispersion in Porous Media , 1963 .

[54]  Vedat Batu,et al.  A generalized two-dimensional analytical solution for hydrodynamic dispersion in bounded media with the first-type boundary condition at the source , 1989 .

[55]  M. Saar,et al.  A dimensionless number describing the effects of recharge and geometry on discharge from simple karstic aquifers , 2009 .

[56]  G. Marsily.  Quantitative Hydrogeology: Groundwater Hydrology for Engineers , 1986 .

[57]  Joanne T. Fredrich,et al.  Predicting macroscopic transport properties using microscopic image data , 2005 .

[58]  L. Boltzmann.  Weitere Studien über das Wärmegleichgewicht unter Gasmolekülen , 1970 .

[59]  Zhaoli Guo,et al.  Explicit finite-difference lattice Boltzmann method for curvilinear coordinates. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[60]  A. Howard,et al.  Early Development of Karst Systems: 2. Turbulent Flow , 1995 .