Structure of the turbulent boundary layer over a rod-roughened wall

Turbulent coherent structures near a rod-roughened wall are scrutinized by analyzing instantaneous flow fields obtained from direct numerical simulations (DNSs) of a turbulent boundary layer (TBL). The roughness elements used are periodically arranged two-dimensional spanwise rods, and the roughness height is k/δ = 0.05 where δ is the boundary layer thickness. The Reynolds number based on the momentum thickness is varied in the range Re_θ = 300–1400. The effect of surface roughness is examined by comparing the characteristics of the TBLs over smooth and rough walls. Although introduction of roughness elements onto the smooth wall affects the Reynolds stresses throughout the entire boundary layer when scaled by the friction velocity, the roughness has little effect on the vorticity fluctuations in the outer layer. Pressure-strain tensors of the transport equation for the Reynolds stresses and quadrant analysis disclose that the redistribution of turbulent kinetic energy of the rough wall is similar to that of the smooth wall, and that the roughness has little effect on the relative contributions of ejection and sweep motions in the outer layer. To elucidate the modifications of the near-wall vortical structure induced by surface roughness, we used two-point correlations, joint weighted probability density function, and linear stochastic estimation. Finally, we demonstrate the existence of coherent structures in the instantaneous flow field over the rod-roughened surface.     

Effect of Roughness on Wall-Bounded Turbulence

Direct numerical simulation of turbulent incompressible plane-channel flow between a smooth wall and one covered with regular three-dimensional roughness elements is performed. While the impact of roughness on the mean-velocity profile of turbulent wall layers is well understood, at least qualitatively, the manner in which other features are affected, especially in the outer layer, has been more controversial. We compare results from the smooth- and rough-wall sides of the channel for three different roughness heights of h ⁺= 5.4, 10.8, and 21.6 for Re _τ of 400, to isolate the effects of the roughness on turbulent statistics and the instantaneous turbulence structure at large and small scales. We focus on the interaction between the near-wall and outer-layer regions, in particular the extent to which the near-wall behavior influences the flow further away from the surface. Roughness tends to increase the intensity of the velocity and vorticity fluctuations in the inner layer. In the outer layer, although the roughness alters the velocity fluctuations, the vorticity fluctuations are relatively unaffected. The higher-order moments and the energy budgets demonstrate significant differences between the smooth-wall and rough-wall sides in the processes associated with the wall-normal fluxes of the Reynolds shear stresses and turbulence kinetic energy. The length scales and flow dynamics in the roughness sublayer, the spatially inhomogeneous layer within which the flow is directly influenced by the individual roughness elements, are also examined. Alternative mechanisms involved in producing and maintaining near-wall turbulence in rough-wall boundary layers are also considered. We find that the strength of the inner/outer-layer interactions are greatly affected by the size of the roughness elements.     

Turbulent boundary layer flow with a step change from smooth to rough surface

A direct numerical simulation (DNS) dataset of a turbulent boundary layer (TBL) with a step change from a smooth to a rough surface is analyzed to examine the characteristics of a spatially developing flow. The roughness elements are periodically arranged two-dimensional (2-D) spanwise rods, with the first rod placed 80θ_in downstream from the inlet, where θ_in denotes the inlet momentum thickness. Based on an accurate estimation of relevant parameters, clear evidence for mean flow universality is provided when scaled properly, even for the present roughness configuration, which is believed to have one of the strongest impacts on the flow. Compared to previous studies, it is shown that overshooting behavior is present in the first- and second-order statistics and is locally created either within the cavity or at the leading edge of the roughness depending on the type of statistics and the wall-normal measurement location. Inspection of spatial two-point correlations of the streamwise velocity fluctuations shows a continuous increase of spanwise length scales of structures over the rough wall after the step change at a greater growth rate than that over smooth wall TBL flow. This is expected because spanwise energy spectrum shows presence of much energetic wider structures over the rough wall. Full images of the DNS data are presented to describe not only predominance of hairpin vortices but also a possible spanwise scale growth mechanism via merging over the rough wall.     

