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.     

Particle image velocimetry measurements of a backward-facing step flow

Particle image velocimetry (PIV) measurements were carried out on a backward-facing step flow at a Reynolds number of Re_h=U_Xh/9=4,660 (based on step height and freestream velocity). In-plane velocity, out-of-plane vorticity, Reynolds stress and turbulent kinetic energy production measurements in the x-y and x-z planes of the flow are presented. Proper orthogonal decomposition was performed on both the fluctuating velocity and vorticity fields of the x-y plane PIV data using the method of snapshots. Low-order representations of the instantaneous velocity fields were reconstructed using the velocity modes. These reconstructions provided insight into the contribution that the various length scales make to the spatial distribution of mean and turbulent flow quantities such as Reynolds stress and turbulent kinetic energy production. Large scales are found to contribute to the Reynolds stresses and turbulent kinetic energy production downstream of reattachment, while small scales contribute to the intense Reynolds stresses in the vicinity of reattachment.     

Measuring turbulence energy with PIV in a backward-facing step flow

Turbulence energy is estimated in a backward-facing step flow with three-component (3C, stereo) particle image velocimetry (PIV). Estimates of turbulence energy transport equation for convection, turbulence transport, turbulence production, viscous diffusion, and viscous dissipation in addition to Reynolds stresses are computed directly from PIV data. Almost all the turbulence energy terms in the backward-facing step case can be measured with 3C PIV, except the pressure-transport term, which is obtained by difference of the other turbulence energy terms. The effect of the velocity spatial sampling resolution in derivative estimations is investigated with four two-dimensional PIV measurement sets. This sampling resolution information is used to calibrate the turbulence energies estimated by 3C PIV measurements. The focus of this study is on the separated shear layer of the backward-facing step. The measurements with 3C PIV are carried out in a turbulent water flow at Reynolds number of about 15,000, based on the step height h and the inlet streamwise maximum mean velocity U₀. The expansion ratio (ER) is 1.5. Turbulence energy budget profiles in locations x/h=4, x/h=6, and x/h=10 are compared with DNS data of a turbulent flow. The shapes of profiles agree well with each other. Different ERs between the PIV case (1.5) and the DNS case (1.2) cause higher values for the turbulence energies measured by PIV than the energies by DNS when x/h=10 is approached. PIV results also show that the turbulence energy level in these experiments is generally higher than that of the DNS data.     

Three-dimensional particle-tracking velocimetry measurement of turbulence statistics and energy budget in a backward-facing step flow

Using a three-dimensional (3-D) particle-tracking velocimeter, detailed turbulent flow measurements were made in a plane channel with a one-sided 50% abrupt expansion, which acted as a backward-facing step. The turbulent channel flow reached a fully developed state well upstream of the step. The Reynolds number based on the upstream centerline velocity and the step height H was 5540. With the mean reattachment point located at 6.51H downstream of the step, the measurement region ranged from −2H upstream to 12H downstream of the step. Various turbulent statistics and the energy budget were calculated from numerous instantaneous vector distributions. As in previous experimental investigations, the Reynolds normal and shear stresses had maximum values upstream of the reattachment. The stress anisotropy tensor revealed a peculiar phenomenon near the reattachment wall, wherein the spanwise normal stress was the largest among the three normal stresses. The triple velocity correlations indicated large values in the separating shear layer, and hence the turbulent diffusion was a major term in the energy budget. Comparison was made between the present results and those of the direct numerical simulation (DNS) of Le et al. (1993), and it was found that the mean and fluctuating velocities, the Reynolds shear stress, and the turbulent energy budget were in excellent agreement, although there was a considerable difference in the inflow 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.     

Large eddy simulations on the effect of the irregular roughness shape on turbulent channel flows

Large Eddy Simulations (LES) are carried out to investigate on the mean flow in turbulent channel flows over irregular rough surfaces. Here the attention is focused to selectively investigate on the effect induced by crests or cavities of the roughness. The irregular shape is generated through the super-imposition of sinusoidal functions having random amplitude and four different wave-lengths. The irregular roughness profile is reproduced along the spanwise direction in order to obtain a 2D rough shape. The analysis of the mean velocity profiles shows that roughness crests induce higher effect in the outer-region whereas roughness cavities cause the highest effects in the inner-region with a reduced effect in the external region. The numerical simulations have been carried out at friction Reynolds number Re_τ=395. Similar results have been found for the higher order statistics: turbulence intensities or shear stresses. The analysis of the Reynolds stress anisotropy tensor confirms the existence of specific roles of cavities and crests in the turbulence modulation.     

