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.     

Unsteady reflection of a shock wave from a body with a cylindrical recess

In a number of cases of supersonic flow past bodies with recesses pulsations in the flow arise [1–3]. Experiments [4, 5] indicate that stabilization of the steady supersonic flow past the body with a recess on which a shock wave is incident takes place after a series of oscillations of the bow wave. Numerical calculation of the interaction of a supersonic jet with a cylindrical cavity [6] reveals that damped pressure pulsations arise inside the cavity if the jet is homogeneous, and undamped pulsations it is inhomogeneous. The authors explain the damping of the pulsations by the influence of artificial viscosity. This paper investigates experimentally and theoretically (by numerical methods) the oscillations of the bow shock wave and the parameters of the flow behind it in the case of unsteady reflection of a shock wave from a body with a cylindrical recess turned towards the flow. The problem is posed as follows. A plane shock wave with constant parameters impinges on a cylinder with a cavity. The unsteady flow originating from this interaction is investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 199–202, September–October, 1984.     

Supersonic motion of bodies in a gas with shock waves

The authors consider the problem of supersonic unsteady flow of an inviscid stream containing shock waves round blunt shaped bodies. Various approaches are possible for solving this problem. The parameters in the shock layer on the axis of symmetry have been determined in [1, 2] by using one-dimensional theory. The authors of [3, 4] studied shock wave diffraction on a moving end plane and wedge, respectively, by the through calculation method. This method for studying flow around a wedge with attached shock was also used in [5]. But that study, unlike [4], used self-similar variables, and so was able to obtain a clearer picture of the interaction. The present study gives results of research into the diffraction of a plane shock wave on a body in supersonic motion with the separation of a bow shock. The solution to the problem was based on the grid characteristic method [6], which has been used successfully to solve steady and unsteady problems [7–10]. However a modification of the method was developed in order to improve the calculation of flows with internal discontinuities; this consisted of adopting the velocity of sound and entropy in place of enthalpy and pressure as the unknown thermodynamic parameters. Numerical calculations have shown how effective this procedure is in solving the present problem. The results are given for flow round bodies with spherical and flat (end plane) ends for various different values of the velocities of the bodies and the shock waves intersected by them. The collision and overtaking interactions are considered, and there is a comparison with the experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 141–147, September–October, 1984.     

Shock wave in a gas-liquid medium

The propagation of weak shock waves and the conditions for their existence in a gas-liquid medium are studied in [1]. The article [2] is devoted to an examination of powerful shock waves in liquids containing gas bubbles. The possibility of the existence in such a medium of a shock wave having an oscillatory pressure profile at the front is demonstrated in [3] based on the general results of nonlinear wave dynamics. It is shown in [4, 5] that a shock wave in a gas-liquid mixture actually has a profile having an oscillating pressure. The drawback of [3–5] is the necessity of postulating the existence of the shock waves. This is connected with the absence of a direct calculation of the dissipative effects in the fundamental equations. The present article is devoted to the theoretical and experimental study of the structure of a shock wave in a gas-liquid medium. It is shown, within the framework of a homogeneous biphasic model, that the structure of the shock wave can be studied on the basis of the Burgers-Korteweg-de Vries equation. The results of piezoelectric measurements of the pressure profile along the shock wave front agree qualitatively with the theoretical representations of the structure of the shock wave.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 65–69, May–June, 1973.     

The reflection of a plane shock wave from a body with groove

When a plane shock wave impinges on bodies with grooves and when a supersonic stream of gas flows past such bodies a complicated flow pattern develops. In a number of cases oscillations of the bow wave [1–3] and an anomalous heating of the gas in the groove [4–6] have been observed. Unsteady reflection of shock waves from bodies with grooves and the processes occurring inside the grooves have been investigated comparatively little.Translated from Izvestiya Akademii Nauk SSSR, Hekhanika Zhidkosti 1 Gaza, No. 5, pp. 180–186, September–October, 1935.The authors wish to thank V. I. Ivanov for carrying out the calculations.     

