4.1. Flow visualisations and vortex structures during the transient

Figure 5 shows 3D isosurface plots of low- and high-speed streaks, $\mp u^{\prime \prime }=(u-\langle \overline{u}\rangle )$ , and negative ${\it\lambda}_{2}$ at several instants of the transient flow. Note that, for the sake of clarity of the structures, only a quarter of the streamwise and a fifth of the spanwise extents of the simulated domain are shown. Here, ${\it\lambda}_{2}$ is the second largest eigenvalue of the symmetric tensor $\unicode[STIX]{x1D64E}^{2}+{\it\bf\Omega}^{2}$ , where $\unicode[STIX]{x1D64E}$ and ${\it\bf\Omega}$ are the symmetric and antisymmetric parts of the velocity gradient tensor, ${\rm\nabla}\unicode[STIX]{x1D66A}$ .

The negative isosurface of ${\it\lambda}_{2}$ is used to identify vortex cores (Jeong & Hussain Reference Jeong and Hussain1995). Note that all visualised quantities are normalised by wall units of the initial flow. It is seen that, following the step increase in flow rate, the flow undergoes a distinctive transient phase. During this process, a well-organised roughness-induced, strong, head-up hairpin structure is initially generated around each roughness element. These vortices, which are similar to each other in both size and strength over all elements in the flow domain, are then convected downstream in the flow direction and evolve in the subsequent instants. In the following, the birth and development of such vortices over the roughness elements are explained.

Referring to figure 5, soon after the commencement of the transient ( $t^{\ast }=0.09$ ), a region of strong vortex (shown by ${\it\lambda}_{2}$ isosurfaces) develops around the roughness crest extending along the spanwise ridge line. Later ( $t^{\ast }=0.2$ ), the vortex is split into two parts. One part remains at the location of its birth (that is, at the crest of the roughness element) throughout the rest of the transient process (see later instants). The second part is convected downstream and plays a major role in the transition. The lower parts of this new vortex traverse alongside the two trough lines, and the upper necklace part located behind the roughness crest goes downstream along the ridge line. At $t^{\ast }=0.25$ , a typical primary hairpin vortex has emerged, with two counter-rotating legs alongside the ridge lines and the head in the roughness wake (around the trough). This vortex has travelled half of the roughness wavelength ( $x/{\it\lambda}=0.5$ ) by this instant. At this stage, the vortices are still very close to the surface. By $t^{\ast }\sim 0.4$ , the head of the quasi-streamwise primary vortex has reached the crest of the successive roughness element, and starts to move away from the wall.

Figure 6(a) shows a close-up view of the hairpin structures at this instant. It is also seen that the region of the strong hairpin vortices is associated with the low-pressure region. From approximately $t^{\ast }=0.8$ , a significant burst of the primary hairpins is observed and the vortices are replaced by more random and smaller vortical structures. By $t^{\ast }\sim 4$ , there are no strong structures identified near the roughness and the flow in the roughness sublayer has effectively reached the final flow condition. This trend is consistent with the development of the turbulence quantities (e.g. $U_{m}^{+}$ and $u^{\prime }v^{\prime }$ ) discussed in the following sections, in which, by $t^{\ast }\sim 4$ , the inner layer ( $y^{\prime +}<100$ ) has reached the final flow values.

The formation of streaks is another distinctive feature of the transition process. As shown in figure 5, at a very early stage ( $t^{\ast }=0.25$ ), positive streaks (green structures) are formed, which go along the troughs of the roughness elements. At this time, patches of negative velocity are formed within the primary vortices. At a later time ( $t^{\ast }=0.4$ ), the negative velocity streaks are formed as such patches are connected. This is also the time when the heads of the primary vortices move away from the wall (see side-view plots of figure 6 b–g). This creates a situation where all the primary vortices sit on top of the negative streaks, which is a typical scenario identified in bypass transition (Jacobs & Durbin Reference Jacobs and Durbin2001; Nolan, Walsh & Mceligot Reference Nolan, Walsh and Mceligot2010). In fact the positive and negative streaks are themselves supporting evidence of the existence of a strong counter-rotating vortex (Alfonsi Reference Alfonsi2006). Note that for figure 5 an isosurface level that better depicts the streaky structures during the transient flow was chosen and applied to all instants; as a consequence, the streaky structures at $t^{\ast }=0.09$ , which exist due to initial turbulence, are not visible.

Figure 6. Primary vortex structure and pressure. (a) Close-up view of vortex structure ( $-{\it\lambda}_{2}$ ) and negative pressure at $t^{\ast }=0.4$ . (b–g) Isosurface of ${\it\lambda}_{2}$ of figure 5 shown in side view ( $x$ – $y$ plane): (b)  $t^{\ast }=0.09$ , (c)  $t^{\ast }=0.25$ , (d)  $t^{\ast }=0.42$ , (e)  $t^{\ast }=0.6$ , (f)  $t^{\ast }=0.8$ and (g)  $t^{\ast }=4$ . For clarity, the side-view plots are shown for half a roughness spanwise wavelength. The dashed lines in the side-view plots show the location of $y/k_{t}=3$ .

The existence of standing and hairpin vortices, and the evolution and breakdown of the latter in the transient flow, bears a strong qualitative resemblance to those of the roughness-induced laminar–turbulent transition described by Acarlar & Smith (Reference Acarlar and Smith1987) for a low-speed Blasius boundary layer and by De Tullio et al. (Reference De Tullio, Paredes, Sandham and Theofilis2013) for a supersonic boundary layer for an isolated roughness element. This is particularly interesting considering that, in the present case, (i) the initial flow is a turbulent flow in the transitionally rough regime ( $k_{s}^{+}\sim 15$ ), and (ii) the transient process occurs over a surface that is entirely covered by close-packed distributed roughness elements in both the streamwise and spanwise directions. It is shown later (e.g. refer to figure 9 c) that, in the transient flow considered herein, a boundary layer is formed as a result of the rapid increase of Reynolds number. This creates a condition similar to that in a spatially developing boundary layer, resulting in the generation of vortices and later their breakdown.

Figure 7. Vector plot and contour of ${\it\lambda}_{2}$ (negative values are shaded) at $x/{\it\delta}=1.503$ ( $x/{\it\lambda}=0.005$ ).

