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.     

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.     

Experimental study of parallel and inclined turbulent wall jets

An experimental study of the flow field in a two-dimensional wall jet has been conducted. All measurements were carried out using hot-wire anemometry. The experimental facility has a rectangular slot nozzle of high aspect ratio l/b = 100 (where l and b are the length and height slot, respectively). Mean velocities and Reynolds stresses were determined with three nozzle Reynolds numbers (Re = 1 × 10⁴, 2 × 10⁴ and 3 × 10⁴) and four different inclination angles between the wall and the flow velocity at the nozzle (β = 0°, 10°, 20° and 30°). Results indicate that all wall jets are self-preserving in the developed region. Normal to the wall two regions can be identified: one similar to a plane free jet and the other similar to a boundary layer. Downstream the interaction between these two regions creates a mixed or third region. The logarithmic region increases with the distance from the nozzle and with the Reynolds number. For the inclined wall jet, the spreading rate expressed in terms of jet half-width or maximum velocity decay with respect to the streamwise distance, asymptotes to a linear law. The streamwise locations where the jet becomes self-similar are farther from the exit than in parallel wall jet. The slope of both half-width and maximum velocity decay in the developed region are affected by both wall jet inclination angle and nozzle exit Reynolds number.     

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.     

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.     

Effects of irregular two-dimensional and three-dimensional surface roughness in turbulent channel flows

Wall-resolved Large Eddy Simulation of fully developed turbulent channel flows over two different rough surfaces is performed to investigate on the effects of irregular 2D and 3D roughness on the turbulence. The two geometries are obtained through the superimposition of sinusoidal functions having random amplitudes and different wave lengths. In the 2D configuration the irregular shape in the longitudinal direction is replicated in the transverse one, while in the 3D case the sinusoidal functions are generated both in streamwise and spanwise directions. Both channel walls are roughened in such a way as to obtain surfaces with statistically equivalent roughness height, but different shapes. In order to compare the turbulence properties over the two rough walls and to analyse the differences with a smooth wall, the simulations are performed at the same Reynolds number Re_τ = 395. The same mean roughness height h = 0.05δ (δ the half channel height) is used for the rough walls.The roughness function obtained with the 3D roughness is larger than in the 2D case, although the two walls share the same mean height. Thus, the considered irregular 3D roughness is more effective in reducing the flow velocity with respect to the 2D roughness, coherently with the literature results that identified a clear dependence of the roughness function on the effective slope (see Napoli et al. (2008)), higher in the generated 3D rough wall. The analysis of higher-order statistics shows that the effects of the roughness, independently on its two- or three-dimensional shape, are mainly confined in the inner region, supporting the Townsend’s wall similarity hypothesis. The tendency towards the isotropization is investigated through the ratio between the resolved Reynolds stress components, putting in light that the 3D irregular rough wall induces an higher reduction of the anisotropy, with respect to the 2D case.     

Flow of drag-reducing surfactant solutions in rough pipes

The paper reports the results of experimental study of the flow of hexadecyltrimethylammonium chloride (CTAC) solutions with addition of sodium salicylate (NaSal) in the rough pipes. Measurements were performed in the range of the surfactant concentration from 200 to 400 ppm at a constant molar ratio CTAC/NaSal of 1:2. Five pipes of the relative roughness k/D varying from 1.2 × 10^?2 to 5.6 × 10^?2, obtained by the covering of inner surface of the pipes with glued silicon carbide particles of different size, were studied. The roughness was observed to increase the drag of flow of CTAC/NaSal solutions already at Reynolds numbers higher than 800. With increasing relative roughness k/D, the critical value of Reynolds number, at which the drag reduction disappears, was found to decrease. However, no influence of the roughness on the critical shear stress was noted. The ratio of the critical Reynolds number for rough pipes to that of hydraulically smooth pipes was independent of the surfactant concentration. The degree of drag reduction by the flow of surfactants was greater in rough pipes than in smooth pipes.     

Influence of adverse pressure gradient on rough-wall turbulent flows

