pith. sign in

arxiv: 2604.22150 · v1 · submitted 2026-04-24 · ⚛️ physics.flu-dyn

Waves dictate the yo-yoing decay of a viscoelastic mixing layer

Giulio Foggi Rota , Piyush Garg , Jason Tang , Marco Edoardo Rosti This is my paper

Pith reviewed 2026-05-08 10:23 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords viscoelastic fluidmixing layerwavesyo-yoingenergy budgetnonlinear dispersiondirect numerical simulation
0
0 comments X

The pith

Waves develop in decaying viscoelastic mixing layers and drive oscillatory yo-yoing of the mean flow.

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

The paper establishes that a time-decaying mixing layer in a viscoelastic fluid develops waves that cause the mean flow velocity to oscillate up and down instead of decaying monotonically as it does in Newtonian fluids. This yo-yoing arises because elastic polymers store and release energy while being rotated by the large-scale shear, creating a feedback that injects energy back into the fluid at regular intervals. Direct numerical simulations reveal the phenomenon, and an energy-budget analysis isolates the polymer contribution as the driver. An analytical model then yields exact wave solutions whose nonlinear dispersion relation predicts both the oscillation period and the range of parameters where yo-yoing occurs. The result supplies a simple explanation for recently observed anomalies in unsteady viscoelastic experiments.

Core claim

In a viscoelastic mixing layer the mean flow does not decay steadily; instead it yo-yos because waves propagate through the layer, with elastic polymers injecting energy into the fluid and being rotated by the mean shear in a repeating cycle. Analytical solutions of the linearized problem admit wave modes with a nonlinear dispersion relation that fixes the yo-yoing period and the domain of parameters in which the oscillation persists.

What carries the argument

Analytical wave solutions with nonlinear dispersion that couple polymer stress to mean shear and close the energy budget of the mixing layer.

If this is right

  • The period of mean-flow oscillation is set by the nonlinear dispersion relation and can be predicted without solving the full time-dependent equations.
  • Yo-yoing occurs only inside a bounded region of Reynolds number and Weissenberg number; outside that region the layer decays monotonically.
  • The same polymer-shear feedback mechanism accounts for unsteady anomalies seen in more complex viscoelastic geometries.
  • Energy is periodically transferred from polymers back to the mean flow, violating the monotonic dissipation expected for Newtonian fluids.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar wave-driven oscillations may appear in other canonical unsteady shear flows once polymers are added, such as decaying wakes or jets.
  • The analytical dispersion relation offers a quick way to estimate oscillation periods in industrial mixing processes that involve viscoelastic liquids.
  • If the waves persist at higher Reynolds numbers they could alter mixing efficiency and scalar transport in viscoelastic turbulence.

Load-bearing premise

The viscoelastic constitutive model and the numerical scheme together faithfully represent polymer dynamics without introducing artifacts that artificially sustain the waves.

What would settle it

A direct numerical simulation or experiment in which the same initial mixing layer is evolved with a different polymer model or with polymers removed shows strictly monotonic decay of the mean velocity profile.

Figures

Figures reproduced from arXiv: 2604.22150 by Giulio Foggi Rota, Jason Tang, Marco Edoardo Rosti, Piyush Garg.

Figure 1
Figure 1. Figure 1: FIG. 1 view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 view at source ↗
read the original abstract

We find that waves develop in a time-decaying mixing layer of viscoelastic fluid, leading the mean-flow to yo-yo. This is in sharp contrast with Newtonian fluids, where laminar mixing layers evolve monotonically. We combine direct numerical simulations with a theoretical analysis of the energy budget for the flow to uncover the underlying physical mechanism. The yo-yoing of the mean-flow is shown to be driven by the elastic polymers injecting energy into the fluid and, in turn, being rotated by the large-scale mean shear. We then provide the mathematical model of the problem and solve it analytically, finding wave solutions with non-linear dispersion predicting the period of the yo-yoing and the parameter range where it occurs. As decaying mixing layers are one of the simplest and canonical examples of unsteady flows, the phenomenon identified here explains the anomalies recently observed in experiments of unsteady viscoelastic flows in complex geometries.

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

