Pith. sign in

REVIEW 2 major objections 4 minor 36 references

Unstable mode around the 3D boundary layer flow

T0 review · 2 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper constructs unstable modes for the linearized Navier-Stokes equations around three-dimensional boundary-layer shear flows in the small-viscosity limit, with growth rates of order exp(t/sqrt(nu)).

desk verdict A genuinely new 3D boundary-layer instability mechanism with a clean core proof, but the advertised 'generic' profiles and 'any δ>0' claim outrun what the stated hypotheses actually support. read the letter →

arxiv 2509.06089 v1 pith:4PM2JJHY submitted 2025-09-07 math.AP physics.flu-dyn

Cheng-Jie Liu , Mengjun Ma , Di Wu , Zhu Zhang This is my paper · ORCID
classification math.APphysics.flu-dyn MSC 35Q3076D1076E05
keywords 3DboundarylayerNavier-StokeslinearstabilityTollmien-SchlichtingwavesRayleighequationOrr-Sommerfeldzero-viscositylimitPrandtlexpansionspectralinstability
topics Navier-Stokes Regularity
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

When the streamwise and spanwise velocity profiles are linearly independent near the wall, the authors prove there is a growing mode with rate exp(t/sqrt(nu)), faster than the classical two-dimensional Tollmien-Schlichting instability. The instability arises because the three-dimensional perturbation tilts the base flow into a modified profile U=(sigma/alpha)u_s+(beta/alpha)v_s that is positive, increasing, and convex at the boundary, which destabilizes the Rayleigh equation. The same mode persists in the viscous Orr-Sommerfeld equation through a Rayleigh-Airy iteration, because the viscous boundary layer only shifts the eigenvalue by a small amount. If correct, this means generic high-Reynolds-number three-dimensional boundary layers are analytically unstable, and Prandtl expansion validity in three dimensions cannot be expected without analytic regularity.

What carries the argument

The central object is the modified shear flow U(Z)=lambda_1 u_s(Z)+lambda_2 v_s(Z), with (lambda_1,lambda_2) chosen as in Proposition 2.1 so that, under the independence condition (2.1) and the quantitative bound (2.2), U is positive for Z>0, increasing and convex at the wall. This U acts as the effective advection speed in the Rayleigh and Orr-Sommerfeld equations for the vertical velocity. The proof machinery is the second-order asymptotic expansion of the Rayleigh solution phi_Ray(0;c) (Lemmas 3.4 and 3.7), which yields the unstable inviscid eigenvalue c_Ray approx alpha U_inf^2/U'(0) + i alpha^2 U''(0)U_inf^4 pi/U'(0)^4, followed by a Rayleigh-Airy iteration that constructs the Orr-Somme

What would settle it

Solve the linearized eigenvalue problem (1.4) numerically for the exponential profiles u_s=1-e^{-Z}, v_s=v_infty(1-e^{-kZ}) with k>1 and k|v_infty| large, at nu=10^{-6}. Theorem 1.1 predicts an eigenvalue with c_i ~ alpha^2 U''(0)U_inf^4 pi/U'(0)^4 and growth rate ~nu^{-1/2}; if no unstable eigenvalue appears, or if the growth rate scales differently in nu, the central claim fails. A second check: test profiles violating (2.2); the theorem makes no prediction there, so finding instability would not falsify it but would show the hypotheses are not sharp.

Watch

Extended reading notes

Core claim

Theorem 1.1 states: for any small viscosity 0<nu<<1 and a class of spanwise profiles v_s satisfying the structural conditions (2.1)-(2.2), the linearized Navier-Stokes system around the shear flow (u_s(z/sqrt(nu)), v_s(z/sqrt(nu)), 0) admits a nontrivial solution of the form e^{-i alpha c nu^{-1/2} t} e^{i nu^{-1/2}(sigma x+beta y)} (tilde u, tilde v, tilde w)(z/sqrt(nu)), with amplitudes in H^1, wave speed c=c_r+i c_i with c_r>0, c_i>0, and growth |(u,v,w)| ~ e^{C nu^{-1/2} t}. The eigenvalue is found first for the inviscid Rayleigh equation, where U''(0)>0 supplies a positive imaginary part, and then shown to persist for the Orr-Sommerfeld equation: the viscous boundary layer only shifts t

Load-bearing premise

The construction needs the quantitative structural inequality (2.2), which forces the modified profile U=lambda_1 u_s+lambda_2 v_s to be positive, increasing, and convex at the wall; this is not shown for generic profiles, and for the explicit exponential family it is only verified under a large-spanwise-amplitude condition, so the word 'generic' in the abstract overreaches the hypotheses.

