Pith. sign in

REVIEW 3 major objections 4 minor 38 references

Level crossing instabilities in inviscid isothermal compressible Couette flow

T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

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

Pith's one-line read Repeated level crossings make compressible Couette flow unstable in stripes.

desk verdict A solid analytical core with clean stability theorems and explicit mode formulas, but the headline infinite zebra-stripe pattern rests on a numerical search that lacks convergence checks and is extrapolated from low modes. read the letter →

arxiv 2412.20813 v1 pith:SW6Q47NY submitted 2024-12-30 physics.flu-dyn math-phmath.MP

Govind S. Krishnaswami , Sonakshi Sachdev , Pritish Sinha This is my paper · ORCID
classification physics.flu-dynmath-phmath.MP MSC 76E0576N10 PACS 47.20.Ft47.40.-x
keywords Couetteflowlinearstabilitycompressibleisothermallevelcrossingzebrastripescontinuousspectrumnoninflectionalmodesgrowthrate
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

The paper studies the linear stability of an inviscid, isothermal, compressible shear flow with a linear velocity profile, inviscid Couette flow, in a finite channel. In the incompressible limit this flow has no unstable modes, but once compressible perturbations are allowed, an infinite tower of noninflectional discrete modes appears. The paper argues that as the Mach number $M$ or wavenumber $k$ grows, neighboring purely oscillatory modes repeatedly collide on the imaginary axis, split into complex pairs, and thereby open finite windows of instability. These windows alternate with stable regions and form zebra-like stripes in the $k$-$M$ plane, with growth rates growing like the square root of distance from each transition. The picture matters because it shows that compressibility alone, without viscosity or inflection points, can produce a rich, organized pattern of linear instabilities.

What carries the argument

The central object is the vortex-perturbation equation $\tilde{w}'' - \frac{2i\tilde{k}}{\Gamma}\tilde{w}' - (\tilde{k}^2 + M^2\Gamma^2)\tilde{w} = 0$, where $\Gamma = \tilde{\gamma} + i\tilde{k}\tilde{y}$ is the Doppler-shifted growth rate and $\tilde{w}$ is the vorticity perturbation; equivalently, a second-order equation for the normal velocity with Dirichlet boundary conditions at the channel walls. These equations are transformable into a confluent hypergeometric (Whittaker) equation, which expresses the same spectrum but is not solved analytically. The argument is carried by the merger-demerger of adjacent imaginary eigenvalues as $\tilde{k}$ or $M$ varies, captured near each transition by the biquadratic characteristic equation $(\tilde{\gamma}^2 - \tilde{\gamma}_c^2 + A^2(M_c-M))^2 + 4\tilde{\gamma}_c^2 A^2(M_c-M) = 0$, whose two branches produce the square-root growth law. A $\mathbb{Z}_2\times\mathbb{Z}_2$ symmetry constrains eigenvalues to four-tuples and fixes the allowed eigenfunction forms.

What would settle it

Calculate the full spectrum at larger $M$ and $k$ with an independent high-order method, for example by solving the transcendental equation from the confluent hypergeometric form, and check whether the first instability arch for $\tilde{k}=1$ opens at $M\approx 4.203$ and whether the imaginary and real parts of $\tilde{\gamma}$ follow the square-root law with fitted slope near $0.5$; a missing arch, a spurious arch, or a mode with $|\Im\tilde{\gamma}|<\tilde{k}$ that never crosses would show that the infinite sequence of stripes is not real.

Watch

Extended reading notes

Core claim

For inviscid isothermal compressible Couette flow, the linearized eigenvalue problem reduces to a single second-order ODE for vorticity perturbations with the growth rate appearing nonlinearly. In the small-wavenumber, small-Mach, and high-frequency limits, the equation has an infinite tower of neutrally stable modes labeled by an integer $n$; each has a conjugate partner, and modes come in four-tuples. Beyond those limits, the purely imaginary eigenvalues of neighboring modes meet on the imaginary axis, leave the axis as a complex-conjugate pair, and later merge again, producing an arch-shaped interval of instability. The paper identifies this repeated level-crossing mechanism as the source of an infinite sequence of stability transitions, derives analytic stability theorems and bounds that delineate the stable regions, and fits the numerically computed eigenvalues near each crossing to a canonical biquadratic characteristic equation that gives square-root power-law behavior. A separate continuous spectrum of neutrally stable eigenmodes exists along the imaginary axis, with eigenfunctions that are only piecewise smooth across critical layers.

