A lattice Boltzmann equation for diffusion

The formulation of lattice gas automata (LGA) for given partial differential equations is not straightforward and still requires some sort of magic. Lattice Boltzmann equation (LBE) models are much more flexible than LGA because of the freedom in choosing equilibrium distributions with free parameters which can be set after a multiscale expansion according to certain requirements. Here a LBE is presented for diffusion in an arbitrary number of dimensions. The model is probably the simplest LBE which can be formulated. It is shown that the resulting algorithm with relaxation parameter =1 is identical to an explicit finite-difference (EFD) formulation at its stability limit. Underrelaxation (0<<1) allows stable integration beyond the stability limit of EFD. The time step of the explicit LBE integration is limited by accuracy and not by stability requirements.     

Asymptotic analysis of the lattice Boltzmann equation

In this article we analyze the lattice Boltzmann equation (LBE) by using the asymptotic expansion technique. We first relate the LBE to the finite discrete-velocity model (FDVM) of the Boltzmann equation with the diffusive scaling. The analysis of this model directly leads to the incompressible Navier–Stokes equations, as opposed to the compressible Navier–Stokes equations obtained by the Chapman–Enskog analysis with convective scaling. We also apply the asymptotic analysis directly to the fully discrete LBE, as opposed to the usual practice of analyzing a continuous equation obtained through the Taylor-expansion of the LBE. This leads to a consistency analysis which provides order-by-order information about the numerical solution of the LBE. The asymptotic technique enables us to analyze the structure of the leading order errors and the accuracy of numerically derived quantities, such as vorticity. It also justifies the use of Richardson’s extrapolation method. As an example, a two-dimensional Taylor-vortex flow is used to validate our analysis. The numerical results agree very well with our analytic predictions.     

Generalized Boltzmann equation for lattice gas automata

In this paper a theory is formulated that predicts velocity and spatial correlations between occupation numbers that occur in lattice gas automata violating semi-detailed balance. Starting from a coupled BBGKY hierarchy for then-particle distribution functions, cluster expansion techniques are used to derive approximate kinetic equations. In zeroth approximation the standard nonlnear Boltzmann equation is obtained; the next approximation yields the ring kinetic equation, similar to that for hard-sphere systems, describing the time evolution of pair correlations. The ring equation is solved to determine the (nonvanishing) pair correlation functions in equilibrium for two models that violate semidetailed balance. One is a model of interacting random walkers on a line, the other one is a two-dimensional fluid-type model on a triangular lattice. The numerical predictions agree very well with computer simulations.     

Global solutions of the Boltzmann equation on a lattice

The nonlinear Boltzmann equation with a discretized spatial variable is studied in a Banach space of absolutely integrable functions of the velocity variables. Conservation laws and positivity are utilized to extend weak local solutions to a global solution. This is shown to be a strong solution by analytic semigroup techniques.Supported by National Science Foundation Grant ENG-7515882.     

Thermal equation of state for lattice Boltzmann gases

The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar--Gross--Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar--Gross--Krook equation (LBGKE) induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics.     

The lattice Boltzmann phononic lattice solid

I present a Boltzmann lattice gas-like approach for modeling compressional waves in an inhomogeneous medium as a first step toward developing a method to simulate seismic waves in complex solids. The method is based on modeling particles in a discrete lattice with wavelike characteristics of partial reflection and transmission when passing between links with different properties as well as phononlike interactions (i.e., collisions), with particle speed dependent on link properties. In the macroscopic limit, this approach theoretically yields compressional waves in an inhomogeneous acoustic medium. Numerical experiments verify the method and demonstrate its convergence properties. The lattice Boltzmann phononic lattice solid could be used to study how seismic wave anisotropy and attenuation are related to microfractures, the complex geometry of rock matrices, and their couplings to pore fluids. However, additional particles related to the two transverse phonons must be incorporated to correctly simulate wave phenomena in solids.     

A lattice Boltzmann model for the generalized Burgers-Huxley equation

In this paper we develop a lattice Boltzmann model for the generalized Burgers-Huxley equation (GBHE). By choosing the proper time and space scales and applying the Chapman-Enskog expansion, the governing equation is recovered correctly from the lattice Boltzmann equation, and the local equilibrium distribution functions are obtained. Excellent agreement with the exact solution is observed, and better numerical accuracy is obtained than the available numerical result. The results indicate the present model is satisfactory and efficient. The method can also be applied to the generalized Burgers-Fisher equation and be extended to multidimensional cases.     

Derivation of lattice Boltzmann equation via analytical characteristic integral

A lattice Boltzmann(LB) theory, the analytical characteristic integral(ACI) LB theory, is proposed in this paper.ACI LB theory takes the Bhatnagar–Gross–Krook(BGK)-Boltzmann equation as the exact kinetic equation behind Navier–Stokes continuum and momentum equations and constructs an LB equation by rigorously integrating the BGK-Boltzmann equation along characteristics. It is a general theory, supporting most existing LB equations including the standard lattice BGK(LBGK) equation inherited from lattice-gas automata, whose theoretical foundation had been questioned. ACI LB theory also indicates that the characteristic parameter of an LB equation is collision number, depicting the particle-interaction intensity in the time span of the LB equation, instead of the traditionally assumed relaxation time, and the over-relaxation time problem is merely a manifestation of the temporal evolution of equilibrium distribution along characteristics under high collision number, irrelevant to particle kinetics. In ACI LB theory, the temporal evolution of equilibrium distribution along characteristics is the determinant of LB method accuracy and we numerically prove this.     

Boltzmann equation on a lattice: Existence and uniqueness of solutions