Editorial extensions

If this is right

  • The 3D linearized Navier-Stokes boundary-layer problem has unstable modes with growth rate of order nu^{-1/2}, faster than the two-dimensional Tollmien-Schlichting growth rate.
  • The instability appears already at the Euler level (inviscid Rayleigh equation) and persists for small viscosity, because the viscous boundary layer only shifts the eigenvalue by O(nu^{1/4-}).
  • For any nonzero spanwise amplitude delta, no matter how small, instability occurs; at delta=0 the flow is stable in Gevrey-3/2 for concave profiles, so delta=0 is a bifurcation point.
  • The growing mode has frequency of order nu^{-1/2}, so it lives on the boundary-layer scale and is genuinely three-dimensional: it does not occur when u_s and v_s are linearly dependent.
  • Validity of the 3D Prandtl expansion cannot be expected without analytic regularity, in contrast to the 2D case where Gevrey-3/2 suffices.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A direct numerical check for intermediate spanwise amplitudes, where Lemma 2.2 is silent, would test whether the instability extends beyond the proven parameter range; the theorem itself makes no claim there.
  • If the instability persists at arbitrarily small spanwise amplitude, the bifurcation picture in Remark 1.5 suggests that physically realized 3D boundary layers with weak spanwise drift should show rapid linearized transient growth even when the secondary flow is infinitesimal.
  • The mechanism inverts the usual role of convexity: U''(0)>0 destabilizes rather than stabilizes, which could guide spanwise-forcing design or numerical experiments aimed at suppressing boundary-layer transition.
  • The paper stops at the linear level; an analogous nonlinear instability statement for the 3D Navier-Stokes system would be a natural next step, following the pattern of the 2D Tollmien-Schlichting nonlinear theory.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper studies linearized 3D Navier–Stokes equations around a boundary-layer shear flow u_s(z/√ν), v_s(z/√ν), 0. It seeks normal modes e^{-iαcν^{-1/2}t} e^{iν^{-1/2}(σx+βy)} (ũ,ṽ,w̃)(z/√ν). By decoupling the vertical velocity, the authors reduce the problem to an Orr–Sommerfeld equation for w̃, with an induced base flow U = (σ/α)u_s + (β/α)v_s. Under structural conditions on u_s,v_s guaranteeing U'(0)>0, U''(0)>0, U>0, they construct an inviscid Rayleigh unstable mode via a detailed asymptotic expansion of φ_Ray(0;c), obtaining c_Ray ≈ α U_∞^2/U'(0) + i α^2 U''(0)U_∞^4 π/U'(0)^4. They then use a Rayleigh–Airy iteration to show this instability persists for the viscous Orr–Sommerfeld equation, and recover the tangential velocity components from the vertical one. The main theorem asserts a spectral instability with growth rate e^{t/√ν} for a class of spanwise profiles satisfying the quantitative condition (2.2).

Significance. If the proof is correct, this is a substantial contribution: it gives the first rigorous construction of a 3D boundary-layer instability that is stronger than the classical 2D Tollmien–Schlichting instability, and it indicates that 3D Prandtl expansions fail in general without analytic regularity. The argument is genuinely parameter-free at the level of the dispersion relation: the eigenvalue c_Ray is solved from φ_Ray(0;c)=0, and the asymptotic expansions in Lemmas 3.4 and 3.7 are explicit. The use of Rouché's theorem in Theorem 4.5 is appropriate given the claimed analyticity of φ_Ray(0;c). The paper also gives a self-contained recovery of the tangential velocity components in H^1. The main caveat is that the advertised 'generic' and 'any δ>0' statements go beyond what the written hypotheses and Lemma 2.2 establish.

major comments (2)
  1. [Abstract; Remark 1.5; Proposition 2.1; Lemma 2.2] The abstract says the instability occurs for 'generic boundary layer profiles', and Remark 1.5 says it occurs for any spanwise amplitude δ>0. The theorem actually requires the quantitative structural inequality (2.2), which is not shown to be generic and is not implied by the linear-independence condition (2.1). Lemma 2.2 verifies (2.2) for the exponential family only under k>1 and k|v∞| ≥ max{1, sqrt(2/(k-1))}; the small-amplitude regime is not treated in the paper. Thus 'generic' and 'any δ>0' are unsupported as written. Either prove (2.2) for an open dense set (or for all linearly independent profiles), or revise the abstract and Remark 1.5 to state the actual conditional class.
  2. [§4, definition of H2 and estimates (4.10)–(4.12)] The parameter set H2 is stated as 'α, c_r, c_i ∼ O(1)', but the actual regime is α≪1 and c_i ∼ α^2. This is internally inconsistent: c_i∼O(1) and c_i≪α≪1 cannot both hold. The convergence proof of the Rayleigh–Airy iteration relies on 'Since α, c_i ∼ O(1)' to make the factor |ε|^{1/4} c_i^{-3/2}|log c_i| small. With c_i∼α^2 the factor is ν^{1/8} α^{-13/4}|log α|, so the written proof only gives convergence for ν sufficiently small depending on α. Theorem 4.5 states existence of ν0 for each fixed α, so the argument can be repaired by tracking the α-dependence explicitly, but the current text does not do so. This needs to be fixed in a revision.
