{"id":"9dca637a-e1e8-4f5e-8195-d4ec16a64098","arxiv_id":"2505.21650","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":18,"one_line_summary":"A weakly nonlinear Landau model for compressible Kelvin-Helmholtz instability is proposed, but its coefficients are not derived and the linearized base flow lacks the shear that creates the instability.","lead":"This paper claims to extend Blumen's linear stability theory of compressible shear layers into the weakly nonlinear regime, deriving a Landau-Stuart amplitude equation with Mach-number-dependent coefficients. The derivation omits the shear term in the linearized equations, and the reported growth rates do not follow from the displayed dispersion relation, so the claimed bifurcation predictions are not supported.","discovery_kind":"incremental","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The linearized equations drop the shear term v1 d\\bar{u}/dy that drives Kelvin-Helmholtz instability, leaving only real acoustic waves, so the claimed bifurcations are not about shear layers.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing defect: the first-order linearization in Eqs. (11)–(13) drops the shear term v1 d\\bar{u}/dy that is present in Eq. (5) and that is physically responsible for Kelvin-Helmholtz instability. My independent check of the algebra confirms the omission. With that term absent, the dispersion relation (21) is real-valued and describes only acoustic waves in a uniform stream; there is no mechanism for exponential growth, so σ in Figure 2 is unexplained. Since all derivations of ζ and μ in Table I and the entire bifurcation analysis rely on this flawed linear problem, the paper's central claim—that weakly nonlinear analysis yields alternating supercritical and subcritical Hopf bifurcations in compressible KHI—does not follow from the equations it presents. The paper does include substantial derivational effort in Appendix A, and it correctly states the form of the full momentum equation, but the decisive linearization error means the results are not results about shear layers. I also considered whether the paper might be treating a uniform base flow and calling it a shear layer; even granting that, Figure 3 explicitly contrasts with Blumen's tanh-profile results, and the boundary conditions in Section II specify a sheared base flow, so the internal inconsistency stands. The reader's verdict of REJECT is appropriate; no adjustment is needed. A concrete numerical test of the Rayleigh equation would settle any residual doubt, but the algebraic omission is already decisive.","tokens_in":20444,"tokens_out":3360,"duration_ms":32468,"concrete_test":"Re-derive the O(ε) perturbation equations retaining the shear term v1 d\\bar{u}/dy in the x-momentum equation. For a shear layer with \\bar{u}(y) = tanh(y) (the Blumen profile), substitute normal modes and obtain the compressible Rayleigh-type eigenvalue problem; solve it numerically. Check whether Im(ω) > 0 exists for some kx, ky, and M. If unstable eigenvalues are found, compare their growth rates with Figure 2; if no instability is found, the paper's Figure 2 and all Landau coefficient results are not solutions of the stated equations. Alternatively, verify directly whether the dispersion relation (21) satisfies the full O(ε) equations when d\\bar{u}/dy ≠ 0; analytic substitution will show it does not.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the Landau-Stuart equation (26) with Mach-dependent coefficients captures nonlinear saturation and alternating Hopf bifurcations of compressible shear-layer Kelvin-Helmholtz instability. For this to hold, the underlying linear problem must actually contain the velocity gradient responsible for KHI. The full x-momentum equation (5) includes the term v ∂\\bar{u}/∂y; after linearizing with u = \\bar{u} + εu1, v = εv1, this term becomes v1 d\\bar{u}/dy. The paper's first-order equations (11)–(13) omit this term, effectively setting d\\bar{u}/dy = 0. Equations (15)–(17) then reduce to a constant-coefficient algebraic system whose dispersion relation (21), ω = kx\\bar{u} ± k/M, is purely real. Consequently, the growth rate σ = Im(ω) is identically zero in the model, yet Figure 2 and the surrounding text plot positive and negative σ as a function of M and k—a direct internal contradiction. Because the first-order solution used to construct the second-order nonlinear terms (Appendix A) and to evaluate the Landau coefficient ζ is the acoustic wave from (21), none of the subsequent Landau coefficients or bifurcation diagrams can be properties of shear-layer instability. The claim that the paper generalizes Blumen's shear-layer results is thus unsupported; it is instead a study of amplitude dynamics around a stable acoustic mode with no physical connection to the Kelvin-Helmholtz mechanism.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper claims to extend Blumen's linear stability theory of inviscid compressible shear layers to the weakly nonlinear regime. Using a multiple-scales expansion of the compressible Euler equations, the authors derive a Landau-Stuart equation for the disturbance amplitude, compute Mach-number-dependent Landau and linear growth coefficients, and present phase portraits, bifurcation diagrams, and amplitude evolutions that indicate alternating supercritical and subcritical Hopf bifurcations of the Kelvin-Helmholtz instability. The analysis incorporates finite-boundary conditions and is framed as a generalizable framework for nonlinear dynamics of compressible shear flows.","tokens_in":20940,"tokens_out":5991,"duration_ms":55611,"significance":"If the derivation were correct, the paper would offer a valuable systematic weakly nonlinear framework for compressible shear-layer instability, extending a classical linear result and providing concrete predictions about long-time state selection. The paper includes a physically motivated problem setup and a broad set of figures illustrating dynamical regimes. However, the central derivation is flawed: the first-order linearized equations omit the shear term that defines the Kelvin-Helmholtz mechanism, and the Landau coefficients are not derived from the model. As a result, the claimed bifurcation structure and dynamics are not established for compressible shear layers, and the significance is not realized in the present form.","major_comments":[{"comment":"The first-order x-momentum equation (11) omits the term v1 d\\bar{u}/dy that is present in the full x-momentum equation (5). With this term absent, the linearized problem reduces to a constant-coefficient algebraic system, and the dispersion relation (21), ω = kx\\bar{u} ± sqrt(kx²+ky²)/M, is purely real. Hence Im(ω) = σ = 0 identically, which contradicts Figure 2, where positive and negative growth rates are plotted. Since the velocity gradient is the physical source of Kelvin-Helmholtz instability, the subsequent nonlinear analysis is not a study of shear-layer instability.","section":"Sec. III.A, Eqs. (11)-(13) and (21)"},{"comment":"The coefficients μ and ζ in Eq. (26) are stated to be 'obtained from linear stability analysis' and 'estimated values,' but no derivation is shown. The solvability condition (25) is never evaluated, and the values in Table I cannot be obtained from the dispersion relation (21), which has no imaginary part. The phase portraits and bifurcation diagrams in Figs. 5-9 are therefore consequences of the input coefficients, not of the compressible Euler equations. This is a circular argument.","section":"Sec. III.B and Table I"},{"comment":"The second-order problem is internally inconsistent. The normal mode assumption (A4) treats the second-order amplitudes as constants, yet Eqs. (A6)-(A7) and the resulting ODE (A9) contain y-derivatives such as \\hat{p}'_2 and \\hat{v}'_2. Additionally, the term v1 d\\bar{u}/dy appears at O(ε²) in (A1), whereas expansion of Eq. (5) shows that it belongs at O(ε). These inconsistencies invalidate any Landau coefficient derived from this system.","section":"Appendix A, Eqs. (A1)-(A9)"},{"comment":"The base flow is specified only at the boundaries, U(±h)=±1, with no profile between, and Sec. III.A states that u0=\\bar{u} represents a steady, uniform base flow. If \\bar{u} is uniform, there is no shear and no Kelvin-Helmholtz mechanism; if it is y-dependent, the linearized problem must be solved as an eigenvalue problem in y, not as the algebraic system (15)-(17). The paper does not resolve this ambiguity, which undermines the connection to Blumen's shear-layer analysis.","section":"Sec. II, base flow specification"}],"minor_comments":[{"comment":"The term 'Mach number' is repeatedly misspelled as 'Mech number' (e.g., Fig. 3 caption, Sec. IV).","section":"Throughout"},{"comment":"The 'Neumann condition' is misspelled as 'Newmann condition'.","section":"Sec. II"},{"comment":"The interpolation formula for the critical Mach number has unbalanced parentheses and is difficult to parse.","section":"Eq. (27)"},{"comment":"The caption describes a red dashed curve of maximum growth rate, but with Im(ω)=0 from Eq. (21) there is no growth rate to maximize.","section":"Fig. 2"},{"comment":"Reference 4 contains a typographical error: 'N. . Chaturvedi' has an extra period.","section":"References"},{"comment":"The introduction states that the Landau coefficient is 'computed explicitly as a function of compressibility,' but no explicit expression or derivation is provided anywhere in the manuscript.","section":"Sec. I"}],"recommendation":"reject","confidential_remarks":"The manuscript's central derivation is not valid: the linear problem omits the shear term that drives Kelvin-Helmholtz instability, and the Landau coefficients appear to be prescribed rather than derived. Even a major revision would require redoing the linear stability analysis from scratch, including a proper treatment of the y-dependent base flow, which is beyond the scope of a normal revision. The paper is not suitable for publication in its present form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, bottom line: this paper's stated subject is Kelvin-Helmholtz instability in compressible shear layers, but the linearized equations it solves are for a uniform stream. The shear term v1 d\\bar{u}/dy that appears in the full momentum equation (5) is dropped in the first-order set (11)-(13). The resulting dispersion relation, ω = kx \\bar{u} ± k/M, is purely real. There is no instability in that problem, yet Figure 2 shows growth-rate contours. Those contours cannot come from the equations in the paper.\n\nThe paper does a few things cleanly: the multiple-scales setup is standard, the acoustic-wave dispersion relation is correctly derived, and the text is readable in the early sections. The comparison figures with Blumen are at least an attempt to connect to the literature. But that's where the credit ends.\n\nThe Landau-Stuart equation (26) is the generic amplitude equation, and the coefficients in Table I are simply entered by hand. The paper never derives ζ from the solvability condition it writes down; it just presents values. The resulting bifurcation diagrams are therefore a parameter study of the Landau equation with zeta chosen to alternate sign, not a prediction from the compressible shear-layer problem. The finite-boundary claim is also unsupported: a single plane wave with nonzero ky cannot satisfy v=0 at y=±h, and the paper does not solve a boundary-value problem in y.\n\nThe second-order equations in Appendix A do retain v1 d\\bar{u}/dy, but that doesn't repair the linear operator. If the linear problem is stable, nonlinear terms cannot create a shear instability in a weakly nonlinear expansion around that state.\n\nWho is this for? Nobody, really, in its current form. A reader familiar with stability theory will see the missing shear term in the first two pages. The proper fix, redoing the calculation with a hyperbolic-tangent base profile and a correct eigenfunction problem, is a new paper.\n\nRecommendation: desk reject. It does not deserve referee time because the central derivation is internally inconsistent with the governing equations, and the Mach-dependent bifurcations are fitted inputs, not derived outputs.","headline":"The linearized equations drop the shear term that defines Kelvin-Helmholtz instability, so the paper's central results are about acoustic waves in uniform flow, not shear layers.","tokens_in":21412,"tokens_out":2516,"would_cite":false,"duration_ms":25716,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A weakly nonlinear expansion of the compressible Euler equations yields a Landau-Stuart amplitude equation whose Mach-dependent coefficients switch between saturating and destabilizing regimes for shear-layer instability.","keywords":["weakly nonlinear stability","Landau-Stuart equation","Kelvin-Helmholtz instability","compressible shear layer","Mach number","Hopf bifurcation","multiple-scale expansion"],"falsifier":"Solve the linearized compressible Euler equations retaining the $v_1\\,d\\bar{u}/dy$ term for a hyperbolic-tangent profile $\\bar{u}(y)=\\tanh(y)$ and check the imaginary part of $\\omega$; if no wavenumber $k$ gives $\\Im(\\omega)\\neq 0$, the growth-rate maps in the paper are not Kelvin-Helmholtz growth rates. A complementary check is a direct numerical simulation of a compressible mixing layer at $M=1.1$ versus $M=1.3$: the first should saturate to a finite-amplitude limit cycle, while the second should show runaway amplitude growth.","tokens_in":20280,"feed_emoji":"🌀","tokens_out":10541,"duration_ms":106572,"temperature":0.7,"pith_summary":"This paper tries to show that the nonlinear fate of a compressible shear layer can be read off from a one-equation model: the Landau-Stuart amplitude equation, with Mach number entering through the coefficients. Using a multiple-scales expansion of the compressible Euler equations, it derives that equation and computes how the nonlinear Landau coefficient changes sign as the Mach number increases. The central prediction is alternating supercritical and subcritical Hopf bifurcations, with regimes where disturbances settle into finite-amplitude oscillations alternating with regimes where they blow up, separated by five critical Mach numbers. If correct, the analysis turns a hard fluid-dynamical stability problem into a tractable dynamical-systems problem, giving concrete predictions about saturation levels, limit cycles, and transition thresholds in supersonic shear flows. That matters for predicting unsteadiness in high-speed aerodynamics and astrophysical jets.","feed_headline":"Mach number decides if shear-layer disturbances saturate or blow up","feed_subtitle":"Mach number alternately stabilizes and destabilizes compressible shear layers: five critical thresholds.","key_machinery":"The central object is the Landau-Stuart equation $dA/dT = \\mu A - \\zeta |A|^2 A$ for the complex amplitude $A$ on the slow time scale $T = \\varepsilon^2 t$. The coefficient $\\mu$ is the linear growth rate and $\\zeta$ is the complex Landau coefficient; the real part of $\\zeta$ decides whether nonlinearity saturates the disturbance (supercritical case) or amplifies it (subcritical case), while the imaginary part controls phase rotation. The equation is obtained from a multiple-scales expansion of the compressible Euler equations: first-order normal modes modulate in amplitude, quadratic self-interaction generates second harmonics and mean-flow corrections, and a solvability condition at third order projects the resonant nonlinear terms onto the adjoint eigenfunction. The Mach number enters through every coefficient, which is what lets the paper read off bifurcation type as a function of $M$.","core_discovery":"The paper claims that a weakly nonlinear multiple-scales expansion of the two-dimensional compressible Euler equations, with a slow-time amplitude modulation, yields the Landau-Stuart equation $dA/dT = \\mu A - \\zeta |A|^2 A$ for the complex perturbation amplitude $A$. The linear growth rate $\\mu$ and the complex Landau coefficient $\\zeta$ inherit their dependence on Mach number from the first- and second-order perturbation equations. The substance of the claim is that $\\Re(\\zeta)$ changes sign repeatedly as $M$ increases: negative values give subcritical Hopf bifurcations where finite-amplitude disturbances explode, while positive values give supercritical Hopf bifurcations where disturbances saturate into limit cycles. The paper reports five critical Mach numbers $M_c = 0.2564, 1.0283, 1.2755, 1.3313, 1.5484$ at which the nonlinear stability character switches, and it interprets phase portraits and trajectories as evidence of saturation, spiral decay, and blow-up in the complex amplitude plane. It further claims that the bounded-domain boundary conditions used here lower perturbation amplitudes and confine pressure fluctuations to the shear-layer center, altering the Mach-number dependence relative to the classical unbounded analysis.","pith_inferences":["Restoring the omitted $v_1\\,d\\bar{u}/dy$ term in the first-order equations and recomputing the Landau coefficient would be the most direct test of whether the five critical Mach numbers survive; the alternating bifurcation pattern may change substantially.","Because the coefficient $\\zeta$ is computed at a single wavenumber, a natural extension is to allow spatial modulation of $A$, yielding a complex Ginzburg-Landau equation that would connect the present amplitude dynamics to the streamwise growth of mixing layers.","The same multiple-scales machinery should transfer to stratified or magnetized shear layers, where the control parameter would be a Richardson or Alfvén Mach number instead of $M$, giving cheap predictions for astrophysical jet stability and filament formation."],"forward_implications":["At Mach numbers where $\\Re(\\zeta)>0$, linearly unstable disturbances saturate into finite-amplitude periodic states (supercritical Hopf bifurcation), so the long-time state is a stable vortex street rather than unbounded growth.","At Mach numbers where $\\Re(\\zeta)<0$, the same disturbances blow up without saturation (subcritical Hopf bifurcation), meaning finite-amplitude perturbations can destabilize a linearly stable flow and generate hysteresis.","Five critical Mach numbers, 0.2564, 1.0283, 1.2755, 1.3313, and 1.5484, mark alternating stable and unstable nonlinear regimes in the compressible shear layer.","The finite-domain boundary conditions reduce perturbation amplitudes and localize pressure fluctuations near the layer center, so confinement acts as a stabilizing mechanism.","Linear stability analysis alone is insufficient: the sign of the real part of the Landau coefficient, not just the linear growth rate, determines whether the flow transitions smoothly or abruptly."],"supporting_citations":[{"why":"Supplies the linear inviscid compressible shear-layer framework that the weakly nonlinear expansion extends.","marker":"[13]"},{"why":"Adds complex flow profiles and density stratification to the linear dispersion relation, providing the multi-mode baseline for the nonlinear study.","marker":"[14]"},{"why":"Gives the high-Mach-number asymptotic behavior of the linear modes, informing the regime in which the Landau coefficient is evaluated.","marker":"[15]"},{"why":"Provides an earlier second-order nonlinear spatial stability analysis of compressible mixing layers, serving as a baseline for the weakly nonlinear framework.","marker":"[10]"},{"why":"Demonstrates weakly nonlinear amplitude equations for self-sustained detonations, the methodological analogue for deriving Landau-type saturation.","marker":"[9]"}],"fun_headline_variants":["Mach number flips shear layers between saturation and blow-up","Five critical Mach numbers switch compressible shear stability","Supersonic shear: Mach number alternates stable and explosive regimes","Nonlinear compressible shear: Mach number controls ultimate fate"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The first-order perturbation equations drop the term $v_1\\,d\\bar{u}/dy$ from the momentum balance, so the background velocity is effectively treated as uniform; without that term, the equations cannot represent the velocity-gradient mechanism that drives Kelvin-Helmholtz instability.","fun_headline_variants_meta":{"raw":{"variants":["Mach number flips shear layers between saturation and blow-up","Five critical Mach numbers switch compressible shear stability","Supersonic shear: Mach number alternates stable and explosive regimes","Nonlinear compressible shear: Mach number controls ultimate fate"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000719,"raw_usage":{"total_tokens":3302,"prompt_tokens":1093,"completion_tokens":2209,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":2142}},"tokens_in":709,"tokens_out":2209,"duration_ms":18650,"temperature":1.0,"reasoning_tokens":2142,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:25:24.765574+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the linearized compressible Euler equations retaining the $v_1\\,d\\bar{u}/dy$ term for a hyperbolic-tangent profile $\\bar{u}(y)=\\tanh(y)$ and check the imaginary part of $\\omega$; if no wavenumber $k$ gives $\\Im(\\omega)\\neq 0$, the growth-rate maps in the paper are not Kelvin-Helmholtz growth rates. A complementary check is a direct numerical simulation of a compressible mixing layer at $M=1.1$ versus $M=1.3$: the first should saturate to a finite-amplitude limit cycle, while the second should show runaway amplitude growth.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the linear inviscid compressible shear-layer framework that the weakly nonlinear expansion extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Adds complex flow profiles and density stratification to the linear dispersion relation, providing the multi-mode baseline for the nonlinear study."},{"cited_title":"Ladeinde \\ and\\ author J","cited_arxiv_id":null,"evidence_quote":"Gives the high-Mach-number asymptotic behavior of the linear modes, informing the regime in which the Landau coefficient is evaluated."},{"cited_title":"\\ Chaturvedi , author S","cited_arxiv_id":null,"evidence_quote":"Demonstrates weakly nonlinear amplitude equations for self-sustained detonations, the methodological analogue for deriving Landau-type saturation."}],"review_version":1}