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.     

Properties of the streamwise velocity fluctuations in the inertial layer of turbulent boundary layers and their connection to self-similar mean dynamics

Well-resolved streamwise velocity measurements are used to investigate three measures of self-similarity in the spatial inertial sublayer of turbulent boundary layers. The emergence of self-similarity in the inertial sublayer requires a high Reynolds number, and thus a relatively wide range of ${\delta }^{+}=\delta u_{\tau }/\nu$ $(1400\lesssim {\delta }^{+}\lesssim 20\text{,}000)$ is explored. The measures investigated include the Kullback–Leibler divergence (KLD) used in turbulent flow analysis by Tsuji et al. (2005), the logarithmic decrease of the even statistical moments studied by Meneveau and Marusic (2013), and the diagnostic plot of Alfredsson and Örlü (2010). These measures are compared with the analyses of Fife et al. (2005) that determine and exploit an invariant form of the mean momentum equation. A primary focus is on domain(s) where the self-similar behaviors are analytically predicted and empirically observed. The present findings indicate that the approximately constant KLD and approximately logarithmic moment profiles reside in a region that is interior to the bounds of the self-similar inertial domain associated with the mean momentum equation. Conversely, the bounds of the self-similar region on the diagnostic plot correspond closely to the theoretically estimated bounds. Results are briefly discussed relative to Townsend’s notion of outer layer similarity, and, on the inertial domain, the physical existence of uniform momentum zones segregated by narrow vortical fissures.     

Prediction of resolvent mode shapes in supersonic turbulent boundary layers

This work applies resolvent analysis to compressible zero-pressure-gradient turbulent boundary layers with freestream Mach numbers between 2 and 4, focusing exclusively on large scale motions in the outer region of the boundary layer. We investigate the effects of Mach number on predicted flow structures, and in particular, look at how such effects may be attributed to changes in mean properties. By leveraging the similarity between the compressible and incompressible resolvent operators, we show that the shape of the streamwise velocity and temperature components of resolvent response modes in the compressible regime can be approximated by applying ideas from wavepacket pseudospectral theory to a simple scalar operator. This gives a means of predicting the shape of resolvent mode components for compressible flows without requiring the singular value decompositions of discretized operators. At a Mach number of 2, we find that accurate results are obtained from this approximation when using the compressible mean velocity profile. At Mach numbers of 3 and 4, the quantitative accuracy of these predictions is improved by also considering a local effective Reynolds number based on the local mean density and viscosity.     

Hairpin vortices in the largest scale of turbulent boundary layers

We have conducted direct numerical simulations of a turbulent boundary layer for the momentum-thickness-based Reynolds number ${\mathit{Re}}_{\theta }$ = 180–4600. To extract the largest-scale vortices, we coarse-grain the fluctuating velocity fields by using a Gaussian filter with the filter width comparable to the boundary layer thickness. Most of the largest-scale vortices identified by isosurfaces of the second invariant of the coarse-grained velocity gradient tensor are similar to coherent vortices observed in low-Reynolds-number regions, that is, hairpin vortices or quasi-streamwise vortices inclined to the wall. We also develop a percolation analysis to investigate the threshold-dependence of the isosurfaces and objectively identify the largest-scale hairpin vortices in terms of the coarse-grained vorticity, which leads to the quantitative evidence that they never disappear even in fully developed turbulent regions. Hence, we conclude that hairpin vortices exist in the largest-scale structures irrespective of the Reynolds number.     

Turbulent boundary layers over sparsely-spaced rod-roughened walls

