Hypersonic flow of nonequilibrium air over blunt bodies

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

Investigation of supersonic three-dimensional flow past pointed axisymmetric bodies

The flow pattern near bodies of revolution with very long cylindrical and pointed nose sections is studied in the framework of an ideal gas model by means of a numerical method based on MacCormack's difference scheme. The existence of internal shock waves, oriented in both the longitudinal and the transverse directions, in the shock layer is established. The variation of the aerodynamic coefficients of the configuration with its length, angle of attack, and free stream Mach number is investigated. The calculated and experimental data are compared, and the connection between the flow parameters on the body surface and the position of the separation line of the boundary layer on its lateral face is established. A method of calculating the influence of the boundary layer on the values of the aerodynamic coefficients of bodies of revolution of large aspect ratio at small angles of attack is proposed. Axisymmetric flow near blunt bodies has been studied in detail in [1].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 127–133, September–October, 1986.The author expresses his gratitude to A. N. Pokrovskii for his help in calculating the boundary layer parameters on the surfaces of the considered configurations.     

Nonequilibrium hypersonic flow of gas past a wing of small aspect ratio

It is well known [1] that nonequilibrium physicochemical processes taking place in gases at high temperature influence the gas-dynamic parameters and aerodynamic characteristics of bodies in hypersonic flight. In the present paper, the thin shock layer method [2–4] is used to consider the problem of nonequilibrium hypersonic flow of gas past a wing of small aspect ratio at an angle of attack. It is shown that the flow component of the vorticity is conserved along the streamlines, and this property is exploited to obtain an analytic solution of the equations of the three-dimensional nonequilibrium shock layer. The influence of the disequilibrium on the thickness of the shock layer and the pressure distribution is investigated.     

Determination of the magnitude of the separation criterion for a three-dimensional laminar boundary layer

A separation criterion, i.e., a definite relationship between the external flow and the boundary layer parameters [1], can be used to estimate the possibility of the origination of separation of a two-dimensional boundary layer. A functional form of the separation criterion has also been obtained for a three-dimensional boundary layer [2] on the basis of dimensional analysis. As in the case of the two-dimensional boundary layer, locally self-similar solutions can be used to determine the specific magnitude of the separation criterion as a function of the values of the governing parameters. Locally self-similar solutions of the two-dimensional laminar boundary-layer equations have been found at the separation point for a perfect gas with a linear dependence of the coefficient of viscosity on the temperature (Ω=1) and Prandtl number P=1 [3, 4]. The influence of blowing and suction has been studied for this case [5]. Self-similar solutions have been obtained for Ω=1, P=0.723 for the limit case of hypersonic perfect gas flow [6]. Locally self-similar solutions of the three-dimensional laminar boundary-layer equations at the separation point are presented in [7] for a perfect gas with Ω=1, P=1. There are no such computations for Ω≠1, P≠1; however, the results of computing several examples for a two-dimensional flow [8] show that the influence of the real properties of a gas can be significant and should be taken into account. Self-similar solutions of the two- and three-dimensional boundary-layer equations at the separation point are found in this paper for a perfect gas with a power-law dependence of the viscosity coefficient on the enthalpy (Ω=0.5, 0.75, 1.0) for different values of the Prandtl number (P=0.5, 0.7, 1.0) in a broad range of variation of the external stream velocity (v ₁ ² /2h₁* = 0–0.99) and the temperature of the streamlined surface. Magnitudes of the separation criterion for a laminar boundary layer have been obtained on the basis of these data.     

Plate cascade in subsonic unsteady gas flow

Accounting for fluid compressibility creates serious difficulties in solving the problem of oscillations of a grid of thin, slightly curved profiles in a subsonic stream. The problem has been solved in [1–3] for a widely-spaced cascade without stagger whose profiles oscillate in phase opposition. The phenomenon of aerodynamic (acoustic) resonance, which may arise in a grid in the direction transverse to the stream for definite values of the stream velocity and profile oscillation frequency, was discovered in [2]. An approximate solution of the problem in which account is not taken of the effect of the vortex trails on the gas flow has been obtained in [4]. In [5, 6] Meister studied in the exact linear formulation the problem of unsteady gas motion through an unstaggered cascade of semi-infinite plates. In [7] Meister considered a grid of profiles with finite chords, but the problem solution was not completed. The problem of subsonic gas flow through a staggered lattice whose profiles oscillate following a single law with constant phase shift was solved most completely in the studies of Kurzin [8, 9] using the method of integral equations. A method of solving the problem for the case of arbitrary harmonic oscillation laws for the lattice profiles was indicated in [10]. The results of the calculation of the unsteady aerodynamic forces for the particular case of a plate cascade without stagger are presented in [9,11], and the possibility of the occurrence of aerodynamic resonance in the cascade in the directions transverse to and along the stream is indicated.Another method of solving the problem is given in [12], in which the more general problem of unsteady subsonic gas flow through a three-dimensional cascade of plates is solved. In the present study this method is applied to the solution of the problem of oscillations of staggered plate cascades in a two-dimensional subsonic gas flow. The results are presented of an electronic computer calculation of the unsteady aerodynamic characteristics of the cascade profiles, which show the essential influence of fluid compressibility on these characteristics. In particular, a sharp decrease of the aerodynamic damping in the acoustic resonance regimes is obtained.     