To investigate further the generation and evolution of the strong primary hairpin vortices during the early stage of the transient process, the development of the velocity fields ( $v,w$ ) at several instants is shown in figure 7 in the form of vector plots in a vertical $y$ – $z$ plane across two roughness elements. The plane is chosen at $x/{\it\delta}=1.503$ ( $x/{\it\lambda}=0.005$ ), which is very close to the crest of the forward side of the element. The vector frames are shaded with contours of the negative ${\it\lambda}_{2}$ , which is normalised by $(u_{{\it\tau}_{0}}/{\it\delta})^{2}$ . Very shortly after the commencement of the transient ( $t^{\ast }=0.09$ ), upward activities (although relatively weak) are observed to occupy much of the surface of the roughness, which can be identified from the velocity vectors. The second frame ( $t^{\ast }=0.3$ ) corresponds to an instant close to when the primary hairpin structures are formed. It can be seen from the velocity vectors and the strong negative values of ${\it\lambda}_{2}$ that a pair of counter-rotating vortices is present within each trough. These are signatures of the hairpin legs. The vortices are centred in a low-speed streak and are surrounded by high-speed streaks (not shown here). At $t^{\ast }=0.5$ , the signature of a hairpin head can be observed just above the trough. Here the velocity vectors go outwards at an elevated level and the low-speed streak starts to lift from the wall. The fourth frame ( $t^{\ast }=0.8$ ) shows the vectors for the time when the breakdown of the organised hairpins has been started, and at $t^{\ast }=1.5$ (fifth frame) most such structures near the roughness have already been pumped and broken up. At a later time ( $t^{\ast }\sim 4$ ) the organised structures near the roughness have completely disappeared.

It is useful to note that, during the first stage ( $t^{\ast }<0.8$ ), signatures of the primary hairpin’s heads and legs (exhibited in figure 7 as strong vortical structures) are confined up to $y/k_{t}\sim 3$ above the roughness crest. This is also confirmed by the development of energy spectra discussed in § 4.6, where the strong spikes exhibited at the wavenumbers corresponding to the roughness wavelength and its multiples largely vanish for locations above $y/k_{t}\sim 3$ . In studies for which the flow field is available (e.g. using DNS or particle image velocimetry), the height of the roughness sublayer is normally defined as the distance at which the spatial variation of Reynolds stresses is diminished (Bhaganagar, Kim & Coleman Reference Bhaganagar, Kim and Coleman2004; Hong et al. Reference Hong, Katz and Schultz2011; Lee, Sung & Krogstad Reference Lee, Sung and Krogstad2011). This sublayer has been found to be 2–5 times the roughness height. Another definition of the roughness sublayer thickness used in the literature is the distance beyond which the statistical values of rough and smooth wall start to collapse (Bakken et al. Reference Bakken, Krogstad, Ashrafian and Andersson2005; Schultz & Flack Reference Schultz and Flack2005). Following the above methods, the present steady-flow results show that, for both the low- and high-Reynolds-number cases, the roughness sublayer is approximately $y\sim 3k_{t}$ above the crest.

Overall, the observations from figures 5 to 7 demonstrate that the transient process following the step increase of flow rate is characterised by a rather orderly process of one cycle of the birth of uniformly distributed quasi-streamwise primary hairpin vortices, followed by the evolution and eventual breakdown of these vortices. The number of such hairpin vortices is equal to the number of roughness elements. At a very early stage (say around $t^{\ast }=0.1$ ), a vortex is formed on the crest of each of the roughness elements, and it later splits into two ( $t^{\ast }=0.2$ ): a ‘standing’ vortex, which remains locally and largely inactive during later development, and a ‘primary hairpin’, which is convected downstream and forms a primary hairpin vortex with two counter-rotating legs joining at the head, at ( $t^{\ast }\sim 0.25$ ). Later, the primary vortices are convected downstream and their heads move away from the wall. Alongside this, two sets of streaks are generated within the span of every roughness element. The low-speed streaks sit below the hairpin vortices, which are accompanied by strong ejection activities. At $t^{\ast }\sim 0.8$ , a significant breakup of the vortices has started and eventually the orderly vortical structures are replaced by more chaotic structures that are typical of turbulent flows over a rough wall. It will take some time for the newly generated turbulence to propagate to the core of the flow, as will be demonstrated later.

4.2. Skin friction factor

The skin friction factor coefficient, $C_{f}$ , is defined as

(4.1) $$\begin{eqnarray}C_{f}=\frac{{\it\tau}_{w}}{\frac{1}{2}{\it\rho}U_{b}^{2}},\end{eqnarray}$$

where ${\it\tau}_{w}$ is the wall shear stress, denoted by ${\it\tau}_{ws}$ and ${\it\tau}_{wr}$ for the upper (smooth) and lower (rough) walls, respectively. For the smooth wall, the calculation of ${\it\tau}_{ws}$ is straightforward and can be obtained from the mean velocity profile,

(4.2) $$\begin{eqnarray}{\it\tau}_{ws}={\it\mu}\left.\frac{\partial \langle \overline{u}\rangle }{\partial y}\right|_{y=y_{ws}}.\end{eqnarray}$$

For the rough wall, shear stress contains the contribution of both pressure drag ( $P_{d}$ ) and friction drag ( ${\it\tau}_{f}$ ), and is calculated using ${\it\tau}_{ws}$ and the mean pressure gradient of the channel ( ${\it\Pi}$ in (2.1)) via the expression (Seddighi Reference Seddighi2011)

(4.3) $$\begin{eqnarray}{\it\tau}_{wr}=(P_{d}+{\it\tau}_{f})=\frac{{\it\Pi}S_{e}}{L_{x}}-{\it\tau}_{ws},\end{eqnarray}$$

where $S_{e}$ is the total flow area. In fact, it is the total area minus that associated with the roughness elements; so for the smooth case it is $S_{e}=(2{\it\delta})L_{x}$ . For the rough case, the height of the channel ( $2{\it\delta}$ ) varies slightly due to the roughness and an average is used. Leonardi, Orlandi & Antonia (Reference Leonardi, Orlandi and Antonia2005) proposed a method that is similarly based on the global momentum balance formula for the total pressure drop. The wall shear stress for each wall is determined from the total pressure loss and the zero-crossing in the wall-normal profile of the turbulent shear stress:

(4.4a,b ) $$\begin{eqnarray}{\it\tau}_{wr}=\frac{{\it\Pi}(2{\it\delta})}{(1+d_{1}/d_{2})},\quad {\it\tau}_{ws}=\frac{{\it\Pi}(2{\it\delta})}{(1+d_{2}/d_{1})},\end{eqnarray}$$

where $d_{1}$ and $d_{2}$ are the distances of the zero-crossing from the upper and lower wall, respectively. For the present results, both methods are used and are found to give similar values with a difference of less than 1 %. Friction factors for the top smooth wall obtained for the steady-state rough-case simulations are in close agreement with those from the literature (Monty Reference Monty2005).