Direct numerical simulations (DNSs) of spatially developing turbulent boundary layers (TBLs) over sparsely-spaced two-dimensional (2D) rod-roughened walls were performed. The rod elements were periodically arranged along the streamwise direction with pitches of p_x/k = 8, 16, 32, 64 and 128, where p_x is the streamwise spacing of the rods, and k is the roughness height. The Reynolds number based on the momentum thickness was varied from Re_θ = 300–1400, and the height of the roughness element was k = 1.5θ_in, where θ_in is the momentum thickness at the inlet. The characteristics of the TBLs, such as the friction velocity, mean velocity, and Reynolds stresses over the rod-roughened walls, were examined by varying the spacing of the roughness features (8 ⩽ p_x/k ⩽ 128). The outer-layer similarity between the rough and smooth walls was established for the sparsely-distributed rough walls (p_x/k ⩾ 32) based on the profiles of the Reynolds stresses, whereas those are not for p_x/k = 8 and 16. Inspection of the interaction between outer-layer large-scale motions and near-wall small-scale motions using two-point amplitude modulation (AM) covariance showed that modulation effect of large-scale motions on near-wall small-scale motions was strongly disturbed over the rough wall for p_x/k = 8 and 16. For p_x/k ⩾ 32, the flow that passed through the upstream roughness element transitioned to a smooth wall flow between the consecutive rods. The strong influence of the surface roughness in the outer layer for p_x/k = 8 and 16 was attributed to large-scale erupting motions by the surface roughness, creating both upward shift of the near-wall turbulent energy and active energy production in the outer layer with little influence on the near-wall region.     

A silicon transducer for the determination of wall-pressure fluctuations in turbulent boundary layers

Small and sensitive silicon sensors for turbulent wall-pressure fluctuation measurements have been designed and fabricated using microelectronic technology. For the detection of the pressure fluctuations piezoresistive gauges are placed on a diaphragm and the resistance of these gauges depends on the stresses in the diaphragm. For the determination of the performance of these pressure transducers comparisons with conventional microphones were carried out in a well-defined two-dimensional boundary layer. Power spectra from the silicon pressure transducer revealed a slope in the inertial sublayer corresponding approximately to the 5/3-law of Kolmogorov, and the normalized wall pressure fluctuations agreed well with other direct measurements.     

Effects of the roughness height in turbulent boundary layers over rod- and cuboid-roughened walls

Direct numerical simulations (DNSs) of spatially developing turbulent boundary layers (TBLs) over two-dimensional (2D) rod-roughened walls and three-dimensional (3D) cuboid-roughened walls are conducted to investigate the effects of the roughness height on the flow characteristics in the outer layer. The rod elements are periodically aligned along the downstream direction with a pitch of p_x/θ_in = 12, and the cuboid elements are periodically staggered with a pitch of p_x/θ_in = 12 and p_z/θ_in = 3, where p_x and p_z are correspondingly the streamwise and spanwise pitches of the roughness and θ_in is the momentum thickness at the inlet. The first surface roughness is placed 80θ_in downstream from the inlet, leading to a step change from a smooth to rough surface. The rod and cuboid roughness height (k) is varied in the range of 0.1 ≤ k/θ_in ≤ 1.8 (13 ≤ δ/k ≤ 285), respectively (δ is the boundary layer thickness), and the Reynolds number based on the momentum thickness (θ) is varied in the range of Re_θ = 300 ~ 1400. For each case, the self-preservation form of the velocity-defect and the turbulent Reynolds stresses is achieved along the downstream direction. As the roughness height increases, the roughness function (ΔU⁺) extracted from the mean velocity profiles increases, although the velocity-defect profiles for the rough-wall cases show good agreement with the profile from the smooth-wall case. The magnitude of the Reynolds stresses in the outer layer increases with an increase of k/δ. The outer layer similarity between the flows over the rough and smooth-walls is found when δ/k ≥ 250 and 100 for the 2D rod and 3D cuboid, respectively. The continuous increase of the Reynolds stresses in the outer layer with an increase of k/δ is explained by a large population of very long structures over the rough-wall flows. Because the characteristic width of the structures increases continuously with an increase of k/δ for the rod and cuboid roughness, a wide width of the structures leads to frequent spanwise merging between adjacent structures. The active spanwise merging events with an increase of k/δ increase the streamwise coherence of the structures with the appearance of significant meandering.     

