Stability of spatially developing boundary layers in pressure gradients

A new formulation of the stability of boundary-layer flows in pressure gradients is presented, taking into account the spatial development of the flow and utilizing a special coordinate transformation. The formulation assumes that disturbance wavelength and eigenfunction vary downstream no more rapidly than the boundary-layer thickness, and includes all terms nominally of order R(-1) in the boundary-layer Reynolds number R. In Blasius flow, the present approach is consistent with that of Bertolotti et al. (1992) to O(R(-1)) but simpler (i.e. has fewer terms), and may best be seen as providing a parametric differential equation which can be solved without having to march in space. The computed neutral boundaries depend strongly on distance from the surface, but the one corresponding to the inner maximum of the streamwise velocity perturbation happens to be close to the parallel flow (Orr-Sommerfeld) boundary. For this quantity, solutions for the Falkner-Skan flows show the effects of spatial growth to be striking only in the presence of strong adverse pressure gradients. As a rational analysis to O(R(-1)) demands inclusion of higher-order corrections on the mean flow, an illustrative calculation of one such correction, due to the displacement effect of the boundary layer, is made, and shown to have a significant destabilizing influence on the stability boundary in strong adverse pressure gradients. The effect of non-parallelism on the growth of relatively high frequencies can be significant at low Reynolds numbers, but is marginal in other cases. As an extension of the present approach, a method of dealing with non-similar flows is also presented and illustrated. However, inherent in the transformation underlying the present approach is a lower-order non-parallel theory, which is obtained by dropping all terms of nominal order R(-1) except those required for obtaining the lowest-order solution in the critical and wall layers. It is shown that a reduced Orr-Sommerfeld equation (in transformed coordinates) already contains the major effects of non-parallelism.

Transition to turbulence in the free convection boundary layers on an inclined heated plate

An experimental study has been made of transition to turbulence in the free convective flows on a heated plate. Observations have been made with the plate vertical and inclined at angles up to about 50° to the vertical, both above and below the plate. A fibre anemometer was used to survey the region of intermittent turbulence. Information has thus been obtained about the range of Grashof numbers over which transition takes place. Even when the plate is vertical the region of intermittent turbulence is long. When it is inclined, this region becomes still longer in the flow below the plate as a result of the stabilizing stratification, a Richardson number effect. It is possible to have a whole flow such that it should be described as transitional, not laminar or turbulent. It was noticed that in this flow and the vertical plate one, the velocity during the laminar periods could be either of two characteristic values, one of them close to zero. The behaviour above an inclined plate could be interpreted largely as a trend towards the behaviour described in a preceding paper.

Three-dimensional laminar compressible boundary layers with large injection

The effect of large mass injection on the following three-dimensional laminar compressible boundary-layer flows is investigated by employing the method of matched asymptotic expansions: (i) swirling flow in a laminar compressible boundary layer over an axisymmetric surface with variable cross-section and (ii) laminar compressible boundary-layer flow over a yawed infinite wing in a hypersonic flow. The resulting equations are solved numerically by combining the finite-difference technique with quasi-linearization. An increase in the swirl parameter, the yaw angle or the wall temperature is found to be capable of bringing the viscous layer nearer the surface and reducing the effects of massive blowing.

Effect of longitudinal surface curvature on boundary layers

We consider here the higher order effect of moderate longitudinal surface curvature on steady, two-dimensional, incompressible laminar boundary layers. The basic partial differential equations for the problem, derived by the method of matched asymptotic expansions, are found to possess similarity solutions for a family of surface curvatures and pressure gradients. The similarity equations obtained by this anaylsis have been solved numerically on a computer, and show a definite decrease in skin friction when the surface has convex curvature in all cases including zero pressure gradient. Typical velocity profiles and some relevant boundary-layer characteristics are tabulated, and a critical comparison with previous work is given.

Universality of scaling laws in correlation between velocity and shear stress in turbulent boundary layers

In this study, we analyse simultaneous measurements (at 50 Hz) of velocity at several heights and shear stress at the surface made during the Utah field campaign for the presence of ranges of scales, where distinct scale-to-scale interactions between velocity and shear stress can be identified. We find that our results are similar to those obtained in a previous study [Venugopal et al., 2003] (contrary to the claim in V2003, that the scaling relations might be dependent on Reynolds number) where wind tunnel measurements of velocity and shear stress were analysed. We use a wavelet-based scale-to-scale cross-correlation to detect three ranges of scales of interaction between velocity and shear stress, namely, (a) inertial subrange, where the correlation is negligible; (b) energy production range, where the correlation follows a logarithmic law; and (c) for scales larger than the boundary layer height, the correlation reaches a plateau.

