Pressure drag reduction via imposition of spanwise wall oscillations on a rough wall

Abstract The present study tests the efficacy of the well-known viscous drag reduction strategy of imposing spanwise wall oscillations to reduce pressure drag contributions in transitional and fully rough turbulent wall flow. This is achieved by conducting a series of direct numerical simulations of a turbulent flow over two-dimensional (spanwise-aligned) semi-cylindrical rods, placed periodically along the streamwise direction with varying streamwise spacing. Surface oscillations, imposed at fixed viscous-scaled actuation parameters optimum for smooth wall drag reduction, are found to yield substantial drag reduction ($\gtrsim$25 %) for all the rough wall cases, maintained at matched roughness Reynolds numbers. While the total drag reduction is due to a drop in both viscous and pressure drag in the case of transitionally rough flow (i.e. with large inter-rod spacing), it is associated solely with pressure drag reduction for the fully rough cases (i.e. with small inter-rod spacing), with the latter being reported for the first time. The study finds that pressure drag reduction in all cases is caused by the attenuation of the vortex shedding activity in the roughness wake, in response to wall oscillation frequencies that are of the same order as the vortex shedding frequencies. Contrary to speculations in the literature, this study confirms that the mechanism behind pressure drag reduction, achieved via imposition of spanwise oscillations, is independent of the viscous drag reduction. This mechanism is responsible for weakening of the Reynolds stresses and increase in base pressure in the roughness wake, explaining the pressure drag reduction observed by past studies, across varying roughness heights and geometries.


Introduction and motivation
Wall-bounded turbulent flows are common in a wide range of engineering applications, such as flows over ships and submarines, and in pipes.In such cases, the bounding walls generally have a 'rough' surface topography, which adds a drag penalty in the form of pressure drag.This drag is generated from the local flow separation and subsequent wake formation behind the protruding roughness elements, similar to that noted for bluff bodies (Choi et al. 2008).An increase in the size of surface roughness increases the relative contributions of the pressure drag to the total drag, but reduces the relative viscous drag contributions (Jiménez 2004;Chung et al. 2021).The latter is associated with obliteration of the drag-producing near-wall cycle, by roughness, thereby making pressure drag the dominant contributor to the total drag in the 'fully-rough' regime.These fundamental differences in smooth and rough wall draggenerating mechanisms pose a significant challenge for shipping as well as piping industries, wherein the exposed surfaces degrade from a hydraulically smooth to fully-rough regime in their operation cycle.In such cases, a drag reduction strategy that can attenuate both smooth and fully-rough wall drag-generating mechanisms is required, and this forms the focus of the present study.Here, we will refer to a fully-rough scenario as when the total drag almost entirely comes from the pressure drag.In this respect, it is important to take note of some recent findings (Marusic et al. 2021) based on high friction Reynolds number (Re τ ) experiments, which reveal inertia-dominated scales also make statistically significant contributions to the total skin-friction drag at Re τ ≳ 10 3 .Contributions from these inertial scales, however, are kept negligible in this study by limiting all our simulations to low Re τ (≲ 300).
1.1.Literature on rough wall drag reduction Numerous active and passive drag reduction strategies have been proposed in the past for turbulent wall-bounded flows (Corke & Thomas 2018).However, the majority of these strategies focused on reducing the viscous drag by targeting the near-wall selfsustaining cycle (Kim 2011;Jiménez 2018), rendering them ineffective for attenuating the fully-rough wall drag.Some studies have found large-eddy break-up devices (LEBUs) to be effective in reducing both smooth (Corke et al. 1982) and rough wall skin-friction drag (Bandyopadhyay 1986).However, using LEBUs also adds significant parasitic drag, effectively making the total drag reduction negligible.Interestingly, recent studies by Banchetti et al. (2020), Nguyen et al. (2021) and Garcia et al. (2021) found that the well-known strategy of imposing spanwise oscillations on a turbulent rough wall flow can significantly reduce the pressure drag emerging from the roughness.This is promising since the same drag reduction strategy is also known to substantially reduce the smooth wall (i.e.viscous) drag (Akhavan et al. 1993;Gatti & Quadrio 2013, 2016;Ricco et al. 2021), suggesting continued drag reduction, if implemented on a surface with varying roughness properties over time.Notably, the same strategy has also been found to yield significant drag reduction in the scenario of a transonic flow over a wing (Quadrio et al. 2022), indicating success of this strategy across a broad range of flows.
The drag of a rough wall, as well as its overlying boundary layer properties, are a function of the geometric parameters describing the surface topography (Chung et al. 2021).This has often led to consideration of simplified roughness geometries for a systematic investigation of the influence of these parameters.For instance, past studies have considered square (Leonardi et al. 2003;Lee & Sung 2007) or cylindrical (Furuya et al. 1976;Leonardi et al. 2015) rods spanning the entire width of the testing domain (known as two-dimensional or 2-D roughness), or 3-D sinusoidal (Chan et al. 2015), hemispherical (Wu et al. 2020) or cubic roughness (Coceal et al. 2007;Lee et al. 2011) distributed over the wall.The same philosophy has also been adopted by past studies understanding the efficacy of drag reduction strategies on the rough wall drag (Bandyopadhyay 1986;Banchetti et al. 2020;Nguyen et al. 2021;Garcia et al. 2021).For example, Banchetti et al. (2020) simulated a turbulent channel flow with an isolated 2-D bump on one wall, to generate a pressure drag contribution to the total drag.The drag reduction strategy considered by Banchetti et al. (2020) involved imposition of streamwise travelling waves of spanwise velocity on the bump wall (Quadrio et al. 2009), which was found to enlarge the separation bubble in the bump's wake, while also strongly stabilizing its temporal activity.A pressure drag reduction of ∼10% was noted by Banchetti et al. (2020), alongside some viscous drag reduction, for actuation parameters that previously yielded maximum drag reduction for a smooth wall channel flow.Similarly, Nguyen et al. (2021) simulated a turbulent channel flow with transverse (i.e.2-D) square bars placed on both walls, positioned at a streamwise offset of 39 times the bar height (imposed by periodicity of the computational domain).Their study involved imposing a time-oscillating pressure gradient in the spanwise direction across the entire channel cross-section, which yielded a maximum pressure drag reduction of 22% at an optimum oscillating frequency (alongside some viscous drag reduction).Their observations were consistent with Banchetti et al. (2020), who also observed pressure drag reduction due to a decrease in the pressure upstream of the rough element, along with an increase in the base pressure (i.e.downstream of the bar).Based on a control volume analysis of the change in momentum between consecutive roughness elements, Nguyen et al. (2021) hypothesized the reduction of the Reynolds shear stresses, in the wake of the roughness elements, to be a manifestation of the reduction of pressure as well as viscous drag.A similar weakening of the Reynolds stresses in the roughness wake was noted by (Banchetti et al. 2020), suggesting applicability of the same pressure drag reduction mechanism despite differences in the roughness geometry (2-D bump or square roughness).
While both of these past studies (Banchetti et al. 2020;Nguyen et al. 2021) have significantly improved our understanding, they were confined to a quasi-isolated arrangement of roughness elements, thereby limiting their applicability to transitional rough wall scenarios.Although no detailed discussion on the pressure drag reduction mechanism was presented in either of these studies, it was speculated to be a direct consequence of the viscous drag reduction.For instance, Nguyen et al. (2021) hypothesized that the pressure drag reduction "could be due to the reduced dynamic pressure of the flow encountering the bar".Similarly, Banchetti et al. (2020) mentioned "the pressure distribution is modified by changes in friction", suggesting that pressure drag reduction is dependent on the viscous drag reduction mechanism.The present study tests this hypothesis by imposing spanwise wall oscillations on a fully rough scenario (i.e.having negligible viscous drag contributions), obtained by maintaining small streamwise offsets between subsequent roughness elements (Leonardi et al. 2003).