3 major / 2 minor

Summary. The manuscript reports that waves develop in a time-decaying mixing layer of a viscoelastic fluid, causing the mean flow to exhibit yo-yoing (oscillatory) decay, in contrast to the monotonic evolution seen in Newtonian fluids. Using direct numerical simulations (DNS) combined with energy-budget analysis, the authors identify the mechanism as elastic polymers injecting energy into the fluid while being rotated by the large-scale mean shear. They then derive an analytical model yielding wave solutions with nonlinear dispersion that predict both the yo-yoing period and the parameter range in which the phenomenon occurs. The work is positioned as explaining anomalies in recent experiments on unsteady viscoelastic flows in complex geometries.

Significance. If the central claim holds, the result would be significant for the field of viscoelastic fluid dynamics: it identifies a previously unrecognized wave-mediated mechanism that governs the decay of a canonical unsteady flow (the mixing layer) and links polymer elasticity directly to mean-flow oscillations. The combination of DNS, energy-budget diagnostics, and an analytical model with falsifiable predictions (period and parameter window) is a methodological strength that could enable direct tests against future experiments. This could resolve discrepancies in complex-geometry viscoelastic flows and inform modeling of industrial mixing or polymer processing under unsteady conditions.

major comments (3)
  1. [Numerical Methods / DNS Setup] The specific constitutive equation used for the polymer stress in the DNS (e.g., Oldroyd-B, FENE-P, or Giesekus) is never stated, nor are the values of the Weissenberg number, viscosity ratio, or Reynolds number. This information is load-bearing because the energy-injection mechanism and the existence of the waves depend on the precise form of the stress tensor; without it, neither reproducibility nor the absence of model-specific artifacts can be assessed.
  2. [Analytical Model / Wave Solutions] The analytical wave solutions with nonlinear dispersion (presented as predicting the yo-yoing period) rely on derivation steps and closure assumptions that are not shown in sufficient detail. It is unclear whether the model linearizes the polymer stress or assumes a fixed mean-shear profile; if either assumption breaks once the mean flow begins to oscillate at finite amplitude, the predicted period will not correspond to the DNS. A direct, quantitative comparison of the analytically predicted period versus the DNS-measured period across multiple parameter values is required to substantiate the claim.
  3. [Energy Budget Analysis] The energy-budget analysis asserts that polymers inject energy while being rotated by the mean shear, yet no explicit term-by-term decomposition of the kinetic-energy or elastic-energy equations is provided, nor is any grid-convergence or stress-tensor accuracy check reported. These omissions leave open the possibility that numerical dissipation or discretization artifacts, rather than the physical mechanism, drive the observed waves and yo-yoing.
minor comments (2)
  1. [Figures] Figure captions should explicitly list the dimensionless parameters (Re, Wi, etc.) corresponding to each panel so that the reader can connect the visualized yo-yoing behavior to the analytical predictions.
  2. [Analytical Model] Notation for the wave speed and dispersion relation in the analytical section should be defined once and used consistently; currently some symbols appear without prior definition.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the detailed and constructive review. We address each major comment below and will incorporate the suggested clarifications and additions into the revised manuscript.

