“Ⅳ型”激波干扰中流-热-固耦合问题一体化计算分析

针对高超声速进气道前缘"Ⅳ型"激波干扰产生的气动加热与结构传热多物理场耦合计算问题，发展了一种基于有限体积法的流-热-固一体化计算方法。该方法采用一体化控制方程组统一离散求解外部高速流场与内部结构温度场，规避了传统分区耦合算法在时间域内交替迭代的繁琐数据交换策略。另外，提出一种新的双温阻模型计算流-固交界面的物性参数以保证计算准确性，采用LU-SGS隐式时间迭代和自适应时间步长以提高计算效率。采用经典高超声速二维圆管流-热-固耦合算例对该一体化方法进行验证，计算结果与试验值和参考文献数据吻合较好，证明了该方法的可靠性和正确性。利用一体化方法对高超声速前缘"Ⅳ型"激波干扰流-热-固耦合问题进行定常/非定常计算与分析，给出了温度与热流的时变特性，计算结果表明，激波干扰作用产生的超声速"喷流"不断冲击壁面，使得壁面最大压力系数增大约9倍，壁面最大热流增大约4.7倍，给高速飞行器的热防护设计与选材带来严峻挑战。同时，也表明了一体化计算方法可以较好地用于长航时飞行条件下与复杂飞行环境下的高超声速热防护系统的热环境特性分析与综合性能评估。     

高超声速飞行器流动特征分析

在非流线型构件或突起物的扰动效应、高马赫数和低雷诺数极限效应、低湍流度环境效应和由激波或摩擦导致的气动加热效应等4个方面的影响下,未来高超声速飞行器涉及的流动主要表现出这样的特点:典型流动结构强度高、尺度大,如强激波和厚边界层;局部流动结构数量多;激波、膨胀波和边界层结构之间相互干扰十分严重;转捩、压力脉动和一些流动结构对细微因素非常敏感;压力、摩擦应力和热流峰值现象普遍;升阻比屏障难以突破;流场同时依赖大量无量纲参数和有量纲参数,导致实验模拟难度大。本文在回顾传统高超声速流动主要流动现象的基础上,对上述7个方面涉及的典型流动现象的基础研究现状、问题本质和因果关系进行综合描述,讨论如何更有效地面对基础研究和工程实际问题。该文既可为解决典型流动现象中尚未解决的基础研究提供帮助,也可为如何合理地利用有限的已知知识解决工程应用问题提供指导。     

Numerical simulations of shock/boundary layer interactions using time-dependent modeling techniques: A survey of recent results

Supersonic or hypersonic flows within and around flight vehicles inevitably involve interactions of strong shock waves with boundary layers. Flows within inlet/isolator configurations, and flows induced by control surface deflections are some examples. Such interactions are time dependent in nature and are often subject to low-frequency, large-scale motion that induces local pressure and heating loads. With recent increases in available computer power, it has now become possible to simulate such interactions at experimentally relevant Reynolds numbers using time-dependent techniques, such as direct numerical simulation (DNS), large-eddy simulation (LES), and hybrid large-eddy simulation/Reynolds-averaged Navier–Stokes (LES–RANS) methods. This paper will survey some recent work in this area and will describe insights in shock/boundary layer interaction physics gained by using these high-fidelity methods. Attention will be focused on studies that compare directly with experimental data at the same (or nearly the same) Reynolds number. Challenges in the application of these techniques to even more complicated high-speed flow fields are also outlined.     

面对称高超飞行器几何参数与总体性能相关性研究

气动技术是高超声速飞行器的重要支撑技术,气动布局的选择与确定是影响飞行器总体方案论证的重要因素.以一类面对称滑翔飞行器为研究对象,采用针对高超声速飞行器的快速工程算法,建立参数化外形的气动数据库,实现由气动布局要素到总体性能指标的映射,并通过数据分析方法,如相关性分析、灵敏度分析等,梳理出影响高超声速滑翔飞行器不同总体性能的关键气动外形几何参数,并进行相关性排序,研究结果可指导高超声速滑翔飞行器的布局设计及优化.     

