COIL化学非平衡流动的一种解耦方法研究

提出了一种计算氧碘化学激光器化学非平衡流动的解耦方法.求解二维可压缩非定常N-S方程及组分连续方程,计算氧碘化学激光器中的组分混合反应流场.空间离散格式为AUSM up格式,用四步龙格-库塔方法作显式时间推进.假设混合气体为热力学完全而热值非完全气体,化学反应模型采用有限速率反应模型.使用一种"点隐式"方法来处理化学源项的刚性问题.计算结果表明本文方法可以有效解决计算中的刚性问题,与现有隐式算法相比减少了大量的矩阵运算,公式推导简单,易于实现和改进.     

三维化学非平衡无粘流场的数值模拟

本文用ＮＮＤ格式求解了化学非平衡流场，利用点隐式处理组分连续方程的刚性，采用全流场激波捕捉办法和特征理论，成功地进行了物面边界条件的处理，模拟了高马赫数来流条件下捻激波流场。化学反应采用了五组分五反应模型，建立起来的三维流场计算程序可用于复杂外形化学非平衡无粘流场的数值模拟。     

高超声速热化学非平衡钝体绕流数值模拟

对高超声速热化学非平衡钝体绕流进行研究,发展了高温热化学非平衡条件下的高超声速气动热数值模拟方法。采用二阶精度迎风TVD格式的N-S方程有限体积法多块结构网格求解器,考虑了空气5组分及化学非平衡效应的影响,对组分生成源项及振动能源项采用了点隐式的处理,提高了计算效率。对球锥模型和RAM-CII飞行实验等算例进行了数值气动热模拟,分析了不同非平衡模型在热化学非平衡条件下流场的影响。计算结果表明建立的数值模拟方法具有较高的精度。     

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

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

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

电弧喷射推力器化学非平衡数值模拟

为了准确把握电弧喷射推力器工作过程物理机理与特征,采用化学非平衡模型对不同压缩室直径下低功率氮电弧喷射推力器工作情况进行了数值模拟。模型考虑了工作过程中的分解反应、电离反应和复合反应,化学动力学模型为4组分,4反应的有限速率化学反应模型。采用二阶精度NND格式数值求解耦合电磁源项和化学反应源项的N S方程组,采用有限控制容积积分方法求解由麦克斯韦方程组推导出的电磁场方程。数值模拟的结果揭示了推力器内部电弧能量输入作用和高温电离气体的离解电离状况,给出了不同压缩室直径下推力器的推力、比冲和推进效率。结果分析表明,压缩室直径对推力器性能具有较大影响。     

非平衡流解耦方法及其计算激波诱导燃烧的应用验证

介绍并发展了一种改进的化学非平衡流动解耦方法.将流动控制方程与化学反应生成源项解耦处理,组分对流项的矢通量分裂基于流动方程密度对流项的分裂形式给出.对流项采用五阶WENO格式求解,化学反应源项的刚性采用简化的隐式方法进行处理.将该方法应用于预混H₂/O₂和H₂/Air激波诱导燃烧计算,得到的定常及非定常燃烧过程与实验观测相符,振荡频率与实验测量较吻合,结果表明该方法应用于多组分多步反应流计算是可行的.      

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

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

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

甲烷/空气超声速燃烧流动数值模拟

用流动项与化学反应生成源项解耦处理的化学非平衡流动计算方法，从薄层近似三维N-S方程出发，采用ENO差分格式数值模拟了超声速冲压加速器简化模型中，高速甲烷气流从后体多个喷口射入高超声速空气流形成流场，研究了高温异质气体效应和19组分65反应模型的非平衡效应对冲压加速器表面压力分布的影响。计算表明新型解耦方法适合反应机理复杂的碳氢燃料超燃冲压发动机内部流动模拟，为开发应用软件系统打下基础。     

气体运动论BGK格式的翼型绕流数值模拟

