基于格子Boltzmann方法的超声速非平衡流模拟

基于D1Q4可压缩格子Boltzmann模型,按照流通矢量分裂方法的思路,采用坐标旋转技术构造求解三维带化学反应Navier-Stokes方程对流通量求解器.结合有限体积法求解三维化学非平衡流Navier-Stokes方程,采用时间算子分裂算法解决化学反应刚性问题,数值模拟超声速化学非平衡流的三个经典算例.数值结果表明:在高马赫数下,采用D1Q4可压缩格子Boltzmann模型构造的三维对流通量求解器数值模拟中没有出现非物理解,同时在超声速化学非平衡流场中正确分辨激波、燃烧波等物理现象,精度和分辨率均较高,验证了本文构造的三维对流通量求解器的可靠性,拓宽了D1Q4可压缩格子Boltzmann模型的应用范围,为计算超声速化学非平衡流提供一种新方法.     

Three-dimensional coupled double-distribution-function lattice Boltzmann models for compressible Navier–Stokes equations

Two three-dimensional (3D) lattice Boltzmann models in the framework of coupled double-distribution-function approach for compressible flows, in which specific-heat ratio and Prandtl number can be adjustable, are developed in this paper. The main differences between the two models are discrete equilibrium density and total energy distribution function. One is the D3Q25 model obtained from spherical function, and the other is the D3Q27 standard lattice model obtained from Hermite expansions of the corresponding continuous equilibrium distribution functions. The two models are tested by numerical simulations of some typical compressible flows, and their numerical stability and precision are also analysed. The results indicate that the two models are capable for supersonic flows, while the one from Hermite expansions is not suitable for compressible flows with shock waves.     

A Hybrid Lattice Boltzmann Flux Solver for Simulation of Viscous Compressible Flows

In this paper, a hybrid lattice Boltzmann flux solver (LBFS) is proposed for simulation of viscous compressible flows. In the solver, the finite volume method is applied to solve the Navier-Stokes equations. Different from conventional Navier-Stokes solvers, in this work, the inviscid flux across the cell interface is evaluated by local reconstruction of solution using one-dimensional lattice Boltzmann model, while the viscous flux is still approximated by conventional smooth function approximation. The present work overcomes the two major drawbacks of existing LBFS [28–31], which is used for simulation of inviscid flows. The first one is its ability to simulate viscous flows by including evaluation of viscous flux. The second one is its ability to effectively capture both strong shock waves and thin boundary layers through introduction of a switch function for evaluation of inviscid flux, which takes a value close to zero in the boundary layer and one around the strong shock wave. Numerical experiments demonstrate that the present solver can accurately and effectively simulate hypersonic viscous flows.     

A time-reversal lattice Boltzmann method

In this paper we address the time-reversed simulation of viscous flows by the lattice Boltzmann method (LB). The theoretical derivation of the reversed LB from the Boltzmann equation is detailed, and the method implemented for weakly compressible flows using the D2Q9 scheme. The implementation of boundary conditions is also discussed. The accuracy and stability are illustrated by four test cases, namely the propagation of an acoustic wave in a medium at rest and in an uniform mean flow, the Taylor–Green vortex decay and the vortex pair–wall collision.     

Three-dimensional finite-difference lattice Boltzmann model and its application to inviscid compressible flows with shock waves

In this paper, a three-dimensional (3D) finite-difference lattice Boltzmann model for simulating compressible flows with shock waves is developed in the framework of the double-distribution-function approach. In the model, a density distribution function is adopted to model the flow field, while a total energy distribution function is adopted to model the temperature field. The discrete equilibrium density and total energy distribution functions are derived from the Hermite expansions of the continuous equilibrium distribution functions. The discrete velocity set is obtained by choosing the abscissae of a suitable Gauss–Hermite quadrature with sufficient accuracy. In order to capture the shock waves in compressible flows and improve the numerical accuracy and stability, an implicit–explicit finite-difference numerical technique based on the total variation diminishing flux limitation is introduced to solve the discrete kinetic equations. The model is tested by numerical simulations of some typical compressible flows with shock waves ranging from 1D to 3D. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.     

一个新的用于不可压磁流体动力学的格子波尔兹曼模型

