Outflow of an unexpanded jet into counter supersonic and subsonic flows

The article gives the results of an experimental investigation of the geometric structure of an opposing unexpanded jet. It discusses flow conditions with interaction between the jet and sub- and supersonic flows. It is shown that, with the outflow of an unexpanded jet counter to a supersonic flow, there are unstable flow conditions. For stable flow conditions with one roll, dependences are proposed determining the form of a jet in a supersonic opposing flow. A generalized dependence is obtained for the distribution of the pressure at the surface of a body with a jet, flowing out counter to a subsonic flow. The range of change in the determining parameters are the following: Mach numbers at outlet cross section of nozzle, M _a = 1 and 3; Mach numbers of opposing flow, M = 0.6–0.9 and 2.9; degree of effectiveness of jet, n = p _a /p = 0.5–800 (p _a and p are the static pressures at the outlet cross section of the nozzle and in the opposing flow); the ratios of the specific heat capacities, _a = = 1.4; the drag temperatures of the jet and the flow, T_o = T_oa = 290°K.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 89–96, January–February, 1977.     

Penetration depth of a jet injected into an oncoming supersonic flow

The geometrical characteristics of jets injected through an opening in a flat plate into an oncoming supersonic flow have been studied on a number of occasions [1, 3]. However, the results were analyzed under different suppositions about the important dimensionless parameters. In [1], the degree of underexpansion of the jet, characterized by n = p _a /p, was regarded as decisive; in [3], the experimental points were plotted against the relative dynamic head _a u² _a /(u² ) of the jet. In the present paper, dimensional considerations are used to determine the dimensionless parameters which influence the flow field when an injected jet interacts with an oncoming supersonic gas flow. The influence of these determining dimensionless parameters on the depth of penetration of injected jets into a flow was investigated experimentally. It is shown that the relative depth of penetration is determined basically by the relative specific impulse of the jet, the injection angle, and the shape of the blowing nozzle section.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 151–154, July–August, 1980.     

Supersonic air flow past a sphere with account for vibrational relaxation

A numerical solution is obtained for the problem of air flow past a sphere under conditions when nonequilibrium excitation of the vibrational degrees of freedom of the molecular components takes place in the shock layer. The problem is solved using the method of [1]. In calculating the relaxation rates account was taken of two processes: 1) transition of the molecular translational energy into vibrational energy during collision; 2) exchange of vibrational energy between the air components. Expressions for the relaxation rates were computed in [2]. The solution indicates that in the state far from equilibrium a relaxation layer is formed near the sphere surface. A comparison is made of the calculated values of the shock standoff with the experimental data of [3].Notation uV_max, vV_max velocity components normal and tangential to the sphere surface - V_max maximal velocity - P V _max ² pressure - density - TT temperature - e_viRT vibrational energy of the i-th component per mole (i=–O₂, N₂) - =rb–1 shock wave shape - a _f the frozen speed of sound - HRT/m gas total enthalpy     

Laminar boundary layer on blunt bodies with account for vorticity of the external flow

This paper presents the technique for and results from numerical calculations of the hypersonic laminar boundary layer on blunted cones with account for the vorticity of the external flow caused by the curved bow shock wave. It is assumed that the air in the boundary layer is in the equilibrium dissociated state, but the Prandtl number is assumed constant, =0.72. The calculations were made in the range of velocities 3–8 km/sec, cone half-angles _k=0°–20°. With account for the vortical interaction of the boundary layer with the external flow, the distribution of the thermal flux and friction will depend on the freestream Reynolds number (other conditions being the same). In the calculations the Reynolds number R, calculated from the freestream parameters and the radius of the spherical blunting, varies from 2.5·10³ to 5.10⁴. For the smaller Reynolds numbers the boundary layer thickness on the blunting becomes comparable with the shock standoff, and for R<2.5·10³ it is apparent that we must reconsider the calculation scheme. With R>5·10⁴ for cones which are not very long the vortical interaction becomes relatively unimportant. The results of the calculations are processed in accordance with the similarity criteria for hypersonic viscous gas flow past slender blunted cones [1, 2].     

Flow around stepped bodies with strong compression in a shock layer

