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

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

热化学非平衡来流条件下热化学模型影响研究

为研究热化学非平衡来流条件下热化学模型等计算设定对斜激波压缩流动计算结果的影响,针对尖劈构型和相应的前缘钝化构型的高焓激波风洞实验,采用多种计算设定开展详细的数值模拟研究。计算结果表明,计算采用不同热化学模型,以及来流设定为振动冻结/平衡/非平衡状态,会导致斜激波激波角等参数存在一定差别,其中激波角差别可达约2%。当来流速度一定时,过斜激波后分子内能增量在平动转动能和振动能上的分配方式的差别决定了激波角的差别。前缘钝化情形下,采用不同计算设定所得激波角之间的关系和尖前缘构型的规律一致;但是,采用不同计算设定所得斜激波到壁面距离之间的关系和尖前缘构型的规律有差别,这源于钝化前缘的激波脱体距离的影响。对于自由来流下的斜激波压缩流动问题,若考虑了分子振动能激发但未考虑热力学非平衡(例如热完全气体模型、考虑空气反应的单温度模型等),就斜激波激波角等参数而言,计算误差比量热完全气体模型计算误差更大。     

磁场对高超声速弱电离气体流动的影响

对偶极子磁场作用下的三维钝头体高超声速黏性绕流的化学非平衡流动进行了数值模拟。电导率利用组分公式计算,化学模型为7组元、6反应模型。应用诱导磁场方法将磁场与流场耦合,控制方程的空间离散采用有限体积法,扩散项用中心差分格式,对流项采用二阶迎风格式,时间推进采用3步Runge-Kutta法。计算结果表明,外加偶极子磁场使激波脱体距离增加,壁面摩擦系数和表面热流密度增大。与无磁场作用时相比,在0.153T外加磁场作用下,冻结流中的激波脱体距离增加约3倍,局部壁面摩擦系数最大增加63%,局部表面热流密度最大增加61%;而化学非平衡流中的激波脱体距离增加约0.5倍,局部壁面摩擦系数最大增加47%,局部表面热流密度最大增加31%,且化学非平衡流中激波层内温度的最大值约为冻结流中的64%。     

模型表面温度对光电特性影响的实验研究

对在JF-10氢氧爆轰驱动高焓激波风洞中开展烧蚀材料模型表面温度对光电特性影响的实验方法进行了探索。风洞实验状态的驻室总压约为20MPa,驻室总温约为6000K,自由流速度约为4km／s。实验以锑化铟多单元红外成像系统与电离列阵测试装置为测量手段,用以烧蚀材料为头部的球头钝锥体模型,实验测量激波层中红外辐射的横向剖面分布和壁面附近电离密度的剖面分布,并在烧蚀材料加热和不加热两种情况下进行了对比实验。实验结果表明：烧蚀材料的加热加强了模型头部激波层中的红外辐射,同时也增大了模型表面的电子密度。     

纳秒脉冲等离子体激励器用于圆柱高速流动控制的数值模拟

利用数值模拟,研究了纳秒脉冲介质阻挡放电（NS-DBD）等离子体激励器在圆柱高速流动控制中的应用。首先,研究了单电极NS-DBD等离子体激励器在静止空气中放电后的流场特性。研究表明在介质阻挡放电形成的等离子体区域,有局部能量快速注入,放电结束5 μs后在上极板后端点位置形成了一个局部温度高达900 K的热点,由此引发很强的压力扰动,形成以上极板后端点位置为中心,扩散速度约为声速的半圆形压缩波。在此基础上,通过数值模拟研究了NS-DBD等离子体激励器布置在直径为6 mm的圆柱上,来流马赫数为Ma_∞=4.6时,对圆柱脱体激波的控制作用。研究表明介质阻挡放电形成的半圆形压缩波对于脱体激波有很强的干扰作用,激波距离增加了15.7%,激波强度也有相应的减弱,导致阻力减少了13%。     

均匀磁场中高超声速弱电离气体流动数值模拟

