3.1. Drag reduction in terms of Stokes-flow protrusion heights and observed virtual origins

Luchini et al. (Reference Luchini, Manzo and Pozzi1991) and Luchini (Reference Luchini1996) advanced a model for predicting riblet drag reduction using viscous (Stokes) flow. The model predicts the roughness function $\Delta U^+$, which is related to the drag reduction $DR$, quantified by the fractional change in skin-friction coefficient $C_f$ relative to a baseline smooth wall, $DR\equiv 1-C_f/C_{f0}$, where subscript 0 refers to the smooth wall. The roughness function $\Delta U^+$ is given by $\Delta U^+=U_{\delta 0}^+-U_\delta ^+$ at matched ${Re}_\tau$, where $U_\delta ^+$ is the mean velocity at the edge ($z=\delta$) of the wall-bounded turbulence. Since $C_f\equiv 2/{U_\delta ^+}^2$, we can relate $DR$ and $\Delta U^+$ at matched ${\delta ^+={Re}_\tau }$ by ${DR = 1-(1-\Delta U^+/U_{\delta 0}^+)^{-2} \sim -2\Delta U^+/U_{\delta 0}^+}$ for small drag changes, or equivalently, $DR\sim -(2C_{f0})^{1/2}\Delta U^+$ (Spalart & McLean Reference Spalart and McLean2011; Garcıa-Mayoral, Gómez-de-Segura & Fairhall Reference García-Mayoral, Gómez-de-Segura and Fairhall2019). Figure 2 demonstrates the conversion between $\Delta U^+$ and $DR$ for two values of ${Re}_\tau$.

Figure 2 shows $\Delta U^+$ of conventionally shaped riblets (e.g.   ), after Bechert et al. (Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997), and of asymmetrical () and two-scaled () geometries as a function of the riblet size. The size is characterised by the square root of the riblet groove cross-sectional area, $\ell _g^+$, as defined on the right of figure 2. By using $\ell _g^+$, the optimal (maximum) drag reduction occurs at a shape-independent size of ${\ell _{g,{opt}}^+\approx 10.7\pm 1}$ (García-Mayoral & Jiménez Reference García-Mayoral and Jiménez2011a,Reference García-Mayoral and Jiménezb). For riblets below the optimal size ($\ell _g^+\lesssim 10.7$), $-\Delta U^+$ increases with $\ell _g^+$. To predict this trend, Luchini et al. (Reference Luchini, Manzo and Pozzi1991) and Luchini (Reference Luchini1996) proposed a linear viscous model for $-\Delta U^+$ designed for vanishingly small riblets ($\ell _g^+\rightarrow 0$)

(3.1)\begin{equation} -\Delta U^+\approx\mu_0(h_\parallel^+-h_\perp^+)=\mu_0(h_\parallel{/}\ell_g-h_\perp{/}\ell_g)\ell_g^+, \end{equation}

where $h_\parallel$ and $h_\perp$ are the Stokes-flow protrusion heights to locate the virtual origins perceived by streamwise and spanwise motions (velocities), respectively. The protrusion heights scaled by the riblet size, $h_\parallel /\ell _g$ and $h_\perp /\ell _g$, can be routinely obtained from Stokes-flow calculations for a given riblet shape (Luchini et al. Reference Luchini, Manzo and Pozzi1991; Bechert et al. Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997; Grüneberger & Hage Reference Grüneberger and Hage2011; García-Mayoral & Jiménez Reference García-Mayoral and Jiménez2011a), and the protrusion-height difference (${h_\parallel /\ell _g-h_\perp /\ell _g}$) expresses the effect of riblet shape on the drag performance. The proportionality constant, $\mu _0$ in (3.1), relates the protrusion heights to $-\Delta U^+$. The ansatz of Luchini (Reference Luchini1996) implies $\mu _0=1$. However, the empirical values, $\mu _0=0.785$ (Bechert et al. Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997; Grüneberger & Hage Reference Grüneberger and Hage2011) and $\mu _0=0.66$ (Jiménez Reference Jiménez1994; García-Mayoral & Jiménez Reference García-Mayoral and Jiménez2011b) have also been used. Figure 3(a) shows $\Delta U^+$ compensated by the respective ${h_\parallel /\ell _g-h_\perp /\ell _g}$ and compares these DNS results with (3.1) using the aforementioned values of $\mu _0$. Here, the 30 % scatter in $\mu _0$ from the literature is reproduced by the present DNS data for riblets sizes up to $\ell _g^+\approx 8$. This scatter suggests that $\mu _0<1$ may be based on the $\Delta U^+$ of non-vanishing riblet sizes that departs from the linear trend of (3.1). Towards vanishingly small riblet sizes ($\ell _g^+\approx 5$), discrepancies for different geometries are small and tend towards $\mu _0=1$, consistent with Luchini (Reference Luchini1996). However, we observe typical departures from linearity of (3.1) near the optimum ($\ell _g^+\approx 10.7$). Extrapolating (3.1) to the optimal size overpredicts the drag reduction by up to 40 % (figure 3a). For riblets, which yield small drag changes, such overprediction is especially undesirable. To account for this departure at the optimum ($\ell _{g,{opt}}^+$), García-Mayoral & Jiménez (Reference García-Mayoral and Jiménez2011a) introduced an empirical constant $0.83$ when evaluating $DR$ and (3.1), ${DR_{{max}}=(2C_{f0})^{1/2}0.83\mu _0(h_\parallel /\ell _g-h_\perp /\ell _g)\ell _{g,{opt}}^+}$ ($\bullet$, figure 3a). However, the data in figure 3(a) indicate that this empirical constant has a 20 % discrepancy.

Figure 3. (a) Roughness function $\Delta U^+$ normalised by the protrusion-height difference, ${h_\parallel /\ell _g-h_\perp /\ell _g}$, reproduces a scatter in $\mu _0$ from (3.1), which is consistent with the scatter in $\mu _0$ in the literature (indicated by the dashed grey lines). Black marker (${\bullet }$) in (a) is an approximation of $DR_{{max}}$ by García-Mayoral & Jiménez (Reference García-Mayoral and Jiménez2011a,Reference García-Mayoral and Jiménezb) using $\mu _0=1$. (b) Value of $\Delta U^+$ normalised by the difference of the mean and turbulence virtual origins, $\ell _U/\ell _g-\ell _T/\ell _g$ from (3.2) shows a smaller scatter than the protrusion-height normalisation shown in (a).

Using the present DNS dataset, we will next review and assess the physical idea behind the aforementioned discrepancies. The Stokes-flow protrusion heights ($h_\parallel ^+$ and $h_\perp ^+$) in (3.1) are meant to capture the observed virtual origins of the mean flow ($\ell _U^+$) and turbulence ($\ell _T^+$), respectively. In fact, figure 3(b) shows that it is the difference in the observed virtual origins, $\ell _U^+-\ell _T^+$ (measured a posteriori from DNS data), that better quantifies $-\Delta U^+$. These virtual origins are illustrated in figure 4. Here, the quasi-streamwise vortices, which represent the turbulence, perceive an apparent origin at $\ell _T^+$ below the crest ($z^+=-\ell _T^+$) that depends on the conditions at the riblet crest plane: no transpiration (figure 4b) or crest transpiration (figure 4c). We discuss the difference between these two conditions in § 3.2. The virtual origin of the mean flow is offset by $\ell _U^+$ (figure 4a), which is deeper than $\ell _T^+$ for drag-reducing riblets. As such, the turbulence is effectively pushed away from the wall, which causes the mean flow to slip near the grooves, thereby reducing drag. Luchini (Reference Luchini1996) linked these two virtual origins to the drag change $\Delta U^+$ by noting that the turbulence above riblets does not change relative to a smooth wall, apart from a shift of $\ell _T^+$. If true, the relationship between $\Delta U^+$ and $\ell _U^+$ and $\ell _T^+$ follows from the stress balance (at large ${Re}_\tau$), $\text {d} U^+/\text {d} z^+=1+\overline {u^\prime w^\prime }^+$, which constrains the riblet and smooth-wall mean gradients ($\text {d} U^+/\text {d} z^+$) to match when the turbulence profiles ($\overline {u^\prime w^\prime }$) collapse at all wall-normal locations, i.e. are smooth-wall like. Figure 5(a) shows the Reynolds shear stress profile ($-\overline {u^\prime w^\prime }^+$) against the wall-normal distance from the riblet crest, $z^+$. Here, we see a scatter in $-\overline {u^\prime w^\prime }^+$ between riblet shapes and sizes. In figure 5(b), however, the same data are now shifted by their respective $\ell _T^+$ in the wall-normal direction, which collapses with the smooth-wall profile, consistent with Luchini (Reference Luchini1996). We compute $\ell _T^+$ by the wall-normal shift that optimally collapses with the smooth-wall Reynolds stress profile at heights in the range $4\lesssim z^+\lesssim 6$. We also find that the energy distribution across scales of near-wall turbulence over riblets mimics that above a smooth wall. Figure 6 shows the premultiplied two-dimensional cospectra of Reynolds shear stress at matched height above the virtual origin of turbulence of riblets ($z^++\ell _T^+\approx 5$) compared with that of a smooth wall. Here, we observe a general trend that the near-wall turbulent structures are smooth-wall like for riblet sizes below the optimal ($\ell _g^+\lesssim 10$). The smooth-wall-like behaviour above riblets suggests that quasi-streamwise vortices are maintained without a change in their statistics above the small grooves, except that they perceive an apparent origin at $z^+=-\ell _T^+$ (Gómez-de-Segura et al. Reference Gómez-de-Segura, Fairhall, MacDonald, Chung and Garcıa-Mayoral2018a). Now, in figure 4(a), if we set the origin of the smooth-wall mean profile at $z^+=-\ell _T^+$, then the difference between the origins of the smooth and riblet mean velocity profiles is $\ell _U^+-\ell _T^+$. And since there is a unity (viscous-scaled) mean velocity gradient near the wall, the difference in mean velocity approaching the crest is also $U^+-U_0^+=\ell _U^+-\ell _T^+$, where $U_0^+$ is the mean velocity above the smooth wall. Then from the stress balance, since $\overline {u^\prime w^\prime }^+$ profiles collapse, the gradient $\text {d} U^+/\text {d} z^+$ is the same at every $z^++\ell _T^+$ above the crest, and therefore, the velocity difference near the crest propagates into the overlying flow (Luchini Reference Luchini1996). This gives