minor comments (4)
  1. [Throughout] There are several typos: 'centain' in §1.1, 'suffcient' in Lemma 2.2, 'vetical' in the proof of Theorem 1.1, and inconsistent spelling of 'Rouché'. These should be corrected.
  2. [§1.3 and §4] The phrase 'α, c ∼ O(1)' appears in the roadmap and in H2; it should be replaced by 'α≪1, c_i≪c_r∼α' to match the actual scaling used in Sections 3 and 4.
  3. [Lemma 2.2] The proof of Lemma 2.2 compares with the expression 2(k^3v∞^2 -1)/(k|v∞|(1+k)). This expression is correct only after assuming v∞ is dimensionless and the profiles are as in (2.3); the presentation would be clearer if the nondimensionalization were stated explicitly. In particular, the notation k^3v∞^2 could be misread as a dimensional quantity.
  4. [Theorem 4.5 proof] In the Rouché argument, the bound |φ_Ray(0;c)| ≥ ν^δ/(2U_∞^2) uses |∂_c φ_Ray| = O(α|log c_i| + |c|). Since |c|∼α, the condition α|log α|≪1 is needed; this is true for α≪1 but should be stated explicitly.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the unstable eigenvalue is solved from a dispersion relation; the structural hypotheses are genuine assumptions, not relabeled conclusions.

full rationale

The paper derives, rather than assumes, the instability. The eigenvalue c_Ray is obtained by solving the dispersion relation phi_Ray(0;c)=0 after an explicit asymptotic expansion (Theorem 3.9), and the viscous eigenvalue is obtained by a Rouché argument near c_Ray (Theorem 4.5); no fitted parameter is renamed as a prediction. The base flow U = (sigma/alpha)u_s + (beta/alpha)v_s is constructed from the profiles via Proposition 2.1, but this is a hypothesis selection, not an assumption of the conclusion. Citations to [9,12] are used as a solution toolbox (Rayleigh–Airy iteration) and are external to the authors; citations to the authors' own prior works are background or motivational and are not load-bearing for the main theorem. The known gap between the abstract's 'generic' claim and the quantitative condition (2.2), including the unsupported small-amplitude claim in Remark 1.5, is a correctness/scope concern rather than a circularity: the proof is conditional on (2.2), but the theorem does not assume the existence of the unstable mode. The derivation chain is therefore self-contained.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The central claim rests on the structural profile hypotheses (2.1)-(2.2), the Fourier mode ansatz, and the technical Rayleigh-Airy estimates. There are no invented physical entities. The effective profile U is derived, not postulated.

free parameters (2)
  • alpha (total rescaled wavenumber) = small, fixed as nu -> 0
    The theorem holds for any alpha << 1; the eigenvalue c_Ray is computed as a function of alpha, so alpha is a tunable perturbation parameter, not fitted to data. The growth rate constant scales as alpha^3.
  • lambda_1, lambda_2 (combination weights in U) = explicit formulas from the profile derivatives
    Chosen by Proposition 2.1 to make U positive and convex at the wall; they are not fitted to data but are constructed from the boundary derivatives of u_s and v_s.
assumptions (5)
  • domain assumption The incompressible Navier-Stokes equations in the half-space with no-slip boundary conditions
    System (1.1); the central setup.
  • domain assumption The base flow is a stationary shear flow (u_s(z/sqrt(nu)), v_s(z/sqrt(nu)), 0) maintained by an external force F^nu
    Section 1.2; standard in the spectral instability literature.
  • ad hoc to paper Profiles u_s, v_s satisfy the structural inequalities (2.1)-(2.2), ensuring the induced profile U satisfies (1.9)
    Proposition 2.1; these are restrictive hypotheses that are not shown to cover 'generic' 3D boundary layers.
  • standard math The Rayleigh-Airy iteration converges for the parameter regime H2; bounds in Propositions 3.1 and 4.2
    Sections 3.1 and 4.2; these are proven in the paper but are technical estimates the reader must accept.
  • standard math Rouche's theorem and analyticity of the dispersion relation phi(0;c)
    Theorem 4.5; standard complex analysis.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Unstable mode around the 3D boundary layer flow." pith.science (2026). https://pith.science/paper/4PM2JJHY

