带配平翼钝体高超声速粘性绕流的数值模拟

采用ＮＮＤ差分格式，通过求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了带配平翼钝体的高超声速粘性绕流。文中首先以球头驻点压强和热流为参数，讨论了在差分方程左端采用ＬＵ－ＳＧＳ、ＬＵ－ＡＤＩ和对角化ＡＤＩ三种不同隐式算法收敛效率的异同。然后以带配平翼钝体的高超声速粘性绕流为模型，对所研制的程序进行了计算验证。在计算中采用了代数方法和求解椭圆型方程方法相结合的网格生成技术，针对配平翼外形给出了贴体性、     

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

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

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

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

跨声速钝体绕流数值模拟研究

通过数值求解ＮＳ方程,进行了三方面工作（１）格式和算法对返回舱跨声速流场特性的影响研究;（２）网格形态对捕捉声速弱激波的有效性研究;（３）模拟了无支杆和有支杆返回舱跨声大攻角绕流流场,给出了后球体上的跨声速激波和后支杆上的第二道激波。     

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

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

绕平板上直圆柱体的高超声速湍流分离流

本文给出直圆柱上游中心线上壁面压力和热流特性及其分离激波运动的特征参数。试验条件是:自由流马赫数为7.8和5,单位长度雷诺数为3.5×10~7/米和4.7×10~7/米。结果表明:1.直圆柱上游中心线上的流场是一个有二次分离的流场,平均壁面压力分布呈双峰型,平均热流分布略有不同,出现一近似平台区,但在分离激波上游影响起始区域和直圆柱前缘邻近区域(x/D<0.5),两者变化极其相似。来流马赫数愈高,分离激波在中心线上的上游影响区长度和峰压及峰热值愈大。2.分离激波是极其不稳定的,出现大尺度低频运动,激波流向运动尺度约0.5D,振荡频率主要集中在2.5kHz以下。     

过渡领域高超声速圆柱绕流直接模拟

本文用直接模拟蒙特卡罗方法模拟了再入速度为７．５ｋｍ／ｓ，高度分别为８０、８５、９０ｋｍ，物壁是完全漫反射有限催化壁的圆柱绕流。在８０ｋｍ考察了壁面催化率的影响。     

特型钝体绕流动态特性的低速实验研究

在低速自由射流风洞上 ,利用热线及相关设备对由半球头部、颈部及旋成体构成的特型钝体周围非定常流动的基本特性、联合特性及湍流度分布进行了测试与分析 ,实验工况为迎角α=- 1 5°～ 1 5°、侧滑角 β =0°～ 1 0°。结果表明 ,在上述工况范围内 ,该钝体周围流场中速度脉动的能量分布平坦 ,属宽频带随机信号 ,流场中没有发生明显的流动分离 ,流场动态品质良好 ;凹陷区内气流的脉动以不同的速度向下游空间传播 ;在上述α、β变化范围内 ,钝体颈部凹陷区的最大湍流度高达 1 0 .9～ 2 4.8%。     

高超声速钝双楔绕流流动转捩与分离流动的壁温影响

为研究壁温对吸气式高超声速飞行器进气道转捩流动的影响,选取钝双楔这一典型外形,基于德国Aachen工业大学Thomas与Herbert所开展的双楔高超声速风洞试验,分析了一些已有的计算流体力学(CFD)研究内容,并结合本文不同方法的CFD数值模拟结果,讨论了不同壁面温度对该双楔模型高超声速绕流流动转捩与流动分离的影响。对于双楔模型,流动分离一般发生在拐角附近,由于流动分离旋涡的剪切作用会诱发流动转捩,转捩又会改变流动分离强度、分离涡尺寸,若分离流动存在非定常特征则将导致非定常旋涡运动与流动转捩的复杂相互作用。通过比较已有文献的CFD数值模拟结果与本文计算结果,表明只有按照转捩思路开展的数值模拟才能够反映该风洞试验情况。计算结果与试验数据的比较显示,文献中按照第一压缩面层流与第二压缩面湍流状态计算得到的结果能够在一定程度上与风洞试验数据相符,本文使用MUSCL格式、剪切应力输运(SST)湍流模型与γ-Reθ关联转捩模型这种计算方法,得到的结果与试验数据符合较好,正确地反映了风洞试验情况。分析还表明,在分离流动之前的区域,随着壁面温度的升高,壁面热流会下降,近壁区域黏性系数变大,边界层内速度剖面不饱满,速度边界层较厚,厚的速度边界层容易发生流动分离现象。     

