pith. sign in

arxiv: 2604.12074 · v1 · submitted 2026-04-13 · ⚛️ physics.flu-dyn

On the optimal period of spanwise wall forcing for turbulent drag reduction

Maurizio Quadrio , Federica Gattere , Marco Castelletti , Alessandro Chiarini This is my paper

Pith reviewed 2026-05-10 14:59 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords spanwise wall forcingturbulent drag reductionStokes layerdirect numerical simulationchannel flowwall oscillationsenergy saving
0
0 comments X

The pith

Adding an extra spanwise body force to wall oscillations decouples period from Stokes layer thickness and raises drag reduction by one third while turning net energy savings positive.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines turbulent channel flow under spanwise forcing to find better ways to cut skin-friction drag. Standard wall oscillations create a Stokes layer whose thickness is fixed by the oscillation period, and an optimal period exists but its meaning is unclear. By adding an independent spanwise body force, the thickness and period can be set separately, showing that pure wall oscillation is only one limited way to apply the forcing. This change produces higher drag reduction and positive net energy savings instead of losses. A reader would care because it reframes an established technique as suboptimal and points to more effective actuation methods for reducing energy use in flows.

Core claim

Harmonic wall oscillations generate a periodic transverse Stokes layer whose thickness is determined by the forcing period. Augmenting the wall oscillation with an additional spanwise body force makes the thickness and period independent. For the Reynolds numbers and forcing amplitudes examined, optimal performance occurs at substantially smaller periods and larger thicknesses than the classical Stokes layer, increasing maximum drag reduction by approximately one third and improving maximum net energy saving from negative 35 percent to positive 16 percent.

What carries the argument

Augmented spanwise forcing that combines wall oscillation with an independent spanwise body force, decoupling the Stokes layer thickness from the oscillation period.

Load-bearing premise

The extra spanwise body force can be realized physically at the needed strength without adding large extra energy costs or side effects, and the DNS results at the tested Reynolds numbers extend to other conditions.

What would settle it

Independent simulations or experiments that apply the reported optimal smaller period and larger thickness with the body force but find no increase in drag reduction or no shift to positive net energy saving would show the claim is not correct.

Figures

Figures reproduced from arXiv: 2604.12074 by Alessandro Chiarini, Federica Gattere, Marco Castelletti, Maurizio Quadrio.