以气体运动论BGK(Bhatnagaar-Gross-Krook)格式为基础,通过改进非均匀网格下的通量计算方法,把气体运动论BGK方法应用到二维贴体网格上模拟了可压缩的无黏和黏性流动。在重构阶段,以交界面周围4个网格单元为一组,先把每个网格单元的宏观量(密度、动量和总能量)转换在相对于交界面的坐标系下,然后使用van Leer限制器插值求出交界面的宏观量,在求解阶段再使用BGK理论求出流过交界面的通量。通过上述方法使贴体网格和均匀网格,在计算格式上一致,与其他求解二维贴体网格的气体运动论BGK方法相比,该方法相对简单,计算量减小。结合贴体网格通量计算和气体运动论BGK格式方法理论,分别对亚、跨、超声速下的NACA0012翼型绕流进行了数值模拟和对比研究,所得可压缩无黏流动的数值结果与北大西洋公约组织航天研究与发展咨询组(AGARD)的结果吻合良好。本文还模拟了低雷诺数状态下的可压缩黏性流动,所得结果与文献中的数值计算结果吻合良好,证明了该方法在数值模拟及激波捕捉方面具有鲁棒性,为在复杂流动中的应用提供了一种简单的途径。     

非平衡流计算方法及其模拟激波诱导振荡燃烧

给出一种数值模拟高超声速化学非平衡流动的计算方法，采用流动项与化学反应生成源项解耦处理，即可以有效地解决非平衡流动计算中遇到的刚性问题，还具有公式推导简单和应用方便的优点。从薄层近似N－S方程出发，采用ENO差分格式，数值模拟了超音速H2／Air预混合气体钝体绕流场，较好分辨出流场中激波诱导振荡燃烧非定常流动现象，表明计算方法具有较高的精度。对计算方法的效率和网格数的影响也进行研究。     

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

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

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

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

非结构网格上的TV-van Leer混合迎风格式

为了改进迎风型对流通量格式的激波捕捉计算稳定性,对迎风格式产生数值激波不稳定现象的原因进行分析,并且在分析结果基础上构造新型迎风型格式。为此,首先设计了数值算例,分析非结构网格上迎风格式出现数值激波不稳定现象的主要影响因素。结果表明,在一定条件下,网格和流动状态共同作用导致产生不合理的数值结果,并影响计算的稳定。其中,网格间的动量通量对数值激波不稳定现象影响显著。因此,基于TV格式构造了一种自适应调节动量通量耗散项的混合格式,从而在激波附近引入耗散,抑制扰动的产生和传播。数值计算结果显示新型混合格式在非结构网格离散条件下能够稳定地捕捉激波,同时混合格式保持了TV格式原有的黏性流动计算分辨率。     

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

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

预处理方法在化学非平衡流场数值模拟中的应用

将预处理方法发展应用于化学非平衡流场的数值模拟中。采用有限差分,LU-SGS隐式方法对二维可压化学非平衡控制方程进行了耦合求解,对化学源项和无粘项隐式处理,粘性项则采用显式。通过化学反应剪切层和激波点火的算例表明,所采用的数值方法能够有效求解化学非平衡流场,证明了本算法的可行性,为下一步的工程应用奠定了基础。     

基于粘性非结构网格的三维增升装置N-S方程数值模拟

应用基于粘性非结构网格求解N—S方程方法数值模拟了某型飞机三维增升装置的绕流。在物面附近粘性主导作用区域内，采用三棱柱型非结构网格，在其他流场区域，采用传统四面体非结构网格的混合非结构网格生成技术。湍流模型采用Spalart—Allmaras一方程模型，运用格心格式有限体积法，对控制方程和湍流模型方程进行数值离散。用四步龙格-库塔方法作显式时间推进，采用当地时间步长、隐式残值光顺等加速收敛措施．求解某型飞机全机带增升装置构型的粘性绕流流场，其数值计算与风洞试验结果显示了良好的吻合度。     

二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。     

辐射加热对返回舱气动热环境影响的数值研究

基于高超声速再入飞行器气动热环境预测分析的需要,建立了高温非平衡气体辐射加热对飞行器热环境影响的计算分析手段。采用数值求解化学非平衡N-S方程的方法,对返回舱绕流流场进行模拟,获得高温空气组分质量分数和温度等流场参数分布。基于辐射传输方程,考虑高温气体组分的主要辐射机制,计算分析高温流场气体辐射加热对返回舱热环境的影响。分析表明,在同一飞行弹道条件下,返回舱大底半径尺寸对气动加热的影响较大,在再入热环境严酷区,辐射加热对物面总热流的贡献达30!;产生辐射加热效应的主要机制是高温流场中O和N原子产生连续谱和线状谱以及N2的第一正带系;物面催化效应对辐射加热影响不大。