用格子Boltzmann方法模拟弹性哑铃在牛顿流体中的运动

利用格子Boltzmann方法分别模拟了在同等初始条件下释放的单个带电和不带电弹性哑铃在牛顿流体中的运动．模拟所得的数值结果表明，对于单个不带电的弹性哑铃，当其被置于质心偏离管道中心线的初始位移处水平地释放后，弹性哑铃逐渐向管道的中心线漂移，当其质心漂移到达管道中心线后，弹性哑铃呈现稳定的垂直下落状态，在沉降运动中，弹性哑铃的方向保持水平方向不变．对于同等初始条件下释放的单个带电弹性哑铃，当流体的介电常数小于管道壁的介电常数时，由于静电极化的作用，弹性哑铃在其沉降过程中将受到管道壁的吸引力，该力使得弹性哑铃的运动发生显著的改变，以至于它处于最终垂直下落状态时，其质心尚未到达管道的中心线．     

液滴撞击液膜喷溅过程的LBM模拟

气液两相流动现象广泛存在于自然界和工程应用中,而液滴撞击液膜后产生喷溅的过程则是这一类问题的典型代表.采用格子Boltzmann方法（LBM）对上述过程进行了数值模拟.采用基于“单相”（single-phase）模型的LB两相流方法计算得到了3组Reynolds数（Re）和Weber数（We）组合下液滴撞击液膜后所产生的3种效应.结果表明,随着We数和Re数的不断增大,液滴撞击液膜后将产生铺展、喷溅以及喷溅并伴有小液滴脱落等不同现象,计算所得的结果与实验及理论分析的结果相吻合,表明了LBM研究气液两相流动问题的可行性.     

Burgers方程定态激波解的格子Boltzmann方法模拟

用一维５速格子Ｂｏｌｔｚｍａｎｎ方法推导一维Ｂｕｒｇｅｒｓ方程,采用Ｃｈａｐｍａｎ－Ｅｎｓｋｏｇ展开的局域平衡分布函数形式并由守恒方程确定其展开系数。最后用格子Ｂｏｌｔｚｍａｎｎ方法模拟了对流系数ｖ＝１,扩散系数Ｄ＝０．３３３３的一维Ｂｕｒｇｅｒｓ方程的激波解。计算机数值模拟结果与理论解析解精确吻合。     

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

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

致密多孔介质中气体非达西现象的LBM模拟

针对气体在致密多孔介质中低速渗流时,其渗流规律在渗流曲线的低压段表现出对达西定律线性关系的偏离,存在着非达西现象。采用格子Boltzmann方法,研究气体和多孔介质的特性对气体渗流Klinkenberg效应的影响因素。结果表明:在气体渗流曲线的低压力梯度段,随着气体黏度系数、净围压、渗透率和孔隙率的变小,渗流曲线的非线性临界点向压力梯度增大的方向移动,对达西定律线性关系的偏离更明显。说明在低渗和低压情况下Klinkenberg效应不能被忽略,气体黏度系数和孔隙率对Klinkenberg效应作用有影响;当净围压或渗透率很大时,气体渗流流量和压力梯度符合达西定律线性关系。     

二维对流扩散方程的格子BGK模拟

用D2Q4模型给出二维对流扩散方程的带修正项的BGK型格子Boltzmann方法.数值模拟与理论结果相吻合.     

微通道气体流动的格子-Boltzmann法模拟

格子-Boltzmann方法的迅速发展为计算流体力学提供了一种强有力的工具．结合“贴壁层”的理论对格子—Boltzmann方程在微槽道流动中的应用进行探索，对Kn数跨越速度滑移区和过渡区的若3个微通道中的气体流体进行了模拟，给出了通道内速度、压力、局部阻力系数的变化曲线，并和同条件下用不考虑滑移的计算结果和文献给出的理论结果进行了对比，结果表明和文献给出的理论结果吻合得很好．     

用LBM方法模拟壁面驱动粘性不可压半圆形空腔流

基于格子Boltzmann方法(LBM)数值模拟壁面驱动的粘性不可压半圆形空腔流. 采用具有二阶精度的曲线边界处理方法, 得到了不同雷诺数下的流线图、 涡线图及速度分量沿半圆形中心线的分布. 在小雷诺数的条件下, 流动状态仅由一个涡组成; 随着雷诺数的增加, 出现一个二级涡, 涡的大小与雷诺数有关. 数值结果表明, 格子Boltzmann方法简单有效, 适合处理该问题.     

基于GPGPU的LBM障碍绕流的实时模拟

基于显卡的通用计算(GPGPU)是近年来并行计算和快速绘制的热点.格子Boltzmann方法(LBM)作为流体动力学的新方法,其并行性好,常常用于基于物理的流体模拟,且具有适用于复杂边界障碍的特性,但计算较为复杂.利用GPGPU技术来加速LBM的流体计算模型,构建了基于图形处理器(GPU)的流体计算框架,实现了格子Boltzmann计算的D2Q9和D3Q15模型,并用于实时的障碍绕流模拟.     

等离子熔射射流中粉末飞行过程的数值模拟

在正六边形7-bit格子Boltzmann(LB)方法等离子射流的温度场和速度场的计算模型基础上,采用单个颗粒加速方程,建立了一个随机算法,实现了对粉末颗粒在射流场中运动过程的模拟;计算结果通过动画演示了,粉末飞行过程,结果表明在射流场一定的情况下,减小送粉口径或增加粉末颗粒直径可以提高粉末利用率。但后者会降低粉末飞行速度;选择合适的喷涂距离可以获得更高的粉末飞行速度．     

