The flow gradients in the vicinity of a shock wave for a thermodynamically imperfect gas

Supersonic rotational planar and axisymmetric flows of a non-viscous, non-heat-conductive gas with arbitrary thermodynamic properties in the vicinity of a steady shock wave are studied. The differential equations describing the gas flow upstream and downstream of the discontinuity surface and the dynamic compatibility conditions at this discontinuity are used. The gas flow non-uniformity in the shock vicinity is described by the spatial derivatives of the gasdynamic parameters at a point on the shock surface. The parameters are the gas pressure, density, and velocity vector. The derivatives with respect to the directions of the streamline and normal to it, and of the shock surface and normal to it, are considered. Spatial derivatives of all gasdynamic parameters are expressed through the flow non-isobaric factor along the streamline, the streamline curvature, and the flow vorticity and non-isoenthalpy factors. An algorithm for determining these factors of the gas flow downstream of a shock wave is developed. Example calculations of these factors for imperfect oxygen and thermodynamically perfect gas are presented. The influence coefficients of the upstream flow factors on the downstream flow factors are calculated. The gas flow in the vicinity of the shock is described by the isolines of gasdynamic parameters. Uniform planar and axisymmetric flows at different distances from the axis of symmetry are examined; the isobars, isopycnics, isotachs and isoclines are used to characterize the downstream flow behind a curved shock in an imperfect gas.     

Two-dimensional self-similar flows generated by the interaction between a shock and low-density gas regions

A plane time-dependent flow generated by the interaction between a normal shock and a low-density gas region occupying a quarter of the plane is theoretically investigated. Numerical simulation is performed on the basis of the Euler equations. It is established that after the shock has come in contact with the low-density region two-dimensional self-similar flows of different type can develop. On regular interaction the original shock is refracted on the low-density region with the matching of the accelerated and original shock and the refracted contact discontinuity at a common point. On irregular interaction a complicated flow occurs; it includes curved and oblique shocks, a contact discontinuity with points of inflection, multiple matching points, a high-pressure jet, and a layered vortex. The jet and vortex structures are investigated in detail. The tendency of the gasdynamic structure development with variation in the control parameters of the problem is determined. A simplified, near-analytical technique for estimating the slopes of the main shocks and the gas parameters behind them is proposed.     

Numerical solution of the problem of wave diffraction by a wedge

The systematic development of the theory of shock reflection from a solid wall started in [1]. Regular reflection and a three-shock configuration originating in Mach reflection were considered there under the assumption of homogeneity of the domains between the discontinuities and, therefore, of rectilinearity of these latter. The difficulties of the theoretical study include the essential nonlinearity of the process as well as the instability of the tangential discontinuity originating during Mach reflection. Analytic solutions of the problem in a linear formulation are known for a small wedge angle or a weak wave (see [2–4], for example). The solution in a nonlinear formulation has been carried out numerically in [5, 6] for arbitrary wedge angles and wave intensities. Since the wave was nonstationary, the internal flow configuration is difficult to clarify by means of the constant pressure and density curves presented. A formulation of the problem for the complete system of gasdynamics equations in self-similar variables is given in [7] and a method of solution is proposed but no results are presented. Difficulties with the instability of the contact discontinuity are noted. The problem formulation in this paper is analogous to that proposed in [7]. However, a method of straight-through computation without extraction of the compression shocks in the flow field is selected to compute the discontinuous flows. The shocks and contact discontinuities in such a case are domains with abrupt changes in the gasdynamics parameters. The computations were carried out for a broad range of interaction angles and shock intensities. The results obtained are in good agreement with the analytical solutions and experimental results. Information about the additional rise in reflection pressure after the Mach foot has been obtained during the solution.     

Interaction Between a Slow Magnetohydrodynamic Shock Wave and a Tangential Discontinuity