@misc{pith2026250906089,
  author       = {Pith},
  title        = {Pith review of: Unstable mode around the 3D boundary layer flow},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4PM2JJHY}},
  note         = {Machine review of arXiv:2509.06089}
}
abstract

We study the stability properties of boundary layer-type shear flows for the three-dimensional Navier-Stokes equations in the limit of small viscosity $0<\nu\ll 1$. When the streamwise and spanwise velocity profiles are linearly independent near the boundary, we construct an unstable mode that exhibits rapid growth at the rate of $e^{t/\sqrt{\nu}}$. Our results reveal an analytic instability in the three-dimensional Navier-Stokes equations around generic boundary layer profiles. This instability arises from the interplay between spanwise flow and three-dimensional perturbations, and does not occur in purely two-dimensional flows.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages

  1. [1]

    Bian and E

    D. Bian and E. Grenier,Onset of nonlinear instabilities in monotonic viscous boundary layers, SIAM J. Math. Anal. Vol. 56, Iss. 3 (2024)

    work page 2024

  2. [2]

    Bian and E

    D. Bian and E. Grenier,Long-wave instabilities, Sci. China Math. Vol. 67, (2024), 1761–1776

    work page 2024

  3. [3]

    Q. Chen, D. Wei and Z. Zhang,Transition threshold for the 3D Couette flow in a finite channel, Mem. Amer. Math. Soc., 296 (2024), Paper No. 1478

    work page 2024

  4. [4]

    Q. Chen, D. Wu and Z. Zhang,On theL ∞ stability of Prandtl expansions in Gevrey class, Sci. China Math., 65(2022), 2521-2562

    work page 2022

  5. [5]

    M. Fei, T. Tao and Z. Zhang,On the zero-viscosity limit of the Navier–Stokes equations inR 3 without analyticity, J. Math. Pures. Appl., 112 (2018), 170-229

    work page 2018

  6. [6]

    Gallaire, D

    F. Gallaire, D. G´ erard-Varet, and F. Rousset,Three-dimensional instability of planar flows, Arch. Ration. Mech. Anal., 186 (2007), 423–475

    work page 2007

  7. [7]

    G´ erard-Varet and Y

    D. G´ erard-Varet and Y. Maekawa,Sobolev stability of Prandtl expansions for the steady Navier-Stokes equations, Arch. Ration. Mech. Anal., 223(2019), 1319-1382

    work page 2019

  8. [8]

    G´ erard-Varet and E

    D. G´ erard-Varet and E. Dormy,On the ill-posedness of the Prandtl equation.J. Amer. Math. Soc. 23, (2010) 591–609

    work page 2010

