双组元姿控发动机喷管化学反应流场数值模拟

本文对混合比为０．９、１．０、１．１三种状态下工作的双组元自然推进剂（肼／四氧化二氮）姿控发动机喷管内化学反应流动进行了数值模拟。数值模拟时采用了弱耦合点隐式方法的数值方法及肼／四氧化二氮的十二组分、十三个基元反应的有限速率化学反应模型。得到了三种混合比下反应流及混合比为１．０时冻结流发动机的推力和比推力、喷管中的流动参数及各组分的质量分数。分析表明，数值模拟的结果与理论分析一致，结果可靠。本文工作为姿控发动机的喷管设计提供了理论依据。     

再入气动理论和实验的评述

战略武器的发展,要求进一步提高再入飞行器的突防能力、生存力和命中精度。现代高级再入飞行器主要是高β(弹道系数)再入的弹道式再入飞行器和机动武再入飞行器。本文所讨论的再入问题,就是指高级再入飞行器再入地球的大气层时所遇到的问题,整个再入系统虽然还包括材料、结构、遥测、控制等其它方面。但是再入动力学和气动热力学是再入系统的最重要问题之一,本文着重讨论这方面的现状、存在问题和解决问题的技术途径。 再入气动力学和气动热力学的主要问题是烧蚀防热问题,气动力问题,粒子云侵蚀问题,滚动问题,机动再入问题和再入物理等问题。解决这些问题要采用理论和实验相结合的方法,依靠理论计算、地面模拟试验和飞行试验等手段,通过综合分析提高设计计算的精度,寻找有效的措施保证性能要求。     

高速高机动飞行器的横航向偏离预测判据分析

以提高高速飞行器静稳定设计精度为目标,对横航向小扰动方程解的稳定性进行了讨论,推导出了开环动态偏航稳定判据(CnβDYN)和闭环横向控制偏离判据(LCDP)。进而对翼-体组合体高超声速飞行器做了工程估算,利用所得气动数据对此外形进行了全面的横航向稳定性分析,指出了CnβDYN、LCDP与传统横航向稳定性判据的差异,传统的横航向判据不能计入两者的相互耦合效应,也无法预则进行滚转控制时引发的偏航失稳发散。     

扩散抛物化(DP)NS方程组的意义及其在计算流体力学中的应用

本文首先讨论扩散抛物化（DP）NS方程组的早期研究工作：它的提出、数学性质、意义和在CFD中的应用，然后讨论扩散抛物化理论的一些新发展。这些新发展是对NS方程组数值计算进行物理分析的基础上得到的，其中包括NS方程组差分计算时，粘性剪切流对网格间距和格式精度的要求；粘性项只保留剪切项的广义扩散抛物化（GDP）NS方程组，它的性质和应用。由于高Re数流动在NS方程组的差分计算中，网格Re数彼此相差悬殊的特点，产生了计算离散单元守恒方程组的新的算法思路，即离散流体力学（DFD）算法。在DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD格式的构造、它的优点和应用。并以超声速绕前后台阶流动为算例，来说明GDPNS方程组的用处和DFD算法的优点。DPNS方程组、GDPNS方程组、DFD算法是高智提出的，对这些问题，他的合作者从理论、算法、数值验证和某些应用又作了系统的研究。     

用于二维无粘跨音速流计算的隐式TVD有限体积法(英文)

本文结合Harten的TVD(Total variation diminishing)概念和有限体公式建立一种求解守恒型Euler方程组的有效方法,采用通量分裂的迎风格式对方程组的隐式算子部分离散化,得到系数矩阵是主对角线占优的线性方程组,采用线Gauss-Seidel迭代推进法求解,提高了数值收敛速度,给出了二维跨音速内流的计算结果。     

超声速气流中简化微小液滴横向喷射的汽化过程探索