From Streaks to Spots and on to Turbulence: Exploring the Dynamics of Boundary Layer Transition

“...an eerie type of chaos can lurk just behind a facade of order, and yet deep inside the chaos lurks an even eerier type of order.” Douglas Hofstadter Bypass transition to turbulence in boundary layers is examined using linear theory and direct numerical simulations (DNS). First, the penetration of low-frequency free-stream disturbances into the boundary layer is explained using a model problem with two time scales, namely the shear and wall-normal diffusion. The simple model provides a physical understanding of the phenomenon of shear sheltering. The second stage in bypass transition is the amplification of streaks. Streak detection and tracking algorithms were applied to examine the characteristics of the streak population inside the boundary layer, beneath free-stream turbulence. It is demonstrated that simple statistical averaging masks the wealth of streak amplitudes in transitional flows, and in particular the high-amplitude, relatively rare events that precede the onset of turbulence. The third stage of the transition process, namely the secondary instability of streaks, is examined using secondary instability analysis. It is demonstrated that two types of instability are possible: An outer instability arises near the edge of the boundary layer on the lifted, low-speed streaks. An inner instability also exists, and has the appearance of a near-wall wavepacket. The stability theory is robust, and can predict the particular streaks which are likely to undergo secondary instability and break down in transitional boundary layers beneath free-stream turbulence. Beyond the secondary instability, turbulent spots are tracked in DNS in order to examine their characteristics in the subsequent non-linear stages of transition. At every stage, we compare the findings from linear theory to the empirical observations from direct solutions of the Navier-Stokes equations. The complementarity between the theoretical predictions and the computational experiments is highlighted, and it leads to a detailed view of the mechanics of transition.     

An experimental-based model for prediction of flow noise with drag reducing polymers

The drag reduction characteristics of certain high molecular weight polymers have been studied by various investigators. Because of the polymer’s ability to reduce turbulent shear stress and dependence of the boundary layer wall pressure spectral amplitude on the shear stress, polymer has the potential to suppress noise and vibration caused by the boundary layer unsteady pressures. Compared to its effect on drag reduction, polymer additive effects on turbulent boundary layer (TBL) wall pressure fluctuations have received little attention. Kadykov and Lyamshev [Sov. Phys. Acoust. 16 (1970) 59], Greshilor et al. [Sov. Phys. Acoust. 21 (1975) 247] showed that drag reducing polymer additives do indeed reduce wall pressure fluctuations, but they have not established any scaling relationship which effectively collapse data. Some effort has been made by Timothy et al. [JASA 108 (1) (2000) 71] at Penn State University to develop a scaling relationship for TBL wall pressure fluctuations that are modified by adding drag reducing polymer to pure water flow. This paper presents a theoretical model based on the work of the Timothy et al. team at ARL, Penn State University. Through this model one can estimate, reduction in TBL flow induced noise and vibration for rigid smooth surfaces due to release of drag reducing polymers in boundary layer region. Using this theoretical model, flow noise as experienced by a typical flush mounted hydrophone has been estimated for a smooth wall plate as a function of polymer additive concentration. Effect of non-dimensionalisation of the wall pressure fluctuations frequency spectra with traditional outer, inner and mixed flow variables will also be addressed in the paper. The paper also covers a model based on molecular relaxation time in polymer additives which not only reduce drag but also flow induced noise up to certain polymer concentration.     

Large Eddy Simulations in Turbines: Influence of Roughness and Free-Stream Turbulence

The individual and coupled effects of the incoming free-stream turbulence (FST) and surface roughness on the transition of a separated shear layer over a flat plate is numerically investigated using Large eddy simulation (LES). The upper wall of the test section is inviscid and specifically contoured to impose a streamwise pressure distribution over the flat plate to simulate the suction surface of a low pressure turbine (LPT) blade. The interaction of the streamwise streaks due to FST and roughness with the separated shear layer is captured. The streaks induced by FST are intermittent in nature and the streaks due to roughness are steady for a given topology of the rough surface. Both FST and roughness promoted near-wall mixing. The ‘net effect’ of mixing in the pre-separated region is manifested by a shift the inflection point of the velocity profile towards the wall. This resulted in the upstream shift of the transition point and a significant reduction in the size of the separation bubble. The combined effect of FST and roughness further reduced the size of separation bubble. The streamwise evolution of the boundary layer parameters has been compared for different cases. The potential ‘roughness benefit’ obtained in the case of highly loaded turbine blades in terms of its considerable reduction of profile loss is also shown.     

