真实气体效应下高马赫数内转进气道特性研究

为初步研究高马赫数内转进气道在真实气体效应下的工作特性,首先设计额定工作状态Ma=12的高超声速内转进气道,再结合不同气体模型对其进行数值模拟。研究结果表明:化学非平衡气体在流场结构、工作性能和气动加热方面与热完全气体较为相近,与热化学非平衡气体存在一定差别。离解反应发生在边界层内和低速涡流区内,热化学非平衡气体的离解反应程度比化学非平衡气体大。在隔离段内激波反射处,相比完全气体,化学反应气体的静温降低了2000～2500K。高热流区在上壁面喉道位置与下壁面激波反射点位置附近,温度较高的等温壁面、热化学非平衡气体均可降低壁面热流密度,不同壁面条件对隔离段出口性能参数影响较为明显。真实气体效应、壁面温度对隔离段涡流区的影响较为复杂,有待进一步研究。     

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

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

电导率模拟对高超声速MHD控制影响

高温气体电导率是高超声速电磁流动控制数值模拟最重要的参数之一。针对电导率模拟准确性及其对高超声速磁流体控制影响的问题，考虑高超声速飞行器流场中高温气体热化学非平衡效应，采用三维低磁雷诺数磁流体动力学（MHD）数值模拟方法及程序，结合国内外常见电导率处理方法开展典型状态高超声速MHD控制数值模拟，分析电导率模拟对高超声速磁流体流场分布、气动力/热特性的影响。研究表明：磁控热流减缓效果与电导率呈非线性关系，电导率较大时将出现电导率的磁控热饱和现象，其产生的原因可能与化学反应趋向于平衡态存在一定关系；采用定电导率方法，会人为放大磁场洛仑兹力的磁阻力效果，使阻力系数的预测值偏大；不同电导率模型计算得到的电导率分布差异很大，甚至存在数量级的差别，显著影响了磁流体的控制效果，这与电导率模型的适用范围、参数选取原则存在很大关联；对于含多种离解、电离组分的高温气体流动来说，采用基于多电离组分迁移碰撞的电导率模型（本文模型M8），计算与试验一致性最好。     

LENS风洞试验返回器模型气动热特性模拟

国家数值风洞高超声速流动模拟软件HyFLOW的研制对打破国外同类软件的技术壁垒具有重要意义。与国外DPLR软件进行了对比研究，同时系统介绍了HyFLOW软件求解器的数值方法、物理化学模型以及壁面催化特性计算模型等主要方法，采用典型算例对有限催化模型进行了数值验证，最后基于LENS风洞试验146 mm返回器标模外形开展了高超声速气动热特性数值模拟。研究结果表明，HyFLOW软件在高超声速热化学非平衡流动模拟与评估方面的气动力计算精度高，与国外同类软件DPLR相当，同时其壁面催化条件下的气动热计算精度可靠，可信度高。     

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

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

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

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

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.     

Experimental investigation on aero-heating of rudder shaft within laminar/turbulent hypersonic boundary layers