为研究超声速气流中简化液滴的汽化过程问题,本文分析了两相流计算中已有的两相传热模型,并对简化液滴绕流开展数值计算.在来流Ma≤0.6的条件下,数值计算得到的简化液滴-气流之间的传热速率与已有模型得到的结果相一致,而在来流Ma≥0.9的条件下,数值计算得到的简化液滴-气流之间的传热速率与已有模型得到的结果存在很大偏差.由此建立了考虑简化液滴与气流相对超声速相互作用的两相传热模型.进一步,采用Charles B.Henderson阻力系数关系式与新建立的传热模型,对不同直径简化液滴的运动与汽化开展工程计算.在来流2.7Ma的二维平板超声速流场中选取一个截面,作为气相流场,结果显示,(1)简化液滴与主气流存在相对超声速作用.当简化液滴直径dk≤0.12mm时,作用区域约为0.1m～0.4m,当dk>0.12mm时,作用区域明显增大,(2)简化液滴的穿透尺度不超过0.011m/m(深度/长度),时间尺度约为0.28ms～3ms,(3)简化液滴完成汽化的空间尺度约为0.1m(dk>=0.03mm)、0.45m(dk>=0.05mm)与1.24m(dk>=0.075mm),而当dk>0.09mm时,简化液滴完成汽化的空间尺度则大于1.9m.使用考虑简化液滴与气流相对超声速相互作用的两相传热模型与使用传统的传热模型对简化液滴的运动轨迹没有影响,而对简化液滴的汽化过程有较大影响.     

改进油流显示法研究层流和湍流下绕椭球分离流

用改进的油流显示法研究和对比了层流和湍流状态下绕椭球流动的表面摩擦力线结构 .实验中观察到 ,椭球表面摩擦力线在层流和湍流状态下具有不同的拓扑结构。只看流动对称面一侧 ,在层流状态 ,Re=1 .4× 1 0 6,迎角为 3 0°和 2 0°时都有三条分离线 ;而湍流状态同样实验条件下只观察到两条 ,并且分离线的位置推向下游 ,二次分离线长度也大大缩短。实验还观察了物体表面颗粒状突起对表面摩擦力线的影响 .层流状态下 ,雷诺数为 1 .4× 1 0 6时 ,分离线附近的颗粒状突起会显著影响分离线的形状 ,当雷诺数减为 0 .9× 1 0 6时 ,颗粒对分离线只有微弱的影响 .实验观察到的颗粒影响局限在壁面附近。本文用涡相互作用对之作了解释     

基于简化Boltzmann方程的超声速流稀薄效应分析

从近连续流到近自由分子流,计算分析超声速流场的稀薄效应。模型Boltzmann方程先后经简化分布函数和离散速度坐标变换后,采用一个隐式通量修正二阶迎风TVD格式差分求解。壁面取漫反射气体分子模型。在Knudsen数为0.001,0.01,0.1,1.0,10等5种情形下,数值模拟二维圆柱氩气体绕流,观察到了不同的波系及尾迹结构。总阻力系数计算值同实验数据基本吻合,碰撞项Shakov动力学模型反映的流场细节略好于修正BGK模型。稀薄效应趋于弱化激波等超声速流动结构,物体影响域增大。     

航天飞机的粘性干扰效应——航天飞机的空气动力学问题之二

航天飞机在55公里高度以上飞行时,处于高马赫数,中、低雷诺数的粘性干扰区(粘性干扰参数(?)=0.005～O.08)。传统的无粘流/边界层修正法(3DV)已不适于计算无粘/有粘强干扰的情况,而要采用全粘性激波层的 PNS 法。本文简要的介绍了美国航天飞机研制过程中采用的3DV 和 PNS 的计算方法,对航天飞机简化外形和钝锥体的粘性干扰效应进行了分析和讨论。研究结果表明,PNS 法是解决航天飞机高超声速粘性干扰效应的较为有效的数值模拟方法。     

求解气体动力学方程的推进迭代方法

本文由三部分组成,第一部分考察隐式推进迭代求解气体动力学方程的基本方面;第二部分介绍带压力松弛因子的推进迭代方法;一些数值算例在最后一部分给出。计算表明,此方法比目前一般方法具有更高的收敛速度,对于不同问题具有更强的适应能力。