The self-similar problem of the oblique interaction between a slow MHD shock wave and a tangential discontinuity is solved within the framework of the ideal magnetohydrodynamic model. The constraints on the initial parameters necessary for the existence of a regular solution are found. Various feasible wave flow patterns are found in the steady-state coordinate system moving with the line of intersection of the discontinuities. As distinct from the problems of interaction between fast shock waves and other discontinuities, when the incident shock wave is slow the state ahead of it cannot be given and must to be determined in the process of solving the problem. As an example, a flow in which the slow shock wave incident on the tangential discontinuity is generated by an ideally conducting wedge located in the flow is considered. The basic features of the developing flows are determined.     

Solution of Riemann problem for dusty gas flow

A direct approach is used to solve the Riemann problem for a quasilinear hyperbolic system of equations governing the one dimensional unsteady planar flow of an isentropic, inviscid compressible fluid in the presence of dust particles. The elementary wave solutions of the Riemann problem, that is, shock waves, rarefaction waves and contact discontinuities are derived and their properties are discussed for a dusty gas. The generalised Riemann invariants are used to find the solution between rarefaction wave and the contact discontinuity and also inside rarefaction fan. Unlike the ordinary gasdynamic case, the solution inside the rarefaction waves in dusty gas cannot be obtained directly and explicitly; indeed, it requires an extra iteration procedure. Although the case of dusty gas is more complex than the ordinary gas dynamics case, all the parallel results for compressive waves remain identical. We also compare/contrast the nature of the solution in an ordinary gasdynamics and the dusty gas flow case.     

Oblique collision of a plane-polarized Alfvén discontinuity and a slow shock wave propagating in opposite directions in magnetohydrodynamics

Collision of plane fronts of a plane-polarized Alfvén discontinuity and a slow shock wave propagating in opposite directions at a certain angle is considered within the framework of an ideal magnetohydrodynamic model. The initial state of an infinitely conducting medium at rest with a frozen-in magnetic field is assumed to be given. Calculations are carried out for various values of the shock wave Mach number and the magnetic field strength using a special software which makes it possible to find an exact solution of the Riemann problem of breakdown of a discontinuity between the states downstream of the interacting waves by means of a computer. The wave flow structure is investigated and a bifurcation map of flow restructuring is constructed. Domains of the initial parameters for which the interaction differs qualitatively are distinguished. The parameters of the medium and magnetic field are found as functions of the angle between the colliding discontinuities and the inclination of the magnetic field. The results obtained may be used in investigations of magnetic reconnection.     

Evolution of weak shocks in one dimensional planar and non-planar gasdynamic flows

Asymptotic decay laws for planar and non-planar shock waves and the first order associated discontinuities that catch up with the shock from behind are obtained using four different approximation methods. The singular surface theory is used to derive a pair of transport equations for the shock strength and the associated first order discontinuity, which represents the effect of precursor disturbances that overtake the shock from behind. The asymptotic behaviour of both the discontinuities is completely analysed. It is noticed that the decay of a first order discontinuity is much faster than the decay of the shock; indeed, if the amplitude of the accompanying discontinuity is small then the shock decays faster as compared to the case when the amplitude of the first order discontinuity is finite (not necessarily small). It is shown that for a weak shock, the precursor disturbance evolves like an acceleration wave at the leading order. We show that the asymptotic decay laws for weak shocks and the accompanying first order discontinuity are exactly the ones obtained by using the theory of non-linear geometrical optics, the theory of simple waves using Riemann invariants, and the theory of relatively undistorted waves. It follows that the relatively undistorted wave approximation is a consequence of the simple wave formalism using Riemann invariants.     

Application of unsteady mixing equations to some aerodynamic problems

