Investigation of hypersonic three-dimensional flow past blunt bodies within the framework of the parabolized Navier-Stokes equations

The flow of a homogeneous gas in a three-dimensional hypersonic viscous shock layer, which includes the shock wave structure, is examined within the framework of the parabolic approximation of the Navier-Stokes equations. The Navier-Stokes equations are simplified on the basis of the asymptotic analysis carried out in [1], are written in variables of the Dorodnitsyn type [2] and are solved by the method proposed in [3, 4] extended to the case of three-dimensional flows. The flow at zero angle of attack past elliptic paraboloids, two-sheeted hyperboloids and triaxial ellipsoids is calculated. Some results of investigating the flow past such bodies are presented. Flow past a sphere in the analogous approximation was considered in [5], where a comparison was also made with the solution of the complete Navier-Stokes equations [6, 7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 134–142, July–August, 1987.In conclusion, the authors wish to express their warm thanks to V. V. Lunev and G. A. Tirskii for useful discussion and valuable comments.     

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.     

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.     

Similarity of three-dimensional and axisymmetric chemically-nonequilibrium flows in the neighborhood of a plane of symmetry

Three-dimensional chemically-nonequilibrium flow past blunt bodies in the neighborhood of the plane of symmetry is investigated within the framework of viscous shock layer theory. The similarity of three-dimensional and axisymmetric flows, previously established in [1] for a uniform gas, is extended to chemically-nonequilibrium gas flows. It is shown that the problem of determining the heat fluxes and friction stress in the neighborhood of the line of flow divergence can be reduced to the problem of determining these quantities for the axisymmetric body. The validity of the axisymmetric analogy is verified by carrying out numerical calculations for bodies of different shapes re-entering the earth's atmosphere along a gliding trajectory. Various models of surface catalytic activity are considered. The use of similarity relations makes it possible to apply existing programs for calculating axisymmetric flows to the solution of three-dimensional problems.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 115–120, March–April, 1990.     

Separated flow over a strongly cambered conical wing at small angles of attack

The interference of supersonic flows on the concave surface of conical wings was experimentally investigated in [1] for various values of the camber and angles of attack. In order to establish the detailed structure of the interference flow the laminar flow past a wing model in the form of half the surface of a circular cone with vertex angle 2_k = 34° was numerically modeled within the framework of the quasiconical approximation for the three-dimensional Navier-Stokes equations [2]. Under this assumption, confirmed by analysis of the experimental data [1], it was found that the displacement of the external inviscid flow as a result of intense flow separation beyond the leading edges leads to flow patterns similar to those realized on V wing's with a break in the transverse contour [3]. At nonzero angles of attack weak secondary separation was detected beneath the flattened regions of primary separation located in the shaded parts of the concave surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 130–136, July–August, 1989.     

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.     

Supersonic flow over elongated blunt bodies with a cross section of elliptical shape

Up to now computational algorithms have been developed for, and systematic studies have been made of, supersonic flow over axisymmetric bodies both by a stream of ideal gas and by an air stream with equilibrium and nonequilibrium physicochemical transformations [1–6]. Conical flows around bodies having cross sections of different shapes and in a wide range of angles of attack have been studied in detail [7–11]. With the further development of numerical methods the next problem has become the analysis of supersonic flow over blunt bodies of large elongation having cross sections of sufficiently arbitrary shape. The effects of essentially three-dimensional flow (without planes of symmetry) over bodies whose cross sections represent ellipses with a constant or variable ratio of axes along the length of the body are discussed in the present paper.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6. pp. 155–159, November–December, 1976.     

Numerical modeling of three-dimensional viscous flows using a space-marching method with pressure iterations

A method of calculating steady three-dimensional compressible and incompressible flows on the basis of simplified Navier-Stokes equations is examined. Some calculation results are presented. These relate to the problem of supersonic flow past blunt bodies and to incompressible channel flows.Based on paper read at the Seventh Congress on Theoretical and Applied Mechanics, Moscow, August I991. Presented by V. I. Polezhaew.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 132–147, September–October, 1992.     

Unsteady three-dimensional viscous shock layer in nonuniform gas flow in the neighborhood of a stagnation point

The combined influence of unsteady effects and free-stream nonuniformity on the variation of the flow structure near the stagnation line and the mechanical and thermal surface loads is investigated within the framework of the thin viscous shock layer model with reference to the example of the motion of blunt bodies at constant velocity through a plane temperature inhomogeneity. The dependence of the friction and heat transfer coefficients on the Reynolds number, the shape of the body and the parameters of the temperature inhomogeneity is analyzed. A number of properties of the flow are established on the basis of numerical solutions obtained over a broad range of variation of the governing parameters. By comparing the solutions obtained in the exact formulation with the calculations made in the quasisteady approximation the region of applicability of the latter is determined. In a number of cases of the motion of a body at supersonic speed in nonuniform media it is necessary to take into account the effect of the nonstationarity of the problem on the flow parameters. In particular, as the results of experiments [1] show, at Strouhal numbers of the order of unity the unsteady effects are important in the problem of the motion of a body through a temperature inhomogeneity. In a number of studies the nonstationary effect associated with supersonic motion in nonuniform media has already been investigated theoretically. In [2] the Euler equations were used, while in [3–5] the equations of a viscous shock layer were used; moreover, whereas in [3–4] the solution was limited to the neighborhood of the stagnation line, in [5] it was obtained for the entire forward surface of a sphere. The effect of free-stream nonuniformity on the structure of the viscous shock layer in steady flow past axisymmetric bodies was studied in [6, 7] and for certain particular cases of three-dimensional flow in [8–11].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 175–180, May–June, 1990.     

Multicomponent three-dimensional viscous shock layer on blunt bodies with a catalytic surface at angles of attack and yaw

The three-dimensional supersonic flow of nonequilibrium dissociating air past smooth blunt bodies on whose surface heterogeneous chemical reactions are taking place is investigated within the framework of thin viscous shock layer theory. An economical numerical method of solving the equations with an improved order of approximation with respect to the normal coordinate is employed. This method does not require the preliminary solution of the Stefan-Maxwell relations for the diffusion fluxes and makes it possible to calculate flows that do not possess a plane of symmetry. The effect of the angles of attack and yaw, the catalytic reaction model and a number of other parameters of the problem on the pressure, heat flux and equilibrium surface temperature distributions is analyzed with reference to the example of flow past a triaxial ellipsoid.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 143–150, January–February, 1990.Tha authors are grateful to G. A. Tirskii for useful discussions of their results.     

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.     

Quasi-three-dimensional calculation of flow in blade passages of turbines

The flow in turbomachines is currently calculated either on the basis of a single successive solution of an axisymmetric problem (see, for example, [1-A]) and the problem of flow past cascades of blades in a layer of variable thickness [1, 5], or by solution of a quasi-three-dimensional problem [6–8], or on the basis of three-dimensional models of the motion [9–11]. In this paper, we derive equations of a three-dimensional model of the flow of an ideal incompressible fluid for an arbitrary curvilinear system of coordinates based on averaging the equations of motion in the Gromek–Lamb form in the azimuthal direction; the pulsation terms are taken into account in the equations of the quasi-three-dimensional motion. An algorithm for numerical solution of the problem is described. The results of calculations are given and compared with experimental data for flows in the blade passages of an axial pump and a rotating-blade turbine. The obtained results are analyzed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 69–76, March–April, 1991.I thank A. I. Kuzin and A. V. Gol'din for supplying the results of the experimental investigations.     

Investigation of flow past blunt bodies with allowance for radiation by the large-particle method

Much recent work has been done on developing methods of solving gas-dynamic problems in which radiation plays a part (see, for example, [1–7]). This is because the temperature in the shock layer associated with flight in the atmosphere at hypersonic velocities can reach values exceeding 10⁴ °K. In such a case, heat transfer by radiation can make an important contribution to the total heat transfer. With increasing flight velocities, the importance of radiation in heat transfer increases and then becomes predominant. In the present paper, the large-particle method as developed by Belotserkovskii and Davydov [8] is developed to calculate flows with radiation around blunt bodies, including the case when there is distributed blowing from the surface of the bodies into the shock layer, which simulates ablation of a heat-shielding covering under the influence of strong heating by radiation. The results are given of systematic calculations of flow past blunt bodies of various shapes by a stream of radiating air, and the results are compared with the data of other methods.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 106–112, July–August, 1982.     

Three-dimensional laminar boundary layer on a permeable surface in the neighborhood of a symmetry plane

Analytical and numerical methods are used to investigate a three-dimensional laminar boundary layer near symmetry planes of blunt bodies in supersonic gas flows. In the first approximation of an integral method of successive approximation an analytic solution to the problem is obtained that is valid for an impermeable surface, for small values of the blowing parameter, and arbitrary values of the suction parameter. An asymptotic solution is obtained for large values of the blowing or suction parameters in the case when the velocity vector of the blown gas makes an acute angle with the velocity vector of the external flow on the surface of the body. Some results are given of the numerical solution of the problem for bodies of different shapes and a wide range of angles of attack and blowing and suction parameters. The analytic and numerical solutions are compared and the region of applicability of the analytic expressions is estimated. On the basis of the solutions obtained in the present work and that of other authors, a formula is proposed for calculating the heat fluxes to a perfectly catalytic surface at a symmetry plane of blunt bodies in a supersonic flow of dissociated and ionized air at different angles of attack. Flow near symmetry planes on an impermeable surface or for weak blowing was considered earlier in the framework of the theory of a laminar boundary layer in [1–4]. An asymptotic solution to the equations of a three-dimensional boundary layer in the case of strong normal blowing or suction is given in [5, 6].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 37–48, September–October, 1980.     

Analytic investigation of friction and heat transfer in the neighborhood of a three-dimensional stagnation point at low and moderate Reynolds numbers

A study is made of hypersonic three-dimensional flow of a viscous gas past blunt bodies at low and moderate Reynolds numbers with allowance for the effects of slip and a jump of the temperature across the surface. The equations of the three-dimensional viscous shock layer are solved by an integral method of successive approximation and a finite-difference method in the neighborhood of the stagnation point. In the first approximation of the method an analytic solution to the problem is found. Analysis of the obtained solution leads to the proposal of a simple formula by means of which the calculation of the heat flux to a three-dimensional stagnation point is reduced to the calculation of the heat flux to an axisymmetric stagnation point. A formula for the relative heat flux obtained by generalizing Cheng's well-known formula [1] is given. The accuracy and range of applicability of the obtained expressions are estimated by comparing the analytic and numerical solutions. Three-dimensional problems of the theory of a supersonic viscous shock layer at small Reynolds numbers were considered earlier in [2–5] in a similar formulation but without allowance for the effects of slip.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 143–150, March–April, 1988.     

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.     

Interaction of two supersonic radial gas flows

The flow resulting from the collision of two spherically symmetric supersonic gas streams generated by three-dimensional sources has been studied within the framework of the inviscid perfect gas model. When the characteristics of both sources are the same, the problem reduces to that of the interaction between a spherically symmetric flow and a plane barrier [1, 2]. By means of a certain passage to the limit, the flow from one of the sources can be reduced to a uniform supersonic gas flow. In this case the problem reduces to the problem of uniform gas flow past a source considered in [1, 3, 4]. The resulting flows are investigated with reference to all the parameters characterizing the two sources. The shock wave structure, both in the neighborhood of and at points remote from the axis of symmetry, is studied, together with the distributions of the gas dynamic quantities in the shock layers; certain similarity laws are established. The astrophysical applications of the problem associated, in particular, with a certain x-ray radiation mechanism in binary stellar systems [5] are of particular interest.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 159–165, July–August, 1990.The authors are grateful to V. B. Baranov for his constant interest and to S. A. Zhekov and V. V. Usov for useful discussions.     

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.     

Numerical investigation of the structure of nonequilibrium three-dimensional hypersonic flow past blunt bodies

Three-dimensional dissociating air flow past blunt bodies is investigated within the framework of the parabolized Navier-Stokes equations in the thin layer approximation. Multicomponent diffusion, barodiffusion and homogeneous chemical reactions, including dissociation-recombination and exchange reactions, are taken into account. The boundary conditions are assigned in the free stream and at the surface of the body with allowance for heterogeneous catalytic reactions and slip effects. The problem of flow at zero angle of attack past blunt bodies possessing two planes of symmetry is investigated numerically for flow patterns varying from smeared layer structure to almost ideal flow (Re=50-10⁵). The flow conditions corresponded to the motion of a body with lift along a re-entry trajectory [1]. The contribution of the chemical reactions in the shock wave as compared to the diffusion flux at the edge of the shock wave was estimated. The edge of the shock wave is assumed to correspond to the point at which the density profile has the greatest slope. The influence of slip effects and barodiffusion on the flow characteristics is demonstrated. The results of the calculations are compared with calculations based on the thin viscous shock layer model [2].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 143–150, September–October, 1987.The author wishes to express his thanks to G. A. Tirskii and V. V. Lunev for useful discussions and valuable advice.     

Incompressible fluid flow past a circular cylinder with linear relation between the vorticity and the stream function

Incompressible fluid flow with a linear relationship between the vorticity and the stream function past a circular cylinder is studied.Vortical flows about profiles have been considered in several studies [1–15], but in all these studies with the exception of [15] a constant vorticity was assumed (in [15] an approximate solution is found of the problem of incompressible fluid flow about a Zhukovskii profile with parabolic distribution of the velocities in the approaching stream).A freestream velocity profile similar to that considered below occurs, for example, in a planar jet (laminar or turbulent), in the wake behind a bluff body, in the boundary layer along an infinite plane [4,13], in turbulent jet flows with reverse fluid currents [16]. A similar situation also arises in the flow past an array of cylinders with large spacing which is located in the wake of another array.The author wishes to thank V. E. Davidson for posing the problem and for guidance in its solution.