Figure 8 shows the development of the friction coefficient for the transient flow. Also shown are the results of HS2013 for a channel with all smooth walls (smooth case). It is seen that the development of the friction coefficient for the smooth wall is consistent with the results of HS2013, exhibiting a bypass transition behaviour as discussed in § 1. A large undershoot, until $t^{\ast }\sim 8$ , exhibited by the smooth wall is characteristic of the flow in the pre-transition stage during which the flow exhibits a trend of a laminar friction factor (He & Seddighi Reference He and Seddighi2015). The friction factor for the rough wall, however, shows a trend that is significantly different from that of the smooth wall. The development shows neither a large undershoot at the initial stages, nor a long delay (until $t^{\ast }\sim 16$ ) before reaching the corresponding final flow value. This behaviour is not a surprise, noting that strong hairpin vortices are generated very early in time ( $t^{\ast }\sim 0.25$ ), and the transition is already completed by $t^{\ast }\sim 4$ . The variation of $C_{f}$ is largely dominated by the relaxation of the mean velocity profile, which causes a continual reduction from its initial peak value resulting from the sudden increase in flow rate. The effect of the primary vortex generated at $t^{\ast }\sim 0.25$ is reflected as a kink in the variation of $C_{f}$ (see the inset). In the period of $0.25<t^{\ast }<0.8$ , the profile shows an oscillatory variation with several further kinks. The period corresponds to the time during which the primary vortices are pumping up from the wall, whilst convecting downstream of the flow. The vortices remain in $y/k_{t}<1$ until $t^{\ast }\sim 0.8$ , when the major breakup of the primary vortices has started to occur. It is of particular interest to note that, in contrast to the development of the transient flow over a smooth wall in the corresponding period, the $C_{f}$ – $t^{\ast }$ curve for the rough wall is convex during the early transition period ( $t^{\ast }\sim 1.5$ ). This reflects the fact that the flow during this period is not a simple diffusion process (as in a laminar flow), but is also affected by the generation and evolution of vortices, thereby causing a change in effective diffusivity and a higher friction coefficient.

Figure 8. Development of the skin friction coefficient, $C_{f}$ , for the unsteady case; symbols show the smooth-wall results of He & Seddighi (Reference He and Seddighi2013).

Figure 9. Development of mean velocity profiles at several instants: (a) smooth wall; (b) rough wall. (c) Plot of $\overline{u}^{\wedge }$ (the profiles are shifted by an amount proportional to the time instant): thin lines, smooth wall; thick lines, rough wall. Lines 1 and 2 show the thickness of the boundary layer of $\overline{u}^{\wedge }$ , ${\it\delta}_{99}$ , for the rough and smooth walls, respectively. The crest’s level is shown by the short-dashed line. (d) Plot of contour lines of mean velocity $\overline{u}$ , shaded by $\overline{u^{\prime }v^{\prime }}/u_{{\it\tau}0}^{2}$ , for the $y$ – $z$ plane at $x/{\it\lambda}=0$ and $t^{\ast }=0.8$ .

4.3. Mean velocity, r.m.s. of velocity fluctuations

Figure 9 shows the mean velocity profiles for the unsteady flow at several instants. The velocity and effective distance, normalised by inner scaling, are shown in figure 9(a,b) for the smooth- and rough-wall cases, respectively. The development for the smooth wall shows a trend similar to that of the all-smooth channel of HS2013. After an initial reduction ( $t^{\ast }<0.02$ ), the velocity ( $\langle \overline{u}\rangle ^{+}$ ) increases monotonically until the onset of transition ( $t^{\ast }\sim 8$ ). During this period, the extent of the viscous sublayer (approximated by $\langle \overline{u}\rangle ^{+}=y^{\prime +}$ ) increases with time. After the transition begins, the velocity in the core region progressively reduces, reaching a profile close to that in a steady flow at the end of the transition ( $t^{\ast }\sim 16$ ). Later, whilst the profile in the inner region ( ${\sim}y^{\prime +}<100$ ) remains unchanged, the velocity in the wake region ( ${\sim}y^{\prime +}>100$ ) adjusts itself slightly to reach the corresponding final steady values. The velocity over the rough wall, however, evolves in a different manner. After the initial rapid reduction, the values increase monotonically until $t^{\ast }\sim 4$ . By this time, the inner layer ( ${\sim}y^{\prime +}<100$ ) reaches its final profile, remaining almost unchanged thereafter. As discussed in § 4.2 (figure 8), by $t^{\ast }\sim 4$ , the friction factor has almost reached the value corresponding to the final flow. In the outer layer ( $y^{\prime +}>100$ ), the velocity continues to adjust while oscillating around the final profile until approximately $t^{\ast }=16$ . Interestingly, the above trend of velocity development is similar to that reported by De Tullio et al. (Reference De Tullio, Paredes, Sandham and Theofilis2013, figure 18b) in the wake of an isolated roughness in a supersonic boundary layer where laminar–turbulent transition occurs.

Following HS2013, the time-developing boundary layer can be studied by examining the perturbation velocity,

(4.5) $$\begin{eqnarray}\overline{u}^{\wedge }(y,t^{\ast })=\left[\frac{\langle \overline{u}(y^{\prime },t^{\ast })\rangle -\langle \overline{u}(y^{\prime },0)\rangle }{\langle \overline{u}_{c}(t^{\ast })\rangle -\langle \overline{u}_{c}(0)\rangle }\right],\end{eqnarray}$$

where $y^{\prime }=y+{\it\epsilon}_{d}$ for the rough wall and $y^{\prime }$ is the wall-normal distance from the top wall for the smooth wall; and $\langle \overline{u}_{c}\rangle$ is the mean velocity corresponding to the effective centreline for each wall. Figure 9(c) shows such profiles against elevation of $y^{\prime }$ . For the smooth wall, the perturbation velocity profile is initially uniform along the $y^{\prime }$ axis, exhibiting a ‘plug-flow’ behaviour except in the vicinity of the wall, where a thin boundary layer is formed due to the ‘no-slip’ boundary condition on the wall. The boundary layer gradually develops into the core of the flow with time. The $\overline{u}^{\wedge }$ profiles at the various time instants are practically indistinguishable from those in the corresponding flow in a wholly smooth-wall channel of He & Seddighi (Reference He and Seddighi2015), which in turn agrees closely with the Stokes solution for a laminar flow during the pre-transition period ( $t^{\ast }<8$ ). The latter two data are not shown to avoid over-crowding the figure.