Cumulus clouds and convective boundary layers: a tropical perspective on two turbulent shear flows

This paper is a review prepared for the second Marseille Colloquium on the mechanics of turbulence, held in 2011, 50 years after the first. The review covers recent developments in our understanding of the large-scale dynamics of cumulus cloud flows and of the atmospheric boundary layer in the low-wind convective regime that is often encountered in the tropics. It has recently been shown that a variety of cumulus cloud forms and life cycles can be experimentally realized in the laboratory, with the transient diabatic plume taken as the flow model for a cumulus cloud. The plume is subjected to diabatic heating scaled to be dynamically similar to heat release from phase changes in clouds. The experiments are complemented by exact numerical solutions of the Navier-Stokes-Boussinesq equations for plumes with scaled off-source heating. The results show that the Taylor entrainment coefficient first increases with heating, reaches a positive maximum and then drops rapidly to zero or even negative values. This reduction in entrainment is a consequence of structural changes in the flow, smoothing out the convoluted boundaries in the non-diabatic plume, including the tongues engulfing the ambient flow. This is accompanied by a greater degree of mixedness in the core flow because of lower dilution by the ambient fluid. The cloud forms generated depend strongly on the history of the diabatic heating profile in the vertical direction. The striking effects of heating on the flow are attributable to the operation of the baroclinic torque due to the temperature field. The mean baroclinic torque is shown to peak around a quasi-cylindrical sheet situated midway between the axis of the flow and the edges. This torque is shear-enhancing and folds down the engulfment tongues. The increase in mixedness can be traced to an explosive growth in the enstrophy, triggered by a strong fluctuating baroclinic torque that acts as a source, especially at the higher wave numbers, thus enhancing the mixedness. In convective boundary layers field measurements show that, under conditions prevailing in the tropics, the eddy fluxes of momentum and energy do not follow the Monin-Obukhov similarity. Instead, the eddy momentum flux is found to be linear in the wind speed at low winds; and the eddy heat flux is, to a first approximation, governed by free convection laws, with wind acting as a small perturbation on a regime of free convection. A new boundary layer code, based on heat flux scaling rather than wall-stress scaling, shows promising improvements in predictive skills of a general circulation model.

On the k(1)(-1) scaling in sink-flow turbulent boundary layers

Scaling of the streamwise velocity spectrum phi(11)(k(1)) in the so-called sink-flow turbulent boundary layer is investigated in this work. The present experiments show strong evidence for the k(1)(-1) scaling i.e. phi(11)(k(1)) = Lambda(1)U(tau)(2)k(1)(-1), where k(1)(-1) is the streamwise wavenumber and U-tau is the friction velocity. Interestingly, this k(1)(-1) scaling is observed much farther from the wall and at much lower flow Reynolds number (both differing by almost an order of magnitude) than what the expectations from experiments on a zero-pressure-gradient turbulent boundary layer flow would suggest. Furthermore, the coefficient A(1) in the present sink-flow data is seen to be non-universal, i.e. A(1) varies with height from the wall; the scaling exponent -1 remains universal. Logarithmic variation of the so-called longitudinal structure function, which is the physical-space counterpart of spectral k(1)(-1) scaling, is also seen to be non-universal, consistent with the non-universality of A(1). These observations are to be contrasted with the universal value of A(1) (along with the universal scaling exponent of 1) reported in the literature on zero-pressure-gradient turbulent boundary layers. Theoretical arguments based on dimensional analysis indicate that the presence of a streamwise pressure gradient in sink-flow turbulent boundary layers makes the coefficient A(1) non-universal while leaving the scaling exponent -1 unaffected. This effect of the pressure gradient on the streamwise spectra, as discussed in the present study (experiments as well as theory), is consistent with other recent studies in the literature that are focused on the structural aspects of turbulent boundary layer flows in pressure gradients (Harun etal., J. Flui(d) Mech., vol. 715, 2013, pp. 477-498); the present paper establishes the link between these two. The variability of A(1) accommodated in the present framework serves to clarify the ideas of universality of the k(1)(-1) scaling.

Scaling of pressure spectrum in turbulent boundary layers

