火星探测器进入飞行气动测量方法研究

在火星探测任务的进入飞行过程中,开展气动测量将获取大量宝贵的气动特性数据,一方面可以验证进入探测器气动外形,掌握其准确的气动特性,另一方面,可以积累详实的火星大气数据,为后续任务提高落点精度等飞行性能提供支持. 本文针对火星探测器进入飞行任务,基于弹道重建和嵌入式大气数据传感系统(flushair data system, FADS),提出火星探测器进入飞行过程中的气动测量方法. 通过融合火星进入外测弹道信息,利用输出误差法实现攻角、侧滑角的高精度测量;基于嵌入式大气数据传感系统,利用最小二乘最优估计算法,建立了进入飞行动压的测量方法. 仿真分析表明,气动测量方案精度高,进入动压测量精度优于1%,攻角和侧滑角测量精度相较于当前火星探测器进入气动测量而言,至少提升1 个量级. 研究结果将为火星探测等深空探测任务的进入飞行气动测量提供技术参考.      

高超声速火星进入环境中颗粒运动特性研究

火星大气中会发生不同规模的沙尘暴,大气中蕴含的尘埃颗粒会对高速进入的火星探测器表面造成侵蚀并导致壁面热流增加,给探测器的热防护系统设计带来巨大挑战.文章针对高超声速火星进入环境两相流动问题,基于Euler-Lagrange框架建立非平衡流场与颗粒的单向耦合计算方法,采用模态半径为0.35μm的火星大气颗粒分布模型,研究不同尺寸颗粒在流场中的运动轨迹,获得高温相变模型对颗粒运动的影响以及不同粒径颗粒的撞击能量分布.结果表明,颗粒在高温流场中运动会吸热融化甚至蒸发,高温相变模型导致的颗粒直径减小对小尺寸颗粒运动轨迹有较大影响;当前计算状态下,直径3μm以上的颗粒具有较大的Stokes数且颗粒半径在运动过程中基本保持不变,其运动轨迹受流场影响较小,该尺寸颗粒的撞击分数均达95%以上,是造成壁面撞击的主要颗粒尺寸;撞击能量分数结果表明,直径3～10μm之间的颗粒是撞击能量的主要来源,约占总撞击能量的80%.     

真实气体效应对返回舱气动力特性的影响

采用7 组元6 反应化学动力学模型,通过数值方法研究了真实气体效应对阿波罗(Apollo) 返回舱流场及气动力特性的影响. 并利用典型弹道点的飞行和实验数据验证了化学非平衡流计算程序的可靠性. 计算结果表明:真实气体效应主要发生在物面附近很薄的激波层内,真实气体效应使得激波脱体距离减小;真实气体效应使阻力系数和升力系数增加,且在小攻角时增加幅度最大;真实气体效应产生附加的低头力矩,使压心位置后移. 真实气体效应的影响随着马赫数的增加变化不明显.     

基于多层分块算法的激波干扰流场预测

激波-激波干扰流场预测是超声速乃至高超声速流动中最具挑战性的问题之一. 特别地, 第IV类激波干扰由于其在壁面驻点附近产生极高的热载荷而备受关注. 本文针对圆柱诱导的弓形激波和入射斜激波的干扰问题, 分别基于量热完全气体模型和考虑振动激发的热完全气体模型, 数值求解有黏二维可压缩NS方程, 分析了高温气体效应对激波干扰流场结构, 以及第IV类激波干扰流场状态参数的影响. 接着, 本文基于一种具有广义可分离特性的遗传算法 (多层分块算法), 给出能够预测不同气体模型下第IV类激波干扰流场三波点的坐标位置、超声速射流的几何形状等特征性几何结构的数学模型, 进一步获得高温气体效应对激波干扰类型转变准则影响的定量化评估. 激波干扰类型转变准则面上的多组临界工况的激波干扰流场结构以及壁面压力和壁面热流分布的对比结果表明, 不同气体模型下的激波干扰类型和流场结构差异较为显著, 获得的定量化预测模型对工程中气动热环境的预测具有一定的参考价值.      

泰氟龙烧蚀对黏性激波层电子密度的影响

采用黏性激波层基本方程组对有无泰氟龙烧蚀两种情况下的钝锥体化学非平衡绕流作了数值求解,以研究泰氟龙烧蚀对流场电子密度的影响规律.算例结果表明:泰氟龙烧蚀确有降低流场电子密度的效应.平衡催化壁工况下这种效应的强度,显著大于非催化壁工况下的强度;远下游截面处的这种效应的强度,显著大于驻点区的强度.此外,通过计算结果分析,对形成上述影响规律的原因作了初步讨论.     