(3.2)\begin{equation} -\Delta U^+\approx\ell_U^+-\ell_T^+=(\ell_U/\ell_g-\ell_T/\ell_g)\ell_g^+, \end{equation}

which we observe to be valid for riblets below the optimal size $\ell _g^+<10.7$. Figure 5(c) shows the scatter between riblet shape and sizes in the mean velocity profiles along $z^+$. In figure 5(d), the mean velocity profiles as a function of the turbulence virtual-origin height, $z^++\ell _T^+$, now collapse when shifted downwards by their respective $\ell _U^+-\ell _T^+$, thus validating (3.2). We also illustrate (3.2) in figure 3(b) to contrast with (3.1) in figure 3(a). From figure 3(b) we observe that there is still a scatter near the optimal size for (3.2) due to an onset of departure from smooth-wall-like flows, but this scatter is smaller than that of (3.1) shown in figure 3(a). Furthermore, we observe that the linearity between $\Delta U^+$ and $\ell _g^+$ near the optimum is preserved for (3.2).

Figure 4. Small (viscous-scaled) riblets offset the apparent origins of the mean and turbulent flows downwards by (a) $\ell _U$ and (b,c) $\ell _T$, respectively, relative to the riblet crests ($z=0$). The turbulence (quasi-streamwise vortices) does not change relative to a smooth wall, apart from the $\ell _T$ offset (Luchini Reference Luchini1996). The difference between the offsets quantifies the drag reduction, $-\Delta U^+\approx \ell _U^+-\ell _T^+$ (Luchini Reference Luchini1996; Garcıa-Mayoral et al. Reference García-Mayoral, Gómez-de-Segura and Fairhall2019). (b) Luchini et al. (Reference Luchini, Manzo and Pozzi1991) suggested that the near-wall turbulence is dominated by spanwise motions, which is valid for vanishingly small riblets. (c) Presently, we include the effects of transpiration at the crest plane, which is crucial in setting $\ell _T^+$ (Gómez-de-Segura et al. Reference Gómez-de-Segura, Fairhall, MacDonald, Chung and Garcıa-Mayoral2018a; Ibrahim et al. Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021) for non-vanishing riblets.

Figure 5. (a,b) Reynolds shear stress for small riblets $({\ell _g^+\lesssim 10.7})$ with (a) the wall-normal coordinate origin at the crest, $z^+$, and (b) the wall-normal origin at the turbulence virtual origin, $z^++\ell _T^+$, which displays a smooth-wall-like behaviour. Inset shows the same profiles with linear-scaled wall-normal coordinate.(c,d) Mean velocity profile for all small riblets ($\ell _g^+\lesssim 10$) as a function of (c) the crest-origin height, $z^+$, and (d) turbulence virtual-origin height, $z^++\ell _T^+$. The velocity profiles in (d) are shifted downwards by the respective $\ell _U^+-\ell _T^+$ that collapse perfectly with the smooth-wall profile ($\ell _U^+=\ell _T^+=0$), hence $-\Delta U^+\approx \ell _U^+-\ell _T^+$ (Gómez-de-Segura et al. Reference Gómez-de-Segura, Fairhall, MacDonald, Chung and Garcıa-Mayoral2018a; Garcıa-Mayoral et al. Reference García-Mayoral, Gómez-de-Segura and Fairhall2019; Ibrahim et al. Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021).

Figure 6. Premultiplied two-dimensional cospectra of Reynolds shear stress, ${\kappa _x^+\kappa _y^+E_{uw}^+}$, in the plane, ${z^++\ell _T^+\approx 5}$. The contours are normalised by the matched Reynolds stress at the height of the plane ${\overline {u^\prime w^\prime }^+=\int _{0}^{\infty }\int _{0}^{\infty } E^+_{uw}\,\text {d}\lambda _x^+\,\text {d}\lambda _y^+}$, where $\lambda _x$ and $\lambda _y$ are the streamwise and spanwise wavelengths. Contours for smooth walls (– – –) are compared with contours for riblets (filled). All available $\ell _g^+\lesssim 8$ riblets show smooth-wall-like turbulence, whilst the onset of deviation from smooth walls is observed to occur at $\ell _g^+\approx 10$, except for the two-scale trapezoids () at $\ell _g^+\approx 7$. For large riblets ($\ell _g^+\gtrsim 20$), the co-spectra contours generally break at $\lambda _y^+\approx s^+$ (——), which indicates pinning of turbulent structures by the riblet textures. Boxes near the top delimit the region of Kelvin–Helmholtz rollers (García-Mayoral & Jiménez Reference García-Mayoral and Jiménez2011b), which are present in and (Endrikat et al. Reference Endrikat, Modesti, García-Mayoral, Hutchins and Chung2021a).

The ansatz of Luchini (Reference Luchini1996) suggests that the streamwise and spanwise protrusion heights are identical to the virtual origins of the mean and turbulence, i.e. $\ell _U^+=h_\parallel ^+$ and $\ell _T^+=h_\perp ^+$, recovering (3.1) with $\mu _0=1$. Figure 7(a) demonstrates that $\ell _U^+=h_\parallel ^+$ is corroborated by the present data by observing that the ratio $\ell _U^+/h_\parallel ^+\approx 1$, even for sizes larger than the optimal ($\ell _g^+\lesssim 15$). However, for the turbulence virtual origin, $\ell _T^+\ne h_\perp ^+$, even for riblet sizes near the optimum $5\lesssim \ell _g^+<10.7$ (figure 7b). The turbulence virtual origins $\ell _T^+$ are generally deeper than $h_\perp ^+$ ($\ell _T^+/h_\perp ^+ > 1$), but seem to asymptote to $h_\perp ^+$ towards vanishingly small sizes ($\ell _g^+\rightarrow 0$), consistent with Luchini et al. (Reference Luchini, Manzo and Pozzi1991). We also observe that $\ell _T^+$ deviates less from $h_\perp ^+$ for riblets with narrower (higher height-to-spacing ratio) grooves (e.g. ) for sizes below the optimal ($\ell _g^+\lesssim 10.7$). For the post-optimal riblet sizes ($\ell _g^+\gtrsim 10.7$), a much larger deviation from $h_\perp ^+$ is observed. Note that the value of $\ell _T^+$ for these post-optimal riblets does not accurately represent the turbulence virtual origin. When attempting to collapse the Reynolds stress profile in the range $4\lesssim z^+\lesssim 6$ (to find $\ell _T^+$), the virtual origin of the Reynolds stress is not well defined due to the departure from smooth-wall-like turbulence. For pre-optimal riblets, however, the virtual origin of the Reynolds stress profile is the same as the smooth wall when shifted by $\ell _T^+$, as observed in figure 5(b). Within these pre-optimal riblet sizes, we seek an accurate (viscous) drag model capable of predicting $\ell _T^+$ for all riblet shapes. One reason $\ell _T^+>h_\perp ^+$ is because the calculation of $h_\perp ^+$ neglects transpiration at the riblet crest plane (Gómez-de-Segura et al. Reference Gómez-de-Segura, Fairhall, MacDonald, Chung and Garcıa-Mayoral2018a). Here, transpiration refers to the spanwise-varying wall-normal flow in the cross-plane (predominantly induced by quasi-streamwise vortices), and does not refer to the variations of the streamwise velocity in the streamwise direction (cf. Bottaro Reference Bottaro2019). The slip/transpiration simulations of Habibi Khorasani et al. (Reference Habibi Khorasani, Lācis, Pasche, Rosti and Bagheri2022), which permit transpiration in both of these senses, show that it is indeed transpiration due to spanwise variation of the spanwise velocity that determines the displacement of quasi-streamwise vortices, and thus to the near-wall turbulent mixing and the generation of Reynolds stresses upon which $\ell _T^+$ is defined.

Figure 7. Comparison between (a) the observed mean virtual origin $\ell _U^+$ from the DNS, and Stokes-flow streamwise protrusion height $h_\parallel ^+$ and between (b) the observed turbulent virtual origin from the DNS $\ell _T^+$, and Stokes-flow spanwise protrusion height $h_\perp ^+$, for non-vanishing riblet sizes, $\ell _g^+\gtrsim 5$. The mean virtual origin $\ell _U^+$ is computed by the distance from the crest to the virtual origin of the linearly extrapolated mean velocity profile with a gradient measured locally at $z^+\approx 1$. The turbulence virtual origin $\ell _T^+$ is computed by the wall-normal shift that best collapses the smooth-wall Reynolds stress profile within heights in the range $4\lesssim z^+ \lesssim 6$.

3.2. Transpiration effects at the riblet crest plane