Investigation of thermal protection system by forward-facing cavity and opposing jet combinatorial configuration

 This paper focuses on the usage of the forward-facing cavity and opposing jet combinatorial configuration as the thermal protection system (TPS) for hypersonic vehicles. A hemispherecone nose-tip with the combinatorial configuration is investigated numerically in hypersonic free stream. Some numerical results are validated by experiments. The flow field parameters, aerodynamic force and surface heat flux distribution are obtained. The influence of the opposing jet stagnation pressure on cooling efficiency of the combinatorial TPS is discussed. The detailed numerical results show that the aerodynamic heating is reduced remarkably by the combinatorial system. The recirculation region plays a pivotal role for the reduction of heat flux. The larger the stagnation pressure of opposing jet is, the more the heating reduction is. This kind of combinatorial system is suitable to be the TPS for the high-speed vehicles which need long-range and long time flight.     

吸气式高超声速飞行器多学科动力学建模

高超声速飞行器一体化设计中存在气动/热/推进/结构弹性相互耦合的问题,首先根据飞行器的机体/发动机一体化设计思想构造了二维高超声速飞行器模型,并基于激波/膨胀波原理和动量定理建立了气动力模型,采用Chavez和Schmidt建立的超燃冲压发动机推进系统模型;在飞行器结构方面,引入变截面和变质量分布的自由梁结构模型,并采用Eckert参考焓方法分析的气动加热过程中承力梁不同轴向位置温度随时间变化特征,在此基础上运用模态法计算了燃料消耗和气动加热条件下结构的固有频率和振型特征,获得结构弹性变形的模型;最后建立了考虑热气动弹性和推进系统作用的飞行动力学方程。研究结果表明:质量变化对结构弹性特性影响比较显著,而气动加热的影响主要表现在振动频率方面,且会随着加热过程的持续而逐渐增强;结构变形会改变飞行器静配平状态,特别是在机体质量较大的最初飞行阶段,气动加热会强化结构变形对配平特征的影响;线性化系统的动力学特征分析表明,质量减小和结构变形均会增加短周期模态和振荡模态的不稳定特性,而对高度特性的影响不大,气动加热效应会进一步增加飞行力学和气动弹性的耦合特征,并导致弹性模态的稳定性降低。     

高超声速推进系统用单膨胀斜面喷管型面设计和流场模拟

基于特征线法,并考虑变比热的影响,开展了高超声速飞行器用单膨胀斜面喷管设计。利用CFD数值模拟技术,计算得到了设计状态和沿飞行轨迹其它飞行状态下的单膨胀斜面喷管内外流场和特性。结果表明,在设计状态马赫数5时,基于特征线法得到的单膨胀斜面喷管内流场分布符合设计要求,而在其它较低的飞行马赫数下,单膨胀斜面喷管处于过膨胀状态,并且过膨胀的程度随飞行马赫数的降低而愈加严重。在马赫数2.5时,喷管膨胀面气流已发生明显分离,喷管性能急剧恶化。为了提高低马赫数条件下单膨胀斜面喷管的性能,采用变几何结构(调节下斜板角度)或基于二次流控制的单膨胀斜面喷管是必须的。     

高超声速飞行器宽速域翼型/机翼设计与分析

对于现代高超声速飞行器的设计而言,除了需要保证高超声速的性能外,还必须兼顾满足工程需求的亚跨超声速特性。首先,采用雷诺平均Navier-Stokes（RANS）方程流动求解器,结合基于Kriging模型的代理优化算法,开展了高超声速飞行器宽速域翼型的优化设计研究,设计出了一种下表面具有双"S"形特征的新翼型。综合性能评估结果表明,该翼型相比于常规的高超声速翼型,在跨声速和高超声速下具有更加优良的气动特性;其跨声速状态下的升阻比达到78.9,高超声速状态下的升阻比达到5.94,能够实现宽速域内良好的综合气动性能。其次,开展了仿德国"桑格尔号"（SANGER）空天飞机运载机机翼的气动特性研究,对配置宽速域翼型与常规高超声速翼型的机翼进行了气动力特性综合对比分析。结果表明,配置新翼型的机翼在宽速域范围内整体气动性能更优,说明所设计的宽速域翼型在三维机翼上也具有一定的实用价值。     

高超声速飞行器机体/推进一体化设计的启示

机体/推进一体化设计是吸气式高超声速飞行器的关键技术。飞行器的前体和后体既是主要的气动型面,又是发动机进气道的外压缩型面和尾喷管的膨胀型面,一体化设计直接影响飞行器的气动与发动机性能。本文阐述了吸气式高超声速飞行器的主要特点,梳理了飞行器的推阻匹配、升阻比特性、操稳匹配等主要气动设计问题。通过对国外典型高超声速飞行器机体/推进一体化设计技术的综合分析,总结了前体/进气道、后体/尾喷管、边界层强制转捩装置等关键部件的气动设计方法,获得了有意义的启示,可为后续吸气式高超声速技术研究提供重要参考。     

