Numerical validation and parametric investigation on the cold flow field of a typical cavity-based scramjet combustor

The three-dimensional coupled implicit Reynolds Averaged Navier–Stokes (RANS) equations and the two equation standard k–ε turbulence model has been employed to numerically simulate the cold flow field in a typical cavity-based scramjet combustor. The numerical results show reasonable agreement with the schlieren photograph and the pressure distribution available in the open literature. The pressure distribution after the first pressure rise is under-predicted. There are five shock waves existing in the cold flow field of the referenced combustor. The first and second pressure rises on the upper wall of the combustor are predicted accurately with the medium grid. The other three shock waves occur in the core flow of the combustor. The location of the pressure rise due to these three shock waves could not be predicted accurately due to the presence of recirculation zone downstream of the small step. Further, the effect of length-to-depth ratio of the cavity and the back pressure on the wave structure in the combustor has been investigated. The obtained results show that there is an optimal length-to-depth ratio for the cavity to restrict the movement of the shock wave train in the flow field of the scramjet combustor. The low velocity region in the cavity affects the downstream flow field for low back pressure. The intensity of the shock wave generated at the exit of the isolator depends on the back pressure at the exit of the combustor and this in turn affects the pressure distribution on the upper wall of the combustor.     

On the evolution of weak discontinuities in radiative magnetogasdynamics

A method of wavefront analysis is used to analyze the formation of shock waves in a two-dimensional steady supersonic flow past plane and axisymmetric bodies in radiative magnetogasdynamics. The gas is assumed to be perfectly conducting and to be permeated by a magnetic field orthogonal to the trajectories of the gas particles. The medium is taken to be sufficiently hot for the effects of thermal radiation to be significant, which is treated by the optically thin approximation to the radiative transfer equation. Transport equations, which lead to the determination of the distance at which the first characteristic could intersect with a successive one and also to conditions, which insure that no shock will ever evolve on the wave front, are derived. The effect of upstream flow Mach number and the magnetic field strength on the behavior of the wavefront are examined.     

Hypersonic flows over slender circular cones at small angles of attack

This paper presents a perturbation theory for hypersonic flows past pointed-nose slender bodies of revolution at small angles of attack. The theory presents the counter part of other theories on two-dimensional flow, axisymmetric flow, and flow past delta wings, in the case of bodies of revolution. Further restricting the analysis to Newtonian flow, a straightforward perturbation in the angle of attack is applied to the equations obtained and the resulting equations are solved only for circular cones. A striking feature of this approach is the absence of a vortical layer and a uniformally valid solution at body surface and all over the flowfield. In spite of the yaw angle, conical streamlines at cone surface are predicted which bend around towards the leeward plane. Results obtained for the surface pressure and the shock wave of a circular cone compare very well with other approximate calculations and experiment.     

On flows with vibrational relaxation behind an attached shock wave

The flow conditions are represented by an oblique shock attached to the vertex of a straight edged wedge, placed symmetrically with respect to a uniform supersonic flow in front, so that the flow after the shcok is along the wedge surface. The curvature of the attached shock wave, and the flow variable gradients at the wedge tip caused by non-equilibrium vibrational effects are determined for a steady two-dimensional flow of a vibrationally relaxing gas. A relation between the curvature of the shock and that of the stream line is obtained at the wedge-tip. The variation of entropy and an expression for the vorticity generated by the shock towards down-stream flow are determined. Calculations have been carried out for the shock curvature and flow variable gradients for a wide range of wedge angle and the up-stream flow Mach number.     

Drag reduction research in supersonic flow with opposing jet

A CFD study on drag reduction in supersonic flow with opposing jet has been conducted. Flowfield parameters, reattachment point position and surface pressure distributions are obtained and validated with experiments. From the analysis on the physical mechanism of drag reduction, it shows the phenomenon that, when the opposing jet blows, the high pressure region is located between the bow shock wave and the Mach disk, which makes the nose region much lower pressure. As the pressure ratio increases, the high pressure region is gradually pushed away from the surface. Larger the total pressure ratio is, the lower of the drag coefficient is. To study the effect of the intensity of opposing jet more reasonably, a new parameter R_PA has been introduced by combining the flux and the total pressure ratio. The study shows that the same shock wave position and drag coefficient can be obtained with the same R_PA with different fluxes and the total pressures, which means the new parameter could stand for the intensity of opposing jet and could be used to analyze the influence of opposing jet on flow field and aerodynamic force.     

