A hybrid CFD/characteristics method for fast characterization of hypersonic blunt forebody/inlet flow

A hybrid CFD/characteristic method(CCM) was proposed for fast design and evaluation of hypersonic inlet flow with nose bluntness, which targets the combined advantages of CFD and method of characteristics. Both the accuracy and efficiency of the developed CCM were verified reliably, and it was well demonstrated for the external surfaces design of a hypersonic forebody/inlet with nose bluntness. With the help of CCM method, effects of nose bluntness on forebody shock shapes and the flowfield qualities which dominate inlet performance were examined and analyzed on the two-dimensional and axisymmetric configurations. The results showed that blunt effects of a wedge forebody are more substantial than that of related cone cases. For a conical forebody with a properly blunted nose, a recovery of the shock front back to that of corresponding sharp nose is exhibited, accompanied with a gradually fading out of entropy layer effects. Consequently a simplification is thought to be reasonable for an axisymmetric inlet with a proper compression angle, and a blunt nose of limited radius can be idealized as a sharp nose, as the spillage and flow variations at the entrance are negligible, even though the nose scale increases to 10% cowl lip radius. Whereas for two-dimensional inlets, the blunt effects are substantial since not only the inlet capturing/starting capabilities, but also the flow uniformities are obviously degraded.     

Numerical flow visualization of a Single Expansion Ramp Nozzle with hypersonic external flow

Numerical simulation of scramjet asymmetric nozzle flow is carried out to visualize and investigate the effects of interaction between engine exhaust and hypersonic external flow. The Single Expansion Ramp Nozzle (SERN) configuration studied here consists of flat ramp and a cowl with different combinations of ramp angle and cowl geometry. UsingPARAS 3D, simulations are performed for a free stream Mach number of 6.5 that constitutes the external flow around the vehicle. Appropriate specific heats ratio has been simulated for the jet and free stream flow. External shock wave due to jet plume interaction with free stream flow, the internal barrel shock wave and the shear layer emanating from the cowl trailing edge and sidewalls are well captured. Wall static pressure distribution on the nozzle ramp for different nozzle expansion angles has been computed for both with and without side fence. Axial thrust and normal force have been evaluated by integrating the wall static pressure. Effect of cowl length variation and side fence on the SERN performance has also been studied and found to be quite significant. Based on this study, an optimum ramp angle at which the SERN generates maximum axial thrust is obtained. SERN angle of 20° was found to be optimum when the flight axis coincides with nozzle axis.     

Aerodynamic design of a supersonic three-dimensional inlet

The results of designing and numerical gas-dynamic modeling a supersonic three-dimensional inlet of a new type are considered. A ramp of external compression of this inlet is the V-shaped body forming an initial plane oblique shock wave and a subsequent isentropic compression wave. The inlet incorporates an entrance section of internal compression, where also a plane oblique shock wave and a subsequent isentropic compression wave are formed by a cowl. The designed three-dimensional inlet has small inclination angles of compression surfaces, which ensures its low wave drag. According to the estimates of inlet efficiency in terms of the compression ratio and the total pressure recovery factor, it is close to the optimal two-dimensional shocked inlet of external compression considered by Oswatisch as well as Petrov and Ukhov. The flow in the inlet was computed with the use of the Euler and Navier — Stokes codes provided by the commercial package “FLUENT”. The flow in the inlet throat in the design regime computed under the inviscid flow approximation is uniform. The most substantial effect of the flow viscosity in this regime manifests itself in the interaction of the shock wave from the cowl with the boundary layer on the V-shaped compression body in the inlet internal duct. According to computed data, the boundary layer separation does not occur in this case; however, due to viscosity effects, reflected shock waves are formed here which results in significant deviations of flow structure as compared to the computed inviscid flow.     

不同来流马赫数下激光能量形成高超声速进气道虚拟唇口的特性

用基于三维非定常可压缩雷诺平均Navier-Stokes方程的有限体积法计算了马赫数低于设计值6时一种高超声速进气道的性能参数,发现其性能存在明显下降。为提升进气道性能,将功率为15kW的激光能量注入进气道固体唇口前的流场中,形成虚拟唇口,马赫数为4.5,5.0和5.5时,计算得到来流捕获率分别提高了34%,20.6%和15.6%。绘制了不同马赫数下来流捕获率达到峰值时的流场压强云图,说明了虚拟唇口的特性及形成机制。结果表明:来流马赫数越低,来流捕获率越小,但相对于无能量注入时的来流捕获率的提升程度越明显;在不同来流马赫数条件下,通过改变激光能量引致的激波结构和位置,可达到最优状态,即激波与进气道前缘斜激波相交后的透射波打在进气道肩部位置的状态。     

