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

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

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

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

壁面催化效应对高超声速气动热影响研究

高超声速飞行器壁面催化效应会导致激波层中原子在壁面处复合释热,加剧周围气动热环境。针对高超声速流动壁面催化特性,选择不同飞行马赫数及高度条件,采用完全催化和非催化两种条件对球锥模型壁面热流率进行数值模拟计算,研究壁面催化效应对气动热的影响规律。结果表明,固定飞行高度时,壁面催化效应对气动热的影响随马赫数增加而加强,Ma=25条件下驻点处完全催化与非催化热流比值高达1.92;固定飞行马赫数时,在50km高空以上壁面催化效应对气动热的影响随高度增加而减弱;壁面催化效应不仅会影响壁面附近的流场特性及组分分布状态,而且对整个激波层都有一定的影响作用。     

滑移流区化学非平衡流中气动热环境的数值模拟

针对滑移流区化学非平衡流中的气动热环境预测问题，采用多组分化学非平衡纳维-斯托克斯（N-S）方程对轨道验证飞行器（OREX）进行数值模拟，对比了有限催化和非催化、滑移和无滑移壁面条件下的气动加热，研究了壁面催化和滑移效应对气动热的影响规律，并分析了影响气动热的主要机制。结果表明：壁面催化和滑移效应均会影响流场物理量分布，高度为92.82 km时对激波脱体距离的影响较为明显；随着高度减小，有限催化与非催化壁面间的热流偏差增大，而滑移与无滑移壁面间的热流偏差减小；高度低于92.82 km时有限催化壁面计算得到的驻点热流值与飞行数据吻合较好，偏差均在11%以内；机制分析发现，催化效应对壁面附近的组分分布影响较大，滑移边界条件中的温度跳跃条件对壁面附近的温度分布影响显著。     

化学动力学模型对H2/AIR超燃模拟的影响

为了研究超燃流动中不同化学反应模型的影响，提出流动项与化学反应生成源项解耦处理的化学非平衡流动计算方法，从薄层近似N－S方程出发，采用ENO差分格式，数值模拟了超声速H2／Air 预混合气体中钝体激波诱导的振荡燃烧流场，针对飞行速度大于和接近爆轰速度两种不同流动状态，比较Evans-Schexnayder提出的7组分8反应机理（ES模型），Jachimowski提出的13组分32反应机理（J模型）及其简化为9组分19反应机理（J－1模型）对流场参数的影响，计算表明反应机理引起压力场变化较小，对温度等参数影响较大，对于飞行马赫数低于7的问题，取不包括N原子反应的Jachimowski的9组分19反应机理较合适。     

有化学反应的气体的非平衡数值模拟

利用三温度动态模型,对具有化学反应的强正激波后的高温气体的非平衡流动进行了数值模拟,在较宽的马赫数范围内（6～26）,计算出了各种气动参数及气体组分在非平衡区随时间的变化历程,并与有关文献进行了比较。     

火星进入器高超声速化学非平衡特性数值计算研究

火星探测活动是当前国际深空探测领域的热点,而探测器的进入、下降和着陆阶段的气动特性预测是关系其布局和热防护系统设计的关键问题之一。针对火星进入器高超声速进入过程,发展了火星环境的高超声速化学非平衡流动数值模拟方法,建立了进入段气动特性的有效预测技术。以探路者号火星进入器防热大底为研究对象,采用5组分6化学反应模型开展典型弹道点上的化学反应流动数值模拟。结果表明:CO_2在激波层内发生较大程度的离解反应,火星进入器高超声速化学非平衡特性显著。因化学反应作用强烈,激波层内热传导所致热流量相对较低,而组分扩散造成的热流量占一定比重,计算所得的驻点总加热量与经验公式吻合良好。火星环境的化学非平衡特性的有效预测可为我国未来火星探测任务提供技术支持。     

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

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

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

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

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

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

高温化学非平衡湍流边界层脉动量象限分析

高超声速飞行器在较低空域以极高马赫数飞行时，表面会同时存在湍流与化学非平衡流动，但目前针对此类高温化学非平衡湍流边界层流动特性的研究工作还比较有限，对不同湍流特征的主导流动机制的认识还有待于进一步深入。选取高超声速楔形体头部斜激波后的流动状态，设置3种不同的壁面温度，通过直接数值模拟对比了不同壁温条件下的边界层参数分布特性，并采用象限分解技术分析了边界层不同象限流动事件对雷诺剪切应力、湍流热流、湍流质量扩散的贡献。结果显示：在整个边界层中上抛和下扫运动对雷诺剪切应力的贡献占优；冷壁效应会使得流向和法向湍流热流通量的主导流动事件在温度峰值两侧发生改变。O原子组分流向湍流组分扩散主要受到高质量分数流体慢速运动事件和低质量分数流体快速运动事件的影响，而法向湍流组分扩散则主要受到高质量分数流体向上运动事件和低质量分数流体向下运动事件的影响。     

气动热数值模拟中壁面热流后处理方法

针对气动热数值模拟中壁面热流后处理方法进行了系统的数值研究,提出了一种有别于常规微分法和积分法的壁面热流后处理方法.通过3种算例进行了对比分析,结果表明:即使在流场参数相同情况下,采用不同的壁面热流后处理方法计算的热流结果也会有较大差异,且表现出不同的网格依赖性.微分法对网格的依赖性较大,单纯提高差分精度不一定能够减弱网格依赖性和提高计算精度,还与近壁处网格分布密切相关;积分法从能量平衡积分方程出发,极大地减弱了网格依赖性,具有较高的计算精度.所提出方法具有与积分法相当的网格依赖性和计算精度,尤其在驻点区域,其网格依赖性更弱于积分法,且计算简单,提高了计算效率.      