高超声速飞行器激波控制减阻技术

针对高超声速飞行器巨大的激波阻力,采用数值方法研究了由钝头体、气动杆和侧向喷流构成的组合模型的减阻性能。侧向喷流将弓形激波推离气动杆,组合模型的再附激波明显弱于传统气动杆模型,其阻力系数比气动杆模型低了33.52%,从而验证了本文组合模型优异的减阻效率。进行了组合模型的影响因素分析,随侧向喷流总压比和气动杆的长度的增加,再附激波强度减弱,减阻效率升高,但减阻效率的变化速率逐渐减小。随喷口位置向下游移动,再附激波逐渐增强,减阻效率降低,且减阻效率的变化速率逐渐增加。此外本文还研究了以上参数对流场结构及钝头体压力峰值位置的影响。     

Growth and decay of weak discontinuities in radiative gasdynamics

Growth and decay properties of weak discontinuities headed by wave fronts of arbitrary shape in three dimensions are investigated in a thermally radiating inviscid gas flow. The effects of radiative transfer are treated by the use of a general differential approximation for a grey gas of arbitrary opacity including effects of radiative flux, pressure and energy density. The transport equations representing the rate of change of discontinuities in the normal derivatives of the flow variables are obtained, and it is found that the nonlinearity in the governing equations does not contribute anything to the radiation induced waves. In contrast to the radiation induced waves, the nonlinearity in the governing equations plays an important role in the interplay of damping and steepening tendencies of a modified gasdynamic wave. An explicit criterion for the growth and decay of a modified gasdynamic wave along bicharacteristics curve in the characteristic manifold of the governing differential equations is given and the special reference is made of diverging and converging waves. It is shown that there is a special case of a compressive converging wave for which the stabilizing influences of thermal radiation and the wave front curvature are not strong enough to overcome the tendency of the wave to grow into a shock.     

Numerical investigation on the shock wave transition in a three-dimensional scramjet isolator

The scramjet isolator, which is used to prevent the hypersonic inlet from disturbances that arise from the pressure rise in the scramjet combustor due to the intense turbulent combustion, is one of the most critical components in hypersonic airbreathing propulsion systems. Any engineering error that is possible in the design and manufacturing procedure of the experimental model, and the intense heat release in the scramjet combustor, may cause the performance of the isolator to decrease, leading to its lack of capability in supporting the back pressure. The coupled implicit Reynolds Averaged Navier–Stokes (RANS) equations and the two-equation standard k?ε turbulent model have been employed to numerically simulate the flow fields in a three-dimensional scramjet isolator. The effects of the divergent angle and the back pressure on the shock wave transition and the location of the leading edge of the shock wave train have been estimated and discussed. The obtained results show that the present numerical results are in very good agreement with the available experimental shadow-pictures, and the numerical method is more suitable for capturing the shock wave train and predicting the location of the leading edge of the shock wave train in the scramjet isolator than the present two-dimensional numerical methods. This is due to the small width-to-height ratio of the isolator and the intense three-dimensional flow structures. On increasing the divergent angle of the scramjet isolator, the static pressure along the central symmetrical line of the isolator decreases sharply. This is due to the strong expansion wave generated at the entrance of the isolator, and when the divergent angle of the isolator is sufficiently large, namely 1.5°, a zone of negative pressure is formed just ahead of the leading edge of the shock wave train. At the same time, the shock wave train varies from being oblique to being normal, and then back to oblique. With an increase in the prescribed back pressure at the exit of the scramjet isolator, the leading edge of the shock wave train moves forward towards the entrance of the isolator, and when the back pressure is sufficiently large, unstart conditions in the hypersonic inlet can take place if the shock train reaches the inlet.     

