电弧风洞高焓流场模拟与气动热校核研究

为了在电弧风洞更好地开展高超声速飞行器防热材料或结构热考核试验,基于有限体积离散,建立了针对电弧风洞高焓试验状态的多组分热化学非平衡流场数值模拟方法。针对中国空气动力研究与发展中心20MW电弧风洞不同试验模型的高焓流场进行了模拟,获得了试验状态的流场特性和模型表面热流分布,通过与试验测量值比较验证了计算方法。研究发现在喷管出口和试验模型之间的轴向距离很近的情况下,可以采用喷管流场和模型绕流分开模拟。通过对比数值模拟方法中的热化学模型,表明采用Gupta 7组分和5组分空气化学反应数据获得的模型表面热流非常接近,同时相比单温度模型结果,双温度模型结果与试验测量结果更接近。比较了多个状态条件下计算和试验测量获得的试验模型表面热流,发现二者相差都在15%以内,验证了建立的数值方法模拟该风洞高焓流场的可靠性。     

来流条件对热流组分扩散项影响效应分析

使用数值模拟的方法研究了高超声速热化学非平衡流动中不同来流条件对热流组分扩散项的影响。以RAM-CII飞行试验外形为例,使用AUSM+-up格式耦合LU-SGS方法求解带化学反应源项的多组分NS方程。得到如下结论:在非平衡热环境数值模拟研究中,完全催化壁时组分扩散项热流所占比率表征了壁面催化效应的强弱;在同一高度下,随着马赫数升高热流组分扩散项比例越来越显著,而马赫数相同时,随着高度增高热流组分扩散项所占比例越来越小,主要原因在于当地化学反应进行程度不同。     

热与化学耦合非平衡对热流和电子密度的影响

采用单、双温度模型两种粘性激波层方程组对弹头驻点区非平衡流场作了数值模拟。计算结果表明，在算例条件下热非平衡与化学非平衡的耦合效应对壁面热流和电子密度有明显影响。文中给出了完全催化壁和非催化壁两种情况下以及采用不同振动松弛时间公式时上述影响的变化规律。     

高焓激波风洞喷管流场非平衡特性研究

高焓激波风洞是开展高超声速流动研究的重要地面模拟设备,但其产生的高焓气流在喷管中的膨胀过程是一种典型的热化学非平衡流动,试验段特征参数通过直接实验测量难以完全确定。本文通过求解耦合双温度模型的轴对称Navier-Stokes方程,研究了高焓激波风洞中典型状态下气流的热化学非平衡流动特性,分析了焓值对非平衡特性的影响规律。结果表明,喷管出口自由流均匀区域达到出口截面直径的75%以上,能够为实验提供足够的空间;喷管出口自由流处于热化学非平衡状态,在喷管喉道后约1/5喷管长度处气流即已处于冻结流状态,组分浓度和振动温度随气流流动基本不变;焓值在8.4MJ/kg~19.5MJ/kg之间变化时,非平衡程度随着焓值的增加而增强,但是低焓值时非平衡程度的增强更加剧烈。     

面向飞行试验的多源气动数据智能融合方法

风洞试验和飞行试验是飞行器研制过程中进行气动性能分析与优化设计的重要手段，然而，在高超声速飞行条件下，真实气体效应、黏性干扰效应和尺度效应的复杂变化给气动数据精准预测带来巨大挑战。为了提升天地气动数据一致性，针对某外形飞行试验数据开展了典型对象的天地气动数据融合方法研究。结合数据挖掘的随机森林方法，本文提出了一种面向飞行试验的数据融合框架，通过引入地面风洞试验气动数据，实现了对复杂输入参数的特征分析与特征排序，进一步对不同飞行时刻下飞行试验的气动数据开展了交叉验证。结果表明随机森林的机器学习框架对风洞-飞行试验数据关联具有较好的预测与外推能力，可以有效提升气动数据预测精度，相关研究为复杂环境下气动数据多源融合提供了思路。     

驻点壁面催化速率常数确定的研究