Tangential discontinuities [1] are introduced in solving several transient and steady-state problems of gas dynamics. These discontinuities are unstable [2] as a result of the effects of viscosity and thermal conductivity. Therefore it is advisable to replace the tangential discontinuity by a mixing region and account for its interaction with the inviscid flows, establishing on the boundaries of this region the conditions of vanishing friction stress and equality of the velocity and temperature components to the corresponding velocity and temperature components of the inviscid flows. This formulation improves the accuracy of the solution of such problems by posing them as problems with irregular reflection and intersection of shock waves [1].The consideration of the interaction of unsteady turbulent mixing regions with the inviscid flow also permits the formulation of several problems in which the effects of viscosity lead to complete rearrangement of the flow pattern (the lambda-configuration) with the interaction of the reflected shock wave with the boundary layer in the shock tube [3,4], the formation of zones of developed separation ahead of obstacles, etc.).In this connection, §1 presents an analysis of the self-similar solutions of the unsteady turbulent mixing equations (a corresponding analysis of the laminar mixing equations which coincide with the boundary layer equations is presented in [1]). It is shown that these self-similar solutions describe, along with the several problems noted above, the problems of the formation of steady jets and mixing zones in the base wake.As an example, §2 presents, within the framework of the proposed schematization, an approximate solution of the problem of the interaction of a shock wave reflected from a semi-infinite wall with the boundary layer on a horizontal plate behind the incident shock wave. The results obtained are applied to the analysis of reflection in a shock tube. Computational results are presented which are in qualitative agreement with experiment [3, 4].     

Decay of a plane shock wave in a two-parameter medium with arbitrary equation of state

Special curves, called shock polars, are frequently used to determine the state of the gas behind an oblique shock wave from known parameters of the oncoming flow. For a perfect gas, these curves have been constructed and investigated in detail [1]. However, for the solution of problems associated with gas flow at high velocities and high temperatures it is necessary to use models of gases with complicated equations of state. It is therefore of interest to study the properties of oblique shocks in such media. In the present paper, a study is made of the form of the shock polars for two-parameter media with arbitrary equation of state, these satisfying the conditions of Cemplen's theorem. Some properties of oblique shocks in such media that are new compared with a perfect gas are established. On the basis of the obtained results, the existence of triple configurations in steady supersonic flows obtained by the decay of plane shock waves is considered. It is shown that D'yakov-unstable discontinuities decompose into an oblique shock and a centered rarefaction wave, while spontaneously radiating discontinuities decompose into two shocks or into a shock and a rarefaction wave.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 147–153, November–December, 1982.     

Some geometric acoustic problems of one-dimensional unsteady and two-dimensional steady flows

The behavior of discontinuities (weak shocks) of the parameters of a disturbed flow and their interaction with the discontinuities of the basic flow in the geometric acoustics approximation, when the variation of the intensity of such shocks along the characteristics or the bicharacteristics is described by ordinary differential equations, has been investigated by many authors. Thus, Keller [1] considered the case when the undisturbed flow is three-dimensional and steady, and the external inputs do not depend on the flow parameters. An analogous study was made by Bazer and Fleischman for the MGD isentropic flow of an ideal conducting medium [2], while Lugovtsov [3] studied the three-dimensional steady flow of a gas of finite conductivity for small magnetic Reynolds numbers and no electric field. Several studies (for example, [4]) have considered the behavior of discontinuities of the solutions from the general positions of the theory of hyperbolic systems of quasilinear equations. Finally, the interaction of weak shocks (or the equivalent continuous disturbances) with shock waves was studied in [5–11].In what follows we consider one-dimensional (with plane, cylindrical, and spherical waves) and quasi-one-dimensional unsteady flows, and also plane and axisymmetric steady flows. Two problems are investigated: the variation of the intensity of weak shocks in the presence of inputs which depend on the stream parameters, and the interaction of weak shocks with strong discontinuities which differ from contact (tangential) discontinuities.The thermodynamic properties of the gas are considered arbitrary. We note that the resulting formulas for the interaction coefficients of the weak and strong discontinuities are also valid for nonequilibrium flow.     

Equation of a shock wave front

Second-order differential equations of the hyperbolic type are derived for describing the local law of shock wave propagation. The shock waves are assumed to be two-dimensional unsteady in a stationary gas flow and three-dimensional steady in a supersonic flow. The behavior of the characteristics of these equations is investigated as a function of the governing flow parameters and their relative position with respect to the typical bicharacteristics of the characteristic cone behind the shock is analyzed. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 159–165, May–June, 2000.     

Interaction between a shock wave and a longitudinal low-density gas layer