Gómez-de-Segura et al. (Reference Gómez-de-Segura, Fairhall, MacDonald, Chung and Garcıa-Mayoral2018a) and Ibrahim et al. (Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021) show that the turbulence virtual origin $\ell _T^+$ is not set by the virtual origin for the spanwise velocity alone, but that wall-normal velocity (transpiration) also plays a role. They performed textureless DNS channel simulations with a Robin boundary condition prescribed on a reference plane ($z=0$) for the streamwise (${u\rvert _{z=0}=\ell _x\partial {u}/\partial {z}\rvert _{z=0}}$), spanwise (${v\rvert _{z=0}=\ell _y\partial {v}/\partial {z}\rvert _{z=0}}$) and wall-normal (${w\rvert _{z=0}=\ell _z\partial {w}/\partial {z}\rvert _{z=0}}$) instantaneous velocities, where $\ell _x$, $\ell _y$ and $\ell _z$ are prescribed slip/transpiration lengths. These DNSs are termed slip/transpiration simulations. Physically, the streamwise and spanwise slip lengths, $\ell _x$ and $\ell _y$, correspond to the local equivalent distances below the reference plane at which a no-slip boundary condition is satisfied if the flow field is extrapolated linearly below the reference plane (Lauga & Stone Reference Lauga and Stone2003). For the wall-normal velocity, however, $\ell _z$ does not convey a slip effect but provides a local transpiration at the reference plane (Gómez-de-Segura et al. Reference Gómez-de-Segura, Fairhall, MacDonald, Chung and Garcıa-Mayoral2018a). The slip/transpiration simulations of Gómez-de-Segura et al. (Reference Gómez-de-Segura, Fairhall, MacDonald, Chung and Garcıa-Mayoral2018a) and Ibrahim et al. (Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021) demonstrate that $\ell _T^+$ does not depend on the prescribed streamwise slip length $\ell _x^+$, no matter how large, but only on both $\ell _y^+$ and $\ell _z^+$ (cross-flow slip and transpiration lengths). This dependency suggests that $\ell _T^+$ is the origin perceived by the quasi-streamwise vortices that induce velocities in $y$ and $z$ (cf. Habibi Khorasani et al. Reference Habibi Khorasani, Lācis, Pasche, Rosti and Bagheri2022), as first proposed by Luchini et al. (Reference Luchini, Manzo and Pozzi1991). The independence to $\ell _x^+$, meanwhile, suggests that the origin of the streaks (the streamwise velocity fluctuations) does not influence $\ell _T^+$ (figure 11 of Ibrahim et al. Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021).

From the results of these slip/transpiration simulations, Ibrahim et al. (Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021) proposed an empirical expression for $\ell _T^+$ in the smooth-wall-like regime

(3.3) \begin{equation} {\ell_{T,{fit}}^+\approx \begin{cases} \displaystyle \ell_w^++\frac{(\ell_v^+-\ell_w^+)}{1+(\ell_v^+-\ell_w^+)/5}, & \text{if }\ell_v^+ > \ell_w^+,\\ \ell_v^+ & \text{if }\ell_v^+\leqslant\ell_w^+, \end{cases}} \end{equation}

as a function of the virtual origins of the spanwise and wall-normal root-mean-squared (r.m.s.) velocities, $\ell _v^+$ and $\ell _w^+$, respectively. These virtual origins ($\ell _v^+$ and $\ell _w^+$) are computed by collapsing the near-wall r.m.s. profiles from the DNSs with that of the smooth wall, and are related to (but distinct from) the slip/transpiration lengths ($\ell _y^+$ and $\ell _z^+$). Here, when transpiration is allowed at the boundary-condition plane (non-zero r.m.s. wall-normal velocity, i.e. $\ell _w^+\ne 0$), the quasi-streamwise vortices are able to penetrate closer to the wall, which deepens the turbulence virtual origin $\ell _T^+$ (Gómez-de-Segura et al. Reference Gómez-de-Segura, Fairhall, MacDonald, Chung and Garcıa-Mayoral2018a), as shown by (3.3).

For pre-optimal riblets ($\ell _g^+\lesssim 10.7$), we presently use (3.3) to explain the relationship between $\ell _T^+$ and transpiration (i.e. $\ell _w^+$). Figure 8(a) compares the DNS-obtained $\ell _T^+$ and the empirical $\ell _{T,{fit}}^+$, showing good agreement, except for the two-scale trapezoidal riblets () which exhibit the onset of departure from smooth-wall-like turbulence at $\ell _g^+\approx 7$ (figure 6m). Figure 8(b) demonstrates that all the riblet geometries tested here come close to satisfying $\ell _v^+\approx \ell _w^+$ and $\ell _w^+\ne 0$, the latter showing that crest-transpiration effects are relevant. As such, $\ell _T^+\approx \ell _v^+\approx \ell _w^+$ according to (3.3). This agreement (that $\ell _T^+\approx \ell _v^+\approx \ell _w^+$) is shown graphically in figure 9, where the r.m.s. cross-flow velocity profiles ($v^{\prime +}$ and $w^{\prime +}$) of riblets collapse with those of the smooth wall after a wall-normal shift of $\ell _T^+$ (contrast figure 9a,b). In the case of non-vanishing riblets, the a priori spanwise protrusion height $h_\perp ^+$ differs from the a posteriori spanwise virtual origin $\ell _v^+$ (and from $\ell _T^+$), as uniform $v$ and zero $w$ (which are used to calculate $h_\perp ^+$) cannot adequately represent the finite length scale $v$ and $w$ fluctuations, typical of near-wall streamwise vortices near the crest.

Figure 8. (a) Turbulence virtual origin of riblets, $\ell _T^+$ from DNS agree with $\ell _{T,{fit}}^+$ for $\ell _g^+\lesssim 10.7$ with the exception of  riblets due to an onset of non-smooth-wall-like turbulence at $\ell _g^+\approx 7$ in figure 6.(b) Riblets of sizes $\ell _g^+\lesssim 10.7$ are in the $\ell _v^+\approx \ell _w^+$ regime regardless of shape, and thus, based on the empirical fit (3.3) by Ibrahim et al. (Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021), $\ell _{T,{fit}}^+\approx \ell _v^+\approx \ell _w^+$.

Figure 9. (a,b) Root-mean-squared streamwise ($u^{\prime +}$), spanwise ($v^{\prime +}$) and wall-normal ($w^{\prime +}$) velocity fluctuations and (c,d) streamwise vorticity fluctuation profiles for small riblets ($\ell _g^+\lesssim 10.7$) and smooth wall (S395) from table 2 as a function of (a,c) $z^+$ (crest origin) and (b,d) $z^++\ell _T^+$ (turbulence virtual origin). The collapse of $v^{\prime +}$ and $w^{\prime +}$ profiles in (b) further corroborates (3.3), whilst that in (d) is consistent with the profiles from the slip/transpiration simulations of Gómez-de-Segura & García-Mayoral (Reference Gómez-de-Segura and García-Mayoral2020).

Figures 9(c) and 9(d) show that the r.m.s. streamwise vorticity fluctuation intensities, $\omega _x^{\prime +}$, for riblets also collapse with the smooth wall when shifted by $\ell _T^+$ ($\ell _T^+$ having been determined from the Reynolds shear stress), suggesting a strong correspondence between the streamwise vorticity and Reynolds shear stress. This collapse is also observed in the $\omega _x^{\prime +}$ profiles from slip/transpiration simulations of Gómez-de-Segura & García-Mayoral (Reference Gómez-de-Segura and García-Mayoral2020) when $\ell _v^+\approx \ell _w^+$. For the streamwise r.m.s. profiles $u^{\prime +}$, a shift of $\ell _T^+$ do not result in the collapse with the smooth-wall profile (see different peaks at $z^+\approx 15$ and slopes at $z^+\lesssim 5$ in figure 9b); a similar behaviour was observed in the $u^{\prime +}$ profiles of Ibrahim et al. (Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021). This implies that the streamwise fluctuations do not contribute to setting the turbulence virtual origin, $\ell _T^+$.

This validation of (3.3) with riblets suggests that $\ell _T^+$ can be predicted accurately when incorporating crest-transpiration effects. Presently, we propose a model to quantify $\ell _T^+$ a priori by explicitly modelling an effect of the quasi-streamwise vortices near the wall, bypassing the empirical expression to find $\ell _T^+$. In doing so, we will be able to form a physically justified model of the non-trivial behaviour of the empirical model (3.3).

3.3. Models of quasi-streamwise vortices

Vortical structures that consist of spanwise-varying cross-flow velocities have been previously modelled as a sinusoidal wave with one or a combination of two spanwise wavelengths, $\lambda _y^+$ (Coles Reference Coles1978; Chapman & Kuhn Reference Chapman and Kuhn1986; Pollard et al. Reference Pollard, Savill, Tullis and Wang1994). Pollard et al. (Reference Pollard, Savill, Tullis and Wang1994) proposed a spanwise and time-varying sinusoidal function for the prescribed upper bound streamwise ($u^\prime$), spanwise ($v^\prime$) and wall-normal ($w^\prime$) velocity fluctuations, using two values of $\lambda _y^+$ to model the autonomous behaviour of turbulence. Using these upper boundary conditions, they numerically solved the Navier–Stokes equations in a two-dimensional height-restricted domain ($L_z^+\approx 40$) with no-slip smooth or riblet wall at the bottom. This model is similar to that of Chapman & Kuhn (Reference Chapman and Kuhn1986), where they also solved the Navier–Stokes equation in a two-dimensional height-restricted domain, but only for a smooth wall (see also Minnick Reference Minnick2022). Both models by Pollard et al. (Reference Pollard, Savill, Tullis and Wang1994) and Chapman & Kuhn (Reference Chapman and Kuhn1986) resolve the turbulence profile ($\overline {u^\prime w^\prime }$) for both smooth and riblet walls, which can be used to predict $\ell _T^+$. However, these models require a Navier–Stokes simulation (involving the nonlinear terms). Coles (Reference Coles1978) suggested a simpler viscous sublayer model for the turbulence by approximating the velocity fluctuations using a time-independent, sinusoidal function, where the amplitudes are height dependent, tuned using an equation for the mean sublayer profile by Spalding (Reference Spalding1961). These velocities are approximated using one mode with a spanwise wavelength of $\lambda _y^+=100$, following from the typical spanwise spacing between two adjacent near-wall coherent streaks (Kline et al. Reference Kline, Reynolds, Schraub and Runstadler1967; Smith & Metzler Reference Smith and Metzler1983). Using this viscous sublayer model by Coles (Reference Coles1978), a turbulence profile up to $z^+\approx 15$ can be obtained. The viscous sublayer model by Coles (Reference Coles1978) is straightforward, but does not provide an extension to riblet walls. The model from Pollard et al. (Reference Pollard, Savill, Tullis and Wang1994) is able to include the riblet wall, but an unsteady numerical simulation of the two-dimensional Navier–Stokes equations is required.