Cercignani, Greenberg, and Zweifel proved the existence and uniqueness of solutions of the Boltzmann equation on a toroidal lattice under the assumption that the collision kernel is bounded. We give an alternative, considerably simpler, proof which is based on a fixed point argument.     

The Boltzmann equation

An abstract form of the spatially non-homogeneous Boltzmann equation is derived which includes the usual, more concrete form for any kind of potential, hard or soft, with finite cutoff. It is assumed that the corresponding gas is confined to a bounded domain by some sort of reflection law. The problem then considered is the corresponding initial-boundary value problem, locally in time.Two proofs of existence are given. Both are constructive, and the first, at least, provides two sequences, one converging to the solution from above, the other from below, thus producing, at the same time as existence, approximations to the solution and error bounds for the approximation.The solution is found within a space of functions bounded by a multiple of a Maxwellian, and, in this space, uniqueness is also proved.Research supported, in part, by the National Research Council of Canada (NRC A8560)     

A consistent lattice Boltzmann equation with baroclinic coupling for mixtures

We propose a consistent lattice Boltzmann equation (LBE) with baroclinic coupling between species and mixture dynamics to model the active scalar dynamics in multi-species mixtures. The proposed LBE model is directly derived from the linearized Boltzmann equations for mixtures and it has the following two distinctive features. First, it uses the multiple-relaxation-time collision model so that it has the flexibility of independent Reynolds and Schmidt numbers, and better numerical stability. Second, it satisfies the indifferentiability principle therefore leads to a set of consistent hydrodynamic equations for barycentric velocity for mixtures. The proposed LBE model is validated through simulations of decaying homogeneous isotropic turbulence in three dimensions. We simulate both the active and passive scalar dynamics in decaying turbulence for mixtures. We also compute various statistical quantities and their decay exponents in decaying turbulence. Our results agree well with existing results for both scalar dynamics and decaying turbulence.     

Numerical method based on the lattice Boltzmann model for the Fisher equation

In this paper, a lattice Boltzmann model for the Fisher equation is proposed. First, the Chapman-Enskog expansion and the multiscale time expansion are used to describe higher-order moment of equilibrium distribution functions and a series of partial differential equations in different time scales. Second, the modified partial differential equation of the Fisher equation with the higher-order truncation error is obtained. Third, comparison between numerical results of the lattice Boltzmann models and exact solution is given. The numerical results agree well with the classical ones.     

A fluctuating lattice Boltzmann scheme for the one-dimensional KPZ equation

Based on the well-known mapping between the Burgers equation with noise and the Kardar–Parisi–Zhang (KPZ) equation for fluctuating interfaces, we develop a fluctuating lattice Boltzmann (LB) scheme for growth phenomena, as described by the KPZ formalism. A very simple LB-KPZ scheme is demonstrated in 1+1 spacetime dimensions, and is shown to reproduce the scaling exponents characterizing the growth of one-dimensional fluctuating interfaces.     

The generalized Boltzmann equation

From the BBGKY-hierarchy of the kinetic equations taking into account three possible levels of scales, connected with the mean time of particle interactions, mean time between collisions and the hydrodynamic time, the generalized Boltzmann equation is derived. The generalized H-theorem is proven.     

Multiple-relaxation-time lattice Boltzmann model for the convection and anisotropic diffusion equation

A lattice Boltzmann model with a multiple-relaxation-time (MRT) collision operator for the convection–diffusion equation is presented. The model uses seven discrete velocities in three dimensions (D3Q7 model). The off-diagonal components of the relaxation-time matrix, which originate from the rotation of the principal axes, enable us to take into account full anisotropy of diffusion. An asymptotic analysis of the model equation with boundary rules for the Dirichlet and Neumann-type (specified flux) conditions is carried out to show that the model is first- and second-order accurate in time and space, respectively. The results of the analysis are verified by several numerical examples. It is also shown numerically that the error of the MRT model is less sensitive to the variation of the relaxation-time coefficients than that of the classical BGK model. In addition, an alternative treatment for the Neumann-type boundary condition that improves the accuracy on a curved boundary is presented along with a numerical example of a spherical boundary.     

The lattice Boltzmann advection-diffusion model revisited

Advection-diffusion processes can be simulated by the Lattice Boltzmann method. Two formulations have been proposed in the literature. We show that they are not fully correct (only first order accurate). A new formulation is proposed, which is shown to produce better results, both from the point of view of the Chapman-Enskog expansion or when comparing simulations with an exact time-dependent solution of the advection-diffusion equation.     

The p-q model Boltzmann equation

The “p-q” model earlier introduced by the authors to describe persistent scattering under a scalar Boltzmann equation is here examined in detail. After deriving the scattering kernel and exhibiting its properties we obtain moment and similarity solutions and show how the model effectively parametrizes all intermediate conditions between the extremes of diffusion-like “small-scattering” and the strong-collisional limit of “diffuse-scattering” characteristic of earlier, more restrictive models. Both continuous and discrete-variable versions of the model are discussed and shown to be straightforwardly interrelated. Our derivations, carried out in natural energy-like variables, parallel those given recently by Ernst and Hendriks using transform methods.     

应用格子Boltzmann模型模拟一类二维偏微分方程

针对一类含源的二维非线性偏微分方程, 通过Chapman-Enskog展开技术和多尺度分析提出了带修正项的简单格子Boltzmann模型. 用模型模拟了几类二维偏微分方程, 数值模拟结果与精确解相符合. 成功将格子Boltzmann方法应用到二维偏微分方程的数值求解中. 关键词： 二维非线性偏微分方程 格子Boltzmann模型 Chapman-Enskog多尺度展开