以平衡流动作为热环境估算的依据，提出了用数值求解非平衡Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和实验测量热流值确定模型表面材料催化速率常数的方法。用５组分１７个化学反应Ｄｕｎｎ－Ｋａｎｇ空气化学模型和轴对称热化学非平衡Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，对激波管中球头和平头圆柱模型绕流流场进行了数值模拟，给出了驻点热流随催化速率常数变化的分布，并根据激波管实验测量的热流值确定了表面材料Ｐｔ、ＳｉＯ２、Ｎｉ和某种     

高焓流场球头外形气动热试验研究

再入地球大气层时,飞行器的再入速度极高,面临严重的气动加热问题。为了研究高焓流动导致的热化学非平衡现象,在高焓膨胀风洞FD-14X中开展了球头外形的热流测量试验以及CFD仿真模拟。FD-14X为中国空气动力研究与发展中心新建成的高焓膨胀风洞,速度模拟能力达到第二宇宙速度,总温模拟能力超过10 000 K,能够产生总焓70 MJ/kg的试验气体。试验来流总焓16.9～63.5 MJ/kg,球头直径20～50 mm,流场采用自发光拍照,同时CFD仿真采用Park双温非平衡模型计算球头绕流流场。试验与仿真结果表明：来流总焓大于5 MJ/kg时,球头绕流场存在显著的热化学非平衡现象;304钢模型壁面在来流总焓小于20 MJ/kg时表现为非催化壁面特性,在来流总焓大于30 MJ/kg时表现为催化壁面特性;当球头表面镀氧化锆膜、来流总焓49.5 MJ/kg时,球头壁面表现为非催化壁面特性。     

热化学非平衡来流条件下热化学模型影响研究

为研究热化学非平衡来流条件下热化学模型等计算设定对斜激波压缩流动计算结果的影响,针对尖劈构型和相应的前缘钝化构型的高焓激波风洞实验,采用多种计算设定开展详细的数值模拟研究。计算结果表明,计算采用不同热化学模型,以及来流设定为振动冻结/平衡/非平衡状态,会导致斜激波激波角等参数存在一定差别,其中激波角差别可达约2%。当来流速度一定时,过斜激波后分子内能增量在平动转动能和振动能上的分配方式的差别决定了激波角的差别。前缘钝化情形下,采用不同计算设定所得激波角之间的关系和尖前缘构型的规律一致;但是,采用不同计算设定所得斜激波到壁面距离之间的关系和尖前缘构型的规律有差别,这源于钝化前缘的激波脱体距离的影响。对于自由来流下的斜激波压缩流动问题,若考虑了分子振动能激发但未考虑热力学非平衡(例如热完全气体模型、考虑空气反应的单温度模型等),就斜激波激波角等参数而言,计算误差比量热完全气体模型计算误差更大。     

不同热化学非平衡模型对高超声速喷管流场影响的数值分析

采用轴对称NS方程,数值研究不同热化学非平衡模型对高超声速喷管流场的影响,包括:(1)不同组元数的高温空气模型(5组元、7组元1、1组元)的比较;(2)热力非平衡(双温度)的化学动力过程与热力平衡(单温度)的化学动力过程的比较。计算结果表明,高焓风洞实验条件下喷管流场处于热力和化学都是非平衡的状态。在计算条件下,数值模拟以采用7组元或11组元的热化学非平衡模型为宜。     

气体模型对高超声速再入钝体气动参数计算影响的研究

本文采用完全气体模型、振动激发气体模型、平衡气体模型、一温度非平衡气体模型、两温度非平衡气体模型和三温度非平衡气体模型进行了钝体高超声速绕流流场的数值计算。分析了各种气体模型对激波脱体距离、壁面热流、温度分布和密度分布等的影响。结果表明：采用两温度（或三温度）非平衡气体模型计算的激波脱体距离更接近实验情况；采用平衡气体模型计算的壁面热流量高，一温度非平衡气体模型完全非催化的壁面热流最低，而完全气体模型的壁面热流在两者之间。最后给出了各种气体模型的密度等值线比较图。