This paper presents an experimental investigation of adverse pressure gradient turbulent flow over two rough surfaces and a reference smooth surface. The adverse pressure gradient was produced in an asymmetric diffuser whose opening angle was 3°. The rough surfaces comprised sand grains and gravels of nominal mean diameters of 1.55 mm and 4.22 mm, respectively. The tests were conducted at an approach flow velocity of 0.5 m/s and the momentum thickness Reynolds number varied from 900 to 3000. A particle image velocimetry technique was used for the velocity measurements. Profiles of the mean velocity, turbulent intensities, Reynolds stress ratios, mixing length, eddy viscosity and the production terms were then obtained to document the effects of adverse pressure gradient (APG) on low Reynolds number rough-wall turbulent boundary layers. The results indicate that APG thickens the boundary layer and roughness sublayer. The APG and surface roughness also enhanced the production of turbulence as well as the turbulence level when compared with the smooth-wall data.     

Turbulent natural convection scaling in a vertical channel

Using direct numerical simulation (DNS) data, this study appraises existing scaling laws in literature for turbulent natural convection of air in a differentially heated vertical channel. The present data is validated using past DNS studies, and covers a range of Rayleigh number, Ra between 5.4 × 10⁵ and 2.0 × 10⁷. We then appraise and compare the various scaling laws proposed by Versteegh and Nieuwstadt, 1999, Hölling and Herwig, 2005, Shiri and George, 2008, George and Capp, 1979 with the profiles of the mean temperature defect, mean streamwise velocity, normal velocity fluctuations, temperature fluctuations and Reynolds shear stress. Based on the arguments of an inner (near-wall) and outer (channel centre) region, the data is found to support a minus one-third power law for the mean temperature in an overlap region. Using the inner and outer temperature profiles, an implicit heat transfer equation is obtained and we show that a correction term is non-negligible for the present Ra range when compared with explicit equations found in literature. In addition, we determined that the mean streamwise velocity and normal velocity fluctuations collapse in the inner region when using the outer velocity scale. We also find that the temperature fluctuations scale in inner coordinates, in contrast to the outer scaling behaviour reported in the past. Lastly, we show evidence of an incipient proportional relationship between friction velocity, u_τ, and the outer velocity scale, u_o, with increasing Ra.     

Heat transfer due to impinging co-axial jets and the jets’ fluid flow characteristics

This investigation had multiple goals. One goal was to obtain definitive information about the heat transfer characteristics of co-axial impinging jets, and this was achieved by measurements of the stagnation-point, surface-distribution and average heat transfer coefficients. These results are parameterized by the Reynolds number Re which ranged from 5000 to 25,000, the dimensionless separation distance between the jet exit and the impingement plate H/D (4–12), and the ratio of the inner diameters of the inner and outer pipes d/D (0–0.55). The d/D = 0 case corresponds to a single circular jet. The other major goal of this work was to quantify the velocity field of co-axial free jets (impingement plate removed). The velocity-field study included both measurements of the mean velocity and the turbulence intensity.It was found that the variation of the stagnation-point heat transfer coefficient with d/D attained a maximum at d/D = 0.55. Furthermore, the variation of the local heat transfer coefficient across the impingement surface was more peaked for d/D = 0 and became flatter with decreasing d/D. This suggests that for cooling a broad expanse of surface, co-axial jets of high d/D are preferable. On the other hand, for localized cooling, the single jet (d/D = 0) performed the best. In general, for a given Reynolds number, a co-axial jet yields higher heat transfer coefficients than a single jet. Off-axis velocity peaks were encountered for the jets with d/D = 0.105. The measurements of turbulence intensity yielded values as high as 18%.     

Measurements of mean and turbulence quantities in the curved wake of an airfoil

Measurements of mean and turbulence quantities are presented for a curved wake of an airfoil. The wake is generated by placing a NACA 0012 airfoil of 0.150 m chord length at one chord length upstream of a 90° bend. The bend has a square cross-section of 0.457 m × 0.457 m, a mean radius-to-height ratio of R/H=1.17, and concave and convex radii of curvature 0.764 and 0.307 m, respectively. In addition to streamwise curvature, the wake is subjected to varying streamwise and radial pressure. The measurements were carried out at mainstream air velocities of 10, 15 and 20 m/s. The results are presented for the mean streamwise velocity, five components of turbulence stresses, the calculated wake half-width and the maximum velocity defect. The results showed the formation of an asymmetric wake about the wake centreline, with a larger wake half-width on the inner side. The wake half-width on both inner side and outer side of the wake decrease with mainstream velocity, whereas the maximum velocity defect, turbulence stresses increase with mainstream velocity. The turbulence stresses are enhanced on the inner side but suppressed on the outer side.     