This article discusses plane and axisymmetric flows of a nonviscous ideal gas around bodies of stepped form, forming with a Mach number M= and an adiabatic indexN1. The greatest amount of attention is paid to the case where there is no Newtonian free layer, but the shock layer is detached at great distances from the nose of the body.Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 4, pp. 104–112, July–August, 1973.     

Comparison of Calculated and Experimental Data on Supersonic Flow past a Circular Cylinder

The supersonic perfect-gas flow past a circular cylinder is studied on the basis of a numerical analysis of the time-dependent two-dimensional Reynolds equations using a differential q– turbulence model with reference to the experimental conditions. The calculations are carried out at Reynolds and Mach numbers Re=2× 10⁵ and M=1.1, 1.3, and 1.7 and the experimental investigations at Re=1.62×10⁵–2×10⁵ and Mach numbers on the interval 0.7 M 1.7. The calculated and experimental data on the pressure coefficient distribution over the cylinder surface, the location of the separation point on the surface, and the pressure drag coefficient are compared.     

Effect of longitudinal riblets on the aerodynamic characteristics of a straight wing

The results of balance aerodynamic tests on model straight wings with smooth and ribbed surfaces at an angle of attack =–4°–12°, Mach number M=0.15–0.63, and Reynolds number Re=2.4·10⁶–3.5·10⁶ are discussed. The nondimensional riblet spacings ⁺, which determines the effect of the riblets on the turbulent friction drag, and the effect of riblets on the upper and/or lower surface of a straight wing on its drag, lift, and moment characteristics are estimated.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 33–38, March–April, 1995.     

Experimental investigation of the flow in the mixing layer in the initial section of an underexpanded jet in a parallel stream

The results are given of an experimental investigation of the flow in the initial section of a turbulent underexpanded jet exhausting from a profiled nozzle with Mach number M _a = 2.56 at the exit into a parallel stream with Mach number M = 3.1. Analysis of the results of measurement of the fields of the total head p₀ and the stagnation temperature T₀ in conjunction with results of calculation of a jet of an ideal gas make it possible to construct the velocity profile in the mixing layer of the underexpanded jet in the parallel supersonic flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 161–163, January–February, 1981.     

Experimental investigation of separation of turbulent boundary layer in vicinity of a step

The results of an experimental investigation of the separation of a turbulent boundary layer in the vicinity of a step on a flat plate at M = 2 and 3, and Re = U/v = (26–66)·10⁶ m^–1 are given. The step height was varied from 3 to 16 mm, which corresponded to the range of relative heights 1.1 h/ 7.6, where is the thickness of the boundary layer at the point at which the pressure starts to increase in front of the step. The obtained data for the pressure distribution in front of the step, and on its face and top surface, and the results of probe measurements in the separation and adjacent regions provide a more accurate scheme of the flow. The obtained data are compared with the results of other investigations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 73–80, September–October, 1977.We express our thanks to A. M. Kharitonov for valuable comments made during the discussion of this work, and also to M. A. Gol'dfel'd for kindly providing the experimental data for axisytnmetric steps.     

Supersonic nonequilibrium flow past segmental bodies of a mixture simulating the atmosphere of venus

Calculations of the flow of the mixture 0.94 CO₂+0.05 N₂+0.01 Ar past the forward portion of segmentai bodies are presented. The temperature, pressure, and concentration distributions are given as a function of the pressure ahead of the shock wave and the body velocity. Analysis of the concentration distribution makes it possible to formulate a simplified model for the chemical reaction kinetics in the shock layer that reflects the primary flow characteristics. The density distributions are used to verify the validity of the binary similarity law throughout the shock layer region calculated.The flow of a CO₂+N₂+Ar gas mixture of varying composition past a spherical nose was examined in [1]. The basic flow properties in the shock layer were studied, particularly flow dependence on the free-stream CO₂ and N₂ concentration.New revised data on the properties of the Venusian atmosphere have appeared in the literature [2, 3] One is the dominant CO₂ concentration. This finding permits more rigorous formulation of the problem of blunt body motion in the Venus atmosphere, and attention can be concentrated on revising the CO₂ thermodynamic and kinetic properties that must be used in the calculation.The problem of supersonic nonequilibrium flow past a blunt body is solved within the framework of the problem formulation of [4].Notation V body velocity - shock wave standoff - universal gas constant - ratio of frozen specific heats - hRt/m enthalpy per unit mass undisturbed stream P pressure - density - T temperature - m molecular weight - c_p specific heat at constant pressure - (X) concentration of component X (number of particles in unit mass) - R body radius of curvature at the stagnation point - _j rate of j-th chemical reaction shock layer P V ² pressure - density - TT temperature - mm molecular weight Translated from Izv. AN SSSR. Mekhanika Zhidkosti i Gaza, Vol. 5, No. 2, pp. 67–72, March–April, 1970.The author thanks V. P. Stulov for guidance in this study.     