On the rough wall, initially (say, $t^{\ast }<0.25$ ), the flow above the roughness elements again responds like a ‘plug flow’, i.e. the perturbation flow $\overline{u}^{\wedge }$ is uniform along the $y^{\prime }$ axis. The mean velocity within the roughness varies linearly from the value above the element height to zero at the bottom of the roughness elements. During this period ( $t^{\ast }<0.25$ ), the turbulent stresses ( $u_{rms}^{\prime },v_{rms}^{\prime },w_{rms}^{\prime }$ ) above the roughness crest remain largely unchanged (thin lines with symbols in the insets of figure 10 g,i,k). Consequently, at this stage, the flow above the roughness can be viewed as comprising the initial turbulent flow plus a ‘plug flow’ (i.e. the perturbation flow), the latter of which contains only a mean velocity with no fluctuating components. During the next period ( $0.25<t^{\ast }<0.8$ ), the boundary layer that is initially confined within the roughness develops into the core at a rate that is significantly faster than that of the boundary layer over the smooth wall (see the development of the ${\it\delta}_{99}$ profiles in figure 9 c). This period coincides with the growths and lift-ups of the hairpin vortices from the wall. The boundary layer continues developing after $t^{\ast }=0.8$ , but, interestingly, a plateau develops in the profile just above the crest of the roughness elements. This coincides with the breakdown of primary vortices into conventional turbulence in the same region, which clearly provides strong mixing, leading to a flat velocity profile in that region. A similar behaviour was shown with the results of Acarlar & Smith (Reference Acarlar and Smith1987) and Muppidi & Mahesh (Reference Muppidi and Mahesh2012) for roughness-induced transition over, respectively, isolated and distributed roughness.

The perturbation boundary layer can, in principle, be seen as an internal boundary layer (IBL) developing within the initial wall shear flow in the channel. In this regard, it is noted that the concept of IBLs developing within an external boundary layer has been studied extensively and particularly in the context of meteorology (e.g. Smits & Wood Reference Smits and Wood1985; Garratt Reference Garratt1990; Saito & Pullin Reference Saito and Pullin2014). The IBLs often result from a discontinuity in surface roughness (smooth to rough or vice versa) or surface temperature, with the free-stream velocity remaining unchanged. The perturbation introduced due to such a discontinuity is likely to be much milder than the sudden large increase in the mean velocity concerned here. In addition, previous research in the literature seems to focus on the steady development of the IBL after the initial transition, when the boundary layer growth indicates a turbulent power-law behaviour. In contrast, our main interest is the initial transition. Consequently, direct comparison with IBL results from the literature is not attempted herein.

Figure 9(d) shows the line contour plot of ensemble-averaged mean velocity $\overline{u}$ for the $y$ – $z$ plane at $x/{\it\lambda}=0$ , at $t^{\ast }=0.8$ . The plot is shaded by $\overline{u^{\prime }v^{\prime }}$ defined in § 4.4. A clear hump is seen for the contour line of the mean velocity in the region of the strong turbulent shear stress, which coincides with where the hairpin vortices exist. This behaviour is also similar to that of Drews et al. (Reference Drews, Downs, Doolittle, Goldstein and White2011) for the roughness-induced transition behind the roughness. The authors attributed the hump to the strong shear layer due to the existence of the hairpin vortices.

Figure 10. Development of r.m.s. of fluctuating velocities for smooth wall (a–f) and rough wall (g–l) at several locations: (a,c,e,g,i,k), near wall; (b,d,f,h,j,l), core region.

Figure 10 shows the development of r.m.s. of velocity fluctuations normalised by the initial friction velocity ( $u_{{\it\tau}_{0}}$ ), at several near-wall and core-region locations. Results for the smooth wall shown in figure 10(a–f) exhibit a bypass transition behaviour. In the wall region (figure 10 a,c,e), $w_{rms}^{\prime }$ are almost unchanged and even reduce slightly during pre-transition ( $t^{\ast }\lesssim 8$ ). In contrast, the streamwise component increases progressively owing to the formation of streaks (HS2013). During the transition, $u_{rms}^{\prime }$ continues to increase, reaching a maximum, from which it reduces and approaches final steady values towards the end of transition ( $t^{\ast }\lesssim 16$ ). Meanwhile, the other two components, $v_{rms}^{\prime }$ and $w_{rms}^{\prime }$ , increase monotonically, reaching the corresponding final steady values towards the end of transition. In the core region (figure 10 b,d,f), there are long delays in the responses of all three components. Data for the smooth case of HS2013 are shown (by the symbols) in figure 10(a–f) for direct comparison with the results for the smooth wall from the present study (shown by the lines). It can be seen that there is good agreement between the results of HS2013 and those for the smooth wall of the present simulation. In particular, the energy growth during the pre-transition period and the point of onset of transition are very close in the two cases.

This indicates that the flow over one wall of the channel, whether it is smooth or rough, has little influence on the transient flow behaviour on the opposite wall. That is, the local flows over the two sides of the channel are largely independent. This is in contrast to the results of Orlandi (Reference Orlandi2009) on laminar–turbulent transition, which showed that the transition is affected by the presence of the roughness on the opposite wall. Orlandi’s observation is not difficult to understand. The early transition on the rough surface, which generates strong and violent velocity and pressure fluctuations, can potentially induce earlier transition on the smooth wall. The turbulence generated near the rough wall will propagate towards the smooth wall, and hence increase the free-stream turbulence, which could in turn cause an early transition on the smooth wall. However, from figure 11 it is seen that the new turbulence resulting from the early transition over the rough surface, in the present study, is initially confined to the region close to the wall. It propagates relatively slowly to the core of the flow. The wave front of increased turbulence on the rough wall only reaches the location where it might affect the transition after the onset of transition on the smooth wall. This is likely to be the reason that the transition on the smooth wall in the present transient flow is unaffected.

Figure 10(g–l) shows the development of r.m.s. of turbulent–spatial velocity fluctuations, $u_{s,rms}^{\prime \prime }$ , for the rough wall. Also shown in the plots are the corresponding data of r.m.s. of turbulent velocity fluctuations, $u_{s,rms}^{\prime }$ . The former measures the fluctuating velocities due to turbulence and the spatial variability due to the presence of the roughness elements. The latter accounts for the turbulent fluctuations only. It is shown that, whereas the r.m.s. of velocity fluctuations obtained from the two averaging approaches collapse on top of each other in the core region (figure 10 h,j,l), there are major differences between them, in particular for the streamwise and wall-normal components, in the region close to the roughness (figure 10 g,i,k). These results and those from other locations (not shown) show that the discrepancies vanish at above $y/k_{t}\sim 1.3$ , $y^{\prime +0}\sim 19$ . This indicates that significant spatial flow variations due to the roughness geometry remain very close to the element.

