激波风洞高焓流动及其驻点对流和辐射热流测量

在激波风洞中用氢氧燃烧驱动方法获得了总压１４ＭＰａ,总温高达７２００Ｋ的高超声速高焓平衡流,可以模拟再入飞行速度４至５ｋｍ／ｓ的真实气体效应。本文还介绍了高温气流中驻点对流和辐射传热测量技术及其测量结果。     

烧蚀端头锥模型激波风洞试验研究

介绍了在激波风洞中用压电天平进行1m量级烧蚀端头大模型五分量气动力试验研究情况。试验目的是测量出烧蚀端头模型的Cm0和Cn0,以期为端头烧蚀模型的Cm0和Cn0换算到全尺寸飞行器下提供试验依据。针对设计部门关心的Cm0和Cn0测量,在天平设计上采取了粗短主体结构另配支杆的方式;在风洞试验中采用正负攻角等较为有效的方法。风洞试验主要结果如下:烧蚀变形对锥模型轴向力影响较大,对其余的气动力分量影响值较小。在0°攻角下Cm0为10-4量级,Cn0为10-3量级,配平攻角约为0.028°。由此表明:激波风洞压电天平能够测量该烧蚀端头模型的小Cm0和Cn0。     

管内激波串振荡和壁面脉动压力特性

分析了管内激波串振荡引起的壁面压力脉动特性.采用直联风洞实验方式结合动态压力测量技术获得了矩形管道内不同激波串位置下的壁面压力分布数据,对所得压力脉动幅值、功率谱密度曲线等结果进行分析,探讨激波串流动非定常特点与引起压力脉动的机制.     

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

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

运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

氮气辐射强度的激波管测量与验证

超高速飞行器的热防护设计必须考虑激波层内高温气体发射与吸收的辐射能量,需要有效的辐射加热评估手段。相应飞行条件下的光谱辐射强度地面实验测量是验证数值模型和方法、理解高焓流动的重要手段。基于燃烧驱动激波管,发展辐射强度标定技术,针对富氮气环境,开展高温气体光谱辐射强度的高分辨定量化测试,掌握辐射特征,为数值验证提供基础数据。实验获得了激波速度5.70和6.20km/s条件下的气体光谱辐射强度精细结构,数据表明激波波后的非平衡过程对辐射强度存在很大影响。通过求解耦合化学反应动力学模型的Navier-Stokes方程和辐射特性模型,得到对应实验条件下的流场参数和辐射强度,计算结果和实验数据符合很好,验证了数值模拟方法。     

圆弧激波缓聚点火及后续燃烧传播

采用激波管实验和准一维数值模拟的方法,对预混可燃气体中圆弧汇聚激波的自点火现象及后续燃烧波的传播特性进行研究。其中圆弧汇聚激波由平面运动激波通过精确设计的弧形过渡管段转变得到。研究表明:收缩段中圆弧汇聚激波波后的非均匀梯度环境由激波在平直段、弧形过渡段和扇形收缩段中传播所分别诱导的3个梯度区共同构成。随着圆弧汇聚激波的不断增强,圆弧激波后某处首先形成一个无激波的温和反应区。该反应区逆流锋面的初期运动速度远超Chapman-Jouguet（CJ）爆轰波速,而反应产物区流动则呈现出一定弱爆轰波特征。进一步分析发现,该反应锋面本质上是一种"自发反应波"（spontaneous reaction wave）,而非常规意义上的动力学波,其速度与汇聚激波波后气流点火时间梯度的倒数吻合。而后,反应区的扩张速度很快降至CJ爆轰波速以下,伴随反应锋面附近激波的产生以及激波-火焰复合结构的形成。激波-火焰结构最终加速演变为反向传播的爆轰波。在一定的条件下,由于入射激波转变过程和汇聚所构造的特定点火环境,自发反应波可再次赶超爆轰波,成为新的燃烧波前;而当自发反应波速度再次低于CJ爆轰波速时,它将再次转变为爆轰波;在此过程中,原先的爆轰波阵面蜕变为反应产物中传播的激波。     

