高焓电弧风洞试验热化学非平衡流场数值模拟

针对高焓电弧风洞内部流动的热化学非平衡效应及气体组分和振动能量冻结效应导致的试验数据外推困难问题，基于高焓风洞喷管/试验段/试验模型一体化数值模拟的思路，通过数值求解三维热化学非平衡Navier-Stokes方程，开展了FD-15高焓电弧风洞典型运行状态下流场的数值模拟，与典型试验状态的气动热数据进行了对比验证，研究了试验数据外推飞行条件的方法及有效性问题，分析了提高驻室总压对试验数据外推的影响。研究表明:（1）风洞试验段来流离解度高，热化学非平衡效应及其冻结现象严重；（2）热流校核试验测量数据位于一体化数值模拟的完全催化热流和非催化热流之间，分布合理，验证了计算方法和程序的正确性；（3）试验模型安放位置对模型表面压力和热流存在影响，模型与喷管出口的距离越大，模型表面压力和热流越低；（4）当驻室总压较低时，通过双尺度模拟准则（模拟飞行条件总焓和双尺度参数ρ_∞L）外推热流失效，使用部分模拟准则（模拟飞行条件总焓和驻点压力）外推热流也会出现较大差异，在非催化条件下这一现象更加明显；（5）当驻室总压较高时，使用双尺度模拟准则或部分模拟准则外推飞行条件，产生的热流差异明显减小。     

高超声速飞行器表面温度分布与气动热耦合数值研究

针对高超声速飞行器热防护设计中的高温气体非平衡效应问题和气动热环境精确预测问题,基于流场的非平衡Navier-Stokes方程、表面的能量守恒方程和内部的热传导方程,考虑流场的非平衡效应、表面的热辐射效应、催化效应和烧蚀效应以及热防护层内部的热传导效应,建立了初步的表面温度分布与气动热的耦合计算方法,完善了高超声速飞行器气动物理流场计算软件(AEROPH_Flow)。在表面材料为碳-碳(C-C)条件下,对飞行高度为65km和飞行速度为8,10km/s的半球以及飞行高度为50km和飞行速度为8km/s的球锥模型,开展了表面温度分布与气动热的耦合计算,验证了计算方法和计算软件,分析了表面温度分布对气动热环境的影响。研究结果表明:表面温度分布对气动热的计算结果有较大影响,在气动热环境的预测中,不仅要考虑热化学非平衡效应和表面催化效应的影响,还要考虑表面温度分布的影响,最好是采用表面温度分布与气动热耦合计算的方法,以减小表面温度分布对气动热计算结果的影响。为此,需要发展完善非平衡流场/表面催化和烧蚀/热传导温度场(气/表/固)的计算模型、耦合求解技术和计算软件,实现对高超声速飞行器的真实飞行条件下高温气体非平衡效应和气动热环境的精确模拟。     

火星探测器气动热环境和气动力特性的数值模拟研究

针对探测器进入火星大气过程中的高温气体非平衡效应及其气动热环境和气动力特性的精确预测问题,考虑火星大气高温混合气体的化学反应、束缚电子能激发效应、分子振动能激发效应和非平衡效应,以及不同表面催化特性和表面温度条件的影响,建立了数值模拟热化学非平衡流场的计算方法,并通过与平头圆柱试验模型和MESUR火星探测器的试验与计算数据对比分析,验证了计算方法的可行性与有效性,然后对典型火星探测器的气动热环境和气动力特性进行了数值模拟分析。研究结果表明:1)对于本文的计算状态情况,火星探测器流场中化学非平衡效应非常严重,存在一定的热力学非平衡效应;2)热力学模型、化学反应模型和表面催化特性对气动热计算有着重要影响;3)气体模型对气动力特性的预测影响较大,表面温度条件和表面催化特性对气动力系数的计算影响较小。     

局部催化特性差异对气动热环境影响的计算分析

高温气体非平衡效应及其壁面催化效应对高超声速飞行器气动热环境造成显著影响,是当前高超声速飞行器气动热环境预测和热防护设计的关键问题之一。考虑高温空气离解与电离等化学反应、气体分子热力学激发、流动中的非平衡效应和壁面催化效应,通过数值求解三维热化学非平衡Navier-Stokes方程和壁面处质量、能量平衡关系,完善了高温气体热化学非平衡流场有限催化气动热环境数值计算方法和计算程序,采用典型算例进行了考核验证。在此基础上,开展了不同条件下高超声速飞行器热化学非平衡流场气动热环境数值模拟,分析局部催化特性差异对气动热环境的影响。研究表明：所建立的高超声速飞行器热化学非平衡流场有限催化气动热环境数值计算方法及程序,其数值模拟结果与飞行试验、文献符合;局部催化特性差异会导致热流跳变,其热流跳变量与催化特性差异量、材料分布方式等有关;催化特性差异较大时,局部区域热流可能远远高于飞行器全表面完全催化的热流结果,此时将飞行器在全表面完全催化（FCW）和完全非催化（NCW）条件下的数值模拟结果作为实际飞行过程中表面热流的上、下限这一简化处理方式,是不可取的。     

三维等离子体MHD气动热环境数值模拟

