熵格子Boltzmann方法的重整化群代数湍流模型

重整化群理论所建立的湍流模型能够最大程度地减小模型经验性,因此文章尝试将重整化群代数湍流模型引入到熵格子Boltzmann方法中,建立新型的计算模型以对高雷诺数湍流进行模拟研究。同时为了进行比较研究,还建立了熵格子Boltzmann方法的标准大涡模拟模型。完成了对高雷诺数湍流绕流场的模拟计算。结果表明:所建立的熵格子Boltzmann方法重整化群代数湍流模型能够有效地模拟高雷诺数湍流流动问题;其对紧贴壁面处较小尺度湍涡的模拟结果趋近于大涡模拟的结果;重整化群代数湍流模型在对高雷诺数湍流的模拟中表现出耗散模型的特征。     

结合亚格子模型的Lattice-Boltzmann算法

介绍一种亚格子模型与LB(Lattic-Boltzmann)方法相结合的算法,并用其对亚临界雷诺数情况下的圆柱绕流问题进行数值模拟,得到流场的流线图以及圆柱表面的压力系数。通过与实验值相比较,发现数值结果与实验值是比较吻合的,证明该算法可以用于计算雷诺数较大的流体问题,同时具有LB方法计算效率高的优点。     

结合大涡模拟的格子玻尔兹曼方法模拟高雷诺数流动

为了将格子玻尔兹曼方法（Lattice Boltzmann method）推广到高雷诺数流动的数值模拟上,文章研究了将大涡模拟的Smagorinsky模型及其改进的IR模型应用到LBM方法中并发展了LBM-SM模型和LBM-IR模型,其后以顶盖驱动方腔流动（Re=1000;10000;100000）和后台阶流动（Re=389;1000;10000;100000）为标准进行数值模拟与比较。数值结果说明该方法可以应用于高雷诺数流动的数值模拟,并克服了传统LBM方法在模拟高雷诺数流场时会出现非物理振荡的情况。通过分析对比,不仅得到了两种模型能够稳定计算且不失真的参数范围,并且认为LBM-IR模型有效参数的选取范围更为广泛。     

格子Boltzmann方程模拟高雷诺数三维方腔流

为分析方腔流内部流场的特性和验证格子Boltzmann方法模拟湍流的能力,应用标准Smagorinsky涡粘性模型与多松弛时间格子Boltzmann方程(Multiple Relaxation Time Lattice Boltzmann Equation,MRT-LBE)组合对高雷诺数(Re=10 000)三维方腔流进行数值研究,计算了时间平均量,如速度,均方根速度、雷诺应力以及中心断面(y=W/2)处的流线等高线。模拟结果与已有实验和数值模型结果比较可知,MRT-LBE能够精准地计算剪切驱动方腔内流场的变化。另外,将基于图形处理(graphic processor unit,GPU)的计算统一设备架构(Compute Unified Device Architecture,CUDA)并行技术引入到基于MRT-LBE的Smagorinsky模型以提高计算效率,计算效率提高达200倍。     

疏水表面流体流动特性的格子Boltzmann方法模拟

采用格子Boltzmann方法研究了微形貌对固体表面润湿性的影响,在此基础上进一步模拟了具有微形貌的疏水表面通道内的流体流动,从法向速度、剪应力、滑移速度等角度分析了疏水表面的流场特性,揭示了疏水表面滑移流动的产生机制。结果表明,疏水表面的滑移流动是由低表面能作用和微形貌共同引起的。具有微形貌的疏水表面比光滑疏水表面具有更好的减阻效果,原因在于微形貌能够驻留气体,形成的气液自由剪切面加剧了疏水表面的滑移流动,最大滑移速度可以达到主流平均速度的50%左右。     

基于格子Boltzmann方法的液润表面减阻规律

