碱金属引射对尾流化学动力学和电子密度的影响

用把碱金属（Na）加入到空气尾流中进行数值模拟的方法,研究了碱金属引射对高超声速尾流化学动力学和电子密度的影响。采用的化学反应模型为包括O2、N2、O、N、NO、NO＋、Na 、Na、O2-及e-的10组元模型。计算结果表明,碱金属引射,使尾流电子数密度增加1～3个量级。计算结果与文献中的数据符合较好。     

含碱金属杂质的高超声速烧蚀流场特性数值模拟

针对高超声速飞行器的烧蚀现象,采用Newton迭代的空间推进方法,数值求解化学非平衡的抛物化Navier-Stokes(PNS)方程.以黏胶基碳布-氨酚醛树脂为烧蚀材料,建立了含碱金属杂质的18组分33反应的化学反应模型,计算了小钝锥体的壁面温度,并和现有文献结果进行对比验证.在此基础上,对含碱金属杂质的高超声速烧蚀流场进行了模拟,深入分析了含碱金属杂质的烧蚀反应对流场温度、压力和电子数密度的影响.结果表明:烧蚀反应使得流场温度较非烧蚀情况增加了10%~15%;考虑碱金属杂质Na后,流场电子数密度增大了1~2个数量级.      

高超声速钝体尾流光电特性的计算

本文提出了一个高超声速非烧蚀钝体尾流全流场的简化模型和简化的空气化学模型。该模型用于化学非平衡尾流积分电子密度与NO_2化学发光辐射等参数的计算。对于NO_2化学发光辐射既求得了整个频谱范围内的总辐射,也求得了分谱辐射强度。计算结果与弹道靶测量数据进行了比较,它们之间的一致性良好。     

污染对氢-空气燃烧影响的化学动力学分析

本文作了H_2O、O、H、OH和NO对氢-空气燃烧影响的化学动力学计算.计算范围是:压力0.02～0.2MPa,初始温度:910～2000K,当量油气比为1.0.有限速率化学反应模型包括38个反应式和13种组分的H_2-O_2-N_2系统,模型反应了氢-空气化学反应的当前最新研究成果.计算结果表明,在某些温度范围内,H_2O、O、H、OH和NO对着火时间有很大影响.反应时间并不受初始成份中自由原子和基的影响,但是初始成份中含有H_2O,将缩短反应时间.本研究关系到污染对燃烧的影响以及对于超音速燃烧冲压发动机地面试验的影响.燃烧化学动力学显然是污染影响中的一个重要因素.     

高焓风洞中带翼钝锥体尾流电子密度的测量与分析

报道了在爆轰驱动高焓激波风洞中开展带尾翼钝锥体电子密度测试的相关研究工作进展。试验气流为4km/s,密度为0.001kg/m3。诊断尾翼对尾流的影响时,为不影响流场并获得足够的空间分辨率采用针状静电探针;实验结果给出带尾翼模型对尾流电子密度影响的定量结果及受影响的空间区域。     

小钝头锥无粘平衡、非平衡绕流流场的数值分析

本文给出钝头细长锥无粘、平衡和非平衡高超声速绕流流场的数值分析。对存在化学反应的空气,考虑了13种气体组分和22个化学反应,包括N_2、O_2和NO的振动松弛与离介过程之间的耦合反应。头部亚-跨声速区和后身高超声速区分别用直线法和有限差分法进行计算。在来流条件M_∞=20～26,Re_∞=(p_∞u_∞α)/μ_∞=1.5×10~3～5.5×10~6(α为头部半径)的参数范围内,给出了绕流流场和等离子体鞘层诸特性,如激波形状、电子密度和电子碰撞频率等。     

电离非平衡粘性激波层低雷诺数钝体绕流

本文从化学非平衡粘性激波层控制方程出发,用数值法计算了高超音速钝体绕流头部区流场(包括驻点线和音速线)。考虑了电离空气的七种主要成份(N_2,O_2,N,O,NO,NO~ ,e~-)以及相互之间的七个化学反应方程式。认为化学反应的速率是有限的,非平衡扩散采用双组元扩散模型。对激波滑移、壁面滑移、完全催化壁和非催化壁等不同的边界条件进行了计算,得到了流场参数和电子密度分布情况,并与有关文献进行了比较。     

高超声速再入体烧蚀流场计算分析

从化学非平衡NS方程或其简化形式(粘性激波层方程、PNS方程)出发,求解高超声速再入体的绕流和底部尾流流场。首先讨论了碳基材料烧蚀壁面条件的确定方法,计算了包含和不包含碳-碳烧蚀产物的再入体半球的化学反应绕流流场,并和文献结果进行了对比。在此基础上,计算了两个高度条件下,包含和不包含碳-酚醛烧蚀产物的再入小钝锥的绕流和尾流流场,分析了烧蚀产物对流场电子数密度、温度等流动参数的影响。     

烧蚀产物对空气辐射影响的实验研究

