MULTI-OBJECTIVE OPTIMIZATION AND AERODYNAMIC PERFORMANCE ANALYSIS OF THE UPPER SURFACE FOR HYPERSONIC VEHICLES

为分析小攻角巡航条件下吸气式高超声速飞行器上壁面的变化对其气动性能和容积的影响, 以参数化后的飞行器上壁面对称面型线为设计变量, 在飞行马赫数6.5, 飞行高度27 km, 飞行攻角为4°的条件下, 采用计算流体力学为性能分析工具, Pareto多目标遗传算法为优化设计方法, 开展了二维条件下的升阻比/容积双目标优化设计. 在此基础上, 选择典型的二维优化结果, 重构生成对应的三维构型并进行数值分析, 获得了飞行器气动性能和容积间的相互关系. 结果表明在巡航条件下, 尽管二维/三维条件下飞行器的气动参数数值有较大差别, 但在这2种条件下, 飞行器的升阻比和容积间的关系均近似呈线性反比例关系. 同时, 对于三维构型而言, 在给定容积不变的条件下, 通过改变上壁面对称面型线的形状仅能使升阻比获得较小的增量(约0.36%). 相比之下, 当给定升阻比基本不变的条件下, 飞行器容积可调空间相对较大, 约为1.93%. 此外, 计算结果还表明, 在飞行器的容积基本不变情况下, 通过调节上壁面对称面型线, 可使飞行器的俯仰力矩获得5%左右的调节空间, 且其升阻比基本不变.     

About the formulation,verification and validation of the hypersonic flow solver Eilmer

We describe the formulation of the gas dynamics and high‐temperature thermochemical modules of the Eilmer code, an open‐source Navier–Stokes solver for transient compressible flow in two and three dimensions. The core gas dynamics formulation is based on finite‐volume cells, and the thermochemical effects are handled with specialised updating schemes that are coupled into the overall time‐stepping scheme. Verification of the code is explored via a number of case studies that use analytic and semi‐analytic solutions as comparison. These include both smooth and shocked flows and are used to demonstrate the order of spatial accuracy of the code. Cases include manufactured solutions for rather abstract inviscid and viscous flow, an idealised detonation wave supported by a curved body, and the transient flow of an idealised but high‐performance shock tube. Validation of the inviscid gas dynamics and thermochemical models is then explored using data from a selection of experimental studies. These studies include ballistic range experiments with chemically‐inert noble gases and high‐temperature chemically‐reacting air. These comparisons show that the code performs well and they provide a lesson in considering a range of experimental data rather than relying upon isolated data points for validation. These verification and validation cases are described in full detail and will be useful for other code developers of high‐temperature compressible flow solvers. Copyright © 2013 John Wiley & Sons, Ltd.     

PIV技术在超及高超声速流场测量中的研究进展

本文分析了超声速流场对测量技术的特殊要求, 归纳了目前将粒子影像测速仪(particle image velocimetry, PIV)技术应用于超声速流场的测量时所面临的主要技术难点以及主要的解决方法, 分析了超声速流场中所用PIV粒子的主要要求、粒子特性、投放方法等, 介绍了PIV技术在超声速、高超声速流场测量中最新的国内外进展, 特别是给出了国内外关于高超声速流场中激波/附面层的相互干扰, 以及高超声速飞行器超燃冲压发动机主要部件内流场的PIV试验研究的最新进展.      

钝头体高超声速绕流底部失稳特征数值模拟

利用数值模拟方法对高超声速钝锥及Apollo返回舱底部尾迹流场进行了研究, 分析尾迹流动的失稳过程. 对钝锥模型, 在M_∞=6, Re=1.71× 10⁶(Re以球头半径为参考长度)条件下观察到了底部流动的不稳定性. 不添加任何扰动, 数值模拟首先得到的流动是稳定解, 在底部发展出一个主分离区和一个二次分离区, 流动是轴对称状态. 继续进行计算, 发现二次分离线率先变形, 底部流场发展出非定常周期流动. 对Apollo返回舱模型, 在相同条件下 (Re以前面圆弧半径为参考长度), 数值模拟首先得到的流动同样是稳定解, 出现以二次分离线率先变形为起始的结构失稳, 演化出周期性过程, 但持续时间较短, 很快出现了非周期非对称状态. 研究表明, 高超声速钝锥及Apollo返回舱底部流场均存在不稳定性问题, Apollo返回舱的底部流场更加不稳定.     

制度落差与中国跨国企业海外进入模式关系研究

本文分析了制度落差与中国企业海外市场进入模式之间的关系,着重对比了制度顺差与制度逆差情况下,制度差异对企业海外进入模式决策的影响。研究发现,在正式制度顺差的情况下,东道国与母国间正式制度差异越大,企业越倾向于选择以合资形式建立子公司;在正式制度逆差情况下,东道国与母国间正式制度差异越大,企业越倾向于选择以全资形式建立子公司。同时发现,正式制度落差并不影响非正式制度对进入模式的影响方向,即随着母国与东道国非正式制度差异的增大,企业更倾向于以全资形式建立子公司。     

基于粒子群的高超声速飞行器闭环制导研究

