Progress in the development of a new lattice Boltzmann method. (15th August 2019)

• Record Type: Journal Article
• Authors: So, R.M.C.; Leung, R.C.K.; Kam, E.W.S.; Fu, S.C.
• Is Part Of: Computers & fluids. Volume 190(2019)
• Page Start: 440
• Page End: 469
• Publication Date: 2019-08-15
• Subjects: Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2019.04.009
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: A new LBM that remedies four inadequacies of conventional LBM is formulated. It requires only a relaxation time scale of the flow and fluid properties as inputs. It is solved with finite difference splitting on minimal D2Q9 and D3Q13 lattices. It simulates a wide ranges of flow problems with different fluid media correctly. Its solutions agree excellently with existing experimental and numerical results. Abstract: A new modeled Boltzmann equation (MBE) with four improvements made to conventional MBE is formulated. The first improvement is to include the particle internal rotational degree of freedom in the derivation of a continuous equilibrium velocity distribution function f eq ; thus, rendering the MBE applicable to diatomic gas. The second improvement is made in the expansion assumed for fα eq in the lattice Boltzmann equation (LBE). This expansion is expressed in terms of the particle velocity vector ( ξ ) alone; hence, the LBE is no longer limited by a very low Mach number ( M ) assumption, and it also allows the LBE to correctly satisfy the zero divergence of the velocity field for incompressible flow. The third improvement is made to eliminate the bounce-back rule used to model no-slip wall boundary condition for fα because the rule leads to leakage at solid walls and mass conservation is compromised. The fourth improvement is carried out to render the modeled LBE truly valid for hydrodynamic flow simulation. Thus improved, the new lattice Boltzmann method (LBM) is no longer subject to the M << 1 condition and can be confidently extended to hydrodynamic simulations where the flow is truly incompressible. The fα eq expansion coefficients are found to be functions of the primitive variables, their derivatives and their products. In this LBE, fluid properties are inputs, and boundary values of fα are deduced from the primitive variables at the boundaries. Only a D2Q9 and a D3Q15 lattice model are required for 2-D and 3-D flow simulation, respectively. A finite difference splitting method is used to solve the new LBM; the scheme is labeled FDLBM, and it has been used to simulate widely different laminar flow of gas and liquid; incompressible with/without heat transfer, compressible with/without shocks, aeroacoustics with/without scattering, thermo-aeroacoustics, buoyant and double diffusive flow, and non-Newtonian flow, as well as blood flow in arteries with/without blockage. Accurate results are obtained, and they agree with other finite-difference numerical simulations of the same problems, and experimental and theoretical results whenever available. Furthermore, the aeroacoustics results are in agreement with those obtained from direct aeroacoustics simulations of the same problems.