钝锥高超声速绕流空间分离形态数值模拟及拓扑分析

本文通过数值模拟和理论分析，系统研究了钝锥高超声速有迎角绕流时垂直于物面的截面流态沿轴向的演变过程。得出了高超声速情况下，背风分离区的流态在大多数情况下貌似对称，其实是非对称的这一重要结论。考察了周向网格疏密程度对流动结构分辨率的影响，发现对分离流动的数值模拟，必须在分离区内布置足够的网格才能给出正确的结果。     

高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

网格对高超声速钝头体表面热流数值模拟结果的影响

高超声速钝头体热流的精确预估，对高超声速飞行器热防护的设计具有重要意义，而网格划分尺度对预估飞行器表面热流精度具有重要的影响。应用理论分析及数值模拟方法，提出了一种基于壁面网格雷诺数及基于钝体特征长度的来流雷诺数网格划分方法，并给出了壁面网格雷诺数的取值范围。应用所提出的网格划分方法，对不同来流条件的高超声速半圆柱及球头表面热流进行了数值模拟，结果表明，应用所提出的网格划分方法，及给出的可供参考的壁面网格雷诺数的取值范围划分网格，在满足高超声速钝体表面热流精度要求的同时，有利于网格合理分布并提高数值计算效率。     

基于改进SST模型的分离流动数值模拟

对SST湍流模型进行了改进,并通过对典型分离流动进行数值模拟,来检验和提高模型预测分离流动的能力.基于亚声速分离流动,提出了提高雷诺应力的模拟精度和分离流修正的改进方法,并进行了对比研究得出结论;在亚声速分离流动分析结论基础上,采用了可压缩性修正方法,开展了跨声速、超声速和高超声速激波/边界层干扰分离流动的数值模拟研究.结果表明:提高雷诺应力的模拟精度和采用分离流修正明显地提高了SST湍流模型对分离流动的模拟能力;适当地可压缩性修正对超声速和高超声速分离流动的计算精度有改善作用.     

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

立楔诱导的高超声速流动分离数值分析

高超声速飞行器进气道等关键部件引起的激波与边界层相互作用将导致流动分离,从而改变当地压力分布与局部受热情况,影响飞行稳定性与飞行安全,因此需要对高超声速流动的分离现象进行细致研究。采用高精度5阶特征型WENO格式与3阶TVD型Runge-Kutta方法,求解三维Navier-Stokes方程,对立楔诱导的高超声速激波与边界层相互作用引起的分离流动流场结构进行了细致的数值模拟与分析。结果表明,5阶特征型WENO格式分辨率远高于类TVD格式;Ma=6时得到清晰的激波结构、分离涡结构及其演化过程和壁面极限流线的拓扑结构,证明了WENO格式应用于高超声速分离流动的可行性与高分辨率;对不同来流Mach数的对比证明Mach数的增大抑制流动分离,导致分离涡减小。     

高超声速验证飞行器助推分离段流场数值研究

为了模拟有相对运动的多体非定常绕流，发展了基于弹簧近似的非结构动网格方法及耦合动网格的Euler解算器，作为方法验证，模拟了俯爷振动的NACA0012翼型绕流，计算结果与实验及文献结果非常接近。采用二维近似外形，对高超声速验证习行器助推分离段流场进行了数值研究，得到了不同时刻由于多体相对运动形成的干扰流场结构以及分离过程的气动力参数。     

稀薄流高超声速钝化前缘逆向射流流场特征研究

为准确预测稀薄过渡区逆向射流对稀薄大气的干扰流动特征，本文采用直接蒙特卡洛算法（DSMC），在自由流马赫数为7的情况下，对稀薄流中逆向射流干扰下的高超声速流动进行了数值模拟。计算中考虑了60～90 km高度的典型大气环境，研究了压比和自由流克努森数对压力、剪切力和热流密度的影响规律。结果表明，逆向射流改变了流场结构，且具有明显的热防护性能，但减阻效果不明显。逆向射流对高超声速稀薄来流的影响随压比的增加明显增强，形成的马赫盘大小也随之增加。高度的增加伴随稀薄效应不断增强，自由流克努森数增大，导致弓形激波厚度随高度增加而增加，壁面处的回流区范围不断减小至消失。这项研究初步揭示了稀薄流中逆向射流/高超声速流动相互作用的机制。     

AUSM+-up格式在无网格算法中的推广

将AUSM -up格式推广运用到了无网格算法中.在点云离散的基础上,运用曲面拟合和AUSM -up格式求得数值通量;采用四步龙格-库塔方法进行时间推进,并引入当地时间步长和残值光顺等加速收敛措施.为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型跨声速流动、某三段翼型低速绕流、平行错位NACA0012双翼超声速、高超声速流动进行了数值模拟.     

γ-Re_θ转捩模型在高超声速下的应用及分析

为了初步研究γ-Reθ转捩模型在高超声速领域的适用性,在自行开发的CFD程序中添加了该模型,利用T3系列平板算例对该模型的实现进行了验证;针对高超声速流动,选取压缩面绕流和双圆锥扰流等算例,分别利用该转捩模型、全层流模拟、以及传统的湍流模型进行了数值模拟计算,并与实验结果进行对比。计算结果表明,尽管该转捩模型是在低速流动的基础上发展的,它仍能准确预测高超声速流动下的转捩现象。     

