颗粒流体系统的格子Boltzmann数值方法研究进展

介绍Euler-Lagrange框架下基于格子Boltzmann方法 LBM(Lattice Boltzmann Method)发展的两种不同层次（即不同时-空尺度和精度）的颗粒流体系统离散模拟方法,即格子Boltzmann颗粒解析直接数值模拟（LB-based PR-DNS）方法和格子Boltzmann离散颗粒模拟（LB-based DPS）方法,总结了Euler-Euler框架下基于格子Boltzmann双流体模型（LB-based TFM）方面的探索研究。LB-based PR-DNS方法中颗粒尺寸远大于格子步长,能够直接解析出流体在颗粒表面的流动以及颗粒所受完整的动力学信息;LB-based DPS方法中格子步长远大于颗粒直径,其在计算精度、时间耗费和计算效率之间能达到很好的平衡,可获得流体的宏观平均流动及颗粒的运动轨迹信息。LB-based DNS和DPS是探索颗粒流体系统的有力手段,但LB-based TFM应用于模拟颗粒流体系统仍需进一步探索。     

利用格子Boltzmann方法数值模拟Rayleigh-Benard对流

采用格子Boltzmann方法对较大Rayleigh数范围下的二维Rayleigh-Benard对流进行了模拟研究.引入能量分布函数,利用该能量分布函数与粒子速度分布函数耦合来求解一个热流场,能量分布函数与粒子速度分布函数和Boltzmann方程构成了一个新的双分布格子Boltzmann模型.在考虑密度随温度变化的情况下,进行数值模拟,得到了Rayleigh-Benard对流速度、温度随时间的变化规律、系统的流线和等温线分布及平均Nusselt数与Rayleigh数的之间的关系,与相关文献数据进行了对比,模拟结果非常吻合,证明了改进的双分布格子Boltzmann模型的有效性.     

并列圆柱绕流的格子Boltzmann数值模拟

采用非均匀不可压格子Boltzmann模型对低雷诺数下并列圆柱绕流进行了数值模拟,给出了数值计算结果,分析了间距g对圆柱尾流及升力、阻力的影响,并在此基础上得到了4种尾迹模式.此外,研究了流场的初始扰动对流动分岔现象的影响,发现在适当的扰动下可以很快得到同步同相的尾流.对Re=160和200下圆柱的升、阻力进行了对比,结果表明升力和阻力受间距g的影响大于雷诺数.     

用格子Boltzmann方法数值模拟混合层中旋涡的合并过程

用格子Boltzmann方法计算混合层中的流动问题。在流场的入口处加不同频率、振幅和相位的小扰动,观察混合层中旋涡的演进机理,模拟二维混合层中旋涡合并现象。在基本扰动波的基础上,又加入频率为基本波频率一半的亚谐波,得到了两个涡合并的计算结果,当加入的亚谐波频率为基本波频率的三分之一时,得到了三个涡合并的计算结果。这些计算结果与已有文献的结果基本一致,显示用格子Boltzmann方法模拟混合层问题是可行的。     

Ginzburg—Landau方程的格子Boltzmann解

本文给出了用于求解Ginzburg-Landau方程的格子Boltzmann模型。通过构造平衡态分布函数的矩函数，我们给出了宏观物理量与平衡态分布函数之间的变换。作为一种特殊情况，得到了各向同性的平衡态分布函数，这后，我们证明了该算法收敛到Ginzburg-Landau方程。     

利用格子Boltzmann方法模拟单个气泡在复杂流道内的运动特性

利用格子Boltzmann方法模拟了单个气泡在具有三个半圆形喉部的复杂流道内的上升过程.通过分析气泡运动过程中的形态、运动轨迹及运动速度的变化,研究复杂流道对气泡运动特性的影响.研究结果表明,在上升过程中,由于壁面的影响,气泡的形状发生严重的变形,运动轨迹也发生相应的偏转.通过实验结果的对比,验证了模拟结果的正确性.结果表明格子Boltzmann方法可以用于模拟具有复杂边界的两相流问题.     

等离子熔射粉末颗粒飞行过程格子Boltzmann法仿真

