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方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.     

高超声速流动湍流模式评估

选取了超声速二维膨胀-压缩拐角和高超声速双椭球绕流作为基准流动,考察了几种当前CFD(Computational Fluid Dynamics)工程应用较流行的湍流模式:BL(Baldwin-Lomax)模式、SA(Spalart-Allmaras)模式、k-ω模式以及SST(Shear-Stress Transport)模式,通过数值计算结果和实验结果的对比及分析,对有关的湍流模式进行了评估,得到一些有意义的结论,并指出受流动可压缩效应的影响,计算高超声速湍流壁面热流时要对现有的湍流模式进行相应的可压缩修正.      

高超声速尖锥边界层转捩数值模拟

边界层转捩对高超声速飞行器气动力和热产生重要影响.通过联立求解Favré平均层流脉动能方程与Favré平均Navier-Stokes方程,开展了高马赫数条件下尖锥边界层转捩位置的数值预测研究.将数值计算结果与可用的实验结果进行对比,结果表明增大单位Reynolds数可使转捩提前发生,攻角效应可使迎风面转捩延迟,背风面转捩提前.     

基于熵限制的Baldwin-Lomax湍流模型

针对Baldwin-Lomax(BL)模型在模拟超声速复杂流动中的不足,提出一种基于熵限制的新型修正方式.分别采用原始的BL模型(BL-origin)、修正后的BL模型(BL-entropy)以及Wilcox k-ω两方程模型对超声速平板、超声速二维压缩拐角以及膨胀-压缩拐角3个典型算例进行计算,结果表明BL-entropy 能够较好地获取长度尺度进而得到均匀合理的涡粘性分布.同时,该修正方式简单高效并且适用于其他与超声速边界层有关的模型,具有较大的发展潜力.     

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

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

XY-SAS模型对于分离流动的性能分析

XY-SAS(Xu-Yan-Scale Adaptive Simulation)模型是一种新的基于平均涡量的一方程SAS类模型,不但解决了分离涡模拟DES(Detached Eddy Simulation)的雷诺平均N-S方程RANS(Reynolds Averaged Navier-Stokes)-大涡模拟LES(Large Eddy Simulation)交界面问题,而且计算量远小于早期SAS.通过对XY-SAS模型的机理分析,认为该模型克服了早期SAS模型预测的湍流能谱在高波数衰减不足的缺陷,并在分析的基础上对该模型进行了改进,取消了其限制器.选用平板边界层和方柱绕流作为算例,进一步研究了该模型对于稳态流动和大分离流动的性能,数值模拟结果同实验值符合良好,证明XY-SAS模型自身完备,是预测稳态流动和大尺度分离流动的有效方法,实现了涡粘性系数根据当地流动状态动态地自适应调整和光滑过渡.     

基于非Favre过滤的大涡模拟

将直接过滤的可压缩Navier-Stokes(NS)方程组作为可压缩湍流大涡模拟方程组.方程组中因过滤产生的高阶相关项用Taylor级数展开近似表达,得到可压缩张量扩散模型.与传统的采用质量加权过滤(Favre过滤)的可压缩大涡模拟方程组相比,由于这种方法没有密度信息的损失,所以能反映真实的密度脉动信息;方程所求量为流场物理量的过滤值,物理意义明确,使大涡模拟结果与直接数值模拟或实验结果更具有可比性.用这种大涡模拟方程组对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验、Baldwin-Lomax(BL)模型及Smagorinsky亚格子模型的结果进行了对比.      

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

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

应用修正湍流模型对旋转通道内换热的模拟

针对旋转光滑矩形通道分别应用针对旋转状态修正的 k-ε、标准k-ω以及提出的针对旋转状态修正的k-ω湍流模型进行流动和换热的数值模拟,通过与实验结论的对比,讨论了采用不同湍流模型对计算结果的影响.计算工况为旋转数Ro=0.24,流体进口雷诺数Re=25000.计算结果表明:采用所提出的针对旋转状态修正的k-ω湍流模型的计算结果要比采用针对旋转状态修正的k-ε以及标准k-ω湍流模型的计算结果更接近实验结论.     

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

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

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

考虑禁飞区规避的预测校正再入制导方法

针对升力式高超声速飞行器再入滑翔侧向制导问题,提出了一种考虑禁飞区规避的预测校正制导方法.纵向制导采用落点误差预测与指令校正相结合的方式,不断更新倾侧角的幅值,实时修正轨迹纵程.侧向制导设计了一种倾侧角反转逻辑的切换机制,利用航向角误差走廊和航向角导向区域控制飞行器的侧向运动.CAV-H再入滑翔飞行器制导仿真实例表明,该方法不依赖于标准再入轨迹,能够导引飞行器规避禁飞区约束.Monte Carlo仿真验证表明,在随机初始扰动和误差存在的情况下,该制导方法具有良好的鲁棒性.      

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的可解能力,显著降低了灰区的影响。