Burgers-Korteweg-de Vries复合方程的格子Boltzmann方法模拟

针对Burgers-Korteweg-de Vries(cBKdV)复合方程提出一种格子Boltzmann模型.通过恰当地处理色散项u_xxx并运用Chapman-Enskog展开从格子Boltzmann方程推导出宏观方程,从而得到联系微观量与宏观量的局部平衡分布函数.对不同微分方程进行数值实验,数值解与解析解非常吻合,相比于其它数值结果,该格子Boltzmann模型的数值结果更精确,说明该数值模型的高效性.     

二维空腔黏性流的格子Boltzmann方法模拟

用13速六方格子BhatnagarGrossKrook(缩写为BGK)模型模拟二维空腔黏性流.给出了上边界流体作匀速运动时,具有不同雷诺数的空腔黏性流的流场速度分布情况,模拟了在雷诺数Re=3000时,流场中的涡旋形成过程及流场稳定后,腔内密度、压力和温度的分布情况 关键词：     

基于九速四方格子模型的二维对流扩散方程格子Boltzmann方法模拟

用九速正方格子模型给出了二维对流扩散方程的格子Boltzmann方法。由对流扩散方程的对流系数和扩散系数确定了局域平衡分布函数的系数。计算机模拟结果与理论结果吻合很好。     

二维卡门涡街的格子Boltzmann仿真

格子气自动机和格子Boltzmann方法的迅速发展提供了一类求解流体力学问题的新的方法。本文中，我们介绍了Boltzmann方法，解决了格子气方法中的缺点，通过选择适当平衡分布及其参数，导出了Navier－Stokers方程，并得到了声速和粘性系数。最后在微机上模拟了在无限长平板流动问题及绕单分离板的流动问题，得到了卡门涡街。结果表该模型有格子气方法及其它的数值方法所没有的优点，计算更精确、更直观、更有效。     

一维空间Riesz分数阶对流扩散方程的格子Boltzmann方法

建立格子Boltzmann方法（LBM）的D1Q3演化模型,研究一类Riesz空间分数阶对流扩散方程的数值求解问题。对分数阶微积分算子中的积分项离散化处理,得到逼近的标准对流扩散方程。结合Taylor展式和Chapman-Enskog多尺度展开技术得到模型的各个方向上的平衡态分布函数,通过D1Q3演化模型正确恢复所要求解的宏观方程。数值算例验证该方法的有效性。     

十三点格子Boltzmann模型仿真

格子气和格子Boltzmann方法的迅速发展提供了一类求解流体力学问题的新方法。格子Boltzmann方法在保留了格子气模型优点的同时,克服了它的不足之处。本文讨论了一种三迭加HPP十三点模型,通过选择适当的平衡分布及参数,并用Chapman－Enskog展开和多尺度技术导出了Navier－Stokes方程。在微机上模拟了空腔流的流动问题,并与传统方法的计算结果进行了比较,结果表明该模型能较好的模拟复杂流动现象,并具有较好的工程应用背景。     

格子Boltzmann亚格子模型的研究

为了将格子Boltzmann法应用于大雷诺数流动的模拟,本文将Smagorinsky亚格子模型和LBGK模型相结合,并对该亚格子LBM模型进行了研究。利用该亚格子LBM模型,对二维顶盖驱动流进行了模拟,得到了若干大雷诺数下流线图和方腔中心线上无量纲速度分布。计算结果与基准解进行比较,两者相互吻合。     

用格子Boltzmann模型模拟可压缩完全气体流动

采用一种新的格子Boltzmann模型模拟超音速流动。在这种模型中,粒子的速度不受限制,可以取得很广。而平衡分布函数的支集却相对集中,使模型得以简化。粒子速度的这种自适应特性允许流体以较高的马赫数流动。通过引入粒子的势能使得该模型适用于具有任意比热比的完全气体。利用Chapman-Enskog方法,从BGK型Boltzmann方程推导出Navier-Stokes方程。在六边形网格上模拟了马赫数为3的前台阶绕流,得到了合理的结果。     

