Large eddy simulation (2D) using diffusion–velocity method and vortex‐in‐cell

A large eddy Simulation based on the diffusion‐velocity method and the discrete vortex method is presented. The vorticity‐based and eddy viscosity type subgrid scale model simulating the enstrophy transfer between the large and small scale appears as a convective term in the diffusion‐velocity formulation. The methodology has been tested on a spatially growing mixing layer using the two‐dimensional vortex‐in‐cell method and the Smagorinsky subgrid scale model. The effects on the vorticity contours, momemtum thickness, mean streamwise velocity profiles, root‐mean‐square velocity and vorticity fluctuations and negative cross‐stream correlation are discussed. Comparison is made with experiment and numerical work where diffusion is simulated using random walk. Copyright © 2004 John Wiley & Sons, Ltd.     

Near-wall hot-wire measurements

This work continues the studies of Khoo et al. (Exp. Fluids 29: 448–460, 2001), where experiments were performed in turbulent-channel and flat-plate boundary-layer flows using near-wall hot-wire probes. The probability density function (pdf) of the wall-shear stress and streamwise velocity fluctuations in the viscous sublayer, buffer region and beyond were compared and analyzed. The convective velocity U _c of the streamwise velocity fluctuations in the very near-wall region was obtained using a two-point correlation technique. It was found that in the viscous sublayer, U _c is approximately constant at 13u _τ and 15u _τ , respectively, for the channel and boundary-layer flows. Spectra data for the viscous sublayer are presented for the first time, and the normalized spectral plots for different flow conditions collapse at high frequencies or wavenumbers, thus indicating the possible presence of small-scale universality at different Reynolds numbers. The integral time scale corresponding to the streamwise velocity fluctuations in the viscous sublayer is also presented. Received: 18 October 2000/Accepted: 2 April 2001     

Large eddy simulation (2D) of spatially developing mixing layer using vortex‐in‐cell for flow field and filtered probability density function for scalar field

A large eddy simulation based on filtered vorticity transport equation has been coupled with filtered probability density function transport equation for scalar field, to predict the velocity and passive scalar fields. The filtered vorticity transport has been formulated using diffusion‐velocity method and then solved using the vortex method. The methodology has been tested on a spatially growing mixing layer using the two‐dimensional vortex‐in‐cell method in conjunction with both Smagorinsky and dynamic eddy viscosity subgrid scale models for an anisotropic flow. The transport equation for filtered probability density function is solved using the Lagrangian Monte‐Carlo method. The unresolved subgrid scale convective term in filtered density function transport is modelled using the gradient diffusion model. The unresolved subgrid scale mixing term is modelled using the modified Curl model. The effects of subgrid scale models on the vorticity contours, mean streamwise velocity profiles, root‐mean‐square velocity and vorticity fluctuations profiles and negative cross‐stream correlations are discussed. Also the characteristics of the passive scalar, i.e. mean concentration profiles, root‐mean‐square concentration fluctuations profiles and filtered probability density function are presented and compared with previous experimental and numerical works. The sensitivity of the results to the Schmidt number, constant in mixing frequency and inflow boundary conditions are discussed. Copyright © 2005 John Wiley & Sons, Ltd.     

Turbulent flow in a channel at a low Reynolds number

A two-component laser Doppler velocimeter with high spatial and temporal resolution was used to obtain measurements for fully developed turbulent flow of water through a channel with an aspect ratio of 12 : 1 at Re=5700 (based on the centerline velocity and the half-height of the channel). Statistical quantities that were determined are the mean streamwise velocity, the root-mean-square of the fluctuations of the streamwise and the normal velocities, the Reynolds shear stress and higher order moments. Turbulence production is calculated from these quantities. Turbulence statistics obtained from experiments are compared with results from a direct numerical simulation at the same Reynolds number. The good agreement validates a recent DNS, at Re=5700, which is approximately twice as large as used in most previous studies. Received: 12 May 1997 / Accepted: 8 April 1998     