Figure 1
Figure 1. Figure 1: Drag reduction R versus oscillation period 𝑇 + for the ESL forcing with 𝛿 = 𝛿SL (green) and the oscillating wall (black), at 𝐴 + = 12 and 𝑅𝑒𝜏 = 400. time units excluded from the computation of flow statistics to eliminate the initial transient. The oscillation period is varied within the range 10 ≤ 𝑇 + ≤ 200, while the thickness varies in 2 ≤ 𝛿 + ≤ 20. The forcing amplitude is fixed at 𝐴 + = 12. In total, … view at source ↗
Figure 2
Figure 2. Figure 2: Drag reduction map in the (𝑇, 𝛿) plane. The thick line indicates the 𝛿 = 𝛿SL (𝑇) constraint. The green/black dots identify the optimal point for ESL/SL. On the drag reduction map, R is largest over a relatively broad region, remaining above 30% for 20 ≤ 𝑇 + ≤ 50 and 6 ≤ 𝛿 + ≤ 20, and significant drag reduction is also observed for larger periods when 𝛿 + ≈ 8. Values of 𝛿 + exceeding the classical SL optimu… view at source ↗
Figure 3
Figure 3. Figure 3: Comparison between the optimal ESL (green) and SL (black) spanwise velocity profiles. Left: wall￾normal distribution of the spanwise velocity component; right: wall-normal distribution of the spanwise shear. The dashed lines mark the location of 𝛿𝑜 𝑝𝑡 for the two cases. −40 −40 − 40 −40 −20 −20 −20 −20 −20 −20 0 0 20 40 0 20 40 60 80 100 120 140 160 180 200 0 5 10 15 20 T + δ + −40 −20 0 20 40 ∆ww[%] [PIT… view at source ↗
Figure 4
Figure 4. Figure 4: Map in the (𝑇, 𝛿) plane of the integral change 𝛥𝑤𝑤 of the variance of the spanwise fluctuations with respect to the reference case. The thick line indicates the 𝛿 = 𝛿SL (𝑇) constraint. Note the inverted color scale. the channel of the variance of the spanwise velocity fluctuations around the spatial mean, minus the variance ⟨𝑤 2 ⟩0 of the spanwise velocity fluctuations in the reference flow, and expressed … view at source ↗
Figure 5
Figure 5. Figure 5: Net power saving map in the (𝑇, 𝛿) plane. The thick line indicates the 𝛿 = 𝛿SL (𝑇) constraint. The green/black dots identify the optimal point for ESL/SL. 4. Power budget The qualitative observation that the most effective spanwise forcing lies far from the SL line is further supported by considering the net energy saving potential of the ESL, evaluated under the usual assumption of ideal actuation (Baron … view at source ↗
Figure 6
Figure 6. Figure 6: Control power map in the (𝑇, 𝛿) plane: 𝑃𝑤 (left) and 𝑃𝑓 (right) expressed as percentages of the pumping power 𝑃0. The thick line indicates the 𝛿 = 𝛿SL (𝑇) constraint. dominant contribution. Moreover, although the instantaneous power associated with the forcing terms in (2.2) containing the ESL profile can be large, the periodicity of the ESL implies that their time average vanishes. As a result, the last t… view at source ↗
Figure 7
Figure 7. Figure 7: Snapshots of the field of streamwise velocity fluctuations in a wall-parallel plane at 𝑦 ∗ = 14.5. Top left: reference flow; top right: optimal SL; bottom left: optimal ESL; bottom right: suboptimal ESL. All controlled cases are at the same (zero) phase [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Longitudinal mean velocity profiles (top) and r.m.s. value of velocity fluctuations (bottom): optimal SL and ESL are compared to the reference case and to a suboptimal ESL case at large 𝑇. actuators (Gouder et al. 2013; Gatti et al. 2015), alternating slip/no-slip stripes (Fuaad & Arul Prakash 2019), surface dimples (Gattere et al. 2022), sinusoidal riblets (Peet et al. 2008). However, their design could b… view at source ↗
read the original abstract

Turbulent channel flow controlled by spanwise wall oscillations is studied using direct numerical simulations to improve how spanwise forcing reduces skin-friction drag. Harmonic wall oscillations generate a periodic transverse Stokes layer whose thickness $\delta$ is determined by the forcing period $T$. Although an optimal $T$ that maximizes drag reduction is known to exist, its physical significance remains unclear. To elucidate it, we extend the spanwise Stokes layer by augmenting wall oscillation with an additional spanwise body force. In this formulation, $\delta$ and $T$ become decoupled and can be varied independently. The oscillating wall thus appears as a special and suboptimal case of spanwise forcing. Optimal performance is obtained for substantially smaller $T$ and larger $\delta$ than those of the classical Stokes layer. For the conditions examined, with Reynolds number and forcing amplitude held fixed, the maximum drag reduction increases by approximately one third, while the maximum net energy saving improves markedly from $-35\%$ to $+16\%$. These findings suggest that drag-reduction strategies based on spanwise forcing deserve renewed scrutiny: wall oscillation represents only one possible actuation method, and not necessarily the most effective one.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper uses direct numerical simulations of turbulent channel flow to study spanwise forcing for skin-friction drag reduction. By augmenting harmonic wall oscillations with an additional spanwise body force, the Stokes-layer thickness δ is decoupled from the forcing period T. The authors report that optimal performance occurs at substantially smaller T and larger δ than in the classical Stokes layer. For fixed Reynolds number and forcing amplitude, this yields approximately one-third higher maximum drag reduction and improves the maximum net energy saving from −35% to +16%.

Significance. If the numerical results are robust, the work shows that classical wall oscillation is a suboptimal special case of spanwise forcing and that independent control of δ and T can produce substantially better drag reduction and net energy savings. This is a useful parameter-space exploration via DNS and could motivate new actuation concepts. Credit is due for the clean decoupling approach that isolates the roles of δ and T without introducing fitted parameters.

major comments (2)
  1. [Numerical methods] Numerical methods section: the manuscript provides no grid-resolution details (e.g., Δx+, Δy+, Δz+), time-step size, or statistical-convergence criteria (e.g., averaging time in outer units or uncertainty estimates on drag and power integrals). These quantities are load-bearing for the central claim of a one-third DR increase and sign change in net saving.
  2. [§3–4] Formulation of the body force and net-power calculation (likely §3 and §4): the additional spanwise body force is introduced to decouple δ and T, yet the explicit spatial/temporal profile and the precise expression used for its power cost (volume integral of f·u) are not stated. Because net energy saving is defined as the difference between drag-reduction benefit and total power input, any ambiguity here directly affects the reported improvement from −35% to +16%.
minor comments (2)
  1. [Abstract] Abstract: the phrase “one third” should be replaced by the exact percentage or a reference to the figure/table that quantifies it.
  2. [Figures] Figure captions (e.g., those showing DR and net-saving contours versus T and δ): add a brief note on the number of independent simulations and the range of T and δ scanned to support the location of the reported optimum.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and positive assessment of the work's significance. We address the two major comments point by point below. Both points identify omissions that affect reproducibility and clarity; we will revise the manuscript to incorporate the requested details.

read point-by-point responses

  1. Referee: [Numerical methods] Numerical methods section: the manuscript provides no grid-resolution details (e.g., Δx+, Δy+, Δz+), time-step size, or statistical-convergence criteria (e.g., averaging time in outer units or uncertainty estimates on drag and power integrals). These quantities are load-bearing for the central claim of a one-third DR increase and sign change in net saving.

    Authors: We agree that these details are essential for assessing the robustness of the reported drag-reduction and net-energy-saving results. In the revised manuscript we will expand the numerical-methods section to include the grid spacings in wall units (Δx+, Δy+, Δz+), the time-step size (in viscous units), the total averaging time in outer units, and any convergence checks or uncertainty estimates performed on the drag and power integrals. revision: yes

  2. Referee: [§3–4] Formulation of the body force and net-power calculation (likely §3 and §4): the additional spanwise body force is introduced to decouple δ and T, yet the explicit spatial/temporal profile and the precise expression used for its power cost (volume integral of f·u) are not stated. Because net energy saving is defined as the difference between drag-reduction benefit and total power input, any ambiguity here directly affects the reported improvement from −35% to +16%.

    Authors: We accept that the explicit form of the body force and the precise power-cost integral must be stated unambiguously. In the revised manuscript we will add, in §3, the full spatial and temporal expression for the additional spanwise body force, and in §4 we will write out the volume integral used for its power contribution (∫ f·u dV) together with the corresponding term for the wall-oscillation power. These additions will remove any ambiguity in the net-energy-saving calculation. revision: yes

Circularity Check

0 steps flagged

No significant circularity; central claims rest on independent DNS parameter sweeps

full rationale

The paper reports drag-reduction and net-energy-saving values obtained by running separate direct numerical simulations at fixed Reynolds number and forcing amplitude while independently varying the Stokes-layer thickness δ and period T through an added spanwise body force. No fitted parameters are later renamed as predictions, no self-citations supply load-bearing uniqueness theorems, and the reported optima (smaller T, larger δ) are simply the numerically observed maxima among the simulated cases. The derivation chain therefore contains no self-definitional, fitted-input, or self-citation reductions; the performance numbers are direct outputs of the simulations rather than algebraic rearrangements of the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claim rests on the accuracy of the incompressible DNS and on the assumption that the added body force can be applied without additional modeling errors.

axioms (2)
  • standard math The flow obeys the incompressible Navier-Stokes equations
    Standard governing equations for DNS of low-Mach-number channel flow at the Reynolds numbers considered.
  • domain assumption Time- and spanwise-averaged statistics yield a well-defined mean skin-friction drag
    Required to compute the reported drag-reduction percentages from the simulations.
invented entities (1)
  • additional spanwise body force no independent evidence
    purpose: To control Stokes-layer thickness δ independently of forcing period T
    Modeling device introduced in this work to explore a larger parameter space; no independent physical mechanism or external validation is supplied.

pith-pipeline@v0.9.0 · 5511 in / 1333 out tokens · 64397 ms · 2026-05-10T14:59:15.850924+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year eprint label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence ...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

    work page

  3. [3]

    , Touber, E

    Agostini, L. , Touber, E. & Leschziner, M.A. 2014 Spanwise oscillatory wall motion in channel flow: Drag-reduction mechanisms inferred from DNS-predicted phase-wise property variations at Re =1000 . J. Fluid Mech. 743 , 606--635

    work page 2014

  4. [4]

    & Quadrio, M

    Baron, A. & Quadrio, M. 1996 Turbulent drag reduction by spanwise wall oscillations . Appl. Sci. Res. 55 , 311--326

    work page 1996

  5. [5]

    , Bayoda, K.D

    Benard, N. , Bayoda, K.D. , Coma, M. , Pons Prat, J. , Bonnet, J.P. & Moreau, E. 2021 Turbulent skin friction reduction by spatial wall forcing oscillations with non-thermal plasma. In 55th 3AF International Conference on Applied Aerodynamics \/ . Poitiers (France)

    work page 2021

  6. [6]

    , Xu, C.-X

    Choi, J.-I. , Xu, C.-X. & Sung, H. J. 2002 Drag reduction by spanwise wall oscillation in wall-bounded turbulent flows . AIAA J. 40 (5), 842--850

    work page 2002

  7. [7]

    2002 Near-wall structure of turbulent boundary layer with spanwise-wall oscillation

    Choi, K.-S. 2002 Near-wall structure of turbulent boundary layer with spanwise-wall oscillation . Phys. Fluids 14 (7), 2530--2542

    work page 2002

  8. [8]

    & Karniadakis, G

    Du, Y. & Karniadakis, G. E. 2000 Suppressing wall turbulence by means of a transverse traveling wave . Science 288 , 1230--1234

    work page 2000

  9. [9]

    & Arul Prakash, K

    Fuaad, P.A. & Arul Prakash, K. 2019 Enhanced drag-reduction over superhydrophobic surfaces with sinusoidal textures: A DNS study . Comput. Fluids 181 , 208--223

    work page 2019

  10. [10]

    , Chiarini, A

    Gattere, F. , Chiarini, A. & Quadrio, M. 2022 Dimples for Skin-Friction Drag Reduction : Status and Perspectives . Fluids 7 (7), 240

    work page 2022

  11. [11]

    , Zanolini, M

    Gattere, F. , Zanolini, M. , Gatti, D. , Bernardini, M. & Quadrio, M. 2024 Turbulent drag reduction with streamwise-travelling waves in the compressible regime . J. Fluid Mech. 987 , A30

    work page 2024

  12. [12]

    , G \"u ttler, A

    Gatti, D. , G \"u ttler, A. , Frohnapfel, B. & Tropea, C. 2015 Experimental assessment of spanwise-oscillating dielectric electroactive surfaces for turbulent drag reduction in an air channel flow . Exp. Fluids 56 (5), 1--15

    work page 2015

  13. [13]

    & Quadrio, M

    Gatti, D. & Quadrio, M. 2016 Reynolds-number dependence of turbulent skin-friction drag reduction induced by spanwise forcing . J. Fluid Mech. 802 , 553--58

    work page 2016

  14. [14]

    , Quadrio, M

    Gatti, D. , Quadrio, M. , Chiarini, A. , Gattere, F. & Pirozzoli, S. 2025 Turbulent skin-friction drag reduction via spanwise forcing at high Reynolds number . J. Fluid Mech. 1016

    work page 2025

  15. [15]

    , Potter, M

    Gouder, K. , Potter, M. & Morrison, J.F. 2013 Turbulent friction drag reduction using electroactive polymer and electromagnetically driven surfaces . Exp. Fluids 54 (1), 1441

    work page 2013

  16. [16]

    2013 Near-wall turbulence

    Jim \'e nez, J. 2013 Near-wall turbulence . Phys. Fluids 25 (10), 101302

    work page 2013

  17. [17]

    & Choi, K.-S

    Jukes, T.N. & Choi, K.-S. 2012 Dielectric-barrier-discharge vortex generators: Characterisation and optimisation for flow separation control . Exp Fluids 52 , 329--345

    work page 2012

  18. [18]

    , Mangiavacchi, N

    Jung, W.J. , Mangiavacchi, N. & Akhavan, R. 1992 Suppression of turbulence in wall-bounded flows by high-frequency spanwise oscillations . Phys. Fluids A 4 (8) , 1605--1607

    work page 1992

  19. [19]

    , Scarano, F

    Kempaiah, K.U. , Scarano, F. , Elsinga, G.E. , van Oudheusden , B.W. & Bermel, Leon 2020 3-dimensional particle image velocimetry based evaluation of turbulent skin-friction reduction by spanwise wall oscillation . Phys. Fluids 32 (8), 085111

    work page 2020

  20. [20]

    & Hussain, F

    Kim, J. & Hussain, F. 1993 Propagation velocity of perturbations in turbulent channel flow . Phys. Fluids A 5 (3), 695--706

    work page 1993

  21. [21]

    , Skandaji, L

    Laadhari, F. , Skandaji, L. & Morel, R. 1994 Turbulence reduction in a boundary layer by a local spanwise oscillating surface . Phys. Fluids 6 (10) , 3218--3220

    work page 1994

  22. [22]

    2021 Introducing CPL , arXiv:arXiv: 2021.12143v2

    Luchini, P. 2021 Introducing CPL , arXiv:arXiv: 2021.12143v2

    work page arXiv 2021

  23. [23]

    & Quadrio, M

    Luchini, P. & Quadrio, M. 2006 A low-cost parallel implementation of direct numerical simulation of wall turbulence . J. Comp. Phys. 211 (2), 551--571

    work page 2006

  24. [24]

    , Chandran, D

    Marusic, I. , Chandran, D. , Rouhi, A. , Fu, M.K. , Wine, D. , Holloway, B. , Chung, D. & Smits, A.J. 2021 An energy-efficient pathway to turbulent drag reduction . Nat. Commun. 12 (1), 5805

    work page 2021

  25. [25]

    , Sagaut, P

    Peet, Y. , Sagaut, P. & Charron, Y. 2008 Turbulent Drag Reduction using Sinusoidal Riblets with Triangular Cross-Section . AIAA Pap. 2008-3745

    work page 2008

  26. [26]

    2011 Drag reduction in turbulent boundary layers by in-plane wall motion

    Quadrio, M. 2011 Drag reduction in turbulent boundary layers by in-plane wall motion . Phil. Trans. R. Soc. A 369 (1940), 1428--1442

    work page 2011

  27. [27]

    , Frohnapfel, B

    Quadrio, M. , Frohnapfel, B. & Hasegawa, Y. 2016 Does the choice of the forcing term affect flow statistics in DNS of turbulent channel flow? Eur. J. Mech. B / Fluids 55 , 286--293

    work page 2016

  28. [28]

    & Luchini, P

    Quadrio, M. & Luchini, P. 2003 Integral time-space scales in turbulent wall flows . Phys. Fluids 15 (8), 2219--2227

    work page 2003

  29. [29]

    & Ricco, P

    Quadrio, M. & Ricco, P. 2004 Critical assessment of turbulent drag reduction through spanwise wall oscillation . J. Fluid Mech. 521 , 251--271

    work page 2004

  30. [30]

    & Ricco, P

    Quadrio, M. & Ricco, P. 2011 The laminar generalized Stokes layer and turbulent drag reduction . J. Fluid Mech. 667 , 135--157

    work page 2011

  31. [31]

    , Ricco, P

    Quadrio, M. , Ricco, P. & Viotti, C. 2009 Streamwise-traveling waves of spanwise wall velocity for turbulent drag reduction . J. Fluid Mech. 627 , 161--178

    work page 2009

  32. [32]

    & Sibilla, S

    Quadrio, M. & Sibilla, S. 2000 Numerical simulation of turbulent flow in a pipe oscillating around its axis . J. Fluid Mech. 424 , 217--241

    work page 2000

  33. [33]

    2004 Modification of near-wall turbulence due to spanwise wall oscillations

    Ricco, P. 2004 Modification of near-wall turbulence due to spanwise wall oscillations . J. Turbul. 5 (24)

    work page 2004

  34. [34]

    , Skote, M

    Ricco, P. , Skote, M. & Leschziner, M. A. 2021 A review of turbulent skin-friction drag reduction by near-wall transverse forcing . Prog. Aero. Sci. 123 , 100713

    work page 2021

  35. [35]

    , Corke, T.C

    Thomas, F.O. , Corke, T.C. , Duong, A. , Midya, S. & Yates, K. 2019 Turbulent drag reduction using pulsed- DC plasma actuation . J. Phys. D: Appl. Phys. 52 (43), 434001

    work page 2019

  36. [36]

    & Leschziner, M.A

    Touber, E. & Leschziner, M.A. 2012 Near-wall streak modification by spanwise oscillatory wall motion and drag-reduction mechanisms . J. Fluid Mech. 693 , 150--200

    work page 2012

  37. [37]

    , Bogard, D.G

    Trujillo, S.M. , Bogard, D.G. & Ball, K.S. 1997 Turbulent Boundary Layer Drag Reduction using an Oscillating Wall . AIAA Pap. 97-1870

    work page 1997