非均匀格子Boltzmann方法模拟方柱绕流

应用非均匀格子Boltzmann方法对方柱绕流的三种情况进行详细数值模拟，在第1种情况中，方柱位于流场中央，模拟了卡门涡街现象，给出了斯特鲁哈数随雷诺数变化曲线；在第2种情况中，方柱位于流场壁面，分析了雷诺数对方柱后回流区的影响；在第3种情况中，两方柱并列在流场中央，考察了方柱间距对流场的影响。     

波形板内流场的格子波尔兹曼模拟

利用格子波尔兹曼方法模拟了波形板复杂流道的流场,计算了流场的状态。采用的是十三点格子波尔兹曼模型,边界处理成无滑移、反弹边界。结果表明,该模型与其他数值计算方法相比,能更方便地处理边界,因此能有效、精确和稳定地模拟复杂流动现象。     

正弦形沙丘背风坡回流区特性研究

利用开源计算流体力学软件OpenFOAM的标准k-ε湍流模型,对沙丘表面流场进行数值模拟,获得了距沙丘表面不同高度的流场速度分布曲线,考察了沙丘背风坡回流区特性受沙丘高度和宽度的影响,分析了回流区形成的原因。结果表明:沙丘高度较低时,在背风坡贴近沙丘表面的区域气流速度先降低后增加;随着沙丘高度的增加,背风坡气流速度降低更为迅速,气流不稳定,出现回流区。研究还发现:在沙丘宽度或高度一定的情况下,随着宽高比的增加,回流区尺度减小;回流区尺度不仅取决于宽高比,还取决于沙丘宽度与高度的数值。沙丘背风坡回流区萌生时,贴近沙丘表面出现一系列小涡,减小沙丘宽高比,涡的尺度增大、强度增加,演变为回流区。     

不同形状建筑物周围风环境的研究

运用k-ε双方程湍流模型数值模拟了几种不同形状的建筑物和一组建筑物群周围气流绕流的三维流场．结果表明，当来流的风速相同时，不同形状的建筑物的气流绕流流场是不同的，建筑物的几何结构是影响其周围流场分布的主要因素之一；建筑物群周围的绕流过程比单个建筑物更加复杂，除了具有单个建筑物绕流过程中所产生的流动分离和回流等复杂特征外，还存在建筑物之间的互相影响。     

基于OpenFOAM的正弦形沙丘表面流场特性研究

基于开源计算流体力学软件OpenFOAM,对沙丘表面流场进行数值模拟,将所得到的沙丘表面风速分布与野外观测结果进行对比,验证数值模拟的准确性.在此基础上,获得了距沙丘表面不同高度的速度分布曲线和沙丘表面上的速度矢量场,探讨了数值模拟结果对计算网格的敏感性,考察了沙丘宽高比对流场特性的影响.结果表明:OpenFOAM的标准k-ε模型能够预测低矮正弦形沙丘表面流场的速度转变特性;在贴近沙丘表面的内层,沙丘表面特定高度处的速度分布曲线相对于沙丘轮廓存在相位提前,随着高度的增加,相位提前量减少,速度变化趋于平缓;网格数目对外层速度分布影响较小,但对内层的速度分布则有重要影响;足够多的网格数目是获得精确的沙丘表面流场信息的必要条件,在沙丘流场数值模拟中应当引起重视;保持沙丘宽度不变,随着沙丘宽高比的减小,沙丘表面速度分布曲线的相位提前量变大,背风坡出现回流区,宽高比越小,回流区的长度与强度越大.     

微通道气体流动的格子波尔兹曼模拟

通过隐式格子波尔兹曼方程,并采用壁面平衡边界条件以及二阶关系,模拟了微通道气体流动中的非线性压力和壁面滑移速度,模拟结果与Arkilic的解析结果十分吻合,验证了格子波尔兹曼方法在滑移流区的有效性。     

井底边界条件下喷嘴射流流的数值模拟

用数值模拟的方法研究了符合井底边界条件的轴对称单喷嘴淹没非自由射流流场，给出了不同喷嘴直径、不同喷速、不同喷距下的流场磁量图，并研究了上述工况下的射流轴心速度衰减规律。结果表明，除射流区、撞击区、漫流区、返回区之外，还存在一明显的回流旋涡区，射流出口雷诺数对轴心速度衰减的影响很小，而温流的存在命名得射流轴心速度在靠近井底的区域内受到很大程度的影响。     

井底边界条件下单喷嘴射流流场的数值模拟

用数值模拟的方法研究了符合井底边界条件的轴对称单喷嘴淹没非自由射流流场，给出了不同喷嘴直径、不同喷速、不同喷距下的流场矢量图，并研究了上述工况下的射流轴心速度衰减规律。结果表明，除射流区、撞击区、漫流区、返回区之外，还存在一明显的回流旋涡区，射流出口雷诺数对轴心速度衰减的影响很小，而漫流的存在使得射流轴心速度在靠近井底的区域内受到很大程度的影响．数值模拟结果和有关实验数据有较好的一致性。     

基于格子Boltzmann方法的固定方柱绕流研究

研究了格子Boltzmann方法的基本原理和边界条件处理方法,并利用该方法对不同雷诺数下固定单方柱绕流流场进行分析,探究了方柱产生卡门涡街的临界雷诺数和范围,验证了格子Boltzmann方法边界处理和数值模拟的正确性和便捷性.对雷诺数为200时并列双方柱不同分布间距的流场进行了模拟,结果表明,当柱间距为2倍方柱边长时,流体绕流方柱的涡流彼此影响最为明显.     

小尺度新月形沙丘背风侧流场特性的大涡模拟分析