The turbulent/non-turbulent interface in an adverse pressure gradient turbulent boundary layer

The turbulent/non-turbulent interface (TNTI) in an adverse pressure gradient (APG, β = 1.45) turbulent boundary layer (TBL) is explored here by using direct numerical simulation (DNS) data; β is the Clauser pressure gradient parameter. For comparison, the DNS data for a zero pressure gradient (ZPG) TBL is included. The interface is extracted with an approach based on enstrophy criteria. Depending on the enstrophy, the outer boundary layer flow can be classified into the free stream, boundary layer wake, and intermittent flow regimes. The fractal dimension of the interface is obtained by using the box-counting algorithm, and was found to be constant over a long range of box sizes. The TNTI shows a monofractal behavior. The geometric complexity of a TNTI can be determined in terms of the genus, which is defined as the number of handles in a geometric object. We examine the volume and projection area of the genus of the TNTI to analyze the entrainment process. The geometric complexity of the APG TBL interface and the local entrainment are greater than those of the ZPG TBL, as is evident in the increases in the genus near the interface. The local entrainment velocity is dominantly affected by the viscous diffusion at the interface.     

Direct numerical simulation of stable and unstable turbulent thermal boundary layers

This paper presents direct numerical simulations (DNS) of stable and unstable turbulent thermal boundary layers. Since a buoyancy-affected boundary layer is often encountered in an urban environmental space where stable and unstable stratifications exist, exploring a buoyancy-affected boundary layer is very important to know the transport phenomena of the flow in an urban space. Although actual observation may qualitatively provide the characteristics of these flows, the relevant quantitative turbulent quantities are very difficult to measure. Thus, in order to quantitatively investigate a buoyancy-affected boundary layer in detail, we have here carried out for the first time time- and space-developing DNS of slightly stable and unstable turbulent thermal boundary layers. The DNS results show the quantitative turbulent statistics and structures of stable and unstable thermal boundary layers, in which the characteristic transport phenomena of thermally stratified boundary layers are demonstrated by indicating the budgets of turbulent shear stress and turbulent heat flux. Even though the input of buoyant force is not large, the influence of buoyancy is clearly revealed in both stable and unstable turbulent boundary layers. In particular, it is found that both stable and unstable thermal stratifications caused by the weak buoyant force remarkably alter the structure of near-wall turbulence.     

Bursting frequency in turbulent boundary layers

Employing laser Doppler anemometry and VITA techniques, the bursting frequency in turbulent boundary layers has been measured over the Reynolds-number range 320 to 1470. The result indicates that the mean and non-dimensional bursting frequency scaled with the variables appropriate for the wall region was constant and independent of Reynoids number. When the same data are plotted using the outer variables of boundary layer to normalize the bursting frequency, the non-dimensional frequency increases as the Reynolds number increases. This is in agreement with the results of Blackwelder et al. (1983) who used hot wire anemometry and VITA technique. The project is supported by the National Natural Science Foundation of China     

Development of customized shear layers on smooth and rough surfaces

Shear layers having different structural properties were produced downstream of devices used to artificially thicken turbulent boundary layers. The means used to produce different structural characteristics are described, along with the effects of changes in structure on wall heat transfer. Results from layers developing over smooth and rough surfaces indicate that alterations of artificial thickening device geometry resulted in larger variations in wall heat transfer near smooth surfaces. The most significant of these occured when alterations were made of inner boundary layer regions, where mean velocity shows a logarithmic dependence on distance from the wall. Outer region changes in mean velocity and turbulence profiles resulted in less significant changes in wall heat transfer, particularly in the flows over the rough wall.     

Heat transfer in film-cooled turbulent boundary layers at different blowing ratios as affected by longitudinal vortices