Figure 11. Development of shear stress: (a) smooth wall, profiles of $\langle \overline{u^{\prime }v^{\prime }}\rangle$ at several instants; (b,c) rough wall, profiles of, respectively, $\langle \overline{u^{\prime \prime }v^{\prime \prime }}\rangle$ and $\langle \overline{u^{\prime }v^{\prime }}\rangle$ at several instants; (d) rough wall, time variation of $\langle \overline{u^{\prime \prime }v^{\prime \prime }}\rangle$ and $\langle \overline{u^{\prime }v^{\prime }}\rangle$ at several locations close to the rough wall.

The development of the r.m.s. of velocity fluctuations (both $u_{s,rms}^{\prime }$ and $u_{s,rms}^{\prime \prime }$ ) on the rough wall in the core region (figure 10 h,j,l) is broadly similar to that on the smooth wall discussed above, except that the latter has a longer delay. For the rough wall, the response at the outer edge of the near-wall region (say $y^{\prime +0}=54.6$ ) occurs much earlier and the delay increases significantly from that location to the centre of the flow. The behaviour is very different in the near-wall region. Consider the r.m.s. of turbulent–spatial velocity fluctuations (thick lines) first. All three components begin to respond almost immediately after the commencement of the transient, although there are clear differences between them. The $u_{rms}^{\prime \prime }$ at the crest ( $y^{\prime +0}=6.1$ ) increases rapidly to reach a peak at $t^{\ast }=0.6$ , falling back afterwards, approaching its steady value at around $t^{\ast }=4$ . The $u_{rms}^{\prime \prime }$ at $y^{\prime +0}=11.6$ behaves in a similar way but with a small delay and a lower peak. Further from the wall ( $y^{\prime +0}=23.3$ ), there is an additional delay and no peak is seen. The $v_{rms}^{\prime \prime }$ is similar in terms of the response time, but the peaks at $y^{\prime +0}=6.1$ and 11.6 are much smaller. In addition, there are strong oscillatory variations at the early times. The increase of $w_{rms}^{\prime \prime }$ is much slower and there is no overshooting. Comparing the formation and growth of the streaks shown in figure 5 with the development of the r.m.s. of streamwise velocity fluctuations, it is concluded that the rapid increase and strong overshooting of $u_{rms}^{\prime \prime }$ in the early stage of the transient are largely associated with the formation of strong positive and then negative streaks. The formation of the hairpin vortices also contributes to the development of $u_{rms}^{\prime \prime }$ , but this has lesser influence. In contrast, the developments of $v_{rms}^{\prime \prime }$ and $w_{rms}^{\prime \prime }$ are influenced quite strongly by the formation and evolution of the hairpin structures. In particular, these cause the oscillatory behaviour in $v_{rms}^{\prime \prime }$ and $w_{rms}^{\prime \prime }$ . The early responses of the r.m.s. of turbulent velocity fluctuations (thin lines with symbols) are very different from those described above for the turbulent–spatial velocity fluctuations. All three components appear to respond in a similar way, all with a short delay after the commencement of the transient flow. Most significantly, the early rapid overshooting in $u_{s,rms}^{\prime \prime }$ is replaced by more gradual changes in $u_{s,rms}^{\prime }$ . The differences between them reflect the spatial variation of the velocities caused by the roughness topology. These are very strong at the early stages due to the formation of the streaks and vortices directly associated with the roughness elements.

As discussed above, in contrast to the smooth wall in which the wall-normal and spanwise components of the r.m.s. of velocity fluctuations remain almost unchanged for a long time in the pre-transition stage, for the rough wall these components exhibit a response similar to that of the streamwise component, increasing shortly after the commencement of the transient flow. This behaviour is in agreement with the results of the earlier studies on the rough wall, in which the increased wall-normal fluctuation has been discussed as the component that is strongly associated with the near-wall vortical structures, thereby characterising flow behaviour over the rough wall (Orlandi Reference Orlandi2013). For the transition from laminar to turbulent flow over a rough wall, Orlandi (Reference Orlandi2013) showed that the transition occurs once the r.m.s. of wall-normal velocity fluctuations (normalised by instantaneous friction velocity) at the roughness crest is greater than a threshold value of 0.6.

4.4. Turbulent shear stress and quadrant analysis

Figure 11 shows the development of turbulent shear stress at several instants. The turbulent shear stress for the smooth wall (figure 11 a) remains largely unchanged in the early moments of the pre-transition period ( $t^{\ast }<1.6$ ), then increases slowly until the onset of transition on the smooth wall ( $t^{\ast }\sim 8$ ). The response becomes much faster during the transition ( $8<t^{\ast }<16$ ), approaching the corresponding final steady-flow profile in the wall region; and eventually, it reaches a fully developed shape by $t^{\ast }\sim 20$ . The response starts from the wall region and propagates to the core with time.

Figure 12. Development of ensemble-averaged turbulent shear stress (normalised by $u_{{\it\tau}0}^{2}$ ) for an element of figure 2(b) at the roughness crest, $y=0$ , and at several instants: (a)  $t^{\ast }=0.2$ , (b)  $t^{\ast }=0.3$ , (c)  $t^{\ast }=0.4$ , (d)  $t^{\ast }=0.6$ , (e)  $t^{\ast }=1.5$ , (f)  $t^{\ast }=4$ . Dashed lines show the location of $x/{\it\lambda}=0.35$ .

The spatially averaged turbulent-spatial shear stress ( $\langle \overline{u^{\prime \prime }v^{\prime \prime }}\rangle$ ) in the wall region of the rough wall (figure 11 b) initially increases rapidly, overshoots the final steady values, and reaches a peak at $t^{\ast }\sim 0.3$ (see also history plot of figure 11 d). Then, during $0.3<t^{\ast }<0.8$ , the turbulent shear stress profile oscillates strongly before reducing monotonically with time. By $t^{\ast }\sim 4$ , it is approaching its final state in the region of $y^{\prime +0}<30$ (which corresponds to $y^{\prime +}<90$ in the mean velocity profiles of figure 9). The large overshoot observed in the near-wall region of the profile is associated with the formation of the strong well-organised roughness-induced hairpin vortices observed in figure 5. Such structures are initially significantly different from ‘conventional’ turbulence. Later, at late pre-transition stage, the structures begin to break down, leading to the generation of more conventional turbulence.