[目的]近年来,液润表面(LIS)作为一种新型的减阻表面被提出.它将传统疏水表面微沟槽中残存的气体替换为润滑油,进而提高了减阻效果的稳定性.为了更全面地认识液润表面,研究润滑油溶解性对滑移长度的影响.[方法]基于格子Boltzmann方法(LBM)伪势模型,对液润表面的滑移现象进行数值模拟,研究润滑油溶解密度和外部剪切...     

基于格子Boltzmann方法的船舶风力助推转子绕流特性

风力转子推进装置节能环保,降低推力成本,有利于船舶行业的绿色发展。基于多松弛格子Boltzmann方法,对并列风力助推转子绕流进行数值模拟。首先模拟单圆柱绕流的流动特性,以验证程序的可靠性,重点探究双圆柱的间距比和转速比对圆柱绕流特性的影响。获取圆柱的升阻力系数以及尾流流型,验证临界转速比的存在。结果显示:旋转可以有效地抑制涡的生成和脱落,当转速比达到临界转速比时,漩涡彻底消失,流场变得稳定;时均升力系数的绝对值和时均阻力系数随转速比的增大分别增大和减小。     

浸没边界格子Boltzmann方法的改进及转动圆柱绕流模拟

针对浸没边界格子Boltzmann方法计算效率不足的问题,提出一种改进浸没边界格子Boltzmann方法。通过使用作用力格子Boltzmann模型,并简化流固耦合作用力的计算方法,使计算流程得到简化,同时,对比原方法转动圆柱绕流计算结果,计算时间缩短一半以上,提高了数值计算效率。结合改进算法研究转动圆柱绕流的流场升阻力系数、压力系数、流场速度等流体参数场随转速比的关系,揭示了一定的转动圆柱绕流动力学规律。研究结果表明,改进方法准确可靠。     

格子Boltzmann方法在串列双圆柱绕流数值模拟中的应用研究

基于格子Boltzmann方法,对二维静止串列双圆柱绕流进行了数值模拟,并将模拟结果与已有研究结果进行了对比分析。为提高计算效率和计算精度,采用了多块网格耦合算法,并在圆柱曲边界处采用了较为精确的边界处理方法。提取了圆柱的升阻力系数,讨论了圆柱间距对圆柱受力情况和尾流特征的影响。数值模拟在雷诺数Re=200条件下进行,对两圆柱中心间距为1.5~4.0D(D为圆柱直径)之间的典型间距进行了数值模拟,获得了圆柱升阻力系数以及尾流中涡和流线的变化,验证了临界间距的存在,模拟结果和已有研究结果符合得较好。     

多孔介质内溶解与沉淀过程的格子Boltzmann方法模拟

  目的  腐蚀会严重影响舰船设备的使用寿命,了解腐蚀机理及对腐蚀进行预测能有效地减缓腐蚀对舰船设备的危害。  方法  构建了格子Boltzmann腐蚀模型,该模型可描述包含多相多组分流动与扩散、电化学反应和金属点蚀的稳态全过程。应用此模型,研究浸没于液体腐蚀环境中的金属表面单坑点蚀稳态过程;分析腐蚀化学反应速率、腐蚀溶液扩散系数、腐蚀产物扩散系数对腐蚀程度的影响。通过数值模拟,获得金属表面点蚀坑的形貌变化特征。  结果  结果显示,对于金属表面单坑点蚀的稳态过程,由于钝化膜与金属基体构成了“大阴极,小阳极”的电化学腐蚀体系,使腐蚀稳态点蚀将持续向金属材料基体的纵深发展,并且初生蚀孔除自身发生点蚀外,还将在蚀孔底部产生次生蚀孔。在不同影响因素下,腐蚀程度随腐蚀反应速率和反应物组分扩散系数的增大而增大,而随腐蚀产物扩散系数的增大而减小。  结论  通过该腐蚀模型,可模拟出与真实金属相近的腐蚀形貌变化特征。     

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