用格子Boltzmann模型模拟激波现象

根据微观和宏观之间的质量、动量、能量守恒准则,建立了一个两维的D2Q13格子Boltzmann模型,可从该D2Q13模型出发推导出宏观的流体力学方程,所得动量方程和Navier-Stokes方程相比,在黏性输运项上有明显的改进,用该模型对冲击波在障碍物表面上的折射和反射现象的数值模拟都得到了比较满意的结果,而且数值稳定性也很好.证明了D2Q13模型的适宜性 关键词： Boltzmann模型 分布函数 冲击波 流体力学方程     

A higher order lattice BGK model for simulating some nonlinear partial differential equations

In this paper, we consider a one-dimensional nonlinear partial differential equation that has the form ut + αuux + βunux - γuxx + δuxxx = F(u). A higher order lattice Bhatnager-Gross-Krook (BGK) model with an amending-function is proposed. With the Chapman-Enskog expansion, different kinds of nonlinear partial differential equations are recovered correctly from the continuous Boltzmann equation. The numerical results show that this method is very effective.     

A lattice Boltzmann model with an amending function forsimulating nonlinear partial differential equations

This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form $u_t+\alpha uu_{xx}+\beta u^n u_x+\gamma u_{xxx}+\xi u_{xxxx}=0$. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman--Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.     

A lattice Boltzmann model for the generalized Burgers-Huxley equation

In this paper we develop a lattice Boltzmann model for the generalized Burgers-Huxley equation (GBHE). By choosing the proper time and space scales and applying the Chapman-Enskog expansion, the governing equation is recovered correctly from the lattice Boltzmann equation, and the local equilibrium distribution functions are obtained. Excellent agreement with the exact solution is observed, and better numerical accuracy is obtained than the available numerical result. The results indicate the present model is satisfactory and efficient. The method can also be applied to the generalized Burgers-Fisher equation and be extended to multidimensional cases.     

Lattice Boltzmann method for the generalized Kuramoto-Sivashinsky equation

In this paper, a lattice Boltzmann model with an amending function is proposed for the generalized Kuramoto-Sivashinsky equation that has the form u_t+uu_x+αu_xx+βu_xxx+γu_xxxx=0. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. It is found that the numerical results agree well with the analytical solutions.     

Global solutions of the Boltzmann equation on a lattice

The nonlinear Boltzmann equation with a discretized spatial variable is studied in a Banach space of absolutely integrable functions of the velocity variables. Conservation laws and positivity are utilized to extend weak local solutions to a global solution. This is shown to be a strong solution by analytic semigroup techniques.Supported by National Science Foundation Grant ENG-7515882.     

Thermal equation of state for lattice Boltzmann gases

The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar--Gross--Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar--Gross--Krook equation (LBGKE) induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics.     

声悬浮过程的格子Boltzmann方法研究

采用轴对称多弛豫时间格子Boltzmann(LB)方法,研究了圆柱形封闭谐振腔中圆盘形样品的声悬浮过程.模拟结果表明,(001)模式下谐振腔的共振长度L=0.499λ,在谐振腔中心引入样品后共振漂移量δL≈-0.9,这与线性声学理论计算结果基本相符.声悬浮力的LB模拟过程包含了黏滞性效应和共振漂移效应,所获得的模拟结果与理论公式计算值在量值上一致,而且其在细节上更符合实验现象.此外,LB模拟还揭示出了声悬浮过程中的声压波形畸变、声流和声辐射压等非线性声学效应.     

Local second-order boundary methods for lattice Boltzmann models

A new way to implement solid obstacles in lattice Boltzmann models is presented. The unknown populations at the boundary nodes are derived from the locally known populations with the help of a second-order Chapman-Enskog expansion and Dirichlet boundary conditions with a given momentum. Steady flows near a flat wall, arbitrarily inclined with respect to the lattice links, are then obtained with a third-order error. In particular, Couette and Poiseuille flows are exactly recovered without the Knudsen layers produced for inclined walls by the bounce back condition.