The problem of the interaction between a shock wave and a semi-infinite longitudinal plane layer or a cylindrical channel of finite thickness filled with a low-density gas is studied on the basis the Euler equations. The flow gasdynamics, including qualitatively new, regular and irregular, interaction regimes, are described. New gasdynamic flow elements, such as high-pressure jets with an internal wave structure and layered vortices, are found to exist. It is revealed that the gasdynamic precursor growth is decelerated at long time intervals, due to the flow chocking effect and the vorticity development behind its front.     

Propagation of discontinuities against a static background

The solution of the ideal gasdynamic equations describing propagation of a shock wave initiated, for example, by the motion of a piston against an inhomogeneous static background is considered. The solution is constructed in the form of Taylor series in a special time variable which is equal to zero on the shock wave. In the case of weak shock waves divergence of the series serves as the constraint for such an approach. Then the solution is constructed by linearizing the equations about the solution with a weak discontinuity. In the case of a given background the last solution can be always found exactly by solving successively a set of transport equations, all these equations are reduced to linear ordinary differential equations. The presentation begins from the one-dimensional solutions with plane waves and ends by discussion of spatial problems.     

Impact of the interplanetary magnetic field on thewave flow pattern in the neighborhood of the Earth’s bow shock under sharp variations in the solar wind dynamic pressure

The impact of the interplanetary magnetic field on transformation and disintegration of the Earth’s bow shock into a system of magnetohydrodynamic (MHD) shock waves, rotational discontinuities and rarefaction waves under the action of abrupt variations in the solar wind dynamic pressure is simulated in the three-dimensional non-plane-polarized formulation within the framework of the ideal magnetohydrodynamic model using the solution of the MHD Riemann problem of breakdown of an arbitrary discontinuity. This discontinuity arises when a contact discontinuity, on which the solar wind density increases or decreases suddenly and which travels together with the solar wind, impinges on the Earth’s bow shock and propagates along its surface. The interaction pattern is constructed in the quasisteady- state formulation as a mosaic of exact solutions obtained on computer using an original MHD Riemann solver. The wave flow patterns are found for all elements of the surface of the bow shock as functions of their latitude and longitude for various jumps in the density on the contact discontinuity and characteristics parameters of the solar wind and interplanetary magnetic field at the Earth’s orbit. It is found that when the solar wind dynamic pressure increases, a fast MHD shock wave, that first penetrates into the magnetosheath, is always formed. When the solar wind dynamic pressure decreases, the influence of the interplanetary magnetic field can lead to the development of the leading fast MHD shock wave in certain zones on the surface of the Earth’s bow shock. The solution obtained can be used to interpret measurements on spacecraft in the solar wind at the libration point and in the neighborhood of the Earth’s magnetosphere.     

Optimum overtaking compression shocks with restrictions imposed on the total flow-deflection angle

The problem of optimization of gasdynamic variables behind a system of two steady oblique compression shocks with restrictions imposed on the flow-defection angle is considered. The intervals of input parameters, in which this system turns out to be more effective than one shock, are determined. On the basis of an analysis of the system optimal for the static pressure, the physical meaning of the transition from one type of the reflected discontinuity to another is explained for the problem of interaction of overtaking oblique compression shocks. Baltic State Technical University, St. Petersburg 198005. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 4, pp. 99–108, July–August, 1999.     

Inviscid vortex structures in the shock layers of conical flows around V wings

The applicability of the criteria of existence of inviscid vortex structures (vortex Ferri singularities) is studied in the case in which a contact discontinuity of the corresponding intensity proceeds from the branching point of the λ shock wave configuration accompanying turbulent boundary layer separation under the action of an inner shock incident on the leeward wing panel. The calculated and experimental data are analyzed, in particular, those obtained using the special shadow technique developed for visualizing supersonic conical streams in nonsymmetric, Mach number 3 flow around a wing with zero sweep of the leading edges and the vee angle of 2π /3. The applicability of the criteria of existence of inviscid vortex structures is established for contact discontinuities generated by the λ shock wave configuration accompanying turbulent boundary layer separation realized under the action of a shock wave incident on the leeward wing panel. Thus, it is established that the formation of the vortex Ferri singularities in a shock layer is independent of the reason for the existence of the contact discontinuity and depends only on its intensity.     