Existence of a new type of irregular shock-wave interaction

Several types of plane shock-wave interactions [1–5], realized in various cases of practical importance, are now known. In [7] the possibility of the existence of a new type of shock-wave interaction was demonstrated by numerically solving the axisymmetric boundary-value problem. Here, the corresponding two-dimensional boundary-value problem of the interaction between a shock wave and a plane shock is numerically studied, a theoretical basis is obtained for the region of existence of the shock-wave interaction detected, and its principal properties are investigated. Moscow, Dnepropetovsk. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 181–183, September–October, 1988.     

Experimental investigation of supersonic flow past blunt bodies with strong distributed injection

The entry of bodies into planetary atmospheres at high supersonic velocities is accompanied by intense evaporation of the surface due to radiative heat fluxes. A series of problems involving the conduction of investigations of such kind has been proposed by Petrov. In [1], in particular, the entry of a meteorite into an atmosphere was examined. The gasdynamic aspects of this problem have been approximately simulated by many authors by intense injection of gas in theoretical, e.g., [2–5], and experimental [6, 7] studies. The theoretical studies were based on two-layer [3, 4] or three-layer [5] schemes of gas flow between the shock wave and the surface of the body. The aim of the present work was an experimental investigation of the interaction of injection with a counter supersonic flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 84–95, May–June, 1978.     

Investigation of the properties of a shock wave arising in a supersonic flow of plasma,passing through a transverse magnetic field

In a flow of plasma, set up by an ionizing shock wave and moving through a transverse magnetic field, under definite conditions there arises a gasdynamic shock wave. The appearance of such shock waves has been observed in experimental [1–4] and theoretical [5–7] work, where an investigation was made of the interaction between a plasma and electrical and magnetic fields. The aim of the present work was a determination of the effect of the intensity of the interaction between the plasma and the magnetic field on the velocity of the motion of this shock wave. The investigation was carried out in a magnetohydrogasdynamic unit, described in [8]. The process was recorded by the Töpler method (IAB-451 instrument) through a slit along the axis of the channel, on a film moving in a direction perpendicular to the slit. The calculation of the flow is based on the one-dimensional unsteady-state equations of magnetic gasdynamics. Using a model of the process described in [9], calculations were made for conditions close to those realized experimentally. In addition, a simplified calculation is made of the velocity of the motion of the above shock wave, under the assumption that its front moves at a constant velocity ahead of the region of interaction, while in the region of interaction itself the flow is steady-state.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 86–91, January–February, 1975.     

Transition from regular reflection to mach reflection when a shock wave interacts with a wall in a two-phase gas—liquid medium