Scaling of pressure spectrum in zero-pressure-gradient turbulent boundary layers is discussed. Spatial DNS data of boundary layer at one time instant (Re-theta = 4500) are used for the analysis. It is observed that in the outer regions the pressure spectra tends towards the -7/3 law predicted by Kolmogorov's theory of small-scale turbulence. The slope in the pressure spectra varies from -1 close to the wall to a value close to -7/3 in the outer region. The streamwise velocity spectra also show a -5/3 trend in the outer region of the flow. The exercise carried out to study the amplitude modulation effect of the large scales on the smaller ones in the near-wall region reveals a strong modulation effect for the streamwise velocity, but not for the pressure fluctuations. The skewness of the pressure follows the same trend as the amplitude modulation coefficient, as is the case for the velocity. In the inner region, pressure spectra were seen to collapse better when normalized with the local Reynolds stress (-(u'v') over bar) than when scaled with the local turbulent kinetic energy (q(2) = (u'(2)) over bar + (v'(2)) over bar + (w'(2)) over bar)

Investigation of the separated region ahead of three-dimensional protuberances on plates and cones in hypersonic flows with laminar boundary layers

Heat transfer rate and pressure measurements were made upstream of surface pro-tuberances on a flat plate and a sharp cone subjected to hypersonic flow in a conventional shock tunnel. Heat flux was measured using platinum thin-film sensors deposited on macor substrate and the pressure measurements were made using fast acting piezoelectric sensors. A distinctive hot spot with highest heat flux was obtained near the foot of the protuberance due to heavy vortex activity in the recirculating region. Schlieren flow visualization was used to capture the shock structures and the separation distance ahead of the protrusions was quantitatively measured for varying protuberance heights. A computational analysis was conducted on the flat plate model using commercial computational fluid dynamics software and the obtained trends of heat flux and pressure were compared with the experimental observation. Experiments were also conducted by physically disturbing the laminar boundary layer to check its effect on the magnitude of the hot spot heat flux. In addition to air, argon was also used as test gas so that the Reynolds number can be varied. (C) 2014 AIP Publishing LLC.

PREDICTION OF TRANSITIONAL AND SEPARATED BOUNDARY LAYERS IN A COMPRESSOR CASCADE

Numerical simulations were performed of experiments from a cascade of stator blades at three low Reynolds numbers representative of flight conditions. Solutions were assessed by comparing blade surface pressures, velocity and turbulence intensity along blade normals at several stations along the suction surface and in the wake. At Re = 210,000 and 380,000 the laminar boundary layer over the suction surface separates and reattaches with significant turbulence fluctuations. A new 3-equation transition model, the k-k(L)-omega model, was used to simulate this flow. Predicted locations of the separation bubble, and profiles of velocity and turbulence fluctuations on blade-normal lines at various stations along the blade were found to be quite close to measurements. Suction surface pressure distributions were not as close at the lower Re. The solution with the standard k-omega SST model showed significant differences in all quantities. At Re = 640,000 transition occurs earlier and it is a turbulent boundary layer that separates near the trailing edge. The solution with the Reynolds stress model was found to be quite close to the experiment in the separated region also, unlike the k-omega SST solution. Three-dimensional computations were performed at Re = 380,000 and 640,000. In both cases there were no significant differences between the midspan solution from 3D computations and the 2D solutions. However, the 3D solutions exhibited flow features observed in the experiments the nearly 2D structure of the flow over most of the span at 380,000 and the spanwise growth of corner vortices from the endwall at 640,000.

Instabilities in unsteady boundary layers with reverse flow

Instabilities arising in unsteady boundary layers with reverse flow have been investigated experimentally. Experiments are conducted in a piston driven unsteady water tunnel with a shallow angle diffuser placed in the test section. The ratio of temporal (Pi(t)) to spatial (Pi(x)) component of the pressure gradient can be varied by a controlled motion of the piston. In all the experiments, the piston velocity variation with time is trapezoidal consisting of three phases: constant acceleration from rest, constant velocity and constant deceleration to rest. The adverse pressure gradient (and reverse flow) are due to a combination of spatial deceleration of the free stream in the diffuser and temporal deceleration of the free stream caused by the piston deceleration. The instability is usually initiated with the formation of one or more vortices. The onset of reverse flow in the boundary layer, location and time of formation of the first vortex and the subsequent flow evolution are studied for various values of the ratio Pi(x) (Pi(x) + Pi(t)) for the bottom and the top walls. Instability is due to the inflectional velocity profiles of the unsteady boundary layer. The instability is localized and spreads to the other regions at later times. At higher Reynolds numbers growth rate of instability is higher and localized transition to turbulence is observed. Scalings have been proposed for initial vortex formation time and wavelength of the instability vortices. Initial vortex formation time scales with convective time, delta/Delta U, where S is the boundary layer thickness and Delta U is the difference of maximum and minimum velocities in the boundary layer. Non-dimensional vortex formation time based on convective time scale for the bottom and the top walls are found to be 23 and 30 respectively. Wavelength of instability vortices scales with the time averaged boundary layer thickness. (C) 2015 Elsevier Masson SAS. All rights reserved.