Properties of two- and three-dimensional separation flows at supersonic velocities

The results are given of experimental investigations into three-dimensional separation of a turbulent boundary layer in the neighborhood of oblique shock waves, wedge-shaped obstacles, and sweptback steps at Mach numbers M = 2, 2.25, 2.5, 3, 4 and Reynolds numbers Re = u/v = (30–36)· 10⁶ m^–1. The characteristic regimes of the separated flows are considered. There is a discussion of the results of comparison and generalization of the pressure distribution in the two- and three-dimensional separation regions, and empirical dependences are also given for determining some geometrical parameters of these regions. An analogy is found in the characteristic pressures, and pressure distribution for a number of two- and three-dimensional separation flows, which suggests that one could use some of the known methods of analysis of two-dimensional separation of a turbulent boundary layer to calculate estimates for the three-dimensional case. This is confirmed by a comparison of calculated and experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–50, May–June, 1979.I am grateful to A. M. Kharitonov and V. S. Dem'yanenko for helpful comments made during the work and for discussion of this paper, and also to V. M. Filatov for assistance in some of the experiments.     

Hodograph design of lifting airfoils with high critical mach numbers

We wish to construct airfoils that have the highest free-stream Mach number M for a given set of geometric constraints for which the flow is nowhere supersonic. Nonlifting airfoils which maximize M for a given thickness ratio are known to possess long sonic segments at their critical speed. To construct lifting airfoils, we proceed under the conjecture that the optimal airfoil satisfying a given set of constraints is the one possessing the longest possible arc length of sonic velocity. A boundary-value problem is formulated in the hodograph plane using transonic small-disturbance theory whose solution determines an airfoil with long sonic arcs. For small lift coefficients, the hodograph domain covers two Riemann sheets and a finite-difference method is used to solve the boundary-value problem on this domain. A numerical integration of the solution around the boundary yields an airfoil shape, and three examples are discussed. The performance of these airfoils is compared with standard airfoils having the same lift coefficient and , and it is shown that the calculated airfoils have a 6%–10% increase in critical M .This research was supported by the Air Force Office of Scientific Research under Grant F49620-93-1-0022DEF and by the National Science Foundation under Grant DMS-9157546. Research support for M.C.A.K. was given by the IBM Corporation under a Graduate Research Fellowship.     

Asymptotic analysis of equilibrium states for rotating turbulent flows

The equilibrium states of homogeneous turbulence simultaneously subjected to a mean velocity gradient and a rotation are examined by using asymptotic analysis. The present work is concerned with the asymptotic behavior of quantities such as the turbulent kinetic energy and its dissipation rate associated with the fixed point (/kS)=0, whereS is the shear rate. The classical form of the model transport equation for (Hanjalic and Launder, 1972) is used. The present analysis shows that, asymptotically, the turbulent kinetic energy (a) undergoes a power-law decay with time for (P/)<1, (b) is independent of time for (P/)=1, (c) undergoes a power-law growth with time for 1<(P/)<(C ₂–1), and (d) is represented by an exponential law versus time for (P/)=(C ₂–1)/(C ₁–1) and (/kS)>0 whereP is the production rate. For the commonly used second-order models the equilibrium solutions forP/,II, andIII (whereII andIII are respectively the second and third invariants of the anisotropy tensor) depend on the rotation number when (P/kS)=(/kS)=0. The variation of (P/kS) andII versusR given by the second-order model of Yakhot and Orzag are compared with results of Rapid Distortion Theory corrected for decay (Townsend, 1970).     

Flow visualization of compressible vortex structures using density gradient techniques

