高温空气中主要分子带系的F—C因子

从计入分子振转相互作用的电子态波函数出发，按自编程序计算了高温空气中主要分子带系Ｏ２（Ｓ－Ｒ），Ｎ２（１＾＋），Ｎ２（２＾＋），Ｎ２＾＋（－），ＮＯ（β），ＮＯ（γ）的Ｆ－Ｃ。从而得出：当激波管低压端的空气压强为１３３３Ｐａ时，入射激波及反射激波后高温空气的温度与激波强度的关系；各主要分子带系转动态分布与温度的关系；高温空气中主要分子带系Ｆ－Ｃ因子与激波强度的关系。计算结果甚好，好既适用于高温，也     

用解析离散法构造EWNO—FCT格式

本文利用了解析离散法的导数计算公式构造高精度ＷＥＮＯ格式，并引入ＦＣＴ通量修正技术以消除非物理解振荡，通过激波－边界层干扰，激波－旋涡相干结构的数据模拟证明方法既能够有效地捕捉激波间断，又对边界层分离及旋涡结构有很好的分辩率，是一类有效的高精度计算格式。     

化学非平衡高超声速电离流动数值研究

本文对非平衡化学区应高超声速绕流零攻角轴对称钝头体电离空气粘性流动进行了研究，通过带有化学源项的Ｎａｖｉｃｒ－Ｓｔｏｋｅｓ方程对流动作数值模拟，总的连续方程由单个组元守恒方程代替。高温电离空气中存在７个主要的组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ，和ｅ￣－。化学动力学模型由７个组元和这些组元之间６个可能的反应来描述。发展了一种数值方法以求解这组完全耦合的与时间有关的偏微分方程组。利用一类新的有限差分方法来求解这组以守恒形式表示的控制方程。这是一种以ＴＶＤ格式来离散粘性流通量的方法。在数值实验中采用了ＴＶＤＹｅｅ－Ｈａｒｔｅｎ格式，对全部算例都是用激波捕捉法。关于化学组元，我们试图处理完全催化壁、非催化壁和有限催化壁等三种不同的表面边界层条件。计算是对ＲＡＭＣ－Ⅱ试验条件进行的，得出了７１ｋｍ高空情形的典型流场数值解。对表面压力，Ｓｔａｎｔｏｎ数和热流率给出了化学反应气体结果与完全气体结果的比较。计算的流场电子密度与ＲＡＭＣ－Ⅱ实验数据进行了比较，计算结果与实验数据符合是相当好的。     

激波运动与反射流场的ENO算法

在ＴＶＤ格式的基础上，以修正限制函数的方法，将通常的ＴＶＤ格式改造为性能更为优越的ＥＮＯ格式，并利用这种有限差分格式求解欧拉方程和薄层Ｎ－Ｓ方程。分别数值模拟了激波管中运动激波与三种楔形物体相互作用流场；给出了激波在不同物体上的反射及波涡干扰过程中非定常流场的变化情况；算出的波系结构在与ＣＡＲＤＣ小型激波管中做的相关实验结果一致。     

可压缩粘性流绕翼型的N─S方程自适应网格解

从Ｎ－Ｓ方程出发,采用ＬＵ－ＡＤＩ隐式分解方法求解绕翼型的可压缩粘性流动。湍流模式采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式和ｑ－ω二方程微分模式以研究湍流边界层的非平衡效应;为更好地捕捉激波,对网格进行了自适应调整。通过对ＮＡＣＡ００１２和ＲＡＥ２８２２两种翼型在亚音速和跨音速不同马赫数和雷诺数下进行大量计算,表明该方法对改进湍流计算、提高激波分辨率有较好的效果     

TATB含铝炸药粉反应特性的激波管研究

用自制光谱探测技术和回收技术，在爆炸激波管中研究了铝粉掺入了ＴＡＴＢ炸药粉中的快速反应特性。研究表明：随着铝粉含量Ｎ的增加，铝粉的氧化反应和氢化反应加剧，在Ｎ－３３％时达到最大，继续增加铝粉含量，铝粉反应明显减弱。铝粉的颗粒度为４４μｍ较好。     

用解析离散法构造 WENO-FCT 格式

本文利用解析离散法的导数计算公式构造高精度ＷＥＮＯ格式,并引入ＦＣＴ通量修正技术以消除非物理解振荡。通过激波－边界层干扰,激波－旋涡相干结构的数值模拟证明方法既能够有效地捕捉激波间断,又对边界层分离及旋涡结构有很好的分辨率,是一类有效的高精度计算格式。     

《航空精密制造技术》2001年总目次

☆综述☆迎接新世纪的曙光ＮＯ．１振动主动控制中的三个问题ＮＯ．１高速切削及其关键技术的发展现状ＮＯ．２超声波马达技术及其应用ＮＯ．２站在新的起点上ＮＯ．５☆精密加工☆纳米级显微硬度试验研究ＮＯ．１微镜单元的设计与加工工艺研究ＮＯ．１陶瓷喷涂层精密镜面磨削技术的实验研究ＮＯ．１脆性材料塑性域的超精密加工方法ＮＯ．２接触疲劳磨损试验中球形磨损微粒形成机理研究ＮＯ．２非球面加工先进技术ＮＯ．３薄壁零件的超声振动精密切削研究ＮＯ．４新型精密狭缝的设计ＮＯ．４第二届欧洲精密工程和纳米技术国际会议ＮＯ．４滚…     

