Macroscopic model and truncation error of discrete Boltzmann method
Introduction
As a mesoscopic numerical approach, the Boltzmann method has been proven to be a useful mathematical tool to simulate lots of complicated physical phenomena including fluid flow, interface problems, particle and suspension flows for its simplicity and flexibility in numerical aspects [1], [2], [3], [4], [5], [6]. It is constructed based on kinetically microscopic view to mimic macroscopically physical process in the way that the macroscopic fluid behavior is the natural response of assembled microscopic particles [7], [8]. Meticulous particle distributions in velocity phase space are put forward to represent fluid properties from their collective effects. Basically, only streaming and collision operators for the particle distribution evolution are required in practicing a simulation. The macroscopic dynamics will be recovered by assembling microscopic kinetic equation with appropriate equilibrium particle distribution model and collision relaxation time scale. That is, computational parameters regarding particle distribution and collision relaxation time in the discrete Boltzmann method must be properly set up to recuperate the accordingly macroscopic model.
Chapman–Enskog expansion is a viable technique in relating the discrete Boltzmann method to its macroscopic counterpart. It analyzes the particle distribution evolution with two time scales stemmed from flow convection and diffusion processes. The evolution equation is subsequently expanded with a power series of a small parameter (ε) delegating the time scale ratio. Besides, particle distribution is assumed to approximately differ from its equilibrium state by a power series of this small parameter. Macroscopically equivalent equations will be ensued from agglomerating all components of particle distributions and equalizing various order of this pretentious parameter in the resultant evolution equation [4], [9], [10], [11], [12]. For example, in computational fluid dynamics (CFD), the Euler equations are obtained from the lowest order and the Navier–Stokes equations by putting in additional diffusion terms deduced from the second lowest order term. Furthermore, as suggested by Latt [13], with appropriate correction terms in collision process, it will substantially increase the stability and accuracy in numerical simulations which is revealed in the resulting macroscopic model. Although its effectiveness to attain the macroscopic model can be assured, concerns regarding its meaning and interpretation consequently arise [14], [15], [16], [17]:
- 1.
Execution of the Chapman–Enskog expansion is essentially based on presumption of two time scales exposed in the continuous model. However, this is not a prerequisite requirement to implement the discrete Boltzmann method. That is, from the numerical viewpoint, the discrete Boltzmann method may be regarded as a mere mathematical technique to handle particle streaming and collision operators.
- 2.
The intricate meaning of the introduced small parameter (ε) accounting for the time scale ratio between different physical processes requires rigorous phenomenal explanation. In general, it is interpreted as the Knudsen number to measure the ratio between particle mean-free path and system characteristic size. However, Watari [17] pointed out that it is natural to use time step (Δt) as ε which consequentially leads to an alternative in deriving macroscopic model. Without a stringent evidence on the interpretation of ε, it will become conceptually perplexed in the applications of the Boltzmann method.
- 3.
In a convection–diffusion problem, two time scales are used to depict the characteristic times associated with the convection and diffusion processes, respectively. This small parameter (ε) is also taken up to signify the ratio between these time scales. Nevertheless, it becomes ambiguous to explain its physical meaning if there exists an additional physical process with different characteristic time scale.
- 4.
Although macroscopic model can be acquired with the Chapman–Enskog expansion, it is difficult to categorize the associated truncation error adherent to the discrete Boltzmann method. Higher-order terms such as Burnett and super-Burnett hydrodynamics may violate the basic physics behind the Boltzmann equation [15], [18].
- 5.
Deviation of particle distribution function from its equilibrium state is expressed as a power series of the same small parameter. The first (ε) and second () terms are regarded as those brought about by the convection and diffusion processes, respectively. It is not clearly explained the physical meanings of the higher-order terms in this approximate series.
The objective of present study is to propose an alternative technique to derive the macroscopic model of the discrete Boltzmann method without the presumptions in the Chapman–Enskog expansion. Thanks to the linearity introduced by the BGK approximation [23], original formulation describing particle distribution function evolution can be reorganized to its equivalent form for equilibrium state. That is, analyses will be conducted with evolution of equilibrium distribution which is directly related to macroscopic properties. Taylor series expansion incorporating with fundamental algebraic manipulations on the discrete Boltzmann method can then be employed to deduce its macroscopic counterpart. Without the intrusion of artificially multiple time scales and confounding small parameter on the derivation procedure, a more comprehension on the conceptual insight of the discrete Boltzmann method can then be obtained.
Besides this introductory section, the present paper is organized as follows. Analysis is conducted for the one-dimensional convection–diffusion problem in Section 2. Both Chapman–Enskog and present proposition are investigated to depict their essential ingredients for comparison. The present analysis is verified by a numerical experiment performed in Section 3. It emphasizes on the truncation error predicted by present formulation which is difficult to attain in the Chapman–Enskog expansion. In the last section, conclusive remarks drawn from the present study are exhibited. Extensions to deal with convection–diffusion–reaction system, multi-relaxation times in collision process as well as multi-dimensional Navier–Stokes equations are exposed in Appendix.
Section snippets
One-dimensional convection–diffusion equation
To demonstrate the essential ingredients of present formulation, detailed derivation procedure for the simple one-dimensional linear convection–diffusion problem will be described in this section. Both the Chapman–Enskog expansion and present proposition are displayed for comparison.
One-dimensional case
The linear one-dimensional convection–diffusion problem depicted in Eq. (1) with the following initial and boundary conditions in the range of is solved to validate the present analyses: The corresponding exact solution can be easily acquired: with First, this problem is solved by a standard lattice Boltzmann method (LBM) given in Eq. (2) with D1Q2 equilibrium particle
Conclusions
In this paper, a novel technique to derive the macroscopically continuous model for the discrete Boltzmann method is proposed. It is motivated by the fact that equilibrium particle distribution is more suitable than its non-equilibrium form as a key variable for analysis since it is more directly connected to the macroscopic properties. The original evolution of non-equilibrium distribution can be converted to its equivalent equilibrium form owing to the linearity introduced by the BGK
References (26)
- et al.
The lattice Boltzmann equation method: theoretical interpretation, numerics and implications
Int. J. Multiph. Flow
(2003) - et al.
Viscosity of finite difference lattice Boltzmann models
J. Comput. Phys.
(2003) Chapman–Enskog analysis for finite-volume formulation of lattice Boltzmann equation
Physica A
(2013)- et al.
Truncation error analysis of lattice Boltzmann methods
J. Comput. Phys.
(2004) - et al.
Asymptotic analysis of the lattice Boltzmann equation
J. Comput. Phys.
(2005) - et al.
Comparison of analysis techniques for the lattice Boltzmann method
Comput. Math. Appl.
(2009) Higher-order exponential difference schemes for the computations of the steady convection–diffusion equation
J. Comput. Phys.
(1996)- et al.
An evaluation of upwind and central difference approximations by a study of recirculating flow
Comput. Fluids
(1976) - et al.
The lattice Boltzmann equation: theory and applications
Phys. Rep.
(1992) - et al.
Lattice Boltzmann simulation of solid particle suspended in fluid
J. Stat. Phys.
(1995)
Lattice-Gas Cellular Automata: Simple Model of Complex Hydrodynamics
Lattice Boltzmann method for fluid flow
Annu. Rev. Fluid Mech.
Fully Lagrangian and lattice Boltzmann methods for the advection–diffusion equation
J. Sci. Comput.
Cited by (4)
Derivations of Exact Lattice Boltzmann Evolution Equation
2020, Journal of Shanghai Jiaotong University (Science)Derivation of lattice Boltzmann equation via analytical characteristic integral
2019, Chinese Physics B