Mathematical results are derived for the schlieren and shadowgraph contrast variation due to the refraction of light rays passing through two-dimensional compressible vortices with viscous cores. Both standard and small-disturbance solutions are obtained. It is shown that schlieren and shadowgraph produce substantially different contrast profiles. Further, the shadowgraph contrast variation is shown to be very sensitive to the vortex velocity profile and is also dependent on the location of the peak peripheral velocity (viscous core radius). The computed results are compared to actual contrast measurements made for rotor tip vortices using the shadowgraph flow visualization technique. The work helps to clarify the relationships between the observed contrast and the structure of vortical structures in density gradient based flow visualization experiments.Nomenclature a Unobstructed height of schlieren light source in cutoff plane, m - c Blade chord, m - f Focal length of schlieren focusing mirror, m - C _T Rotor thrust coefficient, T/( ² R ⁴) - I Image screen illumination, Lm/m ² - l Distance from vortex to shadowgraph screen, m - n _b Number of blades - p Pressure,N/m ² - p Ambient pressure, N/m ² - r, , z Cylindrical coordinate system - r _c Vortex core radius, m - Non-dimensional radial coordinate, (r/r _c ) - R Rotor radius, m - Tangential velocity, m/s - Specific heat ratio of air - Circulation (strength of vortex), m ²/s - Non-dimensional quantity, ² 8²p r _c ² - Refractive index of fluid medium - ₀ Refractive index of fluid medium at reference conditions - Gladstone-Dale constant, m ³/kg - Density, kg/m ³ - Density at ambient conditions, kg/m ³ - Non-dimensional density, (/ ) - Rotor solidity, (n _b c/ R) - Rotor rotational frequency, rad/s     

Flow in vicinity of blunt body stagnation point in diverging hypersonic stream

In the hypersonic thin shock layer approximation for a small ratio k of the densities before and after the normal shock wave the solution of [1] for the vicinity of the stagnation point of a smooth blunt body is extended to the case of nonuniform outer flow. It is shown that the effect of this nonuniformity can be taken into account with the aid of the effective shock wave radius of curvature R_*, whose introduction makes it possible to reduce to universal relations the data for different nonuniform outer flows with practically the same similarity criterion k. The results of the study are compared with numerical calculations of highly underexpanded jet flow past a sphere.Notations x, y a curvilinear coordinate system with axes directed respectively along and normal to the body surface with origin at the forward stagnation point - R radius of curvature of the meridional plane of the body surface - uV, vV., , p V ² respectively the velocity projections on the x, y axes, density, and pressure - and V freestream density and velocity The indices =0 and=1 apply to plane and axisymmetric flows Izv. AN SSSR, Mekhanika Zhidkosti i Gaza, Vol. 5, No. 3, pp. 102–105, 1970.     

Interaction between a distributed source and the free surface of a viscous fluid

A self-similar solution of the Navier-Stokes equations describing steady-state axisymmetric viscous incompressible fluid flow in a half-space is investigated. The motion is induced by sources or sinks distributed over a vertical axis with a constant density. The horizontal plane bounding the fluid is a free surface. It is found that in the presence of sources a solution of the above type exists and is unique for any value of the Reynolds numberR > 0, but in the case of sinks only on the interval –2 R < 0. The results of calculating the self-similar solutions are presented. The asymptotics of the solutions are found asR 0 andR .Novosibirsk. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 53–65, March–April, 1996.     

The finite element solutions of laminar flow and combined convection of air in a staggered or an in-line tube-bank

A finite element method is used to solve the full Navier-Stokes and energy equations for the problems of laminar combined convection from three isothermal heat horizontal cylinders in staggered tube-bank and four isothermal heat horizontal cylinders in in-line tube-bank. The variations of surface shear stress, pressure and Nusselt number are obtained over the entire cylinder surface including the zone beyond the separation point. The predicted values of total, pressure and friction drag coefficients, average Nusselt number and the plots of velocity flow fields and isotherms are also presented.

---

Die Finite-Elemente-Lösung von laminarer Strömung und kombinierter Konvektion von Luft in einer versetzten oder fluchtenden Rohranordnung

