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

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

高温高速稀薄流的DSMC算法与流场传热分析

在非结构网格下,研究了热力学碰撞传能的6种方式和化学反应的8种类型,详细给出了上述碰撞传能和化学反应碰撞类型的子程序计算框图.在可变硬球(VHS)分子模型、Borgnakke-Larsen唯象模型、Bird的化学反应几率模型以及壁面CLL(Cercignani-Lampis-Lord)反射模型的基础上,用Fortran语言编制了能够模拟内能松弛、热力学非平衡和化学非平衡的稀薄气体DSMC(direct si mulation Monte-Carlo)源程序,在地球大气层和火星大气层中完成了Ballute减速装置的8个工况计算(其中Knudsen数从0.05变到30.0,飞行Mach数从26.3变到11.2),并与NASA Langley研究中心2007年发表的计算结果作了比较,本文的结果令人满意.另外,这里用|T-T_v|/T去分析热力学非平衡,用Damk hler数去分析化学反应非平衡,用Stan-ton数去分析飞行器壁面的传热效果,所有这些分析对高温部件的热防护设计十分有益.      

辐射加热对返回舱气动热环境影响的数值研究

基于高超声速再入飞行器气动热环境预测分析的需要,建立了高温非平衡气体辐射加热对飞行器热环境影响的计算分析手段。采用数值求解化学非平衡N-S方程的方法,对返回舱绕流流场进行模拟,获得高温空气组分质量分数和温度等流场参数分布。基于辐射传输方程,考虑高温气体组分的主要辐射机制,计算分析高温流场气体辐射加热对返回舱热环境的影响。分析表明,在同一飞行弹道条件下,返回舱大底半径尺寸对气动加热的影响较大,在再入热环境严酷区,辐射加热对物面总热流的贡献达30!;产生辐射加热效应的主要机制是高温流场中O和N原子产生连续谱和线状谱以及N2的第一正带系;物面催化效应对辐射加热影响不大。     

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

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

高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。     

高超声速钝头双锥体高温流场辐射研究

研究再入高温绕流及尾流流场的辐射特性对高超声速飞行器突防和识别技术的发展具有重要意义。本文基于Navier-Stokes方程和辐射传输方程,建立了沿任意视角方向的再入体全目标区流场(包括绕流及尾迹流场)辐射特性计算方法,计算分析了在0.4~12.0μm光谱范围内钝头双锥体高温化学反应流场辐射在不同视角(0°~180°)方向以及不同再入高度(40~60km)下的变化特征。     

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

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

高速高温流场电子能非平衡的数值模拟

采用描述电子能非平衡的三温度模型,结合11组分空气的化学反应模型,对多种高速高温热化学非平衡流场开展数值模拟,并与描述电子能平衡的两温度模型结果进行对比,研究电子能非平衡对高超声速流场特性的影响。圆球弹道靶试验算例表明电子能非平衡不影响激波脱体距离。RAM-C II飞行器的4个飞行工况算例表明,尽管两温度和三温度模型结果存在差异,但二者电子数密度分布的趋势和量级接近,均可与飞行试验数据保持一致,其中三温度模型的预测效果更好。FIRE II飞行器极高温流场模拟结果显示,电子能非平衡几乎不影响飞行器表面的对流传热。     

高超声速飞行器高温流场数值模拟面临的问题

随着高超声速飞行器目标光辐射和电磁散射特性研究的发展和深入,高温流场特性日益引起人们的关注。由于高温流场特性研究中涉及到非常多的复杂气动现象,如气动加热、烧蚀、辐射、燃烧、化学反应以及湍流等,因此其数值模拟面临着诸多挑战。这里基于连续流计算流体力学(CFD)技术和稀薄气体蒙特卡罗直接仿真(DSMC)方法,从化学物理模型建模、方法稳定性与数值求解效率出发,分析了高超声速飞行器外部绕流、尾迹和发动机喷焰三方面的流场特性数值模拟在不同弹道、热防护手段和飞行流域环境下所面临的问题。在此基础上提出了数值求解技术和化学物理模型建模今后需要发展的方向,为有效提高高超声速高温流场特性数值模拟效率、增加流场特性预测精度提供了指导,从而为研究流场对高超声速飞行器目标光辐射和电磁散射特性影响提供有效的基础数据。     

高超音速飞行器前缘驻点气动加热的计算

介绍了以高超音速(马赫数不小于6．0)在地球大气层中飞行的飞行器外表面的附面层与飞行器壁面间热交换的物理特性及高超音速附面层的数学模型，给出了该飞行条件下轴对称形式的飞行器前缘驻点空气动力加热的计算方法。     

Numerical simulation of compression corner flows at Mach number 9

A hypersonic vehicle encounters a wide range of conditions during its complete flight regime. These flight conditions may vary from low to high Mach numbers with varying angles of attack. The near-wall viscous dissipation associated with flows at combined high Mach and Reynolds numbers leads to significant wall heat transfer rates and shear stresses. The shock wave/boundary-layer interaction results in a flow separation region, which commonly augments total pressure losses in the flow and lowers the efficiency of aerodynamic control surfaces such as fins installed on a vehicle. The standard turbulence models, when used to resolve such flows, result in incorrect separation bubble size for large separated flows. Therefore, it results in an inaccurate aerodynamic load, such as the wall pressures, skin friction distribution, and heat transfer rate. In previous studies, the application of the shock-unsteadiness correction to the standard two-equation k-ω turbulence model improved the separation bubble size leading to an accurate pressure prediction and shock definition with the assumption of constant Prandtl number. In the present work, the new shock-unsteadiness modification to the k-ω turbulence model is applied to the hypersonic compression corner flows. This new model with variable Prandtl number is based on the model parameter, which depends upon the local density ratio. The computed wall pressures, heat flux and flow field are compared to the experimental data. A parametric study is carried out by varying compression deflection angles, free stream Reynolds number and wall temperatures to compute the flow field and wall data accurately, particularly in the shock boundary layer interaction region. The new shock-unsteadiness modified k-ω model with variable Prandtl number shows an accurate prediction of initial pressure rise location, pressure distribution in the plateau region and heat flux in comparison to the standard k-ω model.     

