一类非线性偏微分方程的四阶格子Boltzmann模型

通过Chapman-Enskog展开技术和多尺度分析,建立了一种新的D1Q4带修正项的四阶格子Boltzmann模型,一类非线性偏微分方程从连续的Boltzmann方程得到正确恢复．统一了KdV和Burgers等已知方程类型的格子BGK模型,还首次给出了组合KdV-Burgers,广义Burgers—Huxley等方程...     

Lattice Boltzmann solver of Rossler equation

We proposed a lattice Boltzmann model for the Rossler equation. Using a method of multiscales in the lattice Boltzmann model, we get the diffusion reaction as a special case. If the diffusion effect disappeared, we can obtain the lattice Boltzmann solution of the Rossler equation on the mesescopic scale. The numerical results show the method can be used to simulate Rossler equation.     

Numerical simulation for the solitary wave of Zakharov–Kuznetsov equation based on lattice Boltzmann method

The Zakharov–Kuznetsov equation is considered, which is an equation describing two dimensional weakly nonlinear ion-acoustic waves in plasma. We focus on using the lattice Boltzmann method to study the Zakharov–Kuznetsov equation. A lattice Boltzmann model is constructed. In numerical experiments, the propagation of the single solitary wave and the collision of double solitary waves are simulated. The results with different parameters are investigated and compared.     

A higher-order moment method of the lattice Boltzmann model for the conservation law equation

In this paper, we proposed a higher-order moment method in the lattice Boltzmann model for the conservation law equation. In contrast to the lattice Bhatnagar–Gross–Krook (BGK) model, the higher-order moment method has a wide flexibility to select equilibrium distribution function. This method is based on so-called a series of partial differential equations obtained by using multi-scale technique and Chapman–Enskog expansion. According to Hirt’s heuristic stability theory, the stability of the scheme can be controlled by modulating some special moments to design the third-order dispersion term and the fourth-order dissipation term. As results, the conservation law equation is recovered with higher-order truncation error. The numerical examples show the higher-order moment method can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the conservation law equation.     

Lattice evolution solution for the nonlinear Poisson–Boltzmann equation in confined domains

The lattice evolution method for solving the nonlinear Poisson–Boltzmann equation in confined domain is developed by introducing the second-order accurate Dirichlet and Neumann boundary implements, which are consistent with the non-slip model in lattice Boltzmann method for fluid flows. The lattice evolution method is validated by comparing with various analytical solutions and shows superior to the classical numerical solvers of the nonlinear Poisson equations with Neumann boundary conditions. The accuracy and stability of the method are discussed. This lattice evolution nonlinear Poisson–Boltzmann solver is suitable for efficient parallel computing, complex geometry conditions, and easy extension to three-dimensional cases.     

Lattice Boltzmann model for two-dimensional generalized sine-Gordon equation

The nonlinear sine-Gordon equation arises in various problems in science and engineering. In this paper, we propose a numerical model based on lattice Boltmann method to obtain the numerical solutions of two-dimensional generalized sine-Gordon equation, including damped and undamped sine-Gordon equation. By choosing properly the conservation condition between the macroscopic quantity $u_t$ and the distribution functions and applying the Chapman-Enskog expansion, the governing equation is recovered correctly from the lattice Boltzmann equation. Moreover, the local equilibrium distribution function is obtained. The numerical results of the first three examples agree well with the analytic solutions, which indicates the lattice Boltzmann model is satisfactory and efficient. Numerical solutions for cases involving the most known from the bibliography line and ring solitons are given. Numerical experiments also show that the present scheme has a good long-time numerical behavior for the generalized sine-Gordon equation. Moreover, the model can also be applied to other two-dimensional nonlinear wave equations, such as nonlinear hyperbolic telegraph equation and Klein-Gordon equation.     

General propagation lattice Boltzmann model for variable-coefficient non-isospectral KdV equation

In this paper, a general propagation lattice Boltzmann model for variable-coefficient non-isospectral Korteweg–de Vries (vc-nKdV) equation, which can describe the interfacial waves in a two layer liquid and Alfvén waves in a collisionless plasma, is proposed by selecting appropriate equilibrium distribution function and adding the compensate function. The Chapman–Enskog analysis shows that the vc-nKdV equation can be recovered correctly from the present model. Numerical simulation for the non-propagating one soliton of this equation in different situations is conducted as validation. It is found that the numerical results match well with the analytical solutions, which demonstrates that the current general propagation lattice Boltzmann model is a satisfactory and efficient method, and could be more stable and accurate than the standard lattice Bhatnagar–Gross–Krook model.     