Drawing inspiration from separation control strategies
Besides being an artefact of the viscous drag-reduced flow upstream of the roughness element, it is plausible that the pressure drag reduction could be a direct consequence of the spanwise oscillations influencing the separating shear layer or roughness wake (Choi et al. 2008;Yakeno et al. 2015).This motivates some discussion on the findings of previous studies that have investigated the flow physics and control of separated shear layers.For both streamlined bodies (for example wings at a non-zero angle of attack) and bluff-bodies subjected to separation/pressure drag, one can find several past studies discussing active (Seifert & Pack 2002;Post & Corke 2006;Yakeno et al. 2015;Cho et al. 2016;Brackston et al. 2016) and passive methods (Lin et al. 1990;Son et al. 2011) for reducing pressure drag via delay in flow separation.The control mechanism used in all these studies aims to energize the separated shear layer (or separating boundary layer) to promote flow reattachment, thereby reducing the separation bubble (i.e.wake) size, which subsequently reduces the pressure drag.Based on the observations of Banchetti et al. (2020) and Nguyen et al. (2021), however, this is not the case when spanwise oscillations are imposed on the flow over 2-D roughness elements (they found separation bubbles to enlarge).Energizing the boundary layer via momentum injection also enhances the shear layer activity (and consequently the turbulent stresses) in the separated flow, which is also opposite to what is observed in case of rough wall drag reduction via spanwise oscillations (Banchetti et al. 2020;Nguyen et al. 2021).
Another way to reduce pressure drag, besides reducing the separation bubble size, is by weakening/attenuating the vortex shedding activity in the bluff body wake (Choi et al. 2008), which increases its base pressure and consequently reduces pressure drag.Interestingly, this control mechanism also reduces the Reynolds stresses in the bluff body wake (Desai et al. 2020;Chopra & Mittal 2022), which matches with observations made by Banchetti et al. (2020) and Nguyen et al. (2021) in the roughness wake, for drag-reduced cases.Several past studies have demonstrated weakening of the shedding activity by enforcing a mismatch in the phase of the separating 2-D shear layer from the bluff body, along its spanwise direction.For example, Bearman & Owen (1998) demonstrated suppression of vortex shedding activity by introduction of spanwise waviness on the front stagnation face of a 2-D rectangular bar.The same concept was later successfully extended to a circular cylinder by Owen et al. (2001), by changing the straight axis of the cylinder to a sinusoidal one.Darekar & Sherwin (2001) explained that such a 3-D forcing, distorts the quasi 2-D separated shear layers, making them less susceptible to roll up into a Karman vortex street, thereby suppressing the vortex shedding activity.Their explanation was reinforced later by Hwang et al. (2013) who performed linear stability analysis on two-dimensional wakes, thereby revealing that the phenomena is not dependent on the 'geometry' of the bluff body.The same suppression has also been achieved in a cylinder wake via active control techniques, such as distributed forcing (Kim & Choi 2005), which imposes blowing and suction on the top and bottom surfaces of the cylinder, with sinusoidally varying intensity along the span.Considering this background, one would expect imposition of time-varying spanwise oscillations, on the 2-D roughness elements immersed in a wall-bounded flow (which essentially act as bluff bodies), to also disrupt the two-dimensionality of the separating shear layer in the roughness wake (for any roughness geometry).

Present contributions
The present study draws inspiration from the previous works of Leonardi et al. (2003Leonardi et al. ( , 2015) ) and Nguyen et al. (2021) and reports direct numerical simulations of a turbulent flow in an open channel configuration, with 2-D semi-cylindrical rods fixed on the bottom wall in a streamwise periodic manner.Individual simulations are conducted for different streamwise offsets between rods, thereby enabling investigation of the drag reduction mechanisms for scenarios of high pressure drag domination (small offsets), as well as those with near equal pressure and viscous drag contributions (large offsets).For convenience, we will henceforth refer to the arrangement of semicylindrical elements as a 'rough' wall, but we note that generalizing the present observations to practical roughness geometries (for example, 3-D, sandpaper, etc.) would require further studies.Here, we impose time-periodic wall oscillations in the spanwise direction to achieve drag reduction (Akhavan et al. 1993), which is tested here for the first time on fully-rough cases.Discussion in §1.2 suggests that weakening Focus on Fluids articles must not exceed this page length of vortex shedding activity could be the plausible mechanism behind pressure drag reduction noted previously for transitional rough wall cases (Banchetti et al. 2020;Nguyen et al. 2021).This hypothesis will be tested rigorously in the present study.