Roughness effects on turbulent plane wall jets in an open channel

This paper reports the effects of surface roughness on the mean flow characteristics for a turbulent plane wall jet created in an open channel. The velocity measurements were obtained using a laser Doppler anemometer over smooth and transitionally rough surfaces. The power law proposed by George et al. (2000) was used to determine the friction velocity. Both conventional scaling and the momentum–viscosity scaling proposed by Narasimha et al. (1973) were used to analyze the streamwise evolution of the flow. The results show that surface roughness increases the skin friction coefficient and the inner layer thickness, but the jet half-width is nearly independent of surface roughness.     

Laminar,transitional and turbulent annular flow of drag-reducing polymer solutions

Mean and rms axial velocity-profile data obtained using laser Doppler anemometry are presented together with pressure-drop data for the flow through a concentric annulus (radius ratio κ = 0.506) of a Newtonian (a glycerine–water mixture) and non-Newtonian fluids—a semi-rigid shear-thinning polymer (a xanthan gum) and a polymer known to exhibit a yield stress (carbopol). A wider range of Reynolds numbers for the transitional flow regime is observed for the more shear-thinning fluids. In marked contrast to the Newtonian fluid, the higher shear stress on the inner wall compared to the outer wall does not lead to earlier transition for the non-Newtonian fluids where more turbulent activity is observed in the outer wall region. The mean axial velocity profiles show a slight shift (～5%) of the location of the maximum velocity towards the outer pipe wall within the transitional regime only for the Newtonian fluid.     

The effects of wall roughness on the separated flow over a smoothly contoured ramp

Experimental measurements address the effects on a turbulent boundary layer of wall roughness on a flat plate and a ramp that produces a separation bubble over the ramp trailing edge. A fully rough flow condition is achieved on the upstream flat plate. The main effect of the wall roughness on the outer layer turbulence on a flat plate is to change the friction velocity. The separation region is substantially larger for the rough-wall case. The rough-wall boundary layer turbulence is less sensitive to the onset of an adverse pressure gradient over the ramp, producing substantially smaller Reynolds stress peaks in upstream flat-plate, wall-unit coordinates.     

Turbulence modulation in a particle-laden flow over a hemisphere-roughened wall

Turbulence modulation by the inertia particles in a spatially developing turbulent boundary layer flow over a hemisphere-roughened wall was investigated using the direct numerical simulation method. The Eulerian and Lagrangian approaches were used for the gas- and particle-phases, respectively. An immersed boundary method was employed to resolve the hemispherical roughness element. The hemispheres were staggered in the downstream direction and arranged periodically in the streamwise and spanwise directions with spacing of p_x/d= 4 and p_z/d= 2 (where p_x and p_z are the streamwise and spanwise spacing of the hemispheres, and d is the diameter). The effects of particles on the turbulent coherent structures, turbulent statistics and quadrant events were analyzed. The results show that the addition of particles significantly damps the vortices structures and increases the length scales of streak structures. Compared with the particle-laden flow over the smooth wall, the existence of the wall roughness decreases the mean streamwise velocity in the near wall region, and makes the peaks of Reynolds stresses profiles shift up. In addition, the existence of particles also increases the percentage contributions to Reynolds shear stress from the Q4 events, however, decreases the percentage contributions from other quadrant events.     

Large Eddy Simulations of turbulent particle-laden channel flow