基于双分布格子Boltzmann模型,建立了适合于流体流动和换热的热格子Boltzmann模型.温度分布函数中采用D2Q9离散速度模型.以热格子Boltzmann模型,模拟了方腔内自然对流的形成及其演化,通过与相关文献的计算结果对比可以发现,热格子Boltzmann模型在处理流体流动与传热方面存在着独特的优点,文中建立的数值模拟计算方法和程序是切实有效的.     

浸没边界格子Boltzmann方法的改进及转动圆柱绕流模拟

  目的  为探究有限长圆柱绕流的流动机理和特性,  方法  采用大涡模拟（LES）数值模型并结合涡识别方法,对三维有限长圆柱绕流进行数值模拟,并对有限长圆柱绕流进行验证和分析。  结果  模拟结果表明,有限长圆柱回流区相对较短,自由端的下洗作用会扰乱卡门涡街,导致阻力系数损失;相对于固定壁面,自由端面对顺流向速度影响更大;自由端面\     

基于Boltzmann方法的手动可调螺距桨敞水特性

针对船舶长期服役后出现的船、机、桨不匹配问题,采用Boltzmann方法建立手动可调螺距螺旋桨的数学模型,对手动可调螺距螺旋桨进行敞水特性计算,结果表明:手动可调螺距螺旋桨随着螺距角的增大,其推力、转矩逐渐增大,当进速较小时,桨叶螺距角小的效率高,进速较大时桨叶螺距角大的效率高,并且桨叶螺距角有一定的工作范围,另外随着...     

基于格子Boltzmann方法的船舶风力助推转子绕流特性

风力转子推进装置节能环保,降低推力成本,有利于船舶行业的绿色发展。基于多松弛格子Boltzmann方法,对并列风力助推转子绕流进行数值模拟。首先模拟单圆柱绕流的流动特性,以验证程序的可靠性,重点探究双圆柱的间距比和转速比对圆柱绕流特性的影响。获取圆柱的升阻力系数以及尾流流型,验证临界转速比的存在。结果显示:旋转可以有效地抑制涡的生成和脱落,当转速比达到临界转速比时,漩涡彻底消失,流场变得稳定;时均升力系数的绝对值和时均阻力系数随转速比的增大分别增大和减小。     

平板壁面湍流脉动压力及其波数-频率谱的大涡模拟计算分析研究

壁面湍流脉动压力是重要的流噪声声源,对壁面湍流脉动压力及其波数-频率谱进行数值计算是流声耦合领域的重要课题,开展相应的研究十分必要。文章采用大涡模拟方法（LES）结合动态亚格子涡模型（DSL）与千万量级的精细网格,对平板壁面湍流脉动压力及其波数-频率谱进行了数值计算,并与试验结果进行了对比分析,验证了数值计算方法的可靠性。首先,介绍了大涡模拟的物理内涵与基本方程,给出了所采用亚格子涡模型的表达式。其次,描述了Abraham试验中矩形试验段的几何特征,给出了网格的剖分形式,并给出了相应的离散求解数值方法以及边界条件的设置。再次,探讨了湍流脉动压力变化规律及其相似律,基于Fourier变换计算得到了湍流脉动压力波数-频率谱,并详细讨论了壁面湍流脉动压力及其波数-频率谱计算值与试验值之间的差异,进行了定量与定性的验证分析,结果表明,计算结果与试验结果吻合良好,计算方法合理可靠,为今后复杂几何模型壁面湍流脉动压力及其波数-频率谱的计算研究工作奠定了基础。最后,基于试验和计算结果,比较分析了常用波数-频率谱理论模型,为波数-频率谱的工程应用提供了参考。     

基于LBM方法的圆盘等速入水空泡的数值模拟

采用单相格子Boltzmann模型进行了物体垂直入水问题的研究.该模型适用于大密度比的气液两相流动,其特点是忽略了系统中气相对液相的动力学影响.模拟了圆盘垂直等速入水过程,给出了入水过程中空泡与自由面的变化规律,研究了圆盘入水空泡相对闭合深度与Froude数之间的关系.数值模拟结果与相关的实验结果吻合良好,验证了格子Boltzmann方法(LBM:Lattice Boltzmann Method)模拟物体入水问题的可行性.     