Flow configuration and simulation strategy
As shown in figure 1a, we consider an open-channel flow of an incompressible fluid of density ρ and kinematic viscosity ν past a rough bottom wall.The latter is composed of a series of plane, smooth wall regions parallel to the x-y plane, in between semi-cylindrical rough elements of radius k, that extend across the span and are spaced at a centre-to-centre distance Λ in the streamwise direction.The flow depth H is defined as the wall-normal distance between the bottom plane wall and the upper free-shear boundary, with roughness height fixed at k = 0.1H for all but one rough wall scenarios (k = 0.2H is considered for one case; see table 1).Throughout this manuscript, we use x, y and z as the streamwise, spanwise and wall-normal directions respectively, with ( U , V , W ) representing instantaneous velocities in these directions.Here, we chose the semi-cylindrical roughness geometry, over the commonly researched square bar geometry (Bandyopadhyay 1986;Leonardi et al. 2003;Lee & Sung 2007;Nguyen et al. 2021), to avoid a geometrically imposed flow separation location on the individual roughness elements.While the presence of sharp corners, on the upstream and downstream faces of the square bars, inherently 'locks' the flow separation location, the same is not observed when a rounded bar geometry is considered (Deshpande et al. 2017;Alam 2022).In this way, the choice of a semi-cylindrical roughness element permits an unambiguous investigation of the effect of spanwise wall oscillations on the roughness wake (i.e.pressure drag).Two sets of direct numerical simulations of the above configuration were performed, and are referred to as the static and actuated cases.In the static cases, the bottom rough wall is kept stationary, while in the actuated cases the entire bottom wall (including the rough elements) oscillates in the spanwise direction at a prescribed oscillation amplitude (A) and frequency (1/T osc ), expressed as V w (t) = −A sin (2πt/T osc ), where t is time.The no-slip boundary condition is imposed at the bottom wall, whereas the standard free-slip boundary condition mimics the free- surface at the top boundary.All the simulations are performed in a computational domain that is periodic in the streamwise and spanwise directions and has a size of L x × L y × L z = 6H × 3H × H, consistent with the domain size considered by previous studies investigating flow over 2-D roughness (Miyake et al. 2001;Nagano et al. 2004;Leonardi et al. 2015).Notably, Lozano-Durán & Jiménez (2014) found this box size to be adequate for computing accurate one-point statistics, even in case of a smooth wall channel flow.The present simulations are carried out using the open-source spectral element solver NEKTAR++ (Cantwell et al. 2015).To exploit the spanwise homogeneity of the configuration, the domain is discretized by twodimensional N elem spectral elements in the x − z plane, combined with Fourier series expansion in the homogeneous y direction along 192 grid points.Within each spectral element, velocity and pressure fields are represented by Lagrange polynomials of degree 7 through the Gauss-Lobatto-Legendre points, giving viscous-scaled nominal grid resolutions, ∆x + ≈ 1.4-5.0 and ∆z + ≈ 0.6-7.3.The spanwise grid resolution is ∆y + ≈ 4.7.More details of the adopted numerical parameters, and evaluation of the adequacy of the grid resolution and domain size, are provided in appendix 1.
The simulation campaign was as follows: initially, eight static cases at varying 6 ≲ Λ/k ≲ 60 were performed at a target friction Reynolds number, Re τ s = U τ s H/ν ≈ 300.Seven of these eight cases had roughness height k = 0.1H while another had k = 0.2H (for which Λ = 10k).Here, U τ s is the friction velocity calculated based on the momentum balance of the static cases, with U where C D s is the mean total drag coefficient for the static case).Detailed information on calculating C D s is given later in this section.As a consequence of varying Λ, and in order to ensure a matched Re τ s , the value of the bulk Reynolds number, Re b = U b H/ν varies amongst the different static cases (here, U b ≡ q f /H is the bulk velocity and q f is the flow rate).However, this was required to fix the roughness Reynolds number, k + = kU τ s /ν = 30 for all static rough wall cases corresponding to k = 0.1H, since it can dictate the characteristics of the separating shear layer/wake of the roughness element.It is worth noting here that the present roughness height is a significant fraction of the characteristic outer-scale (H), owing to which outer-layer similarity of the velocity statistics can't be expected (Jiménez 2004).However, our k value is smaller than or equal to the roughness height in past studies considering 2-D roughness elements (Leonardi et al. 2003(Leonardi et al. , 2015;;Lee & Sung 2007;Nguyen et al. 2021), and has negligible blockage effect on the bulk flow characteristics.
Next, starting from an established flow field of the static cases, their actuated counterparts were simulated such that the same flow rate q f (and Re b ) was maintained for the static and actuated cases at matched Λ/k.Here, for all (except one) actuated cases, we consider the optimum wall actuation parameters noted by Gatti & Quadrio (2016) for a smooth wall channel flow, i.e.A + = A/U τ s = 12 and By imposing oscillations at a fixed T + osc and A + , on rough wall cases having constant roughness Reynolds number (k + = 30), we can unambiguously investigate effect of spanwise wall oscillations on the roughness wake and pressure drag at various Λ/k (i.e. from transitional to fully-rough scenarios).Majority analysis in this manuscript will be based on k = 0.1H (i.e., k + = 30) cases, with one case considered for k = 0.2H to demonstrate the generalizability of present conclusions.Apart from these rough wall simulations, reference static and actuated cases with a smooth bottom wall were also run at matched Re τ s , domain size and oscillation parameter values.All these simulations have been summarized in table 1.
Considering the imposed streamwise periodicity from roughness elements arranged on the bottom wall, as well as the very low Re τ s (which leads to negligible largescale/inertial contributions), the streamwise extent of the present computational domain (L x = 6H) does not act as a limitation for the majority of the rough wall cases (Λ ≤ 20k).In the case of Λ/k = 60, however, our numerical domain only had a single roughness element, owing to which flow statistics may not be completely de-correlated from the effects of the computational domain.Further, based on the previous investigations of Lozano-Durán & Jiménez (2014) and Gatti & Quadrio (2016), the present computational domain size may not be sufficient to accurately estimate two-point or higher-ordered statistics for the smooth wall case.However, we consider results from these two cases (smooth wall and Λ = 60k) solely for completeness and so their uncertainties don't influence the conclusions of our study.