Design of supersonic three-dimensional inlets using two-dimensional isentropic compression flow

The design of supersonic three-dimensional inlets using the V-shaped body forming a two-dimensional flow including an initial oblique shock wave and a subsequent isentropic compression wave is considered. Such a flow appears attractive for inlets design due to a possibility of obtaining high compression levels of external flow over the inlet ramp with small losses of the total pressure. Numerical computations of the flows around the designed configurations were carried out in design and off-design regimes using Euler code. The flow structure was identified, the aerodynamic characteristics of the inlets were determined. The investigation covers the range of supersonic speeds corresponding to the freestream Mach numbers M_∞= 1.8−2.5.     

Design of cylindrical air inlets for high flight speeds from a combination of plane flows

In the present paper, we discuss a procedure for designing of cylindrical air inlets for high flight speeds with the use of V-shaped bodies for forming a plane flow with an initial oblique compression shock. In design regime, characteristics of such air inlets can be obtained by means of simple calculations performed in a broad range of governing parameters. The difference between the performance characteristics of a typical cylindrical inlet in design and off-design flow conditions was elucidated with the help of 3D numerical calculations.     

高超声速进气道飞行器一体化设计技术的发展

对高超声速进气道与飞行器一体化设计技术和发展进行了研究,包含轴对称压缩、二维压缩、侧压以及内压缩进气道在高超声速飞行器上的典型布局设计方案.对高超声速可调进气道类型进行了概述,基于轴对称类型调节、二元平面类型调节以及三维内转调节进气道的典型案例给出了其各自的设计特点,并进一步对宽域飞行和组合动力飞行器采用的多通道可调节...     

高超声速平板边界层/圆柱粗糙元非定常干扰

以高超声速表面湍流控制为应用背景,平板/粗糙元干扰流动为模型,采用大涡模拟方法研究粗糙元流场干扰作用机理.分析粗糙元外形特征对于流动稳定性影响,给出其引起的流动表面参数的变化规律.结果显示超声速边界层在粗糙元作用下产生强逆压梯度并发生分离,粗糙元高度对高位自由剪切层失稳有明显影响,低粗糙元干扰下游流动稳定性,而高粗糙元剪切层发生流向失稳,形成涡串结构;同时粗糙元干扰导致下游摩阻和热流系数较平板略低,可能应用在进气道降热和减阻中.     

转捩与湍流对激波边界层干扰及底部流动结构的影响

为研究转捩与湍流对激波边界层干扰及底部流动结构的影响,文章选取了二维与三维高超声速双斜面进气道模型与大钝头着陆器模型,并使用γ-Re_θ转捩模型开展数值模拟研究.研究表明,对于二维进气道模型,随着前缘钝度的增加,激波边界层干扰位置前移,分离区变大,与层流流动情况相比,有转捩流动发生时,激波边界层干扰位置后移,同时分离流动强度变弱,分离区缩小;对于三维进气道模型,其拐角附近的分离流动呈现明显的三维特征,转捩流动也存在三维流动结构,与静风洞状态相比,噪音风洞状态下,有转捩流动发生,对壁面热流影响较大,对激波系影响很小.对于着陆器模型,底部流动发生转捩,使得底部流动由不稳定非定常的流动结构变为稳定定常的流动结构,这有益于姿态控制设计.      

隔离段激波串流场特征的试验研究进展

高超声速推进技术是国际前沿研究, 其中双模态超燃冲压发动机的发展受到极大关注. 作为超燃冲压发动机的重要部件, 隔离段对发动机的性能和高超声速飞行的实现至关重要, 其中所涉及的流动机理问题也极为复杂. 自从高超声速飞行的概念被提出和论证以来, 相关的理论、试验和仿真研究不断取得进展, 但是对其中的机理问题研究仍有待进一步深入. 本文将从试验研究的角度回顾并综述近年来超燃冲压发动机隔离段的研究进展, 结合精细流动测试技术(Nano-tracer Planar Laser Scattering, NPLS)的发展分析了隔离段流场特征, 包括了激波串流场复杂的三维时空结构特点、湍流特性、非线性迟滞运动、不启动流场特征以及激波前缘检测等. 从风洞设备、隔离段设计、测试技术等方面对隔离段的试验研究进行了分类比较和论述, 对今后隔离段试验研究提出了建议.     