Lattice Boltzmann 3D flow simulations on a metallic foam

In this paper a lattice Boltzmann method (LBM) is used to simulate isothermal incompressible flow in a RCM-NCX-1116 metallic foam. The computational technique is a multiple relaxation time (MRT) lattice Boltzmann equation model. Computer aided X-ray micro-tomography is used to obtain 3D images of the metallic foam, providing the geometry and information required for LB simulations of a single phase flow.Pressure drops are computed and successfully compared to experimental measures and correlated with Ergun’s equation. Invariance of Ergun’s parameters A and B with the sampling rate of the images is observed.     

Convergence of lattice Boltzmann methods for Navier–Stokes flows in periodic and bounded domains

Combining an asymptotic analysis of the lattice Boltzmann method with a stability estimate, we are able to prove some convergence results which establish a strict relation to the incompressible Navier–Stokes equation. The proof applies to the lattice Boltzmann method in the case of periodic domains and for specific bounded domains if the Dirichlet boundary condition is realized with the bounce back rule.     

Application of lattice Boltzmann method for incompressible viscous flows

Because of the presence of corner eddies that change in number and pattern the lid-driven cavity problem has been found suitable to study various aspects of the performance of solution algorithms for incompressible viscous flows. It retains all the difficult flow physics and is characterized by a large primary eddy at the centre and secondary eddies located near the cavity corners. In this work, lid-driven cavity flow is simulated by lattice Boltzmann method with single-relaxation-time and it is compared with those by lattice Boltzmann method with multi-relaxation-time and finite difference method. The effects of the Reynolds number on the size, centre position and number of vortices are studied in detail together with the flow pattern in the cavity. The close agreement of the results bears testimony to the validity of this relatively new approach. However lattice Boltzmann method with multi-relaxation-time model is seen to remove the difficulties faces by the lattice Boltzmann method with single-relaxation-time at higher Reynolds numbers.     

均匀格子气Boltzmann方程的稳定性研究

１．引言格子气的基本方程是在几何空间、速度空间和时间上都是离散的Ｂｏｌｔｚｍａｎｎ方程（Ｂ方程）．这是一个有限差分方程．在离散速度气体运动论中［１］,Ｂ方程在速度空间上是离散的,在几何空间和时间上是连续的．这是一个偏微分方程．人们对离散速度气体Ｂ方程的稳定性和渐近特性的研究已经取得了很多结果．Ｍａａｓｓ~［２］通过构造Ｌｙaｐｕｎｏｖ函数族,在分布函数在空间上均匀的条件下,证明了平衡分布的渐近稳定性．信息函数Ｈ是该函数族的一员．Ｂｅｌｌｏｍｏｅｔａｌ~［３］．采用小扰动线性化方法在初值距离平衡解足…     

Lattice Boltzmann simulation for high-speed compressible viscous flows with a boundary layer

In this study, the lattice Boltzmann method is employed for simulating high-speed compressible viscous flows with a boundary layer. The coupled double-distribution-function lattice Boltzmann method proposed by Li et al. (2007) is employed because of its good numerical stability and non-free-parameter feature. The non-uniform mesh construction near the wall boundary in fine grids is combined with an appropriate wall boundary treatment for the finite difference method in order to obtain accurate spatial resolution in the boundary layer problem. Three typical problems in high-speed viscous flows are solved in the lattice Boltzmann simulation, i.e., the compressible boundary layer problem, shock wave problem, and shock boundary layer interaction problem. In addition, in-depth comparisons are made with the non-oscillatory and non-free-parameter dissipation (NND) scheme and second order upwind scheme in the present lattice Boltzmann model. Our simulation results indicate the great potential of the lattice Boltzmann method for simulating high-speed compressible viscous flows with a boundary layer. Further research is needed (e.g., better numerical models and appropriate finite difference schemes) because the lattice Boltzmann method is still immature for high-speed compressible viscous flow applications.     

Lattice Boltzmann model for the bimolecular autocatalytic reaction–diffusion equation

