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.     

Theory of vortical helical ideal gas flows in laval nozzles

An asymptotic solution is found for the direct problem of the motion of an arbitrarily vortical helical ideal gas flow in a nozzle. The solution is constructed in the form of double series in powers of parameters characterizing the curvature of the nozzle wall at the critical section and the intensity of stream vorticity. The solution obtained is compared with available theoretical results of other authors. In particular, it is shown that it permits extension of the known Hall result for the untwisted flow in the transonic domain [1]. The behavior of the sonic line as a function of the vorticity distribution and the radius of curvature of the nozzle wall is analyzed. Spiral flows in nozzles have been investigated by analytic methods in [2–5] in a one-dimensional formulation and under the assumption of weak vorticity. Such flows have been studied by numerical methods in a quasi-one-dimensional approximation in [6, 7]. An analogous problem has recently been solved in an exact formulation by the relaxation method [8, 9]. A number of important nonuniform effects for practice have consequently been clarified and the boundedness of the analytical approach used in [2–7] is shown.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 126–137, March–April, 1978.The authors are grateful to A. N. Kraiko for discussing the research and for valuable remarks.     

Calculation of interaction of a turbulent near-wake behind a step with a supersonic jet

Existing computational methods [1–5] do not enable one to calculate complex flows behind steps, accounting for nonuniformity of the incident supersonic flow and the effect of compression and expansion waves arriving in the near-wake region. For example, computational methods based on the methods of [1] or [2] are used mainly in uniform supersonic flow ahead of the base edge and, for the most part, cannot be used to calculate flow in annular nozzles with irregular conditions. An exception is reference [6], which investigated flow in an annular nozzle behind a cylindrical center-body. The present paper suggests a method, based on references [7, 8] for calculating the base pressure behind two-dimensional and three-dimensional steps, washed by a supersonic jet.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 43–51, November– December 1977.     

Draft of a nozzle with the circulating outflow of a vortical flow of gas

The presence of circulation in an outflowing gas leads to a change in the working parameters of a nozzle. The question of the mass flow rate and the draft of a nozzle without a diffusor (a point) for twisted flows has been studied theoretically and experimentally [1–6]. The use of nozzles with a supersonic part introduces a considerable degree of complication into the method for the analytical calculation of the draft characteristics and the program for their experimental investigation. In [2, 7], a theory of a nozzle is formulated for a model of a potential circulating flow of gas; in [5, 8], an electronic computer was used to solve the complete system of the equations of gasdynamics for the motion of a rotating flow along a nozzle; in [7, 9], an investigation was made of a variational problem of the shaping of a diffusor for a circulation flow. The calculation of the draft, carried out in the above-mentioned communications (with the exception of [2], in which a study was made of a partial model of an eddyless rotational motion), is bound up with labor-consuming computer calculations. In the present article, in a development of [3, 6], a quasi-one-dimensional theory of a supersonic nozzle for a vortical flow of gas is formulated and verified experimentally.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 145–149, September–October, 1975.     

Computation of base flow characteristics with uniform blowing

Blowing at bluff body base was considered under different conditions and for small amount of blowing this problem was solved using dividing streamline model [1]. The effect of supersonic blowing on the flow characteristics of the external supersonic stream was studied in [2–4]. The procedure and results of the solution to the problem of subsonic blowing of a homogeneous fluid at the base of a body in supersonic flow are discussed in this paper. Analysis of experimental results (see, e.g., [5]) shows that within a certain range of blowing rate the pressure distribution along the viscous region differs very little from the pressure in the free stream ahead of the base section. In this range the flow in the blown subsonic jet and in the mixing zones can be described approximately by slender channel flow. This approximation is used in the computation of nozzle flows with smooth wall inclination [6, 7]. On the other hand, boundary layer equations are used to compute separated stationary flows with developed recirculation regions [8] in order to describe the flow at the throat of the wake. The presence of blowing has significant effect on the flow structure in the base region. An increasing blowing rate reduces the size of the recirculation region [9] and increases base pressure. This leads to a widening of the flow region at the throat, usually described by boundary-layer approximations. At a certain blowing rate the recirculation region completely disappears which makes it possible to use boundary-layer equations to describe the flow in the entire viscous region in the immediate neighborhood of the base section.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 76–81, January–February, 1984.     

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.     

