Subharmonic instability of large-scale wavy structures in two-dimensional channels
Pith reviewed 2026-05-09 23:24 UTC · model grok-4.3
The pith
At high Reynolds numbers, large-scale wavy structures in two-dimensional channels become unstable through a subharmonic torsional mode that grows and splits them into multiple wave trains.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
At Re = 200000 the large-scale wavy structures extracted from DNS via singular value decomposition support a subharmonic torsional instability mode with growth rate λ_r = 0.18. This mode produces a half-wavelength shift, deforms the primary waves, splits them into multiple wave trains, and drives opposite-phase evolution. The same structures remain linearly stable at Re = 3000, consistent with the laminar-like flow field observed in that regime. The subharmonic mode thereby offers a distinct pathway for turbulence generation in two-dimensional channels beyond the transitional regime near Re ≈ 10000.
What carries the argument
Floquet-based secondary instability analysis applied to singular-value-decomposition-extracted large-scale wavy structures, which identifies the subharmonic torsional mode responsible for the half-wavelength shift and positive growth.
If this is right
- The instability is absent below the transitional regime and appears only at sufficiently high Reynolds numbers such as 200000.
- The unstable mode deforms the wavy structures and splits them into multiple wave trains that evolve out of phase.
- This subharmonic process supplies a concrete mechanism for breaking down large-scale coherence and generating turbulence in two-dimensional channels.
- The quantified growth rate of the mode reaches λ_r = 0.18 at Re = 200000.
Where Pith is reading between the lines
- The same subharmonic torsional mechanism may operate in other two-dimensional systems that develop large-scale coherent structures through inverse cascades, such as geophysical or plasma flows.
- Controlled introduction of half-wavelength perturbations in high-Re channel simulations could directly measure the predicted growth rate and splitting behavior.
- The instability may couple with or be modulated by the small-scale fluctuations that the SVD filtering removes, suggesting a possible interaction between scales in the full turbulent state.
Load-bearing premise
The singular value decomposition accurately isolates the dominant large-scale wavy structures so that the subsequent Floquet analysis captures the true instability mechanism without significant contamination from unresolved small-scale fluctuations.
What would settle it
A high-Re direct numerical simulation that imposes a half-wavelength perturbation on the extracted wavy structures and observes no exponential growth at a rate near 0.18, or no subsequent splitting into out-of-phase wave trains, would falsify the instability claim.
Figures
read the original abstract
A particular interest on two-dimensional turbulence is the inverse energy cascade from small to large sales, which leads to an energy condensation accompanied by the formation of large-scale vortical structures. Indeed, such a phenomenon is observed in the two-dimensional channel (2DCH) with large Reynolds numbers, where prominent large-scale wavy structures play a central role in the momentum and energy transfer across the inhomogeneous wall-normal direction \citep{Falkovich2018}. Yet, the instability of these wavy structures remains poorly understood, and it is unknown whether they have the capacity to generate turbulence. To address this, we first conduct the direct numerical simulation (DNS) of Navier-Stokes equations for 2DCH, then extract the large-scale wavy structures through the singular value decomposition, and finally perform a Floquet-based secondary instability analysis. Two bulk Reynolds numbers are examined in particular, i.e. $Re = 3000$ and $Re = 200000$, which lie on opposite sides of the transitional regime near $Re \approx 10000$ and cover the previously reported simulation domain. At $Re = 3000$, the large-scale wavy structure is found to be linearly stable, consistent with the laminar-like DNS flow field. However, at $Re = 200000$, a subharmonic torsional mode is identified, which leads to a definite growth rate ($\lambda_r = 0.18$) for the wavy structures with a half wave-length shift. Temporal reconstruction shows that this unstable mode deforms and splits into multiple wave trains and evolves in the opposite phase. Compared to the TS (Tollmien-Schlichting) wave of laminar flow, the subharmonic mode found here offers a novel understanding for the generation of turbulence in larger Reynolds number two-dimensional channels.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript performs DNS of 2D channel flow at Re=3000 and Re=200000, extracts the dominant large-scale wavy structures via SVD, and subjects them to Floquet secondary instability analysis. At the lower Re the extracted structure is linearly stable, consistent with the laminar-like DNS; at the higher Re a subharmonic torsional mode is reported with growth rate λ_r=0.18. Temporal reconstruction of the unstable eigenmode shows deformation, splitting into multiple wave trains, and opposite-phase evolution, which the authors propose as a route to turbulence generation in high-Re 2D channels.
Significance. If the SVD mode can be shown to constitute a sufficiently clean base flow, the identification of a subharmonic torsional instability supplies a concrete, falsifiable mechanism linking the inverse-cascade condensate to the onset of small-scale turbulence. The pipeline (DNS → SVD → Floquet) is standard and reproducible in principle; the concrete growth-rate value and the half-wavelength shift provide clear targets for future verification.
major comments (2)
- [§3.2 and §4.1] §3.2 and §4.1: the claim that the leading SVD mode furnishes a clean, steady base flow for the Floquet operator is load-bearing for the reported λ_r=0.18. At Re=200000 the inverse cascade populates a broad range of scales; without reported energy spectra of the retained SVD mode, residual of the steady Navier-Stokes equations on that mode, or comparison against a filtered base flow, it is impossible to rule out contamination by unresolved fluctuations that would alter the linear operator and the extracted growth rate.
- [§4.2] §4.2, Floquet dispersion relation: the subharmonic character and the precise value λ_r=0.18 are presented without accompanying resolution or convergence tests (grid spacing, number of retained SVD modes, Floquet truncation). Because the low-Re case is stable, it does not probe the same scale-separation issue; an independent check at Re=200000 is required to establish that the growth rate is not an artifact of the numerical representation of the base flow.
minor comments (2)
- [Abstract] Abstract, line 3: 'large sales' is a typographical error for 'large scales'.
- [§2] The manuscript would benefit from a brief statement of the numerical scheme, time-stepping method, and domain size used in the DNS, even if only by reference to a prior publication.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback and for recognizing the potential significance of the subharmonic torsional instability as a mechanism linking the inverse-cascade condensate to turbulence onset. We address each major comment below and will revise the manuscript accordingly to strengthen the supporting evidence.
read point-by-point responses
-
Referee: [§3.2 and §4.1] the claim that the leading SVD mode furnishes a clean, steady base flow for the Floquet operator is load-bearing for the reported λ_r=0.18. At Re=200000 the inverse cascade populates a broad range of scales; without reported energy spectra of the retained SVD mode, residual of the steady Navier-Stokes equations on that mode, or comparison against a filtered base flow, it is impossible to rule out contamination by unresolved fluctuations that would alter the linear operator and the extracted growth rate.
Authors: We agree that explicit verification of the base-flow cleanliness is essential. The SVD isolates the energetically dominant structure by construction, but to rule out contamination we will add to the revised manuscript: the kinetic-energy spectrum of the leading SVD mode (demonstrating its dominance over higher modes), the L2 residual of the steady Navier-Stokes operator evaluated on this mode, and a direct comparison of the SVD mode against a spectrally filtered DNS snapshot. These diagnostics will quantify any residual fluctuations and confirm that they do not materially affect the extracted growth rate. revision: yes
-
Referee: [§4.2] §4.2, Floquet dispersion relation: the subharmonic character and the precise value λ_r=0.18 are presented without accompanying resolution or convergence tests (grid spacing, number of retained SVD modes, Floquet truncation). Because the low-Re case is stable, it does not probe the same scale-separation issue; an independent check at Re=200000 is required to establish that the growth rate is not an artifact of the numerical representation of the base flow.
Authors: We concur that dedicated convergence tests at the high-Re case are required. In the revised manuscript we will include a new subsection (or appendix) presenting systematic checks at Re=200000: (i) variation of DNS grid spacing while keeping the SVD base flow fixed, (ii) sensitivity to the number of retained SVD modes used to construct the base flow, and (iii) convergence of the Floquet eigenvalue with respect to the number of retained Fourier modes in the secondary-instability expansion. These tests will demonstrate that both the subharmonic character and the value λ_r=0.18 remain robust. revision: yes
Circularity Check
No circularity: DNS-SVD-Floquet chain is independent of its outputs
full rationale
The derivation proceeds as DNS of the 2D channel flow at fixed Re, followed by SVD extraction of the dominant coherent structure to serve as a base flow, followed by a standard Floquet secondary-instability eigenvalue problem whose growth rate is computed from that base flow. The reported λ_r = 0.18 is an eigenvalue of the linearized operator on the extracted field; it is not fitted to itself, not defined by the instability result, and not justified by any self-citation chain. No step renames a known empirical pattern or imports a uniqueness theorem from the authors' prior work. The procedure is therefore self-contained against external benchmarks and receives the default non-circularity finding.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math The flow obeys the incompressible two-dimensional Navier-Stokes equations.
- domain assumption Singular value decomposition isolates the dominant large-scale wavy structures.
Reference graph
Works this paper leans on
-
[1]
Exact third-order structure functions for two-dimensional turbulence , journal =
Xie, Jin-Han and B. Exact third-order structure functions for two-dimensional turbulence , journal =. 2018 , doi =
work page 2018
-
[2]
Journal of Fluid Mechanics , volume =
Xie, Jin-Han and Huang, Shi-Di , title =. Journal of Fluid Mechanics , volume =. 2022 , doi =
work page 2022
-
[3]
Physical Review Letters , volume =
Zhu, Hang-Yu and Xie, Jin-Han and Xia, Ke-Qing , title =. Physical Review Letters , volume =. 2023 , doi =
work page 2023
- [4]
-
[5]
Wu, X. and Dong, M. , title =. Journal of Fluid Mechanics , volume =. 2016 , doi =
work page 2016
-
[6]
K. Effects of streamwise-elongated and spanwise-periodic surface roughness elements on boundary-layer instability , journal =. 2020 , doi =
work page 2020
-
[7]
Dong, M. and Liu, Y. and Wu, X. , title =. Journal of Fluid Mechanics , volume =. 2020 , doi =
work page 2020
-
[8]
Dong, M. and Zhao, L. , title =. Journal of Fluid Mechanics , volume =. 2021 , doi =
work page 2021
-
[9]
Dong, M. and Zhang, M. , title =. Journal of Fluid Mechanics , volume =. 2022 , doi =
work page 2022
-
[10]
Dong, M. and Wan, D. , title =. Journal of Fluid Mechanics , volume =. 2025 , doi =
work page 2025
-
[11]
Chen, X. and Huang, G. L. and Lee, C. B. , title =. Journal of Fluid Mechanics , volume =. 2019 , doi =
work page 2019
-
[12]
Chen, X. and Dong, S. and Tu, G. and Yuan, X. and Chen, J. , title =. Journal of Fluid Mechanics , volume =. 2022 , doi =
work page 2022
-
[13]
Kellay, H. and Wu, X.-l. and Goldburg, W. I. , title =. Phys. Rev. Lett. , volume =. 1995 , doi =
work page 1995
-
[14]
Kellay, H. and Goldburg, W. , title =. Rep. Prog. Phys. , volume =. 2002 , doi =
work page 2002
-
[15]
Tran, T. and Chakraborty, P. and Guttenberg, N. and Prescott, A. and Kellay, H. and Goldburg, W. and Goldenfeld, N. and Gioia, G. , title =. Nat. Phys. , volume =. 2010 , doi =
work page 2010
-
[16]
Kellay, H. and Tran, T. and Goldburg, W. and Goldenfeld, N. and Gioia, G. and Chakraborty, P. , title =. Phys. Rev. Lett. , volume =. 2012 , doi =
work page 2012
-
[17]
Vilquin, A. and Jagielka, J. and Djambov, S. and Herouard, H. and Fischer, P. and Bruneau, C.-H. and Chakraborty, P. and Gioia, G. and Kellay, H. , title =. Sci. Adv. , volume =. 2021 , doi =
work page 2021
-
[18]
L'vov, V. S. and Procaccia, I. and Rudenko, O. , title =. Phys. Rev. E , volume =. 2009 , doi =
work page 2009
-
[19]
Samanta, D. and Ingremeau, F. and Cerbus, R. and Tran, T. and Goldburg, W. I. and Chakraborty, P. and Kellay, H. , title =. Phys. Rev. Lett. , volume =. 2014 , doi =
work page 2014
-
[20]
Markeviciute, V. K. and Kerswell, R. R. , title =. J. Fluid Mech. , volume =. 2021 , doi =
work page 2021
-
[21]
Falkovich, G. and Vladimirova, N. , title =. Phys. Rev. Lett. , volume =. 2018 , doi =
work page 2018
-
[22]
Macroscopic effects of the spectral structure in turbulent flows , author=. Nature Physics , volume=. 2010 , publisher=
work page 2010
-
[23]
Physical review letters , volume=
Testing a missing spectral link in turbulence , author=. Physical review letters , volume=. 2012 , publisher=
work page 2012
-
[24]
Inertial ranges in two-dimensional turbulence , author=
-
[25]
Computation of the energy spectrum in homogeneous two-dimensional turbulence , author=. Physics of Fluids , volume=. 1969 , publisher=
work page 1969
-
[26]
International Journal of Computational Fluid Dynamics , volume=
A low-communication-overhead parallel DNS method for the 3D incompressible wall turbulence , author=. International Journal of Computational Fluid Dynamics , volume=. 2021 , publisher=
work page 2021
-
[27]
Stability and transition in shear flows , author=. 2012 , publisher=
work page 2012
- [28]
-
[29]
Schlichting, H. , journal =. Zur Enstehung der Turbulenz bei der Plattenströmung , url =
-
[30]
Landau, Lev Davidovich. On the Problem of Turbulence. Compt. Rend. Acad. Sci. URSS. 1944. doi:10.1016/b978-0-08-010586-4.50057-2
-
[31]
Herbert, Thorwald , title =. Physics of Fluids , volume =. 1983 , month =
work page 1983
-
[32]
Stuart, J. T. , title=. Journal of Fluid Mechanics , volume=. 1958 , doi=
work page 1958
-
[33]
The basic behaviour in plane Poiseuille flow , volume=
On the non-linear mechanics of wave disturbances in stable and unstable parallel flows Part 1. The basic behaviour in plane Poiseuille flow , volume=. Journal of Fluid Mechanics , author=. 1960 , pages=. doi:10.1017/S002211206000116X , number=
-
[34]
The development of a solution for plane Poiseuille flow and for plane Couette flow , volume=
On the non-linear mechanics of wave disturbances in stable and unstable parallel flows Part 2. The development of a solution for plane Poiseuille flow and for plane Couette flow , volume=. Journal of Fluid Mechanics , author=. 1960 , pages=. doi:10.1017/S0022112060001171 , number=
-
[35]
Nagata, M. , year=. Three-dimensional finite-amplitude solutions in plane Couette flow: bifurcation from infinity , volume=. doi:10.1017/S0022112090000829 , journal=
-
[36]
Nagata, M. , year=. Bifurcations in Couette flow between almost corotating cylinders , volume=. doi:10.1017/S0022112086000605 , journal=
-
[37]
Journal of Fluid Mechanics , author=
Accurate solution of the Orr–Sommerfeld stability equation , volume=. Journal of Fluid Mechanics , author=. 1971 , pages=. doi:10.1017/S0022112071002842 , number=
-
[38]
In AGARD Laminar-Turbulent Transition 10 p (SEE N78-14316 05-34 , year=
Finite amplitude stability of plane parallel flows , author=. In AGARD Laminar-Turbulent Transition 10 p (SEE N78-14316 05-34 , year=
-
[39]
Ehrenstein, U. and Koch, W. , year=. Three-dimensional wavelike equilibrium states in plane Poiseuille flow , volume=. doi:10.1017/S0022112091002653 , journal=
-
[40]
Coherent structures in wall-bounded turbulence , volume=
Jiménez, Javier , year=. Coherent structures in wall-bounded turbulence , volume=. doi:10.1017/jfm.2018.144 , journal=
-
[41]
Annual Review of Fluid Mechanics , volume=
Exact coherent states and the nonlinear dynamics of wall-bounded turbulent flows , author=. Annual Review of Fluid Mechanics , volume=. 2021 , publisher=
work page 2021
-
[42]
International journal for numerical methods in fluids , volume=
Onset of transition in boundary layers , author=. International journal for numerical methods in fluids , volume=. 1988 , publisher=
work page 1988
-
[43]
Annual review of fluid mechanics , volume=
Secondary instability of boundary layers , author=. Annual review of fluid mechanics , volume=
-
[44]
H. B. Squire , journal =. On the Stability for Three-Dimensional Disturbances of Viscous Fluid Flow between Parallel Walls , urldate =
-
[45]
Journal of Fluid Mechanics , volume=
Secondary instability of wall-bounded shear flows , author=. Journal of Fluid Mechanics , volume=. 1983 , publisher=
work page 1983
-
[46]
Transition to turbulence in two-dimensional Poiseuille flow , volume=
Jiménez, Javier , year=. Transition to turbulence in two-dimensional Poiseuille flow , volume=. doi:10.1017/S0022112090001008 , journal=
-
[47]
Annales scientifiques de l'École Normale Supérieure, Serie 2 , volume=
Sur les équations différentielles linéaires à coefficients périodiques , author=. Annales scientifiques de l'École Normale Supérieure, Serie 2 , volume=
-
[48]
Applied Mathematics and Mechanics (English Edition) , volume =
Han, Yufeng and Liu, Jianxin and Luo, Jisheng , title =. Applied Mathematics and Mechanics (English Edition) , volume =. 2018 , doi =
work page 2018
-
[49]
Secondary subharmonic instability of hypersonic boundary layer in thermochemical equilibrium over a flat plate , author=. Physics of Fluids , volume=. 2021 , publisher=
work page 2021
-
[50]
The approximation of one matrix by another of lower rank , author=. Psychometrika , volume=. 1936 , publisher=
work page 1936
- [51]
-
[52]
Journal of Fluid Mechanics , volume=
Transition to turbulence in plane Poiseuille and plane Couette flow , author=. Journal of Fluid Mechanics , volume=. 1980 , publisher=
work page 1980
-
[53]
Reports on Progress in Physics , volume =
Two-dimensional turbulence , author =. Reports on Progress in Physics , volume =
-
[54]
Annual review of fluid mechanics , volume=
Two-dimensional turbulence , author=. Annual review of fluid mechanics , volume=. 2012 , publisher=
work page 2012
-
[55]
Journal of Fluid Mechanics , volume=
Bifurcations of two-dimensional channel flows , author=. Journal of Fluid Mechanics , volume=. 1986 , publisher=
work page 1986
-
[56]
Physical Review Fluids , volume=
Unbounded two-dimensional wall turbulence induced by inverse cascade , author=. Physical Review Fluids , volume=. 2024 , publisher=
work page 2024
-
[57]
Annual Review of Fluid Mechanics , volume=
Nonlinear nonmodal stability theory , author=. Annual Review of Fluid Mechanics , volume=. 2018 , publisher=
work page 2018
-
[58]
Journal of Fluid Mechanics , volume=
The structure of the vortices in freely decaying two-dimensional turbulence , author=. Journal of Fluid Mechanics , volume=. 1996 , publisher=
work page 1996
-
[59]
Journal of Fluid Mechanics , volume=
Dipoles and streams in two-dimensional turbulence , author=. Journal of Fluid Mechanics , volume=. 2020 , publisher=
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.