We carry out direct numerical simulations of turbulent flow over riblets, streamwise- aligned grooves that are designed to reduce drag by modifying the near-wall flow. Twenty riblet geometries and sizes are considered, namely symmetric triangular with tip angle, and, asymmetric triangular, blade and trapezoidal. To save on computational cost, simulations are performed using the minimal-channel flow configuration. With this unprecedented breadth of high-fidelity flow data near the wall, we are able to obtain more general insights into the flow physics of riblets. As observed by García-Mayoral & Jiménez (J. Fluid Mech., vol. 678, 2011, pp. 317-347), we confirm that the drag-change curves of all the present groove geometries better collapse when reported with the viscous-scaled square root of the groove area, rather than the riblet spacing. Using a two-dimensional generalization of the Fukagata-Iwamoto-Kasagi identity in difference form we isolate the different drag-change contributions. We show that the drag increase associated with dispersive stresses carried by secondary flows can be as important as the one associated with the turbulent stresses and the pre-eminence of dispersive stresses can be estimated by the groove width at the riblet mean height.@en

Riblets reduce skin-friction drag until their viscous-scaled size becomes large enough for turbulence to approach the wall, leading to the breakdown of drag-reduction. In order to investigate inertial-flow mechanisms that are responsible for the breakdown, we employ the minimal-span channel concept for cost-efficient direct numerical simulation (DNS) of rough-wall flows (MacDonald et al. in J Fluid Mech 816: 5–42, 2017). This allows us to investigate six different riblet shapes and various viscous-scaled sizes for a total of 21 configurations. We verify that the small numerical domains capture all relevant physics by varying the box size and by comparing to reference data from full-span channel flow. Specifically, we find that, close to the wall in the spectral region occupied by drag-increasing Kelvin–Helmholtz rollers (García-Mayoral and Jiménez in J Fluid Mech 678: 317–347, 2011), the energy-difference relative to smooth-wall flow is not affected by the narrow domain, even though these structures have large spanwise extents. This allows us to evaluate the influence of the Kelvin–Helmholtz instability by comparing fluctuations of wall-normal and streamwise velocity, pressure and a passive scalar over riblets of different shapes and viscous-scaled sizes to those over a smooth wall. We observe that triangular riblets with a tip angle α= 30 ∘ and blades appear to support the instability, whereas triangular riblets with α= 60 ∘–90 ∘ and trapezoidal riblets with α= 30 ∘ show little to no evidence of Kelvin–Helmholtz rollers.@en

Riblets reduce skin-friction drag until their viscous-scaled size becomes large enough for turbulence to approach the wall, leading to the breakdown of drag-reduction. In order to investigate inertial-flow mechanisms that are responsible for the breakdown, we employ the minimal-span channel concept for cost-efficient direct numerical simulation (DNS) of rough-wall flows (MacDonald et al. in J Fluid Mech 816: 5–42, 2017). This allows us to investigate six different riblet shapes and various viscous-scaled sizes for a total of 21 configurations. We verify that the small numerical domains capture all relevant physics by varying the box size and by comparing to reference data from full-span channel flow. Specifically, we find that, close to the wall in the spectral region occupied by drag-increasing Kelvin–Helmholtz rollers (García-Mayoral and Jiménez in J Fluid Mech 678: 317–347, 2011), the energy-difference relative to smooth-wall flow is not affected by the narrow domain, even though these structures have large spanwise extents. This allows us to evaluate the influence of the Kelvin–Helmholtz instability by comparing fluctuations of wall-normal and streamwise velocity, pressure and a passive scalar over riblets of different shapes and viscous-scaled sizes to those over a smooth wall. We observe that triangular riblets with a tip angle α= 30 ∘ and blades appear to support the instability, whereas triangular riblets with α= 60 ∘–90 ∘ and trapezoidal riblets with α= 30 ∘ show little to no evidence of Kelvin–Helmholtz rollers.@en

Riblets reduce skin-friction drag until their viscous-scaled size becomes large enough for turbulence to approach the wall, leading to the breakdown of drag-reduction. In order to investigate inertial-flow mechanisms that are responsible for the breakdown, we employ the minimal-span channel concept for cost-efficient direct numerical simulation (DNS) of rough-wall flows (MacDonald et al. in J Fluid Mech 816: 5–42, 2017). This allows us to investigate six different riblet shapes and various viscous-scaled sizes for a total of 21 configurations. We verify that the small numerical domains capture all relevant physics by varying the box size and by comparing to reference data from full-span channel flow. Specifically, we find that, close to the wall in the spectral region occupied by drag-increasing Kelvin–Helmholtz rollers (García-Mayoral and Jiménez in J Fluid Mech 678: 317–347, 2011), the energy-difference relative to smooth-wall flow is not affected by the narrow domain, even though these structures have large spanwise extents. This allows us to evaluate the influence of the Kelvin–Helmholtz instability by comparing fluctuations of wall-normal and streamwise velocity, pressure and a passive scalar over riblets of different shapes and viscous-scaled sizes to those over a smooth wall. We observe that triangular riblets with a tip angle α= 30 ∘ and blades appear to support the instability, whereas triangular riblets with α= 60 ∘–90 ∘ and trapezoidal riblets with α= 30 ∘ show little to no evidence of Kelvin–Helmholtz rollers.@en

We study the flow above non-optimal riblets, specifically large drag-increasing and two-scale trapezoidal riblets. In order to reach large Reynolds numbers and large scale separation while retaining access to flow details, we employ a combination of boundary-layer hot-wire measurements and direct numerical simulation (DNS) in minimal-span channels. Although the outer Reynolds numbers differ, we observe fair agreement between experiments and DNS at matched viscous-friction-scaled riblet spacings in the overlapping physical and spectral regions, providing confidence that both data sets are valid. We find that hot-wire velocity spectra above very large riblets with are depleted of near-wall energy at scales that are (much) greater than. Large-scale energy likely bypasses the turbulence cascade and is transferred directly to secondary flows of size, which we observe to grow in strength with increasing riblet size. Furthermore, the present very large riblets reduce the von Kármán constant of the spanwise uniform mean velocity in a logarithmic layer and, thus, reduce the accuracy of the roughness-function concept, which we link to the near-wall damping of large flow structures. Half-height riblets in the groove, which we use as a model of imperfectly repeated (spanwise-varying) riblets, impede in-groove turbulence. We show how to scale the drag optimum of imperfectly repeated riblets based on representative measurements of the true geometry by solving inexpensive Poisson equations.@en