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.     

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 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.     

A novel lattice Boltzmann model for the Poisson equation

In this paper, a novel lattice Boltzmann model is proposed to solve the Poisson equation through modifying equilibrium distribution function. Compared with previous models, which can be viewed as the solvers to diffusion equation, the present model is a genuine solver to the Poisson equation, and the transient term derived by previous models is eliminated. Numerical solutions agree well with analytical solutions, which indicates the potential of the present model for solving the Poisson equation.     

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.     

A lattice Boltzmann model for the eddy–stream equations in two-dimensional incompressible flows

A lattice Boltzmann model for two-dimensional incompressible flows with eddy–stream equations is proposed. By using two kinds of distribution functions and employing several higher-order moments of equilibrium distribution functions, the eddy equation and stream function equation with the second-order truncation error are obtained. In the numerical examples, we compared the numerical results of this scheme with those obtained by other classical method. The numerical results agree well with the classical ones.     

A new lattice Boltzmann model for the laplace equation

In this paper, a new lattice Boltzmann equation which is independent of time is proposed. Based on the new lattice Boltzmann equation, some steady problems can be modeled by the lattice Boltzmann method. In the further study, the Laplace equation is investigated with the method of the higher-order moment of equilibrium distribution functions and a series of partial differential equations in different space scales. The numerical results show that the new method is effective.     

A lattice Boltzmann model for blood flows

A lattice Boltzmann model for blood flows is proposed. The lattice Boltzmann Bi-viscosity constitutive relations and control dynamics equations of blood flow are presented. A non-equilibrium phase is added to the equilibrium distribution function in order to adjust the viscosity coefficient. By comparison with the rheology models, we find that the lattice Boltzmann Bi-viscosity model is more suitable to study blood flow problems. To demonstrate the potential of this approach and its suitability for the application, based on this validate model, as examples, the blood flow inside the stenotic artery is investigated.     

A lattice Boltzmann model for Maxwell’s equations

In this paper, a lattice Boltzmann model for the Maxwell’s equations is proposed by taking separate sets of distribution functions for the electric and magnetic fields, and a lattice Boltzmann model for the Maxwell vorticity equations with third order truncation error is proposed by using the higher-order moment method. At the same time, the expressions of the equilibrium distribution function and the stability conditions for this model are given. As numerical examples, some classical electromagnetic phenomena, such as the electric and magnetic fields around a line current source, the electric field and equipotential lines around an electrostatic dipole, the electric and magnetic fields around oscillating dipoles are given. These numerical results agree well with classical ones.     

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.     

用格子Boltzmann方法模拟非线性热传导方程

对具有非线性源项和非线性扩散项的热传导方程建立格子Boltzmann求解模型.在演化方程中增加了两个关于源项分布函数的微分算子,对演化方程实施Chapman-Enskog展开.通过对演化方程的进一步改进,恢复出具有高阶截断误差的宏观方程.对不同参数选取下的非线性热传导方程进行了数值模拟,数值解与精确解吻合得很好.该模型也可以用于同类型的其他偏微分方程的数值计算中.     

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

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

Equilibration Rate for the Linear Inhomogeneous Relaxation-Time Boltzmann Equation for Charged Particles

Abstract

We study the long-time behavior of a linear inhomogeneous Boltzmann equation. The collision operator is modeled by a simple relaxation towards the Maxwellian distribution with zero mean and fixed lattice temperature. Particles are moving under the action of an external potential that confines particles, i.e., there exists a unique stationary probability density. Convergence rate towards global equilibrium is explicitly measured based on the entropy dissipation method and apriori time independent estimates on the solutions. We are able to prove that this convergence is faster than any algebraic time function, but we cannot achieve exponential convergence.     

A lattice Boltzmann model for the Fokker-Planck equation

In this paper, a lattice Boltzmann model is presented for solving one and two-dimensional Fokker-Planck equations with variable coefficients. In particular, it is efficient to simulate one-dimensional stochastic processes governed by the Fokker-Planck equation. Numerical results agree well with the exact solutions, which indicates that the proposed model is suitable for solving the Fokker-Planck 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.     

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

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

The linearized Boltzmann equation: Sound-wave propagation in a rarefied gas

An analytical version of the discrete-ordinates method (the ADO method) is used to establish concise and particularly accurate solutions to the problem of sound-wave propagation in a rarefied gas. The analysis and the numerical work are based on a rigorous form of the linearized Boltzmann equation (for rigid-sphere interactions), and in contrast to many other works formulated (for an infinite medium) without a boundary condition, the solution reported here satisfies a boundary condition that models a diffusely-reflecting vibrating plate. In addition and in order to investigate the effect of kinetic models, solutions are developed for the BGK model, the S model, the Gross-Jackson model, as well as for the (newly defined) MRS model and the CES model. While the developed numerical results are compared to available experimental data, emphasis in this work is placed on the solutions of the problem of sound-wave propagation as described by the linearized Boltzmann equation and the five considered kinetic models. Received: November 22, 2004; revised: February 24, 2005     

KINETIC FLUX VECTOR SPLITTING FOR THE EULER EQUATIONS WITH GENERAL PRESSURE LAWS

This paper attempts to develop kinetic flux vector splitting(KFVS)for the Euler equa-tions with general pressure laws.It is well known that the gas distribution function forthe local equilibrium state plays an important role in the construction of the gas-kineticschemes.To recover the Euler equations with a general equation of state(EOS),a newlocal equilibrium distribution is introduced with two parameters of temperature approx-imation decided uniquely by macroscopic variables.Utilizing the well-known connectionthat the Euler equations of motion are the moments of the Boltzmann equation wheneverthe velocity distribution function is a local equilibrium state,a class of high resolutionMUSCL-type KFVS schemes are presented to approximate the Euler equations of gas dy-namics with a general EOS.The schemes are finally applied to several test problems for ageneral EOS.In comparison with the exact solutions,our schemes give correct location andmore accurate resolution of discontinuities.The extension of our idea to multidimensionalcase is natural.     

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.     

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.