A study is made of the transition from regular reflection to Mach reflection when a plane moderately strong or weak shock wave interacts with a wall in a two-phase gas—liquid medium. An equilibrium model that differs from the model of Parkin et al. [1] by the introduction of the adiabatic velocity of sound is used to investigate shock wave reflection in the complete range of gas concentrations. For the reflection of weak shock waves, nonlinear asymptotic expansions [2] are used. In the limiting cases, the results agree with those already known for single-phase media [2, 3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 190–192, September–October, 1983.     

Flow in the neighborhood of the stagnation line on a blunt body traveling at supersonic speed through a zone with elevated temperature and vibrational energy of the molecules

The unsteady flow in the neighborhood of the stagnation line on a sphere traveling at supersonic speed through a plane layer of diatomic gas with elevated temperature and nonequilibrium excitation of the molecular vibrations is investigated. (The source of the inhomogeneity could be a gas discharge [1].) The problem is solved using the viscous shock layer model which makes it possible to take molecular transport processes into account and analyze the unsteady heat transfer. Such flows were previously calculated in [2] within the framework of the inviscid gas model.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i GazNo. 3, pp. 183–185, May–June, 1990.     

Shock strengthening in a wedge-shaped cavity

The unsteady problem of the entry of a shock wave of arbitrary intensity into a wedge-shaped cavity is examined. An exact solution of the non-linear problem of reflection of a plane wave from a nonplanar wall is found for certain cavity angles. Numerical wave focusing calculations are carried out for arbitrary cavity angles. A single scaling law is obtained for gas flows with waves of moderate and high intensity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 123–129, September–October, 1987.     

Hypersonic flow-past of plane blunt bodies by a nonviscous radiating gas

An analytic solution is obtained in the work in a Newtonian approximation [1] for the flow-past problem for a plane blunt body by a steady-state uniform hypersonic inviscous space-radiating gas flow. The hypersonic flow-past problem for axisymmetrical blunt bodies by a nonviscous space-radiating gas has been previously considered [2–4]. In this case a satisfactory solution of the problem was obtained even in a zero-th approximation by decomposing the unknown values in terms of a parameter equal to the ratio of gas densities before and after passage of the shock wave. The solution of the problem in a zero-th approximation with respect to in the case of flow-past of plane blunt bodies does not turn out to be satisfactory, since the departure of the shock and the radiant flux to the body as gas flows into the shock layer turns out to be strongly overstated under nearly adiabatic conditions. Freeman [5] demonstrated that results may be significantly improved for flow-past of a plane blunt body by a nonradiating gas if a more precise expression is used for the tangential velocity component expressed in a new approximation with respect to the parameter . This refinement is applied in this work for solving the flow-past problem for a plane blunt body by a space-radiating gas. The distribution of the gasdynamic parameters in the shock layer, the departure of the shock wave, and the radiant heat flux to the surface of the body are found. The solution obtained is analyzed in detail for the example of flow-past regarding a circular cylinder.Translated from Zhurnal Prikladnoi Mekhanikii Tekhnicheskoi Fiziki, No. 3, 68–73, May–June, 1975.     

Regular reflection of a magnetohydrodynamic shock wave from a conducting wall

In two-dimensional supersonic gasdynamics, one of the classical steady-state problems, which include shock waves and other discontinuities, is the problem concerning the oblique reflection of a shock wave from a plane wall. It is well known [1–3] that two types of reflection are possible: regular and Mach. The problem concerning the regular reflection of a magnetohydrodynamic shock wave from an infinitely conducting plane wall is considered here within the scope of ideal magnetohydrodynamics [4]. It is supposed that the magnetic field, normal to the wall, is not equal to zero. The solution of the problem is constructed for incident waves of different types (fast and slow). It is found that, depending on the initial data, the solution can have a qualitatively different nature. In contrast from gasdynamics, the incident wave is reflected in the form of two waves, which can be centered rarefaction waves. A similar problem for the special case of the magnetic field parallel to the flow was considered earlier in [5, 6]. The normal component of the magnetic field at the wall was equated to zero, the solution was constructed only for the case of incidence of a fast shock wave, and the flow pattern is similar in form to that of gasdynamics. The solution of the problem concerning the reflection of a shock wave constructed in this paper is necessary for the interpretation of experiments in shock tubes [7–10].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 102–109, May–June, 1977.The author thanks A. A. Barmin, A. G. Kulikovskii, and G. A. Lyubimov for useful discussion of the results obtained.     

Solution of the direct problem of the mixed flow of a conducting gas in a laval nozzle with a meridional magnetic field

We consider the direct problem in the theory of the axisymmetric Laval nozzle (including sonic transition) for the steady flow of an inviscid and nonheat-conducting gas of finite electrical conductivity. The problem is solved by numerical integration of the equations of unsteady gas flow using an explicit difference scheme that was proposed by Godunov [1,2], and was used to calculate steady and unsteady flows of a nonconducting gas in nozzles by Ivanov and Kraiko [3]. The subsonic and the supersonic flows of a conducting gas in an axisymmetric channel when there is no external electric field, the magnetic field is meridional, and the magnetic Reynolds numbers are small have previously been completely investigated. Thus, Kheins, Ioller and Élers [4] investigated experimentally and theoretically the flow of a conducting gas in a cylindrical pipe when there is interaction between the flow and the magnetic field of a loop current that is coaxial with the pipe. Two different approaches were used in the theoretical analysis in [4]: linearization with respect to the parameter S of the magnetogasdynamic interaction and numerical calculation by the method of characteristics. The first approach was used for weakly perturbed subsonic and supersonic flows and the solutions obtained in analytic form hold only for small S. This is the approach used by Bam-Zelikovich [5] to investigate subsonic and supersonic jet flows through a current loop. The numerical calculations of supersonic flows in a cylindrical pipe in [4] were restricted to comparatively small values of S since, as S increases, shock waves and subsonic waves appear in the flow. Katskova and Chushkin [6] used the method of characteristics to calculate the flow of the type in the supersonic part of an axisymmetric nozzle with a point of inflection. The flow at the entrance to the section of the nozzle under consideration was supersonic and uniform, while the magnetic field was assumed to be constant and parallel to the axis of symmetry. The plane case was also studied in [6]. The solution of the direct problem is the subject of a paper by Brushlinskii, Gerlakh, and Morozov [7], who considered the flow of an electrically conducting gas between two coaxial electrodes of given shape. There was no applied magnetic field, and the induced magnetic field was in the direction perpendicular to the meridional plane. The problem was solved numerically in [7] using a standard process. However, the boundary conditions adopted, which were chosen largely to simplify the calculations, and the accuracy achieved only allowed the authors [7] to make reliable judgments about the qualitative features of the flow. Recently, in addition to [7], several papers have been published [8–10] in which the authors used a similar approach to solve the direct problem in the theory of the Laval nozzle (in the case of a nonconducting gas).Translated from Izvestiya Akademiya Nauk SSSR, Mekhanika Zhidkosti i Gaza., No. 5, pp. 14–20, September–October, 1971.In conclusion the author wishes to thank M. Ya. Ivanov, who kindly made available his program for calculating the flow of a conducting gas, and also A. B. Vatazhin and A. N. Kraiko for useful advice.     

Study of unsteady flow past bodies behind shock waves

Several theoretical and experimental studies have been devoted to the problem of the nonstationary action of the stream behind a shock wave on bodies of varied shape. In particular, in [1], the pressure and density are calculated for flow about bodies of the more typical shapes in the initial stage of the process. The basic relations which accompany the interaction of shock waves are considered in [2, 3]. The analysis of the phenomena of diffraction of shock waves on the sphere, cylinder, and cone is presented in [4]. Problems of unsteady flow about a wing are examined in [5, 6]. A detailed review of the foreign studies on unsteady flow is given in [7]. Of great practical interest is the question of the time for flow formation and the magnitudes of the unsteady loads during this period. Experimental investigations have been made recently [8, 9] in which some criteria are presented for estimating the bow shock formation time for supersonic flow about the sphere and cylinder with flat blunting. However the question of the formation time of the stationary pressure on the body surface is not referred to in these studies and no relationship is shown between the transient position of the reflected wave and the corresponding unsteady pressure on the surface. Moreover, in [8] the dimensionless time criterion is determined very approximately, independently of the Mach number of the shock wave. The present study was undertaken with the object of determining the basic criteria which characterize unsteady flow about bodies behind a plane shock wave which has time-independent parameters, and clarification of the shock wave reflected from the body and the pressure on the surface of the body during the transient period. The most typical body shapes were studied: 1) a cylinder with flat face aligned with the stream; 2) a spherically-blunted cylinder; and 3) a cylinder transverse to the stream. The experiments were conducted in a conventional shock tube using the single-diaphragm scheme. The measurements of the pressure on the models and the velocity of the incident shock wave were made using the technique analogous to that of [10, 11]. A highspeed movie camera was used to record the pattern of the wave diffraction on the body. The Mach number of the incident shock wave varied in the range from M=1.5 to M≈6.0, which corresponded to a range of Mach numbers M_∞ of the stream behind the shock wave from 0.6 to 2.1. The calculations of the required gas dynamic parameters for high temperatures were made with account for equilibrium dissociation of the air on the basis of the data of [10, 12, 13]. The magnitude of the relative maximal shock wave standoff Δ at the stagnation point obtained in the present experiments was compared with the values of Δ from other studies. In the case of the flat-blunted cylinder it was in good agreement with the results of [8–14], and in the case of the spherically-blunted cylinder and the transverse cylinder it was in agreement with the results of [15].     