This paper scrutinises the Large Eddy Simulation (LES) approach to simulate the behaviour of inter-acting particles in a turbulent channel flow. A series of simulations that are fully (four-way), two-way and one-way coupled are performed in order to investigate the importance of the individual physical phenomena occurring in particle-laden flows. Moreover, the soft sphere and hard sphere models, which describe the interaction between colliding particles, are compared with each other and the drawbacks and advantages of each algorithm are discussed. Different models to describe the sub-grid scale stresses with LES are compared. Finally, simulations accounting for the rough walls of the channel are compared to simulations with smooth walls. The results of the simulations are discussed with the aid of the experimental data of Kussin J. and Sommerfeld M., 2002, Experimental studies on particle behaviour and turbulence modification in horizontal channel flow with different wall roughness, Exp. in Fluids, 33, pp. 143–159 of Reynolds number 42,000 based on the full channel height. The simulations are carried out in a three-dimensional domain of 0.175 m × 0.035 m  × 0.035 m where the direction of gravity is perpendicular to the flow. The simulation results demonstrate that rough walls and inter-particle collisions have an important effect in redistributing the particles across the channel, even for very dilute flows. A new roughness model is proposed which takes into account the fact that a collision in the soft sphere model is fully resolved and it is shown that the new model is in very good agreement with the available experimental data.     

Outer flow scaling of smooth and rough wall turbulent boundary layers

Mean velocity profiles in a zero pressure gradient turbulent boundary layer were measured on a hydraulically smooth surface and three different rough surfaces created from sand paper, perforated plate, and woven wire mesh. The physical size and geometry of the roughness elements were chosen to encompass both transitionally and fully rough flow regimes. The mean velocity profiles were measured using a Pitot tube in a subsonic wind tunnel, for Reynolds numbers (based on momentum thickness) ranging from 3,730 to 12,260. Three different outer velocity scales were used to analyze the defect profile. The results show that application of a so called mixed outer scale causes the velocity profile in the outer region to collapse onto the same curve for different Reynolds numbers and roughness conditions. Although the mixed scale collapses defect profiles on different surfaces, the effect of surface roughness is still observed in the outer region.     

The effect of sidewalls on rectangular jets

An experimental study is presented regarding the influence of sidewalls on the turbulent free jet flow issuing from a smoothly contracting rectangular nozzle of aspect ratio 15. “Sidewalls” are two parallel plates, flush with each of the slots’ short sides, practically establishing bounding walls extending the nozzle sidewalls in the downstream direction. Measurements of the streamwise and lateral velocity mean and turbulent characteristics have been accomplished, with an x-sensor hot wire anemometer, up to an axial distance of 35 nozzle widths, for jets with identical inlet conditions with and without sidewalls. Centreline measurements for both configurations have been collected for three Reynolds numbers, Re_D = 10,000, 20,000 and 30,000. For Re_D = 20,000 measurements in the transverse direction were collected at 13 different downstream locations in the range, x = 0–35 nozzle widths, and in the spanwise direction at three different downstream locations, x = 2, 6 and 25 nozzle widths.Results indicate that, the two jet configurations (with and without sidewalls) produce statistically different flow fields. Sidewalls do not lead to the production of a 2D flow field as undulations in the spanwise mean velocity distribution indicate. They do increase the two-dimensionality of the jet increasing the longevity of 2D spanwise rollers structures formed in the initial stages of entrainment, which are responsible for the convection of longitudinal momentum towards the outer field, establishing larger streamwise mean velocities at the jet edges. In the near field, up to 25 nozzle widths, lower outward lateral velocities in the presence of the sidewalls are held responsible for the decrease of turbulent terms including rms of velocity fluctuations and Reynolds stresses. Skewness factors increase monotonically across the shear layers from negative values to positive forming sharp peaks at the outer edges of the jet, illustrative of the presence of well defined 2D roller structures in the jet with sidewalls.     

On the impacts of coarse-scale models of realistic roughness on a forward-facing step turbulent flow

The present work explores the impacts of the coarse-scale models of realistic roughness on the turbulent boundary layers over forward-facing steps. The surface topographies of different scale resolutions were obtained from a novel multi-resolution analysis using discrete wavelet transform. PIV measurements are performed in the streamwise–wall-normal (x–y) planes at two different spanwise positions in turbulent boundary layers at Re_h = 3450 and δ/h = 8, where h is the mean step height and δ is the incoming boundary layer thickness. It was observed that large-scale but low-amplitude roughness scales had small effects on the forward-facing step turbulent flow. For the higher-resolution model of the roughness, the turbulence characteristics within 2h downstream of the steps are observed to be distinct from those over the original realistic rough step at a measurement position where the roughness profile possesses a positive slope immediately after the step’s front. On the other hand, much smaller differences exist in the flow characteristics at the other measurement position whose roughness profile possesses a negative slope following the step’s front.