壁面催化条件对热环境预测的影响

基于高超声速飞行器热流精确预测的需要,文章采用数值模拟的方法研究了高超声速热化学非平衡流动中壁面催化条件对热环境的影响.针对ram-cⅡ气动外形,使用ausm+up格式耦合lu-sgs方法求解带化学反应源项的多组分ns方程.得到如下结论:有限催化壁结果介于完全催化壁以及非催化壁之间:壁面催化特性的差异主要体现在热流的组分扩散项部分,因此可以通过评估完全催化壁时组分扩散项所占比率大小来预测壁面热流准确值.     

高超声速喷管非平衡黏性流动的数值研究

分别采用平衡气体模型、化学非平衡气体模型和热化学非平衡气体模型,通过求解轴对称Navier-Stokes方程组,数值模拟了法国Marseille高焓激波风洞锥型喷管(H₀=10.3MJ/kg),分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、振动温度、马赫数和组分质量数在轴对称线上的分布与喷管出口附近的速度和温度沿径向分布等结果。计算结果表明:化学反应速度加快,会导致喷管出口温度增加,振动能的冻结会导致喷管出口温度降低。     

碳/酚醛燃气舵热化学烧蚀过程数值研究

为了研究碳/酚醛燃气舵热化学烧蚀问题,在Fluent平台上采用UDF二次开发方法,对碳/酚醛燃气舵二维非定常流固热耦合过程进行了数值仿真研究。对几何建模、边界条件、热解气体逸出过程以及热化学烧蚀导致的边界退移等问题进行了详尽的描述,并选取了合适的计算模型。对不同舵偏角下燃气舵温度分布、热解气体逸出通量以及热化学烧蚀等问题进行了分析研究。结果表明:燃气舵表面化学反应以及热解气体的逸出过程能够有效降低燃气舵壁面以及内部温度,下降温度最高值约为325K;燃气舵前缘一直是热化学烧蚀最严重区域,达到0.5mm,迎风面烧蚀量次之,而背风面几乎不存在热化学烧蚀现象;迎风面热化学烧蚀会随着舵偏角的增大而变得更为严重。     

钝体标模高焓风洞试验和飞行试验相关性的数值分析

在高焓风洞喷管膨胀流动中，会出现组分和振动能量的非平衡冻结现象，给试验数据的分析、外推和使用带来困难。笔者选择有飞行试验数据的钝锥体ELECTRE作为高焓风洞试验的标模，用热化学非平衡Navier-Stokes软件，计算了飞行条件和相应的考虎组分和振动能量的非平衡冻结效应试验条件的模型绕流流场，用双尺度参数ρL和Stanton数，分析试验条件下的热流数据外推飞行条件的问题。研究结果说明：在模型头部区域，保持总焓和双尺度参数ρL不变，热流数据从试验条件外推到飞行条件是可行的；在模型尾部区域，试验条件和习行条件的Stanton数有较大差别，用双尺度参数ρL把热流数据从试验条件外推到飞行条件有较大误差。最后提出了用CFD设计高焓风洞试验条件的思路，并识别真实气体效应显著改变热流分布的高焓风试验能力区域。     

稀薄流热化学非平衡效应的DSMC研究

以钝体绕流驻点线流场为例，应用直接仿真蒙特卡罗法，模拟了五组元混合气体（Ｎ２，Ｏ２，Ｎ，Ｏ，ＮＯ）中的碰撞能松弛，离散反应，交换反应及复合反应，给出了碰撞分子各自由度间的能量交换模式及各类化学反应的选择模式。ＤＳＭＣ仿真结果与宏观粘性激波层方程计算结果，在飞行高度较低时基本相符，表明采用的热化学非平衡非模式是令人满意的，随着飞行高度增高稀薄气体效应增大，二者在激波附近有显著差异，表明高空飞行必须考     

电导率模拟对高超声速MHD控制影响