A DPIV study of a starting flow downstream of a backward-facing step

 In this paper an experimental investigation of a starting water flow downstream of a backward-facing step is described. The Reynolds number of the asymptotic steady flow is Re≈4300 based on the step height of s=2 cm and the free stream velocity of U=21.4 cm/s. Velocity measurements were performed with video-based DPIV (Digital Particle Image Velocimetry) at a sampling frequency of 25 Hz. The main purpose of this study is to reveal the temporal development of global structures which could not be analyzed with single-point probes. It was found that at initialization of the flow a regular vorticity street is formed, which collapses at a normalized time t ^* =U t/s≈17 due to vorticity interactions. After this the flow is dominated by complicated vorticity roll-up and shedding dynamics in the recirculation region. The starting phase is terminated for t ^* >40. Prior to the collapse of the vorticity street values of 9 times the steady state asymptotic wall normal stress and of twice the steady state negative wall shear stress were observed. The early increasing slope of the reattachment length is constant over a time of approximately t ^* =8. The collapse of the vorticity street and the vorticity interactions thereafter cause fluctuations both in the velocity in the free shear layer and in the reattachment length. The fully developed flow has a dominant frequency corresponding to a Strouhal number St=fs/U≈0.04. Received: 20 September 1996/Accepted: 1 April 1997     

Reynolds-number-dependence of the maximum in the streamwise velocity fluctuations in wall turbulence

A survey is made of the standard deviation of the streamwise velocity fluctuations in near-wall turbulence and in particular of the Reynolds-number-dependency of its peak value. The following canonical flow geometries are considered: an incompressible turbulent boundary layer under zero pressure gradient, a fully developed two-dimensional channel and a cylindrical pipe flow. Data were collected from 47 independent experimental and numerical studies, which cover a Reynolds number range of R _θ=U _∞ θ/v=300−20,920 for the boundary layer with θ the momentum thickness and R ⁺=u _*R/v=100-4,300 for the internal flows with R the pipe radius or the channel half-width. It is found that the peak value of the rms-value normalised by the friction velocity, u _*, is within statistical errors independent of the Reynolds number. The most probable value for this parameter was found to be 2.71±0.14 and 2.70±0.09 for the case of a boundary layer and an internal flow, respectively. The present survey also includes some data of the streamwise velocity fluctuations measured over a riblet surface. We find no significant difference in magnitude of the normalised peak value between the riblet and smooth surfaces and this property of the normalised peak value may for instance be exploited to estimate the wall shear stress from the streamwise velocity fluctuations. We also consider the skewness of the streamwise velocity fluctuations and find its value to be close to zero at the position where the variance has its peak value. This is explained with help of the equations of the third-order moment of velocity fluctuations. These results for the peak value of the rms of the streamwise velocity fluctuations and also the coincidence of this peak with the zero value of the third moment can be interpreted as confirmation of local equilibrium in the near-wall layer, which is the basis of inner-layer scaling. Furthermore, these results can be also used as a requirement which turbulence models for the second and triple velocity correlations should satisfy. The authors are indebted to Prof. P. Bradshaw for making available his list of references on this topic and for his remarks on “active” and “inactive” motions. We also gratefully acknowledge discussions with Prof. I. Castro regarding the value of σ _u ⁺ above rough walls.     

Streamwise evolution of an inclined cylinder wake