绕凸台的高超声速分离流动研究

本文研究高超声速流动绕三维低凸台的流动特性。凸台高度与边界层厚度之比介于0.5～0.8,凸台周边倾角介于14°～45°。在高超声速风洞中,来流马赫数为5,单位雷诺数为2.6～6.0×10~7/米。实验过程中测量了模型中心线及特殊部位之表面压强分布,根据纹影记录及表面流动显示确定激波与边界层干扰流场特性及分离区的变化。发现分离区在中心线上距凸台最远,而在凸台肩部距周边沿法线方向最近,确定了中心线上最大压强比的位置。     

高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。     

超音速钝体无粘绕流的计算

研究用三维特征线相容性关系的迎风差分格式,数值模拟超音速和高超音速钝体无粘绕流流场。用钝体反方法计算头激波及其波后的流场,用轴对称流特征线法确定初值面上的流动参数。算例的计算结果与实验数据比较吻合     

飞机横侧动稳定性及反应的奇异摄动理论

应用奇异摄动法分析横侧小扰动运动的动稳定性问题,将四阶问题分解为小根、大根和复根模态,它对应于经典理论中的螺旋、滚转和荷兰滚模态,并给出了动态反应过程的解。     

飞船外形高超音速流场N-S方程数值解

用二阶迎风TVD格式求解N-S访方程,对简化的飞船外形高超音速粘性绕流进行了数值模拟。给出了物面压力分布及流场结构,与实验及其他计算结果符合较好。计算表明,肩部圆弧半径对后体流场影响较小,但严重影响肩部附近热流分布。当肩部圆弧半径是头部半径的20％时,肩部热流低干驻点处热流;半径为5％时,肩部热流超过驻点处热流30％。     

钝缘舵高超音速湍流分离特性

给出由半圆柱前缘舵诱导的高超音速湍流分离的实验结果。实验气流Ｍａｃｈ数为７．８,单位长度Ｒｅ数为３．５×１０７ｍ－１。结果表明：钝缘舵诱导的湍流分离极不稳定,分离激波出现大尺度低频振荡,使壁面压力和热流率无量纲标准偏差在主分离线附近达最大值。Ｍａｃｈ数愈高,最大无量纲标准偏差值越大。在前缘区前缘直径是控制分离流场尺度和平均壁面压力、热流率分布的主要参数     

近高超声速流中旋成体俯仰阻尼导数的计算

本文用改进的二阶激波-膨胀法导出了旋成体法向力静导数的二阶Eggers公式,并计算了俯仰尼阻导数。它与修正的Newton理论相比,由于在压力分布的计算中考虑了流线上前方气动参数的影响,理论上更加完善。计算结果与实验数据比较表明,这是一种比较有效的工程计算方法。     

广义奇异摄动降阶方法在控制器降阶中的应用

 首先推导了广义奇异摄动降阶和奇异摄动降阶的 Schur方法,解决了均衡变换可能出现病态的问题,并且算法具有较好的鲁棒性。然后,将其应用于控制器降阶,并同其它方法比较,认为原系统带宽附近的控制器降阶误差对降阶效果影响较大,而广义奇异摄动降阶法通过参数α的选择,能取得较好的降阶效果。可供工程设计参考。     

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

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

缝隙大小对高超音速分离流的影响

本文提出了一些缝隙大小对高超音速分离流影响的实验结果。结果表明,缝隙大小明显地影响分离区范围,缝隙加宽使分离区缩小;缝隙存在也改变了分离状态下控制翼面的压力分布与传热分布。     

钝锥有攻角分离流动的数值模拟及其分析

本文定性分析了开式分离的性状,并对钝锥有攻角超声速绕流的开式分离作了数值模拟。分析指出,开式分离可能存在两种形态,第一种分离线的起点为正常点,第二种分离线的起始为鞍、结点(包括螺旋点)的组合。对于第一种形态,分离线的起点是横向分离的起始点,除分离线外,分离面上的流线不是从分离线的起点发出的。对文中计算的情况,流动属第一种开式分离。计算证实了定性分析的结论。计算和分析均指出,对第一种开式分离,在分离的起始区域,分离流面尚未卷曲,但在下游,则变成卷曲面。文中还研究了围绕物体的流管在分离诱导下的变形情况。