Experimental Investigation of Three-Dimensional Vortex Structures Downstream of Vortex Generators Over a Backward-Facing Step

An experimental investigation of vortex generators has been carried out in turbulent backward-facing step (BFS) flow. The Reynolds number, based on a freestream velocity U₀ = 10 m/s and a step height h = 30 mm, was Re_h = 2.0 × 10⁴. Low-profile wedge-type vortex generators (VGs) were implemented on the horizontal surface upstream of the step. High-resolution planar particle image velocimetry (2D-2C PIV) was used to measure the separated shear layer, recirculation region and reattachment area downstream of the BFS in a single field of view. Besides, time-resolved tomographic particle image velocimetry (TR-Tomo-PIV) was also employed to measure the flow flied of the turbulent shear layer downstream of the BFS within a three-dimensional volume of 50 × 50 × 10 mm³ at a sampling frequency of 1 kHz. The flow control result shows that time-averaged reattachment length downstream of the BFS is reduced by 29.1 % due to the application of the VGs. Meanwhile, the Reynolds shear stress downstream of the VGs is considerably increased. Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) have been applied to the 3D velocity vector fields to analyze the complex vortex structures in the spatial and temporal approaches, respectively. A coherent bandwidth of Strouhal number 0.3 < St_h < 0.6 is found in the VG-induced vortices, and moreover, Λ-shaped three-dimensional vortex structures at St_h = 0.37 are revealed in the energy and dynamic approaches complementarily.     

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.     

Characterization of Turbulent Structures in a Transonic Backward-Facing Step Flow

A transonic backward-facing step flow, at a free stream Mach number of 0.8 and a Reynolds number of 1.86 × 10⁵ with respect to the step height, was investigated experimentally by means of planar and stereo Particle Image Velocimetry (PIV) measurements for multiple fields of view. The primary aim of this analysis is to examine whether the large temporal variations of the reattachment location is associated with the presence of large scale coherent flow structures. The mean flow reattaches ≈6.1±0.2 times the step height downstream of the step. This value fluctuates temporally as much as ±3 step heights. Measurements of the wake flow in horizontal planes show that the strong variations of the reattachment length are associated with spanwise variations of the streamwise velocity. Two-point correlations revealed large–scale coherent regions with a length of up to 7 step heights and a dominant spanwise wave-length of 1.5…2.5 step heights. Furthermore, close to the step large structures are found, which span more than 5 step heights in spanwise direction. The Reynolds stress distribution of the separated region strongly suggests that the initial streamwise momentum is transferred to the vertical component as well as to the spanwise component in comparable portions by the deformation of the initial Kelvin-Helmholtz vortices and the generation of secondary ones. As a result, the separated shear layer is characterized by eddies of various sizes and orientations. The mean flow field only shows the primary separation bubble and a secondary recirculation region. No stationary streamwise vortices could be found for the tested Reynolds number.     

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.     

DNS and LES of Passively Controlled Turbulent Backward-Facing Step Flow

A passive control approach (no external energy input) for an unsteady separated flow case was investigated numerically. A surface-mounted control fence was positioned upstream of a backward-facing step, and as an oncoming flow a thin and fully developed turbulent boundary layer with a thickness of δ/h = 0.8 was used. The objective of the passive control was to enhance the entrainment rate of the shear layer bounding the separation zone behind the step, thereby reducing the mean reattachment length,〈 X _r0 〉. Direct Numerical Simulations (DNS) and Large-Eddy Simulations (LES) at Re _h = 3000 (based on the step height, h, and the free stream velocity, U _∞) were carried out for the uncontrolled and the controlled flow case. The LES results were in good agreement with the DNS reference solutions. Adaptively controlled feedback simulations showed that a certain minimum distance between the step edge and the upstream position of the control fence is required to achieve a maximum reduction of the reattachment length.     

