Influence of turbulent Prandtl number on heat transfer of a flat plate

In computations involving heat transfer in turbulent flow past bodies it is necessary to assume turbulent Prandtl number distribution across the boundary layer. A review and comparison of results obtained by different authors are given, e.g., in [1–5]. Unfortunately, the existing data are so contradictory that, at present, it does not appear to be possible to establish reliably a function that determines turbulent Prandtl number distribution across the boundary layer. The absence of sufficiently reliable and general results on the distribution of turbulent Prandtl number led to the result that in the majority of studies conducted in earlier years its value was assumed a constant and either close to or equal to one. The effect of turbulent Prandtl number on the intensity of heat transfer from a flat plate is numerically investigated in the present paper. The thermal, turbulent boundary layer equation is integrated for this purpose at different values of turbulent Prandtl number and results are compared with experimental data. Results from [6], where the thermal boundary layer was numerically integrated with Pr_t=1 and compared with experimental data, were used for comparison in the present paper. The same numerical integration procedure as in [6] was used here.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 81–85, July–August, 1984.     

Use of turbulent energy balance equation in the theory of turbulent flows along walls

The turbulent flow of an incompressible fluid is considered in a plane channel, a circular tube, and the boundary layer on a flat plate. The system of equations describing the motion of the fluid consists of the Reynolds equations and the mean kinetic energy balance equation for turbulent fluctuations. On the basis of an analysis of experimental data, hypotheses are formulated with respect to the eddy kinematic viscosity and lengthl entering into the expression for specific dissipation of turbulent energy into heat. It is assumed that in the central (outer) region of the flow in a channel, andl are constants, and expressions are taken for them which are used for a free boundary layer; near the walll varies linearly and almost linearly. Results of calculations of the turbulent energy distribution, the mean velocity, and the drag coefficient are in good agreement with the existing experimental data. The values of two empirical coefficients, which enter into the system of equations as the result of the hypotheses, are close to those obtained for a free boundary layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 25–33, May–June, 1973.     

Boundary layer in the vicinity of the stagnation point of a body of axial symmetry in a turbulent stream

The flow in the boundary layer in the vicinity of the stagnation point of a flat plate is examined. The outer stream consists of turbulent flow of the jet type, directed normally to the plate. Assumptions concerning the connection between the pulsations in velocity and temperature in the boundary layer and the average parameters chosen on the basis of experimental data made it possible to obtain an isomorphic solution of the boundary layer equations. Equations are obtained for the friction and heat transfer at the wall in the region of gradient flow taking into account the effect of the turbulence of the impinging stream. It is shown that the friction at the wall is insensitive to the turbulence of the impinging stream, while the heat transfer is significantly increased with an increase in the pulsations of the outer flow. These properties are confirmed by the results of experimental studies [1–4].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 83–87, September–October, 1973.     

Characteristics of an unstably localized disturbance in a compressed boundary layer

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

Numerical investigation of a spatially developing turbulent natural convection boundary layer along a vertical heated plate

Large-eddy simulation (LES) on a spatially developing natural convection boundary layer along a vertical heated plate was conducted. The heat transfer rate, friction velocity, mean velocity and temperature, and second-order turbulent properties both in the wall-normal and the stream-wise direction showed reasonable agreement with the findings of past experiments. The spectrum of velocity and temperature fluctuation showed a -2/3-power decay slope and -2-power decay slope respectively. Quadrant analysis revealed the inclination on Q1 and Q3 in the Reynolds stress and turbulent heat flux, changing their contribution along the distance from the plate surface. Following the convention, we defined the threshold region where the stream-wise mean velocity takes local maximum, the inner layer which is closer to the plate than the threshold region, the outer layer which is farther to the plate than the threshold region. The space correlation of stream-wise velocity tilted the head toward the wall in the propagating direction in the outer layer; on the other hand, the correlated motion had little inclination in the threshold region. The time history of the second invariant of gradient tensor Q revealed that the vortex strength oscillates both in the inner and the outer layers in between the laminar and the transition region. In the turbulent region, the vortex was often dominant in the outer layer. Instantaneous three-dimensional visualization of Q revealed the existence of high-speed fluid parcels associated with arch-shape vortices. These results were considered as an intrinsic structure in the outer layer, which is symmetrical to the structure of canonical smooth/rough wall bounded layer flow in forced convection.     

DNS and scaling laws from new symmetry groups of ZPG turbulent boundary layer flow