过渡领域高超声速圆柱绕流直接模拟

本文用直接模拟蒙特卡罗方法模拟了再入速度为７．５ｋｍ／ｓ，高度分别为８０、８５、９０ｋｍ，物壁是完全漫反射有限催化壁的圆柱绕流。在８０ｋｍ考察了壁面催化率的影响。     

细长体大攻角绕流的数值模拟

用数值模拟方法研究了细长体大攻角的绕流特性（求解在周向３６０°范围内进行）。数值模拟的出发方程是三维可压流的全Ｎ－Ｓ方程,格式为二阶迎风ＴＶＤ格式。通过数值模拟研究了附近,攻角１９°、３２．５°和５０°状态下,细长体绕流场中的激波、分离流动及表面压力分布等重要特征。     

垂直于物面的横截面上流态的拓扑

对于旋成体有迎角绕流，本文研究了垂直于物面的横截面上的流态结果。利用解析方法，给出了物面轮廓线上、迎风和背风对称线上以及全横截面上截面流线方程奇点的分布规律以及奇点沿流动方向的演变规律，从而可定性地确定截面流态的结构。通过求解ＮＳ方程数值模拟了钝锥有迎角高超声速绕流流场，计算获得的垂直于物面的横截面流线的拓扑结果和理论分析完全一致。     

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

A scaling procedure for measuring thermal structural vibration generated by wall pressure fluctuation

This paper attempts to develop a scaling procedure to measure structural vibration caused simultaneously by wall pressure fluctuations and the thermal load of hypersonic flow by a wind tunnel test. However, simulating the effect of thermal load is difficult with a scaled model in a wind tunnel due to the nonlinear effect of thermal load on a structure. In this work, the temperature variation of a structure is proposed to indicate the nonlinear effect of the thermal load,which provides a means to simulate both the thermal load and wall pressure fluctuations of a hypersonic Turbulent Boundary Layer(TBL) in a wind tunnel test. To validate the scaling procedure,both numerical computations and measurements are performed in this work. Theoretical results show that the scaling procedure can also be adapted to the buckling temperature of a structure even though the scaling procedure is derived from a reference temperature below the critical temperature of the structure. For the measurement, wall pressure fluctuations and thermal environment are simulated by creating hypersonic flow in a wind tunnel. Some encouraging results demonstrate the effectiveness of the scaling procedure for assessing structural vibration generated by hypersonic flow. The scaling procedure developed in this study will provide theoretical support to develop a new measurement technology to evaluate vibration of aircraft due to hypersonic flow.     

非对称腔体合成双射流矢量特性数值研究

利用合成射流全流场计算X-L模型,对非对称腔体合成双射流激励器的射流矢量特性进行了数值模拟,分析了激励器不同腔体体积比、不同振动膜速度、不同频率下的合成双射流流场。研究结果表明,当激励器两腔体不对称时,所形成的合成双射流不再垂直向下游流动,而是会发生偏转,即射流具有矢量特性,偏转角度的大小可以通过激励器工作条件进行调节。其机理是振动膜对大小不同的两腔体的相对压缩量不同,使得两出口处所形成射流的能量和低压区不再对称,因此导致合成双射流向低能量(大腔体)一侧偏转。振动膜振幅和频率的提高对低压区大小分别起到增加和减弱的作用,所以偏转角度又随着振动膜最大速度的增加而增大,随工作频率的增加而减小。     

推力耦合的高超声速飞行器气动伺服弹性研究

对于采用吸气式超燃冲压发动机的高超声速飞行器,其发动机推力可能与机身弹性发生耦合影响,从而引起所谓的推力耦合气动伺服弹性(ASE)问题。为对其耦合原理及影响进行研究,以简化的飞行器纵向模型为对象,考虑结构弹性、非定常气动力、冲压发动机以及控制系统之间的相互耦合作用,建立了推力耦合的高超声速飞行器气动伺服弹性问题的一般建模框架和分析流程。采用牛顿冲击理论计算高超声速非定常气动力,基于准一维流动假设分析发动机性能。算例结果表明,考虑发动机推力的耦合影响后,飞行器的短周期特性和气动伺服弹性特性均有明显改变,气动伺服弹性稳定裕度下降可达16%,应当引起飞行控制系统设计部门的重视。     

一种前体加宽型高超声速进气道试验方案研究

 根据矩形截面高超声速进气道前体的流动特征,对一种前体加宽型高超声速进气道试验方案开展了数值仿真及高焓风洞试验研究。首先,对不同前体宽度的高超声速进气道开展了三维数值仿真研究,结果显示:随着前体宽度的增加,进气道的流量系数和静压比逐渐增加,而总压恢复系数和隔离段出口马赫数逐渐减小,表现为先急后缓,且当来流马赫数和来流攻角变化时依旧保持上述变化规律。其次,对前体加宽型高超声速进气道试验方案开展了高焓风洞试验研究,结果表明:加宽前体可有效地提高进气道的流量系数,较为真实地反映此类进气道的流动特征,试验结果与数值仿真结果吻合较好。考虑到进气道性能参数随前体宽度变化规律表现为先急后缓,建议在试验条件下前体宽度比取0.5~0.8之间较为适宜。     

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.     

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

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