高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。     

扰动形状对钝头体非对称流动的影响

通过在钝头体头部施加人工扰动块可以得到确定的大攻角下的非对称背涡结构。为了研究扰动块形状对非对称背涡结构的影响,本文在攻角50°、雷诺数Re_D=1.54×105的条件下,利用数值模拟对周向角90°、子午角10°的扰动位置的半球形、D型及方形3种扰动块形状分别进行了研究。研究发现在同一扰动位置,半球形扰动主控下的背涡结构为右涡型,而D型扰动和方形扰动主控下的背涡结构呈现左涡型,且方形扰动主控下的背涡结构的非对称性弱于其他2种扰动主控的非对称背涡。通过分析发现扰动块所引起的微流动直接影响钝头体非对称背涡结构。因此为了更精准地通过施加人工扰动得到确定的非对称背涡结构,应尽量选择形状简单、表面平滑过渡的扰动块形状。      

化学非平衡边界层绕流计算

本文介绍用化学非平衡边界层方程组,数值计算高超声速钝头体绕流流场、组元浓度和电子密度分布情况。计算中考虑了电离空气的七种主要成分以及相互间的化学反应,比较了变熵、普朗特数及Lewis数的影响,讨论了完全催化壁和非催化壁两种情况,计算结果与文献及试验进行了比较,结果令人满意。     

熵层对高超声速二维钝楔气动参数的影响

采用薄激波层理论和流管质量守恒相结合的方法,分析了高超声速二维钝楔边界层外缘熵的分布规律,研究了熵层对边界层外缘密度、马赫数以及壁面气动加热等气动参数的影响.结果表明:边界层外缘熵在倒圆-肩部区下降最为剧烈,熵层对气动参数的影响在高超声速下不可忽略,特别是使转捩区和湍流区的气动加热增加约53.6%.因此,将表面流态控制在层流模式对高超声速飞行器热防护具有重要意义.      

前缘钝化和化学非平衡效应对斜爆震发动机进气道特性的影响

为了探究前缘钝化、化学非平衡效应对斜爆震发动机进气道的工作性能及出口温度边界层分布的影响,采用热完全气体、化学非平衡气体两种模型对不同顶板、唇口前缘钝化半径下斜爆震发动机进气道进行数值模拟,结果表明：相比基准进气道,钝化后进气道上壁面温度边界层厚度较大,下壁面温度边界层厚度较小;在化学非平衡气体模型下,顶板前缘钝化半径(R1)>2mm时进气道的顶板附近和分离区内离解反应较为明显,唇口前缘钝化半径(R2)>2mm时进气道的唇罩、唇口板附近离解反应较为明显;当钝化半径≥4mm时,两种气体模型下进气道出口总压恢复系数和静温的相对变化量绝对值大于0.5%,有必要考虑化学非平衡气体效应对进气道出口性能的影响。     

高焓风洞中钝体近尾流红外辐射测试技术

介绍在JF-10氢氧爆轰驱动高焓激波风洞中开展再人流场红外辐射测量的实验技术。风洞试验状态的驻室总压为19．6MPa，驻室总温为7920K。实验以球头钝锥体为试验模型，测量其近尾流红外辐射能量通量的横向分布。测量采用插入式锑化铟多元红外成像系统，波段范围为2．27-6．0μm。试验数据呈现明显的规律性。试验结果表明：利用这一测量技术能够提供高焓条件下有较高空间分辨率的、较为准确的红外实验数据。     

高焓风洞中带翼钝锥体尾流电子密度的测量与分析

报道了在爆轰驱动高焓激波风洞中开展带尾翼钝锥体电子密度测试的相关研究工作进展.试验气流为4km/s,密度为0.001kg/m3.诊断尾翼对尾流的影响时,为不影响流场并获得足够的空间分辨率采用针状静电探针;实验结果给出带尾翼模型对尾流电子密度影响的定量结果及受影响的空间区域.     

带气膜冷却结构的高超声速平板不同前缘形状下表面传热特性研究

当带红外成像制导系统的飞行器在稠密大气层内做高超声速飞行时,必须采取主动冷却方式防止严重气动加热造成的窗口材料热畸变以及复杂流场造成的气动光学畸变。本文根据成像窗口周围流动具有受高超声速钝头体绕流和气膜冷却结构(即背面为空腔的超声速后台阶)共同作用的特点,在 KD-01高超声速炮风洞中开展了带气膜冷却结构的高超声速平板在不同前缘形状下表面传热特性的试验研究,测量了 Ma8来流条件下喷缝下游表面传热系数,试验获得了2种前缘形状的带气膜冷却结构的高超声速平板喷缝周围瞬态流场 NPLS 图像。通过分析试验数据,得出以下结论:对于带气膜冷却结构(气膜不工作状态)的高超声速平板,模型前缘的形状对喷缝下游区域的表面热流整体分布有明显影响,在钝前缘情形下,表面热流分布接近相同前缘形状的平板边界层为层流状态时的表面热流分布;在尖前缘情形下,表面热流分布则表现出从层流边界层状态向充分发展湍流边界层状态变化的特性;喷缝下游分离和再附区表面传热特性和超声速后台阶流动类似,取决于喷缝上缘处边界层相对厚度。     

高分辨率格式与钝舵绕流数值模拟

 采用Beam-Warming/TVD格式，Roe二阶、三阶以及四阶通量差分分裂格式数值计算雷诺平均Navier-Stokes方程，湍流附加粘性采用Baldwin-Lomax模型计算，针对超声速绕钝圆舵流动进行了数值模拟比较研究，计算给出了流场的详细结构，物面压力分布与实验值符合良好，尤其是Roe三阶格式给出了较好的计算结果，两个超声速回流区以及湍流边界层分离也得到了很好的数值模拟。     

超声速平面钝体流动的快速解法

本文介绍了一种平面流动的快速欧拉方程解法。该方法将原参数非定常欧拉方程组重新组合成以广义称曼变量表示的欧拉方程组，再用二点二步迎风格式离散解。针对钝体流动，本文先建立了动网格下的方程，构造了动网格下的算法。提出了一种简单的激波处理方法。计算结果表明，该方法速度快，稳定性好，对初场不敏感。