Load-bearing premise

The claim that instabilities repeat indefinitely rests on a numerical eigenvalue search that is assumed to track the true spectrum of a non-self-adjoint operator without missing or spurious modes, extrapolating from the first few arches to all modes and all parameter ranges.

Editorial extensions

If this is right

  • For small $k$ or small $M$ the flow is guaranteed neutrally stable; the paper proves $M<1/2$ is sufficient for stability at all wavenumbers.
  • Each mode $n$ undergoes an infinite sequence of stable-unstable transitions as $M$ or $k$ grows, so the number of participating modes increases while the maximum growth rate in each arch decreases.
  • The first instability always involves the ground mode colliding with its conjugate, and for the ground mode the resulting windows alternate between $(0^*,0)$ and $(0,1)$ mergers.
  • Near every transition, the real and imaginary parts of the growth rate scale as the square root of distance in $M$ or $k$, with a coefficient fixed by the biquadratic ansatz.
  • The continuous spectrum fills the imaginary interval $-\tilde{k}<\Im\tilde{\gamma}<\tilde{k}$ and overlaps the discrete spectrum in the level-crossing regime, but its eigenfunctions are not smooth across critical layers.

Reading between the lines

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

  • If the infinite-tower extrapolation is correct, a similar zebra-stripe pattern should appear for every mode, shifted toward larger $k$ and $M$; locating the first unstable band for a high mode $n$ would test the extrapolation.
  • Because the discrete operator is nonnormal, the paper's exponential-mode analysis does not rule out finite-time transient growth that could blur the sharp linear stability boundaries in numerical or experimental settings.
  • The level-crossing mechanism is not tied to isothermal thermodynamics or two-dimensional perturbations: the same merger-demerger route should appear in adiabatic and three-dimensional perturbations, so the stripes may be a general organizing feature of compressible shear-flow stability.
  • A direct probe of the $k$-$M$ plane by measuring the response to monochromatic perturbations could look for the predicted alternating stable and unstable bands in an experiment or simulation.
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

3 major / 4 minor

Summary. The paper studies the linear stability of inviscid, isothermal, compressible Couette flow in a finite channel. It derives a compressible analogue of Rayleigh's equation, reduces it to second-order ODEs for the vorticity and velocity perturbations, and establishes symmetry properties of the spectrum. For uniform-vorticity Couette flow it proves several stability criteria: a critical-layer condition for instability, a semicircle-type bound, lower bounds on |Im(gamma)| for neutral modes, and a condition for neutral modes to possess critical layers. It then constructs analytic infinite towers of neutral modes in the small-wavenumber, small-Mach-number, and high-frequency limits. Beyond these regimes, a numerical Fredholm-alternative search is used to find complex growth rates, showing repeated level crossings between adjacent modes that open and close instability arches and produce a claimed zebra-stripe pattern of instability regions in the (k,M) plane. A canonical square-root power law near crossings is proposed and fit to the numerics. A continuous spectrum of non-smooth neutral modes is also described.

Significance. If the numerical results are correct, the paper would provide a striking new picture of noninflectional compressional instabilities in inviscid Couette flow, complementing earlier work by Glatzel and Renardy and extending the classical Rayleigh/Howard framework to isothermal compressible flows. The analytic parts are the strongest contribution: the integral identities in Section 4.2 and the asymptotic towers in Section 5 are clean, internally consistent, and match the numerics in overlapping regimes. The paper also gives a useful semicircle theorem and explicit neutral-stability bounds. However, the central quantitative claims about infinite level-crossing sequences and zebra stripes currently rest on an unvalidated numerical search and on extrapolation from a handful of low-lying arches; no code, data, convergence study, or independent verification is provided. The canonical square-root law is an ansatz fitted to the same numerical spectra it aims to explain. With additional validation, the paper would be a solid contribution; in its present form the main instability-pattern claims remain insufficiently supported.