Use of direct variational methods to optimize axisymmetric Laval nozzles in the case of equilibrium and nonequilibrium two-phase flows

In the construction of the optimal profile of a Laval nozzle when there are subsonic regions in the flow, the use of effective methods such as the general method of Lagrangian multipliers [1] becomes very difficult. In the present paper, direct variational methods are therefore used. For nozzles, these methods were used for the first time to profile the supersonic parts of nozzles in the case of nonequilibrium two-phase flows by Dritov and Tishin [2]. For equilibrium flows, they have been used to optimize supersonic nozzles [3, 4] and in the construction of a profile of the subsonic part of a nozzle ensuring parallel sonic flow in the minimal section of the nozzle [3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 181–183, January–February, 1982.I thank A. N. Kraiko for a number of helpful comments in a discussion of the formulation of the problem.     

Asymptotic analysis of transonic flows

In this paper we derive the equations of the second and third approximations for the stream function of two-dimensional and axisymmetric potential transonic flow of an inviscid gas and find their particular solutions corresponding to certain transonic flows.A similar study concerning the second approximation of subsonic and supersonic flow was made by Van Dyke [1] and Hayes [2]. The second approximation for the velocity potential of transonic flow has been examined in detail by Hayes [3]. Euvrard [4, 5] has investigated the asymptotic behavior of transonic flow far from a body, while Fal'kovich, Chernov, and Gorskii [6] have studied the flow in a nozzle throat.The transonic asymptotic analysis for the stream function is presented in this paper.     

Flow around a sphere by a two-phase supersonic jet

Results are presented of a calculation of the flow around a sphere of a two-phase supersonic jet, discharging into a vacuum. Calculations were performed by the determination method with use of a difference grid constructed on the basis of characteristic ratios [1], The parameters of the unperturbed jet were determined with the two-velocity and two-temperature model of mutually penetrating flows of continuous media (gas and particles) [2, 3] by the network method [4]. In calculating the flow around the sphere, as in [5–7], it was assumed that the particles do not affect the gas flow in the shock layer. An analysis of the effect of particles on gasdynamic parameters in a shock layer was performed in [8].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 171–176, November–December, 1978.The authors are grateful to A. N. Nikulin for providing the program for calculation of flow about a blunt body by a uniform supersonic flow.     

Quasi-one-dimensional approximation in two-dimensional problems of gas dynamics

A useful means of constructing approximate flow models is the hydraulic (for two-dimensional problems quasi-one-dimensional) approach, based on averaging the initial nonuniform flows over some direction or cross section [1]. In this case, at the expense of a rougher model it is possible to reduce the dimensionality of the problem. Here, this approach is extended to unsteady two-dimensional gas-dynamic processes; certain problems (flow around a cone or a blunt body, jet flows) are considered in the framework of the quasi-one-dimensional model obtained, and results are compared with the solutions of the corresponding two-dimensional problems.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 136–143, March–April, 1989.     

Variational problem of optimum side wall design for a narrow three-dimensional nozzle

The optimum design of the side walls of the supersonic section of a three-dimensional nozzle with two planes of symmetry is considered in the narrow channel model approximation, which reduces three-dimensional to two-dimensional flow. This nozzle realizes maximum thrust for given sonic or supersonic inlet flow, upper and lower walls, maximum permissible length and pressure outside the nozzle. In general, an approximate solution of the variational problem can be obtained by the indeterminate control contour method [1]. For nozzles with nonexpanding end sections of the upper and lower walls this is a rigorous solution. Numerical algorithms, based on the method of characteristics, for constructing the optimum, side walls and calculating the flow in narrow channels are developed in the formulation adopted using the optimality conditions found, which generalize the wellknown conditions for plane and axisymmetric configurations [1]. In addition, the three-dimensional supersonic flow in the nozzles thus designed has been calculated in accordance with a shock-capturing marching scheme [2], which for the uniform grids employed in the calculations gives a second-order approximation. A rather complex relation is established between the thrust of the optimum configurations constructed and the shape of their inlet cross sections.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.2, pp. 102–112, March–April, 1992.The authors are grateful to L. E. Sternin for drawing their attention to the problem and to V. A. Vostretsova for assisting with the work.     

Calculation of unsteady supersonic flow past a two-dimensional cascade of plates ween vorticity inhomogeneities of the flow act on it

In the framework of the linear theory of small perturbations the problem of unsteady subsonic flow past a two-dimensional cascade of plates has been considered in a number of papers. Thus, the unsteady aerodynamic characteristics of a cascade of vibrating plates were calculated in [1] by the method of integral equations, while the same method was used in [2, 3] to calculate the sound fields that are excited when sound waves Coming from outside or vorticity inhomogeneities of the oncoming flow act on the cascade. The problem of a two-dimensional cascade of vibrating plates in a supersonic flow was solved in [4, 5]. In [4] the solution was constructed on the basis of the well-known solution of the problem of vibrations of a single plate, while in [5] a variant of the method of integral equations was used which differed slightly from the usual formulation of this method [1–3]. The approach proposed in [5] is used below to calculate the unsteady flow past a two-dimensional cascade of plates in the case when vorticity inhomogeneities of a supersonic oncoming flow act on it. Equations are obtained for the strength of the unsteady pressure jumps arising in such a flow and the vortex wakes shed from the trailing edges of the plates. Examples of the calculations illustrating the accuracy of the method and its possibilities are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp, 152–160, May–June, 1986.     

Slender axisymmetric cavities in a supersonic flow

Slender axisymmetric cavities in a subsonic flow of compressible fluid were investigated in [1–4]. In [5] a finite-difference method was used to calculate the drag coefficient of a circular cone, near which the shape of the cavity was determined for subsonic, transonic, and supersonic water flows; however, in the supersonic case the entire shape of the cavity was not determined. Here, on the basis of slender body theory an integrodifferential equation is obtained for the profile of the cavity in a supersonic flow. The dependence of the cavity elongation on the cavitation number and the Mach number is determined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 179–181, January–February, 1989.     

Numerical study of steady supersonic separated flow past spatial lifting systems

At the present, spatial lifting systems are usually calculated numerically using linear approximation. However, the practical application of such systems at moderate and large angles of incidence requires new approaches that allow for various nonlinear effects such as large disturbances, flow separation, and jumps in entropy across shock waves. The existing investigations [3, 4] generally cover only simple systems (bodies of revolution, wings, and so on). Here, a numerical method is proposed for investigating supersonic flows past complicated spatial systems. The method extends and continues the well-known methods widely used to solve analogous problems in subsonic aerodynamics [5, 6]. Some examples of the computation of the aerodynamic parameters for flows past wings and spatial lifting systems are also given.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 142–148, May–June, 1993.     

Characteristics of an unstably localized disturbance in a compressed boundary layer

Experiments have demonstrated [1] that the transition of streamline-type flow into turbulent flow in a boundary layer occurs as a result of the formation and development of turbulent spots apparently arising from small natural disturbances. A study of the nonlinear evolution and interaction of localized disturbances requires knowledge of their characteristics to a linear approximation [2]. In the current work, results are presented of calculations of such characteristics for the first two unstable modes in a supersonic boundary layer on a two-dimensional plate (M = 4.5, T_w = 4.44).Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 50–53, January–February, 1976.     

Determination of optimum nozzle parameters in a gasdynamic laser

In connection with the use of supersonic nozzles to create lasers, the question arises of the optimum parameters of the nozzle and the gas mixture from the aspect of obtaining the greatest population inversion of the energy levels of internal degrees of freedom of molecules of the working gas and the greatest output power of the lasers. A rather complete concept of the kinetic processes taking place during the escape of a relaxing gas mixture containing carbon dioxide through a supersonic nozzle has now been developed on the basis of calculated and experimental data. In [1–4] the problems of optimization of the parameters of a CO₂-N₂-H₂O-He mixture and of the shape of the nozzle were set up and solved in a one-dimensional steady-state formulation. The influence of the two-dimensionality of the stream in an optimum nozzle on the laser characteristics is studied in the present report. The method of through calculation suggested in [5] is used to calculate the two-dimensional flow of a relaxing gas.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 23–26, September–October, 1978.     

“Lateral” interaction of a supersonic underexpanded ideal-gas jet with surfaces of different shape

A numerical investigation is made of the interaction of an underexpanded jet of an inviscid and nonheat-conducting gas issuing from an axisymmetric conical nozzle with plane, cylindrical, and spherical surfaces. It is assumed that the flow turning angle for flow about a barrier is smaller than the critical angle, and subsonic regions are absent in the flow field studied. The effect of the characteristic parameters (Mach number at the nozzle exit, jet underexpansion) on the flow pattern and jet forces is analyzed. The results of numerical calculations are compared to the results of approximate theories and experimental data. A theoretical solution of the problem of the effect of a supersonic jet on a surface of given shape, even in the approximation of an inviscid, nonheat-conducting gas, is quite difficult. A reason for this is that the flow region contains shock waves interacting with each other, contact discontinuities, and zones of mixed sub-and supersonic flow. As far as is known to the authors, the results obtained for three-dimensional problems for the interaction of supersonic jets with each other or with barriers are primarily experimental (for example, [1–6]). A numerical analysis of the interaction of axisymmetric ideal-gas jets was carried out in [7–10]. In [7] a three-dimensional form of the method of characteristics was used to calculate the initial interaction region for two supersonic cylindrical jets (with Mach number M=10) intersecting at an angle of 60. The interaction of several jets has been considered in [8, 9], where the solution was obtained according to the Lax—Wendroff method without elimination of the discontinuity lines of flow parameters. In [10] the lateral interaction of axisymmetric supersonic jets with each other and with a plate is investigated by means of a straight-through calculationTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 3–8, November–December, 1974.The authors thank A. N. Kraiko for useful discussions of the results, and A. L. Isakov and É. N. Gasparyan for kindly providing the experimental data.     

Nozzle thrust increase by flow rotation

Until recently, in constructing optimum supersonic axisymmetric nozzles, gas motions without rotation about the axis have been considered. This question has been studied rather completely for equilibrium gas flows in [1–6]. Avoidance of flow rotation is a limitation which may reduce the nozzle thrust, other conditions being equal. It is easy to show that for zero nozzle length or for a length that permits obtaining uniform flow at the exist that maximum thrust is provided by uniform flow parallel to the nozzle axis. If the length limitation does not make it possible to do this, then the use of freedom in the gas rotation may increase the nozzle thrust. Using a very simple example it is shown that this possibility is realized.     

Conservation of the resultant force vector in the plane of the angle of attack for slender star-shaped bodies in supersonic flows

At high supersonic flight speeds bodies with a star-shaped transverse and power-law longitudinal contour are optimal from the standpoint of wave drag [1–3]. In most of the subsequent experimental [4–6] and theoretical [6–9] studies only conical star-shaped bodies have been considered. For these bodies in certain flow regimes ascent of the Ferri point has been noted [10]. In [11] the boundary-value problem for elongated star-shaped bodies with a power-law longitudinal contour was solved for the case of supersonic flow. The present paper deals with the flow past these bodies at an angle of attack. It is found that for arbitrary star-shaped bodies with any longitudinal (in particular, conical) profile the aerodynamic forces can be reduced to a wave drag and a lift force, the lateral force on these bodies being equal to zero for any position of the transverse contour.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 135–141, November–December, 1989.     

Acoustic,entropy, and vorticity perturbations in a channel of variable cross section

A study is made of the problem of the propagation of infinitesimally small perturbations in a gas stream moving in a channel of variable cross section when the flow cannot be regarded as isentropic and irrotational. The solution is found in the framework of the linear theory of the flow of an ideal gas and the quasi-one-dimensional hydraulic approximation for the steady regime. For irrotational and isentropic perturbations in a nozzle, this problem was considered in [1–4]. In [1], the problem is generalized to take into account entropy perturbations in the nozzle for the case of longitudinal oscillations. The present paper treats arbitrary modes in a nozzle and takes into account not only entropy but also vorticity perturbations in the moving stream. For each of the three perturbation types — acoustic, entropy, and vorticity — the solutions are expanded in series in cylindrical functions. It is shown that in the considered approximation each oscillation mode can be analyzed independently of the others. In the special case of flow in a Laval nozzle, the concept of impedance (admittance), which is widely used in acoustics, is generalized to take into account entropy and vorticity perturbations. The contribution to the flow dynamics of the acoustic, entropy, and vorticity perturbations is estimated numerically for longitudinal and transverse modes.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 91–98, January–February, 1982.