电磁流动控制技术是一个多学科交叉融合的重要研究方向,在高超声速飞行器气动特性优化、气动热环境减缓、边界层转捩和等离子体分布等流动控制方面显示出广阔的应用前景。考虑高超声速飞行器绕流流场中发生的离解、复合、电离和置换等化学反应,气体分子振动能激发以及化学非平衡效应,耦合电磁场作用并基于低磁雷诺数假设,通过数值模拟求解三维非平衡Navier-Stokes流场控制方程和Maxwell电磁场控制方程,建立磁场与三维等离子体流场耦合数值模拟方法及程序,采用典型算例进行考核。在此基础上,开展不同条件下磁场对再入三维等离子体流场以及气动热环境影响分析。研究表明:建立的高超声速飞行器的等离子体流场与磁场耦合计算方法及程序,其数值模拟结果与文献符合,外加磁场使飞行器头部弓形激波外推,磁场强度越强,激波面外推距离越大;不同磁场强度环境下,流场中温度峰值大小略有变化,变化幅度较小;磁场对绝大部分区域的热流有减缓作用,作用的大小与飞行高度、马赫数以及磁场的配置紧密相关;当前的计算条件下,飞行的高度越高,磁场的作用越明显。     

高温气体效应对高超声速磁流体控制的影响

高温气体效应会严重影响高温气体流场的流动特性，进而影响高超声速磁流体控制效率。基于低磁雷诺数假设，通过耦合求解带电磁源项的三维Navier-Stokes流场控制方程和电场泊松方程，开展完全气体模型、平衡气体模型、化学非平衡气体模型、热化学非平衡气体模型等条件下的高超声速磁流体控制数值模拟，分析气体模型对磁流体控制的影响，研究高温气体各种非平衡效应及焦耳热振动能量配比等对高超声速磁流体控制的影响规律。研究表明：化学非平衡效应对高超声速磁流体控制影响显著，采用化学非平衡气体模型模拟得到的磁控增阻特性介于完全气体模型和平衡气体模型之间，平衡气体和完全气体模型磁控热流变化的定性规律，与非平衡气体模型模拟结果差异很大；热力学非平衡效应对高超声速磁流体控制的影响，与焦耳热振动能量作用比率紧密相关，随该配比增大，磁场增阻效果由67%降到约12%；高温气体效应会极大地降低磁控增阻效果，会明显地增强部分表面区域的磁控热流减缓效果，要准确数值模拟高超声速磁流体控制，必须有效地考虑化学和热力学非平衡效应，同时选用接近实际情况的焦耳热振动能量配比。     

火星探测器表面材料催化特性数值模拟研究

针对火星探测器高超声速进入过程中的表面材料催化作用及其对气动热环境影响这一问题，建立了变壁面温度的火星大气表面材料催化作用模型，并基于火星大气物理化学模型和求解三维热化学非平衡N-S方程的数值方法，对典型火星探测器防热大底进行了数值模拟，获得了不同催化特性下的高超声速非平衡流场和气动热数据，分析了表面材料催化特性对气动热环境影响的规律性。研究结果表明：表面催化特性对壁面附近组分分布影响很大，催化反应进程主要受O原子浓度限制；有限催化热流随催化效率增大而增大，完全催化峰值热流比催化效率为1的有限催化峰值热流高25%~64%；表面温度随催化特性的变化规律与热流变化规律类似。有限催化模型能根据表面材料的催化特性精细化预测表面热流和温度，为防热设计提供更精确合理的参考标准。     

化学模型对数值模拟等离子体流动的影响研究

针对典型再入飞行试验条件,基于求解化学非平衡N-S方程的数值方法,对再入体非平衡等离子体绕流进行数值模拟,比较分析了两种典型化学模型对电子数密度数值模拟结果的影响。研究发现,7组分化学模型的电子数密度模拟结果具有较好的一致性,而11组分化学模型的数值模拟结果则存在较大差异,产生差异的主要机制在于不同模型中氧原子的电离效应,总体上Park模型的数值模拟结果与飞行测量结果具有较好的一致性。     

磁控热防护系统在天地往返运载器上的应用仿真

磁控热防护技术在高超声速领域显现出广泛的应用前景。考虑高超声速流动磁流体力学控制涉及的等离子体生成机制、多电离组分导电机理以及电磁流动能量/动量输运机制,通过耦合求解电磁场泊松方程和带电磁源项的高温热化学非平衡流动控制方程组,搭建了高超声速磁控热防护数值模拟平台。结合美国航天飞机"哥伦比亚"号（OV-102）近似外形和5种磁场配置方案,较为系统地开展了磁控热防护系统在高超声速"滑翔返回式"天地往返运载器上的应用仿真研究。结果表明：搭建的磁控热防护仿真平台具备偶极子磁场、均匀磁场、螺线管磁场及多个磁场组合条件下复杂外形飞行器气动热环境数值模拟能力,其校验结果与文献或飞行试验数据符合较好;采用合适的磁场配置能有效降低航天飞机的表面热流,显著改善了航天飞机的气动热环境,典型状态的表面热流下降25%以上;局部磁场方向与流动方向的夹角,在一定程度上决定了洛伦兹力的强度和方向,对磁控效果的影响明显。