Development of a turbulent boundary layer after a step from smooth to rough surface

The flow developing downstream of a step change from smooth to rough surface condition is studied in the light of Townsend’s wall similarity hypothesis. Previous studies seem to support the hypothesis for channel and pipe flows, but there are considerable controversies about its application to boundary layers and in particular to surface roughness formed by spanwise bars. It has been suggested that this controversy arises from insufficient separation of scales between the boundary layer thickness and the roughness length scale. An experimental investigation has therefore been undertaken where the flow evolves from a fully developed smooth wall boundary layer at high Reynolds numbers over a step in surface roughness (Re _θ = 13,400 at the step). The flow is mapped through the development of the internal layer until the flow is fully developed over the rough wall. The internal layer is found to grow as δ ∼ X ^0.73, and after about 15 boundary layer thicknesses at the step, the internal layer has reached the outer edge of the incoming layer. At the last rough wall measurement station, the Reynolds number has grown to Re _θ ≈ 32,600 and the ratio of boundary layer to roughness length scales is δ/k ≈ 140. The outer layer differences between the smooth and the rough wall data were found to be sufficiently small to conclude that for this setup the Townsend’s wall similarity hypothesis appears to hold.     

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.     

Turbulence in rough-wall boundary layers: universality issues

Wind tunnel measurements of turbulent boundary layers over three-dimensional rough surfaces have been carried out to determine the critical roughness height beyond which the roughness affects the turbulence characteristics of the entire boundary layer. Experiments were performed on three types of surfaces, consisting of an urban type surface with square random height elements, a diamond-pattern wire mesh and a sand-paper type grit. The measurements were carried out over a momentum thickness Reynolds number (Re _θ) range of 1,300–28,000 using two-component Laser Doppler anemometry (LDA) and hot-wire anemometry (HWA). A wide range of the ratio of roughness element height h to boundary layer thickness δ was covered (0.04 ￡ h/d ￡ 0.400.04 \leq h/\delta \leq 0.40). The results confirm that the mean profiles for all the surfaces collapse well in velocity defect form up to surprisingly large values of h/δ, perhaps as large as 0.2, but with a somewhat larger outer layer wake strength than for smooth-wall flows, as previously found. At lower h/δ, at least up to 0.15, the Reynolds stresses for all surfaces show good agreement throughout the boundary layer, collapsing with smooth-wall results outside the near-wall region. With increasing h/δ, however, the turbulence above the near-wall region is gradually modified until the entire flow is affected. Quadrant analysis confirms that changes in the rough-wall boundary layers certainly exist but are confined to the near-wall region at low h/δ; for h/δ beyond about 0.2 the quadrant events show that the structural changes extend throughout much of the boundary layer. Taken together, the data suggest that above h/δ ≈ 0.15, the details of the roughness have a weak effect on how quickly (with rising h/δ) the turbulence structure in the outer flow ceases to conform to the classical boundary layer behaviour. The present results provide support for Townsend’s wall similarity hypothesis at low h/δ and also suggest that a single critical roughness height beyond which it fails does not exist. For fully rough flows, the data also confirm that mean flow and turbulence quantities are essentially independent of Re _θ; all the Reynolds stresses match those of smooth-wall flows at very high Re _θ. Nonetheless, there is a noticeable increase in stress contributions from strong sweep events in the near-wall region, even at quite low h/δ.     

Large-scale spanwise periodicity in a turbulent boundary layer induced by highly ordered and directional surface roughness