高速双锥绕流中热化学与输运模型影响研究

激波与边界层之间相互作用是高超声速飞行中的常见现象,对飞行器气动性能与飞行安全至关重要.对于高焓来流,流场中通常存在复杂的物理化学现象,此时准确模拟流场中激波边界层相互作用的难度大,相关物理化学建模仍有待进一步考察和研究.本文针对最近文献中纯净空气高超声速双锥绕流实验开展数值研究,分别研究了不同热化学模型与输运模型对壁面压力与热流的影响.热力学模型包括完全气体、热力学平衡和非平衡模型,化学模型包括冻结和非平衡化学模型,输运模型包括经典的Wilke/Blottner/Eucken模型与更加复杂的Gupta/SCEBD模型,以及考虑壁面催化/非催化影响的模型.计算了6个不同算例,涵盖了低焓至高焓来流等不同工况.壁面压力与热流的数值计算结果与实验结果符合较好;对于低焓来流,计算结果主要受到分子内能分布的影响,输运模型对计算结果的影响不大;对于高焓来流,一方面计算结果受到化学反应与壁面催化的影响较大,另一方面不同输运模型对计算结果的影响也更加明显.      

高速双锥绕流中热化学与输运模型影响研究

激波与边界层之间相互作用是高超声速飞行中的常见现象,对飞行器气动性能与飞行安全至关重要.对于高焓来流,流场中通常存在复杂的物理化学现象,此时准确模拟流场中激波边界层相互作用的难度大,相关物理化学建模仍有待进一步考察和研究.本文针对最近文献中纯净空气高超声速双锥绕流实验开展数值研究,分别研究了不同热化学模型与输运模型对壁面压力与热流的影响.热力学模型包括完全气体、热力学平衡和非平衡模型,化学模型包括冻结和非平衡化学模型,输运模型包括经典的Wilke/Blottner/Eucken模型与更加复杂的Gupta/SCEBD模型,以及考虑壁面催化/非催化影响的模型.计算了6个不同算例,涵盖了低焓至高焓来流等不同工况.壁面压力与热流的数值计算结果与实验结果符合较好;对于低焓来流,计算结果主要受到分子内能分布的影响,输运模型对计算结果的影响不大;对于高焓来流,一方面计算结果受到化学反应与壁面催化的影响较大,另一方面不同输运模型对计算结果的影响也更加明显.     

高超声速非平衡流动-辐射特性数值模拟研究

飞行器高超声速飞行过程中所承受对流加热和辐射加热可能具有相当的量级,因此合理准确预测气动加热需要将二者进行综合考虑.文章发展了具有非玻尔兹曼电子能级分布和振动能级分布的高温空气碰撞辐射模型,并耦合一维激波后流动方程计算不同飞行条件下激波后的非平衡流动特性,采用逐线辐射输运模型计算获得激波后非平衡辐射特性、辐射强度和辐射输运通量,深入比较分析了不同飞行高度和马赫数对非平衡流动和辐射输运过程的影响.计算结果表明对于高空高马赫飞行条件,其波后流动存在显著的热力学非平衡、化学非平衡和能级非平衡特征,在近激波区域高振动能级和原子高束缚电子激发态明显低于玻尔兹曼分布.在高空高马赫条件下真空紫外辐射占据主导地位,主要是由高能原子束缚-束缚跃迁造成的.随着高度和马赫数的下降,激波层内气体解离和电离程度降低,原子辐射贡献下降,分子辐射贡献增加,导致红外、可见光和紫外波段的辐射输运增强,真空紫外辐射输运过程减弱.     

高温激波管化学非平衡流动数值模拟研究

为了预测氢氧定容燃烧驱动的高温激波管性能,需要准确分析激波管非定常化学非平衡流动过程.本文在破膜前的驱动段定容燃烧以及破膜后的化学非平衡流动数值模拟中,引入双时间步长方法,发展高温激波管化学非平衡流动数值模拟方法,该方法在时间上具有二阶精度.计算结果与目前存在的激波管流动解析解以及零维化学反应系统的数值解进行了比较,吻合较好.对于典型高温激波管状态,采用有限体积方法离散准一维流动Euler控制方程,并通过将流动过程和化学反应动力学过程耦合求解,获得了激波管内部的化学非平衡流动特征.     

探测器对超音速刚性盘-缝-带型降落伞系统的影响