Zusammenfassung Eine Methode der finiten Elemente wird zur Lösung der vollständigen Navier-Stokes- und der Energiegleichung für die Probleme der laminaren kombinierten Konvektion an drei isothermen geheizten horizontalen Zylindern in versetzter Rohranordnung sowie für vier isotherme geheizte horizontale Zylinder in fluchtender Anordnung verwendet.Die Veränderung der Wandschubspannung, des Druckes und der Nusselt-Zahl werden für die gesamte Zylinderoberfläche, einschließlich des Bereiches nach dem Ablösepunkt, bestimmt. Die Werte des gesamten Widerstandsbeiwertes aufgrund von Druck und Reibung, die durchschnittliche Nusselt-Zahl und die Diagramme des Geschwindigkeitsfeldes und der Isothermen werden ebenfalls aufgezeigt.

---

Nomenclature C specifie heat - C _D total drag coefficient - C _f friction drag coefficient - C _p pressure drag coefficient - D diameter of cylinder,L=2R ₀ - G, g gravitational acceleration - Gr Grashof number, g(T_w–T )D ³/v ² - h local heat transfer coefficient - K thermal conductivity - L spacing between the centers of cylinder - M _l shape function - N _i shape function - Nu, local and average Nusselt numbers - P dimensionless pressure, p^*/u ² - p ^*,p pressure, free stream pressure - Pe Peclet number,RePr - Pr Prandtl number, c/K - Ra Rayleigh number,Gr Pr - Re Reynolds number,Du /v - R ₀ radius of cylinder - T temperature - T _w temperature on cylinder surface with fixed value - T free stream temperature - v dimensionless x-direction component of velocity,v ^*/u - u ^* x-direction component of velocity - u free stream velocity - v dimensionless Y-direction component of velocity,v ^*/u - v ^* Y-direction component of velocity - X x-direction axis - x dimensionless x-direction coordinate,x ^*/D - x^* x-direction coordinate - Y Y-direction axis - y dimensionless Y-direction coordinate,y ^*/D - y ^* Y-direction coordinate Greek symbols coefficient of volumetric thermal expansion - plane angle - dynamic viscosity - kinematic viscosity, / - density of fluid - _w dimensionless surface shear stress, _* ^w /u ² - skw/* surface shear stress - dimensionless temperature,      

Scalar Minimizers with Fractal Singular Sets

Lack of regularity of local minimizers for convex functionals with non-standard growth conditions is considered. It is shown that for every >0 there exists a function aC() such that the functional admits a local minimizer uW^1,p() whose set of non-Lebesgue points is a closed set with dim()>N–p–, and where 1<p<N<N+<q<+.     

Investigation of Dynamic Stability Regimes for Axisymmetric Bodies in Supersonic Flow

The problem of supersonic flow past a slender blunt cone with allowance for the reverse boundary-layer effect on the outer flow is solved with the aim of studying the influence of the boundary layer on the damping coefficient of axisymmetric body oscillations. It is assumed that the body executes plane angular, both low-amplitude and low-velocity, oscillations about a center of rotation. A modified version of the method [1] is applied for calculating the time-dependent flow past a body with the viscosity effect taken into account. The high accuracy of the flow parameter determination provided by this technique is confirmed by wind- tunnel experiments on a large-scale cone model (L1 m) at Mach numbers M=4 and 6. The agreement between the calculated and measured data forms the basis for the numerical investigation of the blunt-cone damping coefficient over a wide range of freestream Mach (M=4–20) and Reynolds (Re _L =10⁶–10⁸) numbers. At moderate freestream Mach numbers (M=4 and 6) an appreciable Re _L effect on the damping coefficient was not detected. However, on the hypersonic range this effect manifests itself more strongly, especially when there is gas injection into the boundary layer from the vehicle surface.     

Effect of the flow regime in an underexpanded jet issuing into a subsonic counterstream on the pressure distribution on a body

The pressure distribution on a cone with a half-angle =75°, from which a single central underexpanded jet issues into a subsonic counterstream, has been experimentally investigated. The effect of the flow regime in the jet on the pressure distribution is demonstrated. Generalized relations for the pressure on the body are obtained for various jet-flow momentum ratios J and flow Mach numbers M = 0.35–0.9; the Mach number M_a at the exit of the conical nozzle with half-angle _a=10° was equal to 2.9. The working medium of the jet and the flow was air with stagnation temperatures T_0a = T₀ 260–265°K. The ratio of the nozzle outlet radius to the radius of the maximum cross section of the cone _a/R_M=0.1.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 182–185, January-February, 1988.