Direct Numerical Simulation of an Accelerating Channel Flow

Direct Numerical Simulation of a linearly accelerating channel flow starting from an initially statistically steady turbulent flow has been performed. It is shown that the response of the accelerating flow is fundamentally the same as that of the step-change transient flow described in He and Seddighi (J Fluid Mech 715:60–102, 2013). The flow structure again behaves like a boundary layer bypass transition undergoing three distinct phases, namely, (i) initially (pre-transition), the flow is laminar-like and the pre-existing turbulent structures are modulated resulting in elongated streaks leading to a strong and continuous increase in the streamwise fluctuating velocity but little changes in the other two components; (ii) it then undergoes transition when isolated turbulent spots are generated which spread and merge with each other, and (iii) they eventually cover the entire surface of the wall when the flow is fully turbulent. The similarity between the turbulence responses in the two flows is significant noting the contrasting features of the two types of mean flow unsteadiness: in the step-change flow, a sharp boundary layer is resulted in nearly instantly on the wall which closely resembles the spatially developing boundary layer, whereas the linear flow acceleration causes a continuing change of velocity gradient adjacent to the wall which propagates into the flow field with time, resulting in a gradually-developing boundary layer. There are, however, quantitative differences in the detailed behavior of the two flows and especially the transition is much delayed in the accelerating flow. It is also shown that the late pre-transition and early transition stages in both flows are characterised by significantly increased inwards sweep events in the wall region and ejection events in the outer layer. The flatness of the wall-normal velocity increases markedly near the wall around the time of onset of transition as a consequence of the huge intermittency of the velocity fluctuations. That is, there are long periods of quiescent flow coupled with occasional turbulent bursts.     

Effects of two-dimensional V-shaped grooves on turbulent channel flow

Detailed Laser Doppler velocimeter (LDV) measurements have been carried out in a turbulent rectangular channel flow with one rough wall. The roughness elements of two-dimensional spanwise 120° V-shaped grooves are periodically arranged with different depths and pitches. The Reynolds number based on the centerline velocity, and the channel half height ranges from 2,740 to 20,000. The comparisons of turbulence statistics over smooth and rough walls indicate that the present roughness leads to a significant change in the turbulence both in the inner and in the outer flow. Particularly, the distribution density of the grooves is a key parameter to evaluate the effect of roughness. The low-Reynolds-number dependence of turbulence statistics is also observed. The rough walls with the same pitch-to-depth ratio exhibit the equivalent roughness function under the corresponding Reynolds numbers. The disagreement of velocity defect profiles between smooth and rough walls challenges the defect universal law. The variations of the turbulence stresses and Reynolds shear stress decomposition in the outer layer suggest that the turbulent motions may be modified by the present grooves. The importance of sweep events for the present groove-roughened walls is reflected by the differences in relative contribution to Reynolds shear stress from each quadrant and the higher-order moments over smooth and rough walls.     

Large eddy simulation on the effect of free-stream turbulence on bypass transition