管内正激波/湍流边界层干扰的数值模拟

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。计算与实验结果的比较表明，方法可准确地预测激波结构，激波与边界层干扰区流动特征，激波位置、激波前马赫数和壁面压力分布等。     

激波风洞中三角翼体模型的大攻角气动力测量和局部方法的应用

在激波管风洞中，用模型自由飞方法测量了两种不同迎风面三角翼体模型的大攻角气动力。应用局部方法以一种迎风面外形的实验结果为依据，估算了另一迎风面外形的的气动特性，结果与实验值吻合较好。另外，本文还进一步估算了航天飞机外形的气动特性，结果与实验值吻合较好，另外，本文还进一步估算了航天飞机外形的气动特性，与美国ＣＡＬＳＰＡＮ公司４８英寸（１．２米）激波风洞中的实验值相比，结果也令人满意，从而表明，局部方     

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

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用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算法与流场传热分析

在非结构网格下,研究了热力学碰撞传能的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数去分析飞行器壁面的传热效果,所有这些分析对高温部件的热防护设计十分有益.      

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.     

NNW-HyFLOW高超声速流动模拟软件框架设计

NNW-HyFLOW在国家数值风洞（NNW）工程支持下,由风雷开源软件提供基础框架支撑,将建设成为基于结构/非结构混合网格,面向高超声速应用领域的国产自主工业CFD软件,具备高温气体热化学非平衡效应模拟及其相关气动力、气动热和气动物理特性计算分析等主要特色功能。本文从框架设计、数据结构、耦合方法、并行计算方法以及接口设计等方面对软件设计思想和框架特点进行了介绍,给出了求解器采用的理论模型、核心数值算法及其实现方法,结合HEG风洞试验、RAM-C飞行试验、Electre飞行试验以及航天飞机OV102飞行试验等典型算例开展了数值模拟。研究表明：NNW-HyFLOW具有底层代码复用、功能兼容性好、扩展能力强和接口灵活等优点,其当前测试版本已经具备了较好的高超声速非平衡流动数值模拟能力,在热化学非平衡效应及其影响的气动力特性、气动热环境和等离子体分布特性预测与评估方面,具有较高的数值计算精准度,初步满足了高超声速复杂飞行器高温非平衡流动数值模拟的需求。     

不同热化学非平衡模型对高超声速喷管流场影响的数值分析

采用轴对称NS方程,数值研究不同热化学非平衡模型对高超声速喷管流场的影响,包括:(1)不同组元数的高温空气模型(5组元、7组元1、1组元)的比较;(2)热力非平衡(双温度)的化学动力过程与热力平衡(单温度)的化学动力过程的比较。计算结果表明,高焓风洞实验条件下喷管流场处于热力和化学都是非平衡的状态。在计算条件下,数值模拟以采用7组元或11组元的热化学非平衡模型为宜。     

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

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

姿控发动机内流场及高空羽流流场的DSMC一体化数值模拟

通过引入“平衡抽样方法”有效解决了DSMC求解近连续流区所遇到的碰撞顶计算效率问题;对于过渡流及自由分子流区,构造了适合于DSM仿真的分子碰撞传能模型并配置了一种光滑过渡以减少误差的权函数。在此基础上首次实现了对姿控发动机喷管内流场及包括倒流区在内的高空弱流流场的DSMC方法一体化数值模拟。     

Numerical analysis of hypersonic thermochemical non-equilibrium environment for an entry configuration in ionized flow

A theoretical methodology for thermochemical non-equilibrium flow combing with the HLLC (Harten-Lax-van Leer Contact) scheme was applied to study the hypersonic thermochemical non-equilibrium environment of an entry configuration in ionized flow. A two-temperature controlling model was utilized and the Gupta’s 11 species (N₂, O₂, NO, O, N, NO⁺, N₂⁺, O₂⁺, N⁺, O⁺, e^?) thermochemical non-equilibrium model was taken. Firstly, numerical calculations of hypersonic thermochemical non-equilibrium environments for different aerodynamic shapes were carried out to verify the reliability of the method above. Then, the method was used to research the effects of ionization and wall catalysis on the hypersonic thermochemical non-equilibrium environment of the entry configuration in ionized flow. The shock stand-off distance can be reduced by thermochemical reactions but doesn’t continue to decrease significantly when ionization occurs. The shock stand-off distance calculated by the 11 species model is 4.2% smaller than that calculated by the 5 species (N₂, O₂, NO, O, N) thermochemical non-equilibrium model without considering ionization. Ionization reduces wall heat flux but increases wall pressure a little. The effect of ionization on aerothermal loads is greater than that of aerodynamic loads. The thermochemical reactions of electrons and ions catalyzed at the wall increase wall heat flux significantly but make a small change in wall pressure. The maximum wall heat flux obtained by only considering the electrons and ions catalyzed at the partially catalytic wall condition is 11.8% less than that calculated at the super-catalytic wall condition.     

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

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[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.