局域能量注入对高超声速进气道流动影响研究

研究了局域能量脉冲注入条件下高超声速进气道流场的扰动情况。数值模拟采用三维雷诺平均N-S方程,分别利用ROE格式和二阶中心格式对对流通量和粘性通量进行离散处理;用高斯-赛德尔隐式格式对方程进行时间推进求解,采用k-ε两方程模型用于湍流的数值模拟。开展了高超风洞平板流场能量注入实验,获取了高速纹影图像,并对实验结果与计算结果进行比较。研究表明,能量注入产生的冲击波能与高超流动产生的斜激波发生强烈干扰,脉冲能量的引入可能引起高超声速进气道流量俘获率产生剧烈震荡,从而导致进气道流场的性能急剧下降。     

Gas-dynamic problems in off-design operation of supersonic inlets (<Emphasis Type="Italic">review</Emphasis>)

Modern concepts of operation of supersonic inlets of high-velocity air-breathing engines are analyzed. It is demonstrated that the flow in the engine duct becomes extremely complicated in off-design modes of inlet operation, which can lead to unpredictable consequences, in particular, to inlet unstart. The term “inlet unstart” is considered in the present paper as a synonym of the absence of theoretical understanding and prediction of gas-dynamic phenomena. Various approaches are proposed to ensure self-regulation of the inlet-combustor system for air-breathing engines. Possible directions of further research are indicated for the purpose of stable operation of inlets in a wide range of flight conditions.     

标模外形化学非平衡流场数值模拟

利用三维并行计算代码求解Navier-Stokes方程,数值模拟标模(ELECTRE)化学非平衡绕流,研究真实气体效应对标模气动热特性的影响,反应模型为Dunn和Kang的7组元7反应化学动力学模型.利用典型弹道点的飞行试验数据验证化学非平衡流计算程序的可靠性.在此基础上,研究不同壁面催化条件下攻角和高度变化对热流的影响.计算表明:真实气体效应主要发生在物面附近很薄的激波层内,并使激波脱体距离减小;完全催化壁驻点热流值高于非催化壁热流值;随着攻角增大,热流分布差异明显,而且攻角越大时,物面电子数密度越小;飞行高度越高,O₂和N₂离解程度越低,驻点热流越低.     

On the question of starting conditions for frontal axisymmetric inlets tested in hot-shot wind tunnels

The work presents the results of an analysis of starting conditions for some frontal axisymmetric inlets of internal compression tested at freestream Mach numbers М = 3?8.4 in the hot-shot wind tunnels based at Khristianovich Institute of Theoretical and Applied Mechanics (ITAM). The results of these inlets test are compared with the data of numerical computations of inviscid, laminar, and turbulent flows carried out by the pseudo-unsteady method. There were determined the inlet throat areas limiting either with regard to the inlet starting or with regard to providing the maximally possible degree of geometric compression of the inlet-captured supersonic airstream at its deceleration in the already started inlet. Reshaping of computed flow patterns in the inlets depending on the variation of the minimal cross section of the inlet internal duct is analyzed.     

Acoustics of an obstacle inside a reactive silencer

The transmission and reflection of sound in a cylindrical duct containing several discontinuities is investigated. A building-block method, which gives the transmission and reflection matrices for a complex system from those of the parts, is applied to bifurcations, sudden area changes with or without extended inlets, and spherical obstacles (which may be lossy). In some cases the solution can be interpreted in terms of multiple reflections. When the lengths between the discontinuities are small it is important to include also non-propagating modes, and this is especially true when the sudden area change is obtained from the area change with extended inlet in the limit of vanishing inlet. For an expansion chamber (a portion of the duct with larger radius) with or without an obstacle and with or without inlets a number of numerical results with variation in frequency are presented. Numerical results for the various building-block elements of the expansion chamber are also considered.     

气动光学效应湍流脉动模型系数修正

以光学窗口的气动光学效应研究为背景,研究影响光学传输预测的流场湍流脉动量预测.选取平板、压缩拐角、凹腔流动模型,采用大涡模拟（LES）方法研究超声速湍流脉动工程预测模型的系数修正.结果表明：LES方法能够获得无激波干扰、强激波干扰及底部大分离条件下,湍流密度脉动的定量分布,并据此给出脉动工程模型的系数修正,结果已经应用于型号飞行器的光学窗口气动光学效应预测.     