绝大多数现有的格子波尔兹曼磁流体动力学模型其实是用可压缩方法来模拟不可压磁流体。而这些可压缩效应在数值模拟中往往会带来意想不到的误差。在这篇文章中，我们提出了一个全新的可用于的不可压格子波尔兹曼磁流体动力学模型，并且进行了哈特曼流的数值模拟。模拟结果与哈特曼流的解析解非常吻合。这个方法需要一个假设条件来消除误差。我们做了大量的数值试验，并且与Dellar教授的模型进行了详细的分析与比较。     

A multiple-relaxation-time lattice Boltzmann method for high-speed compressible flows

This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time(SRT) lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations,it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.     

Consistent lattice Boltzmann method

Lack of energy conservation in lattice Boltzmann models leads to unrealistically high values of the bulk viscosity. For this reason, the lattice Boltzmann method remains a computational tool rather than a model of a fluid. A novel lattice Boltzmann model with energy conservation is derived from Boltzmann's kinetic theory. Simulations demonstrate that the new lattice Boltzmann model is the valid approximation of the Boltzmann equation for weakly compressible flows and microflows.     

Lattice Bhatnagar-Gross-Krook Simulations in 2-D Incompressible Magnetohydrodynamics

Lattice Boltzmann Method is recently developed within numerical schemes for simulating a variety of physical systems. In this paper a new lattice Bhatnagar-Gross-Krook (LBGK) model for two-dimensional incompressible magnetohydrodynamics (IMHD) is presented. The model is an extension of a hydrodynamics lattice BGK model with 9 velocities on a square lattice, resulting in a model with 17 velocities. Most of the existing LBGK models for MHD can be viewed as compressible schemes to simulate incompressible flows. The compressible effect might lead to some undesirable errors in numerical simulations. In our model the compressible effect has been overcome successfully. The model is then applied to the Hartmann flow, giving reasonable results.     

Lattice Bhatnagar-Gross-Krook Simulations in 2-D Incompressible Magnetohydrodynamics

Lattice Boltzmann Method is recently developed within numerical schemes for simulating a variety of physical systems. In this paper a new lattice.Bhatnagar-Gross-Krook （LBGK） model for two-dimensional incompressible magnetohydrodynamics （IMHD） is presented. The model is an extension of a hydrodynamics lattice BGK model with 9 velocities on a square lattice, resulting in a model with 17 velocities. Most of the existing LBGK models for MHD can be viewed as compressible schemes to simulate incompressible flows. The compressible effect might lead to some undesirable errors in numerical simulations. In our model the compressible effect has been overcome successfully. The model is then applied to the Hartmann flow, giving reasonable results.     

Nonexistence of H Theorem for Some Lattice Boltzmann Models

In this paper, we provide a set of sufficient conditions under which a lattice Boltzmann model does not admit an H theorem. By verifying the conditions, we prove that a number of existing lattice Boltzmann models does not admit an H theorem. These models include D2Q6, D2Q9 and D3Q15 athermal models, and D2Q16 and D3Q40 thermal (energy-conserving) models. The proof does not require the equilibria to be polynomials.     

Non-hydrodynamic modes and general equations of state in lattice Boltzmann equations

The lattice Boltzmann equation is commonly used to simulate fluids with isothermal equations of state in a weakly compressible limit, and intended to approximate solutions of the incompressible Navier–Stokes equations. Due to symmetry requirements there are usually more degrees of freedom in the equilibrium distributions than there are constraints imposed by the need to recover the Navier–Stokes equations in a slowly varying limit. We construct equilibria for general barotropic fluids, where pressure depends only upon density, using the two-dimensional, nine velocity (D2Q9) and one-dimensional, five velocity (D1Q5) lattices, showing that one otherwise arbitrary function in the equilibria must be chosen to suppress instability.     

一个三维多松弛时间全速域格子Boltzmann模型

构建一个既适用于低速不可压流体又适用于高速可压缩流体的三维自由参数多松弛时间格子Boltzmann模型.模型中,根据SO（3）群的不可约表述基函数构造转化矩阵,根据恢复可压Navier-Stokes方程的需要选取非守恒矩平衡值.通过von Neumann稳定性分析模型参数对数值稳定性的影响,并给出建议选择范围.模型经过基准问题的验证,模拟结果与解析解及其它数值结果符合较好.     

