Numerical investigation of the effect of rotational relaxation rate on the nitrogen shock wave structure

The nitrogen rotational relaxation process is numerically investigated on the basis of the Monte-Carlo method and the calculation of the classic three-dimensional equations of motion of rigid-rotator molecules. The dependence of the relaxation parameterZ _R on the initial values of the translational and rotational temperaturesT _t andT _R and the asymmetry parameter in the molecular interaction potential is obtained. The problem of the shock wave structure is solved using the Pullin phenomenological model whose parameters are determined on the basis of the solution of the relaxation problem with exact calculation of the molecular interaction. The dependence of the flow fields and the relative thickness of the density profile on the quantityZ _R is obtained and the calculated data are compared with the results of experiments. It is shown that the data are in better agreement for larger values of the asymmetry parameter.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 159–167, July–August, 1995.     

Experimental and numerical studies of rotational relaxation behind a strong shock wave in air

This paper describes the experimental and numerical investigations of unknown characteristics of the rotational nonequilibrium phenomena behind a strong shock wave in air. Experiments were carried out using a piston-driven shock tube with helium as driving gas and air as driven (test) gas, operated as a two-stage shock tube. In the experiments, emission spectra of NO were measured to evaluate the rotational temperature behind a strong shock wave. The numerical calculations use the computational code for the thermal and chemical nonequilibrium flow behind a strong shock wave developed by the present author's group, where 11 chemical species (N, O, NO, N, O, N, O, NO, N, O, e) and 48 chemical reactions of high-temperature air are considered. The thermal nonequilibrium is expressed by introducing an 8 temperature model composed of translational temperature, rotational and vibrational temperatures for N, O, NO, and electron temperature. The coupling of a rotation, vibration and dissociation (CRVD) model was incorporated to take sufficiently into account the rotational nonequilibrium. The calculations were conducted for the same conditions as the experimental ones. From the calculated flow properties, emission spectra were re-constructed using the code for computing spectra of high temperature air “SPRADIAN”. Furthermore, rotational and vibrational temperatures of NO (0,1) were determined from a curve fitting method and compared with the computed results. Received 12 September 2001 / Accepted 18 February 2002     

Criterial analysis of the effect of vibrations of an airfoil surface element on the transonic flow structure

A criterial analysis of the effect of forced vibrations of airfoil surface elements on the shock wave structure of a transonic flow around the airfoil is performed. The parameter responsible for regimes of interaction of vibrationally moving zones of the airfoil with the closing shock wave is determined. The influence of this parameter on the wave drag of the airfoil is studied.     

Stability of the Hypersonic Shock Layer on a Flat Plate

The stability of hypersonic viscous gas flow in a shock layer in the neighborhood of a flat plate is considered. The stability of the velocity, temperature, density, and pressure profiles calculated on the basis of the complete viscous shock layer equations is investigated within the framework of the linear stability theory with allowance for the shock wave relations. The calculated perturbation growth rates and phase velocities are compared with the experimental data obtained by means of electron-beam fluorescence.     

Characteristic features of boundary-layers in vibration–dissociation non-equilibrium

The interaction between vibrational excitation and chemical reactions in high-temperature flows is examined through an example of a typical boundary-layer, i.e. the thermal layer generated by a shock wave reflecting at the end-wall of a shock-tube. During the development of this boundary-layer in space and time, particular features are pointed out, like the existence of a freezing zone close to the wall, responsible of dissociation rate constants remaining significant in this zone where the translational temperature is decreasing. The catalytic role of the wall is also examined and its influence is important on the vibrational temperature profiles and therefore on the rate constants themselves.     

Sound absorption in a shock wave

It was shown in [1–4] that the reflection of a sound wave or its transmission through a shock front should be accompanied by attenuation or intensification of the wave is regarded as a discontinuity. In accordance with current representations [5, 6], a shock wave includes a viscous shock and a lengthy relaxation zone. Equilibrium is established with respect to translational and rotational degrees of freedom in the viscous shock and with respect to internal degrees of freedom in the relaxation zone. The result of the interaction of the shock and sound waves is determined by the relationship between the length of the sound wave and the width of the shock wave.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 90–94, May–June, 1986.     

Shock wave interaction in a dusty gas and the appearance of fully dispersed waves

A problem of regular (symmetric and asymmetric) interaction of plane shock waves in a steady-state dusty-gas flow is considered. The possibility of the formation of wave structures is revealed, in which either all or some of the incident or reflected waves degenerate into fully dispersed waves, i.e. zones in which the parameters of both phases vary continuously. Using the Rankine-Hugoniot relations for a one-velocity “effective-gas” model, the ranges of nondimensional governing parameters (the Mach number, the angles between the incident waves and the free stream, the phase specific-heat ratio, and the particle mass concentration) are found, which correspond to different wave configurations. In the framework of a two-fluid dusty-gas model, the flow structure in the region of symmetric interaction of the shocks is calculated numerically for typical configurations containing fully dispersed waves. The flow in the region of a normal fully dispersed wave is also calculated. Good agreement between the calculated wave structure and the data known in the literature is obtained. A range of governing parameters in which the carrier-phase temperature has a local maximum inside the wave structure is found.     