The streamwise evolution of an inclined circular cylinder wake was investigated by measuring all three velocity and vorticity components using an eight-hotwire vorticity probe in a wind tunnel at a Reynolds number Re_d of 7,200 based on free stream velocity (U _∞) and cylinder diameter (d). The measurements were conducted at four different inclination angles (α), namely 0°, 15°, 30°, and 45° and at three downstream locations, i.e., x/d = 10, 20, and 40 from the cylinder. At x/d = 10, the effects of α on the three coherent vorticity components are negligibly small for α ≤ 15°. When α increases further to 45°, the maximum of coherent spanwise vorticity reduces by about 50%, while that of the streamwise vorticity increases by about 70%. Similar results are found at x/d = 20, indicating the impaired spanwise vortices and the enhancement of the three-dimensionality of the wake with increasing α. The streamwise decay rate of the coherent spanwise vorticity is smaller for a larger α. This is because the streamwise spacing between the spanwise vortices is bigger for a larger α, resulting in a weak interaction between the vortices and hence slower decaying rate in the streamwise direction. For all tested α, the coherent contribution to [`(v<sup>2</sup>)] \overline{{v^{2}}} is remarkable at x/d = 10 and 20 and significantly larger than that to [`(u²)] \overline{{u^{2}}} and [`(w<sup>2</sup>)]. \overline{{w^{2}}}. This contribution to all three Reynolds normal stresses becomes negligibly small at x/d = 40. The coherent contribution to [`(u²)] \overline{{u^{2}}} and [`(v²)] \overline{{v^{2}}} decays slower as moving downstream for a larger α, consistent with the slow decay of the coherent spanwise vorticity for a larger α.     

On the friction drag reduction mechanism of streamwise wall fluctuations

Understanding how to decrease the friction drag exerted by a fluid on a solid surface is becoming increasingly important to address key societal challenges, such as decreasing the carbon footprint of transport. Well-established techniques are not yet available for friction drag reduction. Direct numerical simulation results obtained by Józsa et al. (2019) previously indicated that a passive compliant wall can decrease friction drag by sustaining the drag reduction mechanism of an active control strategy. The proposed compliant wall is driven by wall shear stress fluctuations and responds with streamwise wall velocity fluctuations. The present study aims to clarify the underlying physical mechanism enabling the drag reduction of these active and passive control techniques. Analysis of turbulence statistics and flow fields reveals that both compliant wall and active control amplify streamwise velocity streaks in the viscous sublayer. By doing so, these control methods counteract dominant spanwise vorticity fluctuations in the near-wall region. The lowered vorticity fluctuations lead to an overall weakening of vortical structures which then mitigates momentum transfer and results in lower friction drag. These results might underpin the further development and practical implementation of these control strategies.     

Experimental measurement of flow past cavities of different shapes

Experiments are carried out in order to investigate the flow structure past a rectangular, triangular and semi-circular cavity of length-to-depth ratio of 2:1 using the Particle Image Velocimetry (PIV) technique. The experiments are performed in a large scale water channel with three different upstream velocities resulting in Reynolds numbers of 1230, 1460 and 1700, based on inflow momentum thickness, for each cavity type. Contours of constant averaged streamwise and transverse components of velocity, contours of constant averaged vorticity, Reynolds stress and streamline plots for each cavity type for the aforementioned three Reynolds numbers are presented. In addition, streamwise velocity, Reynolds stress and turbulence intensity values are compared for all cavity types. Effect of cavity shape on flow structure within the cavity is discussed in detail. Moreover, spectrum of instantaneous streamwise velocity fluctuations in shear layer near the downstream of the leading corner and the upstream of the trailing corner of the cavities are obtained and it was found that no organized oscillations are present in the flow; rectangular and triangular cavities have the largest amplitudes while semi-circular cavity has the smallest. Calculated turbulence intensities also reveal that the maximum turbulence intensities occur at cavity lid in the centerline section and rectangular and triangular cavities have larger turbulence intensity compared to semi-circular cavity.     

Construction of three-dimensional images of flow structure via particle tracking techniques

Construction of three-dimensional images of flow structure, based on the quantitative velocity field, is assessed for cases where experimental data are obtained using particle tracking technique. The experimental data are in the form of contiguous planes of particle images. These contiguous data planes are assumed to correspond to successive spatial realizations in steady flow, or to phase-referenced realizations in an unsteady flow.Given the particle images on contiguous planes, the in-plane velocity fields are determined. Then, the out-of-plane velocity field is obtained using a spectral interpolation method. Application of this method allows, in principle, construction of the three-dimensional vorticity field and the streamline patterns.A critical assessment is made of the uncertainties arising from the in-plane interpolation of the velocity field obtained from particle tracking and the evaluation of the out-of-plane velocity component. The consequences of such uncertainties on the reconstructed vorticity distributions and streamline patterns are addressed for two basic types of vortex flows: a columnar vortex, for which the streamlines are not closed and are spatially periodic in the streamwise direction; and for a spherical (Hill's) vortex exhibiting closed streamline patterns, and no spatial periodicity.     

Breakdown of the Reynolds Analogy in a Stagnation Region Under Inflow Disturbances

A systematic analysis is performed for the Reynolds analogy breakdown at stagnation-region flow and heat transfer in the presence of inflow disturbances. The Reynolds analogy breakdown between momentum and energy transfers in a stagnation region is scrutinized by varying the Reynolds number (5000≤Re≤20000), the amplitude (0.00075≤A≤0.003) and the length scale (λ/δ=10.6). A spanwise sinusoidal variation is superimposed on the velocity component normal to the wall. Self-similarity solutions are obtained with trigonometric series expansions. The Reynolds analogy criterion demonstrates that the rate of change of skin friction is different from that of wall heat transfer. Different evolutions of the rates of skin friction and wall heat transfer are due to the difference between 〈s'v'〉 and 〈v'T'〉. An in-depth analysis on 〈s'v'〉 and 〈 v'T'〉 is performed by analysis using disturbance correlations based on the fluctuating velocity transport equations in vorticity form. It is found that the pressure fluctuations, the wall blocking and the Lamb vectors are responsible for the breakdown of the Reynolds analogy. A direct comparison is made between momentum and energy balances associated with the three responsible mechanisms. A common finding is that their profiles are changed significantly at a location where the evolution of the streamwise vortex is strong. Received 12 May 2000 and accepted 6 March 2001     

Solid/free-surface juncture boundary layer and wake

 The Reynolds-averaged flow for a solid/free-surface juncture boundary layer and wake is documented. The three mean-velocity components and five of the Reynolds stresses are measured for a surface-piercing flat plate in a towing tank using a laser-Doppler velocimeter system for both boundary-layer and wake planes in regions close to the free surface. The experimental method is described, including the foil-plate model, laser-Doppler velocimeter system, conditions, and uncertainty analysis. The underlying flow data is in excellent agreement with benchmark data. Inner (near the plate and wake centerplane and below the free surface) and outer (near the free surface) regions of high streamwise vorticity of opposite sign are observed, which transport, respectively, high mean velocity and low turbulence from the outer to the inner and low mean velocity and high turbulence from the inner to the outer portions of the boundary layer and wake. For the wake, the inner region of vorticity is relatively weak. The physical mechanism for the streamwise vorticity is analyzed with regard to the Reynolds-averaged streamwise vorticity equation. The anisotropy of the crossplane normal Reynolds stresses closely correlates with the vorticity and, additionally, indicates similarity, i.e., its nature is such that it only depends on the proximity to the plate and free surface boundaries or wake centerplane symmetry plane. Free-surface effects on the Reynolds stresses are analyzed with regard to the behavior close to the free surface of the turbulent kinetic energy and the normal components of the anisotropy tensor and the anisotropy invariants. Close to the free surface, the turbulent kinetic energy is nearly constant and increases for the inner and outer portions, respectively, of the boundary layer and wake and the normal components of the anisotropy tensor and the anisotropy invariants roughly correspond to the limiting values for two-component turbulence. The similarities and differences between the present results and analysis with those from related studies are discussed. The data and analysis should have practical application with regard to the development of turbulence models for computational fluid dynamics methods for the Reynolds-averaged Navier–Stokes equations. Received: 27 May 1997/Accepted: 1 August 1997     

Turbulent structure during transition to self-similarity in a round jet

 The developing turbulent region of a round jet was investigated using an improved implementation of digital particle image velocimetry (DPIV). The two-dimensional flow field in planes normal and parallel to the axial velocity was measured at locations between 15 and 30 diameters downstream, for two Reynolds numbers of 5500 and 16,000. The study consisted of instantaneous snapshots of the velocity and vorticity fields as well as measurements of velocity correlations up to third order. In this regime, the Reynolds number had a significant effect on both the instantaneous flow structure and the profiles of mean velocity across the jet. Coherent streamwise structures were present in the jet core for the lower Reynolds number. Additional structures whose evolution was governed by time scales two orders of magnitude larger than the convective scale inside the jet were observed in the entrainment field. The velocity correlations provided further support for the validity of DPIV turbulence measurements. The data was consistent with the equations of motion and momentum was conserved. DPIV measurements of turbulent kinetic energy components agreed with the hot-wire measurements of previous studies. Received: 27 November 1996/Accepted: 14 July 1997     

Experimental study of flow past a square cylinder at high Reynolds numbers

 Measurements of two-components of velocity in the wake of a square cylinder using a hot-wire anemometer are reported. Two Reynolds numbers, namely 8700 and 17,625, have been considered. The measurements were carried out in a low-speed, low-turbulence wind tunnel. Benchmark experiments at much lower Reynolds numbers show good agreement between the present experiments and published results. At higher Reynolds numbers, the experimental data reveal anticipated trends in terms of wake recovery and turbulence decay. Both velocity and velocity fluctuations show symmetry about the wake axis. The experimental data have been compared with the large eddy simulation (LES) calculation reported by Wang et al. [University of Illinois at Urbana – Champaign (1996) Report CFD 96-03] and LDV measurements of Lyn et al. [J Fluid Mech (1995) 304: 285–319]. The agreement among the three sets is generally acceptable in terms of the time-averaged velocity components, but not the velocity fluctuations. The turbulence fluctuations in the present experiments are seen to be lower than in the referred work. The differences have been traced to factors such as the aspect ratio, blockage ratio and upstream turbulence. Experiments with increased upstream turbulence did show a reduction in the discrepancy between the present experiments and the published data. An assessment of the experimental data in terms of physical mechanisms revealed that (a) streamwise normal stresses were correlated with the vortex centers, and (b) the turbulence kinetic energy profiles are similar to the turbulence shear stress. Spectral analysis of the velocity signals was carried out in the present work. Energy transfer from the mean flow to the streamwise velocity fluctuation was confirmed in the near wake. A redistribution of the kinetic energy between the streamwise and transverse components of velocity over a longer distance downstream was subsequently observed. Received: 17 May 1999/Accepted: 29 December 1999     

Experimental investigation of the unsteady structure of a transitional plane wall jet

 A laminar wall jet undergoing transition is investigated using the particle image velocimetry (PIV) technique. The plane wall jet is issued from a rectangular channel, with the jet-exit velocity profile being parabolic. The Reynolds number, based on the exit mean velocity and the channel width, is 1450. To aid the understanding of the global flow features, laser-sheet/smoke flow visualizations are performed along streamwise, spanwise, and cross-stream directions. Surface pressure measurements are made to correlate the instantaneous vorticity distribution with the surface pressure fluctuations. The instantaneous velocity and vorticity field measurements provide the basis for understanding the formation of the inner-region vortex and the subsequent interactions between the outer-region (free-shear-layer region) and inner-region (boundary-layer region) vortical structures. Results show that under the influence of the free-shear-layer vortex, the local boundary layer becomes detached from the surface and inviscidly unstable, and a vortex is formed in the inner region. Once this vortex has formed, the free-shear-layer vortex and the inner-region vortex form a vortex couple and convect downstream. The mutual interactions between these inner- and outer-region vortical structures dominate the transition process. Farther downstream, the emergence of the three-dimensional structure in the free shear layer initiates complete breakdown of the flow. Received: 8 November 1995/Accepted: 6 November 1996     

Modeling of drag reduction in turbulent channel flow with hydrophobic walls by FVM method and weakly-compressible flow equations

In this paper the effects of hydrophobic wall on skin-friction drag in the channel flow are investigated through large eddy simulation on the basis of weaklycompressible flow equations with the MacCormack's scheme on collocated mesh in the FVM framework. The slip length model is adopted to describe the behavior of the slip velocities in the streamwise and spanwise directions at the interface between the hydrophobic wall and turbulent channel flow. Simulation results are presented by analyzing flow behaviors over hydrophobic wall with the Smagorinky subgrid-scale model and a dynamic model on computational meshes of different resolutions. Comparison and analysis are made on the distributions of timeaveraged velocity, velocity fluctuations, Reynolds stress as well as the skin-friction drag. Excellent agreement between the present study and previous results demonstrates the accuracy of the simple classical second-order scheme in representing turbulent vertox near hydrophobic wall. In addition, the relation of drag reduction efficiency versus time-averaged slip velocity is established. It is also foundthat the decrease of velocity gradient in the close wall region is responsible for the drag reduction. Considering its advantages of high calculation precision and efficiency, the present method has good prospect in its application to practical projects.     

The influence of temperature on the measurements of Reynolds stresses in shear free turbulence near a wall

 Temperature changes have a significant influence on the measurements of Reynolds stresses in turbulent boundary layers. As compared to the spanwise velocity fluctuations the streamwise turbulence intensity is especially sensitive to temperature deviations. Although this is a general statement its importance is clearly elucidated in a shear-free turbulence near a solid wall, since the mixing due to turbulence production is minimized in this flow. A consequence of temperature influence on hot-wire measurements is that frictional heating from the wall has produced contradictory results in different experiments on shear-free turbulence. In the current paper, measurements of streamwise and spanwise turbulence intensities have been conducted at different wall temperatures, thereby simulating the contradictory results mentioned above. A simple model has been developed showing that the turbulence intensities are affected by both the rms. value of the temperature fluctuations and the correlation between fluctuating temperature and velocity. These correlations are measured and the developed model is used to explain deviations in earlier measurements on shear-free turbulence. Moreover, the individual magnitudes of the two correlations in the temperature correction are estimated and their individual importance is discussed. Received: 17 February 1997 / Accepted: 18 January 1998     

The compressible mixing layer: an LES study

This article employs LES to simulate temporal mixing layers with Mach numbers ranging from M _c  = 0.3 to M _c  = 1.2. A form of approximate deconvolution together with a dynamic Smagorinsky subgrid model are employed as subgrid models. A large computational domain is used along with relatively good resolution. The LES results regarding growth rate, turbulence levels, turbulence anisotropy, and pressure–strain correlation show excellent agreement with those available from previous experimental and DNS results of the same flow configuration, underlining the effectiveness and accuracy of properly conducted LES. Coherent structures during the transitional stage change from spanwise aligned rollers to streamwise-aligned thinner vortices at high Mach number. In the quasi-self-similar turbulent stage, the resolved-scale vorticity is more isotropic at higher M _c , and its vertical correlation length scale is smaller. The ratio of the vertical integral length scale of streamwise velocity fluctuation to a characteristic isotropic estimate is found to decrease with increasing M _c . Thus, compressibility leads to increased spatial decorrelation of turbulence which is one reason for the reduction in pressure–strain correlation with increasing M _c . The balance of the resolved-scale fluctuating vorticity is examined, and it is observed that the linear production by mean shear becomes less important compared to nonlinear vortex stretching at high M _c . A spectral decomposition of the pressure fluctuations into low- and intermediate-to-high-wave numbers is performed. The low-wave number part of the pressure field is found not to correlate with the strain field, although it does have a significant contribution to the r.m.s of the fluctuating pressure. As a consequence, the pressure–strain correlation can be analyzed using a simplified Green’s function for the Poisson equation as is demonstrated here using the LES data.     

Study of Density Effects in Turbulent Buoyant Jets Using Large-Eddy Simulation

Large-Eddy simulations (LES) of spatially evolving turbulent buoyant round jets have been carried out with two different density ratios. The numerical method used is based on a low-Mach-number version of the Navier–Stokes equations for weakly compressible flow using a second-order centre-difference scheme for spatial discretization in Cartesian coordinates and an Adams–Bashforth scheme for temporal discretization. The simulations reproduce the typical temporal and spatial development of turbulent buoyant jets. The near-field dynamic phenomenon of puffing associated with the formation of large vortex structures near the plume base with a varicose mode of instability and the far-field random motions of small-scale eddies are well captured. The pulsation frequencies of the buoyant plumes compare reasonably well with the experimental results of Cetegen (1997) under different density ratios, and the underlying mechanism of the pulsation instability is analysed by examining the vorticity transport equation where it is found that the baroclinic torque, buoyancy force and volumetric expansion are the dominant terms. The roll-up of the vortices is broken down by a secondary instability mechanism which leads to strong turbulent mixing and a subsequent jet spreading. The transition from laminar to turbulence occurs at around four diameters when random disturbances with a 5% level of forcing are imposed to a top-hat velocity profile at the inflow plane and the transition from jet-like to plume-like behaviour occurs further downstream. The energy-spectrum for the temperature fluctuations show both −5/3 and −3 power laws, characteristic of buoyancy-dominated flows. Comparisons are conducted between LES results and experimental measurements, and good agreement has been achieved for the mean and turbulence quantities. The decay of the centreline mean velocity is proportional to x ^−1/3 in the plume-like region consistent with the experimental observation, but is different from the x ⁻¹ law for a non-buoyant jet, where x is the streamwise location. The distributions of the mean velocity, temperature and their fluctuations in the near-field strongly depend upon the ratio of the ambient density to plume density ρ_a/ρ₀. The increase of ρ_a/ρ₀ under buoyancy forcing causes an increase in the self-similar turbulent intensities and turbulent fluxes and an increase in the spatial growth rate. Budgets of the mean momentum, energy, temperature variance and turbulent kinetic energy are analysed and it is found that the production of turbulence kinetic energy by buoyancy relative to the production by shear is increased with the increase of ρ_a/ρ₀. Received 16 June 2000 and accepted 26 June 2001     

Evolution of isolated streamwise vortices in the late stages of boundary layer transition

Experiments were conducted with two, smooth hills, lying well within the boundary layer over a flat plate mounted in a wind tunnel. One hill was shallow, with peak height 1.5 mm and width 50 mm; the other, steep, 3 mm high and 30 mm wide. Since the hills occupied one-half of the tunnel span, streamwise vorticity formed near the hills’ edge. At a freestream speed of 3.5 m/s, streaks formed with inflectional wall-normal and spanwise velocity profiles but without effecting transition. Transition, observed at 7.5 m/s, took different routes with the two hills. With the steep hill, streamwise velocity signals exhibited the passage of a wave packet which intensified before breakdown to turbulence. With the shallow hill there was a broad range of frequencies present immediately downstream of the hill. These fluctuations grew continuously and transition occurred within a shorter distance. Since the size of the streamwise vorticity generated at the hill edge is of the order of the hill height, the shallow hill generates vorticity closer to the wall and supports an earlier transition, whereas the steep hill creates a thicker vortex and associated streaks which exhibit oscillations due to their own instability as an additional precursor stage before transition.