基于动态Smagorinsky-MRT-LBM的圆柱绕流湍流研究

采用多松弛格子玻尔兹曼与大涡模拟的动态Smagorinsky亚格子涡粘模型相结合的方法对圆柱绕流湍流进行模拟,研究300至10 000不同雷诺数下的单圆柱绕流湍流涡的周期性变化。研究结果表明:计算所得的斯考特数值和阻力系数值与文献中实验结果以及有限元仿真结果基本一致,说明该方法对于不同雷诺数的湍流模拟是合理的。     

近平板圆柱绕流的边界元分析

应用边界元方法,对低雷诺数粘性不可压流体的半无界近平板圆柱绕流问题进行数值模拟,得到圆柱表面的压力、面力分布规律,给出数值模拟的圆柱绕流的注场,数值计算与理论分析和实验结果比较完全一致。     

正三角形布置下三圆柱绕流的数值模拟

利用有限体积法和非结构性网格对正三角形排列下的三圆柱绕流进行了数值模拟.取雷诺数Re=100,并考虑4个不同的间距比(1.5,3,5,7),着重研究了间距比对绕流流场特性的影响,结果表明:在小间距比的情形,下游圆柱发生的偏流现象为单稳态偏流;在大间距比的情形,尾流效应使得上游圆柱脱落的剪切层在下游圆柱间上下波动.     

基于动态亚格子模型的方柱绕流大涡模拟

目的 对高雷诺数下的方柱绕流进行数值模拟,得到方柱下流场时间平均结果以及湍流流场的流场特征,用于验证亚格子模型对大涡模拟(LES)计算精确性的影响.方法 基于有限体积法(FVM)的流体计算开源软件OpenFOAM,采用两种不同动态亚格子模型下的大涡模拟方法,模拟高雷诺数(Re =2.2 ×104)下的方柱绕流.结果 通过数值计算得到流场完整的平均和脉动结构,并得到阻力系数、升力系数、斯托拉哈数等重要的气动力参数,经对比与试验实际情况吻合.结论 高雷诺数复杂流场下亚格子应力的封闭形式对大涡模拟的计算结果有重要的影响,采用动态模拟技术一方程模型计算结果优于常规的Smagorinsky模型,在高雷诺数复杂流场下的计算是高效的,为结构抗风的数值研究提供了参考.     

高雷诺数下并列双圆柱绕流的大涡模拟

为澄清并列双圆柱结构发生偏向流现象的流场机理,采用大涡模拟(LES)方法,在高雷诺数下(Re=1.4×10~5)研究了并列双圆柱的气动性能及其流场特性随圆柱间距比P/D(P为圆心间距,D为圆柱直径)的变化规律,重点探讨了小间距并列双圆柱的偏向流现象及其对圆柱气动性能的作用机理.研究结果表明:大涡模拟方法得到的气动力结果与文献风洞试验值吻合良好;随着并列双圆柱间距的增大,绕流场会呈现单一钝体、偏向流和平行涡街等多种流态结构;当P/D=1.1时,绕流场会间歇性地出现单一钝体和偏向流流态,两种流态的气动性能和流场特性有很大差异,圆柱的气动力会随时间发生剧烈变化,呈现非稳态特征;当P/D=1.2～1.5时,绕流场呈现偏向流流态,两个圆柱的气动力和尾流呈现不对称现象,偏向流的偏转方向会出现间歇性地变化,尾流涡脱强度弱,气动力脉动小;当P/D=2～4时,绕流场总体呈现平行涡街流态,尾流涡脱强度强,气动力脉动大,气动干扰减弱.     

二维格子Boltzmann方法圆柱绕流湍流模拟

采用多松弛格子Boltzmann方法的大涡模拟技术对二维圆柱绕流问题进行了数值模拟.运用二阶插值反弹格式处理圆柱边界,同时结合Wall-Adapting Local Eddy-viscosity模型(WALE)对壁面附近湍流粘性进行修正,测算了其升、阻力系数和涡脱落频率.结果显示:在雷诺数小于300时,升、阻力系数以及涡脱落频率均跟实验值吻合较好;对于较高Re数问题,其三维效应不能忽略,模拟结果跟实验值有所偏差,但通过引入WALE模型进行壁面修正,该方法较其他二维方法在估算阻力系数上精度有明显提高.     

粘性流体圆柱绕流的有限元模拟

对粘性流体圆柱绕流进行了数值模拟。根据三阶Taylor展开弱解的概念,在二阶分裂步Taylor-Galerkin有限元法基础上,采用时间推进和张量分析的方法推导了对流扩散通用微分方程的ETG有限元离散格式。对Reynolds数为32、102和250的粘性流体圆柱红流流场进行了数值仿真,数值模拟结果与实验结果吻合良好。表明隐含流线迎风耗散作用并具有三阶部分展开的Taylor-Galerkin有限元法稳定性好,精度较高,可用于计算包含对流扩散耦合传递动量的绕流问题及类似的流动问题。     

圆柱绕流涡致振动的平面湍流数值模拟

为了研究工程中存在的湍流涡致振动问题,应用平面湍流k-ε模型和重整化群k-ε模型,结合分块耦合方法及任意拉格朗日欧拉法（ALE）数值模拟了圆柱绕流涡致振动,求解基于非交错网格系统, 其中包含满足连续性方程的压力泊松方程解法.计算了湍动能、湍流耗散率、湍流黏性系数及柱面涡量的变化情况.结果表明,当雷诺数为5 000～2×105时计算所得的阻力系数与实验数据吻合良好.研究发现,在涡致振动情况下的湍流耗散率比固定圆柱涡脱落情况下小得多,湍流黏性系数与雷诺数呈对数律增长.     

有限长波浪形圆柱绕流数值模拟

为研究一端固定在壁面上,另一端为自由端的有限长波浪形圆柱的减阻抑振效果,对有限长波浪形圆柱进行了计算.首先,采用大涡模拟数值模型的方法,对雷诺数为3 900时的有限长直圆柱和不同波长、波幅组合后的12种有限长波浪形圆柱进行了计算;其次,针对计算结果进行后处理,得到并比较了不同组合形式下的有限长波浪形圆柱的升、阻力系数大...     

流体横掠圆柱体传热特性数值模拟的流态模型选择

采用层流模型、标准k-ε紊流模型、RNG k-ε紊流模型,对亚临界雷诺数条件下(Re=300~300 000)的二维圆柱扰流进行了数值模拟与仿真,得到了圆柱壁面平均努塞尔数随雷诺数变化的规律。并得出结论:当Re≤3 900时,采用层流模型较合适;当3 900≤Re≤22 000范围内时,采用标准k-ε紊流模型较合适;当22 000≤Re≤300 000时,采用RNG k-ε紊流模型较合适。