Fine structure of interplanetary shock fronts from the data of the BMSW instrument of the PLASMA-F experiment

The paper is concerned with studying the thickness of fronts of 38 interplanetary shocks detected by the BMSW instrument, which is a part of the scientific payload of the SPEKTR-R spacecraft, which was launched into a highly elliptical orbit in 2011. The main parameters of the interplanetary shocks have been calculated as follows: the ratio of thermal pressure to magnetic pressure before the front β, the angle between the shock front normal and the undisturbed magnetic field θ_Bn, the ratio of the shock propagation velocity to the magnetosonic velocity in the undisturbed region M_ms, and the shock front velocity relative to the Earth. It has been shown that the front thickness determined from the plasma parameters approximately matches the front thickness obtained from the magnetic field measurements and lies between 0.5 and 5 proton inertial lengths. In some events, the oscillations have been observed (upstream and downstream of the shock) in plasma parameters and in the magnetic field data. The length has been found to be between 0.5 and 6 proton inertial lengths for the preceding oscillations and between 0.5 and 10 proton inertial lengths for the following oscillations. The average value of the proton inertial length is 62 km.     

Statistical properties of the fluxes of energetic ions and solar wind in the region of intersection of shock wave fronts in the distant heliosphere

The character of statistical distributions of the intensity of energetic charged particles, solar wind flux, and the interplanetary magnetic field strength is analyzed using the data obtained by the Voyager 1 and Voyager 2 spacecraft in the distant heliosphere. A comparison of the distributions in the region of crossings of shock wave fronts in 1991 and in 2004 is carried out, and their similarities and differences are discussed.     

Thin body motion under free surface with formation of final length cavity

The two-dimensional problem of high velocity moving slender body interaction with fluid free surface was regarded. Gravity was neglected as compared with fluid inertia. The problem is relevant to underwater motion of a bullet or shell, which could take place in orbital collisions of particles with fluid-filled containments. The problem is also relevant to evaluating wave-breakers effect on the streaming flow. The solution was obtained for the case of final length cavity formation behind the body. The solution allows determining free surface shape, drag and lift forces. In the limiting cases of relatively small depth or high compressibility, the obtained solution permitted analytical approximation.     

激波干扰对发汗冷却影响的数值模拟研究

基于宏观表征体元(REV)的数值模拟方法开展了激波干扰对异质发汗冷却影响的数值模拟研究,获得了外部激波干扰与引射气体边界层耦合相互作用流场特征.研究结果表明,不同冷却介质对于冷却效率有显著的影响,冷却介质比热容越大,相同注入率条件下的冷却效果更好;入射激波干扰会显著影响多孔材料表面的压力分布,使得多孔材料内部冷却介质会...     

高超声速化学平衡流辐射流场数值解

本文用数值方法研究了高超声速大钝头体绕流辐射流场。采用粘性激波层控制方程和平衡化学反应气体模型。辐射部份包括连续辐射和线辐射。采用切平板辐射输运模型解决了空间推进的问题，计算结果表明：辐射传热与飞行器再入飞行速度，高度及头部半径密切相关，在计算的条件下，当再入速度小于７．６公里／秒时，辐射传热的影响可以忽略，当再入速度为１１公里／秒左右时，辐射传热可以与对流传热相比拟。当再入速度更大时，辐射传热远     

激波在狭缝中绕射过程研究

针对固体火箭发动机药柱存在裂纹的情况，对运动激波在狭缝中绕射传播的非定常流场进行了计算，得到了清晰的流场变化图谱，从流场图中可分辨出激波在狭缝中的传播过程。对不同深宽比的狭缝进行研究，发现流场结构存在明显差别，在深宽比较小时狭缝内压力波振荡较小，当深宽比超过一定值时狭缝内的压力波产生剧烈振荡。研究结果表明，该方法可很好地描述激波在狭缝中的传播过程。     

超声速进气道喘振的机理研究

