A Boltzmann-based mesoscopic model for contaminant transport in flow systems

Abstract The objective of this paper is to demonstrate the formulation of a numerical model for mass transport based on the Bhatnagar–Gross–Krook (BGK) Boltzmann equation. To this end, the classical chemical transport equation is derived as the zeroth moment of the BGK Boltzmann differential equation. The relationship between the mass transport equation and the BGK Boltzmann equation allows an alternative approach to numerical modeling of mass transport, wherein mass fluxes are formulated indirectly from the zeroth moment of a difference model for the BGK Boltzmann equation rather than directly from the transport equation. In particular, a second-order numerical solution for the transport equation based on the discrete BGK Boltzmann equation is developed. The numerical discretization of the first-order BGK Boltzmann differential equation is straightforward and leads to diffusion effects being accounted for algebraically rather than through a second-order Fickian term. The resultant model satisfies the entropy condition, thus preventing the emergence of non-physically realizable solutions including oscillations in the vicinity of the front. Integration of the BGK Boltzmann difference equation into the particle velocity space provides the mass fluxes from the control volume and thus the difference equation for mass concentration. The difference model is a local approximation and thus may be easily included in a parallel model or in accounting for complex geometry. Numerical tests for a range of advection–diffusion transport problems, including one- and two-dimensional pure advection transport and advection–diffusion transport show the accuracy of the proposed model in comparison to analytical solutions and solutions obtained by other schemes.

[1]  C. Chu.  Kinetic‐Theoretic Description of the Formation of a Shock Wave , 1965 .

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

[3]  Mohamed Salah Ghidaoui,et al.  Computation of Open Channel Flow Using a Gas-Kinetic Based Scheme , 1998 .

[4]  Rolf D. Reitz,et al.  One-dimensional compressible gas dynamics calculations using the Boltzmann equation , 1981 .

[5]  Daniel H. Rothman,et al.  Cellular‐automaton fluids: A model for flow in porous media , 1988 .

[6]  S. Zaleski,et al.  Lattice Boltzmann model of immiscible fluids. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[7]  Yue-Hong Qian,et al.  Simulating thermohydrodynamics with lattice BGK models , 1993 .

[8]  William G. Gray,et al.  A Boltzmann based model for open channel flows , 2001 .

[9]  Y. Pomeau,et al.  Lattice-gas automata for the Navier-Stokes equation. , 1986, Physical review letters.

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

[11]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[12]  Mohamed Salah Ghidaoui,et al.  Low-Speed Flow Simulation by the Gas-Kinetic Scheme , 1999 .

[13]  Kun Xu.  A Gas-Kinetic Scheme for the Euler Equations with Heat Transfer , 1999, SIAM J. Sci. Comput..

[14]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[15]  Kun Xu,et al.  BGK-Based Scheme for Multicomponent Flow Calculations , 1997 .

[16]  Kun Xu,et al.  Analysis and implementation of the gas-kinetic BGK scheme for computational gas dynamics , 1997 .

[17]  Antony Jameson,et al.  Gas-kinetic finite volume methods, flux-vector splitting, and artificial diffusion , 1995 .

[18]  B. Leer,et al.  Flux-vector splitting for the Euler equations , 1997 .

[19]  D. Martínez,et al.  Lattice Boltzmann magnetohydrodynamics , 1994, comp-gas/9401002.

[20]  Philip L. Roe,et al.  The use of the Riemann problem in finite difference schemes , 1989 .

[21]  G. Lebon,et al.  Extended irreversible thermodynamics , 1993 .

[22]  J. Steger,et al.  Flux vector splitting of the inviscid gasdynamic equations with application to finite-difference methods , 1981 .

[23]  T. Teichmann,et al.  Introduction to physical gas dynamics , 1965 .

[24]  Antony Jameson,et al.  BGK-Based Schemes for the Simulation of Compressible Flow , 1996 .

[25]  J. Quirk.  A Contribution to the Great Riemann Solver Debate , 1994 .

[26]  B. P. Leonard,et al.  The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection , 1991 .

[27]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[28]  C. Hirsch.  Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows , 1990 .

[29]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[30]  Takaji Inamuro,et al.  Lattice Boltzmann simulation of flows in a three‐dimensional porous structure , 1999 .

[31]  H. Ted Davis,et al.  Application of the lattice‐Boltzmann method to study flow and dispersion in channels with and without expansion and contraction geometry , 1999 .

[32]  S. Osher,et al.  Stable and entropy satisfying approximations for transonic flow calculations , 1980 .

[33]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

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