The Lie group, or symmetry approach, developed by Oberlack (see e.g. Oberlack [26] and references therein) is used to derive new scaling laws for various quantities of a zero pressure gradient turbulent boundary layer flow. The approach unifies and extends the work done by Oberlack for the mean velocity of stationary parallel turbulent shear flows. From the two-point correlation (TPC) equations the knowledge of the symmetries allows us to derive a variety of invariant solutions (scaling laws) for turbulent flows, one of which is the new exponential mean velocity profile that is found in the mid-wake region of flat-plate boundary layers. Further, a third scaling group was found in the TPC equations for the one-dimensional turbulent boundary layer. This is in contrast to the Navier–Stokes and Euler equations, which have one and two scaling groups, respectively. The present focus is on the exponential law in the outer region of turbulent boundary layer corresponding new scaling laws for one- and two-point correlation functions. A direct numerical simulation (DNS) of a flat plate turbulent boundary layer with zero pressure gradient was performed at two different Reynolds numbers Re=750,2240. The Navier–Stokes equations were numerically solved using a spectral method with up to 140 million grid points. The results of the numerical simulations are compared with the new scaling laws. TPC functions are presented. The numerical simulation shows good agreement with the theoretical results, however only for a limited range of applicability. PACS 02.20.-a, 47.11.+j, 47.27.Nz, 47.27.Eq     

Turbulent mixing of relaxing gases in a supersonic nozzle

Processes of turbulent mixing of hot nitrogen and carbon dioxide gas (CO₂) in a supersonic nozzle are investigated in connection with gasdynamic lasers with selective thermal excitation. The actual flow pattern in the nozzle and cavity is simulated by the system of equations of a turbulent boundary layer. The results of calculations of the gain profile, the CO₂ concentration, and the laser power agree satisfactorily with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 156–160, March–April, 1981.We are very grateful to B. S. Aleksandrov and V. K. Pozdyshev for numerous helpful discussions.     

DNS of a spatially developing turbulent boundary layer with passive scalar transport

A direct numerical simulation (DNS) of a spatially developing turbulent boundary layer over a flat plate under zero pressure gradient (ZPG) has been carried out. The evolution of several passive scalars with both isoscalar and isoflux wall boundary condition are computed during the simulation. The Navier–Stokes equations as well as the scalar transport equation are solved using a fully spectral method. The highest Reynolds number based on the free-stream velocity U_∞ and momentum thickness θ is Re_θ=830, and the molecular Prandtl numbers are 0.2, 0.71 and 2. To the authors’ knowledge, this Reynolds number is to date the highest with such a variety of scalars. A large number of turbulence statistics for both flow and scalar fields are obtained and compared when possible to existing experimental and numerical simulations at comparable Reynolds number. The main focus of the present paper is on the statistical behaviour of the scalars in the outer region of the boundary layer, distinctly different from the channel-flow simulations. Agreements as well as discrepancies are discussed while the influence of the molecular Prandtl number and wall boundary conditions is also highlighted. A Pr scaling for various quantities is proposed in outer scalings. In addition, spanwise two-point correlation and instantaneous fields are employed to investigate the near-wall streak spacing and the coherence between the velocity and the scalar fields. Probability density functions (PDF) and joint probability density functions (JPDF) are shown to identify the intermittency both near the wall and in the outer region of the boundary layer. The present simulation data will be available online for the research community.     

Heat transfer at and near the stagnation point in turbulent flow over bodies

Results are presented of experimental investigations of heat transfer in the neighborhood of the stagnation point in flow of a turbulent gas over bodies. It is assumed that the outer flow is capable of rendering the boundary layer turbulent over the whole body surface, i.e., the hypothesis is invoked that there is a turbulent stagnation point. Using the method of integral relations [1] and the flat plate heat-transfer law, transformed in such a way as to satisfy the heat-transfer conditions at the stagnation point, simple formulas have been obtained for calculating the heat flux.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 177–181, July–August, 1975.     

Nonmonotonic relationship between the Reynolds number and the lift coefficient of a plate in a hypersonic rarefied gas flow

The aerodynamic coefficients of a plate in a hypersonic stream at low Reynolds numbers are studied over a wide range of similarity parameters. The dependence of the lift coefficientC _Y on the tangential force coefficient, the finite Mach number at the outer edge of the boundary layer and the velocity-slip and temperature-jump boundary conditions is taken into consideration. The nonmonotonic character of the relationship betweenC _Y and the Reynolds number, revealed previously in experiments, is explained within the framework of the viscous hypersonic interaction model.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 186–189, January–February, 1996.     

Bulk flow parameters in turbulent wall boundary layers near separation