These aforementioned models motivate our present modelling approach in § 4. We briefly discuss the present approach with the aid of an instantaneous wall-normal velocity field in the cross-plane in figure 10. In this region, a single isolated quasi-streamwise vortex induces an ejection event (Robinson Reference Robinson1991), which lifts up a streak, that is spaced approximately 100 viscous units from the neighbouring streak (Kline et al. Reference Kline, Reynolds, Schraub and Runstadler1967; Smith & Metzler Reference Smith and Metzler1983). These quasi-streamwise vortices are invariant across the outer-flow geometry and Reynolds number (Schlatter et al. Reference Schlatter, Li, Örlü, Hussain and Henningson2014). Each vortex has an average diameter of $d^+\approx 20\unicode{x2013}40$ (Kim, Moin & Moser Reference Kim, Moin and Moser1987; Robinson Reference Robinson1991) and is essentially streamwise aligned with a length of 300 viscous units (Jeong et al. Reference Jeong, Hussain, Schoppa and Kim1997). The cores of these vortices are, on average, located 20 viscous units above a smooth wall (Kim et al. Reference Kim, Moin and Moser1987). This average is based on the location of the local peak in the r.m.s. streamwise vorticity profile (figure 9c,d).

Figure 10. Instantaneous wall-normal velocity field in the $yz$-plane of the smooth-wall minimal-channel case (S395). Mean flow is out of the page. Near the wall ($10\lesssim z^+\lesssim 20$), the intense wall-normal-velocity regions (dark blue and dark red) are induced by quasi-streamwise vortices (circular arrows). Grey arrows indicate the cross-flow velocity field, where ejection events (which lift up streaks) can be identified and observed to have a spacing of about 100 viscous units in the span (Kline et al. Reference Kline, Reynolds, Schraub and Runstadler1967; Smith & Metzler Reference Smith and Metzler1983). On average, the core diameter of these vortices are 30 viscous units (dashed circular region of intense streamwise vorticity) with the centre (‘$+$’ symbols) located 20 viscous units above the wall. The dashed rectangular area represents the region of interest for the viscous vortex model, where the upper boundary is exposed to an overlying vortex that induces wall-normal (transpiration) and spanwise fluctuations.

The scales of the quasi-streamwise vortices inform the region of interest for the present model. Here, we employ a two-dimensional domain with a smooth or riblet wall at the bottom, following Pollard et al. (Reference Pollard, Savill, Tullis and Wang1994), but we solve the (viscous) Stokes equations, further restricting the domain height to the near-wall region, $z^+\approx 12$ (see shaded region, figure 10). Since the $\ell _T^+$ shift is near the riblet wall, we do not need an explicit model for quasi-streamwise vortices (including their autonomous behaviour) farther away from the riblet wall. Instead, restricting the domain height to the near-wall region allows us to model the viscous response below these vortices. As such, the boundary condition at $z^+\approx 12$ of our present model is an ensemble-averaged flow induced by quasi-streamwise vortices, characterised by spanwise and wall-normal velocities varying sinusoidally in the spanwise direction and modelled using a single wavelength, $\lambda _y^+$, similar to Coles (Reference Coles1978). Jiménez, del Álamo & Flores (Reference Jiménez, del Álamo and Flores2004) noted that the spanwise wavelength of a vortex is $\lambda _y^+\approx 50$, as they observed the peak in the wall-normal velocity spectral density at $\lambda _y^+\approx 50$ and $\lambda _x^+\approx 300$, aligning closely with the signature of a quasi-streamwise vortex (Kim et al. Reference Kim, Moin and Moser1987; Robinson Reference Robinson1991; Jeong et al. Reference Jeong, Hussain, Schoppa and Kim1997; Schlatter et al. Reference Schlatter, Li, Örlü, Hussain and Henningson2014). In Appendix B, we perform a parameter calibration and sensitivity study to further inform our choice of $\lambda _y^+=50$. We further show in § 4 that a single $\lambda _y^+$ from the full turbulent flow signal is sufficient to determine $\ell _T^+$. We neglect the streamwise velocity $u$ because it is not responsible in setting $\ell _T^+$ (Ibrahim et al. Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021). Since $u$ is absent, the present model does not provide a turbulence profile $\overline {u^\prime w^\prime }$. Instead, we calculate $\ell _T^+$ from the model by collapsing the r.m.s. streamwise vorticity profiles with a smooth-wall reference case. Indeed, the r.m.s. vorticity profiles also shift with the turbulence above riblets as shown in figure 9(d).

The spirit of the model is similar to the stretched spiral vortex by Lundgren (Reference Lundgren1982), where a simple flow structure allows us to calculate properties of the turbulent flow (cf. Pullin & Saffman Reference Pullin and Saffman1993). In our present model, we allow a simple vortical structure to interact with non-vanishing riblet textures in order to yield $\ell _T^+$. The present approach departs from homogenisation techniques, which require the texture sizes to be vanishingly small compared with the turbulent scales. It also bypasses resolving the autonomous, nonlinear behaviour of turbulence over riblets (Pollard et al. Reference Pollard, Savill, Tullis and Wang1994; Minnick Reference Minnick2022), by predicting an averaged near-wall response induced by a quasi-streamwise vortex.

4.1. Formulation and methodology

This section includes a brief overview of the formulation and methodology of the viscous vortex model. We discuss further details of the model in §§ 4.2–4.3.

Figure 11 shows the viscous vortex domain used to model the flow below a quasi-streamwise vortex (representing turbulence). Here, we model a two-dimensional flow in the cross ($yz$) plane, i.e. $\partial /\partial x=0$, given that the quasi-streamwise vortices are relatively long, ${\approx }300$ (Jeong et al. Reference Jeong, Hussain, Schoppa and Kim1997; Schlatter et al. Reference Schlatter, Li, Örlü, Hussain and Henningson2014) compared with the optimal riblet spacing ($s^+\approx 15$, table 2). The model is an ensemble-averaged flow structure below a quasi-streamwise vortex (figure 10). The strength of this averaged structure is obtained by averaging the cross-flow Fourier amplitudes from smooth-wall DNSs (see § 4.2). Nonlinearity is essential in setting the strength of the near-wall structures (streaks and vortices) of the autonomous self-sustaining process in the buffer region ${10\lesssim z^+\lesssim 30}$ (e.g. Jiménez Reference Jiménez2018). However, by modelling only the flow very near the wall below the self-sustaining process in the buffer layer, we bypass directly contending with the nonlinearity. Instead, we only use a product of this nonlinear process, namely the strength of the vortical structures, as model input, which we measure once in a smooth-wall simulation since this part is identical to that of small riblets (e.g. figure 6). Thus, the modelled flow is governed by the steady Stokes equations

(4.1a–c)\begin{align} -\frac{1}{\rho}\frac{\partial p}{\partial y} + \nu\left({\frac{\partial^2 v}{\partial y^2}+\frac{\partial^2 v}{\partial z^2}}\right) = 0, \quad -\frac{1}{\rho}\frac{\partial p}{\partial z} + \nu\left({\frac{\partial^2 w}{\partial y^2}+\frac{\partial^2 w}{\partial z^2}}\right) = 0, \quad \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0, \end{align}

which are derived by neglecting nonlinearity and ensemble averaging the $y$ and $z$ components of (2.1). We neglect the $x$ component of (2.1) because $\ell _T^+$ only depends on the cross-flow velocities (Ibrahim et al. Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021). Up to here, this is the same set-up as that for calculating $h_\perp ^+$ (Luchini et al. Reference Luchini, Manzo and Pozzi1991), that is meant to be $\ell _T^+$. However, instead of prescribing spanwise-homogeneous, or infinitely large scale velocities at the top (Luchini et al. Reference Luchini, Manzo and Pozzi1991), we will instead prescribe the top velocities with a given scale (wavelength) and intensities (amplitudes)

(4.2a,b)\begin{equation} w(y,z=h) =-A\sin(\beta y), \quad v(y,z=h) =-B\sin(\beta y +\varPhi), \end{equation}

as illustrated at the top of figure 11. Here, $A$ and $B$ are the prescribed amplitudes, $\varPhi$ is the prescribed phase difference, $h$ is the prescribed height of the model domain measured from the riblet crest and $\beta$ is the period that is related to the prescribed wavelength, $\lambda _y$ (${\beta =2{\rm \pi} /\lambda _y}$). The top boundary condition (4.2) models the action below a vortex in the cross-plane, similar to the formulation of the sublayer model by Coles (Reference Coles1978). Furthermore, transpiration effects at the crest are readily incorporated through (4.2a), which is crucial to accurately predicting $\ell _T^+$ (Ibrahim et al. Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021). Using the modelled cross-flow velocities at the upper boundary (4.2) and a no-slip smooth or riblet wall at the bottom, as well as periodic boundary conditions in the span, we can solve for $v$ and $w$ using (4.1).

