复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

基于Shan-Chen模型的格子Boltzmann方法在微流动模拟研究中的应用

对格子Boltzmann方法的本质及Shan-Chen模型的核心机制进行了全面阐述, 并从应用实例角度对基于Shan-Chen模型的格子Boltzmann方法在微流动模拟方面的有效性、适应性进行了详细分析. 结果表明, Shan-Chen模型易于耦合微观条件下占主导作用的微观力, 拓宽了格子Boltzmann方法在微流动模拟方面的应用. 同时, Shan-Chen模型在润湿性边界条件表征方面的优势, 使得这种方法在微结构表面的滑移效应模拟方面具有很好的应用前景.     

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

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

基于格子Boltzmann方法的电流变液动力学建模及数值模拟

电流变悬浮液内部结构对外电场的快速响应发生在指定的控制空间中,在这一特指的时间和空间尺度上电流变悬浮液的物理行为特征主要为剪切速率低和流动阻尼大,即Mach和Reynolds数一般不大,可以视为微尺度流动来加以研究。针对这一流动特征,基于介观动理论的格子Boltzmann方法,建立了电流变悬浮液两相流动的离散颗粒运动模型,通过该模型进行了动力学模拟,结果表明,该模型解决了分子动力学模型难以描述的因颗粒运动造成局部流场流变特性改变的难题,以及该流场双向耦合过程中对颗粒运动的影响。     

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

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

粘弹流体轴对称流动的格子Boltzmann方法

基于BGK碰撞模型，通过在迁移方程中引入作用力项，建立了粘弹流体的轴对称格子Boltzmann模型．通过Chapman-Enskog展开，获得了准确的柱坐标下轴对称宏观流动方程．采用双分布函数对运动方程和本构方程进行迭代求解，模拟分析了粘弹流体管道流动，获得了流场中的速度和构型张量的分布，通过与解析解进行比较，验证了模型的准确性．研究了作为粘弹流体流动基准问题的收敛流动，对涡旋位置进行了定量分析，将回转长度的计算结果与有限体积法进行了比较，两种数值结果十分吻合．研究结果表明，模型能够准确表征粘弹流体的轴对称流动，具有较广阔的应用前景．     

Ginzburg—Landau方程的格子Boltzmann解

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

微尺度气体流动

了解微尺度气体流动特点是微机电系统设计和优化的基础.有关的研究可以上溯到20世纪初Knudsen的平面槽道流动质量流量的测量和Millikan的小球阻力系数的测量,实验结果揭示了稀薄气体效应即尺度效应对气体运动的重要影响.由于流动特征长度很小,微尺度气流经常处于滑流区甚至过渡领域,流动的相似参数为Knudsen数和Mach数.因此可以考虑利用相似准则,通过增大几何尺寸、减小压力的途径,解决微机电系统实验观测遇到的困难.为解决直接模拟Monte Carlo方法分析微机电系统中低速稀薄气流遇到的统计涨落困难,我们提出了信息保存法(IP),该方法能够有效克服统计散布,并已成功用于多种微尺度气流.      

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

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

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

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

LB方法应用实例——水电站水锤计算

建立了模拟水电站水锤的一维ＬＢ模型．通过实例计算与特征线法进行了比较，结果表明ＬＢ方法效果很好，将之引入水电站流场模拟是可行的，值得进一步探讨．     

流向振荡圆柱绕流的格子Boltzmann方法模拟

用一种新近发展起来的格子Boltzmann方法(LBM)在相对较小的雷诺数(Re \le 200)条件下模拟了不可压缩的流向振荡圆柱绕流问题, 考查了涡脱落模态和升阻力特性. 通过模拟, 在近尾流区发现了实验研究中已经发现的对称/反对称的涡脱落模态, 包括有些传统数值方法未发现的模态. 研究了频率锁定区域的范围及其与振幅的关系, 发现振幅越大, 发生锁定的频率区域越宽. 此外还对升阻力进行了定量意义的模拟,研究了振荡频率和振幅与升阻力的关系.      

基于热格子Boltzmann方法的自然对流数值模拟

