Shocked Flows Induced by Supersonic Projectiles Moving in Tubes

A numerical study on shocked flows induced by a supersonic projectile moving in tubes is described in this paper. The dispersion-controlled scheme was adopted to solve the Euler equations implemented with moving boundary conditions. Four test cases were carried out in the present study: the first two cases are for validation of numerical algorithms and verification of moving boundary conditions, and the last two cases are for investigation into wave dynamic processes induced by the projectile moving at Mach numbers of M-p = 2.0 and 2.4, respectively, in a short time duration after the projectile was released from a shock tube into a big chamber. It was found that complex shock phenomena exist in the shocked flow, resulting from shock-wave/projectile interaction, shock-wave focusing, shock-wave reflection and shock-wave/contact-surface interactions, from which turbulence and vortices may be generated. This is a fundamental study on complex shock phenomena, and is also a useful investigation for understanding on shocked flows in the ram accelerator that may provide a highly efficient facility for launching hypersonic projectiles.

准定常强激波马赫反射波形的数值研究

应用DCD频散控制激波捕捉格式,求解二维、多组分、带有化学反应的Euler方程组,数值模拟了准定常强激波的马赫反射问题。研究结果表明:与经典马赫反射理论相比,在强激波条件下,激波诱导的气体分子振动激发和化学反应使马赫反射的三波点轨迹角变小、马赫杆高度变低、楔顶附体激波倾角变小;马赫杆的相对突出量随入射激波马赫数和楔角的增大而增大,而气体分子的振动、离解等真实气体效应能进一步加剧马赫杆的向前突出。

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.

Dns Of Compressible Turbulent Boundary Layer Around A Sharp Cone

Direct numerical simulation of the turbulent boundary layer over a sharp cone with 20 degrees cone angle (or 10 degrees half-cone angle) is performed by using the mixed seventh-order up-wind biased finite difference scheme and sixth-order central difference scheme. The free stream Mach number is 0.7 and free stream unit Reynolds number is 250000/inch. The characteristics of transition and turbulence of the sharp cone boundary layer are compared with those of the flat plate boundary layer. Statistics of fully developed turbulent flow agree well with the experimental and theoretical data for the turbulent flat-plate boundary layer flow. The near wall streak-like structure is shown and the average space between streaks (normalized by the local wall unit) keeps approximately invariable at different streamwise locations. The turbulent energy equation in the cylindrical coordinate is given and turbulent energy budget is studied. The computed results show that the effect of circumferential curvature on turbulence characteristics is not obvious.

Experimental investigation on the cavity-based scramjet model

The present work focused on improving the engine performance with different fuel equivalence ratios and fuel injections. A scramjet model with strut/cavity integrated configurations was tested under Mach 5.8 flows. The results showed that the strut may sreve as an effective tool in a kerosene-fueled scramjet. The integration of strut/cavities also had great effect on stablizing the combustion in a wide range of fuel equivalence ratio. The one-sdimensional analysis method was used to analyze the main characteristics of the model. The two-stage fuel injection should have better performance in increasing the chemical reaction rate in the first cavity region.

Analysis and application of ellipticity of stability equations on fluid mechanics