read point-by-point responses

  1. Referee: [Numerical Methods / DNS Setup] The specific constitutive equation used for the polymer stress in the DNS (e.g., Oldroyd-B, FENE-P, or Giesekus) is never stated, nor are the values of the Weissenberg number, viscosity ratio, or Reynolds number. This information is load-bearing because the energy-injection mechanism and the existence of the waves depend on the precise form of the stress tensor; without it, neither reproducibility nor the absence of model-specific artifacts can be assessed.

    Authors: We agree that these details are essential for reproducibility. The DNS were performed with the Oldroyd-B constitutive model at Weissenberg number Wi = 8, viscosity ratio β = 0.6, and Reynolds number Re = 1200. In the revised manuscript we will add an explicit subsection in the numerical methods section stating the constitutive equation, all parameter values, and their ranges, together with a brief justification of the model choice. revision: yes

  2. Referee: [Analytical Model / Wave Solutions] The analytical wave solutions with nonlinear dispersion (presented as predicting the yo-yoing period) rely on derivation steps and closure assumptions that are not shown in sufficient detail. It is unclear whether the model linearizes the polymer stress or assumes a fixed mean-shear profile; if either assumption breaks once the mean flow begins to oscillate at finite amplitude, the predicted period will not correspond to the DNS. A direct, quantitative comparison of the analytically predicted period versus the DNS-measured period across multiple parameter values is required to substantiate the claim.

    Authors: We will expand the derivation of the analytical model in a new appendix, explicitly showing every step and stating the closure assumptions (linearized polymer stress about the instantaneous mean shear, with the mean profile allowed to evolve slowly). We will also add a direct quantitative comparison: a table and accompanying figure reporting the analytically predicted oscillation periods against the periods extracted from DNS for at least four distinct parameter combinations, thereby testing the validity of the assumptions at finite amplitude. revision: yes

  3. Referee: [Energy Budget Analysis] The energy-budget analysis asserts that polymers inject energy while being rotated by the mean shear, yet no explicit term-by-term decomposition of the kinetic-energy or elastic-energy equations is provided, nor is any grid-convergence or stress-tensor accuracy check reported. These omissions leave open the possibility that numerical dissipation or discretization artifacts, rather than the physical mechanism, drive the observed waves and yo-yoing.

    Authors: We will include a new appendix that presents the complete term-by-term decomposition of both the kinetic-energy and elastic-energy budgets, with each contribution evaluated from the DNS fields. In addition, we will report grid-convergence results (comparing three successively refined meshes) and verification checks on the divergence-free condition of the stress tensor to confirm that the observed waves and energy injection are not numerical artifacts. revision: yes

Circularity Check

0 steps flagged

No significant circularity; analytical wave solutions are independent of simulation outputs.

full rationale

The paper reports DNS observations of yo-yoing in viscoelastic mixing layers, followed by an energy-budget analysis to identify the polymer-shear mechanism, then introduces a separate mathematical model that is solved analytically to obtain wave solutions with nonlinear dispersion. These solutions are presented as predicting the observed period and the parameter window in which yo-yoing occurs. No quoted derivation step shows that the dispersion relation, wave amplitude, or selected parameters are obtained by fitting to the yo-yoing period or mean-flow oscillation extracted from the same DNS runs; the analytical treatment is constructed from the constitutive equations and linearized or weakly nonlinear assumptions rather than by construction from the simulation data. Self-citations, if present, are not load-bearing for the central prediction. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard viscoelastic fluid modeling assumptions and the validity of energy budget analysis in DNS, but the abstract does not specify free parameters, additional axioms, or new entities.

axioms (1)
  • domain assumption Viscoelastic fluid behavior is captured by a standard constitutive model involving elastic polymers that can store and release energy.
    Implied by references to elastic polymers, energy injection, and rotation by mean shear in the mechanism description.

pith-pipeline@v0.9.0 · 5453 in / 1274 out tokens · 69155 ms · 2026-05-08T10:23:33.706248+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

