Asymptotic solution to the problem of hypersonic flow over bodies in the neighborhood of the separation point of a thin shock layer

A study is made of the problem of hypersonic flow of an inviscid perfect gas over a convex body with continuously varying curvature. The solution is sought in the framework of the asymptotic theory of a strongly compressed gas [1–4] in the limit M when the specific heat ratio tends to 1. Under these assumptions, the disturbed flow is situated in a thin shock layer between the body and the shock wave. At the point where the pressure found by the Newton-Buseman formula vanishes there is separation of the flow and formation of a free layer next to the shock wave [1–4]. The singularity of the asymptotic expansions with respect to the parameter ₁ = ( –1)/( + 1) associated with separation of the strongly compressed layer has been investigated previously by various methods [3–9]. Local solutions to the problem valid in the neighborhood of the singularity have been obtained for some simple bodies [3–7]. Other solutions [7, 9] eliminate the singularity but do not give the transition solution entirely. In the present paper, an asymptotic solution describing the transition from the attached to the free layer is constructed for a fairly large class of flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 99–105, January–February, 1982.     

Influence of vibration-dissociation coupling on heat transfer and hypersonic drag

A numerical investigation is carried out within the framework of the multicomponent total viscous shock layer model [6, 7], according to which when Re 100 the flow near a blunt body can be divided into a shock wave zone and a viscous shock layer. At the inner edge of the shock wave the generalized Rankine-Hugoniot relations are imposed, and in the shock layer the complete system of viscous shock layer equations is solved with allowance for vibrational relaxation and nonequilibrium dissociation and ionization reactions.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 141–151, May–June, 1990.     

Initial stage of reflection of a plane shock wave from a cylinder,sphere, and ellipsoid of revolution

The stage of regular reflection of a plane shock wave from a blunt body (cylinder, sphere, and ellipsoid of revolution) is considered. At the point of intersection of the reflected shock wave and the surface of the body, analytic expressions are found for the derivative of the Mach number of the wave with respect to the time, the curvature of the wave, the normal derivatives of the density and the pressure, and the derivative of the Mach number along the wave front. It is shown that the flow has a singularity at = _* < _** (s** is the limiting angle [1] of regular reflection of a shock wave from a rigid surface). The distribution of the parameters in the region between the reflected shock wave and the surface of the body is found up to terms of third order in the time. The density distribution behind the reflected shock wave was measured experimentally, and also the shape of the reflected wave at different instants of time.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 94–100, March–April, 1980.     

Stabilization of plasma flow in a transverse magnetic field

Results are given of a theoretical and experimental investigation of the intensive interaction between a plasma flow and a transverse magnetic field. The calculation is made for problems formulated so as to approximate the conditions realized experimentally. The experiment is carried out in a magneto-hydrodynamic (MHD) channel with segmented electrodes (altogether, a total of 10 pairs of electrodes). The electrode length in the direction of the flow is 1 cm, and the interelectrode gap is 0.5 cm. The leading edge of the first electrode pair is at x = 0. The region of interaction (the region of flow) for 10 pairs of electrodes is of length 14.5 cm. An intense shock wave S propagates through argon with an initial temperature T_o = 293 °K and pressure p_o = 10 mm Hg. The front S moves with constant velocity in the region x < 0 and at time t = 0 is at x = 0. The flow parameters behind the incident shock wave are determined from conservation laws at its front in terms of the gas parameters preceding the wave and the wave velocity W_S. The parameters of the flow entering the interaction region are as follows: temperature T ₀ ¹ = 10,000 °K, pressure P ₀ ¹ = 1.5 atm, conduction ₀ ¹ = 3000 ^–1·m^–1, velocity of flow u ₀ ¹ = 3000 m·sec^–1, velocity of sounda ₀ ¹ = 1600 m·sec^–1, degree of ionization = 2%, 0.4. The induction of the transverse magnetic field B = [0, B_y(x), 0] is determined only by the external source. Induced magnetic fields are neglected, since the magnetic Reynolds number Re_m 0.1. It is assumed that the current j = (0, 0, jz) induced in the plasma is removed using the segmented-electrode system of resistance R_e. The internal plasma resistance is R_i = h(A)^–1 (h = 7.2 cm is the channel height; A = 7 cm² is the electrode surface area). From the investigation of the intensive interaction between the plasma flow and the transverse magnetic field in [1–6] it is possible to establish the place x* and time t* of formation of the shock discontinuity formed by the action of ponderomotive forces (the retardation wave R_T), its velocity W_T, and also the changes in its shape in the course of its formation. Two methods are used for the calculation. The characteristic method is used when there are no discontinuities in the flow. When a shock wave R_T is formed, a system of nonsteady one-dimensional equations of magnetohydrodynamics describing the interaction between the ionized gas and the magnetic field is solved numerically using an implicit homogeneous conservative difference scheme for the continuous calculation of shock waves with artificial viscosity [2].Translated from Izvestiya Akademiya Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 112–118, September–October, 1977.     

