两相格子Boltzmann方法的液滴冲击流动液膜数值模拟

基于最新的两相格子Boltzmann模型(Lee模型)在数值方法上进行了改进,进一步提高了模型在大密度比下的数值稳定性。基于改进后的模型首先研究了二维液滴冲击静止液膜,计算得到的溅射根部铺展半径与实验结论一致吻合,从而验证了改进后模型的正确性。进一步数值模拟了液滴冲击流动液膜的溅射过程,考虑了不同的液膜液滴速度比和液膜相对厚度对流动过程的影响。计算结果与冲击静止液膜的计算结果和实验结论进行了对比,总结了溅射过程中溅射铺展半径和溅射高度随时间的变化规律,论述了流动现象产生的内在机理。     

基于格子Boltzmann方法的空泡溃灭过程模拟

针对自由边界的三维空泡溃灭问题,以格子Boltzmann方法为流场核心求解器,结合伪势多相流模型与Carnahan-Starling(C-S)状态方程,开展了定压边界条件下,不同温度、初始溃灭半径的三维空泡溃灭过程的模拟研究．通过模拟给出了空泡溃灭过程中流场的两相界面图,并定量分析了三维空泡溃灭时间随着半径变化的关系．结果表明：在定压强边界下,空泡的初始半径与总溃灭时间呈线性相关,且不同温度、不同初始半径的空泡溃灭具有相似的过程;相同温度条件下,空泡溃灭速度随时间不断增加,空泡初始半径越大,空泡溃灭速率最大值出现的时间越晚;对于相同半径条件下,温度比越大,空泡溃灭最大速率出现的时间越晚．     

Burgers方程的格子Boltzmann方法模拟

用一维格子Boltzmann方法构造O(ε^4)的Burgers方程模型．格子Boltzmann方法的数值模拟结果与具有特定边界条件的Burgers方程的解析解精确吻合．     

用格子Boltzmann方法模拟MKDV方程的行波解

用精确到0(ε^4)的5速格子Boltzmann模型较系统地对MKDV方程：ut θu^2ux βuxxx=0的行波解进行了模拟，并将模拟的数值结果与相应的理论结果进行了比较，模拟值与理论值吻合很好。     

非挥发性溶质溶液两相问题的格子Boltzmann方法

溶液的气液两相体系广泛存在于工程问题中。由于其多组分和多相态以及组分之间相互作用的复杂性,用传统数值方法进行研究较为困难。该文提出了一种适用于非挥发性物质稀溶液的气液两相问题的格子Boltzmann方法。在Shan-Chen多组分格子Boltzmann模型的基础上,引入溶质与溶剂粒子间的作用力,建立了合理的作用力模型,完成了对溶液气液两相体系的数值模拟。模拟出了溶质的非挥发性以及表面非活性物质溶液的表面性质。表面吸附及表面张力的计算结果与相应理论公式符合很好。     

基于CUDA的格子Boltzmann数值模拟加速实现

针对近年来利用CUDA技术在个人计算机显卡的GPU上实现LBM并行加速计算的研究越来越多,但对在GPU中使用不同GPU存储器进行计算的具体实现算法以及其对计算性能的影响分析研究甚少,文章实现了在GPU中使用不同存储器进行IBM并行计算,给出了具体的实现算法,并以平面Poiseuille流为算例,在普通个人计算机上,分别使用NVIDIA GeForce GTS 450 GPU和Intel Core i5-760 4核CPU进行计算.结果表明,两者计算结果吻合得很好,最高获得了约107倍的加速比,验证了在GPU上进行LBM并行计算的可行性以及加速性能,为在低成本的个人计算机上高效率地解决计算流体力学中的复杂计算问题提供了一种非常有效的途径.     

带自由面流体运动的单相格子Boltzmann方法模拟

采用Thurey提出的单相格子Boltzmann方法模拟了带自由面流体的运动.它在标准的单相模型的基础上将每个格子标记为气体、液体或自由面.格子的类型通过计算格子间的质量通量来重新初始化.与多相模型相比,它不需要求解气相格子的运动,提高了计算效率.针对此方法,给出了一种来流和出流边界条件的提法.最后,计算了来流被竖直圆柱阻挡,破碎成两股后又接合起来运动的动态过程.结果表明,此方法能描述来流经圆柱后高度降低,并形成涡等现象.     

热格子Boltzmann方法分析及应用

格子Boltzmann方法（lattice Boltzmann method,LBM）是一种基于气体动理论的介观计算方法,其物理背景清晰、边界处理简单,已成功应用于等温（或无热）流动中.简要介绍现有的几种热格子Boltzmann模型,并运用几种热格子模型求解热Couette流、方腔自然对流等典型算例,对比不同热格子模型的数值稳定性、准确性、模型的计算效率等.将两种热格子模型用于多孔介质内的流动与传热问题中,对比热格子模型在处理复杂结构时的数值特性.     

用格子Boltzmann方法模拟正弦曲线底边方腔内流动

采用格子Boltzmann方法数值模拟正弦曲线底边方腔内流动, 分析了格子Boltzmann方法处理曲线边界的特性.在曲线边界的处理中采用二阶精度的曲线边界处理方法, 将反弹格式和内插法相结合, 计算方法可靠、准确且易于执行.计算了流线图、等涡线图和涡心位置, 并分析了流场随Re数的变化.     

沸腾过程的格子Boltzmann方法模拟