Figure 11. The viscous vortex model. Stokes flow (4.1) is solved in the grey-coloured two-dimensional domain with the given spanwise and wall-normal velocity boundary condition at the top ($z=h$), periodic sides and a no-slip smooth or riblet wall at the bottom. For the riblet surface, a uniform average of the solutions between $\varDelta =0$ and $\varDelta =s$ along the coordinate $\eta = y + \varDelta$ is performed to incorporate the unpinned, smooth-wall-like character of quasi-streamwise vortices. The model turbulence virtual origin, $\ell _{T,{VV}}$, is measured by the negative offset of the height of the domain $h$ measured from the crest plane such that the r.m.s. streamwise vorticity matches with that of the equivalent smooth wall at the local minimum (i.e. $z^++\ell _{T,{VV}}^+\approx 5$).

Even though the Stokes equation (4.1) is formally in the limit of zero Reynolds number, we retain the effect of finite friction Reynolds number, ${Re}_\tau$, in the Stokes calculation through, for example, the riblet spacing, $s^+$

(4.3)\begin{equation} s^+=\lambda_y^+{(s/\lambda_y)}_{{VV}}, \end{equation}

where $\lambda _y^+$ is associated with the quasi-streamwise vortex, and the ratio of riblet size to the wavelength, ${(s/\lambda _y)}_{{VV}}$, is prescribed in the Stokes calculation. Through this ratio, ${(s/\lambda _y)}_{{VV}}$, the model captures the non-vanishing riblet sizes relative to the near-wall turbulence (represented by the averaged quasi-streamwise vortex), with the latter characterised by a given universal scale, $\lambda _y^+$. From here on, sizes represented in the viscous vortex model are quantified in viscous units ($+$), using a similar calculation as per (4.3).

The prescribed parameters ($A$, $B$ and $\varPhi$) in (4.2) are obtained through a Fourier analysis of the DNS cross-flow velocities for a smooth and riblet wall, which we outline in detail in § 4.2. As the turbulence above small ($\ell _g^+\lesssim 10.7$) riblets remains smooth-wall like, we use the same $A$, $B$ and $\varPhi$ to model the flow above both smooth and riblet walls. These parameters vary with height above the wall, $z^+$, and wavelength, $\lambda _y^+$. Once a height and wavelength are chosen, these set the height and span of the model domain, respectively (figure 11). The domain height should be within $10\lesssim H^+\lesssim 15$ such that the modelled flow is below the average height of the vortex centres (Kim et al. Reference Kim, Moin and Moser1987). This height also provides enough room for the flow to recover from the riblet wall boundary condition. The domain span is restricted by positive-integer ($n$) multiples of $\lambda _y^+$, i.e. $L_y^+=n\lambda _y^+$, to ensure the spanwise periodicity of the top boundary condition (4.2).

Next, we specify the spanwise location of the riblet textures relative to the quasi-streamwise vortex, measured by $\varDelta$, illustrated at the bottom of figure 11. As we will show from the DNS in § 4.3, the quasi-streamwise vortices above smooth-wall-like riblets ($\ell _g^+\lesssim 10$) have an equal probability of residing at any spanwise location above the riblet wall, which indicates that these vortices are unpinned, like those above a smooth wall. We model this unpinned character of these vortices by a uniform average of the Stokes field across spanwise shifts of the riblet wall in the range $0\leqslant \varDelta \leqslant s$.

Using the above procedure, the modelled flow fields above smooth and riblet walls can be obtained. These flow fields are used to calculate $\ell _{T,{VV}}^+$, which is the turbulence virtual origin predicted by the viscous vortex model. Here, $\ell _{T,{VV}}^+$ is found by iterating the negative offset of the riblet domain height $h^+$ measured from the crest plane, such that the modelled r.m.s. streamwise vorticity profile, $\omega _x^{\prime }$, matches with that of the equivalent smooth wall at the local minimum (i.e. $z^++\ell _{T,{VV}}^+\approx 5$). As we will show in § 4.3, once $\omega _x^\prime$ at the local minimum is matched, the entire profile above the riblet crests ($z^+>0$) collapses with the equivalent modelled smooth wall profile. This is analogous to the collapse of the r.m.s. vorticity profiles measured from DNSs of riblet and smooth wall, as shown in figure 9(d), which reinforces that the modelled turbulence above riblets remains smooth-wall like. The output of the viscous vortex model, $\ell _{T,{VV}}^+$, for a range of riblet shape and sizes are reported and compared with $\ell _T^+$ measured from DNS in § 4.4.

4.2. Smooth-wall-like turbulence

This section outlines the process to obtain the magnitudes of $A$, $B$ and $\varPhi$ in (4.2), through a Fourier analysis of the DNS cross-flow velocities for smooth and riblet walls. We also show that the DNSs of riblets are not required for the present model because $A$, $B$ and $\varPhi$ are the same as that for smooth walls at matched height above the turbulence virtual origin.

Figure 12 illustrates the process to obtain samples of the Fourier coefficients of the wall-normal velocity from the minimal-channel smooth-wall DNS (S395). We perform the same analysis for the spanwise velocity. Here, the goal is to extract the parameters for a single spatially localised near-wall quasi-streamwise vortex using a windowed Fourier analysis. At each cross-plane (at fixed $x$ and $t$), we divide the spanwise and wall-normal ($v,w$) velocity fields into overlapping segments of spanwise length, $L_{ys}^+\approx 150$. Figure 12(a) shows a segment for $0\lesssim y^+ \lesssim 150$, in which a vortex can be identified at $y^+\approx 60$. We also take segments at spanwise increments of 50 units, $50\lesssim y^+ \lesssim 200$ and $100\lesssim y^+\lesssim 250$ (both are not shown in figure 12a) to ensure we capture a vortex at the centre of the segment. For cases where the DNS channel spans are wider (i.e. SF395 and S1000 in table 2), the fields are segmented with the same $L_{ys}^+\approx 150$ and spanwise increments, resulting in more segments per $yz$-field. Segmenting the cross-plane is crucial as these vortices are usually isolated in the cross-plane (Robinson Reference Robinson1991; Jiménez et al. Reference Jiménez, del Álamo and Flores2004; Jiménez Reference Jiménez2018). The Fourier amplitudes from unwindowed samples contain distortions resulting from the average effect of multiple vortices (figure 10). The difference between the Fourier amplitudes obtained from segmented and unsegmented signals is discussed in detail in Appendix B.2. In figure 12(b), we multiply each segmented field by the normalised Hann window function, $\varOmega (y^+)=({2/3})^{1/2}[1-\cos (2{\rm \pi} y^+/L_{{ys}}^+)]$, so that the fields are periodic in $y$, and we can take the discrete Fourier transform at each height $z^+$

(4.4)\begin{align} \left\{{\hat{v}^+,\hat{w}^+}\right\}(\lambda_{y,m}^+={2{\rm \pi}}/{\kappa^+_m},z^+) = \frac{1}{N_{y{s}}}\sum_{n=0}^{N_{y{s}}-1}\left[{\left\{{{v}^+_n,{w}^+_n}\right\}(y^+_n,z^+)\varOmega(y^+_n)}\right] \exp\left({-\text{i}y_n^+\kappa^+_m}\right)\!, \end{align}

where $\kappa _m=2{\rm \pi} m/L_{{ys}}$, $y_n=n L_{{ys}}/N_{y{s}}$ and $N_{y{s}}$ is the number of spanwise discrete points at a fixed height in a segment. In physical space, the velocities for each positive and finite wavelength $\lambda _y^+\equiv 2{\rm \pi} /\beta ^+$ can be written in terms of the Fourier coefficients as

(4.5)\begin{equation} \left\{{v^+_{\lambda_y^+},w^+_{\lambda_y^+}}\right\}({y}^+{,{z}^+}) = 2\left\{{|\hat{v}^+|,|\hat{w}^+|}\right\}\sin\left({{\beta^+}y^++ \angle{\left\{{\hat{v}^+,\hat{w}^+}\right\}} + {\rm \pi}/2}\right)\!, \end{equation}

for a fixed height $z^+$ and for each segment, where $|(\hat {\cdot})|\equiv [{\textrm {Re}\{(\hat{\cdot})\}^2 + \textrm {Im}\{(\hat{\cdot})\}^2}]^{1/2}$ and ${\angle {(\hat{\cdot})}}$ is the argument of ${{(\hat{\cdot})}}$ from the positive $\lambda _y^+$. The factor of 2 in (4.5) is included in the amplitude to combine both positive and negative $\lambda _y^+$. We illustrate (4.5) in figure 12(c) for a single $\lambda _y^+\approx 50$, extracted from the full flow field in figure 12(b). We shift the spanwise coordinate ${y}^+=\tilde {y}^+ - (\angle {\hat {w}^+}+{\rm \pi} /2)/\beta ^+$ so that these sinusoids resemble the upper boundary velocities of the viscous vortex model in (4.2)

(4.6a,b)\begin{align} w^+_{\lambda_y^+}({\tilde{y}}^+,z^+) = 2|\hat{w}^+|\sin\left({{\beta^+}\tilde{y}^+}\right), \quad v^+_{\lambda_y^+}(\tilde{y}^+,z^+) = 2|\hat{v}^+|\sin\left({{\beta^+}\tilde{y}^++ \angle{\hat{v}^+}-\angle{\hat{w}^+}}\right)\!. \end{align}

We extract each component of the signal, namely the amplitude ($a_j,b_j$) and phase difference ($\varphi _j$), given by

(4.7a,b)\begin{equation} b_j(\lambda_y^+,z^+) \equiv 2|\hat{v}^+|, \quad a_j(\lambda_y^+,z^+) \equiv 2|\hat{w}^+|, \quad \varphi_j(\lambda_y^+,z^+) \equiv \angle{\hat{v}}-\angle{\hat{w}}, \end{equation}