Experiments were conducted in a turbulent boundary layer near separation along a flat plate. The pressure gradient in flow direction was varied such that three significant boundary layer configurations could be maintained. The flow in the test section thus had simultaneously a region of favourable pressure gradient, a region of strong adverse pressure gradient with boundary layer separation and a region of reattached boundary layer. Specially designed fine probes facilitated the measurements of skin friction and velocity distribution very close to the wall. Bulk flow parameters such as skin friction coefficient C _f, Reynold's number Re_δ2 and shape factors H and G, which are significant characteristics of wall boundary layers were evaluated. The dependence of these parameters on the Reynolds number and along the test section was explored and the values were compared with other empirical and analytical formulae known in the literature.     

Cone in a supersonic flow near a surface with a turbulent boundary layer

The published information about the interaction of incident shocks and a turbulent boundary layer relate to cases of a thin boundary layer ( 1–3 mm) on a flat surface. The present study relates to supersonic flow with Mach number M = 3 and stagnation pressure p₀=1.2 MPa past cones near a surface with a thick boundary layer formed on a plate abutting the lower edge of a plane nozzle.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 177–180, July–August, 1991.     

Quantitative visualization of compressible turbulent shear flows using condensate-enhanced Rayleigh scattering

This paper describes several flow visualization experiments carried out in Mach 3 and Mach 8 turbulent shear flows. The experimental technique was based on laser scattering from particles of H₂O or CO₂ condensate that form in the wind tunnel nozzle expansion process. The condensate particles vaporize extremely rapidly on entering the relatively hot fluid within a turbulent structure, so that a sharp vaporization interface marks the outer edge of the rotational shear layer fluid. Calculations indicate that the observed thin interface corresponds to a particle size of 10 nm or less, which is consistent with optical measurements, and that particles of this size track the fluid motions well. Further, calculations and experiments show that the freestream concentration of condensate required for flow visualization has only a small effect on the wind tunnel pressure distribution. Statistics based on the image data were compared to corresponding results from probe measurements and agreement was obtained in statistical measures of speed, scale, and orientation of the large-scale structures in the shear layer turbulence. The condensate-enhanced Rayleigh scattering technique is judged to be a useful tool for quantitative studies of shear layer structure, particularly for identifying the instantaneous boundary layer edge and for extracting comparative information on the large-scale structures represented there.     

Three-Dimensional Simulation of Turbulent Hot-Jet Ignition for Air-CH<Subscript>4</Subscript>-H<Subscript>2</Subscript> Deflagration in a Confined Volume

This work describes essential aspects of the ignition and deflagration process initiated by the injection of a hot transient gas jet into a narrowly confined volume containing air-CH₄-H₂ mixture. Driven by the pressure difference between a prechamber and a long narrow constant-volume-combustion (CVC) chamber, the developing jet or puff involves complex processes of turbulent jet penetration and evolution of multi-scale vortices in the shear layer, jet tip, and adjacent confined spaces. The CVC chamber contains stoichiometric mixtures of air with gaseous fuel initially at atmospheric conditions. Fuel reactivity is varied using two different CH₄/H₂ blends. Jet momentum is varied using different pre-chamber pressures at jet initiation. The jet initiation and the subsequent ignition events generate pressure waves that interact with the mixing region and the propagating flame, depositing baroclinic vorticity. Transient three-dimensional flow simulations with detailed chemical kinetics are used to model CVC mixture ignition. Pre-ignition gas properties are then examined to develop and verify criteria to predict ignition delay time using lower-cost non-reacting flow simulations for this particular case of study.     

Ignition of a gas in a boundary layer at a heated plate

The problem of the ignition of a moving homogeneous gaseous combustible mixture in a boundary layer along a heated flat semiinfinite plate is one of the main problems of the ignition of a combustible mixture in a flow (for example, [1]). The formulation of the problem includes the two-dimensional equations of motion and the equations of the transfer of heat and of the reacting substance, written taking a chemical reaction into consideration, as well as boundary conditions, and should lead to determination of the steady-state fields of the concentration and the temperature and, by the same token, of the position of the combustion zone. Different approximate numerical solutions of the problem were analyzed in [1–5]. One of the most important characteristics of the process is the length of the ignition, i.e., the distance from the edge of the plate to the point at which, thanks to the intrinsic chemical heat evolution in the gas, the heat flux from the plate to the gas becomes equal to zero. In the present work, for the case of large values of the activation energy of the chemical reaction and a sufficiently great temperature difference between the wall and the flow, an approximate expression is obtained for the length of the ignition.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 142–148, September–October, 1977.The authors thank V. M. Shevtsov for his aid in making the calculations.     

Numerical prediction of turbulent boundary layer noise from a sharp-edged flat plate