Definition of flow statistics
For the purpose of analyzing flow properties in this manuscript, we define a spanwise and time-averaging operator ⟨φ⟩ of any instantaneous flow quantity φ, that exploits the inherent periodicity of length Λ along x (figure 1a) following: where T obs is the simulation steady-state interval and N cyl = L x /Λ is the number of rough elements or streamwise periodic units.In (2.1), x spans from the centre of one rough element (x = 0) to that of the next downstream (x = Λ), which we refer to here as the streamwise periodic unit of the computational domain (0 ≤ x/Λ ≤ 1). Figure 1(a) highlights this periodic unit in light blue shading for reference.
In the presence of the rough elements, the x-component of the instantaneous hydrodynamic force acting upon the bottom wall (i.e. the total drag, D), comprises both viscous ( V ) and pressure ( P ) force contributions.The streamwise variation of the time-and span-averaged total drag force (per unit streamwise width) is given as: where n is the unit normal and p is the hydrostatic pressure (i.e.different from pressure force, P ).We can further integrate (2.2) in the streamwise direction to obtain the total hydrodynamic force (i.e. the total drag, D) acting on the rough wall as: where dS is the differential surface length along the bottom boundary.Note that the pressure drag force, ⟨P ⟩(x) = 0 over the plane sections in between the semicylindrical rough elements (figure 1a), while V is the sum of viscous force contributions from both the plane and semi-cylindrical sections.Following (2.3), the total drag coefficient, are respectively the coefficients of pressure and viscous force contributions to C D .Similarly, the streamwise-resolved viscous force coefficient, To differentiate the force coefficients associated with the static and actuated cases, subscripts 's' and 'a' are used respectively, which are tabulated in table 1 for all the simulation cases.Here, with Re b kept constant, the friction velocities are different for the static (U τ s ) and actuated (U τ a ) cases, which are determined directly from the mean streamwise momentum balance by following: This yields the definition of the total drag reduction, which is: DR% reported in table 1.For example, DR% for the present smooth wall case is 33.9% which, despite the marginal computational box size (Gatti & Quadrio 2016), lies within the predictions made by Gatti & Quadrio (2013, 2016) at this Re τ s .