气体微尺度流动的格子Boltzmann模拟

采用反弹-镜面反射(BSR)组合边界格式的格子Boltzmann方法(LBM)对气体微尺度流动进行了模拟计算，模拟结果能够与已有的研究结果较好吻合.分析了BSR组合边界格式中比例系数的选取与Knudsen数的关系，研究结果表明，当Kn∈(002，016)时，比例系数随Knudsen数的增加呈现先减小后增大的趋势，当Knudsen数约在0055~0072范围内时，比例系数有最小值.在此基础上，对BSR组合边界格式中的比例系数进行了修正，明显改善了数值计算的结果，为微尺度气体流动的LBM模拟提供了参考.     

圆柱绕流的LBM模拟

Lattice Boltzmann Method(LBM)是一种近年来发展的一种数值方法。它具有并行效率高,边界处理简单的特点。本文采用一种能对曲线边界进行较好处理的方法,用LBM对Re=100圆柱绕流进行了计算,计算结果和经典结果一致。进一步,对柱群间复杂流场做了模拟,结果表明,此方法在处理复杂边界是有效的,并且具有较好的并行效率。     

用Lattice Boltzmann方法计算矩形柱的绕流问题

采用二维9速度Lattice Boltzmann模型, 数值模拟流体流过矩形截面柱体的绕流问题. 在边界处, 采用二阶精度的插值边界处理方法, 计算了流动的Strouhal数及柱体受到的升力和阻力系数, 给出了流场的流线和等涡线图. 使来流方向与矩形柱的长边方向平行, 计算结果表明, 改变矩形的长/宽比, 流场的Strouhal数随长/宽比呈线性变化.     

用格子Boltzmann方法计算椭圆柱绕流的阻力

用格子Boltzmann方法模拟椭圆柱绕流, 研究椭圆柱形状对阻力的影响. 对圆柱绕流问题进行了数值模拟, 阻力系数的数值计算结果与相关文献数值相符. 计算了当Re=200, 椭圆柱纵轴长度不变、 横轴长度逐渐变大时几种不同形状的椭圆柱绕流, 并用插值方法处理了曲线边界, 用动量转换法计算了曲线边界受力. 计算得到了不同形状椭圆柱绕流的流线、 涡线以及阻力系数随横轴/纵轴长度比的变化趋势. 通过分析流线和涡线的变化, 给出了阻力变化的机理.     

正型排布四圆柱绕流流动模式的格子Boltzmann模拟

利用格子Boltzmann数值方法模拟了绕正型排布的四个圆柱的二维层流流动,研究了雷诺数Re=70时,16种不同圆柱节距比下的流动模式,提出了无因次数α,β对流动模式进行分类,α为上游圆柱与下游圆柱的关联程度,β为尾涡形态的关联程度。结果表明α主要受横向节距比的影响,当α<0.4时为稳定屏蔽流,当0.4≤α<1.0时为摆动屏蔽流,当α≥1.0时为涡脱落流态;而β主要受纵向节距比的影响,当β<0.26时为单体模式,当0.26≤β<1.0时为干扰模式,当β≥1.0时为涡分离模式。本研究为多柱绕流流态的定量划分进行了探索。     

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

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

基于Lattice Boltzmann方法的方柱绕流模拟

从分子动力论出发,探讨了模拟不可压缩粘性流动的Lattice Boltzmann（LB）方法．根据粒子平衡分布函数和非平衡反弹思想,提出了高效的出口边界条件处理格式．结合区域分裂技术设计了适合集群计算机运算的LB并行算法,并用C＋＋程序语言开发了LB并行计算程序．通过对不同雷诺数时的方柱绕流模拟,验证了本文方法的有效性,同时详细分析了雷诺数对方柱绕流特性的影响．     

基于格子玻尔兹曼方法的自仿射粗糙裂隙渗流模拟

基于格子玻尔兹曼方法模拟研究了自仿射粗糙裂隙的渗流特征.通过连续随机叠加技术,生成了不同粗糙指数的自仿射粗糙裂隙.为验证格子玻尔兹曼方法,将格子玻尔兹曼结果与泊肃叶流解析解对比.通过格子玻尔兹曼方法与自仿射粗糙裂隙技术的结合,成功模拟了粗糙裂隙内部的流场分布,并分析了粗糙裂隙内部的渗流路径与水力等效隙宽的变化.结果表明,不均匀的隙宽分布对裂隙的流场分布影响较大;尽管两组裂隙内部的Re一致,但其最大流速有一定差异.     

光滑壁面明渠水流的三维格子玻耳兹曼模拟

明渠水流与泥沙运动、河床演变等问题密切相关,明渠水流特征的研究对于揭示明渠水流作用下的泥沙运动机理具有十分重要的意义.为此.应用新型数值模拟方法三维格子玻耳兹曼方法对明渠均匀流进行了研究.首先介绍了格子玻耳兹曼方法的理论基础、边界条件处理以及流动的驱动方式,接着分别模拟了明渠层流和光滑壁面紊流的近底水流运动.结果表明,得到的明渠层流流速和切应力分布与解析解吻合较好,同时较好地反映了光滑壁面明渠紊流的流速分区结构、紊动特性和阻力规律,从而说明格子玻耳兹曼方法在明渠水流模拟中具有较好的适用性,为进一步采用该方法研究明渠流近底泥沙颗粒运动机理奠定了基础.