An efficient hybrid uncorrelated wall plane waves–boundary element method (UWPW-BEM) technique is proposed to predict the flow-induced noise from a structure in low Mach number turbulent flow. Reynolds-averaged Navier-Stokes equations are used to estimate the turbulent boundary layer parameters such as convective velocity, boundary layer thickness, and wall shear stress over the surface of the structure. The spectrum of the wall pressure fluctuations is evaluated from the turbulent boundary layer parameters and by using semi-empirical models from literature. The wall pressure field underneath the turbulent boundary layer is synthesized by realizations of uncorrelated wall plane waves (UWPW). An acoustic BEM solver is then employed to compute the acoustic pressure scattered by the structure from the synthesized wall pressure field. Finally, the acoustic response of the structure in turbulent flow is obtained as an ensemble average of the acoustic pressures due to all realizations of uncorrelated plane waves. To demonstrate the hybrid UWPW-BEM approach, the self-noise generated by a flat plate in turbulent flow with Reynolds number based on chord Re_c = 4.9 × 10⁵ is predicted. The results are compared with those obtained from a large eddy simulation (LES)-BEM technique as well as with experimental data from literature.     

The turbulent/non-turbulent interface at the outer boundary of a self-similar turbulent jet

The flow at the outer boundary of a submerged self-similar turbulent jet at Re=2᎒³ is investigated experimentally by means of combined particle image velocimetry (PIV) laser-induced fluorescence (LIF) measurements. The jet fluid contains a fluorescent dye so that the LIF data can be used to discriminate between the jet fluid and the ambient fluid. The axial velocity, Reynolds stress, and vorticity are determined relative to the jet boundary. The results are compared against the conventional profiles, and the results of a direct numerical simulation of the turbulent far-wake behind a flat plate. The results show a sharp transition between rotational and irrotational fluid at the fluid interface, and the existence of a layer of irrotational velocity fluctuations outside the turbulent region.     

Effect of nonisothermality of the surface of a thin profile on the stability of a laminar boundary layer

Investigations of the stability of a subsonic laminar boundary layer have shown that, other things being equal, the stability of the laminar flow is considerably improved by cooling the entire surface of the body to a constant temperature T_w=const lower than the temperature of the free stream [1–3]. This is attributable to an increase in the critical Reynolds number of loss of stability and a decrease in the range of unstable perturbations of the Tollmien-Schlichting wave type when the surface is cooled. Recently, in the course of investigating the stability of laminar flow over a flat plate it was found [4, 5] that a similar improvement in flow stability can be achieved by raising the temperature of a small part of the surface near the leading edge of the plate. In this study we examine the possibility of delaying the transition to turbulent flow by creating a nonuniform temperature distribution along the surface of thin profiles, where the development of an adverse pressure gradient in the outer flow has a destabilizing effect on the boundary layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 36–42, September–October, 1986.In conclusion, the authors wish to thank M. N. Kogan for useful discussions of their results.     

Statistics of Particle Accumulation in Spatially Developing Turbulent Boundary Layers

We present the results of a Direct Numerical Simulation of a particle-laden spatially developing turbulent boundary layer up to Re _θ ?=?2500. Two main features differentiate the behavior of inertial particles in a zero-pressure-gradient turbulent boundary layer from the more commonly studied case of a parallel channel flow. The first is the variation along the streamwise direction of the local dimensionless parameters defining the fluid-particle interactions. The second is the coexistence of an irrotational free-stream and a near-wall rotational turbulent flow. As concerns the first issue, an inner and an outer Stokes number can be defined using inner and outer flow units. The inner Stokes number governs the near-wall behavior similarly to the case of channel flow. To understand the effect of a laminar-turbulent interface, we examine the behavior of particles initially released in the free stream and show that they present a distinct behavior with respect to those directly injected inside the boundary layer. A region of minimum concentration occurs inside the turbulent boundary layer at about one displacement thickness from the wall. Its formation is due to the competition between two transport mechanisms: a relatively slow turbulent diffusion towards the buffer layer and a fast turbophoretic drift towards the wall.     

Near-wall flow in a three-dimensional boundary layer on the endwall of a 30° bend

 Turbulence measurements are reported on the three-dimensional turbulent boundary layer along the centerline of the flat endwall in a 30° bend. Profiles of mean velocities and Reynolds stresses were obtained down to y ⁺≈2 for the mean flow and y ⁺≈8 for the turbulent stresses. Mean velocity data collapsed well on a simple law-of-the-wall based on the magnitude of the resultant velocity. The turbulence intensity and turbulent shear stress magnitude both increased with increased three-dimensionality. The ratio of these two quantities, the a ₁ structure parameter, decreased in the central regions of the boundary layer and showed profile similarity for y ⁺<50. The shear stress vector angle lagged behind the velocity gradient vector angle in the outer region of the boundary layer, however there was an indication that the shear stress vector tends to lead the velocity gradient vector close to the wall. Received: 16 July 1996/Accepted: 14 July 1997