The development of $\langle \overline{u^{\prime }v^{\prime }}\rangle$ is shown in figure 11(c). In comparison with the spatially averaged turbulent-spatial shear stress (figure 11 b), the wall-region response of $\langle \overline{u^{\prime }v^{\prime }}\rangle$ is much slower and the large overshooting over the final steady values at the early moments of the transition shown in $\langle \overline{u^{\prime \prime }v^{\prime \prime }}\rangle$ (figure 11 b) are not seen here. It starts to respond at $y^{\prime +0}\sim 5$ (just below the crest), at $t^{\ast }\sim 0.25$ , developing a peak that increases and moves towards the core, reaching $y^{\prime +0}\sim 15$ at $t^{\ast }\sim 0.4$ . Between $0.4<t^{\ast }<1.5$ , the peak value increases while its $y$ location remains largely unchanged. Thereafter, the peak value reduces slightly but the values increase further away from the wall ( $y^{\prime +0}\sim 20$ ) until $t^{\ast }\sim 4$ , when, similar to $\langle \overline{u^{\prime \prime }v^{\prime \prime }}\rangle$ , the near-wall values approach the corresponding final steady-state flow, and the core region values reach the corresponding final-state values with further delays.

Figure 11(d) shows the time variation of both turbulent and turbulent-spatial shear stresses at several locations close to the rough wall. Strong oscillations in $\langle \overline{u^{\prime \prime }v^{\prime \prime }}\rangle$ can be seen clearly at the crest during the initial stages (until $t^{\ast }\sim 0.8$ ). Such strong oscillations, which are similar to those observed for $v_{rms}^{\prime \prime }$ and, to a lesser extent, for $w_{rms}^{\prime \prime }$ in figure 10(g,i,k), are absent in $\langle \overline{u^{\prime }v^{\prime }}\rangle$ . The spatially averaged turbulent-spatial shear stress, $\langle \overline{u^{\prime \prime }v^{\prime \prime }}\rangle$ , in the locations very close to the roughness, is the sum of the spatially averaged turbulent shear stress, $\langle \overline{u^{\prime }v^{\prime }}\rangle$ , and a so-called dispersive (or form-induced) stress, which results from the roughness geometry. Whereas a significant discrepancy is exhibited between $\langle \overline{u^{\prime }v^{\prime }}\rangle$ and $\langle \overline{u^{\prime \prime }v^{\prime \prime }}\rangle$ for $y/k_{t}=0,0.6$ , in the early stages of transient flow ( $t^{\ast }<0.8$ ), the values collapse for the corresponding profile at $y/k_{t}=1.8$ . In fact, an inspection of the results at all locations (not shown) indicates that, similar to the r.m.s. of velocity fluctuations, the discrepancy between the two stresses disappears for locations above $y/k_{t}\sim 1.3$ , indicating a negligible effect of dispersive stress in that region.

Figure 13. Ensemble-averaged turbulent shear stress, and sweep and ejection events during the transient flow for the $y$ – $z$ plane at $x/{\it\lambda}=0.35$ . Plots are shaded by: (a–d)  $\overline{u^{\prime }v^{\prime }}^{+0}$ ; (e–h)  $\overline{u^{\prime }v^{\prime }}_{(Q2,h4)}^{+0}$ (strong ejections); (i–l)  $\overline{u^{\prime }v^{\prime }}_{(Q4,h4)}^{+0}$ (strong sweep). Contour lines of $\overline{u^{\prime }v^{\prime }}^{+0}$ are also shown in each plot (negative values shown by dashed lines): (a,e,i)  $t^{\ast }=0.3$ , (b,f,j)  $t^{\ast }=0.4$ , (c,g,k)  $t^{\ast }=0.6$ , (d,h,l)  $t^{\ast }=4$ .

The influence of the primary vortex is further investigated in figure 12, which shows the contours of the ensemble-averaged $\overline{u^{\prime }v^{\prime }}$ over a horizontal plane at the crest. The footprint is a whole element as shown in figure 2(b). At $t^{\ast }=0.2$ , an area of high turbulent shear stress (the dark patch) is seen just behind the roughness crest, indicating the initiation of the primary vortex at this location. This strong turbulent shear stress region is then convected downstream. Consistent with the vortex and vector plots (figures 5–7), the signature of the head of a hairpin can be seen in the frame at $t^{\ast }=0.3$ , in the sense that there is a high value region of turbulent shear stress which coincides with the hairpin location detected at this instant. At $t^{\ast }=0.4$ , the large dark patch has disappeared, replaced by two long thinner regions. Note the periodic nature of the plots: the two dark patches in quadrants 2 and 3 can be seen as the continuation of those in quadrants 1 and 4. These thin regions can be associated with the legs of the counter-rotating vortex (see also figure 13). This is consistent with figure 6(b–g), which shows that the hairpin heads are ejected away from this elevation by this time, leaving behind the hairpin legs. At $t^{\ast }=0.6$ , the thinner dark patches are, similarly, associated with the remaining parts of trailing hairpin legs that are about to be ejected away from the crest. At $t^{\ast }=1.5$ , the primary vortex has left this elevation and the shear stress distribution is closely similar to that of the flow at $t^{\ast }=4$ . In addition, the distribution of the steady flow ( $t^{\ast }=27.2$ ; not shown), which is similar to that at $t^{\ast }=4$ , is consistent with the experimental results of Hong et al. (Reference Hong, Katz and Schultz2011) and Talapatra & Katz (Reference Talapatra and Katz2013), which show that the forward side of the element is a region of strong positive shear stress.

Quadrant analysis is a useful tool for studying the contribution of various flow events to the total turbulence production. Using the definition of the hyperbolic hole introduced by Lu & Willmarth (Reference Lu and Willmarth1973), the contribution of each quadrant to $\overline{u^{\prime }v^{\prime }}$ is defined as

(4.6) $$\begin{eqnarray}(\overline{u^{\prime }v^{\prime }})_{Q}=\lim _{x\rightarrow \infty }\frac{1}{T}\int _{0}^{T}(u^{\prime }v^{\prime })I(t)\,\text{d}t,\end{eqnarray}$$

where $I(t)$ is an indicator function defined so that

