提高k-ωSST模型对翼型失速特性的模拟能力

采用CFD软件fluent,k-ω SST模型对S825翼型进行了二维数值模拟.针对k-ω SST模型对翼型分离失速特性模拟不准的情况,对分离区域内湍流的强非平衡输运特性进行了分析.研究表明:原始模型由于不能准确模拟分离区内湍流的强非平衡输运特性,导致其对翼型失速特性模拟失效;提出了通过修正模型系数a₁及β_∞^*,提高k-ω SST模型预测非平衡湍流输运特性的方法,从而提高对翼型失速特性的模拟精度.     

超声速复杂流动中湍流模型的性能评估

对8°,16°,20°,24°超声速二维压缩拐角进行了数值模拟,系统评估了SA(Spalart-Allmaras), k-ω, SST(Shear Stress Transport)3种工程常用湍流模型对激波/边界层干扰复杂流动的模拟适用能力.对于"弱"干扰、小分离,工程常用湍流模型能够较准确地预测其压力、摩擦力、热流分布,而对于"强"干扰、大分离则会造成较大的偏差.同时,随着分离区的增大,各湍流模型的适用性能差别更加明显,数值方法也有一定的影响.      

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA (Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.     

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-Stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA(Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.     

用GAO-YONG湍流模型计算激波/边界层干扰

采用GAO-YONG可压缩湍流方程组,模拟了平板激波/湍流边界层干扰现象.运用SIMPLE算法求解方程组,并分别采用三阶精度的QUICK格式和中心格式离散对流项和扩散项.计算结果较好预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征: 分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波.对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布和速度型与实验值比较吻合很好.结果表明GAO-YONG可压缩湍流方程组能够高精度模拟平板激波/湍流边界层干扰流动.      

使用GAO-YONG湍流方程组对扩压器流动的计算

采用基于同位网格的SIMPLE方法求解GAO-YONG不可压湍流方程组,对二维扩压器流动进行了数值模拟.通过界面速度动量插值法解决了压力锯齿波问题,并且采用了正交贴体网格和二阶QUICK格式离散对流项以提高计算精度.与实验结果以及BL(Baldw-Lomax)模型计算结果的比较表明,不需要任何经验系数及壁面函数的GAO-YONG不可压湍流模型方程组能够对有压力梯度的湍流流动做出很好的预测.计算发现,在机械能方程中引入平均流压力梯度的作用,对GAO-YONG湍流方程组正确模拟逆压力梯度流动起到了关键作用.      

基于非湍流脉动动能方程的高超声速转捩预测

通过求解基于稳定性理论的非湍流脉动动能输运方程,预测了来流Ma为5.91的裙锥绕流流动的层流-湍流转捩点,并采用SST(Shear-Stress Transport)湍流模型结合代数间歇因子模型对转捩流场进行了数值模拟.结果表明非湍流脉动动能输运方程模型可以很好地捕捉逆压梯度对转捩前期边界层内的不稳定波动频率的影响.等温壁面条件下计算得到的转捩位置与静风洞实验结果基本一致.壁面冷却促进第2模态不稳定波动的增长,使转捩提前发生.代数间歇因子模型模拟高马赫数流动的转捩区长度较短,温度峰值偏低.     

PANS方法在双圆柱绕流数值模拟中的性能分析

PANS方法以传统的雷诺平均Navier-Stokes(RANS)方程为模板,引入控制参数对原方程进行修改,是一类能够有效模拟分离湍流流动的混合方法.采用基于Menter 剪切应力输运(SST)湍流模型的PANS方法,选取双圆柱绕流作为研究对象, 分析双圆柱流动结构,验证方法的准确性与可靠性,考察PANS方法中模型参数选取对于计算结果的影响,评估不同类型PANS方法的性能.研究表明:模化湍动能比例这一模型参数对计算结果有显著影响,且计算域内统一的模型参数取值难以处理复杂流动问题.在此基础上使方法中模化湍动能比例取值随流场物理信息和网格尺度变化,所得结果与SST DES计算结果及实验结果吻合较好,表明该方法适用于复杂湍流数值模拟.      

基于支板燃烧室的喷管化学非平衡效应

采用有限体积法全隐式格式和代用燃料C₁₂H₂₃的10组元13步化学反应Arrhenius有限速率模型研究煤油燃料超燃冲压发动机单边膨胀喷管(SERN, Single Expansion Ramp Nozzle)内的化学非平衡流动,通过建立支板燃烧室——喷管模型有效解决了单边膨胀喷管模型的"入口薄层"问题.计算结果表明,整个单边膨胀喷管内,流动呈现化学非平衡效应,喷管入口附近区域尤其显著;非平衡流动喷管性能明显高于冻结流动,随发动机当量混合比 ε增加,非平衡流动的喷管推力系数和升力系数相对冻结流动的百分比增量 δ不断升高,当 ε=0.8时,推力系数百分比增量 δ_CF达到9.41%,升力系数百分比增量 δ_CY达到16.39%,化学非平衡效应对煤油燃料超燃冲压发动机尾喷管性能的影响不可忽略.     

XY-SAS模型中不同网格尺度限制器的影响分析

XY-SAS(Xu & Yan Scale Adaptive Simulation)模型通过在SA(Spalart-Allmaras)模型耗散项的湍流尺度中引入基于平均涡量的冯卡门长度尺度,在不稳定流动区域具有类似LES(Large Eddy Simulation)的计算性能,且不存在DES(Detached Eddy Simulation)方法中RANS(Reynolds Averaged Navier-Stokes)与LES的交界面问题.由于SAS类模型预测的湍流能谱存在高波数衰减不足的缺陷,需要进一步在湍流尺度中引入网格尺度作为限制器.采用不同形式的网格尺度限制器构造了3种XY-SAS模型,对Re=3 900的圆柱绕流进行数值模拟,并与DES及实验数据进行比较.结果表明网格尺度限制器对XY-SAS模型的计算性能有较大影响,以最小网格间距作为高波数衰减限制器是较为合理的选择,此时模型对于大分离流动能够达到DES的模拟能力.     

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

SST-DES在小分离流动数值模拟中的改进

在基于剪切应力运输（SST,Shear Stress Transport）分离涡模拟（DES,Detached Eddy Simulation）应对小分离流动的数值模拟中,上游雷诺平均方法（RANS,Reynolds Averaged Navier Stokes）区域内较大的涡黏性通过对流扩散方程传递到下游大涡模拟方法（LES,Large Eddy Simulation）区域,导致自由剪切层解析湍流解迟迟得不到发展,剪切层失稳延后,出现强烈的灰区效应。为此将SST-DES中体现LES求解的类亚格子模式重新构造成标准k方程亚格子模式,防止多变量在RANS和LES交界面处对湍流黏性系数的混合影响,最终得到了一种改进的SST-DES方法。采用该方法对AS239翼型小分离数值模拟并与DES和延迟分离涡模拟（DDES,Delayed Detached Eddy Simulation）结果进行对比,可以看出该方法在壁面边界层保留了SST-DDES特性的同时,在尾迹区大大提高了LES的可解能力,显著降低了灰区的影响。     

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

近壁区Reynolds应力输运方程的理论模型

在湍流边界层的近壁区,采用对称共振三波的理论模型描述相干结构,根据理论模型对Reynolds应力输运方程逐项进行计算和分析,结果与直接数值模拟符合很好.这不仅在理论上有益于对湍流物理机制的了解,而且展现了一种可能性,即根据相干结构的理论知识来改造湍流模型,使之具有更清晰的物理内涵以提高近壁区的预估精度.     

高超声速热流计算湍流模型性能评估

采用计算流体力学方法,通过对高超声速来流下双椭球热流计算,综合分析了影响壁面热流预测精度的因素.针对工程上4种常用的湍流模型:BL(Baldwin-Lomax)模型、SA(Spalart-Allmaras)模型、k-ω模型及SST(Shear-Stress Transport)模型,评估了其在高超声速热流预测中的性能,获得了湍流模型下壁面法向网格雷诺数对热流计算的影响规律.此外,对高马赫数下不同后掠角钝舵模型进行热流计算,获得了后掠角对钝舵前缘热流峰值的影响规律.研究表明,两方程湍流模型相比一方程和零方程模型更适合高超声速气动热计算,其中SST模型对热流预测精度更高,总体性能更为优异;钝舵驻点线最大热流随后掠角增大呈非单调变化,在后掠角约为22°时热流达到峰值.      

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

合成射流控制NACA0015翼型大攻角流动分离

为了研究合成射流激励器处于NACA0015翼型回流区时对其分离流动的控制,采用商用计算流体力学软件Fluent 6.1求解Reynolds平均Navier-Stokes方程,通过对翼型气动力特性、脱落漩涡结构以及射流孔口附近流动结构的分析,揭示了合成射流处于分离区时对边界层控制的机理.结果表明:当合成射流孔口处于回流区时仍可有效推迟翼面边界层分离点,缩小回流区范围,从而有效提高翼型的升力.当射流方向垂直于壁面,无量纲频率以及吹气速度比都等于1时,翼型平均升力系数提高40%左右.