应用γ-Reθ湍流模型和PSE方法研究高超声速流动

采用考虑转捩的γ-Reθ湍流模型研究了高超声速复杂流动,结合PSE(parabolic stability equa-tion)稳定性分析方法和e-N方法,通过给出边界层最不稳定点的参数,用以控制γ-Reθ的转捩经验关联式,减少对经验参数的依赖.计算了锥体高超声速流场,得到的壁面换热量和边界层转捩位置与实验吻合较好.计算了平板激波/边界层干扰和高超声速拐角压缩流场,得到了准确的波系结构、压力分布和分离区的大小.结果表明,该方法与传统湍流模型相比具有较大的优越性.     

基于γ-Reθt转捩模型的翼型气动特性计算精度改进研究

将结合SSTk-ω湍流模型的γ-Reθt转捩模型与二维雷诺平均NS（RANS）方程数值求解相耦合,发展了考虑转捩影响的翼型气动特性计算方法。得到的计算结果与没有考虑转捩影响的计算结果的比较表明,与实验更加吻合,计算精度得到了明显提高。     

高超声速锥柱裙模型边界层转捩的弹道靶实验

为研究高超声速边界层转捩现象、给边界层计算提供可靠的对比数据,在中国空气动力研究与发展中心超高速弹道靶上开展了锥柱裙模型高超声速边界层转捩的自由飞实验。所采用的锥柱裙模型全长105mm,飞行速度1.94km/s(犕犪=5.65),单位雷诺数4.32×107~1.20×108 m-1。使用激光阴影成像技术,获得了锥柱裙模型边界层转捩和湍流边界层发展的图像,测得的湍流边界层厚度在0.6~2.2mm 之间,湍流涡的流向尺寸与边界层厚度的比值介于0.3~0.8之间且沿流向呈下降趋势。实验结果表明:弹道靶实验能够获得给定飞行环境下的高超声速边界层转捩图像,从图像中可以清晰判断转捩位置或区域、测量边界层厚度和分析湍流涡的尺寸。     

转捩模式与转捩准则预测高超声速边界层流动

对原始的k-ω-γ转捩模式和"层流+转捩准则"模型进行了改进,在2种方法中分别增加了横流模态时间尺度和横流转捩准则用于预测横流失稳诱导转捩。通过对网格预处理可并行计算获得边界层外缘信息以及边界层内横流速度。采用不同雷诺数条件下的0°攻角尖锥以及HIFiRE-5外形对2种方法预测高超声速边界层转捩的性能进行了比较分析。研究结果表明,2种方法均能正确反映高超声速边界层转捩起始位置和转捩区长度随雷诺数的变化趋势,但不能捕捉转捩区热流峰值;"层流+转捩准则"模型计算得到的传热系数在全湍流区较k-ω-γ转捩模式偏高。对于同时存在流向不稳定和横流不稳定的HIFiRE-5外形,改进的k-ω-γ转捩模式和改进的"层流+转捩准则"模型相比于原始的模型均能更加准确地预测中心线两侧横流失稳诱导形成的转捩;对于中心线附近因速度剖面拐点引起的边界层转捩,"层流+转捩准则"模型由于与边界层厚度相关,预测得到的转捩位置较试验结果靠前,k-ω-γ转捩模式与试验结果吻合很好。     

Menter转捩模型在高超声速流动模拟中的改进及验证

通过压力梯度参数和湍流普朗特数的修正,对Menter等构造的低速经验关系式转捩模型进行了改进,将基于局部信息的Menter转捩模型用于高超声速流动转捩数值模拟。在多个高超声速尖锥转捩流动模拟中,对改进后模型的可靠性进行了检验,算例考核了不同壁面温度、来流湍流度、流向压力梯度等多种流动条件。结果表明,数值计算和风洞试验的壁面温度、温度恢复因子和热流等符合得较好,修正后模型较好地模拟了转捩起始位置和转捩区长度,在高超声速边界层转捩预测中具有一定潜力。     

横流不稳定性转捩预测模型