Minimizing Effects of Reflected Shocks in Ballistic Ranges

Shock waves originating from hypersonic projectiles flown in ballistic ranges and subsequently reflected from the range walls have been found to produce undesirable effects on the measurements of projectile wakes. This paper reports on the mechanisms through which reflected shock waves perturb the wake, the design, and installation of full-scale fiberglas wedge shock attenuation treatments, the design of equipment for minimum shock reflection, and the results of measurements of the performance of the treatment.     

进气道内激波/边界层干扰及控制研究进展

进气道作为高速航空航天飞行器的重要气动部件,对飞行器的气动力特性、结构重量、隐身性能等有显著影响。激波/边界层干扰现象是高速进气道内普遍存在的一类流动现象,对进气道的性能有突出的影响。发生于进气道内的激波/边界层干扰现象主要可分为正激波/边界层干扰、斜激波/边界层干扰以及三维激波/边界层干扰几类,由于受到侧壁壁面和进气道内背景波系的影响,这些干扰现象偏离了传统基于简化模型的研究结果,具有显著的耦合干扰特征,干扰区间内三维特征明显。概述了发生于进气道内的激波/边界层干扰特性及相关研究进展,并对目前进气道内激波/边界层干扰现象的控制方法进行了总结。     

涡波效应宽速域气动外形设计

拓展了密切锥乘波体设计方法的应用,推导了设计方法中激波出口型线、流线追踪起始线与平面形状轮廓线之间的几何关系,建立了定平面乘波体设计方法。通过定制乘波体的平面形状引入涡效应,提出涡波效应宽速域气动布局的概念,即在高超声速状态下使用激波效应、在低速状态下使用漩涡效应提升布局的总体性能。以双后掠布局为例,使用CFD方法评估其高速和低速状态的气动性能,与带锥体的平板进行对比,分析了升阻比、升力系数以及流场特性,初步给出了非线性增升效果。计算结果表明：当前定平面乘波体布局在低速状态和高超声速状态均具有较好的气动性能,弥补了传统乘波体的性能缺陷,为宽速域气动布局的设计提供了新的思路。     

Recent advances in the aerothermodynamics of spiked hypersonic vehicles

Among a variety of design requirements, reducing the drag and aeroheating on hypersonic vehicles is the most crucial one. Unfortunately, these two objectives are often conflicting. On one hand, sharp slender forebodies design reduces the drag and ensures longer ranges and more economic flights. However, they are more vulnerable to aerodynamic heating. On the other hand, blunt forebodies produce more drag, however, they are preferred as far as aeroheating is concerned. In addition, in the context of hypersonic vehicles, blunt geometries are preferred over slender ones for practical implications such as higher volumetric efficiency, better accommodation of crew or on-board equipment.In principle, a blunt vehicle flying at hypersonic speeds generates a strong bow shock wave ahead of its nose, which is responsible for the high drag and aeroheating levels. There have been a number of efforts devoted towards reducing both the drag and the aeroheating by modifying the flowfield ahead of the vehicle's nose. Of these techniques, using spikes is the simplest and the most reliable technique. A spike is simply a slender rod attached to the stagnation point of the vehicle's nose. The spike replaces the strong bow shock with a system of weaker shocks along with creating a zone of recirculating flow ahead of the forebody thus reducing both drag and aeroheating.Since their introduction to the high-speed vehicles domain in the late 1940s, spikes have been extensively studied using both experimental facilities and numerical simulation techniques. The present paper is devoted to surveying these studies and illustrating the contributions of the authors in this field. The paper also raises some of the areas in the field that need further investigations.     

超/高超声速飞行器动态稳定性导数极快速预测方法

飞行器设计早期阶段需要预测大量工况下的动导数。本文发展了一种面向超/高超声速飞行器的动导数极快速预测方法：首先基于当地流活塞理论,将飞行器进行小幅非定常运动所受到的气动力分为受自由来流引起的无附加扰动项以及受物面变形或运动引起的附加扰动项;通过当地表面斜度法、激波后等熵关系求解物面当地流动参数,进而结合非定常运动规律求出飞行器所受非定常气动力;再采用待定系数法对非定常气动力进行提取、辨识,最终得到超/高超声速飞行器动导数。该方法克服了传统方法对CFD流场参数的依赖和耦合,具有极高的计算效率;同时典型算例验证表明,该方法在超声速、高超声速工况下都能够很好预测动导数变化趋势。将该方法应用于复杂外形飞行器动导数预测,并讨论了与CFD方法的误差来源。本文方法可作为高速飞行器总体设计阶段布局选型的工具。     