By using characteristic analysis of the linear and nonlinear parabolic stability equations (PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub-characteristic analysis of PSE, the linear PSE are proved to be elliptical and hyperbolic-parabolic for velocity U, in subsonic and supersonic, respectively; the nonlinear PSE are proved to be elliptical and hyperbolic-parabolic for relocity U + u in subsonic and supersonic, respectively. The methods are gained that the remained ellipticity is removed from the PSE by characteristic and sub-characteristic theories, the results for the linear PSE are consistent with the known results, and the influence of the Mach number is also given out. At the same time, the methods of removing the remained ellipticity are further obtained from the nonlinear PSE.

Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms = 1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock with the material interface, and effect of initial perturbation modes on R-M instability are investigated numerically. It is noted that the shock refraction is a main physical mechanism of the initial phase changing of the material surface. The multiple interactions of the reflected shock from the origin with the interface and the R-M instability near the material interface are the reason for formation of the spike-bubble structures. Different viscosities lead to different spike-bubble structure characteristics. The vortex pairing phenomenon is found in the initial double mode simulation. The mode interaction is the main factor of small structures production near the interface.

Effects of entry conditions on cracked kerosene-fueled supersonic combustor performance

Supersonic combustion of thermally cracked kerosene was experimentally investigated in two model supersonic combustors with different entry cross-section areas. Effects of entry static pressure, entry Mach number, combustor entry geometry, and injection scheme on combustor performance were systematically investigated and discussed based on the measured static pressure distribution and specific thrust increment due to combustion. In addition, the methodology for characterizing flow rate and composition of cracked kerosene was detailed. Using a pulsed Schlieren system, the interaction of supercritical and cracked kerosene jet plumes with a Mach 2.5 crossflow was also visualized at different injection temperatures. The present experimental results suggest that the use of a higher combustor entry Mach number as well as a larger combustor duct height would suppress the boundary layer separation near the combustor entrance and avoid the problem of inlet un- start.

Waverider configurations derived from general conical flowfields

A method based on the computational fluid dynamics (CFD) is presented for a flexible waverider's design. The generating bodies of this method could be any cones. In addition, either the leading edge or the profile of the scramjet's inlet is used as the waverider's definition curve, parameterized by the quadric function, the sigmoid function or the B-spline function. Furthermore, several numerical examples are carried out to validate the method and the relevant codes. The CFD results of the configurations show that all the designs are successful. Moreover, primary suggestions are proposed for practical design by comparing the geometrical and aerodynamic performances of the cone-derived waveriders at Mach 6.

A novel oblique detonation structure and its stability

Oblique detonation structures induced by the wedge in the supersonic combustible gas mixtures are simulated numerically. The results show that the stationary oblique detonation structures are influenced by the gas flow Mach number, and a novel critical oblique detonation structure, which is characterized by a more complicated wave system, appears in the low Mach number cases. By introducing the inflow disturbance, its nonstationary evolution process is illustrated and its stability is verified.

利用图象分析法求湍流剪切层的光学密度分布曲线

本文介绍利用一个Microneye Bullet摄象头,一台ApplellPlus微计算机及一台EpsonRX-80打印机从干涉条纹读出信息并由此求出剪切层中光学密度的方法,讨论了在编写计算机程序时考虑过的问题,从一张由湍流不均匀剪切流的实验中取得的Much-Zehnder干涉照片画出流场中的光学密度分布以及剪切层的扩展情况。最后还用同样方法画出激光通过剪切层后相位角偏差的均方根变化曲线。

Numerical investigation of dynamic stall of an oscillating aerofoil

The flow structure around an NACA 0012 aerofoil oscillating in pitch around the quarter-chord is numerically investigated by solving the two-dimensional compressible N-S equations using a special matrix-splitting scheme. This scheme is of second-order accuracy in time and space and is computationally more efficient than the conventional flux-splitting scheme. A 'rigid' C-grid with 149 x 51 points is used for the computation of unsteady flow. The freestream Mach number varies from 0.2 to 0.6 and the Reynolds number from 5000 to 20,000. The reduced frequency equals 0.25-0.5. The basic flow structure of dynamic stall is described and the Reynolds number effect on dynamic stall is briefly discussed. The influence of the compressibility on dynamic stall is analysed in detail. Numerical results show that there is a significant influence of the compressibility on the formation and convection of the dynamic stall vortex. There is a certain influence of the Reynolds number on the flow structure. The average convection velocity of the dynamic stall vortex is approximately 0.348 times the freestream velocity.

Numerical-simulation of low supersonic-flow around a sphere

The controlled equations defined in a physical plane are changed into those in a computational plane with coordinate transformations suitable for different Mach number M(infinity). The computational area is limited in the body surface and in the vicinities of detached shock wave and sonic line. Thus the area can be greatly cut down when the shock wave moves away from the body surface as M(infinity) --> 1. Highly accurate, total variation diminishing (TVD) finite-difference schemes are used to calculate the low supersonic flowfield around a sphere. The stand-off distance, location of sonic line, etc. are well comparable with experimental data. The long pending problem concerning a flow passing a sphere at 1.3 greater-than-or-equal-to M(infinity) > 1 has been settled, and some new results on M(infinity) = 1.05 have been presented.

High-speed streams from coronal holes and the accelerating mechanism of the solar-wind

In this paper, the general Mach number equation is derived, and the influence of typical energy forms in the solar wind is analysed in detail. It shows that the accelerating process of the solar wind is influenced critically by the form of heating in the corona, and that the transonic mechanism is mainly the result of the adjustment of the variation of the crosssection of flowing tubes and the heat source term.The accelerating mechanism for both the high-speed stream from the coronal hole and the normal solar wind is similar. But, the temperature is low in the lower level of the coronal hole and more heat energy supply in the outside is required, hence the high speed of the solar wind; while the case with the ordinary coronal region is just the opposite, and the velocity of the solar wind is therefore lower. The accelerating process for various typical parameters is calculated, and it is found that the high-speed stream may reach 800 km/sec.

爆轰驱动激波管缝合激波马赫数计算

采用数值求解的方法得出了氢氧爆轰驱动激波管的缝合状态参数。以空气为试验气体时，缝合激波马赫数随着H2摩尔浓度的增加而增加，H_2摩尔浓度达到90％左右时达到最大。当缝合马赫数较高时，需要考虑高温真实气体效应的影响，缝合激波马赫数较理想气体的高。以氢空气混合气为试验气体时，缝合激波马赫数较以空气为试验气体的小。通过调整驱动气体与被驱动气体的初始参数，可以得到即能恰好消除Taylor波又能缝合的运行状态。