The aero-heating of the rudder shaft region of a hypersonic vehicle is very harsh, as the peak heat flux in this region can be even higher than that at the stagnation point. Therefore, studying the aero-heating of the rudder shaft is of great significance for designing the thermal protection system of the hypersonic vehicle. In the wind tunnel test of the aero-heating effect, we find that with the increase of the angle of attack of the lifting body model, the increasement of the heat flux of the rudder shaft is larger under laminar flow conditions than that under turbulent flow conditions. To understand this, we design a wind tunnel experiment to study the effect of laminar/turbulent hypersonic boundary layers on the heat flux of the rudder shaft under the same wind tunnel freestream conditions. The experiment is carried out in the ?2 m shock tunnel(FD-14 A) affiliated to the China Aerodynamics Research and Development Center(CARDC). The laminar boundary layer on the model is triggered to a turbulent one by using vortex generators, which are 2 mm-high diamonds. The aero-heating of the rudder shaft(with the rudder) and the protuberance(without the rudder) are studied in both hypersonic laminar and turbulent boundary layers under the same freestream condition. The nominal Mach numbers are 10 and 12, and the unit Reynolds numbers are2.4 × 10~6 m~(-1) and 2.1 × 10~6 m-1. The angle of attack of the model is 20°, and the deflection angle of the rudder and the protuberance is 10°. The heat flux on the model surface is measured by thin film heat flux sensors, and the heat flux distribution along the center line of the lifting body model suggests that forced transition is achieved in the upstream of the rudder. The test results of the rudder shaft and the protuberance show that the heat flux of the rudder shaft is lower in the turbulent flow than that in the laminar flow, but the heat flux of the protuberance is the other way around,i.e., lower in the laminar flow than in the turbulent flow. The wind tunnel test results is also validated by numerical simulations. Our analysis suggests that this phenomenon is due to the difference of boundary layer velocities caused by different thickness of boundary layer between laminar and turbulent flows, as well as the restricted flow within the rudder gap. When the turbulent boundary layer is more than three times thicker than that of the laminar boundary layer, the heat flux of the rudder shaft under the laminar flow condition is higher than that under the turbulent flow condition. Discovery of this phenomenon has great importance for guiding the design of the thermal protection system for the rudder shaft of hypersonic vehicles.     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

Energisation of O+ and O+ 2 Ions at Mars: An Analysis of a 3-D Quasi-Neutral Hybrid Model Simulation

We have studied the loss of O⁺ and O⁺ ₂ ions at Mars with a numerical model. In our quasi-neutral hybrid model ions (H⁺, He⁺⁺, O⁺, O⁺ ₂) are treated as particles while electrons form a massless charge-neutralising fluid. The employed model version does not include the Martian magnetic field resulting from the crustal magnetic anomalies. In this study we focus the Martian nightside where the ASPERA instrument on the Phobos-2 spacecraft and recently the ASPERA-3 instruments on the Mars Express spacecraft have measured the proprieties of escaping atomic and molecular ions, in particular O⁺ and O⁺ ₂ ions. We study the ion velocity distribution and how the escaping planetary ions are distributed in the tail. We also create similar types of energy-spectrograms from the simulation as were obtained from ASPERA-3 ion measurements. We found that the properties of the simulated escaping planetary ions have many qualitative and quantitative similarities with the observations made by ASPERA instruments. The general agreement with the observations suggest that acceleration of the planetary ions by the convective electric field associated with the flowing plasma is the key acceleration mechanism for the escaping ions observed at Mars.     

磁控热防护系统在天地往返运载器上的应用仿真

磁控热防护技术在高超声速领域显现出广泛的应用前景。考虑高超声速流动磁流体力学控制涉及的等离子体生成机制、多电离组分导电机理以及电磁流动能量/动量输运机制,通过耦合求解电磁场泊松方程和带电磁源项的高温热化学非平衡流动控制方程组,搭建了高超声速磁控热防护数值模拟平台。结合美国航天飞机"哥伦比亚"号（OV-102）近似外形和5种磁场配置方案,较为系统地开展了磁控热防护系统在高超声速"滑翔返回式"天地往返运载器上的应用仿真研究。结果表明：搭建的磁控热防护仿真平台具备偶极子磁场、均匀磁场、螺线管磁场及多个磁场组合条件下复杂外形飞行器气动热环境数值模拟能力,其校验结果与文献或飞行试验数据符合较好;采用合适的磁场配置能有效降低航天飞机的表面热流,显著改善了航天飞机的气动热环境,典型状态的表面热流下降25%以上;局部磁场方向与流动方向的夹角,在一定程度上决定了洛伦兹力的强度和方向,对磁控效果的影响明显。     

A quasi-one-dimensional model for hypersonic reactive flow along the stagnation streamline