The effect of free-stream turbulence (FST) on bypass transition in a zero-pressure-gradient boundary layer is investigated by means of Large Eddy Simulation (LES). The broadband turbulent inflow is synthesized to validate the feasibility of LES. Both a zero-thickness plate and one with super-ellipse leading-edge are addressed. The calculated Reynolds-averaged fields are compared with experimental data and decent agreement is achieved. Instantaneous fields show the instability occurs in the lifted low-speed streaks similar to earlier DNS results, which can be ascribed to outer mode. Various inflows with bi-/tri-mode interaction are specified to analyze effects of particular frequency mode on the instability pattern and multifarious transition or non-transition scenarios are obtained. Outer instability is observed in the cases with one low-frequency mode and one high-frequency mode inflow as reported by Zaki and Durbin (2005), and with one more high-frequency mode appended. Inner instability is observed in the case with a low-frequency dominant inflow, while the high-frequency mode is indispensable to induce the secondary instability. Furthermore, the results show that the transition onset is highly sensitive to low-frequency mode while the transition rate is highly sensitive to high-frequency mode. Finally, the formational frequency of turbulent spot (FFTS) is counted and the frequency of laminar streaks is demonstrated by spectral analysis.     

Adverse pressure gradient turbulent flows over rough walls

An experimental study of a fully developed turbulent channel flow and an adverse pressure gradient (APG) turbulent channel flow over smooth and rough walls has been performed using a particle image velocimetry (PIV) technique. The rough walls comprised two-dimensional square ribs of nominal height, k = 3 mm and pitch, p = 2k, 4k and 8k. It was observed that rib roughness enhanced the drag characteristics, and the degree of enhancement increased with increasing pitch. Similarly, rib roughness significantly increased the level of turbulence production, Reynolds stresses and wall-normal transport of turbulence kinetic energy and Reynolds shear stress well beyond the roughness sublayer. On the contrary, the distributions of the eddy viscosity, mixing length and streamwise transport of turbulence kinetic energy and Reynolds shear stress were reduced by wall roughness, especially in the outer layer. Adverse pressure gradient produced a further reduction in the mean velocity (in comparison to the results obtained in the parallel section) but increased the wall-normal extent across which the mean flow above the ribs is spatially inhomogeneous in the streamwise direction. APG also reinforced wall roughness in augmenting the equivalent sand grain roughness height. The combination of wall roughness and APG significantly increased turbulence production and Reynolds stresses except in the immediate vicinity of the rough walls. The transport velocities of the turbulence kinetic energy and Reynolds shear stress were also augmented by APG across most part of the rough-wall boundary layer. Further, APG enhanced the distributions of the eddy viscosity across most of the boundary layer but reduced the mixing length outside the roughness sublayer.     

A plane turbulent wall jet on a fully rough surface

The mean velocity field and skin friction characteristics of a plane turbulent wall jet on a smooth and a fully rough surface were studied using Particle Image Velocimetry. The Reynolds number based on the slot height and the exit velocity of the jet was Re = 13,400 and the nominal size of the roughness was k = 0.44 mm. For this Reynolds number and size of roughness element, the flow was in the fully rough regime. The surface roughness results in a distinct change in the shape of the mean velocity profile when scaled in outer coordinates, i.e. using the maximum velocity and outer half-width as the relevant velocity and length scales, respectively. Using inner coordinates, the mean velocity in the lower region of the inner layer was consistent with a logarithmic profile which characterizes the overlap region of a turbulent boundary layer; for the rough wall case, the velocity profile was shifted downward due to the enhanced wall shear stress. For the fully rough flow, the decay rate of the maximum velocity of the wall jet is increased, and the skin friction coefficient is much larger than for the smooth wall case. The inner layer is also thicker for the rough wall case. The effects of surface roughness were observed to penetrate into the outer layer and slightly enhance the spread rate for the outer half-width, which was not observed in most other studies of transitionally rough wall jet flows.     

单一倾斜热线探头在三维湍流场中的应用

研究的目的是确立用单一倾斜热线探头进行三维湍流测试的基础技术．在阐明基本原理，基础方程式的基础上，对比实验结果，证明应用所阐明的方法、用单一倾斜热线探头进行三维湍流测试是切实可行的．对边界层壁面附近的测试结果表明考虑平均速度和变动速度的高次相关项后可以得到更加精确的实验结果．     

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.     

Görtler vortices promoted by wall roughness