Numerical simulation of hypersonic flow around a sharp cone

A hypersonic viscous flow in the shock layer near sharp cones is considered. The profiles of density and velocity, the slopes of the shock wave, and the pressures and heat fluxes calculated on the basis of the full viscous shock layer equations are compared with available experimental and theoretical data.     

Computational study of curved shock–vortex interaction

The interaction between a curved shock wave and a compressible vortex is numerically investigated. The investigation concentrates on the local deformation of the shock structure due to the shock–vortex interaction. The essentially non‐oscillatory (ENO) scheme is used to solve the unsteady two‐dimensional Euler equations. A curved shock wave is obtained by the diffraction of an initially planar shock wave around a right‐angled corner and then allowed to interact with a strong compressible vortex superimposed on the flow. The same vortex affects the shock wave differently depending on the placement of the vortex because of the varying strength of the shock wave. This effect could range from a non‐symmetric deformation of the shock wave to a local disruption in the shock structure depending on the strength of the shock wave in the interaction region. This process leading to a local disruption in the shock structure is analyzed in detail. It is shown that such a disruption in the shock structure can be predicted by simple one‐dimensional considerations. Copyright © 1999 John Wiley & Sons, Ltd.     

Aspects of self-similar diffraction gas flows and their numerical simulation

A large number of papers has been devoted to the investigation of the interaction of a plane shock wave with bodies of various geometric shapes, and they have been generalized and classified for a stationary body in [1, 2]. Separate results of experimental and theoretical investigations of the interaction of a shock wave with a wedge, cone, sphere, and cylinder moving with supersonic velocities are contained in [3–9]. Analysis of the available results shows that the features of the unsteady gas flows formed in this case largely depend on the nature of the boundary-value problem that arises for the system of differential gas dynamic equations. The question of the wave structure of the unsteady gas flow and the accuracy of the obtained solution is central to the numerical investigation of the present class of problems. The most characteristic types of unsteady self-similar gas flows that arise on the interaction of a plane shock wave with bodies of a wedge or convex corner type are calculated on the basis of an explicit numerical continuous calculation method of the second order of accuracy. The accuracy of the numerical solutions is discussed on the basis of a comparison with the experimental data. The case of the interaction of a shock wave with the rarefaction wave that arises in a supersonic flow past a convex corner is considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 146–152, July–August, 1986.     

Aerodynamic focusing of inertial particles in the shock-wave intersection region

The flow structure in the region of intersection of two stationary plane shock waves in a dusty gas is investigated. The regular symmetric and asymmetric shock interaction and the symmetric Mach interaction are considered. For a small particle mass concentration, on the basis of numerical calculations using the full Lagrangian approach it is shown that behind the shock wave intersection point a long thin region is formed, in which the particle trajectories intersect and the particle concentration sharply increases. A parametric study of the particle concentration distributions in this region is performed and the range of governing parameters on which the particle focusing effect is maximal is found.     

Development of the Richtmyer-Meshkov instability during interaction of the diffusion mixing layer of two gases with transient and reflected shock waves

The evolution of the initially disturbed mixing layer of two gases of different densities under the action of an incident shock wave, shock waves reflected from the end face, compression and expansion waves is studied in a two-dimensional unsteady approximation on the basis of the previously formulated mathematical model of mechanics of a two-velocity two-temperature mixture of gases. Problems of wave interaction with a sinusoidally disturbed diffusion layer are solved numerically. It is shown that the calculated width of the mixing region is in good agreement with experimental data.     

A three-dimensional numerical study on instability of sinusoidal flame induced by multiple shock waves

The instabilities of a three-dimensional sinusoidally premixed flame induced by an incident shock wave with Mach = 1.7 and its reshock waves were studied by using the Navier–Stokes(NS) equations with a single-step chemical reaction and a high resolution, 9th-order weighted essentially non-oscillatory scheme. The computational results were validated by the grid independence test and the experimental results in the literature. The computational results show that after the passage of incident shock wave the flame interface develops in symmetric structure accompanied by large-scale transverse vortex structures. After the interactions by successive reshock waves, the flame interface is gradually destabilized and broken up, and the large-scale vortex structures are gradually transformed into small-scale vortex structures. The small-scale vortices tend to be isotropic later.The results also reveal that the evolution of the flame interface is affected by both mixing process and chemical reaction. In order to identify the relationship between the mixing and the chemical reaction, a dimensionless parameter, η, that is defined as the ratio of mixing time scale to chemical reaction time scale, is introduced. It is found that at each interaction stage the effect of chemical reaction is enhanced with time.The enhanced effect of chemical reaction at the interaction stage by incident shock wave is greater than that at the interaction stages by reshock waves. The result suggests that the parameter η can reasonably character the features of flame interface development induced by the multiple shock waves.     

关于超音速流中运动激波若干问题的研究

