Smooth and rough turbulent boundary layers at high Reynolds number

A 2-D turbulent boundary layer experiment with zero pressure gradient (ZPG) has been carried out over a rough and a smooth surface using two cross hot-wire probes. Wind tunnel speeds of 10 m/s and 20 m/s were set up in order to investigate the effects of the upstream conditions and the Reynolds number on the downstream flow. For a given set of upstream conditions, such as the wind tunnel speed, trip wire size and location, the three components of the velocity field were measured from about 14 m from the inlet of the wind tunnel to 30 m downstream. This experiment is unique because it achieves Reynolds numbers as high as R120,000, for which measurements of the mean velocity are reported. It is shown that by fixing the upstream conditions, the mean deficit profiles collapse with the freestream velocity, , but to different curves depending on the upstream conditions and surface roughness. Moreover, the effects of the upstream conditions, the Reynolds number, and roughness are completely removed from the outer flow when the mean deficit profiles are normalized by the Zagarola/Smits scaling, . Consequently, the true asymptotic profile in the turbulent boundary layer is found in ZPG flow regardless of the range of Reynolds number, surface conditions and initial conditions.     

Upstream roughness and Reynolds number effects on turbulent flow structure over forward facing step

This paper presents results of experiments conducted to investigate the effects of Reynolds number and upstream wall roughness on the turbulence structure in the recirculation and recovery regions of a smooth forward facing step. A reference smooth upstream wall and a rough upstream wall made from sand gains were studied. For the smooth upstream wall, experiments were conducted at Reynolds number based on the freestream velocity and step height (h), Re_h = 4940, 8400 and 8650. The rough wall experiments was performed at Re_h = 5100, 8200 and 8600 to closely match the corresponding Re_h experiment over the smooth wall. The reattachment lengths in the smooth wall experiments were L_r/h ≈ 2.2, but upstream roughness significantly reduced these values to L_r/h ≈ 1.3. The integral scales within the recirculation bubbles were independent of upstream roughness and Reynolds number; however, upstream roughness significantly increased the spatial coherence and integral scales outside the recirculation bubbles and in the recovery region. Irrespective of the upstream wall condition, the redeveloping boundary layer recovered at 25h from reattachment.     

Computation of high Reynolds number flow around a circular cylinder with surface roughness

Incompressible high-Reynolds-number flows around a circular cylinder are analyzed by direct integration of the Navier-Stokes equations using finite-difference method. A generalized coordinate system is used so that a sufficient number of grid points are distributed in the boundary layer and the wake. A numerical scheme which suppresses non-linear instability for calculations of high-Reynolds-number flows is developed. The computation of an impulsively started flow at Re = 1200 is compared with corresponding experimental observations, and excellent agreements are obtained.A series of computations are carried out on the flow around a circular cylinder with surface roughness. The height of the roughness in these computations is 0.5% of the diameter. The range of Reynolds numbers is from 10³ to 10⁵; no turbulence model is employed. Sharp reduction of drag coefficient is observed near Re = 2 × 10⁴, which indicates that the critical Reynolds number is captured in the present computation.     

Some Reynolds number effects on two- and three-dimensional turbulent boundary layers

Experimental data for a two-dimensional (2-D) turbulent boundary layer (TBL) flow and a three-dimensional (3-D) pressure-driven TBL flow outside of a wing/body junction were obtained for an approach Reynolds number based on momentum thickness of Re _θ =23,200. The wing shape had a 3:2 elliptical nose, NACA 0020 profiled tail, and was mounted on a flat wall. Some Reynolds number effects are examined using fine spatial resolution (Δy ⁺=1.8) three-velocity-component laser-Doppler velocimeter measurements of mean velocities and Reynolds stresses at nine stations for Re _θ =23,200 and previously reported data for a much thinner boundary layer at Re _θ =5,940 for the same wing shape. In the 3-D boundary layers, while the stress profiles vary considerably along the flow due to deceleration, acceleration, and skewing, profiles of the parameter correlate well and over available Reynolds numbers. The measured static pressure variations on the flat wall are similar for the two Reynolds numbers, so the vorticity flux and the measured mean velocities scaled on wall variables agree closely near the wall. The stresses vary similarly for both cases, but with higher values in the outer region of the higher Re _θ case. The outer layer turbulence in the thicker high Reynolds number case behaves similarly to a rapid distortion of the flow, since stream-wise vortical effects from the wall have not diffused completely through the boundary layer at all measurement stations. Received: 9 June 2000/Accepted: 26 January 2001     