Effect of pressure gradient of friction and heat transfer in a dusty boundary layer

Approximate analytic expressions for the local friction and heat transfer coefficients in a dusty laminar boundary layer are obtained and tested in the case of an incompressible carrier phase, power-law variation of the external gas flow velocity and small velocity and temperature phase disequilibrium. These expressions supplement the numerical analysis of the dusty boundary layer on a blunt body [1, 2] and the asymptotic calculation of the friction and heat transfer in a quasiequilibrium dusty gas boundary layer on a plate [3]. The combined effect of dustiness and pressure gradient on the friction and heat transfer coefficients is discussed. The results obtained can be used for the practical calculation of the friction and heat transfer in a quasiequilibrium dusty laminar boundary layer and for interpreting the corresponding experimental data. Tomsk. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 105–108, September–October, 1988.     

Oscillations of a plate cascade in a transonic gas flow

The transonic unsteady flow of a gas through a cascade of thin, slightly curved plates is quite complex and has received little study. The main difficulties are associated with the nonlinear dependence of the aerodynamic characteristics on the plate thickness. In [1] it is shown that, for a single thin plate performing high-frequency oscillations in a transonic gas stream, the variation of the unsteady aerodynamic characteristics with plate thickness may be neglected. For a plate cascade, the flow pattern is complicated by the aerodynamic interference between the plates, which may depend significantly on their shape. A solution of the problem of transonic flow past a cascade without account for the plate thickness has been obtained by Hamamoto [2].The objective of the present study is the clarification of the dependence of the aerodynamic characteristics of a plate cascade on plate thickness in transonic unsteady flow regimes. The nonlinear equation for the velocity potential is linearized under the assumption that the motionless plate causes significantly greater disturbances in the stream than those due to the oscillations. A similar linearization was carried out for a single plate in [3]. The aerodynamic interference between the plates is determined by the method presented in [4]. As an example, the aerodynamic forces acting on a plate oscillating in a duct and in a free jet are calculated.     

Laminar boundary layer in flow of a nonequilibrium ionized gas

The motion of a hypersonic body is accompanied by an increase in the gas temperature in the boundary layer up to tens of thousands of degrees, which causes the gas to ionize. Under these conditions there are problems in calculating coefficients of viscosity, diffusion, and heat conduction. Investigations have shown that it is invalid to extrapolate the widely used approximations for transport coefficients in the high temperature region [1–3]. This paper considers the laminar boundary layer in the vicinity of the stagnation point of a blunt body in a stream of monatomic nonequilibrium ionized gas. The main thrust is a more accurate calculation of transport coefficients and an investigation of their effect on profiles of the gasdynamic parameters. A specific calculation is performed for argon by way of example.     

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.     

Method of mean-mass parameters for a three-dimensional boundary layer in a flow with vorticity

When a gas flows with hypersonic velocity over a slender blunt body, the bow shock induces large entropy gradients and vorticity near the wall in the disturbed flow region (in the high-entropy layer) [1]. The boundary layer on the body develops in an essentially inhomogeneous inviscid flow, so that it is necessary to take into account the difference between the values of the gas parameters on the outer edge of the boundary layer and their values on the wall in the inviscid flow. This vortex interaction is usually accompanied by a growth in the frictional stress and heat flux at the wall [2, 3]. In three-dimensional flows in which the spreading of the gas on the windward sections of the body causes the high-entropy layer to become narrower, the vortex interaction can be expected to be particularly important. The first investigations in this direction [4–6] studied the attachment lines of a three-dimensional boundary layer. The method proposed in the present paper for calculating the heat transfer generalizes the approach realized in [5] for the attachment lines and makes it possible to take into account this effect on the complete surface of a blunt body for three-dimensional laminar, transition, or turbulent flow regime in the boundary layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 80–87, January–February, 1981.     