A lattice Boltzmann model for the bimolecular autocatalytic reaction–diffusion equation is proposed. By using multi-scale technique and the Chapman–Enskog expansion on complex lattice Boltzmann equation, we obtain a series of complex partial differential equations, complex equilibrium distribution function and its complex moments. Then, the complex reaction–diffusion equation is recovered with higher-order accuracy of the truncation error. This equation can be used to describe the bimolecular autocatalytic reaction–diffusion systems, in which a rich variety of behaviors have been observed. Based on this model, the Fitzhugh–Nagumo model and the Gray–Scott model are simulated. The comparisons between the LBM results and the Alternative Direction Implicit results are given in detail. The numerical examples show that assumptions of source term can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the complex reaction–diffusion equation.     

Numerical solutions of two-dimensional Burgers’ equations by lattice Boltzmann method

In this paper, the two-dimensional Burgers’ equations with two variables are solved numerically by the lattice Boltzmann method. The lattice Bhatnagar–Gross–Krook model we used can recover the macroscopic equation with the second order accuracy. Numerical solutions for various values of Reynolds number, computational domain, initial and boundary conditions are calculated and validated against exact solutions or other published results. It is concluded that the proposed method performs well.     

General propagation lattice Boltzmann model for a variable-coefficient compound KdV-Burgers equation

In this paper, a general propagation lattice Boltzmann model for a variable-coefficient compound Korteweg-de Vries-Burgers (vc-cKdVB) equation is investigated through selecting equilibrium distribution function and adding a compensation function, which can provide some more realistic models than their constant-coefficient counterparts in fluids or plasmas. Chapman–Enskog analysis shows that the vc-gKdVB equation can be recovered correctly from the present model. Numerical simulations in different situations of this equation are conducted, including the propagation and interaction of the bell-type, kink-type and periodic-depression solitons and the evolution of the shock-wave solutions. It is found that the numerical results match well with the analytical solutions, which demonstrates that the current lattice Boltzmann model is a satisfactory and efficient algorithm. In addition, it is also shown the present model could be more stable and more accurate than the standard lattice Bhatnagar–Gross–Krook model through adjusting the two free parameters introduced into the propagation step.     

Lattice Boltzmann model for two-dimensional unsteady Burgers’ equation

In this paper, a special lattice Boltzmann model is proposed to simulate two-dimensional unsteady Burgers’ equation. The maximum principle and the stability are proved. The model has been verified by several test examples. Excellent agreement is obtained between numerical predictions and exact solutions. The cases of steep oblique shock waves are solved and compared with the two-point compact scheme results. The study indicates that lattice Boltzmann model is highly stable and efficient even for the problems with severe gradients.     

Lattice Boltzmann model for a generalized Gardner equation with time-dependent variable coefficients

In this paper, lattice Boltzmann model for a generalized Gardner equation with time-dependent variable coefficients, which can provide some more realistic models than their constant-coefficient counterparts, is derived through selecting equilibrium distribution function and adding the compensate function, appropriately. Effects and approximate value range of the free parameters, which are introduced to adjust the single relaxation time and equilibrium distribution function, are discussed in detail, as well as the impact of the lattice space step and velocity. Numerical simulations in different situations of this equation are conducted, including the propagation and interaction of the solitons, the evolution of the non-propagating soliton and the propagation of the double-pole solutions. It is found that the numerical results match well with the analytical solutions, which demonstrates that the current lattice Boltzmann model is a satisfactory and efficient algorithm.     

Boltzmann方程的精确解

本文建立了非线性Boltzmann输运方程与线性AKNS方程之间的直接联系,从而以不同于Chapman,Enskog和Grad的方式,使Boltzmann方程化归为Dirac方程的求解。在没有其他外场作用的情形中,Boltzmann方程的精确解可以用逆散射方法来求得。     

A finite difference interpretation of the lattice Boltzmann method

Compared to conventional techniques in computational fluid dynamics, the lattice Boltzmann method (LBM) seems to be a completely different approach to solve the incompressible Navier–Stokes equation. The aim of this article is to correct this impression by showing the close relation of LBM to two standard methods: relaxation schemes and explicit finite difference discretizations. As a side effect, new starting points for a discretization of the incompressible Navier–Stokes equation are obtained. © John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:383 402, 2001     

Lattice Boltzmann Constraints for Standard Optimization Problems

This paper presents a simultaneous analysis and design (SAND) approach for the control of incompressible, stationary fluid problems ruled by the lattice Boltzmann equation. Therefore the discretized lattice Boltzmann fix point problem is used as constraint for the optimization problem. This allows for the use of standard non-linear programming techniques. To validate the proposed method a topology optimization problem is solved with the NLP solver WORHP at moderate Reynolds numbers (Re < 1000) for a 2-dimensional flow domain. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)