采用7组元化学模型并应用组分公式计算电导率,通过求解黏性MHD(magneto-hydro-dy-namics)方程组,研究了不同强度均匀磁场对三维钝头体高超声速绕流化学非平衡流动的影响.结果表明,随着磁场的增强,激波脱体距离逐渐增加;总阻力系数和壁面温度逐渐减小.在B_y=0.03 T磁场作用下,与无磁场的结果相比,化学非平衡流中的激波脱体距离增加约7%,总阻力系数减小约5%,局部壁面温度最大降低74%;而冻结流中的激波脱体距离增加约43%,总阻力系数减小约6.9%,局部壁面温度最大降低18%.      

运动激波绕尖劈流动的研究

在激波管中进行双尖劈二维外形的激波绕流数值计算和实验研究,入射激波,Ｍ＝１．８,用激光全息双爆光技术定量测定各瞬时的密度场。计算采用欧拉方程和有限体积法进行离散并采用高精度的ＴＶＤ差分格式,计算与实验两者之间的比较,表明本文采用的计算方法对于解决尖点的绕流十分奏效,根据计算所获得的运动激波系,能够判断局部区域实验测量的密度场梯度方向,从而使密度场的定量测定获可靠度的结果。     

准定常强激波马赫反射波形的数值研究

应用DCD频散控制激波捕捉格式,求解二维、多组分、带有化学反应的Euler方程组,数值模拟了准定常强激波的马赫反射问题.研究结果表明:与经典马赫反射理论相比,在强激波条件下,激波诱导的气体分子振动激发和化学反应使马赫反射的三波点轨迹角变小、马赫杆高度变低、楔顶附体激波倾角变小;马赫杆的相对突出量随入射激波马赫数和楔角的增大而增大,而气体分子的振动、离解等真实气体效应能进一步加剧马赫杆的向前突出.     

用热流探针测量激波速度

激波速度测量是激波管和激波风洞运行状态的一个重要参数,压电传感器或光学方法测速系统成本高,而传统电离探针在激波马赫数较低、波后温度达不到空气电离程度的情况下无法满足实验要求。提出了一种使用同轴热电偶作为测速探针来测激波速度的方法,弥补了电离探针在激波马赫数较低时的不足。介绍了同轴热电偶探针测速原理,并设计了测量激波速度的系统电路。通过信号放大电路锁定激波冲激信号,触发脉冲信号发生电路,实现了一种单通道、多测点的激波风洞测速系统。分别开展以温度与热流为触发信号的风洞实验,结果表明只有使用热流信号才能满足激波测速的时间要求。     

强激波背后发光形态观测及其结果的图像处理

本研究是用实验观测的方法判明低密度空气中发生的强激波的构造并进而分析其特性。为此,首先利用自由活塞二段隔膜激波管获得强激波,并以新型高速照相机拍摄激波背后的发光域,然后以SPICCA图像处理机对拍摄的照片进行等辉度线和浓度变换等处理。研究结果表明:(1)激波背后的发光受其速度影响很大,当激波速度U_s超过11.0km/s时,发光强度将出现两次峰值;(2)空气中氧元素对激波背后的发光情况有重要影响;(3)初压不同时,发光形态不同,当初压增高时,发光时间和强度都增加。这些结论为从理论上分析与研究强激波背后发光的形态与特性提供了可贵的参考资料。     

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

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

十一组元高温平衡、非平衡空气光谱吸收系数计算

飞行器以很高的马赫数再入大气层时 ,头部激波层、尾迹的气体辐射产生的紫外、可见、红外特征信号 ,是地面监测和反导制导系统探测、识别的主要依据 ;同时 ,头部气体层热辐射也是再入体壁面加热热流的重要源项。激波层内高温气体吸收系数是目标紫外、可见光、红外辐射特性计算、热防护计算和气动流场辐射场耦合计算的基本参数。本文分析了再入过程中高超声速稀薄气体流的非平衡现象 ,采用三温度模型 (电子能温度 Te、振动能温度 TV、转动和平动能温度 TR)来表征激波层内非平衡态的气体组份各个能级占有数分布 ,并由原子分子辐射理论直接计算 N2 ,N+ 2 ,N,N+ ,O2 ,O+ 2 ,O,O+ ,NO,NO+和 e-等 1 1种主要空气组元各种能级跃迁对辐射的贡献 ,最终得到平稀和非平衡高温空气吸收系数计算模型。计算结果可见光区与试验结果符合较好 ,紫外、红外区结果稍差     