A spatial laminar boundary layer on a streamline in conical external flow of a uniform gas in the absence of sources and sinks

Theoretical study of a three-dimensional laminar boundary layer is a complex problem, but it can be substantially simplified in certain particular cases and even reduced to the solution of ordinary differential equations.One such particular case is the flow of a compressible gas on a streamline in conical external flow. The case is of considerable practical importance because the local heat fluxes may take extremal values on such lines.Such flow, except for the conical case, has been examined [1–4], and an approximate method has been given [1] on the basis of integral relationships and a special form for the approximating functions. A numerical solution has been given [2, 3] for such flow around an infinite cylinder. It was assumed in [1–3] that the Prandtl number and the specific heats were constant, and that the dynamic viscosity was proportional to temperature. Heat transfer has been examined [4] near a cylinder exposed to a flow of dissociated air.Here we give results from numerical solution of a system of ordinary differential equations for the flow of a compressible gas in a laminar boundary layer on streamlines in conical external flow, with or without influx or withdrawal of a homogeneous gas. It is assumed that the gas is perfect and that the dynamic viscosity has a power-law temperature dependence.     

Effects of second-approximation of boundary layer theory in three-dimensional flow over bodies with a large aspect ratio at small angles of incidence

The influence of the second-approximation effects of laminar boundary layer theory on the heat transfer in three-dimensional hypersonic flow over blunt cones with a large aspect ratio is investigated numerically.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 57–64, March–April, 1995.     

Two-phase couette flow

A study is made of plane laminar Couette flow, in which foreign particles are injected through the upper boundary. The effect of the particles on friction and heat transfer is analyzed on the basis of the equations of two-fluid theory. A two-phase boundary layer on a plate has been considered in [1, 2] with the effect of the particles on the gas flow field neglected. A solution has been obtained in [3] for a laminar boundary layer on a plate with allowance for the dynamic and thermal effects of the particles on the gas parameters. There are also solutions for the case of the impulsive motion of a plate in a two-phase medium [4–6], and local rotation of the particles is taken into account in [5, 6]. The simplest model accounting for the effect of the particles on friction and heat transfer for the general case, when the particles are not in equilibrium with the gas at the outer edge of the boundary layer, is Couette flow. This type of flow with particle injection and a fixed surface has been considered in [7] under the assumptions of constant gas viscosity and the simplest drag and heat-transfer law. A solution for an accelerated Couette flow without particle injection and with a wall has been obtained in [6]. In the present paper fairly general assumptions are used to obtain a numerical solution of the problem of two-phase Couette flow with particle injection, and simple formulas useful for estimating the effect of the particles on friction and heat transfer are also obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–46, May–June, 1976.     

Lattice of plates in an unsteady supersonic flow

The supersonic unsteady flow of a gas past a lattice of thin, slightly curved profiles has been investigated in several studies. The paper [1] is devoted to an evaluation of the effect of wind tunnel walls on the unsteady aerodynamic characteristics of a profile, and [2] investigates the effects of the boundaries of a free jet. These cases are equivalent respectively to the anti-phase and in-phase oscillations of the profiles of an unstaggered grid. In [3] consideration is given to a more general case of gas flow past a profile in a wind tunnel with perforated walls. Flow past a lattice of profiles with stagger is studied in [4], where the magnitude of the stagger angle is limited by the condition that the lattice leading edge is located in the undisturbed stream.In the present paper we present a method of calculation of the unsteady supersonic flow of a gas past a lattice of profiles with arbitrary stagger. As an example the results are presented of the calculation of the aerodynamic forces and moments acting on an oscillating profile in a wind tunnel with solid walls and in a free jet.     

Strong viscous interaction on slender three-dimensional bodies

Hypersonic flow of a viscous gas past axisymmetric power-law bodies in the regime of strong interaction of the laminar boundary layer with the outer inviscid flow was studied in [1, 2]. In this paper the results obtained in [1, 2] are extended to the case of flow past a slender three-dimentional body.     

3D Steady and Unsteady Bifurcations in a Shock-wave/Laminar Boundary Layer Interaction: A Numerical Study

