Determination of soil hydraulic conductivity with the lattice Boltzmann method and soil thin-section technique

A pore-scale modelling of soil hydraulic conductivity using the lattice Boltzmann equation method and thin-section technique is presented in this paper. Two-dimensional thin sections taken from soil samples were used to measure soil pore geometry, and three-dimensional soil structures were then reconstructed based on the thin-section measurements and an assumption that the soil structure can be fully characterised by its porosity and variogram. The hydraulic conductivity of the reconstructed three-dimensional soil structures was calculated from a lattice Boltzmann simulation of water flow in the structures, in which the water-solid interface was treated as a non-slip boundary (zero-velocity boundary) and solved by the bounce-back method. To improve the accuracy of the bounce-back method, the particle distribution functions were located at the centre of the cube that was fully occupied either by water or by solid. The simulated hydraulic conductivity was compared with measured hydraulic conductivity and the results showed good agreement. We also analysed the probability distribution of the simulated water velocity and the results indicated that the transverse velocity components had a non-Gaussian symmetric distribution, while the longitudinal velocity component had a skewed-forward distribution. Both distributions had a marked peak for velocity close to zero, indicating a significant portion of stagnant water.

[1]  E. Fitzpatrick,et al.  Micromorphology of Soils , 1984 .

[2]  B. R. Sehgal,et al.  Evaluation of the Darcy's law performance for two-fluid flow hydrodynamics in a particle debris bed using a lattice-Boltzmann model , 2000 .

[3]  J. Hopmans,et al.  Pore‐scale measurements of solute breakthrough using microfocus X‐ray computed tomography , 2000 .

[4]  Monica Moroni,et al.  Three‐dimensional particle tracking velocimetry studies of the transition from pore dispersion to Fickian Dispersion for homogeneous porous media , 2001 .

[5]  Qisu Zou,et al.  A improved incompressible lattice Boltzmann model for time-independent flows , 1995 .

[6]  Yoshito Nakashima,et al.  Estimate of transport properties of porous media by microfocus X‐ray computed tomography and random walk simulation , 2002 .

[7]  Takaji Inamuro,et al.  A NON-SLIP BOUNDARY CONDITION FOR LATTICE BOLTZMANN SIMULATIONS , 1995, comp-gas/9508002.

[8]  A. Hoekstra,et al.  Permeability of Three-Dimensional Random Fiber Webs , 1998 .

[9]  R. Maier,et al.  Pore-Scale Flow and Dispersion , 1998 .

[10]  D. Martínez,et al.  On boundary conditions in lattice Boltzmann methods , 1996 .

[11]  M. Císlerová,et al.  CT derived porosity distribution and flow domains , 2002 .

[12]  Koo-Chul Lee,et al.  Nonlinear correction to Darcy's law for a flow through periodic arrays of elliptic cylinders , 2001 .

[13]  Yakov A. Pachepsky,et al.  Fractal parameters of pore surfaces as derived from micromorphological data: effect of long-term management practices , 1996 .

[14]  Chen,et al.  Simulation of multicomponent fluids in complex three-dimensional geometries by the lattice Boltzmann method. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[15]  Michael A. Celia,et al.  Recent advances in pore scale models for multiphase flow in porous media , 1995 .

[16]  C. Pan,et al.  Pore-scale modeling of saturated permeabilities in random sphere packings. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  David R. Noble,et al.  Direct assessment of lattice Boltzmann hydrodynamics and boundary conditions for recirculating flows , 1995 .

[18]  Shiyi Chen,et al.  LATTICE BOLTZMANN METHOD FOR FLUID FLOWS , 2001 .

[19]  N. Bird,et al.  Simulation of soil water retention using random fractal networks , 1997 .

[20]  E. Moreau,et al.  Pore networks in an Italian Vertisol: quantitative characterisation by two dimensional image analysis , 1996 .

[21]  Mark L. Rivers,et al.  Using X-ray computed tomography in hydrology: systems, resolutions, and limitations , 2002 .

[22]  J. Crawford,et al.  Effect of sampling volume on the measurement of soil physical properties: simulation with x-ray tomography data , 2002 .

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

[24]  J. Buick,et al.  Gravity in a lattice Boltzmann model , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[25]  W. E. Soll,et al.  Pore scale study of flow in porous media: Scale dependency, REV, and statistical REV , 2000 .

[26]  Albert J. Valocchi,et al.  Pore‐scale modeling of dissolution from variably distributed nonaqueous phase liquid blobs , 2001 .

[27]  L. E. Scriven,et al.  NMR imaging of velocity profiles and velocity distributions in bead packs , 1996 .

[28]  Cathy J. Ridgway,et al.  Void space modeling of mercury intrusion hysteresis in sandstone, paper coating, and other porous media , 1995 .

[29]  Matthaeus,et al.  Recovery of the Navier-Stokes equations using a lattice-gas Boltzmann method. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[30]  John W. Crawford,et al.  On boundary conditions in the lattice Boltzmann model for advection and anisotropic dispersion equation , 2002 .

[31]  Ladd,et al.  Lattice-boltzmann model with sub-grid-scale boundary conditions , 2000, Physical review letters.

[32]  L. Luo,et al.  Analytic solutions of simple flows and analysis of nonslip boundary conditions for the lattice Boltzmann BGK model , 1997 .

[33]  J. Crawford.  The relationship between structure and the hydraulic conductivity of soil , 1994 .

[34]  U. Pietrzyk,et al.  Flow visualization in porous media via Positron Emission Tomography , 1998 .

[35]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[36]  Alex B. McBratney,et al.  Soil Mass, Surface, and Spectral Fractal Dimensions Estimated from Thin Section Photographs , 1996 .

[37]  Dominique d'Humières,et al.  Multireflection boundary conditions for lattice Boltzmann models. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  Wei Shyy,et al.  Regular Article: An Accurate Curved Boundary Treatment in the Lattice Boltzmann Method , 1999 .

[39]  Hans-Jörg Vogel,et al.  A new approach for determining effective soil hydraulic functions , 1998 .

[40]  P. Lallemand,et al.  Lattice Boltzmann method for moving boundaries , 2003 .

[41]  Peter M. A. Sloot,et al.  Implementation Aspects of 3D Lattice-BGK: Boundaries, Accuracy, and a New Fast Relaxation Method , 1999 .

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

[43]  I M Young,et al.  Interactions and Self-Organization in the Soil-Microbe Complex , 2004, Science.

[44]  D. d'Humières,et al.  Local second-order boundary methods for lattice Boltzmann models , 1996 .

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

[46]  P. Viotti,et al.  Lagrangian Analysis of Nonreactive Pollutant Dispersion in Porous Media by Means of the Particle Image Velocimetry Technique , 1996 .

[47]  Muhammad Sahimi,et al.  Pore network modelling of two-phase flow in porous rock: the effect of correlated heterogeneity , 2001 .

[48]  P. Worsfold,et al.  Simulation of water retention and hydraulic conductivity in soil using a three‐dimensional network , 2000 .

[49]  Hans-Jörg Vogel,et al.  Quantitative morphology and network representation of soil pore structure , 2001 .

[50]  John W. Crawford,et al.  An Efficient Markov Chain Model for the Simulation of Heterogeneous Soil Structure , 2004 .