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

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

基于格子Boltzmann方法的两液滴相溶模拟

原多组分Shan-Chen模型在模拟气液流动时受到密度比的限制,经过测试,该模型在模拟气液流动时最大密度比为50左右。本文提出了利用改进后的多组分Shan-Chen模型模拟密度比为175的两液滴相溶,模拟结果与实验结果相符。     

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

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

基于格子Boltzmann方法的储层岩石油水两相分离数值模拟

针对岩石物理试验中ffl现的孔隙流体(油水)两相分离现象,应用格子Boltzmann(LB)]方法中的两相不相溶流体的伪势模型,对油水界面动力学行为进行微观数值模拟,分析多孔介质中两相流动的微观特征,并从理论上给出两相不相溶流体界面张力因子Gf,值的确定方法.模拟由于表面张力造成的油水两相分离现象,在此基础上研究润湿性对真实储层岩心孔隙流体两相分离的影响,并实现全程动态可视化.研究表明,用LB方法进行储层岩石油水两相分离简便易行、形象直观,是研究流体分离规律和特点的重要评价方法.     

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

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

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

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

Burgers方程的格子Boltzmann方法模拟

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

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

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

二维Selkov反应扩散模型的格子Boltzmann方法模拟

建立了二维Selkov反应扩散模型的格子Boltzmann方法,进行计算机模拟演化,模拟出Selkov反应的均匀定态、极限环、图灵斑图、环状波等空间结构,表明在模拟求解流速为零情况下的反应扩散方程时,这种方法是现实可靠的．     

二维反应扩散方程的格子Boltzmann方法模拟

建立了二维反应的散方程五速四方格子模型的格子Boltzann方法，进行计算机模拟演化，分析模拟结果，并与SChen利用七速六方格子模型得到的结果进行比较，得出在模拟求解流速为零情况下的反应扩散方程时，这两种格子模型是一致的结论。     

基于格子Boltzmann方法的量子等离子体离子声波的数值模拟

用格子Boltzmann方法(LBM)研究由KPB(Kadomtsev-Petviashvili-Burgers)方程控制的量子等离子体中离子声波的数值模拟问题, 得到了KPB方程中离子声波的传播演化结果.     

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

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

基于格子Boltzmann方法的量子等离子体离子声波的数值模拟

用格子Boltzmann方法(LBM)研究由KPB(Kadomtsev-Petviashvili-Burgers)方程控制的量子等离子体中离子声波的数值模拟问题, 得到了KPB方程中离子声波的传播演化结果.     

用格子Boltzmann方法研究交通流BML模型

用格子玻尔兹曼研究二维决定论文通流BML模型，得出模型满足的玻尔兹曼方程，并对其进行了计算机模拟，获得了与用CA方法研究该模型相类似的结论，提供了格子玻尔兹曼方法用于研究交通流的可行性。     

用格子Boltzmann方法模拟椭圆柱体在牛顿流体中的二维运动

用格子Boltzmann方法建立了椭圆柱体的二维运动模型，利用所建立的模型对椭圆柱体在牛顿流体中的二维运动进行了数值模拟，并通过对照圆柱体存相同条件下的运动，深入讨论了椭圆柱体二维运动的特征，得到了一些有意义的结果，所建立的模型，可以推广到对其他形状柱体二维运动的模拟．     

热格子Boltzmann方法分析及应用

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

格子Boltzmann方法模拟磁场作用下的融化传热过程

采用格子Boltzmann方法对磁场作用下方腔内相变材料的融化传热及液相流体流动过程进行数值模拟,重点分析磁场倾角对融化相变传热过程的影响,为实现相变传热的主动控制技术提供参考.研究表明磁场可以控制融化相变传热过程,当磁场倾角大于180°时,磁场加速度方向指向热壁面方向,磁场强化换热,热壁面平均Nu_h数增加,融化速率加快;当磁场倾角小于90°、磁场加速度方向指向冷壁面侧时,磁场弱化换热,热壁面平均Nu_h数下降,融化速率降低;磁场倾角为231°时,磁场强化传热效果最佳.     

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

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

格子Boltzmann方法模拟狭窄颈动脉修复前后的流场

采用格子Boltzmann方法模拟真实人体颈动脉血管狭窄前后脉动流场的速度、压强以及壁面剪切应力分布,给出脉动流影响下流体分离区可能出现的位置以及壁面剪切应力的分布情况,分析脉动血液流在二维不对称动脉分叉血管中容易发生动脉粥样硬化的动力学机制,并对模拟结果进行定性分析,与已有结论进行对比,为血管壁病变和动脉粥样硬化形成机制提供有用信息.     

基于格子Boltzmann方程的大涡模拟对湍流时空关联性的研究

将格子Boltzmann方程和大涡模拟(LBE-LES)相结合,提出适应于格子Boltzmann方法(LBM)的涡黏性亚格子尺度模型,开展均匀各向同性湍流时空关联性的研究.采用D3Q19格式计算湍流的三维能谱、湍动能耗散率和其它高阶统计量,与实验和直接数值模拟结果的比较表明,该模型比传统涡黏模型有明显改进.考察了不同亚格子模型预测湍流频率波数能量谱的能力,结果表明,尺度涡产生的横扫作用是造成小尺度涡时间去关联的主要因素,不同波数的频率能量谱之间有一定的相似性,横扫速度是描述湍流频率波数能量谱的特征量.