Effects of embedded longitudinal vortices on heat transfer in film-cooled turbulent boundary layers at different blowing ratios are discussed. These results were obtained in boundary layers at free-stream velocities of 10 and 15 m/s. Film coolant was injected from a single row of holes at blowing ratios of 0.47–1.26. A single longitudinal vortex was induced upstream of the film-cooling holes using a half-delta wing attached to the wind tunnel floor. Heat transfer measurements were made downstream of injection using a constant heat flux surface with 126 thermocouples for surface temperature measurements. For all blowing ratios examined, the embedded vortices cause significant alterations to wall heat transfer and to film cooling distributions. Measurrments of mean temperatures and mean velocities in spanwise planes show that high wall heat transfer regions are associated with regions of high near-wall longitudinal velocity where very little film coolant is present. In addition to high heat transfer regions associated with the vortex downwash, there are also secondary heat transfer peaks. These secondary peaks develop due to shear layer mixing and interaction between the vortex and cooling jets and become higher in magnitude and more persistent with downstream distance as the blowing ratio increases from 0.47 to 1.26.     

Numerical prediction of turbulent thermocapillary convection in superposed fluid layers with a free interface

The present paper introduces a new numerical method for predicting the characteristics of thermocapillary turbulent convection in a differentially-heated rectangular cavity with two superposed and immiscible fluid layers. The unsteady Reynolds form of the Navier–Stokes equations and energy equation are solved by using the control volume approach on a staggered grid system using SIMPLE algorithm. The turbulence quantities are predicted by applying the standard k–ε turbulence model. The level set formulation is applied for predicting the topological changes of the interface separating the two fluid layers and to provide an accurate and robust modeling of the interfacial normal and tangential stresses. The computational results obtained showed good agreement when compared with the previous experimental, numerical and analytical benchmark data for different validation cases in both laminar and turbulent regimes. The present numerical method is then applied to predict the velocity and temperature distribution in two immiscible liquid layers with undeformable interface for a wide range of Marangoni numbers. The laminar-turbulent transition is demonstrated by obtaining the turbulence features at high interfacial temperature gradient which is characterized by high Marangoni number. The effect of increasing Marangoni number on the interface dynamics in turbulent regime is also investigated.     

Measurements of wall shear stress of a turbulent boundary layer using a micro-shear-stress imaging chip

Measurements of wall shear-stress streaks of a turbulent boundary layer in the channel flow were carried out using a MEMS-based micro-shear-stress imaging chip, which contains about 100 sensors. The chip is designed and fabricated by surface micromachining technology. One arrray of 25 micro-shear-stress sensors in the chip that covers a length of 7.5 mm is used to measure the instantaneous spanwise distribution of the surface shear stress. The statistics of high shear-stress streaks were established. Based on the measurement, the physical quantities associated with the high shear-stress streaks, such as their length, width and peak shear-stress level, were obtained. We found out that a high correlation exists between the peak shear-stress level and front-end shear-stress slope of a high shear-stress streak. This important property is currently being applied to the deisgn of a real-time flow control logic.     

On near-wall turbulence-generating events in a turbulent boundary layer on a riblet surface

The boundary layer over a drag reducing riblet surface is investigated using hot-wire anemometry and flow visualisation. The concept of a riblet sublayer is introduced, and a definition is proposed in terms of a region of reduced turbulence energy production formed near the wall by the addition of riblets. The hot wire records are examined using a modified form of quadrant analysis, and results obtained over plain and riblet surfaces are compared. Close to the wall, the addition of riblets produces a marked reduction in the occurrence of ejection (2nd quadrant) events. A corresponding increase in the incidence of sweep (4th quadrant) events is accompanied by the development of a strong tendency toward a preferred event duration, and a preferred interval between events. These changes diminish rapidly with distance from the surface, becoming almost undetectable beyondy ⁺=40. They are discussed in the light of flow visualisation results, and interpreted in terms of mechanisms associated with the interaction between the riblets and the inner boundary layer flow structures. A conceptual model of the flow mechanisms in the riblet sublayer is proposed.