where $j$ is the sample index within the range $1\leqslant j\leqslant N_{{seg}}N_xN_t$. The upper bound (maximum number of samples) of $j$ is the product of the number of segments for each $yz$-field ($N_{{seg}}$), the number of $yz$-fields in the streamwise direction ($N_x$) and the number of time snapshots ($N_t$).

Figure 12. Illustration of the process to extract the Fourier coefficients from both spanwise and wall-normal instantaneous velocities in the cross-plane. Mean flow is out of the page. Instantaneous wall-normal velocity field in the $yz$-plane in figure 10 is segmented into overlapping segments of length, $L_{{ys}}^+\approx 150$: $0\lesssim y^+\lesssim 150$ (shown in a), $50\lesssim y^+\lesssim 200$ and $100\lesssim y^+\lesssim 250$ (the last two not shown in the figure). (b) Segments are multiplied by a normalised Hann window function, ${\varOmega (y^+)=({2/3})^{1/2}[1-\cos (2{\rm \pi} y^+/L_{{ys}}^+)]}$ so that a Fourier transform can be applied. (c) Resulting field for one wavelength, $\lambda _y^+\approx 50$, where samples of the Fourier amplitudes as a function of $z^+$ can be obtained.

Figures 13(a)–13(c) depict the p.d.f. of the samples (4.7) obtained from the smooth-wall DNS cases. These include the full-channel and minimal-channel smooth walls at ${Re_\tau =395}$, as well as a minimal-channel smooth-wall case at ${Re}_\tau =1000$. The samples in these p.d.f.s correspond to a height of $z^+\approx 12$, i.e. below the average vortex core height of $z^+\approx 20$ (Kim et al. Reference Kim, Moin and Moser1987) and a sinusoidal wavelength of $\lambda _y^+\approx 50$ to represent a quasi-streamwise vortex (Jiménez et al. Reference Jiménez, del Álamo and Flores2004). At this specific height and wavelength, the p.d.f.s for the smooth-wall DNS cases collapse, indicating that the variables $A$, $B$ and $\varPhi$ are independent of the channel span and friction Reynolds number ${Re}_\tau$. By averaging the amplitude and phase difference across all samples $j$, we obtain $A=0.2$, $B=0.28$ and $\varPhi =0.31{\rm \pi}$. A value $\varPhi <0.5{\rm \pi}$ indicates that, on average, the quasi-streamwise vortices are tilted with respect to the streamwise direction (cf. Jeong et al. Reference Jeong, Hussain, Schoppa and Kim1997; Schlatter et al. Reference Schlatter, Li, Örlü, Hussain and Henningson2014). Our modelling choices of $z^+=12$ and $\lambda _y^+=50$ (and corresponding values of $A$, $B$ and $\varPhi$) are additionally supported by a parameter calibration and sensitivity study outlined in Appendix B.

Figure 13. Probability density functions (p.d.f.s) of the cross-flow parameters extracted from DNS (4.7a,b). These parameters correspond to $\lambda _y^+\approx 50$ at a height of $z^++\ell _T^+\approx 12$. The p.d.f.s of the smooth-wall minimal-channel case at ${Re}_\tau =395$ are represented by the vertical bars in (a–f), compared with (a–c) the full-channel ${Re}_\tau =395$ smooth-wall case (——) and the minimal-channel ${Re}_\tau =1000$ smooth-wall case (——, grey), as well as (d–f) the blade () riblet cases with $\ell _g^+\approx 5$, 10 and 16 (dark blue to light blue). The red vertical lines represent the mean values of these sample p.d.f.s, which are prescribed in the viscous vortex model, along with $\lambda _y^+=50$ and $H^+=12$.

As demonstrated in figures 13(d)–13( f), the values of $A$, $B$ and $\varPhi$ obtained from minimal-channel DNSs of blade riblets () with sizes of $\ell _g^+\lesssim 10$ (i.e. cases BL08 and BL16 from table 2) are consistent with those from smooth-wall DNSs. For these riblets ($\ell _g^+<10$), the p.d.f.s collapse when the parameters are extracted at equivalent heights above the turbulence virtual origin, $z^++\ell _T^+\approx 12$. This collapse further shows that the turbulence above riblets behaves like that above a smooth wall, and demonstrates that only smooth-wall DNS data are required to obtain $A$, $B$ and $\varPhi$. The p.d.f.s for large, non-smooth-wall-like riblets ($\ell _g^+\approx 16$) in figures 13(d)–13( f) appear similar to the smooth-wall equivalents, but exhibit noticeably lower peaks in $b_j$ and $\varphi _j$, suggesting that $A$, $B$, and $\varPhi$ at $\lambda _y^+\approx 50$ and $z^++\ell _T^+\approx 12$ are not overly sensitive to the drag-increasing mechanisms associated with these riblets, in contrast to the full Reynolds stress cospectra at $z^++\ell _T^+\approx 5$ (figure 6e).

The above observation suggests that the information from one wavelength is sufficient to compute $\ell _T^+$. Moreover, this wavelength ($\lambda _y^+=50$) is also of the same order as the riblet sizes ($10\lesssim s^+\lesssim 15$). Thus, these parameters for obtaining $\ell _T^+$ appear robust and do not require the stricter scale separation between the riblet textures and turbulence imposed on homogenisation techniques.

4.3. The unpinned, smooth-wall-like quasi-streamwise vortices above riblets

In this section, we justify from DNS that the quasi-streamwise vortices above riblets are unpinned like those above a smooth wall. We incorporate this character of vortices in the viscous vortex model by a uniform average of the solutions across spanwise shifts $\varDelta$ of the riblet texture (figure 11).

We have previously shown, from the DNSs, that the flow above small riblets ($\ell _g^+\lesssim 10.7$) is smooth-wall like in terms of the collapse of the Reynolds stress profiles (figure 5b) and the collapse of Fourier amplitudes and phases (figures 6 and 13). Thus, we expect the quasi-streamwise vortices (representing turbulence) above riblet walls to behave as if above featureless smooth walls, which have no preferential spanwise position relative to the riblets. Figure 14 shows that these vortices above either a smooth wall or small riblets are unpinned, through a uniform p.d.f. of the DNS cross-flow velocity phase shifts at wavelengths $\lambda _y^+$ similar to the average diameter of a vortex. Here, we quantify the relative spanwise positions of the vortex by phase shifts of the spanwise and wall-normal velocities, $\phi _v\equiv \angle {\hat {v}}$ and $\phi _w\equiv \angle {\hat {w}}$, respectively. These phase shifts are found from the argument of the discrete Fourier transform of the velocity signals from (4.4). Samples of the phase shifts are taken at a height of 12 units above the turbulence virtual origin, $z^++\ell _T^+\approx 12$, from the DNSs of smooth walls (figure 14a, f) and blade () riblets (figure 14b–e,g–j). For the wavelengths considered ($10\leqslant \lambda _y^+\leqslant 60$), the uniform p.d.f.s (white regions) for the smooth and $\ell _g^+\approx 5$ blade riblets (figure 14a,b, f,g) demonstrate that the flow signals shift between $-{\rm \pi}$ and ${\rm \pi}$ at equal probability, indicative of unconstrained turbulent motions for these scales. This further demonstrates that the turbulent flow is smooth-wall like above riblets of this size. For $8\lesssim \ell _g^+\lesssim 16$ (figure 14c,d,h,i), we observe that wavelengths close to the riblet spacing, $\lambda _y^+\approx s^+$ (——) are pinned by the riblet textures, indicated by the varying p.d.f. along $-{\rm \pi}$ to ${\rm \pi}$. As these pinned wavelengths increase with riblet spacing towards the average diameters of quasi-streamwise vortices ($d^+\approx 25$, Kim et al. Reference Kim, Moin and Moser1987; Robinson Reference Robinson1991; Jeong et al. Reference Jeong, Hussain, Schoppa and Kim1997), the textures begin to interfere with these vortex scales, disrupting the smooth-wall-like behaviour of the cross-flow velocity fields. A small interference at this scale is seen for the $w$ phase shifts $\phi _w$ of $\ell _g^+\approx 10$ riblets ($\lambda _y^+\approx 2s^+\approx 30$, figure 14i), indicating that pinning may be a mechanism resulting in the onset of departure from smooth-wall-like flows at $\ell _g^+\approx 10$ (figure 6). For the riblets of size $\ell _g^+\approx 16$ (figure 14e, f), the peaks observed in the p.d.f.s of $\varphi _v$ and $\varphi _w$ for $\lambda _y^+\approx 25$ indicate that the scales of quasi-streamwise vortices are now pinned by the textures (Choi et al. Reference Choi, Moin and Kim1993; Lee & Lee Reference Lee and Lee2001). The pinning of vortices was thought to stabilise the streamwise streaks (Goldstein et al. Reference Goldstein, Handler and Sirovich1995; Goldstein & Tuan Reference Goldstein and Tuan1998), which thereby explained how smaller riblets reduced drag. However, figure 14(e,j) shows that vortex pinning is only active for the larger (post-optimal) riblets ($\ell _g^+\gtrsim 10$), suggesting that pinning is actually associated with drag increase (see also § 5.1 of García-Mayoral & Jiménez Reference García-Mayoral and Jiménez2011b).

Figure 14. Probability density functions of the phase shifts of $v$ and $w$ fluctuations ($\phi _v$ and $\phi _w$, respectively) as a function of wavelength, $\lambda _y^+$, at $z^++\ell _T^+\approx 12$. These shifts quantify the spanwise location of the flow structures of various sizes $\lambda _y^+$ in the $yz$-plane. The cases considered are (a, f) full-span smooth wall (SF395), and blade () riblets of sizes (b,g) $\ell _g^+\approx 5$, (c,h) $\ell _g^+\approx 8$, (d,i) $\ell _g^+\approx 10$ and (e,j) $\ell _g^+\approx 16$. For each $\lambda _y^+$, white regions indicate a p.d.f. of 0.5 where phase shifts are of equal probability, indicating the turbulent scales are unpinned by the wall. Regions bounded by red contours are where the p.d.f. is above 0.6, and blue contours are where the p.d.f. is less than 0.4, which shows that the turbulent scales are pinned by the riblet textures, particularly when the wavelength is approximately the size of the riblet spacing, $\lambda _y^+\approx s^+$, indicated by a solid black line (——).