高温气体电导率是高超声速电磁流动控制数值模拟最重要的参数之一。针对电导率模拟准确性及其对高超声速磁流体控制影响的问题，考虑高超声速飞行器流场中高温气体热化学非平衡效应，采用三维低磁雷诺数磁流体动力学（MHD）数值模拟方法及程序，结合国内外常见电导率处理方法开展典型状态高超声速MHD控制数值模拟，分析电导率模拟对高超声速磁流体流场分布、气动力/热特性的影响。研究表明：磁控热流减缓效果与电导率呈非线性关系，电导率较大时将出现电导率的磁控热饱和现象，其产生的原因可能与化学反应趋向于平衡态存在一定关系；采用定电导率方法，会人为放大磁场洛仑兹力的磁阻力效果，使阻力系数的预测值偏大；不同电导率模型计算得到的电导率分布差异很大，甚至存在数量级的差别，显著影响了磁流体的控制效果，这与电导率模型的适用范围、参数选取原则存在很大关联；对于含多种离解、电离组分的高温气体流动来说，采用基于多电离组分迁移碰撞的电导率模型（本文模型M8），计算与试验一致性最好。     

热化学非平衡流场辐射特性数值计算

设计了耦合辐射的三维热化学非平衡流场计算方法,以适用于任意形状的网格。采用Jameson有限体积法求解耦合辐射源项的三维N-S方程。化学模型包含11个组元,20个化学反应。辐射源项通过直接求解辐射输运方程RTE获得。在空间和方向上分别离散后,利用有限体积法求解不同方向上的辐射输运方程,精确计算了波长为100-1500nm所有主要组元的吸收系数,逐条计算了目前已观测到的空气组元的原子谱线,而对分子的带状谱采用了带模型计算,并最终得到了驻点处辐射密度随波长变化的分布情况。     

高超声速流动壁面催化复合气动加热特性

针对高超声速流动壁面催化特性,计算了不同壁面催化复合系数条件下的球锥驻点热环境。引入了经验证的数值求解Navier Stokes方程的方法,在不同壁温500K～2500K的条件下分别分析了O 2和N 2气体在壁面处的催化复合气动加热特性,得到如下结论：（1） 原子复合放热将提高近壁面温度梯度,改变近壁面组分分布;原子复合放热一部分加热飞行器形成组分扩散热流,一部分加热近壁气体提高近壁温度梯度。（2） 在壁面催化复合系数较小时,原子复合放热主要转化为组分扩散加热,对于不同壁面温度,壁面催化复合系数α<0.1时,单一气体反应组分扩散热流小于总热流的20%。     

高超声速三维热化学非平衡流场的数值计算对比研究

选取国内外典型的高超声速绕流实验,对比研究了不同化学反应模型(Dunn-Kang模型、Park85模型、Park93模型、Park2001模型、Gupta模型)对高超流动计算结果的影响。重点采用Park双温度模型开展热化学非平衡效应的模拟研究,并同时与5组元单温度模型的计算结果进行了对比。研究表明:Park85-7组元双温度模型与实验结果吻合最好,结果最为可靠;而Park85-5组元单温度模型更适用于工程计算。     

高马赫数下横/逆向喷流干扰流场数值研究

横向喷流和逆向喷流广泛用于高超声速飞行器气动力与气动热控制。采用格心型非结构有限体积法求解基于三温度热化学非平衡模型的全Navier-Stokes方程,对高空、高马赫数来流条件下二维圆柱状构型飞行器的喷流干扰流场进行数值模拟,研究了仅存在横向或逆向喷流以及横/逆向喷流同时存在时的复杂流场结构以及喷流降低热流、减阻、改善升力的具体效果。通过控制变量的方法,探究了不同参数（马赫数、静压）的喷流对流场结构及飞行器的气动力、气动热的影响规律。结果表明：在一定条件下,当逆向喷流与横向喷流同时存在时,下游的横向喷流可以影响到上游的逆向喷流流场结构;逆向喷流可以显著减小高超飞行器阻力,并降低头部壁面热流峰值,而横向喷流对高超飞行器的升力特性有一定提高;在横向喷流已用于飞行器姿态控制的情况下,一定条件下可以同时使用逆向喷流,既可以减阻、又可以降低热流峰值,还可以提升升力。