Numerical simulation of hypersonic viscous flow around space shuttle with method of diffusion analogy

Hypersonic viscous flow around a space shuttle with M(infinity) = 7, Re = 148000 and angle of attack alpha = 5-degrees is simulated numerically with the special Jacobian matrix splitting technique and simplified diffusion analogy method. With the simplified diffusion analogy method the efficiency of computation and resolution of the shock can be improved.

Receptivity to free-stream disturbance waves for blunt cone axial symmetry hypersonic boundary layer

Based on high-order compact upwind scheme, a high-order shock-fitting finite difference scheme is studied to simulate the generation of boundary layer disturbance waves due to free-stream waves. Both steady and unsteady flow solutions of the receptivity problem are obtained by resolving the full Navier-Stokes equations. The interactions of bow-shock and free-stream disturbance are researched. Direct numerical simulation (DNS) of receptivity to free-stream disturbances for blunt cone hypersonic boundary layers is performed.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.

Shape Optimization for Hypersonic Arc-Wing Missiles

National Natural Science Foundation of China (NO.90916013)

Experimental Studies on Scramjet Tested in a Freejet Facility

Two different type scramjet models with side-wall compression and top-wall compression inlets have been tested in HPTF (Hypersonic Propulsion Test Facility) under the experimental conditions of Mach number 5.8, total temperature 1700K, total pressure 4.5MPa and mass flow rate 3.5kg/s. The liquid kerosene was used as main fuel for the scramjets. In order to get fast ignition in the combustor, a small amount of hydrogen was used as a pilot. A strut with alternative tail was employed for increasing the compression ratio and for mixing enhancement in the side-wall compression case. Recessed cavities were used as a flameholder for combustion stability. The combustion efficiency was estimated by one dimensional theory. The uniformity of the facility nozzle flow was verified by a scanning pitot rake. The experimental results showed that the kerosene fuel was successfully ignited and stable combustion was achieved for both scramjet models. However the thrusts were still less than the model drags due to the low combustion efficiencies.

Development And Application Of Flexible Substrate Sensors In Instantaneous Heat Flux Measurement

A new type of sensor with the flexible substrate is introduced. It is applicable in measuring instantaneous heat flux on the model surface in a hypersonic shock tunnel. The working principle, structure and manufacture process of the sensor are presented. The substrate thickness and the dynamic response parameter of the sensor are calculated. Because this sensor was successfully used in measuring the instantaneous heat flux on the surface of a flat plate in a detonation-driven shock tunnel, it may be effective in measuring instantaneous heat flux on the model surface.

乘波构形和乘波飞行器研究综述

乘波构形的特点是高升阻比，下表面上的流动是均匀的，因此是推进系统/机身一体化设计的理想候选构形．乘波飞行器是源于乘波构形的高超音速飞行器，利用了乘波构形的高升阻比，并可为吸气发动机提供已知的均匀流场．本文比较全面地总结了乘波构形的生成方法和乘波飞行器的设计方法，介绍了乘波构形的优化方法及影响因素，给出了优化的乘波构形，并介绍了乘波飞行器的研究进展，提出了今后的研究重点．

扩压段对高超声速推进风洞起动的影响

高超声速推进风洞是进行超燃冲压发动机模型地面模拟的重要试验设备.其中,扩压段的设计非常重要,它不是孤立的,而与主流系统和引射系统密切相关.通过对不同几何形状的扩压段在不同安装位置进行调试,研究在风洞起动过程中整个系统流动状态的演变规律,探索扩压段的优化设计.实验表明,扩压段应安装在距主喷管较近的位置,并具有与主喷管相匹配的入口面积和形状.扩压段的形状应设计成能使主气流通过一系列斜激波串减速的通道.扩压段应具有一定的长度以保证主气流减速到一定程度.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

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.

Supersonic dusty-gas flows with Knudsen effect in interphase momentum exchange

On the basis of the two-continuum model of dilute gas-solid suspensions, the dynamic behavior of inertial particles in supersonic dusty-gas flows past a blunt body is studied for moderate Reynolds numbers, when the Knudsen effect in the interphase momentum exchange is significant. The limits of the inertial particle deposition regime in the space of governing parameters are found numerically under the assumption of the slip and free-molecule flow regimes around particles. As a model problem, the flow structure is obtained for a supersonic dusty-gas point-source flow colliding with a hypersonic flow of pure gas. The calculations performed using the full Lagrangian approach for the near-symmetry-axis region and the free-molecular flow regime around the particles reveal a multi-layer structure of the dispersed-phase density with a sharp accumulation of the particles in some thin regions between the bow and termination shock waves.

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.

Diffusion analogy and shock capturing for solving aerodynamic equations

Improving the resolution of the shock is one of the most important subjects in computational aerodynamics. In this paper the behaviour of the solutions near the shock is discussed and the reason of the oscillation production is investigated heuristically. According to the differential approximation of the difference scheme the so-called diffusion analogy equation and the diffusion analogy coefficient are defined. Four methods for improving the resolution of the shock are presented using the concept of diffusion analogy.

高超声速边界层流动转捩研究

针对高超声速飞行器前体预压缩性所需求的气动构型,开展了具有一级压缩效果的压缩面边界层的流动稳定性分析.采用有限体积法数值求解NS方程组得到基本流场,应用当地局部平行流假设和线性稳定性理论求解了扰动波参数的特征值问题.分析了来流马赫数M∞=6情况下二维扰动波的演化规律,并进一步关联扰动的空间放大率,结合EN方法进行了转捩预测.

扫描式总压耙在高超声速推进风洞流场测量中的应用

在风洞实验中,为了保证实验结果的可靠性,首先需要了解流场的品质.笔者自行设计研制了用于高超声速推进风洞流场测量的带有水冷装置的可移动式扫描总压耙.对于出口截面为300mm×187mm的风洞喷管,通过计算机程序控制,可在3s时间内实现全截面间歇式或连续式扫描,最大移动速度可达250mm/s,而且定位准确.通过扫描结果,分析了流场压力均匀性、稳定性以及实验结果的可重复性,同时还给出了风洞喷管出口截面的总压与马赫数等值线图.从而为超燃冲压模型发动机实验提供参考数据.