major comments (3)
  1. [Appendix B, Eq. (104)] The entire level-crossing scenario of Section 6 rests on the numerical Fredholm-alternative search, but the manuscript gives no convergence study, no grid-refinement test, and no independent check of the criterion that an eigenvalue is found when the norm of the solution 'tends to become unbounded.' The operator T_gamma is non-self-adjoint, and precisely in the level-crossing regime |Im(gamma)| < k the leading coefficient k^2 + M^2 Gamma^2 in (104) can have real zeros inside the channel, so the Fredholm hypotheses are not verified where the method is used. Spectral pollution is a known hazard for such problems. I request a convergence study, error estimates, a comparison with an independent method such as a spectral collocation scheme or the hypergeometric transcendental equation (108), and release of the code and data used to produce Figures 3-7.
  2. [Section 6.3, Fig. 5] The claim of an 'infinite sequence' of alternating instability stripes is extrapolated from the first few arches of the n = 0 mode, computed on only four parameter slices: k = 1 and k = 1.6 for variable M, and M = 4.75 and M = 7.6 for variable k. Figure 5 is itself a sketch using solid segments from those slices, with no computations shown for higher arches, higher modes, or larger parameter values. The paper explicitly states in Section 2 that the authors 'hope to prove the key features of our numerical spectra elsewhere,' which confirms that the unbounded repetition is currently a conjecture. This point is load-bearing for the abstract's claims of an infinite sequence of stability transitions and a zebra-stripe pattern; it should either be proven or clearly presented as a conjecture.
  3. [Section 6.4, Eqs. (86)-(92)] The canonical square-root law is proposed as an ansatz rather than derived from the eigenvalue problem. The coefficients A and B in (86)-(89) and the interpolating arch parameters in (91) are fitted to the same numerical spectra that the ansatz is supposed to explain, so the agreement in Figures 6 and 7 is a consistency check rather than independent confirmation. In particular, the assumption that only the 2x2 block of the two crossing modes contributes to the characteristic equation near a crossing is introduced without justification, and the biquadratic form (90) is not derived from the operator Q. The manuscript should either derive the 2x2 reduction from the resolvent structure or demonstrate that the fitted A-values are parameter-free predictions across multiple crossings and parameter values.
minor comments (4)
  1. [Appendix B title] The heading contains a typo: 'F redholm alternative' should read 'Fredholm alternative.'
  2. [Section 7] The regularity class 'C2' in the text and the remark should be written as C^2, and the sentence 'which is a regularity class' should be clarified: 'C^2' is standard notation, but 'C2' is ambiguous.
  3. [Fig. 5 caption] The caption refers to solid and dashed boundary curves, but the figure as printed does not clearly distinguish them; a legend or explicit labels for the lines M = 1/2 and the curve M^2 = 1/4 + 1/(8k^2) would improve readability.
  4. [Section 6.4] The notation Mc1, Mc2, and gamma_mid in Eq. (91) is introduced only in the surrounding prose; defining all symbols in the displayed equation itself would make the interpolation formula easier to follow.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; the stability analysis is derived from the Euler equations, and the level-crossing instabilities come from an independent numerical search.

full rationale

The eigenvalue problem is derived from the isothermal Euler equations and impenetrable boundary conditions (Eqs. 17, 38, 42), and the stability theorems in Sec. 4.2 follow from integral identities applied to those equations. The infinite tower of modes in Sec. 5 follows from asymptotic solutions of the ODE in the small-k, small-M, and large-|Im gamma| limits, not from the numerical results. The level-crossing instabilities in Sec. 6 are found by a Fredholm-alternative search, which is an independent numerical method, and the paper explicitly notes that the results are numerical spectra rather than proven theorems ('Elsewhere, we hope to prove the key features of our numerical spectra', Sec. 2). The canonical square-root form (86) is introduced as a hypothesis and then tested by fitting the numerically determined eigenvalues; Eq. (92) computes the coefficient A from the same numerically determined arch and compares it with local log-log intercepts. This is an internal consistency check of the assumed functional form rather than an independent prediction, and the central claims of the paper do not rest on this comparison. No load-bearing self-citations or definitional reductions were found; the extrapolation from a few computed arches to an infinite zebra-stripe pattern is an evidentiary limitation, not a circularity.