燃烧室中辐射热流分布的蒙特卡罗计算

本文用蒙特卡罗法计算了燃烧室壁面辐射热流密度分布。计算中考虑了燃气参数如:压力、温度、燃气组分及燃气的吸收和发射特性沿燃烧室轴向、径向和周向的变化。详细地说明了计算过程并进行了算例计算。计算所得热流密度值与分布规律是合理的。      

心形孔气膜冷却特性的数值模拟

为进一步提高航空发动机热端部件的冷却效率,提出了心形气膜冷却孔结构,利用数值模拟分析心形孔的流场特性和冷却特性,并通过与常规圆形孔计算结果的对比,揭示心形气膜孔强化冷却的物理机制.计算结果表明:与圆形孔相比,心形孔能有效抑制反向旋转涡对的生成,冷却气流的贴壁效果得到明显提高,同时心形孔的扩展出口结构使得冷却气流在展向上的分布更为均匀,展向平均气膜冷却效率得到显著提高;在吹风比为0.5~2.0内,心形孔的全局平均冷却效率相对于圆形孔分别提高了70.93%,246.94%,598.9%和879.07%;从热流比分布来看,心形孔在吹风比为1.5下的热流比值最低,表征在吹风比为1.5下心形孔对壁面的保护效果最好.      

超燃冲压发动机流量匹配机理

采用基于集总参数方程的燃烧室性能计算模型,辅以临界面积法,应用于隔离段和燃烧室的一维流场计算,实现了双模态超燃冲压发动机各种模态下的隔离段和燃烧室的流量匹配计算,并分析了流量匹配工作机理.结果表明:在未分离超燃模态与分离超燃模态下,增加燃烧室供油流量,隔离段和燃烧室流量匹配是通过流道中的喉道处马赫数降低来实现的;在跨燃模态与亚燃模态下,增加燃烧室供油流量,流量匹配主要是通过提高燃烧室流道中热力喉道处的总压来实现的.      

Measurement of mass diffusion coefficients of O_2 in aviation fuel through digital holographic interferometry

The experimental apparatus to measure the mass diffusion coefficients of O2 in aviation fuel was constructed based on the digital holographic interferometry method. The theory of mass diffusion coefficient and interference image processing were introduced in detail. The accuracy of the experiment was verified by measuring the mass diffusion coefficient of 0.33 mol/L KCl in aqueous solution at 298.15 K. The mass diffusion coefficients of O_2 in RP3 and RP5 aviation fuels were measured at temperature from 278.15 K to 333.15 K, and the Arrhenius equation was employed to fit the experimental data. In terms of the Stokes-Einstein equation, the viscosities of these two aviation fuels were tested to estimate the correlation among mass diffusion coefficient, viscosity and temperature. A uniform polynomial calculation correlation was proposed to predict the mass diffusion coefficients of O_2 in both RP3 and RP5 aviation fuels, and its accuracy is considerably higher than that of the Stokes-Einstein equation.     

WZ-5燃烧半水煤气并注水蒸汽的试验

提供了常压下,在 WZ- 5燃烧试验台上,采用中热值半水煤气为燃料的,向燃烧室喷注水蒸汽的各种方案的试验研究。研究了烧中热值燃料的燃烧性能以及在燃烧室中注入水蒸汽增加工质及增大燃气比热,在燃烧室高温区中注入水蒸汽,降低温度抑制 NO_x 的生成。      

中空纤维膜机载制氮装置的数学建模分析

建立了考虑浓差极化现象的微分方程数学模型,并用正交配置法求解,对中空纤维膜机载制氮装置(OBIGGS)进行了分析,并对部分状态点的计算结果进行了实验验证.结果表明:富氮气体(NEA)中氧气的质量分数随着进气温度的升高而降低,在达到最小值后又呈上升趋势;在进气与排气压力差保持不变的情况下,随着中空纤维膜排气压力的下降,中空纤维膜的富氮气体质量流量逐渐增加;随着中空纤维膜丝长度的增加,丝内气体的压降和富氮气体质量流量均有所增加;中空纤维膜空气分离制得的富氮气体质量流量越大,则所需进气的质量流量越大,且富氮气体氧气质量分数越高.      

中间层厚度对陶瓷扩散焊接头应力分布影响的数值分析

采用有限元计算方法，对平面应变、平面应力两种应力状态下含中间层的陶瓷接头界面上应力分布规律进行数值分析。得到了接头界面上正应力和剪应力沿板定的分布曲线。结果表明，中间层材料厚度小时，两种应力状态下的正应力和剪应力在接头界面上滑板定的分布基本均匀。中间层材料厚度增大，平面应变状态下试件中部的正应力增大并出现峰值；平面应力状态下在试件自由瑞边缘正应力增大并出现峰值。两种应力状态下界面上的剪应力由板中央到板边缘逐渐增大，界面上边缘处最大剪应力的值也随中间层厚度的增加而增大。因此，在保证接头接合良好的前提下，应采用小厚度的中间层。     

1.5级涡轮实验台前腔燃气入侵实验

在1.5级涡轮实验台上,针对不同转速和不同封严质量流量,分别采用稳态压力测量、瞬态压力测量和二氧化碳体积分数法对燃气入侵现象进行了实验研究,以确定不同工况下主流的动静叶相互作用对燃气入侵的影响以及入侵到盘腔的燃气在腔内与冷气的掺混过程.结果表明:静叶后的时均压力在周向有明显的周期分布,部分区域封严环外的压力高于腔内压力,且随着封严质量流量增加,这一区域逐渐缩小;动叶扫掠带来同频率的压力波动,对于燃气入侵的发生有着重要的影响;通过二氧化碳体积分数实验获得了最小封严质量流量,并得到在封严质量流量不足情况下燃气沿着静盘侧入侵盘腔的结果.