Results and discussions
We begin by comparing the present static rough wall statistics with previously published data (Leonardi et al. 2015).For this purpose, results presented in figures 2(a), 3(a,c) and 4(a,b) will be discussed first, before proceeding towards investigating the change in statistics owing to drag reduction.Figure 2(a) shows fractions of the total drag coefficient (C D s ), coming from the pressure drag (C p s ) and viscous drag (C v s ) contributions, for cases corresponding to both k = 0.1H and 0.2H.For cases of Λ/k ≲ 10, C p s is essentially the sole contributor to C D s , with relative contributions from C v s gradually increasing (and those from C p s decreasing) with increasing Λ/k ≳ 10.At Λ/k = 60, notably, both contribute nearly equally to C D s .Further, for Λ/k = 10, the pressure and viscous contributions to the total drag are found to be similar irrespective of k = 0.1H or 0.2H.This nature of variation of C p s and C v s , with Λ/k, is consistent with Leonardi et al. (2015), suggesting the mean flow physics observed by the latter can be extended to the present study.This is despite Leonardi et al. (2015) employing a circular rod geometry for the rough elements, and conducting simulations for different Λ/k at a constant bulk flow Reynolds number, both of which are different from the approach adopted in the present study.Leonardi  2015) also noted that the variation of pressure and viscous drag with Λ/k did not change significantly on increasing the flow Reynolds number (at least in the low Reynolds number regime).This establishes that a similar variation in pressure drag domination can be expected on changing the streamwise offsets between 2-D rough elements (for 0.1H ≤ k ≤ 0.2H), irrespective of differences in roughness geometry, flow Reynolds number and simulation approach (i.e.constant friction versus bulk Reynolds numbers).Figure 2(a), hence, demonstrates that the present rough wall cases are valid for testing the efficacy of spanwise wall oscillations on varying degrees of pressure drag domination.Based on our definition given in §1, we henceforth refer to the Λ/k ≲ 10 cases as fully-rough, while those corresponding to 10 < Λ/k ≤ 60 referred to as transitionally rough.
The flow physics associated with the variation of C v s and C p s can be understood based on figures 3(a,c) and 4(a,b).The former depicts the normalized Reynolds shear stresses (⟨u ′ w ′ ⟩/U b 2 ) and mean flow streamlines, while the latter shows the timeand span-averaged viscous force coefficients (⟨C v ⟩) and normalized surface-pressure (2⟨p⟩/ρU 2 b ), plotted across the streamwise periodic unit (0 ≤ x/Λ ≤ 1) for Λ = 10k and 20k with k = 0.1H.These two cases were selected as representatives of the fully and transitionally rough scenarios, respectively.Here, u ′ and w ′ are respectively the streamwise and wall-normal velocity fluctuations obtained via u ′ = U − ⟨U ⟩.It is evident from figure 4(a) that ⟨C v ⟩ ≲ 0 across the plane region between the rough elements (1 ≤ x/k ≤ 9), for Λ = 10k.On the other hand, for Λ = 20k, ⟨C v ⟩ < 0 only immediately downstream of the rough element (1 ≤ x/k ≤ 9), after which ⟨C v ⟩ > 0 for 7 ≤ x/k ≤ 19 (figure 4b).These trends are consistent with observations from the mean flow field in figure 3(a), wherein the wake/re-circulation region can be noted across the majority of the plane region between subsequent rough elements, in case of small streamwise offsets (Λ ≲ 10k).Qualitatively similar mean flow fields and surface statistics (⟨C v ⟩, ⟨p⟩) were noted in case of k = 0.1H and 0.2H for Λ = 10k (not shown here for brevity).For larger Λ/k (figure 3c), however, the separated shear layer from upstream roughness reattaches in the plane region between these elements (x ∼ 7k), leading to an extended region of wall-attached flow impinging onto the downstream roughness element (resulting in ⟨C v ⟩ > 0).While the variation of ⟨C v ⟩(x) in the plane region is opposite for the transitional and fully rough cases (figure 4), the nature of variation of the surface-pressure is very similar, with ⟨p⟩(x) much higher on the upstream (relative to the downstream) section of the rough element, which generates the pressure drag.Notably, both of the present results plotted in figures 3(a,c) and 4(a,b) are consistent with the observations made previously by Leonardi et al. (2015), explaining the similar trends in pressure and viscous drag noted in figure 2(a).
With the static cases demonstrated to be consistent with the literature, we next investigate the drag reduction achieved on imposing spanwise wall oscillations for the same rough wall cases.Figure 2(b) plots DR%, on imposition of wall-oscillations at A + = 12 and T + osc = 100, for various Λ/k.The plot suggests substantial DR% is achieved for all wall-oscillation cases at T + osc = 100 corresponding to k = 0.1H.The DR% associated with k = 0.2H, however, is not as high as noted for k = 0.1H at Λ = 10k, and this difference is likely owing to the same spanwise forcing magnitude (A + = 12) imposed on flows with significantly different pressure drag magnitudes (nearly double for k = 0.2H than 0.1H; see table 1).The most interesting observation is the maximum DR (≈ 34%) noted at Λ/k = 7.5, which is at par to that noted via pure viscous drag reduction for a smooth wall (for the same A + and T + osc ).This is a novel result and is particularly notable considering our findings in 2(a), which indicate negligible viscous drag contributions at Λ/k = 7.5.It suggests this DR% is solely from pressure drag reduction, which has to be driven by an independent mechanism.To quantify what fraction of the total DR% is achieved by decreasing the pressure and viscous drag, figure 2  contribute to the total DR% significantly.This can be understood by comparing the mean flow field (figures 3c,d), as well as ⟨C v ⟩(x) and ⟨p⟩(x) between static and actuated cases (figure 4b), for Λ = 20k.Consistent with Banchetti et al. (2020), the pressure drag reduction can be associated with extension of the wake size downstream of the roughness elements (refer to streamlines in figures 3c,d), while viscous drag reduction can be associated with the weakening of the near-wall turbulence in the wall attached flow between the elements (refer to Reynolds shear stresses at z < 0.1H).This is supported quantitatively by figure 4(b), where ⟨C v ⟩ immediately upstream of the rough element (7 ≲ x/k ≲ 19) reduces, compared to the static case.This is accompanied by reduction in ⟨C v ⟩ in the upstream section of the semi-cylindrical surface, combination of which leads to a considerable viscous drag reduction noted in figure 2(b).Similarly, 2⟨p⟩/ρU 2 b can be noted to decrease in the upstream section of the rough element (7 ≲ x/k ≲ 19.5; figure 4b), while it increases downstream of the roughness (0 ≲ x/k ≲ 4).This explains the drop in pressure drag in figure 2(b), which Nguyen et al. (2021) had associated with the decrease in the dynamic pressure in 7 ≲ x/k ≲ 19.Consistent with the literature (Banchetti et al. 2020;Nguyen et al. 2021), the drag reduced flow for Λ = 20k also corresponds with weakening of the Reynolds shear stresses in the roughness wake (z ≈ 0.1H; figures 3c,d).This reaffirms the weakening of the Reynolds shear stresses on pressure drag reduction, irrespective of differences in the roughness geometries (semi-cylindrical, cylindrical or square).
As Λ/k decreases along 60 ≥ Λ/k ≥ 10, the streamwise extent of the wall attached flow between rough elements reduces (figures 3a,c), which reduces the contribution of viscous drag reduction to the total DR% (figure 2b).For 5 ≲ Λ/k ≲ 10, owing to the absence of flow reattachment/wall-attached flow between rough elements (i.e.⟨C v ⟩(x) ≲ 0; Leonardi et al. 2015), almost all the drag reduction is associated with a decrease in the pressure drag.This is observed irrespective of the variation in roughness height, i.e. for both k = 0.1H and 0.2H.The dominant pressure drag reduction can be confirmed by comparing ⟨C v ⟩ and ⟨p⟩ plotted for the static and actuated cases (T + osc = 100) at Λ = 10k and k = 0.1H (figure 4a), where one can note the drop in ⟨p⟩ upstream of the roughness (5 ≲ x/k ≲ 9.5), and increase in ⟨p⟩ downstream of the roughness (0 ≲ x/k ≲ 4), while there is negligible variation in ⟨C v ⟩ across the plane region (1 ≲ x/k ≲ 9).The same can be understood qualitatively by comparing the mean flow fields of the static and actuated cases depicted for Λ = 10k, in figures 3(a,b).In the static scenario, a very small region of wall-attached mean flow can be noted between the two roughness elements, which however completely vanishes on imposition of wall oscillations.The fact that significant pressure drag reduction is noted without any changes in the viscous drag (for 6 ≲ Λ ≲ 10k) suggests existence of a unique pressure drag reduction mechanism.This observation is valid for both k = 0.1H and 0.2H (figure 2b), and confirms that the same pressure drag reduction mechanism would be applicable for varying roughness heights.These results are novel, and confirms that pressure drag reduction in a 2-D rough wall is not dependent on the viscous drag reduction.Another new observation from figures 3(a,b) is that the Reynolds shear stresses are attenuated in the wake of the roughness elements, even when viscous drag reduction is negligible.It suggests that this attenuation of the turbulent stresses, which are noted across 6 ≤ Λ/k ≤ 60 (not shown), is likely associated with the pressure drag reduction mechanism.Despite differences in their roughness geometries compared to the present study, both Banchetti et al. (2020) and Nguyen et al. (2021) noted a similar reduction in the Reynolds stresses and an increase in base pressure in their roughness wake (on imposition of oscillations).These consistent results suggest that the same pressure drag reduction mechanism is applicable for different roughness geometries, which, however can only be confirmed after conducting dedicated simulations (which is beyond the present scope).