Assumptions & free parameters 3 free parameters · 7 assumptions · 0 invented entities

The central derivations rest on standard fluid-dynamic modeling assumptions (inviscid isothermal ideal gas, 2D normal modes, channel walls) and on a numerical search method whose reliability is asserted but not demonstrated. The only fitted quantities are parameters of the phenomenological square-root ansatz, not of the underlying eigenvalue problem.

free parameters (3)
  • A (canonical crossing amplitude) = log A about -1.97 to -2.11 for M=4.203 arch; log A about -1.20 to -1.33 for k-M arch
    Introduced in Eq. (86) as the coefficient of the square-root term; fitted to numerical eigenvalue data in Figs. 6 and 7.
  • B (linear correction coefficient) = not reported
    Added in Eq. (89) to capture the linear drift of Im gamma; fitted to numerical arches.
  • Interpolating arch parameters (Mc1, Mc2, gamma_c1, gamma_c2, gamma_mid) = Mc1=4.203, Mc2=5.4535 for k=1; kc1=0.8755, kc2=1.1415 for M=4.75
    Taken directly from the numerical spectrum to build ansatz (91); these are empirical inputs, not derived.
assumptions (7)
  • domain assumption Inviscid isothermal ideal gas Euler equations with constant temperature
    Governing equations (1)-(3); the whole analysis is restricted to this model.
  • domain assumption Normal mode ansatz with exponential time dependence e^{gamma t}
    Eq. (7); excludes algebraic or transient growth, which the authors note could matter for nonnormal operators.
  • domain assumption Two-dimensional perturbations are the least stable (Squire's theorem extrapolation)
    Section 3.2; the 3D case is argued to be similar but not analyzed.
  • domain assumption Impenetrable channel walls y=+-L with v(+-L)=0
    Sections 3.1 and 4.1; the discrete spectrum depends on these boundary conditions.
  • ad hoc to paper Only the 2x2 block of the two crossing modes contributes to the characteristic equation near a level crossing
    Section 6.4 around Eq. (86); justifies the quadratic and biquadratic ansatz but is not derived from the operator.
  • ad hoc to paper The Fredholm alternative applies to T_gamma and the grid search finds all eigenvalues without spurious results
    Appendix B; no convergence proof or error analysis is given.
  • domain assumption Linear background velocity u=-Omega y (uniform vorticity Couette flow)
    Section 3.1 and Eq. (5); all specific results are for this profile.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Level crossing instabilities in inviscid isothermal compressible Couette flow." pith.science (2026). https://pith.science/paper/SW6Q47NY

@misc{pith2026241220813,
  author       = {Pith},
  title        = {Pith review of: Level crossing instabilities in inviscid isothermal compressible Couette flow},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SW6Q47NY}},
  note         = {Machine review of arXiv:2412.20813}
}
abstract

We study the linear stability of inviscid steady parallel flow of an ideal gas in a channel of finite width. Compressible isothermal two-dimensional monochromatic perturbations are considered. The eigenvalue problem governing density and velocity perturbations is a compressible version of Rayleigh's equation and involves two parameters: a flow Mach number $M$ and the perturbation wavenumber $k$. For an odd background velocity profile, there is a $\mathbb{Z}_2 \times \mathbb{Z}_2$ symmetry and growth rates $\gamma$ come in symmetrically placed 4-tuples in the complex eigenplane. Specializing to uniform background vorticity Couette flow, we find an infinite tower of noninflectional eigenmodes and derive stability theorems and bounds on growth rates. We show that eigenmodes are neutrally stable for small $k$ and small $M$ but that they otherwise display an infinite sequence of stability transitions with increasing $k$ or $M$. Using a search algorithm based on the Fredholm alternative, we find that the transitions are associated to level crossings between neighboring eigenmodes. Repeated level crossings result in windows of instability. For a given eigenmode, they are arranged in a zebra-like striped pattern on the $k$-$M$ plane. A canonical square-root power law form for $\gamma(k,M)$ in the vicinity of a stability transition is identified. In addition to the discrete spectrum, we find a continuous spectrum of eigenmodes that are always neutrally stable but fail to be smooth across critical layers.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