本文利用低密度激波管和辐射计,在不同的温度和密度状态下,就杂质辐射在分谱和累积辐射条件下对高温空气平衡辐射的影响进行的实验研究;给出了在不同条件下的实验结果,并做了简单分析。实验结果表明:少量杂质对高温空气平衡辐射的影响较大,这为计算尾流辐射场时,考虑杂质辐射的影响提供了实验依据。     

化学非平衡边界层绕流计算

本文介绍用化学非平衡边界层方程组,数值计算高超声速钝头体绕流流场、组元浓度和电子密度分布情况。计算中考虑了电离空气的七种主要成分以及相互间的化学反应,比较了变熵、普朗特数及Lewis数的影响,讨论了完全催化壁和非催化壁两种情况,计算结果与文献及试验进行了比较,结果令人满意。     

电弧风洞内球锥全目标绕流场计算分析

针对电弧风洞高焓气体来流条件,通过数值求解三维化学非平衡Navier-Stokes方程,分析无烧蚀球锥全目标绕流场真实气体效应.考虑空气电离化学反应系统,包含7个化学组元和6个反应方程.差分算法是隐式NND格式.依据实验来流条件,计算了不同壁温及不同壁面催化特性下的流场结构和电子数密度,观察了来流非均匀性的影响.结果表明.尾迹区对模拟条件较敏感,其电子数密度随壁温降低而升高,完全催化壁的电子数密度约比非催化壁低近一个量级,实验状态的计算应采用风洞实际来流条件,     

再入飞行器尾迹流场及其雷达散射效应研究

对再入飞行器等离子体尾迹及其雷达散射特性进行了分析、研究和大量的计算。讨论了物形、流场各因素对尾迹雷达散射截面的影响。流场计算使用准一维粘性尾迹方程，以修正基尔方法（多值法）求解，用一阶Ｂｏｒｎ近似完成亚密雷达散射截面（ＲＣＳ）计算。计算中使用８组元混合空气、１４个非平衡化学反应模型，考虑５种不同尺度的小钝头锥形物体，沿再入轨道取６５至３４公里，共１３个高程的飞行条件。通过计算得到了再入体尾迹各流场参数、电子密度分布及湍流亚密尾迹的ＲＣＳ。结果说明再入钝锥细长体粘性尾迹的转捩特性对于等离子体的散射性质具有决定性的作用；再入弹头尾迹等离子体对地面单站雷达发射波的回波主要来源于尾迹湍流亚密的非相干散射；对确定的波长，当环境雷诺数达到临界值之后，可能出现ＲＣＳ的突增现象；不同物形及来流条件造成尾迹转捩位置的改变，从而影响ＲＣＳ的数值及其沿轨道的分布；改变尾迹颈部初值会引起ＲＣＳ值的明显变化。     

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

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

钝体化学非平衡粘性流数值模拟

用有限差分法求解Ｎ Ｓ方程,数值模拟了二维钝楔头部和近尾迹非平衡化学反应流场。模拟条件为Ｍａ∞＝２５,飞行高度ｈ＝７０ｋｍ,半楔角θ＝２０°,钝度比为０６０;气体组分为７个组元,共１１个化学反应式。头部流场结果与已有计算结果符合较好。在Ｍａ∞＝２５条件下,底部近尾迹较大区域内具有较高的温度和较大的电子密度,且下降趋势缓慢。     

高超声速电离绕流的数值模拟

本文从三维N－S方程出发，采用稳式L－U分解算法和7组分15个反应的化学模型，数值模拟高超声速电离空气绕流。首先采用对称TVD格式、AUSMPW＋格式和Van Leer的矢通量分裂格式计算了高超声速球头绕流，并对它们的计算结果做了对比分析。然后用前两种格式，对RAM－C飞行试验模型三个再入高度（81km、71km、61km）的流场进行了数值模拟，计算的流场电子数密度值和试验测量数据符合较好。     

锥体高超声速尾流辐射特性

应用显式有限差分法求解了轴对称尾流边界层方程,求出高超声速锥体尾流化学非平衡流场。在求得温度、速度、密度和粒子浓度场的二维分布的同时,详细求解尾流中各种粒子的辐射强度和辐射量的二维分布。     

电离对高超声速稀薄流飞行器气动热影响

将电离反应模型扩展到（direct simulation Monte Carlo,DSMC）方法中,研究了电离反应效应对高超声速稀薄流飞行器气动热的影响特性.针对稀薄流场中电子出现带来的实际困难,引入“捆绑法”思想处理电子在流场中的运动,并给出了电离反应模型及电离反应处理方法.在以RAM-C Ⅱ飞行器外形为例对增加了电离反应的DSMC代码进行验证的基础上,以“星尘号”探测器外形为研究对象,针对不同飞行高度下5组元混合气体模型（无电离）和11组元混合气体模型（含电离）的化学非平衡流动开展了数值模拟,细致分析和对比了电离反应效应对探测器气动热的影响规律.研究结果表明:采用的电离反应处理方案能够模拟带电离反应的高超声速化学非平衡稀薄流动.在飞行高度为60km时电离反应对探测器气动热的影响最为强烈,使探测器的驻点热流密度降低了5.12%,电离反应对探测器气动热的影响随气体稀薄程度增加而减弱.