(4.7) $$\begin{eqnarray}\displaystyle I(t)=\left\{\begin{array}{@{}ll@{}}1,\quad & \text{if }|u^{\prime }v^{\prime }|_{Q}\geqslant hu_{rms}^{\prime }v_{rms}^{\prime },\\ 0,\quad & \text{otherwise}.\end{array}\right. & & \displaystyle\end{eqnarray}$$

The greater the value of $h$ , the stronger the events. Here, we have carried out ensemble-averaged quadrant analysis, and so the integral is carried out over a number of wavelengths in any horizontal plane rather than over time as defined in (4.6).

Figure 13 shows the development of ensemble-averaged (a–d) turbulent shear stress, $\overline{u^{\prime }v^{\prime }}$ , (e–h) strong ejection events, $\overline{u^{\prime }v^{\prime }}_{(Q2,h4)}$ , and (i–l) strong sweep events, $\overline{u^{\prime }v^{\prime }}_{(Q4,h4)}$ , for the $y$ – $z$ plane at $x/{\it\lambda}=0.35$ at several instants. The location of the plane is shown in figure 12 as dashed lines. At $t^{\ast }=4$ when the near-roughness flow has reached the steady-state condition of the final flow, strong turbulent shear stress exhibits above the forward side (red patch), along the ridge line of the roughness element. This distribution is consistent with the experimental results of Hong et al. (Reference Hong, Katz and Schultz2011) for flow structure over a pyramid roughness in the steady-flow condition. The locus of the strong turbulent shear stress for the first three instants is consistent with the high-stress area seen in figure 12, and is associated with primary counter-rotating vortices that exist at these moments (figure 5). This shows that, although vortices of various strengths are formed behind each roughness element, the primary vortex has the most striking impact on the local turbulence structure in the early stages of the transient flow. During this period, the areas of strong ejection events (green areas in figure 13 e–h) are correlated with the most intense turbulent stress regions, whereas the strong sweep events (figure 13 i–l) are not associated with the regions of high turbulent stress. According to the contour levels in these plots (ranging from 0 to 0.3), it can be seen that more than 30 % of the total positive shear stress comes from those events. That is, the dominant contribution to shear stresses in these early stages is from the violent ejection events. This correlation remains true until $t^{\ast }\sim 0.8$ (not shown), when the primary vortex is ejected from the element surface. In contrast, at $t^{\ast }=4$ and in the region close to the rough element, neither are the ejections the dominant events nor is the influence of the strong primary hairpin exhibited. This is in agreement with many studies of steady flow over a rough wall, which show that, in the region very close to the rough element, sweep events make the dominant contribution to turbulent shear stress.

Figure 14. Development of Reynolds stress anisotropy at several instants. (a–f) For a smooth wall: (a)  $t^{\ast }=0$ , (b)  $t^{\ast }=0.17$ , (c)  $t^{\ast }=4$ , (d)  $t^{\ast }=8$ , (e)  $t^{\ast }=16$ , (f)  $t^{\ast }=27.2$ . (g–l) For a rough wall: (g)  $t^{\ast }=0$ , (h)  $t^{\ast }=0.17$ , (i)  $t^{\ast }=0.4$ , (j)  $t^{\ast }=0.8$ , (k)  $t^{\ast }=1.5$ , (l)  $t^{\ast }=4$ ; arrows, from small to big, denote the values at the bottom of the channel, crest and channel centre. (m) Distribution of the invariant function $F$ at the crest at $t^{\ast }=0.4$ . The plot is shaded by value of $F$ , and contour lines of $\overline{u^{\prime }v^{\prime }}$ are also shown.

4.5. Anisotropy

Turbulence anisotropy can be examined by studying anisotropy-invariant maps (AIMs) of the Reynolds stress anisotropy tensor ( $\unicode[STIX]{x1D623}_{ij}$ ), which is defined as

(4.8) $$\begin{eqnarray}\unicode[STIX]{x1D623}_{ij}=\frac{\langle \overline{u_{i}^{\prime }u_{j}^{\prime }}\rangle }{\langle \overline{u_{l}^{\prime }u_{l}^{\prime }}\rangle }-\frac{{\it\delta}_{ij}}{3},\end{eqnarray}$$

where ${\it\delta}_{ij}$ is the Kronecker delta tensor. An AIM proposed by Lumley (Reference Lumley1978) is constructed using the second and third principal invariants, respectively defined as $-\text{II}=(\unicode[STIX]{x1D623}_{ij}\unicode[STIX]{x1D623}_{ji})/2$ and $\text{III}=(\unicode[STIX]{x1D623}_{ij}\unicode[STIX]{x1D623}_{jk}\unicode[STIX]{x1D623}_{ki})/3$ . These parameters, which are independent of coordinate choice for analysing the overall anisotropy, provide information on turbulence streaky structures and are useful in developing turbulence modelling. In addition, since these normalised parameters are not dependent on the friction velocity, they can be of interest in experimental studies where the friction velocity is not available.

Figure 14(a–l) shows the AIM profiles, at several instants, for flow over the smooth and rough walls, respectively. For the smooth wall (figure 14 a–f), soon after the commencement of the transient, the turbulence in the wall region begins shifting towards one-component turbulence (the top right vertex). This trend reflects the presence of elongated steaks in the pre-transition stage. From late pre-transition ( $t^{\ast }\sim 5.5$ ), the values of $\text{II}$ begin to reduce, and the points move along the right-hand boundary, indicating that the turbulence is becoming axisymmetric. At $t^{\ast }=16$ and later, the turbulence state is characteristic of a steady flow. The AIM profiles for the rough wall are shown in figure 14(g–l). Three arrows in each panel, from small to big in size, respectively, indicate the value at the bottom of the channel, at the roughness crest and at the centre of the channel ( ${\it\delta}_{t}$ ). At $t^{\ast }=0$ , the AIM values for the locations within the roughness height mainly follow the top right line, which characterises the two-component turbulence state, indicating a negligible contribution of the wall-normal velocity. The values exhibit the following changes during the transient flow. Soon after the transient begins ( $t^{\ast }\sim 0.17$ ), the values move towards the right boundary line, which represents a rod-like (cigar-shaped) axisymmetric turbulence state. This turbulence state, which is exhibited until the beginning of the second stage ( $t^{\ast }\sim 0.8$ ), indicates that the two turbulent components should be similar. The r.m.s. of turbulent velocity fluctuations of figure 10(g,i,k) confirm that the values of wall-normal and spanwise components are almost the same during this period. From $t^{\ast }\sim 0.8$ , the AIM values move away from the right boundary line towards the top right line. The AIM developments also show that the value of $\text{II}$ for the roughness crest reduces monotonically until $t^{\ast }\sim 1.5$ , from which point it remains almost unchanged for the rest of the transient flow.

One way of examining the overall anisotropy is via the invariant function, $F=27\text{III}+9\text{II}+1$ (Lumley Reference Lumley1978). Here $F=0$ denotes two-component turbulence, while $F=1$ indicates an isotropic state. Figure 14(m) shows the distribution of the invariant function at the roughness crest for $t^{\ast }=0.4$ . Also shown in the plot is the contour lines of the ensemble average of the wall-normal stress, $\overline{u^{\prime }v^{\prime }}$ . It is seen that the locus of the strong wall-normal stress corresponds to the lowest anisotropy. These areas are associated with the hairpin legs shown in figure 12(c).

Figure 15. Development of pre-multiplied spectra for $\langle \overline{u^{\prime }u^{\prime }}\rangle$ in the streamwise direction. (a–f) Spectra at several instants: (a)  $t^{\ast }=0$ , (b)  $t^{\ast }=0.17$ , (c)  $t^{\ast }=0.6$ , (d)  $t^{\ast }=1.5$ , (e)  $t^{\ast }=4$ , (f)  $t^{\ast }=27.2$ . (g) Time history of the first four dominant wavenumbers ( ${\it\lambda}$ , $2{\it\lambda}$ , $3{\it\lambda}$ , $5{\it\lambda}$ ) at the crest of element, $y=0$ . (h) Development of the spectrum of the fundamental wavelength ( ${\it\lambda}$ ).

4.6. Energy spectra

To examine further the influence of the roughness on the development of the turbulence structure, pre-multiplied streamwise energy spectra are presented in figure 15. The energy spectra of the streamwise normal stress, $E_{11}$ , are defined such that

(4.9) $$\begin{eqnarray}\langle \overline{u^{\prime }u^{\prime }}\rangle =\int _{0}^{\infty }E_{11}\,\text{d}k.\end{eqnarray}$$

The main influence of roughness on the streamwise energy spectra (figure 15 a–f) during the early stages of the transient flow is the presence of spikes at a wavenumber of $k_{{\it\lambda}}=2{\rm\pi}/{\it\lambda}=20.94$ (wavenumber corresponding to the roughness wavelength, ${\it\lambda}=0.3$ ) and its multiples. However, it is seen from the figure that, at the initial and final time instants ( $t^{\ast }=0$ , $t^{\ast }=27.2$ ), there is no evidence of such spikes. This suggests that, for a steady flow, the $\langle \overline{u^{\prime }u^{\prime }}\rangle$ spectra do not exhibit any strong signature of the roughness element anywhere above the crest.

Soon after the commencement of the transient ( $t^{\ast }=0.17$ ), whereas the profiles at $y/k_{t}=0.6$ and above remain almost unchanged, a big response is exhibited at the roughness crest. The structures with wavenumbers corresponding to the roughness wavelength of ${\it\lambda}=0.3$ and, to a lesser extent, its subharmonics are clearly shown. At $t^{\ast }=0.6$ , the response at ${\it\lambda}=0.3$ becomes more pronounced, exhibiting a big spike at the roughness wavelength. Meanwhile, the spectra for other wavenumbers have also responded to the transient, but with much smaller strength. These reflect the fact that the primary vortices and other flow structures, which are generated around individual roughness elements, appear periodically along the streamwise direction in the same wavelength as that of the roughness element during these early times, making a huge contribution to this turbulent quantity $\langle \overline{u^{\prime }u^{\prime }}\rangle$ . This is consistent with figure 5, in which two pairs of high- and low-speed streaks and trains of quasi-streamwise vortices appear along every element. Another thing to be noted is that no spikes are detected in the spectrum at $y/k_{t}=5.2$ . This is consistent with the observation from the vector plots (figure 7) that the direct influence of the roughness topography at the early stages is confined to a distance of $y/k_{t}\sim 3$ . It is interesting to note that, at $t^{\ast }=0.6$ and subsequent instants, the profiles at $y/k_{t}=0$ and 0.6 overlap each other. At $t^{\ast }=1.5$ , all spikes have disappeared. It is interesting to see that, at $t^{\ast }=4$ , the energy spectra of structures with $k_{x}>10$ at locations up to $y/k_{t}=5.2$ have almost reached their corresponding steady-flow values ( $t^{\ast }=27.2$ ).

To examine the development of the spectra further, time histories of the pre-multiplied streamwise energy for the first four most dominant wavenumbers at the crest level ( $y/k_{t}=0$ ) are shown in figure 15(g). It is seen that the energy spectrum at $k_{x}=k_{{\it\lambda}}$ increases after the transient flow begins, which continues until $t^{\ast }\sim 0.3$ . It then decreases half-way by $t^{\ast }\sim 0.4$ , when it begins a sharp increase, exhibiting a peak at $t^{\ast }=0.7$ . The peak value of the streamwise spectrum associated with $k_{x}=k_{{\it\lambda}}$ is approximately four times that of the second wavenumber $k_{x}=2k_{{\it\lambda}}$ (i.e. $(E_{k_{{\it\lambda}}}/E_{k_{2{\it\lambda}}})_{x,max}\approx 4$ ). At later times, the spectrum decreases monotonically, approaching the steady value at $t^{\ast }\sim 1.5$ . A close inspection of the vortex (figure 5) and the contour plots of figure 12 ( $t^{\ast }=0.3$ ), show that the initial increase, until $t^{\ast }=0.3$ , is associated with the period during which the well-defined primary hairpins are being formed. The period of $0.3<t^{\ast }<0.4$ corresponds to times when the primary hairpin is just generated within each element and its head passes up the crest of the roughness, while the associated energetic trailing legs have not yet reached this level. The sharp increase at $0.4<t^{\ast }<0.71$ corresponds to the period in which the trailing legs pass over the roughness wavelength, and move upwards towards the crest level of successive roughness. By $t^{\ast }=0.7$ , most parts of the legs of the vortical structures have moved away from the wall, and so a gradual reduction is exhibited at later times, leading to steady-flow values for this location at $t^{\ast }\sim 1.6$ .

Figure 15(h) shows the development of the energy spectra associated with the fundamental wavelength $k_{x}=k_{{\it\lambda}}$ for the locations up to $y/k_{t}=5$ . An area of high-energy spectra is seen for the locations up to $y/k_{t}\sim 1.5$ during $0.3<t^{\ast }<0.8$ . Consistent with the profiles of figure 15(a–f) and visualisations of figure 5, this indicates that the strong influence of the fundamental wavelength is restricted up to $y/k_{t}\sim 3$ and the effect lasts until $t^{\ast }\sim 0.8$ , when the major breakdown of the primary hairpin structures starts to occur.