基于格子Boltzmann方法的航道整治建筑物影响下过饱和总溶解气体输移释放规律

基于双分布函数-浅水Lattice Boltzmann方法,建立了过饱和总溶解气体(TDG)输移模型,研究某电站下游航道整治建筑物对河道内TDG输移规律的影响。结果显示,整治建筑物的修建对主流区流态和TDG变化影响较小;航道整治建筑物对河道过饱和TDG时空分布的影响主要表征于整治建筑物局部区域范围内,整治建筑物坝田回流区流速降低,水体传质减弱,从而导致TDG浓度降低且低浓度区域面积增加,可以为鱼类提供一定避难空间。该研究为水体中TDG输移扩散的研究提供了一种新的方法,同时为改善电站泄水导致TDG过饱和对鱼类的不良影响和建设生态水道提供科学依据。     

基于IVCBC涡方法高雷诺数下串联双圆柱的数值计算模型及其尾流特征研究

针对串联双圆柱的结构特点,利用IVCBC涡方法适合高雷诺数下数值计算的特点,建立了高雷诺数下串联双圆柱绕流的数值计算模型;采用经典算例验证了IVCBC涡方法的收敛性;探索了雷诺数为2.5×104、间隙比分别为1.1,1.25,2,2.5,3,3.25和4的串联双圆柱绕流的尾流特征;清晰地展示了尾流中较小的漩涡的形成、分裂和融合,详细地阐述了串联双圆柱流体特征发生突变的原因;深刻地揭示了串联双圆在高雷诺数下的绕流机理。研究表明:尾流模式与经典的实验尾流模型吻合较好,两圆柱中间的涡对是串联双圆尾流发生突变的主要原因;间隙之间时,间隙上部小漩涡形成、间隙中间流体的振动与下游圆柱表面上涡的脱落是同步的;该流体模式能清晰地展示尾流中较小漩涡,说明与有网格方法比较,该计算模型具有较大的优越性,为进一步研究高雷诺数下串联双圆柱流体力的特征提供了重要的研究工具。     

二维浅水流动的旋涡动力学问题探讨

  目的  高雷诺数下水翼绕流问题因其原理复杂、计算困难,一直是水动力领域关注的热点。基于涡量-流函数的离散涡（DVM）模型是一种无网格方法,可有效克服数值粘性问题。同时,在涡量集中区分辨率高,可有效模拟高雷诺数流场中的旋涡运动。使用一种随机涡方法对高雷诺数下二维刚性水翼涡发放频率进行预报。  方法  以NACA 66系列水翼为例,与试验数据对比验证其对涡发放频率预报的准确性。应用此方法分析高雷诺数下来流速度及初始攻角大小对水翼涡发放频率的影响规律。  结果  结果表明：涡发放频率会随来流速度呈非线性增加;随着水翼攻角的增加,泄涡尺寸和涡强不断增大,泄涡频率降低。  结论  来流速度和水翼初始攻角对二维水翼涡发放频率有显著影响,可为螺旋桨桨叶二维剖面的涡发放机理探究,甚至全桨的振动噪声预报提供参考。     

用基于SST模型的DES方法数值模拟圆柱绕流

分离涡数值模拟方法(DES)在物面附近求解雷诺平均Navier-Stokes方程,在其他区域采用大涡模拟方法.兼具前者计算量小的优点和后者能模拟大分离湍流流动、计算精度高的优势.采用基于剪切应力传输(SST)两方程湍流模型的DES方法对粘性不可压缩流体的圆柱绕流问题进行数值模拟.通过对所得的速度场、压力场、阻力系数Cd、升力系数Cl、斯特罗哈数St等结果的分析及与文献上的实验和计算数据的比较,说明DES方法对于层流及高亚临界雷诺数的圆柱绕流模拟是合理的.