Unsteady viscous shock layer in the case of supersonic motion through an inhomogeneity

The model of a perfect gas is used in a numerical simulation of unsteady effects in a viscous shock layer near the stagnation streamline near the front part of a rotating blunt body in an inhomogeneous external flow with pressure difference P _t8, temperature difference T _t8, and vorticity difference . The evolution of nonlinear disturbances due to the passage of a heated region and a change of the injection regime is followed. A divergence finite-difference scheme of second order of approximation across the shock layer, realized by vector sweeps with allowance for the boundary conditions on the surface of the body and behind the separated bow shock wave, is used.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 138–145, January–February, 1992.     

Unsteady supersonic flow around blunt bodies with detached shock wave

The numerical method of calculating the supersonic three-dimensional flow about blunt bodies with detached shock wave presented in [1–3] is applied to the case of unsteady flow. The formulation of the unsteady problem is analogous to that of [4], which assumes smallness of the unsteady disturbances.The paper presents some results of a study of the unsteady flow about blunt bodies over a wide range of variation of the Mach number M=1.50– and dimensionless oscillation frequency l/V=0–1.0. A comparison is made with the results obtained from the Newton theory.     

Measurement of parameters of a gas suspension behind a shock wave in foam

A study has been made of the propagation of a shock wave in dry polyhedral foam with cell diameter 1 cm. The experiments were made in a shock tube in the range of Mach numbers M < 1.4 of the shock wave. The interaction of the shock wave with the foam was photographed. This established that the destruction of the foam by the shock wave leads to the formation of a gas-droplet flow behind the shock front. To determine the parameters of the suspension, the flow was probed by He-Ne lasers with different radiation wavelengths. The spectral-transparency method was used to find the modal diameter of the droplets of the gas suspension and the volume concentration of the droplets in the flow. The modal diameter of the droplets was 2m, and the volume concentration of the droplets decreased downstream.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 134–141, May–June, 1993.     

Behavior of sonic lines in a shock layer behind a detached shock wave

A numerical investigation is made into the formation of local supersonic zones in the subsonic flow region between a detached shock wave and the surface of the body in the case of supersonic three-dimensional flow over conical bodies with opening angle _k = 120 ° of the cone in the range of Mach numbers M = 2.5–15.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No.4 pp. 143–145, July–August, 1979.We thank G. I. Petrov for suggesting the problem and for helpful advice and O. M. Belotserkovskii for constant interest in the work.     

Investigation of a “hanging” shock wave near a generating point

The hodograph method is used to plot a hanging shock wave in the plane nonequilibrium supersonic flow of an ideal gas. This paper considers the general case of an analytical solution in the plane of the hodograph at the point of generation of the shock wave. A type of limiting line is established which makes it possible to plot a shock wave (it is found that the shock wave may not extend over the whole flow, with a convolution in the physical plane).Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 30–37, November–December, 1971.     

Experimental study of plasma flow in a disk MHD channel under conditions of spontaneous formation of a current layer

The results of an experimental study of the plasma flow in a disk channel under conditions of strong hydromagnetic interaction are presented. It is shown that if the condition Re_mH ₀ ² /80.2 is satisfied for the magnetic Reynolds number at some point of the stream, then a current layer develops at that point characterized by a high electric-current density and high conductivity and temperature. The formation of the current layer leads to strong local retardation of the stream, the appearance of a shock wave, and a number of other nonlinear hydromagnetic phenomena. The experimental results are in agreement with theoretical studies conducted earlier.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 31–37, May–June, 1974.     

Instability of a hypersonic flow over a wedge in the Newtonian approximation

At the present time, there are a number of works in the literature that treat unsteady hypersonic flows in the Newtonian approximation [1–4]. Since the angle of incidence of the shock wave _s coincides in the zero-order approximation with the angle of inclination of the bodys [1], the latter is usually used in the boundary conditions on the shock. However, in the zero-order approximation _b can be used with the same justification. Both approaches are equally justified and give similar results for a steady flow. For unsteady flows the results can differ radically. It will be shown below that for an investigation of a flow over a fixed wedge with constant conditions in the free stream a steady-state pattern is obtained in the first case and a solution growing in time, in the second case.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 158–160, July–August, 1976.     

Ionization and nonequilibrium radiation of air behind strong shock waves