The principal objective of this paper is to study some unsteady characteristics of an interaction between an incident oblique shock wave impinging a laminar boundary layer developing on a plate plane. More precisely, this paper shows that some unsteadiness, in particular the low frequency unsteadiness, originate in a supercritical Hopf bifurcation related to the dynamics of the separated boundary layer and not necessarily to the coherent structures resulting from the turbulent character of the boundary layer crossing the shock wave. Numerical computations of a shock-wave/laminar boundary-layer interaction (SWBLI) have been compared with a classical test case (Degrez test case) and both two-dimensional and three-dimensional (3D) unsteady Navier–Stokes equations are numerically solved with an implicit dual time stepping for the temporal algorithm and high order AUSM+ scheme for the spatial discretization. A parametric study on the oblique shock-wave angle has been performed to characterize the unsteady behaviour onset. Finally, discussions and assumptions are made about the origin of the 3D low frequency unsteadiness.     

Boundary-layer separation on conical bodies

Some characteristics of the variation in the linear dimensions of the flow separation zones on conical bodies with expanding conical skirts and of variation of the pressure within these zones as a function of variation of the Mach number, Reynolds number, and intensity of the disturbance that causes the boundary layer separation are examined. Experiments were conducted in laminar, transitional, and turbulent flows in flow separation regions. The interaction of viscous and nearly inviscid flows is quite common. This phenomenon occurs in flow past a concave corner, when a compression shock impinges on a boundary layer, and in many other cases. The characteristics of this phenomenon in flow about two-dimensional bodies have been investigated experimentally in [1, 2] and other studies. Attempts have been made to analyze the interaction of compression shocks with the boundary layer theoretically. In “free” separated flows, when the points of separation and reattachment of the boundary layer are not fixed (for example, on a flat plate with a long wedge attached to it), theoretical studies are usually made within the framework of the boundary layer theory with use of the approximate integral methods [3, 4]. In this article we examine some results from studies of free separated flows on conical bodies with conical skirts in laminar, transitional, and turbulent flows (Fig. 1).     

Aerodynamic and thermal characteristics of three-dimensional star-shaped bodies in a rarefied gas

The search for shapes of a body optimal with respect to its aerodynamic characteristics has led to the investigation of star-shaped bodies. Numerous theoretical and experimental studies [1–3] have demonstrated the advantages of such bodies compared with the equivalent bodies of revolution for motion in a dense gas in a wide range of supersonic velocities. Since the folded surface of star-shaped bodies appreciably exceeds in area the surface of the compared bodies of revolution, and in a rarefied gas the forces of viscous friction are of the same order as the pressure forces, it is important to investigate whether star-shaped bodies are optimal for flight at hypersonic speeds in rarefied gases at different altitudes. In the present paper, under the assumptions of the “locality hypothesis” [4] expressions are obtained for the main aerodynamic and thermal characteristics of bodies having a front end in the shape of a regular star with n rays joined smoothly to a circular rear end by the method proposed by Gusarov et al. [3]. A detailed analysis of the main characteristics is made and the range of variation of the aspect ratio and the number of rays of the optimal star-shaped configuration are established for different degrees of rarefaction of the atmosphere.     

Boundary layer on a blunt body in a flow of dusty gas

The problem of interaction of gas-dust flows with solid surfaces arose in connection with the study of the motion of aircraft in a dusty atmosphere [1–2], the motion of a gas suspension in power generators, and in a number of other applications [3]. The presence of a disperse admixture may lead to a significant increase in the heat fluxes [4] and to erosion of the surface [5]. These phenomena are due to the joint influence of several factors — the change in the structure of the carrier-phase boundary layer due to the presence of the particles, collisions of the particles with the surface, roughness of the ablating surface, and so forth. This paper continues an investigation begun earlier [6–7] into the influence of particles on the structure of the dynamical and thermal two-phase boundary layer formed around a blunt body in a flow. The model of the dusty gas [8] has an incompressible carrier phase. The method of matched asymptotic expansions [9] is used to obtain the equations of the two-phase boundary layer. In the frame-work of the refined classification made by Stulov [6], it is shown that the form of the boundary layer equations is different in the presence and absence of inertial precipitation of the particles. The equations are solved numerically in the neighborhood of the stagnation point of the blunt body. The temperature and phase velocity distributions in the boundary layer, and also the friction coefficients and the heat transfer of the carrier phase are found for a wide range of the determining parameters. In the case of an admixture of low-inertia particles that are not precipitated on the body, it is shown that even when the mass concentration of the particles in the undisturbed flow is small their accumulation in the boundary layer can lead to a sharp increase in the thermal fluxes at the stagnation point.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 99–107, September–October, 1985.I thank V. P. Strulov for a discussion.