对于高超声速飞行器的研究来说,凭借其超高的飞行速度和飞行高度,能够在短时间内飞行更远的距离。为了达到作战要求,需要规划整个飞行轨迹,而最为重要的便是上升段的制导问题。本文以X-33高超声速飞行器模型为研究对象,提出基于粒子群算法的闭环制导策略,实时修正飞行轨迹,使飞行器最终准确到达目标位置;并对该方法的可靠性进行了仿真验证。仿真结果表明,基于粒子群算法的闭环制导策略优化精度高,物理概念明确,能满足高超声速飞行器上升段的闭环制导需求。     

Flow characteristics of hypersonic inlets with different cowl-lip blunting methods

Under hypersonic flight conditions,the sharp cowl-lip leading edges have to be blunted because of the severe aerodynamic heating.This paper proposes four cowl-lip blunting methods and studies the corresponding flow characteristics and performances of the generic hypersonic inlets by numerical simulation under the design conditions of a flight Mach number of 6 and an altitude of 26 km.The results show that the local shock interference patterns in the vicinity of the blunted cowl-lips have a substantial influence on the flow characteristics of the hypersonic inlets even though the blunting radius is very small,which contribute to a pronounced degradation of the inlet performance.The Equal Length blunting Manner(ELM)is the most optimal in that a nearly even reflection of the ramp shock produces an approximately straight and weak cowl reflection shock.The minimal total pressure loss,the lowest cowl drag,maximum mass-capture and the minimal aeroheating are achieved for the hypersonic inlet.For the other blunting manners,the ramp shock cannot reflect evenly and produces more curved cowl reflection shock.The Type V shock interference pattern occurs for the Cross Section Cutting blunting Manner(CSCM)and the strongest cowl reflection shock gives rise to the largest flow loss and drag.The cowl-lip blunted by the other two blunting manners is subjected to the shock interference pattern that transits with an increase in the blunting radius.Accordingly,the peak heat flux does not fall monotonously with the blunting radius increasing.Moreover,the cowl-lip surface suffers from severe aerothermal load when the shear layer or the supersonic jet impinges on the wall.     

Evaluation of Euler fluxes by a high-order CFD scheme: shock instability

The construction of Euler fluxes is an important step in shock-capturing/upwind schemes. It is well known that unsuitable fluxes are responsible for many shock anomalies, such as the carbuncle phenomenon. Three kinds of flux vector splittings (FVSs) as well as three kinds of flux difference splittings (FDSs) are evaluated for the shock instability by a fifth-order weighted compact nonlinear scheme. The three FVSs are Steger–Warming splitting, van Leer splitting and kinetic flux vector splitting (KFVS). The three FDSs are Roe's splitting, advection upstream splitting method (AUSM) type splitting and Harten–Lax–van Leer (HLL) type splitting. Numerical results indicate that FVSs and high dissipative FDSs undergo a relative lower risk on the shock instability than that of low dissipative FDSs. However, none of the fluxes evaluated in the present study can entirely avoid the shock instability. Generally, the shock instability may be caused by any of the following factors: low dissipation, high Mach number, unsuitable grid distribution, large grid aspect ratio, and the relative shock-internal flow state (or position) between upstream and downstream shock waves. It comes out that the most important factor is the relative shock-internal state. If the shock-internal state is closer to the downstream state, the computation is at higher susceptibility to the shock instability. Wall-normal grid distribution has a greater influence on the shock instability than wall-azimuthal grid distribution because wall-normal grids directly impact on the shock-internal position. High shock intensity poses a high risk on the shock instability, but its influence is not as much as the shock-internal state. Large grid aspect ratio is also a source of the shock instability. Some results of a second-order scheme and a first-order scheme are also given. The comparison between the high-order scheme and the two low-order schemes indicates that high-order schemes are at a higher risk of the shock instability. Adding an entropy fix is very helpful in suppressing the shock instability for the two low-order schemes. When the high-order scheme is used, the entropy fix still works well for Roe's flux, but its effect on the Steger–Warming flux is trivial and not much clear.     

NORMAL GRID RULE FOR HYPERSONIC HEAT FLUX NUMERICAL SIMULATION

采用数值求解雷诺平均Navier-Stokes 方程方法对高超声速钝双锥模型进行了计算,比较了不同法向网格尺度对热流计算的影响,并在此基础上提出了一种适于热流计算的网格准则. 研究表明:在实际外形的热流计算中,使用过密或过粗的法向网格均不利于得到准确的热流数据;采用该文提出的基于当地温度梯度的法向网格尺度准则可以在保证热流计算精度的同时实现较高的收敛速度.     

迎风凹腔与逆向喷流组合热防护系统冷却效果研究

对迎风凹腔与逆向喷流组合热防护系统的冷却效果进行了分析, 研究了相同总压不同流速的逆向喷流对组合结构的流场、气动受力及壁面传热的影响. 通过与相关的实验结果对比, 验证了数值方法的可靠性. 研究发现:该结构能够有效地对飞行器鼻锥表面进行冷却, 引入很小总压的逆向喷流(逆喷总压比 PR=0.1), 组合结构的冷却效果就可以远远优于单一的迎风凹腔; 相同逆向喷流总压下, 逆喷速度越高, 逆喷流量越大, 外壁面的冷却效果越好; 随逆喷流速提高, 气动阻力也进一步减小. 本文研究的组合结构非常适用于