基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

主流攻角对双射流孔气膜冷却特性的影响

应用压力敏感漆技术,在平板上测量了不同主流攻角(i=-30°,-20°,-10°,0°,10°,20°,30°)下双射流孔的气膜冷却效率,并利用计算流体动力学(CFD)计算得到的流场对气膜冷却效率的规律进行了解析。所研究的双射流孔结构的孔间无量纲横向距离为0.5,孔间无量纲流向距离为3;射流与主流密度比为1.0,吹风比分别为0.5、1.0、1.5、2.0。结果表明小的主流攻角(i=-10°,0°,10°)下,流场中存在反肾型涡对或挤压作用,气膜层与壁面贴附良好,气膜冷却效率最高;大正值攻角(i=20°,30°)下,虽然气膜覆盖面积大,但反肾型涡对退化,气膜冷却效率下降;大负值攻角(i=-20°,-30°)下,流场中有肾型涡对,且气膜横向覆盖受限,气膜冷却效率最低。     

燃烧室结构高温变形非接触测量试验研究

为获得地面试验中主动冷却燃烧室在高温状态下的结构总变形和局部应变,基于数字图像相关法开展了非接触测量试验研究。分别采用加速鲁棒特征算法和改进的圆柱曲面数字图像相关法完成了燃烧室结构总变形和局部应变测量及分析。测量结果表明,燃烧室稳定工作条件下,结构轴向总变形约4.8 mm、局部平均应变约0.0049。测量结果与工程粗略估算结果比较吻合,说明所使用的非接触测量方法有效,能够支撑燃烧室热结构设计,测量数据能够用于三维结构强度数值计算的验证。     

低速压气机静叶非定常Clocking效应数值模拟

采用三维粘性非定常数值模拟方法研究了重复级低速压气机内部非定常流动,其中包括静叶Clocking(时序)效应以及尾迹-边界层、尾迹-尾迹和尾迹-泄漏流的相互作用.对静叶时序效应的结果分析显示,当两级静叶相对周向位置改变时最大效率变化只有0.1%,表明低速压气机中静叶时序效应不明显.另外,非定常结果流动图画显示,在非定常环境下上游周期性尾迹的通过对下游叶片边界层、尾迹以及叶尖泄漏流的发展都有一定影响.     

大扩张角涡轮过渡段非定常流动数值仿真

为进一步减小涡轮过渡段流动损失,深入了解涡轮过渡段中的非定常损失机理,开展了大扩张角过渡段研究。在过渡段的非定常流动机理研究中,过渡段进口流场的最显著特点是:转子泄漏涡、通道涡和尾迹。采用数值方法对大扩张角涡轮过渡段进行3维非定常数值仿真。结果表明:支板尾缘部分的静压波动小于支板前缘部分的;高压涡轮静子尾迹被转子切割后进入转子通道中向下游传播并在过渡段内形成尾迹通道,尾迹在过渡段内的时空演化是过渡段内损失的主要来源;过渡段支板表面负荷分布发生明显的周期性变化,支板表面承受较强的非定常力,在过渡段设计中必须考虑。     

螺旋槽端面气膜密封结构高温特性研究

为了研究螺旋槽端面气膜密封结构在高温下的密封性能,建立了高温密封分析数学模型,研究了螺旋槽气膜密封的气膜温度、压力以及端面变形分布规律,在此基础上探讨了螺旋槽气膜密封的热变形机理,并进一步分析了不同密封压力、转速、环境温度下热效应对开启力和泄漏率的影响。结果表明:在高压、高速条件下,热效应使端面形成发散间隙,导致开启力减小,泄漏率增加;在低压、高速条件下,热效应使端面形成收敛间隙,导致开启力及泄漏率增大。对于螺旋槽端面气膜密封结构,环境温度的升高对端面变形的影响不明显,且环境温度从300 K升至550 K时,考虑端面热效应的开启力减小4%,泄漏率减少36%。     

对转螺旋桨流场气动干扰数值模拟

基于结构网格动态面搭接技术,通过求解三维非定常雷诺平均纳维-斯托克斯(RANS)方程,对单独螺旋桨(SRP)滑流流场进行了数值模拟,结果与风洞试验数据吻合良好,在此基础上深入研究了对转螺旋桨（CRP）的非定常滑流效应及前后桨之间的气动干扰现象。计算结果表明:所采用的非定常RANS方法能够良好刻画对转螺旋桨滑流流场的发展变化特征,适用于分析前后桨气动干扰特性。前后桨桨尖涡的复杂相互作用是导致对转螺旋桨气动干扰的直接原因,后桨桨尖涡对于前桨桨尖涡具有一定耗散作用。在气动干扰影响下,对转螺旋桨气动力会产生周期性波动,波动周期与螺旋桨桨叶数相关, 6×6对转螺旋桨在1个旋转周期内产生了12次波动。气流穿过对转螺旋桨会发生两次加速,因此在相同工况下,对转螺旋桨的拉力系数和功率系数相比于单独螺旋桨分别增大约1倍,效率提升约1.5%。      

摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。      

纵向波纹隔热屏气膜冷却特性实验

针对加力燃烧室纵向波纹隔热屏气膜冷却效果开展了细致的实验研究,利用红外热像仪测量了隔热屏壁面的温度分布,分析了隔热屏板型、吹风比、开孔率等参数对气膜冷却效率的影响。实验中板型选取了平板和纵向波纹隔热屏,吹风比变化范围是0.5~3.0,开孔率变化范围是1.4%~3.7%。结果表明:相比于平板隔热屏的气膜冷却效率沿程逐渐增加,纵向波纹隔热屏的气膜冷却效率随波纹板的起伏而起伏且大于平板隔热屏;随着吹风比的增加气膜冷却效率逐渐加大,在吹风比为3.0时达到最大值;气膜冷却效率在波峰处低,波谷处高,整体上随波纹板的起伏而波动,吹风比越小,气膜冷却效率随波纹板的起伏变化越明显;高吹风比(吹风比为2.0~3.0)下,气膜冷却效率沿程变化与增幅较为缓慢;整体上,随着开孔率的增加气膜冷却效率逐渐加大,小开孔率(开孔率为1.4%、2.7%)下的气膜冷却效率相差不大,但在次流背风侧,开孔率小的气膜冷却效率要小于开孔率大的气膜冷却效率。      

高速进气道出口快速撤锥再起动特性研究

为了探寻出口快速撤锥过程对高速进气道起动性能的影响,对二元高速进气道出口快速撤锥过程开展了非定常数值仿真研究。对比了进气道在一个喘振周期中从不同状态点处进行快速撤锥所表现出的起动性能的差异。研究表明,当来流马赫数低于进气道的加速自起动马赫数时,对于高速进气道因出口堵锥节流所引起的不起动流态(喘振),从不同喘振状态开始撤锥,进气道表现出不同的起动特性,而进气道最终能否顺利建立起动流态与快速撤锥时机的选取有关。