Reynolds number dependence of Reynolds and dispersive stresses in turbulent channel flow past irregular near-Gaussian roughness

Direct numerical simulations of fully-developed turbulent channel flow with irregular rough walls have been performed at four friction Reynolds numbers, namely, 180, 240, 360 and 540, yielding data in both the transitionally- and fully-rough regime. The same roughness topography, which was synthesised with an irregular, isotropic and near-Gaussian height distribution, is used in each simulation. Particular attention is directed towards the wall-normal variation of flow statistics in the near-roughness region and the fluid-occupied region beneath the crests, i.e. within the roughness canopy itself. The goal of this study is twofold. (i) Provide a detailed account of first- and second-order double-averaged velocity statistics (including profiles of mean velocity, dispersive stresses, Reynolds stresses, shear stress gradients and an analysis of the mean force balance) with the overall aim of understanding the relative importance of “form-induced” and “turbulence-induced” quantities as a function of the friction Reynolds number. (ii) Investigate the possibility of predicting the levels of streamwise dispersive stress using a phenomenological closure model. Such an approach has been applied successfully in the context of idealised vegetation canopies (Moltchanov & Shavit, 2013, Water Resour. Res., vol. 49, pp. 8222-8233) and is extended here, for the first time, to an irregular rough surface. Overall, the results reveal that strong levels of dispersive stress occur beneath the roughness crests and, for the highest friction Reynolds number considered in this study, show that the magnitude (and gradient) of these “form-induced” stresses exceed their Reynolds stress counterparts. In addition, this study emphasises that the dominant source of spatial heterogeneity within the irregular roughness canopy are “wake-occupied” regions and that a suitable parameterisation of the wake-occupied area is required to obtain an accurate prediction of streamwise dispersive stress.     

Full 3D correlation tensor computed from double field stereoscopic PIV in a high Reynolds number turbulent boundary layer

The turbulence structure near a wall is a very active subject of research and a key to the understanding and modeling of this flow. Many researchers have worked on this subject since the fifties Hama et al. (J Appl Phys 28:388–394, 1957). One way to study this organization consists of computing the spatial two-point correlations. Stanislas et al. (C R Acad Sci Paris 327(2b):55–61, 1999) and Kahler (Exp Fluids 36:114–130, 2004) showed that double spatial correlations can be computed from stereoscopic particle image velocimetry (SPIV) fields and can lead to a better understanding of the turbulent flow organization. The limitation is that the correlation is only computed in the PIV plane. The idea of the present paper is to propose a new method based on a specific stereoscopic PIV experiment that allows the computation of the full 3D spatial correlation tensor. The results obtained are validated by comparison with 2D computation from SPIV. They are in very good agreement with the results of Ganapthisubramani et al. (J Fluid Mech 524:57–80, 2005a).     

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.     

Slip boundary effect on the critical Reynolds number of subcritical transition in channel flow

In this letter, the effect of slip boundary on the origin of subcritical transition in two-dimensional channel flows is studied numerically and theoretically. It is shown that both the positive and the negative slip lengths will increase the critical Reynolds number of localized wave packet and hence postpone the transition. By applying a variable transformation and expanding the variables about a small slip length, it is illustrated that the slip boundary effect only exists in the second and higher order modulations of the no-slip solution, and hence explains the power law found in simulations, i.e. the relative increment of the critical Reynolds number due to the slip boundary is proportional to the square of the slip length.     

Effect of leading-edge tubercles on the flow over low-aspect-ratio wings at low Reynolds number

Two-dimensional time-resolved particle image velocimetry (TR-PIV) and stereographic particle image velocimetry (SPIV) techniques were used to investigate the effect of leading-edge tubercles on the flow over low-aspect-ratio wing models. The angle of attack is fixed at 10°, and the Reynolds number based on chord length is 5.8 × 10³. It is shown that the leading-edge tubercles can effectively mitigate flow separation in the model and also reduce the contribution of wake vortex to the fluctuating energy of flow. Counter-rotating vortex pairs (CVPs) initiated from the peak of leading-edge tubercles can promote nearby momentum exchange, enhance mixing of the flow and increase the energy contained in the boundary layer, which results in resisting the larger adverse pressure gradient. Therefore, it is concluded that CVPs play an important role in mitigating the flow separation for wings with leading-edge tubercles.     

3-D measurements of separated flow over rigid plates with spanwise periodic cambering at low Reynolds number