Numerical experiments are conducted to investigate spatially developing Görtler vortices and the way in which wall roughness promotes their formation and growth. Several different types of walls are examined and their relative merits as vortex promoters assessed. The only disturbances of the flow are due to the rough wall; hence, at each downstream station the local field feels (1) the upstream flow distribution (produced by the upstream wall conditions) and (2) the local forcing at the wall. Rapid vortex formation and growth, like in the case of ribleted walls, can be qualitatively explained by the positive combination of these two effects; when the two influences on the local flow field compete, e.g. for randomly distributed wall roughness, the equations with the boundary conditions filter the disturbances over some streamwise length, function of the roughness amplitude, to create coherent patches of vorticity out of the random noise. These patches can then be amplified by the instability mechanism. If a thin rough strip is aligned along the span of an otherwise smooth wall to trip the boundary layer, the filtering region is shorter and growth of the vortices starts earlier. Also for the case of an isolated three-dimensional hump a rapid disturbance amplification is produced, but in this case the vortices remain confined and a very slow spanwise spreading of the perturbation occurs. In all naturally developing cases, where no specific wavelengths are explicity favored, the average spanwise wavelengths computed are very close to those of largest growth from the linear stability theory.     

An intermittency model for predicting roughness induced transition

An intermittency transport equation for RANS modeling, formulated in local variables, is extended for roughness-induced transition. To predict roughness effects in the fully turbulent boundary layer, published boundary conditions for k and ω are used. They depend on the equivalent sand-grain roughness height, and account for the effective displacement of wall distance origin. Similarly in our approach, wall distance in the transition model for smooth surfaces is modified by an effective origin, which depends on equivalent sand-grain roughness. Flat plate test cases are computed to show that the proposed model is able to predict transition onset in agreement with a data correlation of transition location versus roughness height, Reynolds number, and inlet turbulence intensity. Experimental data for turbine cascades are compared to the predicted results to validate the proposed model.     

Experimental study of boundary-layer transition on an airfoil induced by periodically passing wake

 Hot-wire measurements are performed in boundary-layer flows developing on a NACA 0012 airfoil over which wakes pass periodically. The periodic wakes are generated by rotating circular cylinders clockwise or counterclockwise around the airfoil. The time- and phase-averaged mean streamwise velocities and turbulence fluctuations are measured to investigate the phenomena of wake-induced transition. Especially, the phase-averaged wall shear stresses are evaluated using a computational Preston tube method. The passing wakes significantly change the pressure distribution on the airfoil, which has influence on the transition process of the boundary layer. The orientation of the passing wake alters the pressure distribution in a different manner. Due to the passing wake, the turbulent patches are generated inside the laminar boundary layer on the airfoil, and the boundary layer becomes temporarily transitional. The patches propagate downstream at a speed smaller than the free-stream velocity and merge together further downstream. Relatively high values of phase-averaged turbulence fluctuations in the outer part of the boundary layer indicate the possibility that breakdown occurs in the outer layer away from the wall. It is confirmed that the phase-averaged mean velocity profile has two dips in the outer region of the transitional boundary layer for each passing cycle. Received: 12 February 2001 / Accepted: 6 July 2001 Published online: 23 November 2001     

On Tollmien–Schlichting-like waves in streaky boundary layers

The linear stability of the boundary layer developing on a flat plate in the presence of finite-amplitude, steady and spanwise periodic streamwise streaks is investigated. The streak amplitudes considered here are below the threshold for onset of the inviscid inflectional instability of sinuous perturbations. It is found that, as the amplitude of the streaks is increased, the most unstable viscous waves evolve from two-dimensional Tollmien–Schlichting waves into three-dimensional varicose fundamental modes which compare well with early experimental findings. The analysis of the growth rates of these modes confirms the stabilising effect of the streaks on the viscous instability and that this stabilising effect increases with the streak amplitude. Varicose subharmonic modes are also found to be unstable but they have growth rates which typically are an order of magnitude lower than those of fundamental modes. The perturbation kinetic energy production associated with the spanwise shear of the streaky flow is found to play an essential role in the observed stabilisation. The possible relevance of the streak stabilising role for applications in boundary layer transition delay is discussed.     

