Theoretical investigations on lattice Boltzmann method: an amended MBD and improved LBM

Abstract

This paper presents theoretical investigations of lattice Boltzmann method (LBM) to develop a completed LBM theory. Based on H-theorem with Lagrangian multiplier method, an amended theoretical equilibrium distribution function (EDF) is derived, which modifies the current Maxwell-Boltzmann distribution (MBD) to include the total internal energy as its parameter. This modification allows the three conservation laws derived directly from lattice Boltzmann equation (LBE) without additional small-parameter expansions adopted in references. From this amended theoretical EDF, an improved LBM is developed, in which the total internal energy like the mass density and mean velocity is a new macroscopic variable to be updated for different times and cells during simulations. The developed method provides a means to consider external forces and energy generation sources as generalised forces in LBM simulations. The corresponding model and implementation process of the improved LBM are presented with its performance theoretically investigated. Analytically hand-workable examples are given to illustrate its applications and to confirm its validity. The paper will excite more researchers and scientists of this area to numerically practice the new theory and method dealing with complex physical problems, from which it is expected to further advance LBM benefiting science and engineering.

Application of shifted lattice model to 3D compressible lattice Boltzmann method

“…Since the first paper about lattice Boltzmann method (LBM) appeared in 1988, the LBM [1][2][3][4][5] as a prominent tool in computational fluid dynamics has attracted much attention and become a major research hotspot. Unlike the traditional method based on macroscopic equations, the LBM is based on a mesoscopic kinetic model and obtains macroscopic quantities through statistics of distribution functions.…”

“…f eq eq[3](u, v, w) = f eq eq[1](w, v, u), f eq eq[4](u, v, w) = f eq eq[1](−u, v, w), f eq eq[5](u, v, w) = f eq eq[1](−v, u, w), f eq eq[6](u, v, w) = f eq eq[1](−w, v, u), f eq eq[8](u, v, w) = f eq eq[7](−u, v, w), f eq eq[9](u, v, w) = f eq eq[7](u, −v, w), f eq eq[10](u, v, w) = f eq eq[7](−u, −v, w), f eq eq[11](u, v, w) = f eq eq[7](u, w, v), f eq eq[12](u, v, w) = f eq eq[7](−u, w, v), f eq eq[13](u, v, w) = f eq eq[7](u, −w, v), f eq eq[14](u, v, w) = f eq eq[7](−u, −w, v), f eq eq[15](u, v, w) = f eq eq[7](v, w, u), f eq eq[16](u, v, w) = f eq eq[7](−v, w, u), f eq eq[17](u, v, w) = f eq eq[7](v, −w, u), f eq eq[18](u, v, w) = f eq eq[7](−v, −w, u), f eq eq[20](u, v, w) = f eq eq[19](v, u, w), f eq eq[21](u, v, w) = f eq eq[19](w, v, u), f eq eq[22](u, v, w) = f eq eq[19](−u, v, w), 094701-8…”