Stability and Transition on a Swept Cylinder in a Supersonic Flow

Results of experimental investigations of the evolution of natural disturbances and laminar–turbulent transition in a supersonic boundary layer on the attachment line of a circular cylinder with a sweep angle of 68° and a freestream Mach number M = 2 are presented. The experimental studies are supplemented by calculations of the mean flow and stability characteristics. Flow regimes in the boundary layer on the attachment line are determined by a hotwire technique as functions of the Reynolds number and height of twodimensional roughness elements. The results are compared with NASA (Ames) experiments.     

Laminar boundary layer of infinitely long elliptic cylinders for an arbitrary yaw angle

Results are presented of the calculation of the laminar boundary layer on infinitely long elliptic cylinders in a supersonic perfect gas flow at an arbitrary angle of attack. It was assumed that the Prandtl number is constant and equal to 0.7, the dynamic viscosity coefficient follows a power-law variation ( T^0.76) with temperature, and there is high heat transfer at the body surface (H_1w=0.05).The calculations showed that a change of the body shape—the ellipticity coefficient =b/a—has a significant effect on the nature of the distribution and the magnitude of the local heat flux.In evaluating the thermal fluxes at the blunt leading edges, swept wings are usually considered as infinitely long yawed cylinders. In studying heat transfer at the surface of bodies of small aspect ratio at high angles of attack, wide use is made of the hypothesis of plane sections, when each section, orthogonal to the longitudinal axis of the body, is considered equivalent to a corresponding yawed infinite cylinder.By now quite detailed studies have been made of the behavior of the boundary layer on an infinitely long yawed circular cylinder with both the laminar and turbulent flow regimes for a compressible gas [1, 2]. However, there are no data on the heat transfer at the surface of a yawed infinite cylinder with arbitrary cross section, although the availability of such data is urgently needed, for example, for the proper selection of the form of the leading edges of the swept wing.This article presents the results of the calculation of the characteristics of the laminar boundary layer on the surface of infinite elliptic cylinders in a supersonic perfect gas flow. The calculations were made over a quite wide range of flight Mach number M, yaw angle , and ellipticity factor . The data presented on the distribution of the relative heat flux along the cylinder directrix may be used also for estimating the heat flux with account for the real properties of air if we know the corresponding value of the heat flux in the vicinity of the stagnation line.     

Structure of the disturbed boundary layer near a cold surface

A study is made of a laminar boundary layer near a cold plate in the regime of weak viscous -inviscid interaction in the limit in which the free-stream Mach number tends to infinity and the temperature factor to zero. Local disturbed flow regions formed due to the presence of small elements of roughness on the surface of the plate are investigated. It is found that thick regions of roughness elements induce disturbances of the frictional stresses and the heat flux of the same order as these quantities in the undisturbed boundary layer, while all thin regions of roughness elements induce only small, linear perturbations and, therefore, cannot cause separation of the boundary layer; the different regimes of flow past the roughness elements are described.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 19–27, January–February, 1991.     

Heat transfer in a thin layer of a viscous heavy noncompressible liquid under wavy flow conditions

The article describes a method for calculating the flow of heat through a wavy boundary separating a layer of liquid from a layer of gas, under the assumption that the viscosity and heat-transfer coefficients are constant, and that a constant temperature of the fixed wall and a constant temperature of the gas flow are given. A study is made of the equations of motion and thermal conductivity (without taking the dissipation energy into account) in the approximations of the theory of the boundary layer; the left-hand sides of these equations are replaced by their averaged values over the layer. These equations, after linearization, are used to determine the velocity and temperature distributions. The qualitative aspect of heat transfer in a thin layer of viscous liquid, under regular-wavy flow conditions, is examined. Particular attention is paid to the effect of the surface tension coefficient on the flow of heat through the interface.Notation x, y coordinates of a liquid particle - t time - v and u coordinates of the velocity vector of the liquid - p pressure in the liquid - c_v, , T,, andv heat capacity, thermal conductivity coefficient, temperature, density, and viscosity of the liquid, respectively - g acceleration due to gravity - surface-tension coefficient - c phase velocity of the waves at the interface - T_w wall temperature - h₀ thickness of the liquid layer - u₀ velocity of the liquid over the layer Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 147–151, July–August, 1970.     

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.     