This study proposes a quasi-one-dimensional model to predict the chemical non-equilibrium flow along the stagnation streamline of hypersonic flow past a blunt body. The model solves reduced equations along the stagnation streamline and predicts nearly identical results as the numerical solution of the full-field Navier-Stokes equations. The high efficiency of this model makes it useful to investigate the overall quantitative behavior of related physical-chemical phenomena. In this paper two important properties of hypersonic flow, shock stand-off distance and oxygen disso-ciation, are studied using the quasi-one-dimensional model with the ideal dissociating gas model. It is found that the shock stand-off distance is affected by both chemical and thermal non-equilibrium. The shock stand-off distance will increase when the flow conditions are changed from equilibrium to non-equilibrium, because the average density of the shock-compressed gas will decrease as a result of the increase in translational energy. For oxygen dissociation, the maximum value of its dis-sociation degree along the stagnation line varies with the flight altitude. It is increased at first and decreased thereafter with the altitude, which is due to the combination effect of the equilibrium shift and chemical non-equilibrium relaxation. The overall variation of the maximum dissociation is then plotted in the speed and altitude coordinates as a reference for engineering application.     

Vibrational relaxation of CO2(m,nl, p) in a CO2N2 mixture. Part 2: Application to a one-dimensional problem

In order to quantify the relaxation mechanism of CO₂(m, n^l, p), the vibrational level populations are calculated for a particular test-case: the vibrational relaxation of a CO₂N₂ mixture along the stagnation streamline in a reentry problem. The N₂ species is chosen as a collision partner because it is a component existing in numerous gaseous mixtures (cf. Part. 1). Excitation and deexcitation processes are taken following Nickerson and Herzfeld. The Navier-Stokes code CELHYO for the simulation of hypersonic laminar viscous flows in chemical and thermal nonequilibrium is used with a new one-dimensional approach, reduction of the Navier-Stokes equations along the stagnation streamline. Mass fraction and ‘vibrational temperature’ distributions of every vibrational level, considered as an independent chemical species, are presented for the two different CO₂N₂ mixture compositions. The validity of the usual assumptions for the vibrational mechanism is examined on the basis of the obtained results.     

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

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用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 范围,分别进行无烧蚀的壁面催化与非催化条件和石墨烧蚀条件下的流场计算,分析各类条件下的流场结构、流动分离特性以及流场热化学参数分布特点,研究壁面条件对流动特性的影响。结果表明：流动分离可能性和分离区范围随着压缩拐角斜面倾角增大而增大,随来流马赫数增大而减小;相对于低壁温条件,无烧蚀的辐射平衡壁温和壁面烧蚀条件下流动分离区增大,斜面上压力、摩阻和热流峰值点也有所后移。     

稀薄流航天器鼻锥迎风凹腔气动力和气动热性能研究

本文采用直接模拟蒙特卡罗(DSMC)方法,对高超声速稀薄流中航天器鼻锥迎风凹腔气动力与气动热性能进行了数值研究。得到了鼻锥外壁面、凹腔侧壁面以及凹腔底面的热流密度分布,分析了不同凹腔深宽比对鼻锥冷却效率以及凹腔腔体内气体参数的影响;以深宽比为1的凹腔为基准,研究了凹腔唇口钝化半径对航天器气动热与气动力的影响。数值结果表明,高超声速稀薄流中迎风凹腔能够降低鼻锥外壁面的热流密度;当凹腔深宽比达到1之后,凹腔侧壁面热流变化趋于一致,热流密度最低点的轴向位置不随深宽比改变,且凹腔底部热流很小;凹腔近底部气体均由稀薄流转化为连续流,腔内气体压力不断振荡;唇口钝化没有明显优势,虽然可以降低鼻锥峰值热流,但是会带来严重的气动力性能下降。     

典型尖楔结构燃气流热试验工况确定方法

