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

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

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

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

熵格子Boltzmann方法的亚格子尺度模型

熵格子Boltzmann方法对格子Boltzmann方法的改进使其具有稳定模拟高雷诺数湍流的潜力。尝试将大涡模拟的重整化群亚格子尺度湍流模型引入到熵格子Boltzmann方法中。同时建立基于熵格子Boltzmann方法的大涡模拟重整化群亚格子尺度模型及大涡模拟标准Smagorinsky亚格子尺度模型,进行高雷诺数湍流圆柱绕流场的模拟计算。结果表明:所建立的熵格子Boltzmann方法重整化群亚格子尺度模型能够有效模拟高雷诺数湍流流动,其对高雷诺数湍流的模拟计算结果略优于熵格子Boltzmann方法的大涡模拟标准Smagorinsky亚格子尺度模型。     

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

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

基于PINN方法的不可压流场求解

将残差网络（Res Net）结构引入物理信息神经网络（PINN）,提出一种基于Res Net的PINN方法（Res Net-PINN）。采用该方法对二维不可压圆柱绕流尾迹流场进行重建和预测,结果表明：Res Net-PINN能更准确地重建绕流尾迹的非定常变化规律;在圆柱绕流尾迹的短期预测方面,Res Net-PINN的预测精度和收敛速度相比全连接PINN均能提升3倍左右。Res Net能提高PINN对非定常流场的求解和预测能力,该研究可为采用机器学习方法求解更复杂流动问题提供参考。     

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

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

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

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

舰船开孔流动噪声控制措施研究

基于计算流体软件,采用大涡模拟方法计算了给定来流速度条件下表面开孔的非定常流场,分析了流场的非定常流频谱特性,采用边界元方法模拟了噪声源特征,对各种流动噪声控制措施进行了对比分析。计算结果表明,开孔形式的改变会直接影响流噪声的水平,设计结果可以为舰船流水孔的降噪设计提供参考。     

舰船开孔流动噪声控制措施研究

基于计算流体软件,采用大涡模拟方法计算了给定来流速度条件下表面开孔的非定常流场,分析了流场的非定常流频谱特性,采用边界元方法模拟了噪声源特征,对各种流动噪声控制措施进行了对比分析.计算结果表明,开孔形式的改变会直接影响流噪声的水平,设计结果可以为舰船流水孔的降噪设计提供参考.     

临界雷诺数下圆柱涡激振动抑制措施研究

拖缆、海洋立管等在流体力的作用下,会发生涡激振动,容易使这类海洋结构物受到破坏。本文将这一力学现象简化为绕二维圆柱流动问题,通过离散求解二维非定常平均雷诺应力方程(URANS),分析临界雷诺数时(8.21×104Re1.54×106)的流动特性,从而得到绕圆柱流动的涡脱落过程、涡激力与雷诺数之间的关系。本文计算结果和已发表的结果有良好的一致性。在此基础上,本文提出了一种在圆柱体后方安装分隔板来抑制涡激力的简单措施,并通过数值方法验证其有效性。     

二维圆柱绕流数值模拟

对于椭球体绕流问题，可以考虑利用二维切片法的思想将其剖成一个个圆剖面，研究圆剖面的二维流动问题即二维圆柱绕流问题。采用Fluent软件，分别用多种流动模型在高雷诺数下对定常和非定常状态的圆柱绕流进行了数值模拟，计算得到圆柱绕流的升阻力系数、分离点位置、Strouhal数、流函数等结果，并与文献中的实验结果进行了比较，得到相应比较合理的计算模型。     

三角形布置圆柱体群绕流特性与流动机理研究

基于浸没边界-多松弛-格子玻尔兹曼方法,对Re=100时等边三角形布置三圆柱体群绕流问题进行数值模拟,分析间距比(K_d)与来流角度(α)两个关键参数对流体力系数及流场特性的影响,揭示其流动机理。计算结果表明：不同流场下尾流模式分为单漩涡、单漩涡体向双漩涡过渡、不规则的双漩涡、规则的双漩涡、双漩涡向三漩涡过渡、不规则的三漩涡和规则的三漩涡模式,在α=30°时流场下游涡街更加规则,尾流模式转变更快。三圆柱所受的时均流体力系数受间隙流流速影响较大,各来流角度下随间距比(除小间距比外)的变化趋势基本一致,尾流模式转变会导致三圆柱群的流体力系数均方根值发生显著变化,在K_d≥3.5后上游圆柱流体力系数均方根值有较大起伏。除α=0°工况外,在K_d≥2.5后中游圆柱的涡脱落频率受α和K_d影响较小,数值相对稳定。     