38 extracted references · 37 canonical work pages

  1. [1]

    Drazin P G and Reid W H, Hydrodynamic instability, 2 nd Ed., Cambridge University Press (2004). 2, 3, 7

    work page 2004

  2. [2]

    Part 2 , J

    Blumen W, Drazin P and Billings D, Shear layer instability of an inviscid compressible fluid. Part 2 , J. Fluid Mech., 71(2), 305 (1975). 2

    work page 1975

  3. [3]

    Choudhury S R and Lovelace R V E, On the Kelvin-Helmholtz instabilities of supersonic shear layers, Astrophysical Journal Part 1, 283, 331 (1984). 2

    work page 1984

  4. [4]

    Rayleigh J W S, The Theory of Sound , Dover (1929). 2, 3, 9

    work page 1929

  5. [5]

    Howard L N, Note on a paper of John W Miles , J Fluid Mech., 10, 609 (1961). 2

    work page 1961

  6. [6]

    Drazin P G and Howard L N, Hydrodynamic Stability of Parallel Flow of Inviscid Fluid , Advances in Applied Mechanics, 9, 1 (1966). 2

    work page 1966

  7. [7]

    Note 1115 (1946)

    Lees L and Lin C C, Investigation of stability of laminar boundary layer in a compressible fluid, National Advisory Committee for Aeronautics, Tech. Note 1115 (1946). 2

    work page 1946

  8. [8]

    Fluids, 6, 1042 (1963)

    Eckart C, Extension of Howard’s circle theorem to adiabatic jets , Phys. Fluids, 6, 1042 (1963). 2

    work page 1963