为考察等离子熔射过程中粉末的飞行过程,本文在已开发的正六边形7-b it格子Bo ltzm ann(LB)方法等离子射流的温度场和速度场的计算模型基础上,采用单个颗粒加速方程,建立了一个随机算法,实现了对粉末颗粒在射流场中运动过程的仿真;计算结果通过动画演示了粉末飞行的全过程,表明初始位置越靠近射流场出口中心的粉末颗粒加速越充分,并且在射流场一定的情况下,减小粉末颗粒直径可以提高粉末速度,但会降低粉末利用率。     

基于格子Boltzmann方法的挤出胀大的数值模拟

将单相格子Boltzmann方法(lattice Boltzmann method, LBM)引入到粘弹流体的瞬态挤出胀大的数值模拟中,建立了基于双分布函数的自由面粘弹性流动格子Boltzmann模型．分析得到的流道中流动速度分布和构型张量结果与理论解十分吻合．对粘弹流体瞬态挤出胀大过程进行了模拟,并分析了运动粘度比和剪切速率对挤出胀大率的影响,得到的胀大率结果与理论分析和其它模拟结果基本一致．表明给出的LBM可以捕捉挤出胀大的瞬态效应．     

基于REV尺度格子Boltzmann方法的页岩气流动数值模拟

结合页岩扫描电镜图像,提出页岩气藏物理模型,采用表征单元体积（representative elementaryvolume, REV） 尺度格子Boltzmann 方法,考虑滑脱效应,模拟页岩气在页岩气藏中的流动. 模拟结果表明,页岩气主要沿着天然裂缝窜进,但在有机质和无机质中也存在缓慢的流动,且有机质中的流速要略大于无机质中的流速. 通过改变地层压力,研究地层压力对页岩气渗流特性的影响. 研究结果表明,整个流场的速度和渗透率均随着地层压力的下降而增加.     

用格子Boltzmann方法模拟Belousov-Zhabotinsky反应中的靶型波

构造了用于Belosov-Zhabotinsky反应的格子Boltzmann模型。通过对使用多组分的分布函数满足的格子Boltzmann方程,进行多重尺度Knudsen数展开,得到了模型的平衡态分布函数的各向同性解。作为算例,给出随机初始条件下反应区域内的靶型波的模拟结果,再现了Belousov-Zhabotinsky反应的经典结果。     

Numerical simulation of laminar jet-forced flow using lattice Boltzmann method

In the paper, a numerical study on symmetrical and asymmetrical laminar jet-forced flows is carried out by using a lattice Boltzmann method （LBM） with a special boundary treatment. The simulation results are in very good agreement with the available numerical prediction. It is shown that the LBM is a competitive method for the laminar jet-forced flow in terms of computational efficiency and stability.     

Parallel computation of turbulent flows using moment base lattice Boltzmann method

This paper describes parallel computing approach for simulating turbulent flows using a moment base lattice Boltzmann method. The distribution functions of the lattice Boltzmann method are expressed by corresponding moments. Choosing proper relaxation times for higher order moments, a minimum numerical dissipation is implicitly added to stabilise the method at high Reynolds numbers. Validation of the method is made by computing free decaying periodic turbulent flows and fully developed turbulent channel flows on a GPU platform. Though the present method requires additional work to calculate the higher order moments, it is shown that additional computational cost is negligible in the GPU computing. The numerical results stably obtained for the turbulent flows are in good agreement with those of a pseudo-spectral method and corresponding DNS database.     

基于格子-波尔兹曼法的流体能耗求解

格子-波尔兹曼法是近年来新兴的一种计算流体力学数值方法。随着这种方法的不断发展,人们将它用于流体的仿真、优化等不同场合。与此同时,一些与流场流速和压强相关的物理量(如能耗)的求解也成为关注的焦点。本文介绍了能耗这一流体宏观量的格子-波尔兹曼法求解及其实现。与传统的有限差分法不同,本文在求解有关的速度梯度时使用了格子-波尔兹曼-矩法,这种方法不但能够避免有限差分法在边界处失效的缺点,而且计算简单,算法局部性好,适合大规模并行计算。本文在分析其数值解精度的基础上,使用这种方法进行了以能耗极小为目标的直通道内椭圆挡块的参数优化。这些分析和算例分别定量和定性地说明了本文算法的准确性。     

On numerical diffusion of simplified lattice Boltzmann method

