用LBM计算沟槽形状边界Couette流

采用动量转换方法,应用格子Boltzmann模型计算流体对边界的作用力,用二维六角形的格子Boltzmann模型计算沟槽形状边界的Coutte流动,得出沟槽形状边界的阻力。计算结果显示,沟槽形状边界的阻力系数随Re数的增大而减小,其变化规律与平板边界的情况一致。     

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

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

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

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

圆柱绕流的LBM模拟

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

用13速六方格子BGK模型模拟空腔流

采用１３速六方格子ＢＧＫ模型人腔粘性流，模拟了空腔边界温度高于空腔内部温度时，腔内流体热的输运过程及腔中流体达到稳定后的流场速度、密度及温度的分布情况。     

用格子Boltzmann方法计算剪切流的方柱绕流问题

用格子Boltzmann方法计算来流为水平剪切流的方柱绕流问题, 得到了在不同速度梯度条件下方柱绕流的流线和等涡线图. 发现在圆柱尾部形成两排涡, 当来流速度梯度较大时, 两排涡有很大的不同. 计算结果表明, 用格子Boltzmann方法计算剪切流的方柱绕流问题是可行的, 计算结果与理论分析相符.     

双板驱动矩形空腔STOKES流动的数值模拟

用边界元方法对A =1.5 ,s为任意值时的双板驱动矩形空腔Stokes流动进行数值模拟 .研究了矩形空腔内涡旋生成和发展的规律 .数值计算结果表明 :在s的绝对值不太大的范围内 ,伴随s的增加 ,矩形空腔内涡旋的生成、发展过程可由s的 5个临界点分成 6个阶段 .第一个临界点 ,s =- 8.5 5 ,1个涡旋变成 2个涡旋 .第二个临界点 ,s =- 0 .118,2个涡旋变成 1个涡旋 .第三个临界点 ,s =- 0 .0 0 0 1,1个涡旋变成 3个涡旋 .第四个临界点 ,s =0 .0 5 ,3个涡旋变成 2个涡旋 .第五个临界点 ,s=15 5 ,2个涡旋变成 3个涡旋 .s =- 1时的数值计算结果与理论分析结果以及实验结果比较完全一致 ,令人非常满意 .特别是 ,本文方法适用于不同高宽比和任意速率比的矩形空腔驱动流动问题 ,因此不但具有方法上的理论意义 ,而且数值计算结果还具有工程应用上的实际意义 .     

基于LBM和FVM的方柱绕流特性对比分析

基于粘性流体理论,本文分别采用格子Boltzmann方法（LBM）和有限体积法（FVM）建立了粘性流场中方柱绕流模型,探究LBM在非光滑曲面钝体绕流方面的应用,并结合FVM进行对比分析。在FVM模型中,采用局部加密的方法对钝体边界进行处理,而在LBM模型中,除了传统的Half-way边界处理方法,还结合了拐角边界处理方法。为获得较好的可对比数据,根据已发表文献中的理论及UDF编译码技术分别对两模型的进出口边界条件进行了讨论和设置。对比分析了两模型下的速度云图以及获得的升、阻力系数,Strouhal数。结果发现方柱上游压力不受涡脱落影响,雷诺数对其影响也较小;两种方法下的速度、无量纲参数吻合较好,但两者最适进出口边界不同,且相同条件下,LBM比FVM数值模拟能更快达到稳定状态。 关键词 方柱绕流 格子Boltzmann方法 有限体积法 对比分析     

低Reynolds数下单排方柱绕流尾迹涡型的LBM模拟

针对工程中遇到的周期柱体绕流问题,采用格子Boltzmann方法(LBM)对低Reynolds数下单排周期排列方柱的定常绕流进行数值模拟,详细分析了柱体后面的尾迹涡型随着Reynolds数的变化情况,并通过对不同Reynolds数下尾迹涡型结构的分析,得到了二射流、三射流和四射流发生合并现象的临界Reynolds数,其中二射流出现合并的临界Reynolds数与文献的数值模拟结果一致,三射流和四射流出现合并的临界Reynolds数则偏大,但更接近实验结果。     

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

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

气体微尺度流动的格子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时为涡分离模式。本研究为多柱绕流流态的定量划分进行了探索。     

用于Beltrami流的格子Boltzmann算法

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

基于格子Boltzmann方法的液力

研究了基于格子Boltzmann方法(LBM)的液力变矩器导轮内流场数值计算理论与方法.首先,提出了LBM中处理旋转周期性边界条件的方法.然后,分析了LBM中各项参数对于计算稳定性以及计算效率之间的影响.为了保证粒子迁移与碰撞计算的稳定进行,必须合理选定弛豫时间τ,从而避免在计算过程中平衡态分布函数feq出现负值的情形.此外,在LBM中运用大涡模拟(LES)可以降低计算稳定性对于弛豫时间τ选取的限制,在一定程度上提高计算效率.最后,在开源代码Palabos的基础上进行功能拓展,实现导轮内流场的仿真计算,得到了导轮尾迹区域瞬时非定常流动分布特征.结果表明,LBM与传统计算流体力学(CFD)方法相比,液流对导轮作用力的计算结果在数值上较为接近.然而,使用LBM可以获得详细流场形成过程的时间历程以及局部复杂的流动细节.     

Lattice-Boltzmann方腔模型的CUDA加速实现

对Lattice Boltzmann方法(LBM)在CUDA下的建模和算法进行了研究,使得该方法在GPU下的计算速度得到提升,大大缩短了计算过程的时间消耗。利用非平衡外推边界条件处理,以LBM方法模拟了D2Q9模型的方腔顶盖驱动流动,采用全局内存和纹理内存存储数据,将模型中9个分布函数存储为二维网格,每个网格分配一个线程,每个线程块包括256个线程,多条线程并行计算。在普通个人计算机上,采用NVIDIA GeForce 9600 GT显卡和CUDA,实现了LBM模拟方腔流动,将计算速度提高到CPU的50倍。     

球面上Turing斑图的格子Boltzmann模拟

利用球坐标系下的格子Boltzmann模型, 求解球表面的Gierer-Meinhardt方程, 得到了球面上的Turing斑图. 与经典差分方法相比, 该模型可用于球表面Turing斑图的数值模拟.     

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

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