62 extracted references

  1. [1]

    0 0 . 2 0 . 4 ⟨Cxy⟩x − 14 0 14 y/δ 0 (d) 0 200 400 tU0/δ 0 2 3 4 5 ⟨tr(C)⟩c (c) 0 200 400 tU0/δ 0 10−40 10−20 100 (e) Vδ0/U 3 0 Rδ0/U 3 0 Pδ0/U 3 0 FIG. 1.Y o-yoing of the viscoelastic mixing layer.(a) Evolution of the absolute value of the mean centreline vorticity|⟨ω⟩ c/⟨ω0⟩c|with the non-dimensional timetU 0/δ0 for the viscoelastic (De= 28) and Newtoni...

  2. [2]

    4 ˜kn 0 1 2 ˜ω n (c)

    2 0 . 4 ˜kn 0 1 2 ˜ω n (c)

  3. [3]

    0 0 . 5 1 . 0 β 0 100 200 300E (d) 5 15 40 100 T /τ

  4. [4]

    0 0 . 5 1 . 0 β 0 100 200 300E (e)

  5. [5]

    5 2 10 50 t∗/τ FIG. 3.Viscoelastic waves.Evolution of the mean centreline vorticity|⟨ω⟩ c/⟨ω0⟩c|over timetU 0/δ0 for varying (a) Deborah numbers (De) withRe= 0.2 andβ= 0.9, and (b) viscosity ratios (β) withRe= 0.2 andDe= 28. Continuous lines denote positive values, dashed lines negative values. (c) Dispersion relation for the viscoelastic waves periodical...

  6. [6]

    Unpacking

    of the form ˜u(˜y,˜t) = +∞X n=0 (aneγ+(˜kn)˜t +b neγ−(˜kn)˜t)Bn(˜y),(7) whereB n is a set of sinusoidal base–functions compatible with the boundary conditions, while the coefficients an, bn follow from the initial conditions, projected over the infinite set of discrete modes with wavenumbers ˜kn = (2n+ 1)π/L(further details on the derivation are found in ...

  7. [7]

    R. B. Bird, R. C. Armstrong, and O. Hassager,Dynamics of polymeric liquids. Vol. 1: Fluid mechanics(John Wiley and Sons Inc., New York, NY, 1987)

    1987

  8. [8]

    Balmforth and R

    N. Balmforth and R. Craster, Geomorphological Fluid Mechanics , 34 (2001)

    2001

  9. [9]

    Morozov and S

    A. Morozov and S. E. Spagnolie, Complex fluids in biological systems: experiment, theory, and computation , 3 (2015)

    2015

  10. [10]

    R. G. Larson,Constitutive equations for polymer melts and solutions: Butterworths series in chemical engineering (Butterworth-Heinemann, 2013). 8

    2013

  11. [11]

    S. S. Datta, A. M. Ardekani, P. E. Arratia, A. N. Beris, I. Bischofberger, G. H. McKinley, J. G. Eggers, J. E. L´ opez-Aguilar, S. M. Fielding, A. Frishman,et al., Physical Review Fluids7, 080701 (2022)

    2022

  12. [12]

    Steinberg, Annual Review of Fluid Mechanics53, 27 (2021)

    V. Steinberg, Annual Review of Fluid Mechanics53, 27 (2021)

    2021

  13. [13]

    E. G. Shaqfeh, Annual Review of Fluid Mechanics , 129 (1996)

    1996

  14. [14]

    Bergeron, D

    V. Bergeron, D. Bonn, J. Y. Martin, and L. Vovelle, Nature405, 772 (2000)

    2000

  15. [15]

    Varshney and V

    A. Varshney and V. Steinberg, Nature communications10, 652 (2019)

    2019

  16. [16]

    de Blois, S

    C. de Blois, S. J. Haward, and A. Q. Shen, Physical Review Fluids8, 023301 (2023)

    2023

  17. [17]

    Groisman and V

    A. Groisman and V. Steinberg, Nature405, 53 (2000)

    2000

  18. [18]

    R. K. Singh, P. Perlekar, D. Mitra, and M. E. Rosti, Nature Communications15, 4070 (2024)

    2024

  19. [19]

    L. Pan, A. Morozov, C. Wagner, and P. Arratia, Physical Review Letters110, 174502 (2013)

    2013

  20. [20]

    Lellep, M

    M. Lellep, M. Linkmann, and A. Morozov, Proceedings of the National Academy of Sciences121, e2318851121 (2024)

    2024

  21. [21]

    Foggi Rota, C

    G. Foggi Rota, C. Amor, S. Le Clainche, and M. E. Rosti, Physical Review Fluids9, L122602 (2024)

    2024

  22. [22]

    Boffetta, A

    G. Boffetta, A. Celani, A. Mazzino, A. Puliafito, and M. Vergassola, J. Fluid Mech.523, 161 (2005)

    2005

  23. [23]

    P. Garg, I. Chaudhary, M. Khalid, V. Shankar, and G. Subramanian, Physical Review Letters121, 024502 (2018)

    2018

  24. [24]

    J. Page, Y. Dubief, and R. R. Kerswell, Physical Review Letters125, 154501 (2020)

    2020

  25. [25]

    Morozov, Physical Review Letters129, 017801 (2022)

    A. Morozov, Physical Review Letters129, 017801 (2022)

    2022

  26. [26]

    Lewy and R

    T. Lewy and R. R. Kerswell, Journal of Fluid Mechanics1007, A55 (2025)

    2025

  27. [27]

    Ho and P

    C.-M. Ho and P. Huerre, Annual Review of Fluid Mechanics16, 365 (1984)

    1984

  28. [28]

    Peltier and C

    W. Peltier and C. Caulfield, Annual Review of Fluid Mechanics35, 135 (2003)

    2003

  29. [29]

    Azaiez and G

    J. Azaiez and G. Homsy, Journal of Fluid Mechanics268, 37 (1994)

    1994

  30. [30]

    Azaiez and G

    J. Azaiez and G. Homsy, Journal of non-newtonian fluid mechanics52, 333 (1994)

    1994

  31. [31]

    P. K. Ray and T. A. Zaki, Physics of Fluids26(2014)

    2014

  32. [32]

    Varshney and V

    A. Varshney and V. Steinberg, Physical Review Fluids3, 103303 (2018)

    2018

  33. [33]

    N. K. Jha and V. Steinberg, Proceedings of the National Academy of Sciences118, e2105211118 (2021)

    2021

  34. [34]

    R. A. Lopez de la Cruz, S. J. Haward, and A. Q. Shen, PNAS nexus4, pgae571 (2025)

    2025

  35. [35]

    Yu and N

    Z. Yu and N. Phan-Thien, Journal of Fluid Mechanics500, 29 (2004)

    2004

  36. [36]

    Hinch and O

    J. Hinch and O. Harlen, Journal of Non-Newtonian Fluid Mechanics298, 104668 (2021)

    2021

  37. [37]

    Hillebrand, R

    F. Hillebrand, R. J. Hill, M. Davoodi, S. J. Haward, A. Q. Shen, R. J. Poole, and S. Varchanis, Phys. Rev. Lett. 135, 024002 (2025)

    2025

  38. [38]

    Foggi Rota, R

    G. Foggi Rota, R. K. Singh, A. Chiarini, C. Amor, G. Soligo, D. Mitra, and M. E. Rosti, Int. J. Multiph. Flow 197, 105630 (2026)

    2026

  39. [39]

    V. F. Zaitsev and A. D. Polyanin,Handbook of exact solutions for ordinary differential equations(Chapman and Hall/CRC, 2002)

    2002

  40. [40]

    Achenbach,Wave propagation in elastic solids(Elsevier, 2012)

    J. Achenbach,Wave propagation in elastic solids(Elsevier, 2012)

    2012

  41. [41]

    Caulfield, Annual Review of Fluid Mechanics53, 113 (2021)

    C. Caulfield, Annual Review of Fluid Mechanics53, 113 (2021)

    2021

  42. [42]

    Gilpin, M

    W. Gilpin, M. S. Bull, and M. Prakash, Nature Reviews Physics2, 74 (2020)

    2020

  43. [43]

    A. E. Hosoi, Physical Review Fluids4, 110508 (2019)

    2019

  44. [44]

    M. E. Cates and S. M. Fielding, Advances in Physics55, 799 (2006)

    2006

  45. [45]

    Fardin, D

    M.-A. Fardin, D. Lopez, J. Croso, G. Gr´ egoire, O. Cardoso, G. H. McKinley, and S. Lerouge, Physical review letters104, 178303 (2010)

    2010

  46. [46]

    J. E. Butler and B. Snook, Annual Review of Fluid Mechanics50, 299 (2018)

    2018

  47. [47]

    Guazzelli and J

    E. Guazzelli and J. F. Morris,A physical introduction to suspension dynamics, Vol. 45 (Cambridge University Press, 2011)

    2011

  48. [48]

    A. D. Rey and M. M. Denn, Annual Review of Fluid Mechanics34, 233 (2002)

    2002

  49. [49]

    D. L. Koch and G. Subramanian, Annual Review of Fluid Mechanics43, 637 (2011)

    2011

  50. [50]

    J. S. Guasto, R. Rusconi, and R. Stocker, Annual Review of Fluid Mechanics44, 373 (2012)

    2012

  51. [51]

    Lauga, Annual Review of Fluid Mechanics48, 105 (2016)

    E. Lauga, Annual Review of Fluid Mechanics48, 105 (2016)

    2016

  52. [52]

    J. M. Ottino and S. Wiggins, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences362, 923 (2004)

    2004

  53. [53]

    C. A. Browne and S. S. Datta, Proceedings of the National Academy of Sciences121, e2320962121 (2024)

    2024

  54. [54]

    Sasmal, Journal of Non-Newtonian Fluid Mechanics , 105393 (2025)

    C. Sasmal, Journal of Non-Newtonian Fluid Mechanics , 105393 (2025)

    2025

  55. [55]

    M. E. Rosti, Fluid Dyn. Res.58, 021401 (2026)

    2026

  56. [56]

    M. E. Rosti, P. Perlekar, and D. Mitra, Science Advances9, eadd3831 (2023)

    2023

  57. [57]

    M. S. Abdelgawad, I. Cannon, and M. E. Rosti, Nature Physics19, 1059 (2023)

    2023

  58. [58]

    Sugiyama, S

    K. Sugiyama, S. Ii, S. Takeuchi, S. Takagi, and Y. Matsumoto, Journal of Computational Physics230, 596 (2011)

    2011

  59. [59]

    Zayernouri and M

    M. Zayernouri and M. Metzger, Physics of Fluids23(2011). 9 METHODS The numerical simulations are performed using the in-house solverFujin[49], which has been exten- sively validated including for viscoelastic fluids [12, 50, 51] (also seehttps://www.oist.jp/research/ research-units/cffu/fujin). The code uses a staggered, uniform Cartesian grid with a seco...

    2011

  60. [60]

    0 1 . 4 1 . 8 ⟨Cxx⟩x − 14 0 14 y/δ 0 tU0/δ 0 ≈ 70 tU0/δ 0 ≈ 73 tU0/δ 0 ≈ 76 tU0/δ 0 ≈ 78 tU0/δ 0 ≈ 81 tU0/δ 0 ≈ 84 0 1 2 ⟨Cyy⟩x − 14 0 14 FIG. S3: The spatial variation of⟨C xx⟩x versusy/δ 0 at various timestU 0/δ0 near the overturning time; the inset shows the unchanging⟨C yy ⟩x. The parameters are the same as the main paper, i.e.Re= 0.2,De= 28,β= 0.9. 3...

  61. [61]

    I. C. Christov, Mechanics Research Communications37, 717 (2010). 6

    2010

  62. [62]

    0000 0 . 0001 0 . 0002 ⟨Cxx⟩x +1 − 14 0 14 y/δ 0 − 2 0 2 4 ⟨Cxy⟩x × 10− 6 − 14 0 14 tU0/δ 0 ≈ 280 tU0/δ 0 ≈ 283 tU0/δ 0 ≈ 286 tU0/δ 0 ≈ 288 tU0/δ 0 ≈ 291 tU0/δ 0 ≈ 294 FIG. S7: Left : spatial variation of⟨C xx⟩x versusy/δ 0 at various times neartU 0/δ0 ∼293.5; the parameters are the same as the main paper, i.e.Re= 0.2,De= 28,β= 0.9. Right : same for⟨C xy⟩...