研究了在无粘完全气体流中的运动激波 ,讨论了激波运动速度D和来流速度U对激波后气流参数的影响 ,包括对激波后的总焓比值和总压比值以及对流转角的影响。计算结果表明它们不同于通常静止激波下所得到的结果。该内容涉及到超音速射流与障碍物或空腔体相互作用时出现的失稳状态下激波的振动和空腔体底部的反常加热问题。     

Simulated configurational temperature of particles and a model of constitutive relations of rapid-intermediate-dense granular flow based on generalized granular temperature

In gas–solid flat-base spout bed with a jet, the flow of particles must go through an intermediate regime where both kinetic/collisional and frictional contributions play a role. In this paper, the statistical framework is proposed to define the generalized granular temperature which sums up the configurational temperature and translational granular temperature. The configurational temperature, translational and rotational granular temperatures of particles are simulated by means of CFD-DEM (discrete element method) in a 3D flat-base spout bed with a jet. The configurational temperatures of particles are calculated from instantaneous overlaps of particles. The translational and rotational granular temperatures of particles are calculated from instantaneous translational and angular velocities of particles. Roughly, the simulated translational and rotational granular temperatures increase, reach maximum, and then decrease with the increase of solids volume fractions. However, the configurational temperature increases with the increase of solids volume fractions. At high solid volume fraction, the predicted configurational temperatures are larger than the translational and rotational granular temperatures, indicating that the rate of energy dissipation do contributes by contact deformation of elastic particles. The generalized granular temperature is proposed to show the relation between the variance of the fluctuation velocity of deformation and the variance of the translational fluctuation velocity of particles. The constitutive relations of particle pressure, viscosity, granular conductivity of fluctuating energy and energy dissipation in rapid-intermediate-dense granular flows are correlated to the generalized granular temperature. The variations of particle pressure, shear viscosity, energy dissipation and granular conductivity are analyzed on the basis of generalized granular temperature in a flat-base spout bed with a jet. The axial velocities of particles predicted by a gas–solid two-fluid model of rapid-intermediate-dense granular flows agree with experimental results in a spout bed.     

Momentum and Energy Transfer in a Shock Wave

The macroparameter profiles in a strong shock wave propagating in a single-component monatomic gas are investigated. The interaction between the molecules is described by the variable-diameter sphere model. Qualitative information concerning the shock wave parameters is obtained by direct statistical simulation (Monte-Carlo) method. Using the data obtained, simple approximate gas-velocity dependences of the stress and the heat flow are found and the linear Newton and Fourier relations are generalized.     

一种冲击波作用下结构毁伤算法研究

针对爆炸冲击波与建筑物结构相互作用过程,分析了冲击波与结构碎块作用机理,发展了一种能够模拟建筑物结构破坏及冲击波传播过程的计算模型和方法。采用建筑物结构工程毁伤载荷作为判据,处理结构在冲击波作用下的破坏问题;利用流固耦合界面算法处理结构运动引起的泄压效应,利用“虚拟网格通气技术”处理结构碎块对冲击波的阻碍作用,模拟了冲击波作用下典型建筑物的毁伤过程及冲击波传播过程。结果表明,该模型在模拟冲击波与结构的作用过程中,压力计算结果与非结构动网格模拟结果符合较好;在典型建筑物毁伤过程的数值模拟中,计算得到的建筑物毁伤效果和冲击波超压分布与建筑物物理毁伤特点符合。     

Slip and Macroparameter Jumps of a Polyatomic Gas with Rotational Degrees of Freedom on a Weakly Curved Spherical Surface

The problem of the slip of a temperature-inhomogeneous polyatomic gas along a spherical surface of small curvature is solved. The solution is obtained using the half-space moment method on the basis of a previously proposed model kinetic equation which takes into account the rotational degrees of freedom of the polyatomic gas. Both the first- and second-order (in the Knudsen number) slip coefficients and the polyatomic gas macroparameter jump coefficients on the phase interface are obtained. These coefficients are given as functions of the tangential momentum accommodation coefficients, the translational and rotational energy accommodation coefficients, and the Prandtl number. The kinetic coefficients are calculated for certain polyatomic gases.     

Radiative cooling of plasma behind a reflected shock front

Calculation of gas flow in a shock tube on the basis of ideal theory [1] leads to results that differ from the real picture. In particular, the calculated velocity of the reflected shock wave exceeds the experimentally measured velocity [2] by about 20%. The calculated parameters of shock-heated gas agree well with the experimental results only directly behind the shock front [3]. The present paper reports a theoretical and experimental investigation of the variation of the plasma parameters behind the front of a reflected shock wave in argon. A picture of the gas-dynamic processes taking place after reflection of the incident shock wave by the end of the shock tube is determined. A method is developed for approximate analytic calculation, this making it possible to determine not only the parameters of the gas directly behind the front of the reflected shock wave for different positions of the wave relative to the end of the shock tube but also the variation of these parameters in other regions behind the reflected shock wave. The calculation takes into account the influence of the boundary layer and radiative cooling in the approximation of a low degree of ionization of the plasma and persistence of equilibrium conditions in the entire region behind the reflected shock wave. The experimental and theoretical profiles of the radiation behind the reflected shock wave are compared.