In this article, we analyze the numerical diffusion in the recently developed simplified lattice Boltzmann method (SLBM) and propose amending strategies towards lower numerical diffusion. It is noted that, in the original SLBM, the intermediate flow properties are utilized to evaluate the nonequilibrium distribution function, which may bring in excessive numerical diffusion. In the revised scheme, this evaluation strategy is nurtured by using the corrected flow properties to calculate the nonequilibrium distribution function. In the meantime, the numerically evaluated nonequilibrium distribution function only approximately fulfills the conservation relationship in the second order of accuracy. Although such approximation does not violate the global order of accuracy, offsetting the extra error would contribute to reducing the numerical diffusion. After implementing the proposed amending strategies, the revised SLBM (RSLBM) is validated through three numerical examples. The results indicate that RSLBM bears comparable order of accuracy as the original SLBM but shows lower numerical error on the same mesh size. And the reduced numerical error facilitates recovery of delicate flow structures. The proposed RSLBM can be flexibly implemented on nonuniform or body-fitted meshes, and in three-dimensional simulations.     

Numerical simulations of droplet formation in a cross‐junction microchannel by the lattice Boltzmann method

An immiscible liquid–liquid multiphase flow in a cross‐junction microchannel was numerically studied using the lattice Boltzmann method. An improved, immiscible lattice BGK model was proposed by introducing surface tension force based on the continuum surface force (CSF) method. Recoloring step was replaced by the anti‐diffusion scheme in the mixed region to reduce the side‐effect and control the thickness of the interface. The present method was tested by the simulation of a static bubble. Laplace's law and spurious velocities were examined. The results show that our model is more advantageous for simulations of immiscible fluids than the existing immiscible lattice BGK models. Computational results of multiphase flow in a cross‐junction microchannel were obtained and analyzed based on dimensionless numbers. It is found that the flow pattern is decided mostly by the capillary number at a small inlet flux. However, at the same capillary number, a large inlet flux will lead to much smaller droplet generation. For this case, the flow is determined by both the capillary number and the Weber number. Copyright © 2007 John Wiley & Sons, Ltd.     

New lattice Boltzmann model for the compressible Navier-Stokes equations

We have developed a fundamentally new type of simple lattice Boltzmann (LB) model for the compressible Navier-Stokes equations based on the kinetic system proposed by Sone. The model uses the kinetic equation of free-molecular type in the streaming process and modifies the distribution function to its Chapman-Enskog type at each time step. Compared with the current LB models, the proposed model is superior in the following two points: (i) there are no inherent errors associated with the Knudsen number; (ii) any flow parameters, including three transport coefficients, can be chosen freely according to our convenience. Numerical tests and error estimates confirm the above statements.     

Numerical study of the properties of the central moment lattice Boltzmann method

Central moment lattice Boltzmann method (LBM) is one of the more recent developments among the lattice kinetic schemes for computational fluid dynamics. A key element in this approach is the use of central moments to specify the collision process and forcing, and thereby naturally maintaining Galilean invariance, an important characteristic of fluid flows. When the different central moments are relaxed at different rates like in a standard multiple relaxation time (MRT) formulation based on raw moments, it is endowed with a number of desirable physical and numerical features. Because the collision operator exhibits a cascaded structure, this approach is also known as the cascaded LBM. While the cascaded LBM has been developed sometime ago, a systematic study of its numerical properties, such as the accuracy, grid convergence, and stability for well‐defined canonical problems is lacking, and the present work is intended to fulfill this need. We perform a quantitative study of the performance of the cascaded LBM for a set of benchmark problems of differing complexity, viz., Poiseuille flow, decaying Taylor–Green vortex flow, and lid‐driven cavity flow. We first establish its grid convergence and demonstrate second‐order accuracy under diffusive scaling for both the velocity field and its derivatives, that is, the components of the strain rate tensor, as well. The method is shown to quantitatively reproduce steady/unsteady analytical solutions or other numerical results with excellent accuracy. The cascaded MRT LBM based on the central moments is found to be of similar accuracy when compared with the standard MRT LBM based on the raw moments, when a detailed comparison of the flow fields are made, with both reproducing even the small scale vortical features well. Numerical experiments further demonstrate that the central moment MRT LBM results in significant stability improvements when compared with certain existing collision models at moderate additional computational cost. Copyright © 2015 John Wiley & Sons, Ltd.