Results of an experimental and numerical study of the aerodynamic heating of the undersurface of delta wings with sharp leading edges at mach numbers M∞=6.1 and 8

The heat transfer between a supersonic flow and the undersurface of delta wings with leading-edge sweep angles x=65 and 70° is investigated in a shock tunnel at angles of attack 15°. The supersonic inviscid flow over these wings in regimes in which the bow shock is attached to the sharp leading edges is calculated numerically. The compressible boundary layer problem is solved for the calculated inviscid flow fields in the laminar, transition and turbulent flow zones. The calculations and experimental values of the heat flux on the surface of the wings are compared. The calculations are in satisfactory agreement with the experimental data in the laminar and transition zones, but diverge significantly (by up to 20%) in the turbulent zone.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 183–188, July–August, 1991.The authors wish to thank A. A. Golubinskii for assisting with the solution of the problem of supersonic inviscid gas flow over a wing.     

Exact solution of the problem of convective heat exchange for an elliptic cylinder and a plate in a fluid with small Prandtl number

The heat exchange problem is solved for an elliptic cylinder and a plate in an incompressible fluid flow with small Prandtl number Pr1. For flow along a plate, the values of the Nusselt number Nu obtained by solving the complete energy equation and the heat boundary layer equation are compared.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 26–31, May–June, 1996.     

Numerical analysis of simulation accuracy for hypersonic heat transfer in subsonic jets of dissociated nitrogen