Eulerian formulation of transport equations for three-dimensional shock waves in simple elastic solids

The propagation of a three-dimensional shock wave in an elastic solid is studied. The material is assumed to be a simple elastic solid in which the Cauchy stress depends on the deformation gradient only. It is shown that the growth or decay of a discontinuity ψ depends on (i) an unknown quantity φ^? behind the shock wave, (ii) the two principal curvatures of the shock surface, (iii) the gradient on the shock surface of the shock wave speeds and (iv) the inhomogeneous term which depends on the motion ahead of the shock surface and vanishes when the motion ahead of the shock surface is uniform. If a proper choice is made of the propagation vectorb along which the growth or decay of the discontinuity is measured, the dependence on item (iii) can be avoided. However,b assumes different directions depending on the choice of discontinuity ψ with which one is concerned and the unknown quantity φ^? behind the shock wave on which one chooses to depend. As in the case of one-dimensional shock waves, the growth (or decay) of one discontinuity may not be accompanied by the growth (or decay) of other discontinuities. A universal equation relating the growth or decay of discontinuities in the normal stress, normal velocity and specific volume is also presented.     

基于多层分块算法的激波干扰流场预测

激波-激波干扰流场预测是超声速乃至高超声速流动中最具挑战性的问题之一. 特别地, 第IV类激波干扰由于其在壁面驻点附近产生极高的热载荷而备受关注. 本文针对圆柱诱导的弓形激波和入射斜激波的干扰问题, 分别基于量热完全气体模型和考虑振动激发的热完全气体模型, 数值求解有黏二维可压缩NS方程, 分析了高温气体效应对激波干扰流场结构, 以及第IV类激波干扰流场状态参数的影响. 接着, 本文基于一种具有广义可分离特性的遗传算法 (多层分块算法), 给出能够预测不同气体模型下第IV类激波干扰流场三波点的坐标位置、超声速射流的几何形状等特征性几何结构的数学模型, 进一步获得高温气体效应对激波干扰类型转变准则影响的定量化评估. 激波干扰类型转变准则面上的多组临界工况的激波干扰流场结构以及壁面压力和壁面热流分布的对比结果表明, 不同气体模型下的激波干扰类型和流场结构差异较为显著, 获得的定量化预测模型对工程中气动热环境的预测具有一定的参考价值.      

Behavior of the vorticity vector in supersonic axisymmetric swirl flows behind a detonation wave

The behavior of the vorticity vector on a discontinuity surface arising in a supersonic nonuniform combustible gas flow with the formation of a shock or detonation wave is studied. In the general case, it is a vortex flow with prescribed distributions of parameters. It is demonstrated that the ratio of the tangential component of vorticity to density remains continuous in passing through the discontinuity surface, while the quantities proper become discontinuous. Results calculated for flow vorticity behind a steady-state detonation wave in an axisymmetric supersonic flow of a combustible mixture of gases are presented. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 15–21, November–December, 2007     

Hypersonic flow of nonequilibrium air over blunt bodies

A numerical method is described for computing nonequilibrium three-dimensional supersonic flow of a gas in the shock layer over the forward surface of blunt bodies with discontinuities of shape. The basic idea is to divide the original system of differential equations into two subsystems, which are solved in succession: first for the gasdynamic variables, the velocity components and the pressure, and then for the relaxation parameters and the enthalpy. To calculate the velocity components and the pressure we use the iterative marching method [1, 2] in the form given in [3]. The relaxation equations and the enthalpy equation are integrated numerically along the stream lines. A discussion is given of the effect of nonequilibrium of physical and chemical reactions on the distribution of parameters in the inviscid shock layer and on the aerodynamic coefficients of blunt bodies in hypersonic air flow. The unsteady aerodynamic coefficients are calculated by the curved body method [4]. The computational algorithm takes the form of a program in “ALGOL-60” for the BéSM-6 computer.