在传统冷壁热流模拟方法的基础上,进一步提出以热壁温度及热流密度的时序变化曲线为控制目标的燃气流热试验工况确定方法,即利用壁温控制目标与实测值的偏差对热壁热流控制目标做一定修正,以尽可能消除和弥补前期试验误差,同时利用300K冷壁边界热流密度数据库插值迭代方法,快速确定一定气动热模拟所需燃气流温度,解决了沿飞行轨迹瞬态热试验技术难题之一。利用CFD数值模拟方法,建立了典型尖楔结构高/中温双路燃气流组合热试验300K冷壁边界热流密度数据库,并针对典型尖楔结构沿某飞行轨迹9个典型状态气动热模拟需求,确定相应双路燃气流热模拟参数。相关数值计算结果显示,驻点区域热流密度平均模拟偏差为4.5%,平板区热流密度平均模拟偏差为4.6%,两者最大模拟偏差均不大于8%,满足工程试验精度要求。同时,瞬态热分析结果显示第45s时,距驻点1mm处最大温度梯度达到21K/mm,距驻点10.1mm处最大温度梯度达到18K/mm,满足气动热大温度梯度效应需求。      

Experimental Investigation into Characteristics of Plasma Aerodynamic Actuation Generated by Dielectric Barrier Discharge

This article carries out synthetic measurements and analysis of the characteristics of the asymmetric surface dielectric barrier discharge plasma aerodynamic actuation.The rotational and vibrational temperatures of an N2 (C3Пu) molecule are measured in terms of the optical emission spectra from the N2 second positive system.A simplified collision-radiation model for N2(C) and N2+(B) is established on the basis of the ratio of emission intensity at 391.4 nm to that at 380.5 nm and the ratio of emission intensity at 371.1 nm to that at 380.5 nm for calculating temporal and spatial averaged electron temperatures and densities.Under one atmosphere pressure,the electron temperature and density are on the order of 1.6 eV and 1011cm-3 respectively.The body force induced by the plasma aerodynamic actuation is on the order of tens of mN while the induced flow velocity is around 1.3m/s.Starting vortex is firstly induced by the actuation;then it develops into a near-wall jet,about 70 mm downstream of the actuator.Unsteady plasma aerodynamic actuation might stimulate more vortexes in the flow field.The induced flow direction by nanosecond discharge plasma aerodynamic actuation is not parallel,but vertical to the dielectric layer surface.     

高超声速典型弹道下的壁板热气动弹性动力学分析

对高超声速流中带有热防护系统（TPS）的二维壁板进行了热气动弹性的双向耦合建模与分析,采用三阶活塞理论计算气动力,通过参考焓法获得气动热流,在有限差分法的基础上进行结构热传导计算,拟合了结构材料特性随温度退化的曲线,最后将气动热模块、气动弹性模块进行双向耦合以考虑气动热与结构形变之间的相互反馈,并在2种典型弹道状态下进行热气动弹性响应分析。结果表明,在X-34A的设计弹道下,双向耦合分析会引起更加剧烈的热应力与热弯矩的变化与较长的瞬态混沌过程。在FALCON弹道下,双向耦合得到的结果加热更为剧烈,而温度下降的过程也更快。对比2种弹道发现,长时间的高超声速飞行更容易引发颤振,而机动性较强的弹道面临的主要问题则是屈曲,需要考虑材料的应力及强度特性。同时说明了双向耦合策略对于现代飞行器在弹道状态下工作的热气弹响应分析的必要性。     

Numerical simulation of hypersonic thermochemical nonequilibrium flows using nonlinear coupled constitutive relations

To predict aeroheating performance of hypersonic vehicles accurately in thermochemical nonequilibrium flows accompanied by rarefaction effect, a Nonlinear Coupled Constitutive Relations (NCCR) model coupled with Gupta’s chemical models and Park’s two-temperature model is firstly proposed in this paper. Three typical cases are intensively investigated for further validation, including hypersonic flows over a two-dimensional cylinder, a RAM-C II flight vehicle and a type HTV-2 flight vehicle. The results predicted by NCCR solution, such as heat flux coefficient and electron number densities, are in better agreement with those of direct simulation Monte Carlo or flight data than Navier-Stokes equations, especially in the extremely nonequilibrium regions, which indicates the potential of the newly-developed solution to capture both thermochemical and rarefied nonequilibrium effects. The comparisons between the present solver and NCCR model without a two-temperature model are also conducted to demonstrate the significance of vibrational energy source term in the accurate simulation of high-Mach flows.     

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

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