研究了流体初始马赫数为 2.0 时, 探测器的存在与否对刚性盘-缝-带型降落伞系统气动减速性能以及流场流体结构特性的影响. 对于非定常可压缩流体的数值模拟, 流场采用了三层块结构自适应网格加密技术, 配合混合形式的TCD (tuned center difference)和WENO (weighted essentially non-oscillatory)计算格式以及基于拉伸涡亚格子模型的大涡模拟方法来处理超音速流中的激波以及大尺度湍流旋涡结构等. 结果表明: 无探测器时, 降落伞系统的流场结构稳定, 扰动较小; 有探测器存在时, 探测器后端的湍流尾迹和伞衣内部逆向运动溢出的流体与伞衣前端的弓形激波周期性的相互作用, 使得激波位置发生前移、激波倾角变小, 伞衣内部流场难以达到平衡稳定状态. 这加剧了降落伞系统的气动阻力振荡脉动变化, 降低了降落伞系统气动阻力系数, 同时也使得降落伞系统流场尾迹结构更加复杂.      

Comparison of different models for non-equilibrium CO2 flows in a shock layer near a blunt body

The paper presents results of a numerical simulation of a supersonic two-dimensional (2D) viscous flow containing CO₂ molecules near a spacecraft entering the Mars atmosphere. The gas–dynamic equations in the shock layer are coupled to the equations of non-equilibrium vibrational and chemical kinetics in the five-component mixture CO₂/CO/O₂/C/O. Transport and relaxation processes in the flow are studied on the basis of the rigorous kinetic theory methods; the developed transport algorithms are incorporated in the numerical scheme. The influence of the vibrational excitation of CO₂ and chemical reactions on the gas flow parameters and heat transfer is analyzed. The obtained results are compared with those found using two simplified models based on the two-temperature and one-temperature vibrational distributions in CO₂. The accuracy of the simplified models and the limits of their validity within the shock layer are evaluated. The effect of bulk viscosity in a flow near a re-entry body is discussed. The role of different diffusion processes, chemical reactions, and surface catalytic properties in a flow of the considered mixture in the shock layer is estimated.     

The effect of an unsteady drag force on the structure of a non-equilibrium region behind a shock wave in a gas-particle mixture

For numerical analysis of shock wave propagation in gas-particle mixtures, drag coefficients of a sphere in steady flows are generally used. However, it is shown both experimentally and numerically that a shock loaded solid sphere experiences unsteady drag forces. The paper describes a model of unsteady drag force and its effect on the structure of the non-equilibrium region behind a shock front traveling in a dusty gas. The results are compared with those obtained by using a steady drag coefficient and are discussed. It is demonstrated that the large drag force at the early stage of the interaction between shock-wave induced flow and a solid particle affects the flow structure that is obtained with a steady drag force.      

On the interference of a weak and a strong shock wave in a dissociating nitrogen flow

The influence of the nitrogen dissociation on the interactions due to the interference of two planar shock waves in a hypersonic high enthalpy flow is theoretically investigated for infinite reaction rates. The two limiting cases of infinitely slow and infinitely fast reactions are modelled as a perfect gas and an ideal dissociating gas in chemical equilibrium.To investigate the influence of finite reaction rates on the interactions of shock waves, experiments are performed in the high enthalpy shock tunnel Göttingen (HEG) with a wind tunnel model consisting of a wedge type shock generator and a transversally mounted cylinder. The pressure and heat transfer loads resulting from the shock wave interferences are measured and the flow field is visualized by means of interferograms. The experimental results are compared with the results of a numerical simulation for a dissociating nitrogen flow and with the experimental results for a perfect gas flow.     

Calculation of stagnation streamline quantities in hypersonic blunt body flows

An approximate method for the efficient calculation of stagnation-streamline quantities in hypersonic flows about spheres or cylinders is suggested. Based on the local similarity of the flow field the two-dimensional Navier-Stokes equations are simplified to a one-dimensional approximation for the stagnation streamline. These equations are solved with an implicit finite-volume scheme. Comparisons with fully two–dimensional Euler and Navier–Stokes calculations for flows about spheres are presented, that include perfect gas flows and flows in chemical non-equilibrium. Comparisons with a number of experiments conclude this report. Received 8 May 1996 / Accepted 31 October 1996     

Study on the numerical schemes for hypersonic flow simulation

Hypersonic flow is full of complex physical and chemical processes, hence its investigation needs careful analysis of existing schemes and choosing a suitable scheme or designing a brand new scheme. The present study deals with two numerical schemes Harten, Lax, and van Leer with Contact (HLLC) and advection upstream splitting method (AUSM) to effectively simulate hypersonic flow fields, and accurately predict shock waves with minimal diffusion. In present computations, hypersonic flows have been modeled as a system of hyperbolic equations with one additional equation for non-equilibrium energy and relaxing source terms. Real gas effects, which appear typically in hypersonic flows, have been simulated through energy relaxation method. HLLC and AUSM methods are modified to incorporate the conservation laws for non-equilibrium energy. Numerical implementation have shown that non-equilibrium energy convect with mass, and hence has no bearing on the basic numerical scheme. The numerical simulation carried out shows good comparison with experimental data available in literature. Both numerical schemes have shown identical results at equilibrium. Present study has demonstrated that real gas effects in hypersonic flows can be modeled through energy relaxation method along with either AUSM or HLLC numerical scheme.     