Measurements of the flow field around a flat plate and rigid plates with spanwise periodic cambering were performed using volumetric three-component velocimetry (V3V) at a Reynolds numbers of 28,000 at α=12° where the flow is fully separated. The Reynolds normal and shear stresses, and the streamwise, spanwise and normal components of the vorticity vector are investigated for three-dimensionality. Flow features are discussed in context of the periodic cambering and corresponding aerodynamic force measurements. The periodic cambering results in spanwise variation in the reversed-flow region, Reynolds stresses and spanwise vorticity. These spanwise variations are induced by streamwise and normal vortices of opposite directions of rotation. Moreover, measurements were carried out for the cambered plates at α=8°, where a long separation bubble exists, to further understand the behavior of the streamwise and normal vortices. These vortices become more organized and increase in strength and size at the lower angle of attack. It is also speculated that these vortices contribute to the increase in lift at and beyond the onset of stall angle of attack.     

A short note on the general boundary element method for viscous flows with high Reynolds number

In this paper, the general boundary element method and the parallel computation are employed to solve laminar viscous flows in a driven square cavity, governed by the exact Navier–Stokes equations. Using the solution at Re=0 as the initial approximation, the convergent numerical results for high Reynolds number at Re=7500 are obtained, for the first time, by the boundary element method. This verifies the validity and great potential of the general boundary element method for highly non‐linear problems, which may greatly enlarge application regions of the boundary element method in science and engineering. Copyright © 2003 John Wiley & Sons, Ltd.     

Quantifying the nonlinear mode competition in the flow over an open cavity at medium Reynolds number

Our purpose is to quantify the rate of intermittency of nonlinearly competing modes, in a dominantly mode-switching scenario. What is the rate of presence of each mode? Can they simultaneously appear in, or disappear from the signal? The study is done in the context of open flows, exhibiting self-sustained oscillations, where air is here flowing over an open cavity. Reynolds numbers are of the order of 14,000. Velocity measurements downstream of the cavity are based on a laser Doppler velocimetry technique. We propose two methods to estimate the rate of presence of each mode: one based on a complex demodulation technique, the other on the distribution of the state vectors in the phase portrait of the signal.     

Coherent and turbulent process analysis in the flow past a circular cylinder at high Reynolds number

This article presents an experimental study of the turbulent flow past a circular cylinder at high Reynolds number by means of advanced optical measurements techniques. Following previous studies using standard PIV and stereoscopic PIV (3C PIV), TRPIV and 3C-TRPIV have been employed in low subsonic wind tunnel environment. The database consisting of statistical and time-dependent fields aims at providing a physical analysis of the coherent and turbulent part, as well as a proper basis for validation and improvement of recent turbulence modelling approaches for strongly detached flow at high Reynolds number. As the nonlinear interaction between the coherent and turbulent dynamics have to be taken into account in a model, particular attention is paid to a decomposition of the flow into a coherent and a turbulent part, and to the analysis of their dynamics. This is achieved both using phase averaging and Proper Orthogonal Decomposition. For phase averaging, the two first POD coefficients are used for the evaluation of the vortex shedding phase angle. Furthermore, selected results of a Detached Eddy Simulation which had been validated by means of the experiment, are also presented to contribute to the physical analysis. The present study's experimental data resolved in space and time allow the confirmation of the conditional averaging for the turbulent stresses evaluation, by alleviating their overestimation due to phase jitter that occurs between the trigger signal and the velocity, when the phase angle is determined from a wall pressure signal. A more accurate physical analysis of the flow is achieved, particularly regarding the occurrence of irregular vortex shedding.     

The restoration coefficient and Reynolds analogy in a boundary layer with injection and suction over the entire Prandtl number range

The restoration and Reynolds analogy coefficients are calculated for a laminar self-similar boundary layer on a permeable plate over the entire possible range of variation in the Prandtl number and the injection and suction parameter.     

Development of localized arc filament RF plasma actuators for high-speed and high Reynolds number flow control