Transitionally rough zero pressure gradient turbulent boundary layers

Near-wall measurements are performed to study the effects of surface roughness and viscous shear stresses on the transitionally rough regime (5 < k ⁺ < 70) of a zero pressure gradient turbulent boundary layer. The x-dependence is known from the eleven consecutive measurements in the streamwise direction, which allows for the computation of the streamwise gradients in the boundary layer equations. Thus, the skin friction is computed from the integrated boundary layer equation with errors of 3 and 5% for smooth and rough, respectively. It is found that roughness destroys the viscous layer near the wall, thus, reducing the contribution of the viscous stress in the wall region. As a result, the contribution in the wall shear stress due to form drag increases, while the viscous stress decreases. This yields Reynolds number invariance in the skin friction as k ⁺ increases into the fully rough regime. Furthermore, the roughness at the wall reduces the high peak of the streamwise component of the Reynolds stress in the near-wall region. However, for the Reynolds wall-normal and shear stress components, its contribution is not significantly altered for sand grain roughness.     

DNS and LES of turbulent flow in a closed channel featuring a pattern of hemispherical roughness elements

Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) were performed for fully-developed turbulent flow in channels with smooth walls and walls featuring hemispherical roughness elements at shear Reynolds numbers Re_τ = 180 and 400, with the goal of studying the effect of these roughness elements on the wall-layer structure and on the friction factor. The LES and DNS approaches were verified first by comparison with existing DNS databases for smooth walls. Then, a parametric study for the hemispherical roughness elements was conducted, including the effects of shear Reynolds number, normalized roughness height (k⁺ = 10–20) and relative roughness spacing (s⁺/k⁺ = 2–6). The sensitivity study also included the effect of distribution pattern (regular square lattice vs. random pattern) of the roughness elements on the walls. The hemispherical roughness elements generate turbulence, thus increasing the friction factor with respect to the smooth-wall case, and causing a downward shift in the mean velocity profiles. The simulations revealed that the friction factor decreases with increasing Reynolds number and roughness spacing, and increases strongly with increasing roughness height. The effect of random element distribution on friction factor and mean velocities is however weak. In all cases, there is a clear cut between the inner layer near the wall, which is affected by the presence of the roughness elements, and the outer layer, which remains relatively unaffected. The study reveals that the presence of roughness elements of this shape promotes locally the instantaneous flow motion in the lateral direction in the wall layer, causing a transfer of energy from the streamwise Reynolds stress to the lateral component. The study indicates also that the coherent structures developing in the wall layer are rather similar to the smooth case but are lifted up by almost a constant wall-unit shift y⁺ (∼10–15), which, interestingly, corresponds to the relative roughness k⁺ = 10.     

Vortex reynolds number in turbulent boundary layers

The effects of vortex Reynolds number on the statistics of turbulence in a turbulent boundary layer have been investigated. Vortex Reynolds number is defined as the ratio of circulation around the vortex structure to the fluid viscosity. The vortex structure of the outer region was modeled and a full numerical simulation was then conducted using a high-order spectral method. A unit domain of the outer region of a turbulent boundary layer was assumed to be composed of essentially three elements: a wall, a Blasius mean shear, and an elliptic vortex inclined at 45° to the flow direction. The laminar base-flow Reynolds number is roughly in the same range as that of a turbulent boundary layer based on eddy viscosity, and the vortex-core diameter based on the boundary-layer thickness is nearly the same as the maximum mixing length in a turbulent boundary layer. The computational box size, namely, 500, 150, and 250 wall units in the streamwise, surface-normal, and spanwise directions, respectively, is approximately the same as the measured quasi-periodic spacings of the near-wall turbulence-producing events in a turbulent boundary layer. The effects of vortex Reynolds number and the signs of the circulation on the moments of turbulence were examined. The signs mimic the ejection and sweep types of organized motions of a turbulent boundary layer. A vortex Reynolds number of 200 describes the turbulence moments in the outer layer reasonably well.