氢气泡DPIV技术在圆柱绕流流场测定中的应用

将氢泡示踪法与数字图像分析技术结合,设计了一套粒子图像流场观测装置,对不同雷诺数下的圆柱绕流流场进行了观测。实验中清晰观察到圆柱绕流的卡门涡街现象。采用基于互相关算法的分析软件对图像进行处理,得到流速矢量图、等值线图及涡量图。根据图像及计算结果,分析了一个拟周期内圆柱后涡旋的发展、演化过程。实测结果表明,随着雷诺数的增大,圆柱绕流近尾迹区长度逐渐缩短。将流动特性实测结果与相关文献进行了对比,结果基本一致。     

格子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倍。     

基于边界元方法的气泡脉动诱导壁压特性

针对圆柱壳结构在水下爆炸气泡脉动载荷作用下的诱导壁压特性,采用Geers-Hunter模型得到气泡脉动运动,在势流假设下采用边界元方法建立了水下爆炸气泡脉动载荷壁压计算方法,通过对不同水下爆炸工况进行计算分析发现:圆柱壳表面的气泡诱导壁压有明显的绕射现象,在爆距较小时迎爆面气泡脉动诱导压力峰值要远大于背爆面压力;在圆柱壳轴向上随距中心点距离的增大绕射特性逐渐减小,壁压向自由场压力逼近。所得到的方法和规律可为潜艇抗冲击设计和评估人员提供参考。     

基于DES方法的三维圆柱受迫振动数值模拟

基于OpenFOAM平台对雷诺数3900下的三维固定圆柱绕流和雷诺数30000下的三维圆柱受迫振动进行数值模拟,对湍流的模拟采用SA-DDES离散涡模型.针对三维圆柱绕流问题,采用PISO算法耦合求解速度-压力场,对圆柱绕流的基础参数如St数、平均阻力系数和平均升力系数等与试验进行了比较,结果吻合良好.针对三维圆柱受迫振动问题,采用PIMPLE算法耦合求解速度-压力场,对不同频率下振幅比为0.3的阻力系数、尾涡形态进行了计算,阻力平均系数与实验值相比吻合较好.     

亚临界雷诺数下圆柱和方柱绕流数值模拟

基于RNG k-ε模型,采用有限体积法对亚临界雷诺数条件下(Re=3×102~3×105)的二维单圆柱和单方柱绕流进行数值模拟与仿真。得到了圆柱和方柱绕流阻力系数Cd与Strouhal数随雷诺数的变化规律。同时对雷诺数Re=22 000下的圆柱、方柱绕流相关特性进行详细对比分析。计算结果表明:同一雷诺数下,单圆柱绕流阻力系数Cd较单方柱低,但圆柱的Strouhal数较单方柱则要高。虽然二者边界层分离点不同,但流场的演变与漩涡的脱落具有一定相似特性。     

并列双圆柱绕流中宽窄尾流的识别方法

采用基于边界瞬时涡量守恒(IVCBC)的纯拉格朗日涡方法,结合并列双圆柱的结构特点,建立双圆柱绕流数值计算模型。在间隙率T/D=1.1~2.6时,并列双圆柱绕流会出现宽尾流(WW)和窄尾流(NW)现象。利用两圆柱间隙中心点的诱导速度方向与间隙流的偏转方向同步的特性,提出一种区别宽、窄尾流的数值计算方法;比较间隙流的偏转方向与该间隙中点的速度偏转方向,并进行典型的算例分析。结果表明:该方法能准确区别宽尾流和窄尾流,可为数值研究这一特殊区域的流体特征提供重要的数值计算方法。     

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

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

基于ADINA有限元软件的圆柱绕流三维数值模拟

利用ADINA有限元计算软件建立基于N-S控制方程的圆柱绕流三维数学模型,模拟不可压缩粘性流体中单个固定圆柱的绕流问题,计算不同雷诺数(Re)条件下圆柱绕流的压力系数、升阻力系数和Strouhal数等相关参数.通过与物理模型试验结果的对比对该数学模型进行验证,进一步分析在大雷诺数条件下圆柱绕流的三维效应.