细长体飞行器自由边界热模态试验与仿真

随着耐热承载一体化材料在新型高超声速飞行器上的应用,承力结构的工作温度不断提高,各类热模态特性逐渐得到关注。针对非平面形状的细长体飞行器自由边界条件下的热模态特性开展了研究。通过研究模拟气动加热条件的圆筒形加热笼、模拟自由边界的耐高温柔性支撑、非接触激光多普勒测振、耐高温激振杆激励等试验方法,获得了细长体结构自由边界条件下随温度变化的前3阶模态变化情况。结果表明：对此类薄壁长圆筒类结构,温度升高对模态频率影响可以超过6 Hz。开展有限元仿真,并与试验取得的热模态结果的变化规律进行对比。结果表明：建立考虑温度对结构弹性模量、热应力影响的壳单元模型,能够较好地预测出前3阶模态频率在全部受热时间范围内的最大下降量,可为高超声速飞行器控制系统设计时的拉偏范围提供参考。     

乘波体构型飞行器的高超声速测压实验研究

乘波体构型是高超声速飞行器的重要气动布局之一。对某多目标优化设计的乘波体构型飞行器进行了高超声速测压实验,对其气动性能进行风洞实验验证。实验马赫数M=6和M=7,迎角α=-4°、-2°、0°、2°、4°、6°、8°。结果表明：该乘波体构型各部件气动性能良好。进气道唇口准确捕捉到压缩激波,激波位置与设计吻合。乘波体上表面流向压力变化不大,有利于减小乘波体飞行阻力。下表面经过进气口内压段时压力有明显的增大,后体膨胀效果显著。在设计状态下,该乘波体飞行器整体气动性能良好。     

JF-12激波风洞在火星进入环境下的运行特性

火星探测是目前国际深空探测的热点,而火星着陆是火星探测的关键技术之一,火星着陆器发展中面临的一个严峻的挑战是其气动环境远远不同于地球大气的空气。然而,现阶段大多数地面试验设备都是以空气为试验气体来设计的,而不是火星大气的CO₂。本文利用高温热化学反应流动数值计算技术,对JF-12反向爆轰驱动激波风洞在火星进入环境下（主要气体成分是CO₂）的运行特性进行了计算模拟,通过调整激波管中驱动/被驱动气体的初始参数和高/低压段的截面积比,来模拟其中的波系产生、传播过程以及反射激波与接触面的相互作用机制。研究发现,相同情况下驱动CO₂的缝合激波马赫数要明显高于空气,通过减小驱动气体的声速和低压段的直径,可以在驱动CO₂时获得驻室压力稳定的试验气体。     

乘波体飞行器气动布局优化设计

由于具有高升阻比,乘波体是高超声速巡航飞行器气动布局的首选方案。文章在求解圆锥激波流场精确解的基础上,应用流线追踪方法,建立了乘波体飞行器气动布局的参数化模型。在此基础上,对飞行器的气动力特性进行了估算。最后,以气动布局参数为设计变量,升阻比最大化为设计目标,对乘波体飞行器进行气动布局优化设计,应用改进的粒子群优化算法(Particle Swarm Optimization,PSO),对优化模型进行求解,得到了优化的气动布局设计方案。     

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

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

Mechanisms of plasma actuators for hypersonic flow control

A summary of recent research progress in hypersonic plasma actuators for flow control is attempted. It is found that the most effective plasma actuator is derived from an electromagnetic perturbation and amplifies by a subsequent viscous–inviscid interaction. Computational efforts using drift-diffusion theory and a simple phenomenological plasma model, as well as experiments in a hypersonic plasma channel, have shown the effectiveness of using electro–aerodynamic interaction as a hypersonic flow control mechanism. In principle, the plasma actuator based on magneto–aerodynamic interaction should have an added mechanism in the Lorentz force, making it even more effective as a flow control mechanism. However, this approach also incurs additional challenges and complications due to the Hall effect.

Magneto–aerodynamic interactions have also been demonstrated for separated flow control, albeit in a very limited scope. Numerical simulations based on a simple phenomenological plasma model have shown the feasibility of separated flow suppression in shock-boundary-layer interaction over a compression ramp at a hypersonic flow of Mach 14.1. The control mechanism relies on the Lorentz force to energize the retarded shear layer in the viscous interacting region, but the effectiveness of momentum transfer via inelastic collision requires further validation.