From the analysis above, the quasi-streamwise vortices ($d^+\approx 25$) are unpinned above riblets with sizes below the optimal, $\ell _g^+\lesssim 10.7$. We incorporate this unpinned character of vortices in the viscous vortex model by uniformly averaging the solutions at all spanwise locations of the riblets relative to the vortex. We quantify the relative spanwise position by $\varDelta$ through a related spanwise coordinate $\eta =y + \varDelta$, which has an origin at a riblet crest (figure 11). Here, the coordinate $y$ is fixed to the location of zero wall-normal velocity at the upper boundary, i.e. $w(y=0,h)=0$ (blue line, figure 11), and $\varDelta$ is the horizontal shift of the wall relative to $y$. Hence, increasing $\varDelta$ horizontally shifts the wall whilst fixing the top velocity boundary condition. Figures 15(a)–15(c) show the solution fields for the spanwise and wall-normal velocities and the streamwise vorticity for $\varDelta =0$ (no shift, filled contours) and $\varDelta =s/2$ (half-riblet-spacing shift, solid line contours) for the two-scale trapezoidal riblets with $\ell _g^+\approx 5.5$, or $s^+=10$, where $s^+$ is computed using (4.3). Near the top, the solutions are similar for the various shifts of $\varDelta$ because of the same prescribed top boundary conditions. However, differences are seen near the wall, particularly for the vorticity (figure 15c), where we observe local peaks in vorticity near the crests. These peaks are also observed in the instantaneous flow fields of Lee & Lee (Reference Lee and Lee2001), which were previously characterised as the centres of vortices. However, the Stokes solution in figure 15(c) suggests that these peaks are instead caused by a rapid change in velocity due to the no-slip crests.

Figure 15. The Stokes solution fields of (top to bottom) spanwise velocity ($v$), wall-normal velocity ($w$) and streamwise vorticity ($\omega _x=\partial w/\partial y - \partial v/\partial z$) for the two-scale trapezoidal () riblets with $\ell _g^+\approx 5.5$ ($s^+=10$), where $s^+$ is computed using (4.3). (a–c) Two-scale riblet solutions for $\varDelta =0$ (filled contours, grey riblets) are compared with that for $\varDelta = s/2$ (——, phantom-lined riblets); (d–f) average two-scale riblet solutions for all $0\leqslant \varDelta < s$ (filled contours) show a collapse with the smooth wall (– – –) above the riblet crest plane, after a vertical offset by $\ell _{T,{VV}}^+$.

Using these two linearly independent solutions ($\varDelta =0$ and $\varDelta =s/2$), we can analytically obtain additional solutions for any other $\varDelta$ (see Appendix C for more details). We then uniformly average the solutions in the range $0\leqslant \varDelta < s$, i.e. one riblet period by

(4.8)\begin{equation} \left\{{v,w,\omega_x}\right\}(y,z) = \frac{1}{s}\int_0^{ s}\left\{{v_\varDelta,w_\varDelta,\omega_{x,\varDelta}}\right\}(\varDelta,y,z)\,\text{d} \varDelta, \end{equation}

where $({\cdot })_\varDelta$ denotes the Stokes solution for a given $\varDelta$. The phase-averaged solutions for the same two-scale trapezoidal riblets with $\ell _g^+\approx 5.5$ are illustrated in figures 15(d)–15( f) where each are compared with the equivalent smooth-wall case. Here, the value of $\ell _{T,{VV}}^+$ is determined by negatively offsetting the riblet domain height by a guessed value of $\ell _{T,{VV}}^+$, and then iterating on this value until the r.m.s. streamwise vorticity, $\omega _x^\prime$, matches with that above the equivalent smooth wall at the local minimum (i.e. $z^++\ell _{T,{VV}}^+\approx 5$). After matching $\omega _x^\prime$ at the local minimum, we also observe that the phase-averaged $v$, $w$ and $\omega _x$ collapse everywhere above the riblet crest plane ($z^+>0$, figure 15d–f), which reflects the collapse of the DNS r.m.s. profiles in figure 9(b,d). The collapse also reaffirms the smooth-wall-like flow above riblets, even in regions close to the wall, after accounting for the unpinned character of quasi-streamwise vortices.

4.4. Drag predictions for riblets

In figure 16, we report the turbulence virtual origin calculated by the viscous vortex model, $\ell _{T,{VV}}^+$, for all riblet geometries from table 2, and for blade riblets () by García-Mayoral & Jiménez (Reference García-Mayoral and Jiménez2011a,Reference García-Mayoral and Jiménezb). Figure 16(a) shows the comparison between the modelled $\ell _{T,{VV}}^+$ and the DNS $\ell _T^+$ using a ratio between the two, $\ell _T^+/\ell _{T,{VV}}^+$. Here, we observe that $\ell _{T,{VV}}^+$ differs from $\ell _T^+$ in the range of 1 %–10 % for all present riblet geometries below the optimal size. This difference is noticeably less than $\ell _T^+/h_\perp ^+$ where the error in predicting $\ell _T^+$ becomes larger (up to $40\,\%$ at the optimal) with increasing size, because $h_\perp ^+$ is applicable only for vanishingly small riblets (figure 7b). For larger (post-optimal) riblets, the increased error in $\ell _{T,{VV}}^+$ can be attributed to the departure from a smooth-wall-like turbulence as well as the pinning of quasi-streamwise vortices, both of which are not accounted for in the viscous vortex model. Figures 16(b)–16(i) show the magnitudes of $\ell _T/\ell _g$ from DNS (coloured filled symbols), $\ell _{T,{VV}}/\ell _g$ from the viscous vortex model (open white symbols and black fitted lines), and $h_\perp /\ell _g$ from the protrusion-height model (grey lines) as a function of $\ell _g^+$. For each riblet shape, we compute $\ell _{T,{VV}}/\ell _g$ for several riblet sizes, then use a best-fit line to obtain $\ell _{T,{VV}}/\ell _g$ for a continuous range of riblet sizes. We observe that the viscous vortex model also predicts $\ell _T^+ \rightarrow h_\perp ^+$ for vanishingly small riblets ($\ell _g\rightarrow 0$) (${\square }$, figure 16b), tending towards the same output as the protrusion-height model (Luchini et al. Reference Luchini, Manzo and Pozzi1991; Luchini Reference Luchini1996). Hence, $\ell _{T,{VV}}^+$ for any riblet size $\ell _g^+$ can be approximated as

(4.9)\begin{equation} \ell_{T,{VV}}^+= h_\perp^++ m_T\ell_g^{+2}, \end{equation}

where $m_T$ is the gradient of the (empirical) best-fit line of $\ell _{T,{VV}}/\ell _g$ against $\ell _g^+$ for a given riblet shape. This gradient, $m_T$, can be computed by using at least two data points for a given riblet shape (e.g. figure 16d): (i) $h_\perp /\ell _g$ from the protrusion-height model for $\ell _g^+=0$, and (ii) $\ell _{T,{VV}}/\ell _g$ for one riblet size below the optimum. In figure 16, we observe that $m_T$ decreases with increasing height-to-spacing ratio, especially obvious for the triangular riblets with systematically varied tip angle $\alpha =30^\circ, 60^\circ, 90^\circ$ (figure 16d–f). Here, a larger $m_T$ is associated with a larger $\Delta U^+$ deviation from the linear protrusion-height model (3.1) at matched $\ell _g^+$, as observed in figure 3(a).

Figure 16. (a) Ratios of the observed $\ell _T^+$ from DNS to the calculated $\ell _{T,{VV}}^+$ show that the present model accurately predicts $\ell _T^+$ for sizes below the optimum, $\ell _g^+\lesssim 10.7$. For convenience, the grey region indicates the range of $\ell _T^+/h_\perp ^+$, as in figure 7(b). (b–i) Turbulence virtual origin per unit $\ell _g$, $\ell _T/\ell _g$, as a function of riblet size, $\ell _g^+$. Coloured markers represent data from DNS, open markers are outputs of the viscous vortex model, $\ell _{T,{VV}}/\ell _g$, with a corresponding linear best-fit line (solid black). Solid grey lines indicate the spanwise protrusion height, $h_\perp /\ell _g$.

Figure 16(g) shows two calculations of $\ell _{T,{VV}}/\ell _g$ at matched sizes for the asymmetric triangle riblets (). Due to the asymmetric nature of the riblets, the Stokes solutions may differ depending on the direction of rotation of the modelled vortex. Here, we computed the additional $\ell _{T,{VV}}/\ell _g$ by reversing the direction of the modelled vortex, i.e. changing the signs of (4.2). These two $\ell _{T,{VV}}/\ell _g$ differ only by 2 %, and hence, the direction of rotation of the modelled vortex does not significantly affect $\ell _{T,{VV}}/\ell _g$ for the asymmetric riblets.