对Eggels和Somers提出的热格子Boltzmann格式进行了改进. 在不可压缩流动的假设下,提出了一种新的温度平衡分布函数,可以克服压缩性对温度统计的影响,并且相应地修正了统计宏观温度的方法. Eggels和Somers的方法对速度和温度均采用半步长反弹格式边界条件,适合无滑移的速度边界条件.但是对温度采用该边界条件在物理本质上显得不够准确,所以在边界上对二者统一采取算法既简单又容易实现的非平衡态外推格式,同时可以与Boltzmann格式的整体二阶精度保持一致. 最后,利用改进的热格子Boltzmann方法(TLBM)模拟了Ra=10^6和Pr=0.71(空气)的方腔中的自然对流,模拟得到的流动参数与其它数值方法的结果吻合得很好,表明改进的热格子Boltzmann方法可以有效准确地模拟非等温流动.      

横向振荡圆柱绕流的格子Boltzmann方法模拟

基于格子Boltzmann方法(LBM)对不可压横向振荡圆柱绕流问题进行了数值研究. 与传统的求解宏观的N-S方程的数值方法不同, LBM求解此类问题不需要采用动网格, 而且不需要对网格进行特殊处理, 从而节约了计算成本. 结果显示, 当振荡频率增加到相应的静止圆柱绕流的自然涡脱落频率附近时, 圆柱后最新形成的集中涡距离柱体越来越近, 直到达到一个极限位置. 随后, 集中涡突然转向圆柱体另一侧脱落. 当振荡频率接近于静止圆柱的自然涡脱落频率时, 发生频率同步的现象. 随着振荡频率远离自然涡脱落频率, 同步现象消失. 在几种次谐振荡和超谐振荡下, 尾流区的涡脱落频率仍为相应的静止圆柱绕流的自然涡脱落频率.      

非牛顿流体流动的格子Boltzmann方法研究进展

格子玻尔兹曼方法（lattice Boltzmann method,LBM）能够直接计算局部剪切速率并可以达到二次精度,因此在非牛顿流动数值模拟中展现出一定优势。尽管已证实LBM 对于非牛顿流动的适用性,但是LBM 需要通过即时调节BGK（Bhatnagar-Gross-Krook）碰撞项中的松弛时间来实时反映黏度改变,当松弛时间接近1/2 时,迭代会出现数值不稳定现象。该文对LBM 在非牛顿流体研究中的进展进行了总结,介绍了增加数值稳定性的方法并对结果的精度进行了比较,在此基础上对LBM 在非牛顿研究中的进一步发展进行了展望。     

建筑物周围流场格子Boltzmann仿真

基于格子BGK方程的三维15点格子Boltz-mann模型,对建筑物周围流场进行了数值仿真,对其流场的划分进行了分析,并对不同雷诺数下建筑物周围的流场进行了数值模拟计算,得到了比较满意的结果,为环境质量的评估提供了新的途径.     

基于波尔兹曼方法研究幂律流体的椭圆柱绕流

基于格子波尔兹曼方法(lattice Boltzmann method) 和幂律流体本构方程, 建立二维流动模型, 将充分发展的速度分布与理论解进行对比, 吻合良好. 对幂律流体的圆柱绕流进行模拟, 采用了反弹格式的无滑移圆柱边界, 并使用应力积分法计算阻力系数, 分析了稳态圆柱绕流时, 阻力系数随幂律指数n 以及雷诺数Re 的变化规律. 分析了椭圆横轴/纵轴长度比和幂律指数n, 对压力系数Cp 和黏度系数Cv 的影响. 得到的变化规律与有限元方法规律一致, 验证了格子波尔兹曼模型的可行性.     

格子Boltzmann方法模拟多孔介质惯性流的边界条件改进

格子Boltzmann方法可以有效地模拟水动力学问题,边界处理方法的选择对于可靠的模拟计算至关重要.本文基于多松弛时间格子Boltzmann模型开展了不同边界条件下,周期对称性结构和不规则结构中流体流动模拟,阐述了不同边界条件的精度和适用范围. 此外,引入一种混合式边界处理方法来模拟多孔介质惯性流, 结果表明:对于周期性对称结构流动模拟,体力格式边界条件和压力边界处理方法是等效的,两者都能精确地捕捉流体流动特点; 而对于非周期性不规则结构,两种边界处理方法并不等价,体力格式边界条件只适用于周期性结构;由于广义化周期性边界条件忽略了垂直主流方向上流体与固体格点的碰撞作用,同样不适合处理不规则模型;体力-压力混合式边界格式能够用来模拟周期性或非周期性结构流体流动,在模拟多孔介质流体惯性流时,比压力边界条件有更大的应用优势,可以获得更大的雷诺数且能保证计算的准确性.