椭圆内聚激波面间断化过程的激波动力学分析

针对椭圆内锥几何约束下激波的非均匀汇聚问题,利用高超声速等价原理,将三维定常椭圆内锥激波转化为二维非定常椭圆内收缩运动激波。根据激波动力学原理,发展出一种既能得到非定常激波面演变过程及参数分布,又能沿着激波面追踪扰动传播过程的“波面-扰动追踪法”。该方法不仅具有快速预测激波非均匀汇聚及其演变过程的特点,而且有助于揭示激波面从连续弯曲演变出间断的内在机理。研究表明：初始沿周向强度均匀,而几何形状偏离轴对称的椭圆内聚激波受到自身产生的非均匀“Shock-Compression”扰动,在向中心汇聚的过程中,激波强度非均匀性出现且不断加剧。由于长轴附近的激波面曲率大,激波强度增长得更快。而激波强度的非均匀性会导致扰动的聚集,使得原本连续光滑的激波面出现间断,进而将初始长、短轴附近的激波面分割为强、弱两对激波段。增大长短轴比,椭圆激波的非均匀性演化更快,激波面更早地出现间断。利用“波面-扰动追踪法”对椭圆激波汇聚过程进行分析,为解决三维定常内锥激波的非均匀汇聚问题提供了新的途径。     

超声速飞行器声爆飞行试验与声爆预测评估

声爆问题是超声速民机研制中首要解决的关键问题之一。声爆飞行试验是研究超声速飞行器声爆特性的最直接手段,可为声爆预测方法和低声爆设计技术提供真实可信的验证数据,对新一代低声爆超声速民机设计具有重要意义。中国航空研究院与中国飞行试验研究院发展了基于传感器阵列的地空一体化地面声爆测量技术,开展了国内首次超声速飞机声爆专项测试飞行试验,采集了多组真实大气条件下的声爆实测数据,验证了飞行试验方案的合理性。在飞行航迹正下方测得的声爆波形具有显著相关性,声爆信号头激波峰值相对误差在18%左右,尾激波峰值相对误差在8%左右,声爆持续时间均为0.1 s。对比分析了地面声爆实测数据与数值预测结果,发现：飞行器从测量阵列正上方飞过时,得到的声爆信号基本形态一致、持续时间较接近,声爆信号头激波、机翼前缘激波峰值相对误差小于5%;由于计算模型简化和声爆长距离传播的非线性累积效应等因素,导致声爆信号预测值与实测值在局部特征上有一定差异;后续还需深入研究真实大气环境下的超声速声爆远场传播预测方法。     

CPR方法与高分辨率二阶格式的混合算法

重构修正方法（correction procedure via reconstruction,CPR）具有紧致高效的优点,但对较强激波的捕捉能力还相对较弱,而加权紧致非线性格式（weighted compact nonlinear scheme,WCNS）具有很强的激波捕捉能力。将基于高阶WCNS插值的二阶格式引入到高阶CPR方法中,构造了一种高效高分辨率的混合激波捕捉格式。首先,基于非线性权偏离线性权的程度的激波侦测器侦测出问题单元,并在问题单元附近引入缓冲单元,其余单元则标记为光滑单元。然后,针对问题单元和缓冲单元采用二阶格式计算,光滑单元采用CPR方法计算,构造混合格式。通过对等熵涡问题、含激波的问题以及激波旋涡干扰问题的数值模拟,测试了混合格式的精度、激波捕捉能力和计算效率。数值模拟结果充分说明了该混合格式具有很强的激波捕捉能力,同时在光滑区具有高分辨特性,可以应用于高超声速流动问题的高效数值模拟中。相比于基于高阶WCNS插值的二阶格式,此格式具有更高的计算效率和更高的分辨率。