Mechanism behind pressure drag reduction
Referring to our previous discussion in §1.2, regarding bluff body drag reduction based on increases in base pressure (i.e.⟨p⟩ at x/k ∼ 1) accompanied by weakening of Reynolds stresses, our observations from figures 3 and 4(a,b) point towards weakening of the vortex shedding in roughness wake as the plausible pressure drag reduction mechanism (Choi et al. 2008).This can be investigated by analyzing the surfacepressure fluctuations (p ′ ) on the rough element; past studies on the flow past a circular cylinder (Deshpande et al. 2017;Desai et al. 2020;Chopra & Mittal 2022) have shown that weakening of its shedding activity is reflected by a drop in the lift force fluctuations, as well as surface-pressure fluctuations at the crest of the cylinder (i.e. at θ = π/2, where θ is the azimuthal angle from the front stagnation point).This sensitivity is owing to its proximity to the azimuthal location, from where the shear layer separates and sheds into the wake (0.4π ≤ θ ≤ 0.65π; Desai et al. 2020;Chopra & Mittal 2022).With this background, we consider p ′ (x, y, z, t) = p(x, y, z, t)−⟨p⟩(x) over the surface of the cylinder, where tilde ( ) and prime ( ′ ) respectively denote instantaneous and fluctuating properties.Figure 4(c) plots the normalized rootmean-square of the surface-pressure fluctuations (2 ⟨p ′ 2 ⟩/ρU 2 b ) across the surface of the semi-cylindrical rough element, with 0 ≤ θ ≤ π/2 in blue and π/2 ≤ θ ≤ π in yellow shading.It is evident that pressure fluctuations drop at θ = π/2, on imposition of wall oscillations at T + osc = 100, for both Λ = 10k and 20k and different roughness heights (k = 0.1H and 0.2H).Notably, the drop in pressure fluctuations for k = 0.2H is relatively smaller than that for k = 0.1H at θ = π/2, and this is consistent with the relatively lower DR% noted in case of the former (figure 2b).
The fact that wall-actuations weaken the vortex shedding suggests that the actuation frequency (1/T + osc ) must have a connection with the shedding frequency.Hence, we test the sensitivity of pressure drag reduction to T + osc , by imposing actuations differing by an order of magnitude.Interestingly, T + osc = 10 yields negligible pressure and total drag reduction for Λ = 10k (table 1).Indeed, when ⟨p⟩, ⟨C v ⟩ and ⟨p ′ 2 ⟩ profiles are plotted for the T + osc = 10 actuation case in figure 4(a), the corresponding profiles overlap with those of the static case, suggesting negligible changes in the mean flow physics.We can thus conclude that pressure drag reduction, on imposition of spanwise wall oscillations, is associated with weakening of the vortex shedding activity, and this explains the sensitivity of DR% to T + osc (also see Nguyen et al. 2021).Further information for interested readers, this weakening of the vortex shedding activity has also been confirmed by comparing the mean spanwise vorticity in the roughness wake, between static and actuated cases, but is not shown here for brevity.
In order to confirm the causal relationship that weakening of shedding activity leads to pressure drag reduction, we investigate the spectral distribution of p ′ at θ = π/2, and of the fluctuating pressure-force coefficient (C ′ p ).Here, where dS is the differential streamwise distance across the entire bottom wall.We analyze their energy spectra in the time/frequency domain to relate the spectral distribution to the imposed walloscillation frequency (1/T + osc ).Figures 5(a,b) respectively plot the two spectra in the premultiplied form (i.e.f ϕ p ′ p ′ and f ϕ C ′ p C ′ p ), for all the static rough wall cases (6 ≤ Λ/k ≤ 60) corresponding to k = 0.1H.They are plotted against the frequency scale (f ), which is normalized in viscous units (f ν/U 2 τ s ).Here, considering that the semi-cylindrical rough element is exposed to an incoming turbulent shear flow, we would expect the spectra to be much more broadband than that noted for a cylinder exposed to a free-stream (Darekar & Sherwin 2001).Indeed, on looking at the spectra for all the static cases, we found that both ϕ p ′ p ′ and ϕ C ′ p C ′ p are most energetic across a broadband frequency range, O(10 −3 ) ≲ f ν/U 2 τ s ≲ O(10 −2 ) (shaded in magenta), which is also exhibited by their cross-correlation spectra, ϕ C ′ p p ′ (not shown here to avoid clutter).This close correspondence in spectral energy distribution suggests that the pressure drag (which is represented by ϕ C ′ p C ′ p ) is directly associated with the vortex shedding activity (represented by ϕ p ′ p ′ for θ = π/2).It should be noted here that the correlation between ϕ p ′ p ′ (at θ = π/2) and ϕ C ′ p C ′ p is not obvious, given that the surface pressure at the crest does not contribute to the pressure drag along x.Such a correlation is only possible if the pressure drag is dependant on the vortex shedding phenomena, which is established by comparing the two spectra.Another noteworthy observation from figures 5(a,b) is that the energetic frequency range for both ϕ p ′ p ′ and ϕ C ′ p C ′ p remains consistent despite variation in the pressure drag domination across 6 ≤ Λ/k ≤ 60.It suggests that despite the variation from fully to transitionally rough scenarios on increasing Λ/k, the pressure drag mechanism always corresponds to the same frequency range (for the present case of k = 0.1H).
Next, on considering wall-oscillations imposed at T + osc = 100 (figures 5c-e), the corresponding oscillation frequency lies within the most energetic frequency range of ϕ p ′ p ′ and ϕ C ′ p C ′ p , for both Λ = 10k and 20k corresponding to k = 0.1H and 0.2H (see green arrows in figures 5c-e).One would thus expect the imposed wall oscillations to weaken the vortex shedding phenomena, likely by disturbing the two-dimensionality of the separated shear layer, and consequently reducing the pressure drag (Darekar & Sherwin 2001;Choi et al. 2008).This explains the significant attenuation of ϕ p ′ p ′ and ϕ C ′ p C ′ p , for the actuated cases (T + osc = 100) as compared to the static cases.The same was also noted in case of their cross-correlation spectra, ϕ C ′ p p ′ (not shown here).Hence, the choice of T + osc = 100 lying within the energetic frequency range, for all rough wall cases (figures 5a,b), explains the significant pressure drag reduction noted for these cases.This can be further confirmed by considering wall oscillations at a frequency much higher than the energetic range, for example T + osc = 10 (figure 5c).Indeed, both ϕ p ′ p ′ and ϕ C ′ p C ′ p spectra for T + osc = 10 overlap with the corresponding spectra for the static case at Λ = 10k, since such high oscillation frequency is incapable of disturbing the vortex shedding activity.This explains the negligible DR% corresponding to T + osc = 10.Along the same lines, in case of k = 0.2H (figure 5e), we can note that the attenuation of both the spectra (for actuated case relative to static case) is not as significant as noted for k = 0.1H (figure 5c), which again is consistent with the different magnitudes of DR% for the two cases (figure 2b).The results for k = 0.2H, thus, support the notion that the magnitude of pressure drag reduction is correlated with the intensity with which shedding activity in the roughness wake is weakened (depicted by ϕ p ′ p ′ ), and this is applicable despite the variation in roughness heights.The results also confirm that imposition of wall oscillation frequency, of the order of the vortex shedding frequencies, leads to pressure drag reduction.