Distribution of radiant thermal flux over the surface of a sphere for hypersonic flow of a nonviscous radiating gas past the sphere

In the present study using the Newtonian approximation [1] we obtain an analytical solution to the problem of flow of a steady, uniform, hypersonic, nonviscous, radiating gas past a sphere. The three-dimensional radiative-loss approximation is used. A distribution is found for the gasdynamic parameters in the shock layer, the withdrawal of the shock wave and the radiant thermal flux to the surface of the sphere. The Newtonian approximation was used earlier in [2, 3] to analyze a gas flow with radiation near the critical line. In [2] the radiation field was considered in the differential approximation, with the optical absorption coefficient being assumed constant. In [3] the integrodifferential energy equation with account of radiation was solved numerically for a gray gas. In [4–7] the problem of the flow of a nonviscous, nonheat-conducting gas behind a shock wave with account of radiation was solved numerically. To calculate the radiation field in [4, 7] the three-dimensional radiative-loss approximation was used; in [5, 6] the self-absorption of the gas was taken into account. A comparison of the equations obtained in the present study for radiant flow from radiating air to a sphere with the numerical calculations [4–7] shows them to have satisfactory accuracy.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 44–49, November–December, 1972.In conclusion the author thanks G. A. Tirskii and É. A. Gershbein for discussion and valuable remarks.     