气相爆震波衍射现象的数值研究

对不同条件下气相爆震波在突扩管道和带扩散段管道中的衍射过程进行了数值研究.结果表明:起爆管的直径和临界爆震管直径对爆震波能否成功衍射传入主爆管起到决定性的作用.壁面反射激波与爆震波相互作用在邻近壁面处产生高温、高压区域,有助于减少爆震的诱导时间以及横波间距.扩散段扩散角越小时,反射激波强度越大.高强度的反射激波会将爆震波分为3段,在分界点处产生高温、高压点,对爆震波的衍射传播起到支持作用.      

RP-3航空煤油着火特性的实验

在化学激波管中利用反射激波着火,采用壁面压力与OH自发光作为着火指示信号,测量了着火温度范围为1100~1600K,压力为0.1,0.2,0.3MPa,当量比为0.5,1.0,1.5时RP-3航空煤油/氧气/氩气混合气的着火延迟时间,分析了着火温度、压力以及当量比对混合气着火延迟时间的影响,并拟合得到了不同压力与当量比下混合气着火延迟时间的Arrhenius关系式.结果表明:在不同压力与当量比下,混合气的着火延迟时间的对数与着火温度的倒数呈线性关系,同时,随着着火温度与压力的升高以及混合气当量比的降低,着火延迟时间逐渐缩短.      

爆轰波穿越惰性气团时的透射激波参数分析

爆轰波穿越惰性气团再进入可燃混合气体时,透射激波的参数主要受到三个基本过程的影响:头激波两次穿越惰性气团与可燃气体之间接触面导致的激波强度下降、放热反应中止、爆轰波后稀疏波对激波的衰减作用。放热反应中止是透射激波衰减的主要原因,不同种类的惰性气体对透射激波参数影响不大。     

超声速流中三波结构的参数优化研究

对三波结构的参数进行了探讨，给出了在反射激波和主激波后一些气动力变量的比值存在极值，文中称这种结构为优化结构，当给定气流的马赫数时，理论上就可确定优化三波结构的入射激波强度和存在区域，这可为一些超声速气动装置的设计提供参考依据。     

Influence of shock attenuation on tailored operation in free piston shock tunnels

The free piston shock tunnel is a type of shock tunnel with high performance. For this type of tunnel, the influence mechanism of shock wave attenuation on tailored operation is explored by numerical simulation and theoretical analysis. By introducing the normalized velocity, the simple constraint equation for shock wave under the tailored operation is deduced. Moreover, the real gas effect is also taken into account in this equation. Based on the equation, the tailored operation of shock tunnels can be predicted with very few calculations. The present study shows that the change rate of the thermodynamic state of the gas behind the shock wave is inconsistent with the attenuation rate of the shock wave, which is the fundamental reason why the wind tunnel achieves tailored operation at a lower Mach number of shock waves. This lower Mach number of shock waves differs from the corresponding ideal value by a factor, which is about the square root of shock attenuation rate.     

平板上直立圆柱引起的三维湍流分离流的传热特性

 1.实验目的与方法 航天飞行器再入大气层时,其表面突出物引起的三维湍流分离,使局部热流急剧增加。突出物激波与湍流边界层的干扰流场很复杂,需要较多的测量数据揭示其流动规律。为此,在激波风洞中Ma_1=5.2,Re/L=2.3×10~7m~(-1)的气流条件下,测量了平板上直立圆柱干扰区中柱前缘及其邻近平板上的热流分布,并拍摄了激光纹影