It is shown that at high velocities of shock waves (V 9.5 km/sec) an important factor influencing the rate of ionization is the depletion of the number of excited states of the atoms through de-excitation. In the case of low pressures (p 1 torr) and for a bounded and optically transparent region of gas heated by the shock wave (for example, for the motion of gas in a shock tube or in a shock layer near a blunt body), the effective ionization rate k_f depends on the pressure [1], which leads to violation of the law of binary similarity which holds under these conditions without allowance for de-excitation. On leaving the relaxation zone, the gas arrives at a stationary state with constant parameters differing from those in thermodynamic equilibrium. The electron concentration and also the radiation intensity in the continuum and the lines are lower than the values for thermodynamic equilibrium. These considerations explain the results of known experiments and some new experiments on ionization and radiation of air behind a travelling shock wave.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 105–112, January–February, 1980.     

Supersonic flow past circular cones at large angles of attack

Flow past sharp-nosed circular cones is investigated for a broad range of freestream Mach numbers M>1 and cone half-angles _c at angles of attack from zero to the value at which conical flow breaks down. Several new results are obtained with regard to the position of the Ferri point, the shape of the local supersonic zones and internal shock wave, and the nonmonotonicity of the windward shock slope as a function of the angle of attack. The existence of flow regimes in which the radial velocity on the windward side is directed toward the apex of the cone is demonstrated. The investigation is carried out numerically with relaxation of the solution in a fictitious time coordinate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No, 6, pp. 79–84, November–December, 1973.     

Influence of vorticity on the heat transfer to blunt bodies in hypersonic flight

When blunt bodies are in hypersonic flight, a high-entropy layer of gas with nonzero vorticity is formed near their surface. The transverse gradients of the entropy, density, and gas velocity in the layer are high, which makes it necessary to take into account its absorption by the boundary layer of finite thickness . This vortex interaction is usually accompanied by an increase in the heat flux q and the frictional stress on the wall compared with their values as calculated in accordance with the classical scheme of a thin boundary layer, when the parameters on the outer edge of the boundary layer are set equal to the inviscid parameters on the body. This effect has been investigated on the side surface of slender (with angle 1 to the undisturbed flow) blunt bodies in a hypersonic stream [1–3]. It is shown in the present paper that the effect can have a stronger and even qualitative influence on the flow over blunt bodies with 1 if the radius of curvature R_s of the detached shock wave on the axis is small compared with the midsection radius R of the body. It is shown that the distributions of the heat fluxes with allowance for the vorticity of the inviscid shock layer are similar in the case of slightly blunt (r₀/R 0) cones with half-angles less than a critical _*.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 50–57, March–April, 1981.     

Heat transfer of gas-particle flow in a supersonic convergent-divergent nozzle

Summary Heat flux, wall heat transfer coefficients, and wall pressures are determined for high velocity flow of gas-solid mixtures in a converging-diverging nozzle. Flow separation accompanied with oblique shock formation occurs in the diverging section of the nozzle. The shock strength is reduced upon the addition of solid particles. The wall pressure in the convergent section of the nozzle appears unaffected by the presence of solid particles. In the divergent section, however, the wall pressure is slightly lowered. At the maximum ratio of solid to air flow used in the experiments (3.7) increases in the heat transfer rate of up to 20 and 50 percent are obtained in the convergent and separated (divergent) regions of the nozzle, respectively. Slightly larger increases in the wall heat transfer coefficients are also obtained. It is concluded that the wall heat flux and heat transfer coefficients are influenced strongly by the presence of disturbances upstream of the nozzle inlet.Nomenclature W _a air flow rate - W _s solids flow rate - x axial distance from nozzle entrance - L axial length of nozzle - specific heat ratio of fluid - A _e exit cross section of flow - A _* throat cross section of flow - P ₀ inlet pressure - P _s wall separation pressure - P _a ambient exhaust pressure - shock wave angle - shock wave deflection angle - M ₁ Mach number upstream of shock wave - Mach number normal to shock wave - q heat flux - k _f thermal conductivity of fluid - T _wi inside wall temperature - T _wo outside wall temperature - T _ad adiabatic wall temperature - h wall heat transfer coefficient - C nozzle constant - A local cross section of flow - c _p specific heat of fluid - Pr Prandtl number - viscosity of fluid - r _c throat radius of curvature - factor accounting for variation of and Units absolute temperature °R(ankine) °F+459.7 - conductivity 1 BTU (hr ft °F)^–1 4.137×10^–3 cal (s cm °C)^–1 - specific heat 1 BTU (1b °F)^–1 1 cal (g °C)^–1 - absolute pressure 1 psia 0.0680 atm Supported in part by aid provided by the UCLA Space Science Center (Grant NsG 236-62 Libby).Listed for readers not familiar with the units adopted in this paper (editor).     