Concluding remarks
The present study investigates drag reduction in a simplified (i.e.2-D) rough wall scenario, with the intention to answer the following two open questions: (i) Can the well-known viscous drag reduction strategy of imposing spanwise wall oscillations yield pressure drag reduction in a fully-rough case (i.e. with negligible viscous drag contributions)?(ii) What is the mechanism driving this pressure drag reduction?These are answered here via a set of systematically conducted direct numerical simulations of a turbulent flow over 2-D semi-cylindrical rods, maintained at matched roughness Reynolds numbers (k + ).The study also considers a case for higher roughness Reynolds number, to investigate the applicability of the pressure drag reduction mechanism across different roughness heights.It is found that spanwise wall oscillations are capable of attenuating the pressure drag when the wall oscillation frequency is of the order of frequency of vortex shedding from the rough elements.In all the rough wall cases simulated for the present study, the spectral distribution of the vortex shedding phenomena corresponded to a nominally similar range of viscous-scaled frequencies.This explained the pressure drag reduction on imposition of spanwise wall oscillations for all rough wall cases, given they were imposed at constant viscous-scaled actuation parameters.
Although the present study is limited to semi-cylindrical roughness geometry, the flow physics observed in association with the pressure drag reduction (namely, weakening of Reynolds stresses and increasing base pressure in the roughness wake) was found to be consistent with past observations noted in case of 2-D bumps (Banchetti et al. 2020) and square bar roughness (Nguyen et al. 2021).This supports applicability of the pressure drag reduction mechanism (i.e.weakening of vortex shedding) for varying roughness geometries, roughness heights and transitional/fully rough scenarios.A thorough quantification of this argument, however, would require an extensive simulation campaign, which is beyond the scope of the present study.
A natural follow-up of this work would involve testing the drag reduction strategy on practical surface conditions, such as sand-paper or 3-D roughness geometries, and at high Re τ .In this case, the pressure drag is likely to be contributed from highly 3-D vortical structures shedding at very small viscous-scaled frequencies, thereby requiring significant input power to directly affect the vortex shedding (and associated pressure drag).In this respect, the ability to affect this high-frequency phenomenon, via low-frequency wall oscillations (Marusic et  the inherent triadic interactions between the large-and small-time scales (Deshpande et al. 2023), could offer an energy-efficient pathway for rough wall drag reduction.

Appendix 1: Numerical parameters of the simulations
As discussed in §2, the present simulations are carried out using the open-source spectral element solver NEKTAR++ (Cantwell et al. 2015).In all the production simulations performed, the grid resolution per the cylindrical roughness element was kept identical.Each element was represented with Lagrange polynomials of degree N p = 7 on Gauss-Lobatto-Legendre points.The computational parameters of the simulation cases are listed in table 2. The size of the quadrilateral elements used to discretise the domain varies in space with the smallest element located at the wall next to the cylinder edge while the largest element is located at the top boundary (shown by the arrows in figure 6a).Denoting the streamwise and wall-normal size of the elements as l elem x and l elem z respectively, the streamwise and wall-normal grid resolutions within the smallest element can be estimated as ∆x + min = l elem+ x /(N p − 1) ≈ 1.38 and ∆z + min = l elem+ z /(N p − 1) ≈ 0.624.The corresponding grid resolution within the largest element are ∆x + max ≈ 4.98 and ∆z + max ≈ 7.29 (as shown in table 2 the values are slightly different for case Λ = 10k * in which k = 0.2H).The grid resolution in the homogeneous spanwise direction is kept constant at ∆y + ≈ 4.7.These grid resolutions are more than sufficient to resolve the near-wall turbulent scales of a smooth wall channel flow (Lozano-Durán & Jiménez 2014).To ascertain the adequacy of the grid resolution to resolve the turbulence associated with the vortex shedding downstream of the roughness elements, we have run two additional simulations (cases Λ = 10k f ine in table 2) that are identical to the static and actuated cases with Λ = 10k but at a higher grid resolution (see figure 6b).In figure 6(c), we compare the variation of the instantaneous total drag coefficient (C D ) as a function of time for both the static and actuated cases at the two grid resolutions.The match between the temporal evolution of the C D values from both grid resolutions indicates the adequacy of the adopted grid resolution to capture the dominant pressure drag reduction mechanisms at the actuation frequencies of the present study.Indeed, the mean drag coefficient values C D s = 36.1 × 10 −3 and C D a = 25.0 × 10 −3 for the static and actuated cases of Λ = 10k f ine are very close to the corresponding values of Λ = 10k (cf.table 1).A more quantitative assessment of the grid resolution is provided in figures 6(d,e) that show the wall-normal profiles of the global velocity fluctuation covariances (averaged over the corresponding steady-state interval).Note that the fluctuations include the turbulent, form-induced and wall-oscillation contributions.The very good match observed between the profiles is a further demonstration of the sufficiency of the considered resolution.
In order to test the adequacy of adopted box size to accommodate the dominant turbulent structures, we have computed the streamwise and spanwise pre-multiplied energy spectra of the streamwise velocity fluctuation on an x-y plane located above the crest of the cylinders where the mean turbulent kinetic energy attains a maximum value (within the shear layer region).Figure 7 shows the spectra for both the static and actuated cases of Λ = 10k and Λ = 60k.Although very large box size is required to achieve complete de-correlation in the streamwise direction, it can be seen that

Figure 1 :
Figure 1: (a) Schematics of an open-channel flow of an incompressible fluid over a rough wall, made up of spanwise-aligned semi-cylindrical rough elements.Terminology is defined in §2.The shaded light blue region shows the extent of the domain over what is referred to as the streamwise periodic unit (0 ≤ x/Λ ≤ 1).Dashed green lines are used to indicate the crest of the roughness elements, referred in §3.(b) Depiction of the two-dimensional quadrilateral spectral elements used to discretize the computational domain (shown only over a streamwise periodic unit).

Figure 2 :
Figure2: (a) Fractional contribution to the total drag experienced by a static rough wall (CD s ), from viscous forces (Cv s ; in blue shading) and pressure forces (Cp s ; in red shading), for cases of varying streamwise offsets between rough elements (Λ/k).(b) Total percentage drag reduction (in black diamonds), percent drag reduction due to decrease in pressure drag (in red squares), and due to decrease in viscous drag (in blue circles), on imposition of wall oscillations at T + osc = 100.In (a,b), fully and partially filled symbols correspond to present simulations with k = 0.1H and 0.2H respectively, while empty symbols in (a) correspond to data ofLeonardi et al. (2015).Dash-dotted magenta lines highlight cases: Λ = 10k and 20k analyzed in figures 3-5.
(b) also shows 100 × (C p s -C p a )/C D s and 100 × (C v s -C v a )/C D s for varying Λ/k.For significantly large streamwise offsets (10k < Λ ≤ 60k), which is the scenario investigated by Banchetti et al. (2020) and Nguyen et al. (2021), both viscous drag reduction and pressure drag reduction

Figure 5 :
Figure 5: Premultiplied frequency spectra of the fluctuating surface-pressure (f ϕ p ′ p ′ ; in black) at the crest of the rough element (θ = π/2), and of the fluctuating pressure force coefficient (f ϕ C ′ p C ′ p ; in red), as a function of the viscous-scaled frequency (f ν/U 2 τs ) for various rough wall scenarios.(a) and (b) respectively plot f ϕ p ′ p ′ and f ϕ C ′ p C ′ p for the static rough wall cases at 6 ≤ Λ/k ≤ 60. (c-e) compare both the spectra for the static and actuated cases at (c) Λ = 10k, k = 0.1H, (d) Λ = 20k, k = 0.1H, and (e) Λ = 10k, k = 0.2H.Green arrows are used to indicate (ToscU 2 τs /ν) on the x-axis, and the magenta shaded background nominally represents the most energetic frequency range for both the spectra (10 −2 ≲ f ν/U

Figure 6 :
Figure 6: Depiction of the spectral element resolution in the vicinity of the cylindrical roughness elements for (a) case Λ = 10k and (b) case Λ = 10k f ine .The arrows point to the smallest and largest elements.(c) Time evolution of the instantaneous total drag force coefficient for case Λ = 10k (blue) and Λ = 10k f ine (red).Solid lines show CD of the actuated case while the dashed lines show those of the corresponding static case.(d,e) Wall-normal profiles of the velocity fluctuation covariances (with respect to the global space and time averaged mean velocity U i) for the (d) static and (e) actuated cases corresponding to Λ = 10k (blue) and Λ = 10k f ine (red).

Table 1 :
A summary of direct numerical simulations for static and actuated cases with A + = A/U τs = 12 and oscillation period, T + osc = ToscU 2 τs /ν.Values marked with * correspond to the sole k = 0.2H case, while k = 0.1H for all other cases.

Table 2 :
al. 2021), by leveraging Numerical paramters of the simulations.The computational box size Lx × Ly × Lz = 6H × 3H × H is the same in all the simulations.N cyl = Lx/Λ represents the number of roughness or streamwise periodic units in the computational box.N elem is the number of two-dimensional spectral elements in the x-z plane, while Ny is the number of Fourier modes in the homogeneous spanwise direction.∆x + , ∆y + and ∆z + are representative grid resolutions (normalized with the static viscous length scale, ν/U τs ).T obs is the steady-state simulation interval over which statistics are accumulated after discarding an initial transient.Values marked with * correspond to the sole k = 0.2H case, while k = 0.1H for all other cases.