提出了一种新的描述气液相变过程的格子Boltzmann理论模型,并利用该模型对液化和蒸发过程进行了模拟.对液化过程的计算结果表明,气相的体积分数、液相和气相的密度与理论解的最大相对误差分别为1.18%、0.23%、0.32%,而蒸发过程的计算结果相应的最大相对误差分别为0.08%、0.48%、0.20%,较前人研究的结果更接近理论解.为进一步验证该模型模拟复杂相变问题的可行性,利用该模型对池内沸腾和垂直管道中的流动沸腾现象进行了模拟,观察到气泡随时间变化不断长大,气泡受浮升力的影响脱离壁面,以及气泡在液相中上升时发生碰撞、聚合等现象.     

Simulation of wind gust structure in the atmospheric boundary layer with Lattice Boltzmann Method

Cold fronts occur in northern East Asia during winter and spring.After cold frontal passage,airflow is downward and accompa-nying strong winds fluctuate significantly;this is termed wind gusts.Analysis of observation data shows that wind gust structure has coherent characteristics.This is important for entrainment of spring dust storms into the upper boundary layer,where they are transported great distances.The Lattice Boltzmann Method(LBM) is a computational fluid technique based on the Boltzmann transport equation.The LBM has been used to study complex motion such as turbulence,because it describes motion at the micro level.In this paper,Large eddy simulation is introduced in the LBM,enabling simulation of turbulent flow in the atmospheric boundary layer.The formation and development of wind gusts are simulated,and a coherent structure with a combination of wave and vortex is obtained.This explains the mechanism of soil erosion and sand entrainment by the coherent structure of wind gusts.     

用Lattice Boltzmann方法模拟方柱绕流

用Lattice Boltzmann方法模拟了两种情况下的方柱绕流：方柱位于流场的边壁和方柱位于流场的中央。在第1情况中，分析了流场的流线、速度矢量分布以及方柱后面绕流形成的回流区长度随Reynolds数的变化曲线；在第2种情况中，除了给出稳定时流场的流线、速度矢量分布和压力分布以外，还模拟了Re＝150时方柱后面形成的周期旋涡现象。两种情况下的计算结果与其他文献中的计算数值和试验数据是一致的。     

方法对于大粘度差两相流的模拟

为了模拟具有大粘度差的两相流,对Boltzmann(LB)方法中的伪势多组分模型进行了改进.伪势模型将粒子作用力转化为速度形式,再将其引入离散的LB方程.改进模型把力的添加形式变为将力离散处理后,引入离散的LB方程中,以此提高数值计算的稳定性.该模型能够比原伪势模型得到更小的虚假速度和更大的粘度差.通过对沥青-水的大粘度差两相流体和粘性指进现象的模拟,验证了改进模型的正确性和有效性.该模型具有较好的适用性,它能模拟的两相流体最大粘度差可达到4 510倍.     

三维低雷诺数后台阶流动的Lattice Boltzmann Method模拟

Lattice Boltzmann Method (LBM) 是近年来兴起的新型流体计算方法,该方法已成功应用模拟各工程中的流体流动.本文利用该方法研究在低雷诺数的不同雷诺数下三维后台阶流动的计算,并与前人的研究结果进行了比较,计算表明:与传统的计算流体力学方法相比,该方法具有极易处理复杂的边界条件和卓越的并行计算性能.已证明该方法在模拟恒温、不可压和低雷诺数条件下的流体流动时,有非常好的计算精度.     

波形板内流场的格子波尔兹曼模拟

利用格子波尔兹曼方法模拟了波形板复杂流道的流场,计算了流场的状态。采用的是十三点格子波尔兹曼模型,边界处理成无滑移、反弹边界。结果表明,该模型与其他数值计算方法相比,能更方便地处理边界,因此能有效、精确和稳定地模拟复杂流动现象。     

可压Euler方程的格子Boltzmann方法:二维情形

根据微观和宏观之间的质量、动量和能量守恒原理，构建了一个可用于模拟可压流体的D2V19格子Boltzmann模型，利用Chapman—Enskog分析可从该模型推导出宏观流体力学方程-Euler方程．利用该模型对三种典型激波管进行数值模拟并与Riemann解析解进行比较和分析，发现该模型能较好的捕捉激波，精度令人满意，证明了该模型对可压流体的适用性．该模型适用于气体多方指数为任意值的系统．     

双重网格的Lattice Boltzmann方法

针对均匀网格的Lattice Boltzmann方法（LBM）计算效率较低的不足，引入有限差分方法（FDM）中双重网格算法的基本思想，结合LBM在连续流场上变量也连续分布的特点，推导出双重网格的LBM。这种算法充分利用粗、细两种网格各自不同的优点，加快了收敛速度，提高了流场的计算效率。通过模拟二维Poiseulle流动和后台阶流动两个算例，证明双重网格算法的计算时间明显地少于均匀网格算法的计算时间，计算结果正确可信。     

假塑性流体液滴撞击壁面上的铺展的格子Boltzmann模拟

流体液滴在壁面上的铺展与沉积在工业应用中存在广泛应用,研究非牛顿流体特性对液滴铺展的影响具有重要的实际意义。基于相场法建立了非牛顿幂律流体与牛顿流体的两相流格子Boltzmann模型,引入包括接触角影响的边界条件,模拟了假塑性流体在幂律指数（n）0.5~1.0以及韦伯数（We）5~45时液滴撞击壁面的铺展过程。结果表明：相对于牛顿流体,幂律流体的非牛顿特性会抑制液滴的铺展,且n越小,越有利于液滴沉积。此外,韦伯数越大,液滴越快地达到稳定。     

基于格子波尔兹曼的图像修复模型研究

在研究格子波尔兹曼图像修复模型的基础上,提出了一种改进的模型,该模型通过引入调节函数,使图像在梯度变化较大的区域减弱其扩散强度,在梯度变化较小的区域增强其扩散强度,同时引入调节因子来控制其作用强度。实验结果证明改进后的模型具有更好的修复效果。