Recently developed localized arc filament plasma actuators (LAFPAs) have shown tremendous control authority in high-speed and high Reynolds number flow for mixing enhancement and noise mitigation. Previously, these actuators were powered by a high-voltage pulsed DC plasma generator with low energy coupling efficiency of 5–10%. In the present work, a new custom-designed 8-channel pulsed radio frequency (RF) plasma generator has been developed to power up to 8 plasma actuators operated over a wide range of forcing frequencies (up to 50 kHz) and duty cycles (1–50%), and at high energy coupling efficiency (up to 80–85%). This reduces input electrical power requirements by approximately an order of magnitude, down to 12 W per actuator operating at 10% duty cycle. The new pulsed RF plasma generator is scalable to a system with a large number of channels. Performance of pulsed RF plasma actuators used for flow control was studied in a Mach 0.9 circular jet with a Reynolds number of about 623,000 and compared with that of pulsed DC actuators. Eight actuators were distributed uniformly on the perimeter of a 2.54-cm diameter circular nozzle extension. Both types of actuators coupled approximately the same amount of power to the flow, but with drastically different electrical inputs to the power supplies. Particle image velocimetry measurements showed that jet centerline Mach number decay produced by DC and RF actuators operating at the same forcing frequencies and duty cycles is very similar. At a forcing Strouhal number near 0.3, close to the jet column instability frequency, well-organized periodic structures, with similar patterns and dimensions, were generated in the jets forced by both DC and RF actuators. Far-field acoustic measurements demonstrated similar trends in the overall sound pressure level (OASPL) change produced by both types of actuators, resulting in OASPL reduction up to 1.2–1.5 dB in both cases. We conclude that pulsed RF actuators demonstrate flow control authority similar to pulsed DC actuators, with a significantly reduced power budget.     

Effects of aspect ratio and orientation on the wake characteristics of low Reynolds number flow over a triangular prism

The wake characteristics of unconfined flows over triangular prisms of different aspect ratios have been numerically analysed in the present work. For this purpose, a fixed Cartesian-grid based numerical technique that involves the porous medium approach to mimic the effect of solid blockage has been utilised. Correspondingly, laminar flow simulations ranging from the sub-critical regime (before the onset of vortex shedding) to the super-critical regime have been considered here within the limits of two-dimensionality. In the sub-critical regime, correlations relating the wake bubble length with Reynolds number (Re) have been proposed for various aspect ratios. Also, the effects of aspect ratio and Reynolds Number on the drag force coefficient (C_D) have been characterised for two different geometrical orientations of the prism (base or apex facing the flow). Subsequently, the critical Reynolds number at the onset of vortex shedding has been predicted for each of the aspect ratio considered, by an extrapolation procedure. The unsteady flow characteristics of the super-critical regime are finally highlighted for different aspect ratios and triangular orientations considered in the study.     

Enhancement of finite Reynolds number effects: Inner-outer sublayer interaction in the turbulent boundary layer

Residual Reynolds number effects in the established data for the velocity profile in turbulent boundary layers (and in pipe or channel flows) are found to be remarkably large. We combine two eddy-viscosity models (with overlapping validity in the inertial sublayer) and show (both analytically and numerically) that this enhancement (which involves a viscous correlation length) arises from inner-outer sublayer interaction.     

Simultaneous entrainment of oil and water droplets in high Reynolds number gas turbulence in horizontal pipe flow

We develop a 1D cross sectional concentration profile model for oil and water droplets that coexist in the turbulent gas phase (of Re ∼ 10⁶) in near horizontal stratified pipe flows. Entrainment of the oil and water mixture from a liquid film near the bottom of the pipe into the gas is modeled based on earlier single-fluid entrainment correlations. A Gamma distribution for the droplet sizes based on the breakup of liquid filaments, is adopted. An explicit algebraic–exponential formula for the total concentration profile for either phase can then be derived.     

Shock-vortex shear-layer interaction in the transonic flow around a supercritical airfoil at high Reynolds number in buffet conditions

This paper provides a conceptual analysis and a computational model for how the unsteady ‘buffeting’ phenomenon develops in transonic, low incidence flow around a supercritical aerofoil, the OAT15A, at Reynolds number of 3.3 million. It is shown how a low-frequency buffet mode is amplified in the shock-wave region and then develops upstream and downstream interaction with the alternating von Kármán eddies in the wake past the trailing-edge as well as with the shear-layer, Kelvin–Helmholtz vortices. These interactions are tracked by wavelet analysis, autoregressive (AR) modelling and by Proper Orthogonal Decomposition. The frequency modulation of the trailing-edge instability modes is shown in the spectra and in the wall-pressure fluctuations. The amplitude modulation of the buffet and von Kármán modes has been also quantified by POD analysis. The thinning of the shear layers, both at the outer edge of the turbulent boundary layers and the wake, caused by an ‘eddy-blocking’ mechanism is modelled by stochastic forcing of the turbulent kinetic energy and dissipation, by small-scale straining of the higher-order POD modes. The benefits from thinning the shear-layers by taking into account the interfacial dynamics are clearly shown in the velocity profiles, and wall pressure distribution in comparison with the experimental data.