耦合传热下激波对超声速气膜冷却影响

针对离散孔式超声速平板气膜冷却，在主流区引入楔角形成激波环境，以研究激波与超声速气膜之间的相互作用。通过计算楔角在0°、15°、20°和25°产生的四种激波强度下，超声速气膜与高温壁面的耦合传热。所得结果表明：适当强度的激波能够抑制气膜入射后产生的反向涡旋对，降低主流对气膜的卷吸，增大壁面平均H2摩尔分数并降低壁面温度。对金属层温度场的分析表明，壁面冷却效果随着激波角的增加而先增加后降低，其中楔角为20°时的流场结构最有利于壁面温度保护。小楔角生成的激波在低冷流马赫数下对冷却效果的改善更明显，大楔角则在高冷流马赫数下更明显，热障涂层（TBC）不影响这种变化趋势；激波的存在削弱了TBC的影响范围。可以揭示超声速气膜在耦合传热条件下的传热机理，为超声速气膜冷却的设计提供参考，或为现有超声速气膜冷却结构的优化提供依据。     

超燃冲压发动机唇口气动热计算研究与分析

采用数值模拟方法,对超燃冲压发动机进气道唇口的气动热进行了计算和分析.依据唇口与前体激波所处相对位置的不同,分三种情形进行了计算研究.结果表明,三种情形下,进气道前缘驻点热流均大于头部.通过算例验证,证明该方法和结论的正确性.研究表明,激波干扰以及前体激波压缩后超声速来流使唇口热流增大,热防护和结构设计时应当予以充分考虑.      

A thorough experimental investigation of shock/shock interferences in high Mach number flows

The Types III and IV interference flows, as defined by Edney, and corresponding heat transfer distributions were investigated experimentally. The model consists of a cylindrically blunted plate and a wedge serving as an oblique shock generator. The ‘thin wall’ technique was used for heat transfer measurements on the cylinder surface. These experiments were carried out in the TsAGI short duration wind tunnel UT-1 at Mach numbers 6 and 16 in air and at Mach number 6.6 in carbon dioxide. The Reynolds number based on the plate bluntness diameter was varied in the range 2.2×10⁴ to 1.6×10⁶. Tests of the cylinder alone (without the wedge) at Mach number 6 and for different Reynolds numbers revealed an influence of incoming disturbances on the stagnation line heat transfer. The influence of the impinging shock location on the interference heat transfer was carefully investigated. Systematic calculations of inviscid flow at Mach number 6 were also performed. Estimations of the maximum interference heat transfer rate, based on these calculations and a boundary layer approach, compare well with the data. Influence of the specific heat ratio on the interference flow was studied. These experiments and calculations revealed important features of interference flow patterns and heat transfer distributions.     

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

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

存在烧蚀的高超声速压缩拐角流动数值研究

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用16组元(N2,O 2,NO,N,O,NO +,N+2,O +,N+,CO,CO 2,C,C2,C3,CN,e-)29个反应的非平衡模型,热力非平衡的双温度模型下,不同反应采用不同控制温度。石墨材料表面反应包括碳的氧化反应、碳催化的O 原子复合反应和碳的升华反应。对15°、18°、24°压缩拐角模型,在自由流 Ma =10~30,总焓值6~55 MJ/kg 范围,分别进行无烧蚀的壁面催化与非催化条件和石墨烧蚀条件下的流场计算,分析各类条件下的流场结构、流动分离特性以及流场热化学参数分布特点,研究壁面条件对流动特性的影响。结果表明：流动分离可能性和分离区范围随着压缩拐角斜面倾角增大而增大,随来流马赫数增大而减小;相对于低壁温条件,无烧蚀的辐射平衡壁温和壁面烧蚀条件下流动分离区增大,斜面上压力、摩阻和热流峰值点也有所后移。     

高超声速进气道再起动特性分析

1引言未来高超声速飞行器飞行必然要经历低马赫数飞行过程,因此高超声速进气道同样要面临低马赫数下的不起动问题,那么进气道一旦进入不起动,如何才能再起动?再起动的特征是怎样的呢?常规内压式进气道再起动过程中存在迟滞回路现象,高超声速进气道的再起动过程是否也有相似的现     

高马赫数下横/逆向喷流干扰流场数值研究

横向喷流和逆向喷流广泛用于高超声速飞行器气动力与气动热控制。采用格心型非结构有限体积法求解基于三温度热化学非平衡模型的全Navier-Stokes方程,对高空、高马赫数来流条件下二维圆柱状构型飞行器的喷流干扰流场进行数值模拟,研究了仅存在横向或逆向喷流以及横/逆向喷流同时存在时的复杂流场结构以及喷流降低热流、减阻、改善升力的具体效果。通过控制变量的方法,探究了不同参数（马赫数、静压）的喷流对流场结构及飞行器的气动力、气动热的影响规律。结果表明：在一定条件下,当逆向喷流与横向喷流同时存在时,下游的横向喷流可以影响到上游的逆向喷流流场结构;逆向喷流可以显著减小高超飞行器阻力,并降低头部壁面热流峰值,而横向喷流对高超飞行器的升力特性有一定提高;在横向喷流已用于飞行器姿态控制的情况下,一定条件下可以同时使用逆向喷流,既可以减阻、又可以降低热流峰值,还可以提升升力。