Lifting bodies using the flow behind axisymmetric conical shock waves

One of the methods of designing aircraft with supersonic flight speeds involves solving an inverse problem by means of the well-known flow schemes and the substitution of rigid surfaces for the flow surfaces. Lifting bodies using the flows behind axisymmetric shock waves belong to these configurations. All lifting bodies using the flow behind a conical shock wave can be divided into two types [1]. Bodies whose leading edge passes through the apex of the conical shock wave pertain to the first type and those whose leading edge lies below the apex of the conical shock wave, to the second. For small apex angles of the basic cone at hypersonic flow velocities an approximate solution of the variation problem was obtained, which showed that the lift-drag ratio of lifting bodies of the second type is higher than that of the first [2]. The present paper gives a numerical solution of the problem for flow past lifting bodies of the second type using the flow behind axisymmetric conical shock waves with half-angles of the basic cone _S=9.5 and 18° The upper surfaces of the bodies are formed by intersecting planes parallel to the velocity vector of the oncoming flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 135–138, March–April, 1986.     

Experimental investigation of flow over sharp cones with near-critical and supercritical vertex angle

This paper presents results of an experimental investigation of supersonic flow over sharp cones with near-critical and supercritical semivertex angles. The authors have determined the drag coefficients and the shock position at supersonic flow velocities corresponding to M = 4.0 over a range of cone semivertex angles from 40 to 130 °, at angle of attack = 0. The experimental drag coefficients are compared with available theoretical values, obtained using both exact and approximate methods of calculation. The experimentally obtained position of the attached shock wave is compared with theory, derived by the method of integral relations in the first approximation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 195–198, March–April, 1978.The authors thank G. E. Sidel'nikov for his help in processing the experimental data by means of his computer program.     

Special case of the density distribution behind a nonstationary shock wave

The density distribution behind a nonstationary shock wave for a definite value of the Mach number M_*, which depends on = c_p/c_v, is considered. Use is made of the previously established fact [1] that for M = M_*() there exists a connection between the first and second derivatives of the density along the normal behind the wave. An investigation is made into the density profile in dimensionless variables behind plane, cylindrical, and spherical shock waves in the neighborhood of the shock front. In the first case, if the gas in front of the wave is homogeneous, only two types of density profile are possible (up to small quantities of third order in the coordinate). In the second and third cases, the form of the density distribution also depends on a parameter, the ratio of the first derivative along the normal of the density behind the wave to the radius of curvature of the wave.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 163–167, November–December, 1979.     

Screening of selective radiation in a boundary layer

An analysis of numerous calculations of the flow and radiative-convective heat exchange in a hypersonic shock layer near a blunt body, both at an impermeable surface [1] and in the presence of ablation [1–4], made it possible to establish some relationships connected with the screening in the boundary layer of radiation from the high-temperature part of the shock layer. It was established that the dependence of the screening coefficient of the boundary layer, which consists of the ratio of the integral flux q⁺ _RW of radiation incident on the body to the flux q⁺ _RW arriving at the outer limit of the boundary layer from the high-temperature region behind the shock wave, on the radius has a nonmonotonic character. In this case the absolute value of the minimum _min is determined mainly by the temperature T_S behind the shock wave and by the injection velocityf _W = (v)_W/(v), while its position (R_min) is determined by the pressure P_S in the shock layer. An analysis of the spectral characteristics of the screening of the boundary layer showed that the region of wavelengths > 0.5 m is very important from the point of view of the radiative heating and destruction of the surface. Components which efficiently absorb high-temperature radiation in this spectral interval are absent from the gaseous products of the disintegration of specific heat-protection coatings. An analysis is made permitting an evaluation of the possibility of introducing into the coatings additives possessing those optical and thermochemical properties which would provide for screening of radiation in this part of the spectrum.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 71–77, September–October, 1978.The authors thank A. B. Karasev for advice and a discussion.     

Experimental investigation of the recovery of the total pressure in a hypersonic wind tunnel with conical and profiled nozzles

Measurements have been made of the coefficient of recovery of the total pressure of a gas flow exhausting from axisymmetric and conical profiled hypersonic nozzles into a cylindrical channel of diameter equal to or greater than the nozzle exit and also in the presence of an Eiffel chamber. The experiments were made at Mach numbers M = 4.83–12.4 in the isentropic core. It is shown that the values of differ slightly (by 5%) from the corresponding value for a normal shock wave at the number M determined for a onedimensional flow by the ratio of the area of the cylindrical channel to the area of the critical section of the nozzle.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 170–173, January–February, 1984.