Aeroacoustic Simulations Using Compressible Lattice Boltzmann Method

This paper presents a lattice Boltzmann (LB) method based study aimed at numerical simulation of aeroacoustic phenomenon in flows around a symmetric obstacle. To simulate the compressible flow accurately, a potential energy double-distribution-function (DDF) lattice Boltzmann method is used over the entire computational domain from the near to far fields. The buffer zone and absorbing boundary condition is employed to eliminate the non-physical reflecting. Through the direct numerical simulation, the flow around a circular cylinder at $Re$=150, $M$=0.2 and the flow around a NACA0012 airfoil at $Re$=10000, $M$=0.8, $α$=$0^◦$ are investigated. The generation and propagation of the sound produced by the vortex shedding are reappeared clearly. The obtained results increase our understanding of the characteristic features of the aeroacoustic sound.     

Study of the 1D lattice Boltzmann shallow water equation and its coupling to build a canal network

The D1Q3 lattice Boltzmann (LB) shallow water equation is analyzed in detail and compared with other numerical schemes. Analytical results are derived and used to discuss the accuracy and stability of the model. We show how such D1Q3 LB models for canal reaches may be easily coupled with various hydraulic interconnection structures to build models of complex irrigation networks.     

Quasi-equilibrium lattice Boltzmann method

A general lattice Boltzmann method for simulation of fluids with tailored transport coefficients is presented. It is based on the recently introduced quasi-equilibrium kinetic models, and a general lattice Boltzmann implementation is developed. Lattice Boltzmann models for isothermal binary mixtures with a given Schmidt number, and for a weakly compressible flow with a given Prandtl number are derived and validated.     

Analytical solution of axi-symmetrical lattice Boltzmann model for cylindrical Couette flows

Analytical solution for the axi-symmetrical lattice Boltzmann model is obtained for the low-Mach number cylindrical Couette flows. In the hydrodynamic limit, the present solution is in excellent agreement with the result of the Navier–Stokes equation. Since the kinetic boundary condition is used, the present analytical solution using nine discrete velocities can describe flows with the Knudsen number up to 0.1. Meanwhile, the comparison with the simulation data obtained by the direct simulation Monte Carlo method shows that higher-order lattice Boltzmann models with more discrete velocities are needed for highly rarefied flows.     

Multiple-relaxation-time lattice Boltzmann model for compressible fluids

We present an energy-conserving multiple-relaxation-time finite difference lattice Boltzmann model for compressible flows. The collision step is first calculated in the moment space and then mapped back to the velocity space. The moment space and corresponding transformation matrix are constructed according to the group representation theory. Equilibria of the nonconserved moments are chosen according to the need of recovering compressible Navier-Stokes equations through the Chapman-Enskog expansion. Numerical experiments showed that compressible flows with strong shocks can be well simulated by the present model. The new model works for both low and high speeds compressible flows. It contains more physical information and has better numerical stability and accuracy than its single-relaxation-time version.     

A solution-adaptive lattice Boltzmann method for two-dimensional incompressible viscous flows

A stencil adaptive lattice Boltzmann method (LBM) is developed in this paper. It incorporates the stencil adaptive algorithm developed by Ding and Shu [26] for the solution of Navier–Stokes (N–S) equations into the LBM calculation. Based on the uniform mesh, the stencil adaptive algorithm refines the mesh by two types of 5-points symmetric stencils, which are used in an alternating sequence for increased refinement levels. The two types of symmetric stencils can be easily combined to form a 9-points symmetric structure. Using the one-dimensional second-order interpolation recently developed by Wu and Shu [27] along the straight line and the D2Q9 model, the adaptive LBM calculation can be effectively carried out. Note that the interpolation coefficients are only related to the lattice velocity and stencil size. Hence, the simplicity of LBM is not broken down and the accuracy is maintained. Due to the use of adaptive technique, much less mesh points are required in the simulation as compared to the standard LBM. As a consequence, the computational efficiency is greatly enhanced. The numerical simulation of two dimensional lid-driven cavity flows is carried out. Accurate results and improved efficiency are reached. In addition, the steady and unsteady flows over a circular cylinder are simulated to demonstrate the capability of proposed method for handling problems with curved boundaries. The obtained results compare well with data in the literature.     

Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.