REVIEW 2 major objections 5 minor 34 references
Tollmien-Schlichting waves near neutral stable curve
T0 review · 2 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read This paper proves that Tollmien-Schlichting waves exist throughout a neighborhood of the neutral stability curve, with stable waves below, growing waves inside, and exact neutral modes at the two branch crossings.
desk verdict The neutral-curve construction is real progress, and the one real gap is a missing quantitative lemma on the Tietjens function on the lower branch, which looks fixable. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the modified Rayleigh operator Ray_ĉ[φ] = (u−ĉ)(∂$_Y^{2}$−$α^{2}$)φ − u''φ with ĉ = c + i|ε| $α^{{-3/2}}$, a regularized version of the singular Rayleigh operator used on the neutral curve. The imaginary shift c0 = |ε| $α^{{-3/2}}$ encodes the idea that viscous diffusion regularizes the critical layer and effectively pushes the eigenvalue left in the complex plane. Around this operator the paper sets up a modified Rayleigh-Airy iteration that builds two independent solutions of the Orr-Sommerfeld equation, the slow mode φ_s and fast mode φ_f, and then solves the dispersion relation φ_s(0)/∂_Y φ_s(0) = φ_f(0)/∂_Y φ_f(0). This dispersion relation, valid for ci ≤ 0, is what produces the neutral crossings.
What would settle it
A direct numerical solution of the Orr-Sommerfeld eigenvalue problem for a profile satisfying (1.8), at small viscosity, that computes Im(c(α)) through the whole range ($Aν^{{1/8}}$, $Bν^{{1/12}}$) would falsify the theorem if it did not show two sign changes near those scalings or if the interior growth rate α Im(c) did not scale like $ν^{{1/4}}$.
Extended reading notes
Core claim
The central result, Theorem 1.1, is that for shear profiles satisfying the structural assumptions (1.8), for every sufficiently small viscosity ν and every α in ($Aν^{{1/8}}$, $Bν^{{1/12}}$), there is a complex phase speed c(α) and a $W^{{2,∞}}$ solution of the linearized Navier-Stokes equations in the Tollmien-Schlichting form. The imaginary part of c changes sign twice: it is negative at the lower-branch edge, zero at an intermediate value A_c $ν^{{1/8}}$, positive in the interior, zero again at B_c $ν^{{1/12}}$, and negative above the upper branch; in the interior the growth rate satisfies α Im(c) ∼ $ν^{{1/4}}$. The paper interprets the two zeros as the lower and upper branches of the neutral stability curve, and identifies the lower branch with the competition between viscous diffusion and sublayer instability and the upper branch with the stabilizing curvature u''(0)<0 of the background profile.
Load-bearing premise
Everything rests on the claim that near the neutral curve the Orr-Sommerfeld operator is accurately approximated, in the upper and main decks, by the modified Rayleigh operator with the specific imaginary shift |ε|$α^{{-3/2}}$, an approximation justified heuristically and whose rigorous estimates are quoted from an earlier paper rather than proved here.
Editorial extensions
If this is right
- The lower and upper branches of the neutral curve are located at α ∼ ν^{1/8} and α ∼ ν^{1/12}, and on each branch there is a neutral T-S wave with Im(c)=0.
- For wave numbers strictly inside the two branches the T-S waves are unstable, with growth rate α Im(c) of order ν^{1/4}; outside the branches they are damped.
- The construction admits shear profiles without analyticity; the only structural conditions are the ones in (1.8), and u''(0)<0 is not needed near the lower branch, which allows the Blasius profile there.
- The transition of linear stability near the lower branch is governed by diffusion versus sublayer instability, while near the upper branch it is governed by the curvature of the background flow.
- The neutral stable curve separates stable from unstable Tollmien-Schlichting modes exactly at the predicted scalings, matching the physical neutral curve of Tollmien, Schlichting and Lin.
Reading between the lines
- If the modified-Rayleigh shift |ε|α^{-3/2} is the right regularization generally, the same device should locate neutral curves for other high-Reynolds shear flows, including channel and pipe flows, where the analogous singular Rayleigh problem appears.
- The paper notes that for the Blasius profile the physical upper branch is α_up ∼ ν^{1/20}; extending the transition argument to that branch would require additional estimates beyond the present theorem, and the method developed here is a plausible starting point.
- An immediate testable extension is to compare the predicted dispersion relation, expressed through the Tietjens-type function F(−κη(0)), with direct numerical solution of the full Orr-Sommerfeld equation at finite but small ν; the crossing locations and the ν^{1/4} growth law are quantitative enough to check numerically.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies linear stability of boundary-layer shear flows over a flat plate. After rescaling, Tollmien-Schlichting waves correspond to solutions of the Orr-Sommerfeld equation (1.7). The authors construct the slow and fast modes of the homogeneous Orr-Sommerfeld operator via a modified Rayleigh-Airy iteration, with the key new ingredient being the replacement of the singular Rayleigh operator by Ray_{c+i|ε|α^{-3/2}}. They then solve the dispersion relation (5.3) and obtain Theorem 1.1: for every α in (Aν^{1/8}, Bν^{1/12}) there is a phase speed c(α) and a W^{2,∞} T-S wave, with Im(c) crossing from negative to positive at lower and upper branch points and with α Im(c) ∼ ν^{1/4} in the interior. Much of the technical machinery, especially the Rayleigh and Airy estimates in Sections 2 and 3, is quoted from the authors' earlier paper [27].
Significance. If the theorem is valid, it gives the first rigorous confirmation of the neutral stable curve for this problem, complementing the construction of unstable T-S waves by Grenier, Guo, and Nguyen. The result is substantive: it identifies the asymptotic locations α ∼ ν^{1/8} and α ∼ ν^{1/12}, gives the growth-rate scaling α Im(c) ∼ ν^{1/4}, and removes the analyticity assumption of [15] through a modified Langer transformation. The derivation is not fitted: c0 is a fixed regularization parameter and the sign changes come from the dispersion relation, so the theorem is explicit and falsifiable. The proof is long and depends heavily on estimates quoted from the authors' earlier paper [27], which limits the self-containedness of the manuscript.
major comments (2)
- [Section 5, Proposition 5.1, Case 1, Eq. (5.9)] The proof of the lower-branch crossing is incomplete. After reducing to F(−κη(0)) + O(ν^{1/16}) = u′(0)c/α, the text invokes “the property of F(z) (see Fig 3.2 in [26])” to find z0 ∈ [2, 2.5] with Im F(z0) = 0 and Im F(z) < 0 for z ∈ [2, z0). It then asserts that a perturbation z1 with |Im z1| ≤ Cν^{1/16} yields Im(F(z0+z1)+O(ν^{1/16})) < −(1/8)ν^{1/16}. This requires a quantitative transversality statement, e.g. a lower bound on |Im F′(z0)| or at least a known finite order of vanishing at z0. A textbook figure gives the qualitative sign only; if the derivative at z0 is zero or too small, the O(ν^{1/16}) error can flip the sign and no A0 with Im(c) < 0 follows. Since Theorem 1.1's first bullet depends exactly on this point, Proposition 5.1 needs a precise lemma giving the definition of F and a quantitative lower bound on |Im F| on a ν-independent interval to the left of z0.
- [Sections 2–3 and Theorem 4.1] The construction of the slow and fast modes, and hence the dispersion relation, rests on a large body of estimates quoted without proof from [27]: Proposition 2.1, Lemmas 2.2 and 2.4, Lemmas 3.1–3.4, Propositions 3.5–3.7, Theorems 3.8–3.10, and Lemma A.5. Because [27] is co-authored by two of the present authors, this is not circular in the sense of fitting a target quantity, but it makes the central claim difficult to verify from the manuscript alone. The key assumption behind (1.12) — that Ray_{c+i|ε|α^{-3/2}} approximates the Orr-Sommerfeld operator to the order required by the iteration — is justified only through these quoted estimates. The paper should either state the needed estimates as explicit hypotheses or provide proofs; at minimum it should explicitly verify that the hypotheses of [27], in particular the c0 range |ε|^{2/3} ≪ c0 ≪ |ε|^{1/3}, are satisfied with the choice c0 = |ε|α^{-3/2} throughout the interval (Aν^{1/8}, Bν^{1/12}).
minor comments (5)
- [Throughout] There are several spelling errors: “Naiver” before (1.3), “Nyugen” in the discussion of [15], “exsit” in Theorem 1.1, and “adn” in reference [3].
- [Proposition 5.2] In the displayed partial derivatives, the third line repeats ∂crFi; it should presumably be ∂ciFr, and the fourth line should be ∂ciFi.
- [Eq. (5.9)] The Tietjens function F(z) is not defined in the paper; the proof should provide its formula or a precise reference with a quantitative statement, rather than only citing a figure.
- [Section 4.2] The notation H1, H2 is used both for the parameter sets in (2.5) and for the source-term functions in (4.12), which is confusing but harmless.
- [Theorem 1.1 and Proposition 5.1] The strict inequality A0 > A appears in Theorem 1.1 but only A0 is mentioned in Proposition 5.1; the statements should be made consistent.
Circularity Check
No significant circularity: c(α) and neutral crossings are derived from the dispersion relation; heavy use of [27] is independent proof-based support, and the Tietjens-function step is an external input with an unproved quantitative margin.
full rationale
The paper's central claims are obtained by solving the dispersion relation (5.3)-(5.4), which is derived from the Orr-Sommerfeld equation via slow/fast mode constructions. No target quantity (Im(c)=0, A_c, B_c, or the ν^{1/4} growth rate) is fitted: the signs of Im(c) are read off asymptotic expansions of Airy-function ratios and the explicit estimates (5.8) and (5.10), and the interior growth α Im(c)∼ν^{1/4} follows from the same relation. The modified Rayleigh operator with ĉ=c+i|ε|α^{-3/2} is chosen within an explicitly bounded interval of admissible regularizations; c0 is a technical parameter, not a fitted constant, and the leading-order sign changes are independent of it. The many lemmas imported from [27] are quoted with their assumptions (e.g., H2, |ci|≤C|ε|^{1/3}); these are previously published proof-based estimates for Rayleigh and Airy equations and do not include the neutral-curve existence as an assumption, so citing them is legitimate independent support even though two of the present authors are co-authors. The one notable concern is in Proposition 5.1, Case 1, where the lower-branch sign change is made to rest on the property of the Tietjens function F taken from Fig 3.2 of Lin [26]: the paper asserts a perturbation z1 giving |Im(F(z0+z1)+O(ν^{1/16}))|>ν^{1/16}/8 without proving the required quantitative slope/transversality of Im F at z0. That is a rigor gap in the verification of the lower-branch crossing, not a circular reduction: the Tietjens property is an external classical input, not the theorem's conclusion, and no equation of the paper is equivalent by construction to the asserted Im(c)=0. Hence the derivation chain is not circular, and the appropriate score is 0.
Assumptions & free parameters
free parameters (1)
- c0 = |ε| α^{-3/2} (imaginary shift in modified Rayleigh operator) =
|ε| α^{-3/2}
assumptions (4)
- domain assumption Profile u(Y) satisfies (1.8): u(0)=0, lim u=1, u'(0)>0, u''(0)<0, and exponential decay of (u-1) and derivatives.
- domain assumption Wavenumber and phase-speed regime: α ∈ (Aν^{1/8}, Bν^{1/12}), c_r ∼ α, |c_i| ≪ α (set H2/H3).
- standard math Tietjens function F(z) has a zero of Im F at some z0 ∈ [2,2.5] with Im F<0 on [2,z0), as read from Fig 3.2 of [26].
- ad hoc to paper The modified Rayleigh operator -Ray_{c+i|ε|α^{-3/2}} approximates the Orr-Sommerfeld operator in the upper and main deck on the neutral curve to the order required by the iteration.
Cite this review
Pith. "Pith review of Tollmien-Schlichting waves near neutral stable curve." pith.science (2026). https://pith.science/paper/IY32DIDM
@misc{pith2026250202258,
author = {Pith},
title = {Pith review of: Tollmien-Schlichting waves near neutral stable curve},
year = {2026},
howpublished = {\url{https://pith.science/paper/IY32DIDM}},
note = {Machine review of arXiv:2502.02258}
}
abstract
In this paper, we study the linear stability of boundary layer flows over a flat plate. Tollmien, Schlichting, Lin et al. found that there exists a neutral curve, which consists of two branches: lower branch $\alpha_{low}(Re)$ and upper branch $\alpha_{up}(Re)$. Here, $\alpha$ is the wave number and $Re$ is the Reynolds number. For any $\alpha\in(\alpha_{low},\alpha_{up})$, there exist unstable modes known as Tollmien-Schlichting (T-S) waves to the linearized Navier-Stokes system. These waves play a key role during the early stage of boundary layer transition. In a breakthrough work (Duke math Jour, 165(2016)), Grenier, Guo, and Nguyen provided a rigorous construction of the unstable T-S waves. In this paper, we confirm the existence of the neutral stable curve. To achieve this, we develop a more delicate method for solving the Orr-Sommerfeld equation by borrowing some ideas from the triple-deck theory. This approach allows us to construct the T-S waves in a neighborhood of the neutral curve.
Reference graph
Works this paper leans on
-
[15]
E. Grenier, Y. Guo and T. Nguyen, Spectral instability of characteristic boundary layer flow s, Duke Math. J., 165(2016), 3085-3146
work page 2016
-
[27]
N. Masmoudi, Y. Wang, D. Wu and Z. Zhang, Tollmien-Schlichting waves in the subsonic regime , Proc. Lond. Math. Soc. (3), 128 (2024), Paper No. e12588, 112 pp
work page 2024
-
[26]
C. C. Lin, The theory of hydrodynamic stability, Cambri dge, at the University Prees, 1955
work page 1955
-
[1]
Y. Almog and B. Helffer, On the stability of laminar flows between plates , Arch. Ration. Mech. Anal., 241(2021), 1281-1401
work page 2021
-
[2]
On the stability of symmetric flows in a two-dimensional channel
Y. Almog and B. Helffer, On the stability of symmetric flows in a two-dimensional chan nel, arXiv:2212.12827. 60 QI CHEN, DI WU, AND ZHIFEI ZHANG
-
[3]
D. Bian adn E. Grenier, Onset of nonlinear instabilities in monotonic viscous boun dary layers , SIAM J. Math. Anal., 56 (2024), 3703–3719
work page 2024
-
[4]
Asymptotic behaviour of solutions of linearized Navier Stokes equations in the long waves regime
D. Bian and E. Grenier, Asymptotic behaviour of solutions of linearized Navier Sto kes equations in the long waves regime, arXiv:2312.16938
-
[5]
Instability of shear layers and Prandtl's boundary layers
D. Bian and E. Grenier, Instability of shear layers and Prandtl’s boundary layers , arXiv:2401.15679
Show all 34 references
-
[6]
Q. Chen, D. Wei and Z. Zhang, Linear stability of pipe Poiseuille flow at high Reynolds num ber, Comm. Pure Appl. Math., Comm. Pure Appl. Math., 76 (2023), 1868–1964
2023
-
[7]
Q. Chen, D. Wei and Z. Zhang, Linear inviscid damping and enhanced dissipation for monot one shear flows , Comm. Math. Phys., 400 (2023), 215–276
2023
-
[8]
Q. Chen, D. Wu and Z. Zhang, On the L∞ stability of Prandtl expansions in Gevrey class , Sci. China Math., 65(2022), 2521-2562
2022
-
[9]
Q. Chen, D. Wu and Z. Zhang, On the stability of shear flows of Prandtl typefor the steady N avier-Stokes equations, Sci. China Math., 66(2023), 679-722
2023
-
[10]
P. G. Drain and W. H. Reid, Hydrodynamic Stability, Camb ridge Monogr. Mech. Appl. Math., Cambridge Univ.Press, Cambridge, 1981
1981
-
[11]
Gao and L
C. Gao and L. Zhang, On the steady Prandtl boundary layer expansions , Sci. China Math., 66(2023), 1993- 2020
2023
-
[12]
G´ erard-Varet and Y
D. G´ erard-Varet and Y. Maekawa, Sobolev stability of Prandtl expansions for the steady Navi er-Stokes equa- tions, Arch. Ration. Mech. Anal., 223(2019), 1319-1382
2019
-
[13]
G´ erard-Varet, Y
D. G´ erard-Varet, Y. Maekawa and N. Masmoudi, Gevrey stability of Prandtl expansions for 2-dimensional Navier-Stokes flows , Duke Math. J., 167(2018), 2531-2631
2018
-
[14]
G´ erard-Varet, Y
D. G´ erard-Varet, Y. Maekawa and N. Masmoudi, Optimal Prandtl expansion around concave boundary layer , Anal. PDE, 17 (2024), 3125–3187
2024
-
[16]
Grenier, Y
E. Grenier, Y. Guo and T. Nguyen, Spectral instability of general symmetric shear flows in a tw o-dimensional channel, Adv. Math., 292 (2016), 52-110
2016
-
[17]
Grenier and T, Nguyen, Sublayer of Prandtl boundary layers , Arch
E. Grenier and T, Nguyen, Sublayer of Prandtl boundary layers , Arch. Ration. Mech. Anal., 229 (2018), 1139–1151
2018
-
[18]
Grenier and T
E. Grenier and T. Nguyen, L∞ instability of Prandtl layers , Ann. PDE, 2019, 5: Paper No. 18, 36 pp
2019
-
[19]
Grenier and T
E. Grenier and T. Nguyen, On nonlinear instability of Prandtl’s boundary layers: the case of Rayleigh’s stable shear flows , J. Math. Pures Appl. (9), 184 (2024), 71–90
2024
-
[20]
Guo and S
Y. Guo and S. Iyer, Regularity and expansion for steady Prandtl equations , Comm. Math. Phys., 382 (2021), 1403-1447
2021
-
[21]
Guo and S
Y. Guo and S. Iyer, Validity of steady Prandtl layer expansions , Comm. Pure Appl. Math., 76 (2023), 3150–3232
2023
-
[22]
Guo and T
Y. Guo and T. Nguyen, Prandtl boundary layer expansions of steady Navier-Stokes flows over a moving plate , Ann. PDE, 2017, 3: Paper No. 10, 58 pp
2017
-
[23]
Y. Guo, Y. Wang and Z. Zhang, Dynamic stability for steady Prandtl solutions , Ann. PDE, 2023, 9: Paper No. 16, 33 pp
2023
-
[24]
Iyer, On global-in-x stability of Blasius profiles , Arch
S. Iyer, On global-in-x stability of Blasius profiles , Arch. Ration. Mech. Anal., 237 (2020), 951–998
2020
-
[25]
Iyer and N
S. Iyer and N. Masmoudi, Global-in-x stability of steady Prandtl expansions for 2D N avier-Stokes flow , arXiv:2008.12347
2008 arXiv
-
[28]
Masmoudi, Y
N. Masmoudi, Y. Wang, D. Wu and Z. Zhang, Mack modes in supersonic boundary layer , arXiv:2405.04853
-
[29]
Olver, D
F. Olver, D. Lozier, R. Boisvert and C. Clark, NIST handb ook of mathematical functions, Cambridge Uni- versity Press, New York, 2010
2010
-
[30]
V. A. Romanov, Stability of plane-parallel Couette flow , Funkcional. Anal. i Priloˇ zen, 7 (1973), 62-73
1973
-
[31]
Schlichting, Boundary layer theory, 4th ed., McGraw -Hill Series in Mechanical Engineering, McGraw-Hill, New York, 1960
H. Schlichting, Boundary layer theory, 4th ed., McGraw -Hill Series in Mechanical Engineering, McGraw-Hill, New York, 1960
1960
-
[32]
Wang and Z
Y. Wang and Z. Zhang, Global C ∞ regularity of the steady Prandtl equation with favorable pr essure gradient, Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire, 38 (2021), 1989–2004
2021
-
[33]
Wang and Z
Y. Wang and Z. Zhang, Asymptotic behavior of the steady Prandtl equation , Math. Ann., 387 (2023), 1289–1331. TOLLMIEN-SCHLICHTING W A VES NEAR NEUTRAL STABLE CUR VE 61
2023
-
[34]
Yang and Z
T. Yang and Z. Zhang, Linear instability analysis on compressible Navier-Stoke s equations with strong bound- ary layer , Arch. Ration. Mech. Anal., 247 (2023), Paper No. 83, 53 pp. (Q. Chen) School of Mathematical Sciences, Zhejiang University, Han gzhou, 310058, P. R. China ...
2023
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.