低温烧蚀材料的烧蚀实验研究

本文给出了低温烧蚀材料的高超声速风洞烧蚀实验研究的结果。模型外形是半锥角为9°的球锥体,头部半径为15毫米和25毫米,材料为蜂蜡和樟脑,实验条件是M_∞=6,α=0°、2°、6°和8°,Re_∞为(1.90～4.42)×10~7米~(-1)。实验结果表明:随着Re_∞的增加,可以获得层流、转捩和湍流烧蚀表面形态和端头烧蚀外形;证实了在实验过程的后期出现平衡外形,在平衡外形中存在角点;融熔型材料和升华型材料的烧蚀表面形态差别较大;有攻角和零攻角的烧蚀外形有较大差异。     

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

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

高超声速半球绕流流场电磁散射特性分析

分析再入流场对飞行器再入通信及空间散射特性的影响,研究电磁波与再入等离子体作用机理.采用分段线性电流密度递归卷积时域有限差分方法计算导电金属半球高超声速绕流流场电磁散射特性,分析半球等离子体包覆绕流流场雷达散射截面(RCS)随入射电磁波频率、双站角、飞行马赫数和高度变化特性.计算表明,前向是全方位散射中RCS取得最大值的方向.马赫数Ma≤10时,马赫数及高度变化对半球高超声速绕流流场L,S波段后向RCS和双站RCS影响很小;在L,S波段,绕流流场及半球本体的后向和前向RCS差距较小.马赫数Ma=14,16时,过密等离子体鞘套的形成增大了本体投影面积,在L,S波段,绕流流场前向RCS比本体前向RCS大;绕流流场存在对半球本体双站散射特性影响很大.      

高超声速平头圆柱绕流化学反应流场的DSMC/EPSM混合算法研究

本文发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的EPSM耦合模型，并通过混合参数进行流区的自动识别，将EPSM方法与蒙特卡罗直接模拟方法(DSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的轴对称情况下高超声速平头圆柱绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

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

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

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

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

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

跨声速半球柱绕流流场的激波系结构

本文在1相似文献   

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

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

高温化学非平衡湍流边界层直接数值模拟

高超声速飞行器在特定工况下飞行时，其表面会存在湍流与化学非平衡耦合作用，但考虑化学反应的高温湍流边界层的研究工作还十分有限。选取高超声速楔形体头部斜激波后的流动状态，分别采用热完全气体模型、空气化学反应模型开展直接数值模拟(DNS)研究，对比分析了化学非平衡效应对湍流统计量、湍流脉动量的影响趋势，并研究了高温湍流边界层标度律。计算结果表明，以吸热离解反应为主的化学非平衡效应会使边界层平均温度和脉动温度显著降低，但对平均速度剖面、速度脉动、雷诺应力的影响较小；在高温化学非平衡湍流边界层的对数区，各流动参数的拓展自相似理论的相对标度指数基本符合标度律。     

完全催化壁驻点高超声速流动加热地面模拟方法研究

在地面高焓风洞中准确模拟高超声速头部驻点区化学反应流动加热,从平衡边界层和冻结边界层驻点气动加热公式出发,分析了离解焓对不同催化特性驻点加热的影响、地面风洞非平衡来流条件下头激波后的流场参数变化规律和与天上飞行条件的差异,分析了壁面催化特性对流场参数的影响,建立了地面高焓风洞模拟完全催化壁驻点气动加热的模拟准则:只要在风洞条件下实现对驻点压力ps、驻点速度梯度βe、来流总焓hs的模拟,则能够复现天上状态的气动热载荷。针对不同半径的飞行器头部驻点开展化学非平衡流场数值模拟进行了验证,表明:地面风洞无法完全模拟飞行器头部绕流流场参数,材料的催化特性对壁面附近流场中的N、O组元和壁面热流影响较大;地面风洞采用半径1∶1的模型模拟完全催化壁驻点加热会偏低于天上状态,而根据三参数模拟准则确定的模型尺寸能够同时复现天上状态驻点区的焓、压力、热流,且随着飞行器头部半径增大头部驻点线近壁面附近的温度、组元N、O梯度与天上会趋于一致。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

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

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

不带稳定翼飞船返回舱俯仰动稳定性研究

联盟号飞船返回舱,如果不存在稳定翼,从高超声速到低Mach数时,力矩曲线有由一个平衡迎角(α1),经二个(a1、α'2)到三个(α1、α2、α3)平衡迎角的情况出现.当有一个平衡迎角时,可能有一临界Mach数Mcr,当M∞＞Mcr时,平衡迎角处是动稳定的;M∞＜Mcr,变为动不稳定的,同时出现极限环.也可能不存在这个Mcr,这时随Mach数降低到M∞=M'cr,即2个平衡迎角的情况,在力矩曲线与α轴相切的平衡点α'2,会出现鞍、结点分叉.M∞＜M'cr后,对于存在三个平衡点的情况,在α和α的相图上,它是结点-鞍点-结点的结构.当α1处是稳定结点形态时,只有在很小的扰动情况下,它才能返回α1的状态.如果α1为不稳定结点状态,在α1附近,微小的扰动,使迎角由α1迅速变成α3或者更大.这是混沌的先兆.     

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

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

气体自发光在高超声速流场显示中的应用

提出了利用高焓气体自发光作为高超声速流场显示的方法，介绍了在使用高焓运行的激波风洞中，对二维模型的高超声速绕流流场使用此种方法的初步结果，可观察到二维棱形柱的尾流和马赫波的相交。结果表明此种方法不需外加光源，对于结构限制无法设置透明部件的模型，无疑是简单可行的。     

钝头飞行器大攻角绕流流场计算

本文用有限差分法求解Ｎ－Ｓ方程，对具有大面积底部飞行器的高超声速大攻角绕流流场进行了计算，通过积分物面压力给出了飞行器的气动力特性，并研究了攻角的影响，物面压力系数和气动力系数与实验数据进行了比较，符合较好，结果可应用到返回舱和弹头的再入气动特性研究。     

垂直于物面的横截面上流态的拓扑

对于旋成体有迎角绕流,本文研究了垂直于物面的横截面上的流态结果。利用解析方法,给出了物面轮廓线上、迎风和背风对称线上以及全横截面上截面流线方程奇点的分布规律以及奇点沿流动方向的演变规律,从而可定性地确定截面流态的结构。通过求解ＮＳ方程数值模拟了钝锥有迎角高超声速绕流流场,计算获得的垂直于物面的横截面流线的拓扑结果和理论分析完全一致。     

气体模型对高超声速再入钝体气动参数计算影响的研究

本文采用完全气体模型、振动激发气体模型、平衡气体模型、一温度非平衡气体模型、两温度非平衡气体模型和三温度非平衡气体模型进行了钝体高超声速绕流流场的数值计算。分析了各种气体模型对激波脱体距离、壁面热流、温度分布和密度分布等的影响。结果表明：采用两温度（或三温度）非平衡气体模型计算的激波脱体距离更接近实验情况；采用平衡气体模型计算的壁面热流量高，一温度非平衡气体模型完全非催化的壁面热流最低，而完全气体模型的壁面热流在两者之间。最后给出了各种气体模型的密度等值线比较图。     

高超声速稀薄气流中球阻力系数的一种简单公式

本文提出了以粘性参数K=M_∞/√Re_D为变量的球阻系数的简单关系式。它适用于高超声速稀薄气流的情况,在Re_D=5×10~2～10~5;M_∞=8～25的范围中与实验数据吻合。此外,按此式计算的球阻系数而得到的标定球的动压值与由皮托压力测定所得结果一致。