应用数值模拟方法对中心锥中心进气混压式进气道的喘振现象进行了研究。在数值计算的基础上,根据进气道出口截面每个网格点的压力、密度、速度等参数计算了进气道喘振过程中流量系数和总压恢复系数随时间的变化情况。同时给出了在喘振过程中激波振荡的振幅、频率、对应的波系图案。并根据进气道头部分离涡的发展情况以及进气道内通道中状态参数的变化情况对喘振产生的机理进行了分析,认为进气道头部分离涡对喘振的产生起关键的作用。     

Similarity solutions for oscillating pointed-nose slender axisymmetric bodies — part I: circular cones

Similarity solutions for pointed-nose slender circular cones performing small-amplitude pitching oscillatons in hypersonic flow at zero angles of attack, are obtained in this paper. A three-dimensional unsteady hypersonic flow perturbation theory found recently by the present author for the pitching oscillations of such bodies is further studied here. It is shown that the equations can be reduced to four sets of boundary-value problems on ordinary differential equations by seeking a perturbation solution in the reduced frequency parameter k. Closed-form simple formulas for the stability derivatives are obtained showing the effect of all the parameters included in the problem except k. The obtained results show good agreement with the few available results.     

乘波体设计方法及设计参数分析

乘波体是一类典型的高超声速飞行器外形,它具有较高的升力及升阻比特性。建立了一套基于激波装配法的乘波体设计方法,该方法采用激波装配法获得精确激波位置,同时引入导波体的概念,通过设计导波体的构型来决定激波面的形状,引入导波体构型参数。通过改变导波体俯视外形前后锥角度、前后体长度比、侧视外形表面后锥角度等参数,比较其对最终乘波体气动特性的影响,并分析了这些参数改变对乘波体气动特性影响的内在机制,为工程化乘波体设计提供理论指导。     

Moving of a blunt body through the dense atmosphere under conditions of severe aerodynamic heating and ablation

Consideration is given for radiative-convective heat transfer in the vicinity of stagnation streamline in the bow shock layer of blunt body moving through the dense atmosphere at velocities up to 50 km/sec. It is shown that radiative fluxes on non-ablating body surface are subjected to an approximate similarity relations.The effects of gaseous ablation layer on radiative fluxes reaching the body surface is investigated. The results obtained are used to treat the available data on trajectories and light curves of bright fireballs.The treatment procedure proposed is shown to provide reliable estimates of the major meteorite parameters such as their mass, density and size.     

斜激波串与上游激波干扰特性试验研究

为研究斜激波串在背压条件下前移与上游激波相互干扰的流场结构和运动规律,在来流为马赫数 2.7 的直管道内设计一种等宽度斜楔,采用动态压力测量、高速纹影和粒子图像测速(PIV)技术等手段进行了试验。研究结果表明：内置斜楔在管道内产生入射激波、分离激波、膨胀波、再附激波和激波诱导分离等复杂上游激波流场,在分离区附近形成有顺压梯度和逆压梯度的区域。当增大下游压比时,斜激波串逐渐向上游激波流场移动;经过斜楔产生的分离区时,斜激波串的移动速度急剧提升,同时出现非对称分离偏转方向的切换。对比了三种长度尺寸的等楔角斜楔所产生的上游激波流场的差异性,发现在相同的斜楔前缘起始点和楔角时,随着斜楔长度的增加,上游激波流场中激波诱导的分离尺度逐渐变大。     

多级压缩锥导乘波体设计与分析

提出了多级压缩锥导乘波体的设计方法,该方法应用吻切锥理论和零攻角圆锥绕流基准流场通过流线追踪生成具有多个压缩面的乘波体。对以吸气式冲压发动机为动力的高超声速飞行器,应用多级压缩乘波前体可充分发挥前体的预压缩作用,为进气道的正常工作提供所需的均匀流场。以二级压缩乘波体为例阐述了该设计方法,设计方法通过对二级压缩基准流场进行重构,使其符合Taylor-Maccoll流动模型以获得新的二级压缩基准流场。同时编写设计程序生成了一级、二级和三级压缩乘波体,通过数值模拟结果校验设计方法的正确性,并对其压缩性、升阻比、总压恢复系数等性能进行了对比分析。