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

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

浸没边界-简化热格子Boltzmann方法研究及其应用

针对流固耦合传热问题,本文提出了一种基于浸没边界-简化热格子玻尔兹曼方法(immersed boundary method-simplified thermal lattice Boltzmann method,IB-STLBM)的耦合模型.不同于传统的格子玻尔兹曼方法使用分布函数演化流场和温度场,简化热格子玻尔兹曼方法(simplified thermal lattice Boltzmann method,STLBM)的演化过程不需要依赖分布函数,只涉及平衡态分布函数和非平衡态分布函数,能够直接演化宏观量,极大减小了计算过程中所占用的虚拟内存,简化了边界条件的实现方式,同时具有较高的稳定性.传统的浸没边界法对流场的计算采用欧拉网格,对固体边界采用拉格朗日网格,认为固体边界是对流场产生某种体积力.在应用浸没边界法时,汲取介观的思想,把固体的介入看作是对流场的干扰,打破了固体附近流体介观微团颗粒原始的平衡状态,这种干扰可以看作是在耦合边界上产生的一个非平衡项,可用非平衡态分布函数来表示.基于此,在模型中浸没边界法与简化热格子玻尔兹曼方法更紧密联系在一起,更大程度发挥二者的优点,整个计算过程更加简单直观,符合物理特性.通过对热圆柱绕流和内含热颗粒的封闭方腔自然对流问题的模拟以及对其结果的分析,验证了该算法在求解流固耦合传热问题的有效性和可行性.      

基于格子Boltzmann方法的挤出胀大的数值模拟

将单相格子Boltzmann方法(lattice Boltzmann method, LBM)引入到粘弹流体的瞬态挤出胀大的数值模拟中,建立了基于双分布函数的自由面粘弹性流动格子Boltzmann模型．分析得到的流道中流动速度分布和构型张量结果与理论解十分吻合．对粘弹流体瞬态挤出胀大过程进行了模拟,并分析了运动粘度比和剪切速率对挤出胀大率的影响,得到的胀大率结果与理论分析和其它模拟结果基本一致．表明给出的LBM可以捕捉挤出胀大的瞬态效应．     

热肋片宽度对复杂方腔自然对流影响的数值分析

采用笔者等发展的二阶全展开Euler-Taylor-Galerkin分裂步有限元方法,经过对算法精度验证,在确定不同网格密度对计算结果影响的基础上,确定出计算所需网格密度的必要条件,并在指定的网格密度条件下,分析了热肋片宽度对封闭方腔内典型自然对流流动的影响,流动对应的Pr=0.71,Ra=104。计算结果表明,随热肋...     

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

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

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

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

悬浮颗粒运动的格子Boltzmann数值模拟

将固体颗粒的牛顿力学和格子Boltzmann方法相结合,研究不规则形状悬浮颗粒在流场中的运动。通过受力分析,精确求得其所受合力、合力矩、合力作用中心等。提出了跟随颗粒运动的动网格计算域技术和模拟悬浮颗粒转动运动的局部数组方法及Euler-Lagrange两套坐标技术。通过对椭圆颗粒运动的数值模拟和对照他人对矩形颗粒的研究,分析了其复杂运动规律,并提供了合理的物理解释。结果表明:运用格子Boltzmann方法和上述特殊技术可以得到与有限元方法相同的模拟精度,且具有计算速度快、对复杂形状边界处理方便灵活、程序简单及特别适合大规模并行计算等优点。     

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

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

不同属性肋片对复杂方腔自然对流影响的数值分析

采用二阶全展开Euler-Taylor-Galerkin分裂步有限元方法,在指定的网格密度条件下,在流动对应的普朗特数取为0.71,雷诺数取为104的情况下,数值分析了热肋、冷肋、上绝热肋、下绝热肋等四种不同属性肋片对封闭方腔内典型自然对流流动的影响.计算结果表明,肋片的存在对封闭方腔内的自然对流及相应的传热效率具有较强的影响,对流流动结构以及平均Nusselt数随肋片的属性发生较大的改变.     

并列圆柱绕流的格子Boltzmann数值模拟

采用非均匀不可压格子Boltzmann模型对低雷诺数下并列圆柱绕流进行了数值模拟,给出了数值计算结果,分析了间距g对圆柱尾流及升力、阻力的影响,并在此基础上得到了4种尾迹模式.此外,研究了流场的初始扰动对流动分岔现象的影响,发现在适当的扰动下可以很快得到同步同相的尾流.对Re=160和200下圆柱的升、阻力进行了对比,结果表明升力和阻力受间距g的影响大于雷诺数.     

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

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

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

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

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

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

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

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

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

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

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

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

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

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

THE NUMERICAL STUDY OF LATTICE BOLTZMANN METHOD BASED ON DIFFERENT GRID STRUCTURE

传统的格子波尔兹曼方法(lattice-Boltzmann method, LBM)通常基于标准均匀网格, 这主要取决于速度的空 间离散格式.均匀网格结构的特点, 使LBM在处理具有复杂边界的问题时遇到较大的困难, 从而限制了它的应用.另外, 对于较为复杂的流动, 其流场存在流动变化剧烈和平缓的区域, 在流动变化剧烈的区域, 往往需要足够的网格点才能更好地捕捉到流场信息, 而均匀网格会使得网格数量过多, 这会增加计算量, 但网格数量过少又无法获得必要的流场信息, 使LBM的计算效率降低.为了解决上述问题, 用不同的网格结构, 以顶盖驱动的腔体内流、柱体绕流和翼型绕流为例, 探讨了提高LBM算法的计算效率和适用性问题.     

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

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