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.

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

Effects of Mach number on turbulent separation behaviours induced by blunt fin

An experimental study of the interaction between shock wave and turbulent boundary layer induced by blunt fin has been carried out at M-infinity = 7.8 using oil flow visualization and simultaneous measurements of fluctuating wall pressure and heat transfer. This paper presents the effects of Mach number on turbulent separation behaviours induced by blunt fin.

Novel experimental phenomena of fine-particle fluidized beds

The dynamic response of bed height and concentration waves in liquid-solid fluidized beds to a step change in the fluidization velocity is considered. We experimentally study the liquid-solid fluidized beds, spherical beadings, with sizes ranging from 230 to 270 mesh and the inner diameter of columns made from glass is 2.4 mm. Experimental results find that under certain conditions, fine particles with large Richardson-Zaki exponent n display different dynamic behavior from usual particles with smaller n during expansion and collapse of the fluidized state. (c) 2007 Elsevier Inc. All rights reserved.

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.

Upwind compact scheme with dispersion control for aerodynamic equations

A compact upwind scheme with dispersion control is developed using a dissipation analogy of the dispersion term. The term is important in reducing the unphysical fluctuations in numerical solutions. The scheme depends on three free parameters that may be used to regulate the size of dissipation as well as the size and direction of dispersion. A coefficient to coordinate the dispersion is given. The scheme has high accuracy, the method is simple, and the amount of computation is small. It also has a good capability of capturing shock waves. Numerical experiments are carried out with two-dimensional shock wave reflections and the results are very satisfactory.

Phase relaxation model for 2-phase flows

A phase relaxation model (PRM) for 2-phase flows is presented in this paper on the basis of three principal assumptions. The basic equations for PRM arc derived from the Boltzmann equations for gas-partlcle mixture, The general characteristics and solving process of the PRM's basic equations are also presented and discussed. Many terms in the PRM's basic equations contain a factor ε= ρgρp/ρg+ρp2 which is an intrinsic small parameter for 2-phase mixture, with ρg and ρp being respectively the densities of gas and particle phases.This makes it possible to simplify the computation of the PRM's basic equations. The model is applied to for example, studying file steady propagation of shock waves in gas-particle mixture. The analysis shows that with an increase of shock wave strength the relaxation process behind a gasdynamics shock front becomes a kind of dynamics relaxation instead of the standard exponential relaxation process. A method of determining experimentally the velocity and tem...更多perature relaxation rates (or times) of gas-particle flows is suggested and analyzed.

The effect of strain softening in the matrix material during void growth

Furthermore, the compressed flow driven by the piston is discussed. The consistent solution of gasdynamical equations including solar gravity is obtained for the unsteady and two-dimensional configuration, which is applied to the region between the piston and shock wave. This solution may satisfy the jump conditions of shock wave, which separates the region of compressed flow and quiet corona.

可燃气体中激波聚焦的点火特性

数值模拟了二维平面激波从抛物面上反射在可燃气体中聚焦的过程，研究了形成爆轰波的点火特性,对理想化学当量比氢气／空气混合气体，在初始压强20kPa的条件下，马赫数2．6-2．8的激波聚焦能产生两个点火区：第1个点火区是反射激波会聚引起的，第2个点火区是由入射激波在抛物面上发生马赫反射引起的．这种条件下流场中会出现爆燃转爆轰，起爆点分别分布在管道壁面、抛物反射面和第2点火区附近．起爆机理分别为激波管道壁面反射、点火诱导激波的抛物面反射和点火诱导的激波与第2点火区产生的爆燃波的相互作用,不同的点火和起爆过程导致了不同的流场波系结构，同时影响了爆轰波传播的波动力学过程．

可燃气体中激波与障碍物作用在下游形成爆轰波的数值研究

对平面激波和单个矩形障碍物作用的过程进行了数值模拟,研究了反射产生的上行爆轰波在下游可燃气体中形成爆轰波的过程。数值结果表明,下游爆轰波形成过程主要有2种模式:爆轰波直接绕射和绕射波在上壁面反射,这和已有的实验结果是一致的。通过研究下游爆轰波的形成过程受入射激波马赫数、混合气体的压力和管道尺度的影响,分析了上游爆轰波向下游传播的波动力学过程,讨论了2种形成过程的作用规律和控制因素,阐明了下游爆轰波的形成规律。

在激波波面中氧分子的非平衡离解

对强激波作用下双原子分子振动与离解耦合的非平衡离解过程进行了理论计算.本工作的特点是将计算起点建立在分子基本参数上,采用主方程理论处理振动与离解的耦合,振动跃迁几率用SSH理论计算,在离解限附近考虑多量子数跃迁并计及原子复合的影响.对O2-Ar体系,计算给出了在正激波后O2分子振动能级分布、振动弛豫时间、离解孕育时间、离解产物浓度、离解速率系数等物理量随时间的演化.计算结果分别与Camac和Wray的实验相符.计算显示,在激波作用的后期,有准稳态的振动能级布居分布.计算结果显示,Park模型低估了非平衡离解速率系数,Hansen模型则高估了非平衡离解速率系数.

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Perturbational finite volume scheme for the one-dimensional Navier-Stokes equations

Starting from the second-order finite volume scheme,though numerical value perturbation of the cell facial fluxes, the perturbational finite volume (PFV) scheme of the Navier-Stokes (NS) equations for compressible flow is developed in this paper. The central PFV scheme is used to compute the one-dimensional NS equations with shock wave.Numerical results show that the PFV scheme can obtain essentially non-oscillatory solution.