Show all 38 references
  1. [9]

    de Mecanique, 16, 227 (1977)

    Shivamoggi B K, Inviscid theory of stability of parallel compressible flows , J. de Mecanique, 16, 227 (1977). 2

    1977

  2. [10]

    Shivamoggi B K, Effects of compressibility on stability characteristics of a free shear layer, Z. Angew. Math. Mech. 59, 405 (1979). 2

    1979

  3. [11]

    Fluids, 31, 3239 (1988)

    Djordjevic V D and Redekopp L G, Linear stability analysis of nonhomentropic, inviscid compressible flows, Phys. Fluids, 31, 3239 (1988). 2

    1988

  4. [12]

    Subbiah M and Jain R K, Stability of Compressible Shear Flows , Journal of Mathemat- ical Analysis and Applications, 151(1), 34 (1990). 2

    1990

  5. [13]

    Padmini M and Subbiah M, Stability of Non-Homentropic, Inviscid, Compressible Shear Flows, Journal of Mathematical Analysis and Applications, 241(1), 56 (2000). 2

    2000

  6. [14]

    1-2 (2011)

    Kishan H, Devi G and Dua N K, Linear instability of an inviscid compressible parallel shear flow, Pure and Applied Mathematika Sciences, LXXIII, No. 1-2 (2011). 2

    2011

  7. [15]

    R-709, 3.1 (1984)

    Mack L M, Boundary-Layer Linear Stability Theory , in Special Course on Stability and Transition of Laminar Flow , AGARD Rep. R-709, 3.1 (1984). 2, 3

    1984

  8. [16]

    2, 5, 13

    Shivamoggi B K, Stability of parallel gas flows , Ellis Horwood (1986). 2, 5, 13

    1986

  9. [17]

    Fluids, 8, 143 (1960)

    Case K M, Stability of inviscid plane Couette flow , Phys. Fluids, 8, 143 (1960). 2, 28

    1960

  10. [18]

    Glatzel W, The linear stability of viscous compressible plane Couette flow, J Fluid Mech., 202, 515 (1989). 2

    1989

  11. [19]

    Fluid Mech

    Duck P W, Erlebacher G and Hussaini Y, On the linear stability of compressible plane Couette flow , J. Fluid Mech. 258, 131 (1994). 2

    1994

  12. [20]

    Rayleigh J W S, On the stability, or instability, of certain fluid motions , Proc. Lond. Math. Soc. 9, 57 (1880). 3

  13. [21]

    3, 8, 13, 31

    Glatzel W, Sonic instabilities in supersonic shear flows , Monthly Notices of the Royal Astronomical Society, 231(3), 795 (1988). 3, 8, 13, 31

    1988

  14. [22]

    Shivamoggi B K and Pohl L, Polytropic Gas Effects in Parker’s Solar Wind Model and Coronal Hole Flows , arXiv:2407.06122. 3

    arXiv

  15. [23]

    Renardy M, Instability proof for some transonic problems with resonant mode crossings , Theoret. Comput. Fluid Dynamics, 7, 457 (1995). 3

    1995

  16. [24]

    Trefethen L N, Trefethen A E, Reddy S C and Driscoll T A, Hydrodynamic stability without eigenvalues , Science 261, 5121, 578 (1993). 4, 8 32

    1993

  17. [25]

    PDE, 7, 24 (2021)

    Antonelli P, Dolce M, Marcati P, Linear stability analysis of the homogeneous Couette flow in a 2D isentropic compressible fluid , Ann. PDE, 7, 24 (2021). 4, 8

    2021

  18. [26]

    Hunt J C R and Durbin P A, Perturbed vortical layers and shear sheltering , Fluid Dyn. Res. 24, 375 (1999). 4

    1999

  19. [27]

    Plasma Phys

    Thyagaraja A, Loureiro N and Knight P J, Spectral and evolutionary analysis of advection-diffusion equations and the shear flow paradigm, J. Plasma Phys. 68, 363 (2002). 4

    2002

  20. [28]

    Lagerstrom P A, Laminar Flow Theory, Theory of Laminar Flows - Vol IV, High Speed Aerodynamics and Jet Propulsion, Section B, Oxford University Press (1964). 9

    1964

  21. [29]

    Dyson F J, Stability of an Idealized Atmosphere. II. Zeros of the Confluent Hypergeo- metric Function, Phys. Fluids, 3, 155 (1960). 13

    1960

  22. [30]

    Vulovi´ c V Z and Prange R E, Randomness of a true coin toss , Phys. Rev. A, 33, 576 (1986). 25

    1986

  23. [31]

    Drazin P G, Introduction to hydrodynamic stability, Cambridge University Press (2002). 28

    2002

  24. [32]

    Press (1927)

    Whittaker E T and Watson G N, A course of modern analysis , 4 th Ed., Cambridge Univ. Press (1927). 30

    1927

  25. [33]

    Ramm A G, A Simple Proof of the Fredholm Alternative and a Characterization of the Fredholm Operators, American Mathematical Monthly, 108, 855 (2001). 31

    2001

  26. [34]

    Rudin W, Functional analysis, McGraw Hill (1973). 31

    1973

  27. [35]

    Douglas J, Thyagaraja A and Kim E, A new resolvent technique for calculating linear eigenspectra in tokamaks, in McKenna C (Ed.), 37 th EPS Conference on Plasma Physics 2010, Curran Associates inc., 17–20 (2010). 31

    2010

  28. [36]

    Fusion 55, 053016 (2015)

    Chandra D, Thyagaraja A, Sen A, Ham C J, Hender T C, Hastie R J, Connor J W, Kaw P and Mendonca J, Modelling and analytic studies of sheared flow effects on tearing modes, Nucl. Fusion 55, 053016 (2015). 31

    2015

  29. [37]

    Abramowitz M and Stegun I, Handbook of Mathematical Functions, Dover (1970). 31

    1970

  30. [38]

    Coddington E A and Levinson N, Theory of ordinary differential equations , McGraw Hill (1955). 31 33

    1955

Pith tools

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