Burgers方程的格子Boltzmann方法模拟

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

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

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

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

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

用于求解Poisson方程的格子Boltzmann模型

提出一个求解Poisson方程的格子Boltzmann模型.通过使用Chapman-Enskog展开和多尺度展开得到了在不同时间尺度下的系列偏微分方程及平衡态分布函数和具有三阶截断的误差修正Poisson方程.用该模型计算Kolmogorov流和Green-Taylor涡流,并与解析解进行比较,计算结果表明,数值结果与经典解析结果基本相符.     

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

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

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

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

完全气体格子Boltzmann热模型

为建立一种具有任意比热比的完全气体多速度格子 Boltzm ann热模型 ,引入粒子的势能来调整压能与热力学能的关系 ;利用 Chapm an- Enskog方法从 BGK型的格子Boltzmann方程推导出了 Navier- Stokes方程和能量方程。对一维正弦波形式的能量衰减过程进行了模拟 ,测得的热扩散率与理论预测值相吻合。还模拟了绕加热平板的二维强制热对流问题 ,结果合理     

基于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方法的两液滴相溶模拟

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

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

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

明渠流中分形絮团破裂的格子Boltzmann方法模拟

应用格子Boltzmann方法模拟了黏性泥沙三维分形絮团在明渠水流中的沉降、破裂过程,揭示了絮团在水流剪切作用下的破裂机理．结果表明,水流流速较小时,絮团将不发生破裂而沉积到海床;在水流流速较大的情况下,强度较小的絮团将在水流剪切作用下发生破裂,破裂后的单个泥沙颗粒或小絮团随水流一起运动,不会沉降到床面,而强度较大的絮团在经受近底剪切变形后仍可沉降到床面．絮团沉降破裂过程的格子Boltzmann模拟结果与已有的直观认识一致．     

反应扩散问题的格子波尔兹曼模型

本文利用格子波尔兹曼方程、Chapman—Enskog展开及多尺度技术,建立了二维反应扩散问题的格子波尔兹曼模型,给出了可用于求解反应扩散系统的LBM方法。     

二维直角管道中高Reynolds流动的格子Boltzmann模拟

运用9 Bit的格子Boltzamnn方法对二维直角管道中高Reynolds流动进行数值模拟, 得到了在不同雷诺数下不可压缩流动涡的涡心的位置分布及其演化斑图.     

用于Beltrami流的格子Boltzmann算法

给出Beltrami流动中流体质点位移的格子Boltzmann算法. 引入Beltrami涡流的概念, 并给出涡量和速度的结构. 应用格子Boltzmann算法模拟Beltrami流的Arnold-Beltrami-Chidress(ABC)解.     

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

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

圆柱绕流的LBM模拟

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

Adaptive Lattice Boltzmann Model for Compressible Flows

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｔｈｅＬａｔｔｉｃｅＢｏｌｔｚｍａｎｎ(ＬＢ)ｍｅｔｈｏｄｈａｓａｔｔｒａｃｔｅｄｍｕｃｈａｔｔｅｎｔｉｏｎａｓａｎｏｖｅｌａｌｔｅｒｎａｔｉｖｅｔｏｔｒａｄｉｔｉｏｎａｌｍｅｔｈｏｄｓｆｏｒｎｕｍｅｒｉｃａｌｌｙｓｏｌｖｉｎｇｔｈｅＮＳｅｑｕａｔｉｏｎｓ[1].ＴｈｉｓｍｅｔｈｏｄｏｒｉｇｉｎａｔｅｄｆｒｏｍａＢｏｏｌｅａｎｍｏｄｅｌｋｎｏｗｎａｓＬａｔｔｉｃｅＧａｓＡｕｔｏｍａｔａ(ＬＧＡ)[2].ＴｈｅｓｔａｎｄａｒｄＬＧＡｍｏｄｅｌｉｍｐｏｓｅｓ,ｆ…     

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

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

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

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