为了准确获得新月形沙丘背风侧湍流流场流动特性,采用基于Smagorinsky的亚格子尺度涡黏模型的大涡模拟(LES)方法,对简化了的实验缩比小尺度沙丘模型绕流气相流场进行了数值研究,模拟了相同来流风速、不同沙丘迎风坡度和高度下背风侧风场的湍流流动模式,获得了背风侧区域地表面摩擦系数分布,比较了回流区长度和湍流强度分布.模拟结果表明:LES方法能较好地揭示背风侧湍流流场特性;流场重附距离随沙丘迎风坡度和高度的增加而增大;背风侧回流区内湍流强度总体上比回流区外大,回流区内沙丘坡脚位置及重附点位置处湍流强度最大值出现在贴近地表附近,随沙丘高度的增加而增大;回流区中部湍流强度最大值出现在回流区顶部,随沙丘高度的增加先增大后基本不变.     

用lattice Boltzmann方法模拟矩形竖腔内混合对流

用lattice Boltzmann方法(LBM)模拟了纵横比(高度与宽度之比)大于1的矩形竖腔内同时存在自然对流和强制对流的混合流动时流场与温度场分布.比较了不同雷诺数(Re)、瑞利数(Ra)下矩形腔的流场分布在不同纵横比时的特点,研究了这些参数值对混合对流换热的影响,提出混合对流存在时Re和Ra对流动和换热的影响可由Richardson数(Ri)综合反映,并得到在不同纵横比时平均Nusselt数(Nu)随Ra和Re的变化关系式.与常用的LBGK模型不同,用了LB力矩模型.     

用格子Boltzmann大涡模拟方法计算电磁力作用下的翼型绕流

利用格子Boltzmann大涡模拟(LBM-LES)方法,对较大雷诺数Re =2.4×105下翼型绕流的电磁控制进行数值研究.结果表明,LBM-LES方法计算过程简单,容易并行,适合处理该问题.     

Adaptive Lattice Boltzmann Model for Compressible Flows

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

格子Boltzmann方法中不同边界处理及不同体力项表达式的对比分析

简要介绍了格子Boltzm ann方法以及耦合温度的格子Boltzm ann方法,就LBM采用不同边界处理分别对二维方腔自然对流进行了模拟以及对比。同时,推导出LBM中更为完整的体力项的表达式,并采用简化的体力项与较完整体力项分别进行数值模拟,将结果进行了对比分析。模拟的结果与前人的结果吻合良好。     

多孔介质振荡流格子-Boltzmann模拟

为研究多孔介质内振荡流流动特性,该文将格子-Boltzmann方法对多孔介质内的振荡流过程进行了数值研究.采用四参数方法构造了多孔介质结构模型,对多孔介质内的振荡流过程的模拟结果表明;多孔介质振荡流的最大压差随最大Renoylds数的增加而迅速增加,在Renoylds数小于100的范围内,最大阻力系数与最大Renoylds数呈反比例关系,与文献报道的实验结果一致;在孔隙率小于0.85范围内,压差与流速之间的相位差变化较小,孔隙率继续增加时,相位差随孔隙率的增加而线性降低.     

用晶格Boltzmann方法模拟动脉分叉流场

简要介绍了晶格Boltzmann方法(LBM),它是模拟流体流动以及为复杂物理现象建模的一个新工具,比传统的数值模拟方法有许多独特的优点,特别是在处理复杂边界方面.并应用LBM对动脉分叉的流场在不同雷诺数情况下进行模拟,给出了速度、剪切应力和压力在分叉处的分布,并且分析了分支管内流体分离区的存在.流体在分离区内几乎处于停滞状态,剪切率非常小,容易在此处形成动脉粥硬化.     

基于格子玻尔兹曼方法的二维T型微通道内液滴生成

采用格子玻尔兹曼方法中的伪势模型,模拟了密度比近似相同且不互溶的两液相在T型微通道内产生不同形态的液滴,模拟结果与实验结果相吻合。分析了毛细数、两相流量比及两相的黏度比等因素对生成的液滴尺寸和间距的影响,模拟结果所得各种因素间的数值公式也与实验结果相差无几,连续相的毛细数和两相流量比对液滴的尺寸和间距的影响比较大,两相的黏度比对液滴的影响微小。     

用Lattice Boltzmann方法计算矩形柱的绕流问题

采用二维9速度Lattice Boltzmann模型, 数值模拟流体流过矩形截面柱体的绕流问题. 在边界处, 采用二阶精度的插值边界处理方法, 计算了流动的Strouhal数及柱体受到的升力和阻力系数, 给出了流场的流线和等涡线图. 使来流方向与矩形柱的长边方向平行, 计算结果表明, 改变矩形的长/宽比, 流场的Strouhal数随长/宽比呈线性变化.     

网格因素对格子Boltzmann方法误差的影响

基于外推法的思想,提出了格子Boltzmann方法中的误差计算方法.该方法在粗细两种网格下,采用Filippova提出的分布函数耦合方法推导出误差的计算公式,其主要误差项与格子大小和分布函数中的非平衡态部分成正比.最后用垂直平板障碍绕流作为一个数值计算实例,验证了理论推导与数值试验结果是一致的.     

利用格子Boltzmann方法计算页岩渗透率

针对页岩中孔隙主要分布在纳米尺度的特点,建立Knudsen数修正固体边界并考虑镜面反弹的格子Boltzmann模型。利用二维平板模型,考察Knudsen数对渗透率的影响,得出在一定误差条件下,须考虑Knudsen数修正的最大孔隙宽度。由页岩的SEM扫描图像重构得到三维数字岩心,利用新模型进行模拟流动,并计算得到绝对渗透率。结果表明:利用Knudsen数修正后,通道中间部分流体速度增大,固体边界处速度变小;压力不变,随着孔隙直径变小,Knudsen数增大,渗透率变小。