One-dimensional problems of the flow in a boundary layer of finite thickness on the end face of a model and in a thin viscous shock layer on a sphere are solved numerically for three regimes of subsonic flow past a model with a flat blunt face exposed to subsonic jets of pure dissociated nitrogen in an induction plasmatron [1] (for stagnation pressures of (10⁴–3·10⁴) N/m² and an enthalpy of 2.1·10⁷ m²/sec²) and three regimes of hypersonic flow past spheres with parameters related by the local heat transfer simulation conditions [2, 3]. It is established that given equality of the stagnation pressures, enthalpies and velocity gradients on the outer edges of the boundary layers at the stagnation points on the sphere and the model, for a model of radius R_m=1.5·10^–2 m in a subsonic jet the accuracy of reproduction of the heat transfer to the highly catalytic surface of a sphere in a uniform hypersonic flow is about 3%. For surfaces with a low level of catalytic activity the accuracy of simulation of the nonequilibrium heat transfer is determined by the deviations of the temperatures at the outer edges of the boundary layers on the body and the model. For this case the simulation conditions have the form: dU_e/dx=idem, p₀=idem, T_e=idem. At stagnation pressuresP ₀2·10⁴ N/m² irrespective of the catalycity of the surface the heat flux at the stagnation point and the structure of the boundary layer near the axis of symmetry of models with a flat blunt face of radius R_m1.5·10^–2 m exposed to subsonic nitrogen jets in a plasmatron with a discharge channel radius R_c=3·10^–2 m correspond closely to the case of spheres in hypersonic flows with parameters determined by the simulation conditions [2, 3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 135–143, March–April, 1990.     

Heat transfer enhancement with deformation of the velocity profile

The problem of enhancing the heat transfer in channels and boundary layers by the appropriate deformation of the fluid velocity profile is considered. The resulting additional hydraulic losses, the price of heat transfer enhancement, are determined. The possibilities of controlling heat transfer by redistributing the fluid velocity in channels are demonstrated with reference to flows at low Prandtl numbers. Laminar and turbulent liquid and gas flows with heat transfer in channels and boundary layers are numerically modeled on the basis of modern models of turbulence (flow development in channels with different initial velocity profiles, flows with wall roughness and boundary layer flows with forces acting on the flow to cause deformation of the velocity profile). In all cases it is found that the heat transfer can be enhanced only at the expense of a considerable increase in the hydaulic losses. A class of self-similar thermal problems for flows in plane diffusers is formulated. The eigenfunctions — temperature modes — for various velocity profiles are determined with allowance for the nonuniqueness of the solution of the classical dynamical problem for a plane diffuser and the corresponding heat transfer coefficients are found.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 94–105, May–June, 1993.The authors are grateful to A. Yu. Klimenko for useful discussions.     

Measurements in the separation region of a gas-particle cross flow

This paper reports the results from a study of the effects of particles on the boundary layer characteristics for cross flow over a single circular cylinder in the near separation region and on a cylinder located in the second row of a small mock-up heat exchanger array. Two size ranges of particles (45–85 m and 100–200 m) with two concentrations were used for the single cylinder tests. For the array, the larger particles were used at a single concentration. For the single cylinder, the results with the larger particles at mass loading ratio of 1.0 show an obvious effect on the velocity profile shape and on the rms levels. For reduced concentrations, and for the smaller particles, the effects are less noticeable. For the array, similar trends are observed and a delay in the separation point is reported.     

Interference between a blunt edge shock layer and an impinging inclined shock at low Reynolds numbers

The interference between the shock layer on a cylinder modeling the leading edge of an air intake and an impinging plane inclined shock is investigated experimentally and numerically for a Reynolds number Re₀=32. The low-pressure wind tunnel experiments made it possible to visualize the flow and determine the local heat transfer in the presence of interference. The corresponding flow regimes were calculated numerically within the framework of the system of Navier-Stokes equations by the through-calculation method. The principal properties of the distribution of the flow characteristics for a low value of the Reynolds number were obtained for various types of interference and the differences with respect to the previously investigated interference regimes for high Reynolds numbers were examined.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 166–171, January–February, 1993.     

Optical investigation of the heat transfer from a rotating cylinder in a cross flow

The heat transfer from a rotating cylinder in an air-cross flow was investigated by purely optical measuring techniques. Flow velocities were measured by a two-dimensional LDV both in the vicinity of the cylinder and in the boundary layer. A new optical device based on light-deflection in a temperature field was developed to examine local temperature gradients in the boundary layer of the rotating cylinder. Finally, a Michelson-interferometer was installed to produce real-time pictures of isothermal lines around the heated cylinder. The impact of rotation on flow patterns, boundary layer behaviour and heat transfer could be clearly identified. It appears that the velocity-ratio acts like an independent parameter, in that flow patterns correspond to this dimensionless number. Furthermore, it seems that rotation dominates over cross flow, both fluid-dynamically and thermally above = 2.This work was carried out at the University of the Federal Armed Forces in Munich/Germany.     

Flow of a viscous fluid in an annular gap with intensive blowing from the walls

Self-similar solutions are obtained in [1, 2] to the Navier-Stokes equations in gaps with completely porous boundaries and with Reynolds number tending to infinity. Approximate asymptotic solutions are also known for the Navier-Stokes equations for plane and annular gaps in the neighborhood of the line of spreading of the flow [3, 4]. A number of authors [5–8] have discovered and studied the effect of increase in the stability of a laminar flow regime in channels of the type considered and a significant increase in the Reynolds number of the transition from the laminar regime to the turbulent in comparison with the flow in a pipe with impermeable walls. In the present study a numerical solution is given to the system of Navier-Stokes equations for plane and annular gaps with a single porous boundary in the neighborhood of the line of spreading of the flow on a section in which the values of the local Reynolds number definitely do not exceed the critical values [5–8]. Generalized dependences are obtained for the coefficients of friction and heat transfer on the impermeable boundary. A comparison is made between the solutions so obtained and the exact solutions to the boundary layer equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 21–24, January–February, 1987.     

Receptivity of the inclined-cylinder attachment-line boundary layer to vortex perturbations

The receptivity of the boundary layer in the neighborhood of the attachment line of a cylinder inclined to the flow with respect to periodic vortex perturbations frozen into the stream is investigated. The problem considered simulates the interaction between external turbulence and the leading-edge swept wing boundary layer. It is shown that if the direction of the external perturbation vector is almost parallel to the leading edge, then the external perturbations are considerably strengthened at the outer boundary layer edge. This effect can cause laminar-turbulent transition on the attachment line at subcritical Reynolds numbers.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 72–85. Original Russian Text Copyright © 2004 by Ustinov.     

LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.     

Boundary layer transition on the surface of a delta wing in supersonic flow

Laminar-turbulent transition on the surface of a delta wing has been experimentally investigated in a supersonic wind tunnel at Mach numbers M_t8=3–5. It is shown that when M,=3, Re_L=6.5·10⁶, and =–5.5° much of the upper surface of the wing in the neighborhood of the line of symmetry is occupied by a wedge-shaped region of turbulent flow. In this region the heat fluxes reach the same values as at the heat transfer maxima induced here by separated flows and may exceed the turbulent heat flux level on the windward surface of the wing. Changing the shape of the under surface of the wing from plane to pyramidal leads to acceleration of the boundary layer transition on the under surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 87–92, May–June, 1989.     

Determination of specific thrust in unchoked regimes and design of a separationless convergent nozzle contour

The variation of the specific thrust R_Y on the angle of inclination of the wall is analyzed within the framework of the ideal gas model using the results of specific impulse and flow rate calculations for conical convergent nozzles. It is shown that in unchoked regimes nozzles with different have almost the same values of R_Y for both subcritical and supercritical pressure ratios _c. On the interval _C < 6 typical of convergent nozzles conical convergent nozzles with =30–90° have almost the same value of the specific thrust, maximal relative to the R_Y of nozzles with < 30°. In the presence of viscosity forces local boundary layer separation may occur in the neighborhood of the entrance section of the convergent nozzle. A method of constructing a separationless convergent nozzle contour with enhanced thrust is developed on the basis of a boundary layer separation criterion. The separationless contour is determined for given values of the flow rate, specific heat ratio, Reynolds number, wall temperature and initial boundary layer displacement thickness.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 158–164, January–February, 1990.     

Method of solution of the hypersonic boundary layer equations in the case of strong viscous-inviscid interaction

A method of solving the boundary layer equations is developed taking into account the strong interaction between the boundary layer and the outer hypersonic inviscid flow. The method is aimed at solving problems whose salient feature is the possible upstream propagation of disturbances over distances comparable with the body length. The procedure for fitting a self-consistent contour of the effective body using an artificially formulated boundary value problem for an ordinary second-order differential equation, which lies at the basis of the method, is considered in detail. The method is applied to the problem of flow around a flat plate with roughness in the form of an embankment or a trench; the calculated results are presented.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 81–89, July–August, 1995.     

Estimation of the acoustic emission spectrum of a boundary layer on a rough surface

Under engineering conditions the surfaces over which fluids flow are not usually hydraulically smooth. In this connection it is important to investigate the generation of sound by a turbulent boundary layer on a rough surface. Turbulent flow over a deformed surface creates dipole sources of sound, which may considerably increase the acoustic emission as compared with the quadrupole emission from a boundary layer on a smooth plate [1, 2]. In the case of sandy roughness estimates of the acoustic field are usually based on the energy summation of the fields generated by flow over the individual roughness elements [3, 4]. In this case not easily verifiable assumptions are made concerning the nature of the turbulent flow near the roughness, and the intensity of the emission is found correct to a constant factor subject to experimental determination. In the present study, in order to calculate the acoustic emission of a boundary layer on a surface with sandy roughness, it is proposed to employ the available experimental data on the cross of the surface pressure.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 20–26, September–October, 1985.     

Calculation of the averaged three-dimensional boundary layer at the axisymmetrical boundaries of the interblade channels of turbine machines

The article discusses the flow of a gas at the blade rim of an axial turbine, consisting of an external steady-state continuous flow of an ideal compressible liquid and a three-dimensional turbulent boundary layer of a compressible liquid at the end surfaces of the rim, averaged in a peripheral direction. It presents an example of a calculation of flow in fixed blades, with a different form of the meridional cross section. In a flow through the rim of a turbine machine between the convex and concave surfaces of adjacent blades there arises a transverse gradient of the static pressure. At the end surface in the boundary layer the lines of the flow are shifted toward the convex side of the profile, and a secondary transverse flow of the liquid arises [1–3]. The article discusses the following: an external two-dimensional steady-state adiabatic flow of an ideal compressible liquid at the surface S₂, which can be taken as the mean surface of the interblade channel, with boundary lines at the peripheral and root end surfaces of the rim; a two-dimensional steady-state adiabatic flow of an ideal compressible liquid at the end surfaces of the rim between the convex and concave sides of the profiles [3, 4]; and a three-dimensional turbulent boundary layer, averaged in a peripheral direction at the end surfaces of the blade rim. The averaged boundary layer is calculated along one coordinate line s, and a simplified model of the quasi-three-dimensional flow is used. The coefficients of friction and heat transfer, and the inclination of the bottom flow lines are averaged.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 22–31, May–June, 1975.The author thanks G. Yu. Stepanov for posing the problem and evaluating the results.