Phase step holographic interferometry applied to hypervelocity,non-equilibrium cylinder flow

For blunt bodies the reduced bow shock wave stand-off distance and the shock layer density rise in the stagnation region as compared to ideal gas flow are phenomena caused by dissociative effects. In this work experiments with aR _BODY=45 mm radius cylinder and an aspect ratio ofL/R _BODY 11 are described. The tests were carried out in the High Enthalpy Shock Tunnel in Göttingen (HEG), a free piston driven shock tunnel. Two different test conditions at reservoir enthalpies of around 21 MJ/kg and Mach numbers of around 9 in air and nitrogen were available.Optical measurements with a holographic phase step interferometer to obtain complete flow field density gradients have been carried out. By increasing the signal-to-noise ratio with the technique of phase stepping over the original recording quality, high quality interferograms are obtained. The high spatial resolution of the holograms results in the creation of hologram-schlieren images which are compared to directly recorded laser-schlieren images.Infinite and finite fringe interferograms and the complete density fields for the two free-stream conditions are presented. The stagnation line densities are quantified. The measured results are shown and compared with Navier-Stokes calculations which account for chemical reactions in the flow. The numerical code underpredicts the stand-off distance of the bow shock wave. It is shown that the flow behind the bow shock wave is in non-equilibrium and that it reaches equilibrium before the body for one condition.This work was an offshoot of an ESA research contract monitored by Dr. D. Vennemann, whose support is gratefully acknowledged. This study of stagnating high enthalpy flows has exceedingly benefited from discussions with S. Brück and V. Hannemann. The operation of the large wind tunnel HEG is the result of a team effort. The with to thank the whole team, represented by Dr. W.H. Beck, for keeping the tunnel going.     

The flow gradients in the vicinity of a shock wave for a thermodynamically imperfect gas

Supersonic rotational planar and axisymmetric flows of a non-viscous, non-heat-conductive gas with arbitrary thermodynamic properties in the vicinity of a steady shock wave are studied. The differential equations describing the gas flow upstream and downstream of the discontinuity surface and the dynamic compatibility conditions at this discontinuity are used. The gas flow non-uniformity in the shock vicinity is described by the spatial derivatives of the gasdynamic parameters at a point on the shock surface. The parameters are the gas pressure, density, and velocity vector. The derivatives with respect to the directions of the streamline and normal to it, and of the shock surface and normal to it, are considered. Spatial derivatives of all gasdynamic parameters are expressed through the flow non-isobaric factor along the streamline, the streamline curvature, and the flow vorticity and non-isoenthalpy factors. An algorithm for determining these factors of the gas flow downstream of a shock wave is developed. Example calculations of these factors for imperfect oxygen and thermodynamically perfect gas are presented. The influence coefficients of the upstream flow factors on the downstream flow factors are calculated. The gas flow in the vicinity of the shock is described by the isolines of gasdynamic parameters. Uniform planar and axisymmetric flows at different distances from the axis of symmetry are examined; the isobars, isopycnics, isotachs and isoclines are used to characterize the downstream flow behind a curved shock in an imperfect gas.     

Decay of a plane shock wave in a two-parameter medium with arbitrary equation of state

Special curves, called shock polars, are frequently used to determine the state of the gas behind an oblique shock wave from known parameters of the oncoming flow. For a perfect gas, these curves have been constructed and investigated in detail [1]. However, for the solution of problems associated with gas flow at high velocities and high temperatures it is necessary to use models of gases with complicated equations of state. It is therefore of interest to study the properties of oblique shocks in such media. In the present paper, a study is made of the form of the shock polars for two-parameter media with arbitrary equation of state, these satisfying the conditions of Cemplen's theorem. Some properties of oblique shocks in such media that are new compared with a perfect gas are established. On the basis of the obtained results, the existence of triple configurations in steady supersonic flows obtained by the decay of plane shock waves is considered. It is shown that D'yakov-unstable discontinuities decompose into an oblique shock and a centered rarefaction wave, while spontaneously radiating discontinuities decompose into two shocks or into a shock and a rarefaction wave.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 147–153, November–December, 1982.