面向高超声速飞行器的激波智能预测方法

高超声速飞行器激波位置的准确预测能够有效提升数值模拟的精度和效率。一方面,对高超声速飞行器激波附近网格进行正交和加密处理,可有效提升数值计算精度;另一方面,使用高超声速飞行器激波位置对计算网格进行修正,能够加速CFD计算收敛过程。提出了一种基于机器学习的高超声速飞行器激波智能预测方法,对典型高超声速飞行器外形进行激波位置的高效准确预测。首先,针对典型高超声速飞行器外形和典型飞行状态,使用数值模拟方法获得收敛的流场,并采用基于Mach数等值线的激波提取方法,从流场中判别激波面并提取构成激波面的关键点位置,形成训练数据;然后采用有监督学习算法,学习关键点位置,并利用二次曲线沿流向拟合关键点形成初步的激波线族;最后,基于剖面压力云图,构造基于投影压力图像的智能预测神经网络,对初步形成的激波线族进行修正,并获得三维激波面。大量的实验结果表明,激波预测模型能够对高超声速飞行器激波位置做出准确预测,预测的激波面与CFD数值计算结果中提取的激波面误差在10^-4量级。     

导弹超声速尾退分离干扰特性试验

超/高超声速尾退分离在防热、保形、隐身、多次投放、回收等方面具有明显优势,有望成为高超声速飞行器载荷投放的优选方案。由此面临一类新的多体分离问题：超/高超声速尾退分离问题(aft super/hypersonic ejection separation, ASES)。超/高超声速尾退分离问题本质上是带空腔底部流动与多体分离构成的耦合问题,具有流场结构复杂、气动非定常非线性非对称效应显著的特点。针对超声速尾退分离问题,采用网格测力和轨迹捕获(captive trajectory system, CTS)风洞试验方法探索了尾退分离干扰流场的结构,发现可根据流场结构和舵效变化分为低速-亚声速无激波、高亚声速-跨声速弱激波、超声速激波和准自由流弱干扰4种典型干扰特征,揭示了尾流场影响后不同区域的全弹气动特性和舵效特性以及控制律、攻角、高度和Mach数对分离位移和姿态的影响规律。相关结论将有助于增强对尾退分离问题的认识,对尾退分离技术的工程实践具有参考价值。     

电弧射流激励器工作特性仿真研究

为了在高超声速飞行器减阻中达到更好的减阻效果,设计了一种电弧射流等离子体激励器。采用有限元法求解非线性多物理方程,对此电弧射流等离子体激励器的工作特性进行了数值模拟,得到了激励器内部的电势、压力、温度和速度分布,综合分析了进气口气体速度、放电电流、激励器管道半径对电势、压力、温度和速度分布的影响。获得了全面的影响规律,通过仿真结果还得到:电弧射流等离子体激励器可产生最高温度为8638 K、最高速度为655 m/s的等离子体射流。当电流20 A,进气速度0.5 m/s,管道半径2.5 mm时,所需功率最小;当电流20 A,入口气体流速5 m/s,管道半径2.5 mm时,出口处平均温度最高;当电流20 A,进口气体速度10 m/s,管道半径2.5 mm时,出口处平均速度最大。并对仿真得到的放电电压进行了实验验证,在等离子体参数相似的情况下,实验结果与仿真结果吻合较好。     

进入轨道偏差对火星科学实验室气动力特性的影响

针对火星科学实验室(MSL)高超声速进入过程,利用三维并行程序求解流体力学Navier-Stokes方程,耦合真实气体模型,分析火星大气中真实气体效应对进入器气动力特性的影响量在进入轨道发生偏差时的变化规律.结果表明:对海盗号的计算结果与飞行数据符合很好,验证了火星大气真实气体模型和计算方法;真实气体效应影响下,激波层厚度大为减小,温度下降明显,进入器阻力系数明显增加,升力系数变化不大,俯仰力矩系数增加,基准状态下配平攻角较完全气体减小约2.2°;高度不变,Ma数增加导致阻力系数和俯仰力矩系数增大,配平攻角和完全气体的差值由1.6°增加到2.6°,表明Ma数变大时真实气体效应引起的气动力变化增强;Ma数不变,高度增加略微减弱波后化学反应,对进入器气动力特性基本没有影响.