The effect of converging–diverging riblet-type surface roughness (riblets arranged in a ‘herringbone’ pattern) are investigated experimentally in a zero pressure gradient turbulent boundary layer. For this initial parametric investigation three different parameters of the surface roughness are analysed in detail; the converging–diverging riblet yaw angle α, the streamwise fetch or development length over the rough surface F_x and the viscous-scaled riblet height h⁺. It is observed that this highly directional surface roughness pattern induces a large-scale spanwise periodicity onto the boundary layer, resulting in a pronounced spanwise modification of the boundary layer thickness. Hot-wire measurements reveal that above the diverging region, the local mean velocity increases while the turbulent intensity decreases, resulting in a thinner overall boundary layer thickness in these locations. The opposite situation occurs over the converging region, where the local mean velocity is decreased and the turbulent intensity increases, producing a locally thicker boundary layer. Increasing the converging–diverging angle or the viscous-scaled riblet height results in stronger spanwise perturbations. For the strongest convergent–divergent angle, the spanwise variation of the boundary layer thickness between the diverging and converging region is almost a factor of two. Such a large variation is remarkable considering that the riblet height is only 1% of the unperturbed boundary layer thickness. Increasing the fetch seems to cause the perturbations to grow further from the surface, while the overall strength of the induced high and low speed regions remain relatively unaltered. Further analysis of the pre-multiplied energy spectra suggests that the surface roughness has modified or redistributed the largest scale energetic structures.     

Effects of the separating shear layer on the reattachment flow structure Part 1: Pressure and turbulence quantities

The effect of the separating shear-layer thickness and shape on the structure of the flow in the reattachment region of a backward-facing step is examined using wall static-pressure profiles and turbulence data for a range of Reynolds number (800 < Re _H< 40,000) and upstream boundary-layer thickness (0 < δ/H < 2). The reattachment pressure and the peak pressure in the reattachment zone decrease in a continuous manner as the upstream boundary layer thickens. The thinnest boundary layers follow the correlation of Roshko and Lau. Using the pressure data, correlations are developed which can be used to predict the level of turbulent shear stress in the near-wall region at reattachment, a location in which experimental data are extremely difficult to obtain.     

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.     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.     

Expansion ratio effects on the separated shear layer and reattachment downstream of a backward-facing step

Experiments were carried out to study the behavior of the incompressible turbulent separated shear layer and subsequent reattachment, downstream of a backward-facing step in a channel. The main objective of the study was to determine the effect of the expansion ratio on the development of the mean velocity and turbulence intensity in the shear layer and on the evolution of wall static pressure downstream of the step. The step height-to-upstream channel height ratio was varied between 0.5 and 2.13 while all inlet conditions were kept constant. Both hot-wire anemometry and frequency shifted laser Doppler anemometry were used for the velocity measurements. The Reynolds number based on free stream velocity and channel height upstream of the step was 16,600. The expansion ratio was found to have a particularly strong influence in the development of the turbulent, separated shear layer. Larger step height-to-inlet channel height ratios lead to higher turbulence intensities and faster growth of the unstable shear layer. As a result of this, shorter normalized reattachment lengths occurred with lager expansion ratios. For all the expansion ratios studied, the mean reattachment lenght was uniform along the spanwise direction except very near the side walls.     

Numerical Study of Turbulent Channel Flow over Surfaces with Variable Spanwise Heterogeneities: Topographically-driven Secondary Flows Affect Outer-layer Similarity of Turbulent Length Scales

Wall-bounded turbulent flows over surfaces with spanwise heterogeneous surface roughness – that is, spanwise-adjacent patches of relatively high and low roughness – exhibit mean flow phenomena entirely different to what would otherwise exist in the absence of spanwise heterogeneity. In the outer layer, mean counter-rotating rolls occupy the depth of the flow, and are positioned such that “upwelling” and “downwelling” occurs above the low and high roughness, respectively. It has been comprehensively shown that these secondary flows are Prandtl’s secondary flow of the second kind (Anderson et al., J. Fluid Mech. 768, 316–347 2015). This behaviour indicates that spanwise spacing, s _y , between adjacent patches of high and low roughness is, itself, a problem parameter; in this study, we have systematically assessed how s _y affects turbulence structure in high Reynolds number channel flows via two-point correlations. “High roughness” is imposed with streamwise-aligned pyramid elements with height, h, selected to be ≈ 5% of the channel half height, H. For \(s_{y}/H \gtrsim 1\), we find that the aforementioned domain-scale mean circulations exist and the surface may be regarded as a topography. For s _y /H ? 0.2, turbulence statistics show characteristics very similar to a homogeneous roughness and thus the surface may be regarded as a roughness. For 0.2 ? s _y /H ? 2, the spatial extent of the counter-rotating rolls is controlled by proximity to adjacent rows, and we define such surfaces as being intermediate. We refer to such surfaces as intermediate state.