由于Langtry和Menter提出的γ-Re_θt边界层转捩模型只能预测流向的边界层转捩现象,因此继续改进该转捩模型使其具有横流不稳定性转捩的预测能力显得非常必要。通过对经典Falkner-Skan-Cooke (FSC)三维边界层相似解的理论分析和数值求解,结合Thwaites压力梯度因子与当地后掠角构建的函数关系来求解复杂构型的当地Hartree压力梯度因子β_H以及形状因子H₁₂,采用由试验数据标定的C1准则关系式获得横流转捩位移厚度雷诺数,从而建立能够对复杂构型进行横流不稳定性转捩预测的转捩判据。应用所建模型对30°前缘后掠角的ONERA-M6机翼和变前缘后掠角的DLR-F5机翼以及标准6:1椭球标模进行了横流不稳定转捩数值模拟,计算结果显示转捩位置均与试验数据吻合较好,证明了该模型的合理性和实用性。     

超高负荷低压涡轮非定常特性实验

针对某后加载超高负荷低压涡轮叶型的定常与非定常气动特性进行了实验研究.结果表明:由于上游尾迹的作用,在低来流雷诺数与来流湍流度下,尾迹诱导转捩叠加自然转捩可以抑制分离、降低叶型损失.在一个尾迹通过周期,由于尾迹诱导转捩的作用,附面层分离周期性地放大、缩小.同时在吸力面速度峰值点下游观测到了尾迹被割裂为两部分的现象,形成了一个主尾迹与一个副尾迹.副尾迹同样可以诱导转捩,但强度较低,对分离抑制效果有限.      

基于准则的大展弦比飞翼气动设计

从设计实际出发,为切实提高气动性能,开展了大展弦比飞翼无人机(UAV)的气动设计及分析研究.在设计分析过程中,依据飞翼无人机的特征,提出了气动设计准则;基于设计准则,采用更新设计的策略,结合变可信度数值模拟、代理模型优化方法构建了优化设计框架;针对飞翼无人机开展了参数化表达、无限插值网格自动生成以及多轮更新优化,得到了优化推荐构型;应用γ-Re_θt转捩模型方法对优化构型的气动性能进行了细致地验证分析.研究结果表明:通过气动设计,飞翼无人机设计构型很好地契合了设计准则,其巡航升阻比相比最初的原始构型提高了14%,γ-Re_θt转捩模型能较细致地分析大展弦比飞翼的流动特征.      

飞行高度对高超声速钝锥边界层稳定性及转捩的影响

研究了飞行高度对高超声速钝锥边界层稳定性及转捩的影响。通过求解三维可压缩Navier-Stokes方程计算了来流Ma=6,半锥角为7°的钝锥在飞行高度20~40 km条件下的基本流场,利用线性稳定性理论（LST）研究了飞行高度对钝锥边界层流动稳定性的影响,最后采用e N方法进行了转捩预测。研究发现,随着飞行高度的增加,流向不稳定N s值和横流不稳定N cf 值均减小,由横流不稳定性引起的圆锥表面大部分区域转捩逐渐转变为流向扰动引起迎风面转捩横流扰动引起背风面转捩,继而横流扰动消失,流向不稳定波引起迎风面转捩。     

基于离散伴随的层流翼优化设计方法

层流技术是未来发展绿色航空的核心技术,可处理大规模设计变量的高效、鲁棒的优化设计方法是推动层流技术工程应用转化的关键环节之一。基于高可信度RANS求解器,结合准三维层流边界层方程、Drela-Giles和C1准则,建立了可同时捕捉Tollmien-Schlichting(TS)和Crossflow(CF)不稳定性的转捩预测方法,与典型风洞和飞行试验对比,验证了该方法的可靠性。通过严格物面插值、获得准确的间歇因子函数以剔除气动力计算引入的数值噪声,精确推导了考虑转捩的耦合伴随方程,并结合无矩阵存储技术、链式求导法则、反向混合自动微分及Coupled Krylov(CK)算法发展了耦合伴随方程高效求解方法,最终建立了基于离散伴随的层流翼梯度优化设计方法。针对典型客机翼身组合体构型的气动优化显著推迟了转捩,获得了10.48%的减阻收益。优化结果表明,建立的层流翼梯度优化方法能够有效处理多种转捩机制并存的复杂三维层流翼优化问题。