Show all 36 references
  1. [9]

    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

  2. [10]

    G´ erard-Varet, Y

    D. G´ erard-Varet, Y. Maekawa and N. Masmoudi,Optimal Prandtl expansion around concave boundary layer, Anal. PDE 17 (2024), No. 9, 3125–3187

    2024

  3. [11]

    Grenier,On the nonlinear instability of Euler and Prandtl equations, Comm

    E. Grenier,On the nonlinear instability of Euler and Prandtl equations, Comm. Pure Appl. Math., 53 (2000), 1067–1091

    2000

  4. [12]

    Grenier, Y

    E. Grenier, Y. Guo and T. Nguyen,Spectral instability of characteristic boundary layer flows, Duke Math. J., 165 (2016), 3085-3146

    2016

  5. [13]

    Grenier and T

    E. Grenier and T. Nguyen,L ∞ instability of Prandtl layers, Ann. PDE, 5(2019), Paper No. 18, 36 pp

    2019

  6. [14]

    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

  7. [15]

    Guo and S

    Y. Guo and S. Iyer,Validity of steady Prandtl layer expansions, Comm. Pure Appl. Math., 76 (2023), no. 11, 3150–3232

    2023

  8. [16]

    Iyer and N

    S. Iyer and N. Masmoudi,Global-in-x stability of steady Prandtl expansions for 2D Navier-Stokes flows, arXiv (2020): 2008.12347

    2020 arXiv

  9. [17]

    Kukavica, V

    I. Kukavica, V. Vicol, F. Wang,The inviscid limit for the Navier-Stokes with data only near the boundary. Arch. Ration. Mech. Anal. 237, (2020) 779–827

    2020

  10. [18]

    Li, N, Masmoudi and T

    W.-X. Li, N, Masmoudi and T. Yang,Well-Posedness in Gevrey Function Space for 3D Prandtl Equations without Structural Assumption, Comm. Pure Appl. Math., 75 (2022), no. 8, 1755-1797

    2022

  11. [19]

    Liu , Y.-G

    C.-J. Liu , Y.-G. Wang and T. Yang,A well-posedness theory for the Prandtl equations in three space variables.Adv. Math., 308(2017), 1074–1126

    2017

  12. [20]

    Liu , Y.-G

    C.-J. Liu , Y.-G. Wang and T. Yang,A global existence of weak solutions to the Prandtl equations in three space variables. Discrete Contin. Dyn. Syst. Ser. S, 9(2016), 6, 2011-2029

    2016

  13. [21]

    Liu , Y.-G

    C.-J. Liu , Y.-G. Wang and T. Yang,On the ill-posedness of the Prandtl equations in three space dimen- sions. Arch. Ration. Mech. Anal., 220(2016), 1, 83-108

    2016

  14. [22]

    Maekawa,On the inviscid limit problem of the vorticity equations for viscous incompressible flows in the half-plane, Comm

    Y. Maekawa,On the inviscid limit problem of the vorticity equations for viscous incompressible flows in the half-plane, Comm. Pure. Appl. Math., 67 (2014), 1045–1128

    2014

  15. [23]

    Maekawa and A

    Y. Maekawa and A. Mazzucato: The inviscid limit and boundary layers for Navier-Stokes flows. Handbook of mathematical analysis in mechanics of viscous fluids, 781–828. UNSTABLE MODE AROUND THE 3D BOUNDARY LAYER FLOW 35

  16. [24]

    Masmoudi, Y

    N. Masmoudi, Y. Wang, D. Wu and Z. Zhang,Mack modes in supersonic boundary layer, arXiv:2405.04853

    arXiv

  17. [25]

    Masmoudi, Y

    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

    2024

  18. [26]

    Mazzucato and M

    A. Mazzucato and M. Taylor,Vanishing viscosity limits for a class of circular pipe flows, Comm. Partial. Differential. Equa., 36 (2010), 328-361

    2010

  19. [27]

    Nguyen and T

    T. Nguyen and T. Nguyen,The inviscid limit of Navier-Stokes equations for analytic data on half space. Arch. Ration. Mech. Anal. 230, (2018), 1103-1129

    2018

  20. [28]

    Prandtl,Uber Flussigkeitsbewegung bei sehr kleiner Reibung, Verhandl

    L. Prandtl,Uber Flussigkeitsbewegung bei sehr kleiner Reibung, Verhandl. 3rd Int. Math. Kongr. Heidel- berg (1904), Leipzig (1905)

    1904

  21. [29]

    H. L. Reed and W. S. Saric,Stability of three-dimensional boundary layers, Ann. Rev. Fluid Mech. 1989. 21: 235–84

    1989

  22. [30]

    Sammartino and R.E

    M. Sammartino and R.E. Caflisch,Zero viscosity limit for analytic solutions of the Navier-Stokes equation on a half-space. II. Construction of the Navier-Stokes solution. Comm. Math. Phys., 192 (1998), 463-491

    1998

  23. [31]

    W. S. Saric, H. L. Reed, and E. B. WhiteStability and transition of three-dimensional boundary layers, Ann. Rev. Fluid Mech. (2003) 35: 413–40

    2003

  24. [32]

    W. S. Saric and L. G. Yeates, Generation of crossflow vortices in a three-dimensional flat-plate flow. In: Kozlov, V.V. (Ed.) (1985): Laminar-Turbulent Transition. Springer-Verlag, Berlin, 429–437

    1985

  25. [33]

    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

  26. [34]

    C. Wang, Y. Wang and Z. Zhang,Zero-viscosity limit of the Navier–Stokes equations in the analytic setting, Arch. Ration. Mech. Anal., 224 (2017), 555-595

    2017

  27. [35]

    C. Wang, Y. Wang and Z. Zhang,Zero-viscosity limit of the compressible Navier–Stokes equations in the analytic setting, arXiv (2023): 2305. 09393

    2023

  28. [36]

    Yang and Z

    T. Yang and Z. Zhang,Linear instability analysis on compressible Navier-Stokes equations with strong boundary layer, Arch. Ration. Mech. Anal., 247 (2023), Paper No. 83, 53 pp. (C.-J. Liu)School of Mathematical Sciences, LSC-MOE, CMA-Shanghai and Institute of Nat- ural Science...

    2023

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.