Equation (4.9) is a second-order polynomial fit to the viscous vortex model, distinct from second- or higher-order asymptotic expansions in homogenisation. Higher-order asymptotic expansions formally require a scale ratio of riblets to turbulence approaching zero for convergence. Even the smallest coherent turbulent scales ($\approx$30, such as the diameter of a near-wall vortex core) are of the order of the optimal riblet size, $s^+\approx 15$. At such a finite ratio, adding more terms need not improve the approximation of the expansion (Bender & Orszag Reference Bender and Orszag1978). As such, higher-order homogenisation techniques are formally invalid for finite viscous-scaled riblet textures (Ibrahim et al. Reference Ibrahim, Gómez-de-Segura, Chung and García-Mayoral2021). The present analysis also suggests that formal expansions, if successfully attempted, would inevitably involve tackling the nonlinear interactions of the flow that sets the model parameters ($A$, $B$, $\varPhi$ and $\lambda _y^+$), noting that even the well-known streak spacing has not been calculated by any formal mathematical procedure. In other words, unless introduced explicitly, obtaining these model parameters requires modelling the self-sustaining process. Instead, here, we use the mean information from the nonlinear flow to retain the use of linear viscous equations in the spirit of Luchini et al. (Reference Luchini, Manzo and Pozzi1991). The present approach effectively bypasses the aforementioned challenges by recognising that the nonlinearity is similar regardless of the wall conditions (i.e. smooth-wall like), and thereby, requires us to determine the flow properties only once. Correspondingly, the ability of the present approach to accurately predict $\ell _T^+$ indicates that the essential physics have been captured.

With the accurate prediction of $\ell _T^+$ at non-vanishing riblet sizes, we can now predict the drag reduction, $\Delta U^+$. As ${-\Delta U^+\approx \ell _U^+-\ell _T^+}$ (figure 5d) and $\ell _U^+\approx h_\parallel ^+$ (figure 7a), the drag prediction of the model can be written as

(4.10)\begin{equation} -\Delta U^+_{{VV}} = h_\parallel^+- \ell_{T,{VV}}^+, \end{equation}

by simply replacing $h_\perp ^+$ with $\ell _{T,{VV}}^+$ in (3.1). Similar to the protrusion-height model, (4.10) consists of a priori quantities, which are solely based on a given riblet shape and size below the optimal ($\ell _g^+\lesssim 10.7$). In figure 17, we note that $\Delta U^+_{{VV}}$ for all riblet shapes agree with $\Delta U^+$ measured from the present DNS, and from DNSs and wall-resolved LES of past studies (El-Samni et al. Reference El-Samni, Chun and Yoon2007; García-Mayoral & Jiménez Reference García-Mayoral and Jiménez2011b, Reference García-Mayoral and Jiménez2012; Bannier et al. Reference Bannier, Garnier and Sagaut2015; Li & Liu Reference Li and Liu2019; Malathi et al. Reference Malathi, Nardini, Vaid, Vadlamani and Sandberg2022; Cipelli Reference Cipelli2023), for riblet sizes up to $\ell _g^+\approx 10.7$. However, the experimental drag reduction of triangular and trapezoidal riblets from Bechert et al. (Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997), shown in figure 17(b,e, f,i), generally exhibits lower magnitudes compared with DNS/wall-resolved LES data. Experimental drag-reduction measurements have been converted to $\Delta U^+$ based on matched bulk Reynolds number following from Bechert et al. (Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997), ${{\Delta C_f}/{C_{f0}}\approx {\Delta U^+}/[{(2C_{f0})^{-1/2}+(2\kappa )^{-1}}]}$ using $\kappa =0.4$. Some error will be introduced into the converted $\Delta U^+$ values as the experiments were conducted in an asymmetric channel set-up (one smooth and one riblet wall). Nonetheless, considering figure 17(b), there is a 30 % discrepancy in the optimal performance of trapezoidal riblets () between the experimental and DNS data. A source of this discrepancy is the finite width of the milled riblet tips (e.g. ), with tip bluntness having been shown to decrease drag performance (Walsh Reference Walsh1990; Bechert et al. Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997; García-Mayoral & Jiménez Reference García-Mayoral and Jiménez2011a; Grüneberger & Hage Reference Grüneberger and Hage2011). While Bechert et al. (Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997) did not provide measurements of the width of the trapezoidal (and triangular) riblet tips, the manufacturing tolerance can be estimated based on tip measurements of scalloped riblets (figure 9 of Bechert et al. Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997), as $R/s \approx 0.01$–$0.02$ or $R \approx 0.03$–$0.11$ mm, where $R/s$ represents the ratio of effective tip radius $R$ to the riblet spacing $s$. Assuming a similar manufacturing tolerance for trapezoidal riblets implies $R/s \approx 0.03$ (), or $R=0.13$ mm based on the provided $s=4$ mm spacing (figure 24 of Bechert et al. Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997). With $R/s \approx 0.03$, the viscous vortex model would predict a degradation in optimal performance of ${\approx }15\,\%$ (figure 17b), accounting for half of the observed 30 % discrepancy. Grüneberger & Hage (Reference Grüneberger and Hage2011) repeated measurements using sharper flat-tipped trapezoidal riblets (equivalent to $R/s\approx 0.003\unicode{x2013}0.005$) with $\alpha =45^\circ$ in the same testing facility, and obtained better agreement with both the protrusion-height model and $\Delta U^+_{{VV}}$ (figure 17i), given the smaller $R/s$. A similar level of accuracy is observed when comparing the experimental data of von Deyn et al. (Reference von Deyn, Gatti and Frohnapfel2022), depicted in figure 17(j), whose riblets were manufactured to a tolerance of $R/s\approx 0.002$. Another source of disagreement with Bechert et al. (Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997) is an unmatched bulk Reynolds number. Small changes in the effective channel height have been shown to significantly change $-\Delta U^+$ (von Deyn et al. Reference von Deyn, Gatti and Frohnapfel2022). Adjusting the experimental drag curve of Bechert et al. (Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997) for the trapezoidal riblets () using an effective channel height starting from $z^+=-\ell _{T,VV}^+$ compared with that from $z^+=-0.4k^+$ (figure 24 of Bechert et al. Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997) results in a 5 % increase in $-\Delta U^+_{{opt}}$. Other sources of discrepancy include the finite fetch of the riblet test plate (cf. Garratt Reference Garratt1990; García-Mayoral & Jiménez Reference García-Mayoral and Jiménez2011a), and the unflushed transition from smooth channel to riblet test plates which induces pressure drag (cf. Bechert, Bruse & Hage Reference Bechert, Bruse and Hage2000; Spalart & McLean Reference Spalart and McLean2011; Li et al. Reference Li, de Silva, Rouhi, Baidya, Chung, Marusic and Hutchins2019; Smith, Yagle & McClure Reference Smith, Yagle and McClure2023).

Figure 17. Roughness function $\Delta U^+$ of riblets as a function of riblet groove size, $\ell _g^+$. Roughness functions measured from DNS (coloured symbols) are compared with the linear protrusion-height model of Luchini et al. (Reference Luchini, Manzo and Pozzi1991) (——, grey) and the present drag prediction (4.10) for perfectly sharp riblets (——) and a tip-rounded trapezoidal riblet with $R/s\approx 0.03$ and $k/s\approx 0.4$ (– – –) in (b). Grey markers are LES/DNS results from (b) Bannier et al. (Reference Bannier, Garnier and Sagaut2015), (e) Choi et al. (Reference Choi, Moin and Kim1993), Malathi et al. (Reference Malathi, Nardini, Vaid, Vadlamani and Sandberg2022), Cipelli (Reference Cipelli2023), ( f) Li & Liu (Reference Li and Liu2019), (h) García-Mayoral & Jiménez (Reference García-Mayoral and Jiménez2011b, Reference García-Mayoral and Jiménez2012) and (k) El-Samni et al. (Reference El-Samni, Chun and Yoon2007). White circles (${\circ }$) are experimental $\Delta U^+$ from (b,e, f,k) Bechert et al. (Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997), (i) Grüneberger & Hage (Reference Grüneberger and Hage2011) and (j) von Deyn et al. (Reference von Deyn, Gatti and Frohnapfel2022). Crosses ($\boldsymbol {\times }$) in (i) are from Bechert et al. (Reference Bechert, Bruse, Hage, Van Der Hoeven and Hoppe1997). These experimental $\Delta U^+$ data, except for von Deyn et al. (Reference von Deyn, Gatti and Frohnapfel2022), are obtained by adjusting the percentage drag reduction $DR$ based on the respective flow conditions.

As (4.10) is quadratic in $\ell _g^+$, the departure near the optimum from linearity of the extrapolated protrusion-height model (by up to 40 %, figure 3a) can also be captured. Figure 18 compares (4.10) and the extrapolated protrusion-height model (3.1) against the optimal drag reduction $\Delta U^+_{{opt}}$ obtained from the DNS/wall-resolved LES. Here, the riblet sizes are determined by the optimal drag reduction from the DNS and LES data. In figure 18(a), the overprediction of optimal drag reduction up to 40 % in the extrapolated protrusion-height model is somewhat mitigated by applying an empirical factor of 0.83 (García-Mayoral & Jiménez Reference García-Mayoral and Jiménez2011a,Reference García-Mayoral and Jiménezb), which instead leaves a discrepancy of $\pm 20\,\%$ among different riblet geometries. Note that, in the original work of Luchini et al. (Reference Luchini, Manzo and Pozzi1991), it was never suggested that the protrusion-height model could be used to predict $\Delta U^+_{{opt}}$, but was only intended to be valid for vanishingly small riblets. In figure 18(b), we observe that replacing $h_\perp ^+$ with $\ell _{T,{VV}}^+$ further improves the optimal drag prediction, and alleviates the 40 % overprediction from before.

Figure 18. Optimal riblet drag reduction obtained from DNS, $\Delta U^+_{{opt}}$, compared with (a) an extrapolation of the protrusion-height model to the optimal size, and (b) the present drag-reduction prediction in (4.10). Markers indicate different riblet shapes consistent with figure 17. Riblet sizes are $10\lesssim \ell _g^+\lesssim 11$ depending on the optimal drag reduction available in the dataset.