Calculation of hydrogen-air combustion behind detached shock wave for supersonic flow past a sphere

Many of the published theoretical studies of quasi-one-dimensional flows with combustion have been devoted to combustion in a nozzle, wake, or streamtube behind a normal shock wave [1–6].Recently, considerable interest has developed in the study of two-dimensional problems, specifically, the effective combustion of fuel in a supersonic air stream.In connection with experimental studies of the motion of bodies in combustible gas mixtures using ballistic facilities [7–9], the requirement has arisen for computer calculations of two-dimensional supersonic gas flow past bodies in the presence of combustion.In preceding studies [10–12] the present author has solved the steady-state problem under very simple assumptions concerning the structure of the combustion zone in a detonation wave.In the present paper we obtain a numerical solution of the problem of supersonic hydrogen-air flow past a sphere with account for the nonequilibrium nature of eight chemical reactions. The computations encompass only the subsonic and transonic flow regions.The author thanks G. G. Chernyi for valuable comments during discussion of the article.     

Irregular reflection of a strong shock wave from a thin wedge

The problem of irregular reflection of a strong shock wave from a rigid wall has been studied [1–3] mainly within the framework of the linear theory. It has been found that near the front of a shock wave there exist a region of large gradients of gasdynamic parameters in which the linear theory is no longer valid [4]. In the present paper we consider the nonlinear problem of Mach reflection when there is interaction between a shock wave of high intensity and a thin wedge. The solution of the problem is constructed on the assumption that the ratio of densities along the front of the impinging shock wave is small [5, 6].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 55–61, September–October, 1974.In conclusion, the authors wish to express their gratitude to A. A. Grib for his interest in the subject and to L. A. Rumyantsev for his help in carrying out the calculations.     

Interaction between the boundary layer and a supersonic flow when part of the wall is rapidly heated

There have been many theoretical studies of aspects of the unsteady interaction of an exterior inviscid flow with a boundary layer [1–9]. The mathematical flow models obtained in these studies by the method of matched asymptotic expansions describe a wide range of phenomena observed experimentally. These include boundary layer separation near the hinge of a flap, the flow in the neighborhood of the trailing edge of an oscillating airfoil [1–2], and the development and propagation of perturbations in a boundary layer excited by an oscillating wall or some other way [3–5]. The present paper studies the interaction of an unsteady boundary layer with a supersonic flow when a small part of the surface of a body in the flow is rapidly heated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 66–70, January–February, 1984.     

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.