{"id":"02e03047-813e-4ba4-a24f-b5b2e4e27df6","arxiv_id":"2501.16268","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A rigorous proof that subsonic boundary layer profiles are structurally stable for the steady compressible Navier-Stokes equations, plus a low Mach number limit with a Prandtl layer.","lead":"This paper proves that steady shear-flow boundary layers for the 2D compressible Navier-Stokes equations are stable for the entire subsonic Mach number range, in a Sobolev setting with small viscosity. It also establishes a low Mach number limit that keeps the Prandtl boundary layer visible at leading order.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The theorem's m in (0,1) assumption does not imply uniform positivity of A=1-m^2 U_s^2; the missing upper bound on U_s in (1.3)-(1.4) is a load-bearing gap.","rationale":"The central claim of Theorem 1.1 is that every shear profile satisfying (1.3)-(1.4) is structurally stable for every Mach number m in (0,1). For this to be true, the coefficient A=1-m^2 U_s^2 must remain uniformly positive, because the quasi-compressible reduction (3.3)-(3.5) divides by A and defines the elliptic operators used throughout the linear analysis. The paper asserts this lower bound follows from m<1, but it only follows if sup U_s <= 1/m; with m arbitrarily close to 1, even a small overshoot above 1 destroys it. Since (1.3) only fixes U_s'(0)=1, positivity, and the limit 1, overshooting profiles are admissible. The reader's weakest_assumption identifies exactly this gap. My independent reading confirms it: Section 3.1 uses A >= 1-m^2, and Section 6.1 explicitly invokes 0 <= U_s <= 1 in the Stokes absorption estimate, neither of which is present in the hypotheses. The intended class is almost certainly the standard monotone Prandtl layer with 0 <= U_s <= 1 or, more generally, m sup U_s < 1, so the proof should be repairable by adding such a condition. Because the gap is in the theorem statement rather than in the core strategy, and the proof is substantial and coherent modulo this hypothesis, the reader's conditional verdict is appropriate and no further adjustment is needed.","tokens_in":92266,"tokens_out":5761,"duration_ms":57680,"concrete_test":"Take the explicit admissible profile U_s(Y)=1-exp(-Y)+eps Y^2 exp(-Y) with eps=0.1 and m=0.99, and compute the minimum of A(Y)=1-m^2 U_s(Y)^2 over Y>=0. If the minimum is negative, then A^{-1} in (3.3) and Lambda in (3.5) are not uniformly positive, so the Rayleigh-Airy estimates in Lemmas 3.1-3.3 cannot hold for this profile; this directly contradicts the assertion in Section 3.1 and confirms that an extra bound on U_s must be added to Theorem 1.1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 3.1 the proof defines A(Y)=1-m^2 U_s(Y)^2 and asserts that m in (0,1) gives A(Y) >= 1-m^2 > 0. This is false under the stated assumptions (1.3)-(1.4), which only impose U_s>0, U_s(0)=0, U_s'(0)=1, U_s -> 1, and algebraic decay; U_s may overshoot 1. For example, U_s(Y)=1-exp(-Y)+eps Y^2 exp(-Y) satisfies (1.3)-(1.4) but has sup U_s>1, so for m sufficiently close to 1, A(Y)<0 on an interval. The coefficient A^{-1} enters the density formula (3.3), the modified vorticity operator Lambda in (3.5), and every Rayleigh-Airy estimate in Lemmas 3.1-3.3 and Proposition 3.4; ellipticity is lost exactly where the iteration needs a uniform positive lower bound. Moreover, Section 6.1 explicitly uses the stronger fact 0 <= U_s <= 1 to close the Stokes estimate around (6.33)-(6.34), and that inequality is not among (1.3)-(1.4). Thus Theorem 1.1 as stated overreaches: the proof requires an additional hypothesis such as m sup U_s < 1 or U_s <= 1. This is repairable, and the intended class is almost certainly the standard monotone Prandtl layer, but the stated theorem is not proven for its full hypothesis class.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the two-dimensional steady compressible Navier-Stokes equations in the half-plane, with the horizontal variable on a torus, near a shear boundary-layer profile (1, U_s(Y), 0). Theorem 1.1 claims, for every Mach number m in (0,1) and under the structural conditions (1.3)-(1.4), that small external-force perturbations admit a unique solution with bounds in a weighted high-regularity space and a zero-mass condition. Theorem 1.2 derives a low-Mach-number limit to the steady incompressible Navier-Stokes equations with a Prandtl layer. The proof proceeds by a zero-mode analysis, a Fourier-mode decomposition, a quasi-compressible Orr-Sommerfeld approximation, a quasi-compressible-Stokes iteration, boundary-layer correctors in low/middle/high frequency regimes, and a nonlinear iteration in Sobolev spaces. The paper is long and technically substantial, and the proof strategy is original in combining the authors' earlier quasi-compressible-Stokes framework with the Rayleigh-Airy machinery of Gérard-Varet and Maekawa.","tokens_in":92521,"tokens_out":8537,"duration_ms":85333,"significance":"If the missing hypothesis identified below is added, this would be a significant contribution: it provides uniform-in-Mach-number structural stability of subsonic steady Prandtl-type boundary layers in Sobolev spaces, over the whole subsonic range m in (0,1), and it gives a low-Mach-number limit in the presence of a Prandtl boundary layer. The proof is coherent and the frequency-by-frequency analysis is a genuine technical achievement. The paper does not assume its own conclusion; it builds on independent published results [14,38]. However, the theorem as stated overreaches because a key ellipticity assumption is not implied by (1.3)-(1.4). The gap is local and repairable, but it is load-bearing for every subsequent estimate.","major_comments":[{"comment":"The paper asserts that m in (0,1) implies A(Y)=1-m^2 U_s(Y)^2 satisfies A(Y) >= 1-m^2 > 0. This is false under assumptions (1.3)-(1.4), which impose U_s>0, U_s(0)=0, U_s'(0)=1, U_s->1 as Y->infinity, and algebraic decay, but do not prevent U_s from exceeding 1 or even 1/m. If U_s overshoots, A can vanish or change sign. Since A^{-1} enters the density formula (3.3), the modified vorticity operator Lambda in (3.5), and all subsequent Rayleigh-Airy estimates in Lemmas 3.1-3.3 and Proposition 3.4, the loss of a uniform positive lower bound for A invalidates the claimed solvability of the compressible Orr-Sommerfeld equation. The theorem should either add the hypothesis 0 < U_s <= 1 (or, more generally, m sup U_s < 1) or prove such a bound from the existing assumptions; the latter is not possible as stated.","section":"Section 3.1, Eq. (3.3)-(3.5)"},{"comment":"The closing estimate for the Stokes regularization system explicitly uses the inequality 0 <= U_s <= 1. Specifically, after Eq. (6.33), the proof absorbs the term involving ||sqrt(U_s) rho|| using m<1 and 0<=U_s<=1. This inequality is not among the structural conditions (1.3)-(1.4); it is an unstated additional hypothesis. Consequently Proposition 6.3, the linear stability theorem Theorem 7.5, and the nonlinear Theorem 1.1 all rely on an assumption that is not part of the theorem statement. The same issue appears in Section 4 around Eq. (4.11), where the proof uses 'A^{-1} ~ 1' without a uniform upper and lower bound on U_s. The repair is the same: add a uniform upper bound on U_s, for example 0 < U_s <= 1, or equivalently m sup U_s < 1, to the hypotheses of Theorems 1.1 and 1.2.","section":"Section 6.1, Eqs. (6.33)-(6.34)"},{"comment":"Several technical lemmas are stated and then deferred to prior work, especially [14, Propositions 5.1, 6.1, 7.11] and [38]. Some of these are classical incompressible results, but here they are applied to compressible operators containing A(Y) and m-dependent coefficients. Because the missing U_s bound affects the validity of these cited results, the authors should either give proofs of these lemmas or provide a precise verification that the hypotheses of the cited propositions are satisfied under the repaired assumptions. This is important for verification but is secondary to the main hypothesis gap.","section":"Lemmas 3.1, 3.2, 5.1, 5.4"}],"minor_comments":[{"comment":"There are numerous typographical and rendering errors, such as 'Foturnately', 'middfle', 'asscciated', 'suffices', and stray '/greaterorsimilar' symbols in Section 3.4 and elsewhere. These should be cleaned up.","section":"Throughout"},{"comment":"The remark says L can be large if the amplitude of the boundary-layer profile is small. Since the proof repeatedly requires L to be small, this remark should be made precise or moved to a heuristic comment, otherwise it may confuse the reader about the quantitative role of L.","section":"Remark 1.1(v)"},{"comment":"The phrase 'the first result concerning the low Mach number limit in the presence of Prandtl boundary layers' is a strong claim. Please check the existing literature carefully and, if this is indeed the first steady result, state the comparison class explicitly.","section":"Theorem 1.2"}],"recommendation":"major_revision","confidential_remarks":"The central mathematical architecture is sound and the paper is likely acceptable after the authors add the missing U_s <= 1 (or m sup U_s < 1) hypothesis and adjust the statements of Theorems 1.1 and 1.2 accordingly. I do not see circularity or novelty concerns. The main work for the revision is to repair the theorem statement and to verify the deferred technical lemmas under the repaired assumptions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead the whole thing carefully. The headline is that this is a real step forward: the first nonlinear Sobolev-space structural stability theorem for compressible Prandtl-type boundary layers across the full subsonic range, plus a low-Mach-number limit that keeps the Prandtl layer. That is genuinely new relative to what I know: previous compressible results either needed analytic regularity, covered only m<1/sqrt(3), or were linear instability analyses. The quasi-compressible-Stokes iteration from [38] is extended here to the steady problem with a clean frequency decomposition, and the higher-order nonlinear estimates in Section 8 are the kind of laborious, substantive work that makes the paper worth referee time. The zero-mode analysis and the mass constraint are handled properly.\n\nNow the soft spots, in order of seriousness.\n\nThe load-bearing gap is real and the stress-test note is right. In Section 3.1, A(Y)=1-m^2 U_s(Y)^2 is asserted to satisfy A>=1-m^2 from m in (0,1) alone. That requires U_s(Y)<=1. Conditions (1.3)-(1.4) only give U_s>0, U_s(0)=0, U_s'(0)=1, decay to 1, and algebraic decay of derivatives; they do not rule out overshoot above 1. For m close to 1, A can vanish or change sign, and then the Orr-Sommerfeld operator loses ellipticity at exactly the point where every Rayleigh-Airy estimate needs a uniform positive lower bound. The paper even uses 0<=U_s<=1 explicitly in Section 6.1 around (6.33), so the proof already assumes more than the theorem states. This is repairable: add the hypothesis m sup U_s < 1 or U_s<=1, and the intended class is almost certainly the standard monotone Prandtl layer. But Theorem 1.1 as written overreaches its hypotheses.\n\nMinor but worth noting: several lemmas (3.1, 3.2, 5.1, 5.4) defer proofs to [14] or prior work. That is not a flaw by itself, but the verification burden is nontrivial because the compressible coefficients change the operators. The high-frequency regime is solved by energy methods and is less delicate, so the gap is isolated in the low/middle frequency analysis.\n\nWho is this for: anyone working on boundary layer stability, compressible Navier-Stokes asymptotics, or low-Mach limits. A serious referee should engage with it, but should require the missing hypothesis and a check that every use of A^{-1} carries the corrected assumption.\n\nRecommendation: send to peer review, conditional on the authors fixing the U_s bound and restating Theorems 1.1-1.2 accordingly. The core machinery holds up; the stated theorem needs tightening.","headline":"A substantial nonlinear stability result for subsonic compressible boundary layers, likely correct in intent but stated with a repairable gap: the uniform subsonic condition needs U_s <= 1, which is not in (1.3)-(1.4).","tokens_in":93161,"tokens_out":837,"would_cite":false,"duration_ms":13715,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76N20","76N10","35Q35","35B35"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper establishes uniform-in-Mach-number structural stability of steady compressible boundary layers in the entire subsonic regime, and derives the first low-Mach Prandtl-layer limit.","keywords":["subsonic compressible Navier-Stokes","boundary layer stability","Prandtl boundary layer","Orr-Sommerfeld equation","quasi-compressible-Stokes iteration","low Mach number limit","shear flow","Rayleigh-Airy iteration"],"falsifier":"Let $m=0.95$ and $U_s(Y)=1-e^{-Y}+\\frac12 Y^2e^{-Y}\\sin(10Y)$, which satisfies $U_s(0)=0$, $U_s'(0)=1$, positivity on $\\mathbb R_+$, and the algebraic decay (1.4). Near $Y=2$ this profile exceeds $1/m\\approx1.053$, so $A(Y)<0$ on an interval; the Rayleigh and Airy estimates in Section 3, which all use $A^{-1}$ as an elliptic weight, then fail. Evaluating $A$ on any admissible profile with $\\sup U_s>1/m$ is therefore a direct check of whether Theorem 1.1 needs an extra hypothesis.","tokens_in":91960,"feed_emoji":"🧮","tokens_out":9112,"duration_ms":85999,"temperature":0.7,"pith_summary":"This paper proves that a compressible shear-flow boundary layer of the form $(\\rho_s,u_s)=(1,U_s(y/\\sqrt\\nu),0)$ is structurally stable for the full subsonic range of Mach numbers $m\\in(0,1)$, provided the external-force perturbation is small with respect to $\\nu^{9/8+}$ and the spatial period $L$ is short. Such stability was previously established only in restricted settings, and the incompressible steady boundary-layer theory is recovered as the $m\\to0$ limit. The authors establish existence, uniqueness, and uniform-in-$m$ estimates for the steady two-dimensional compressible Navier-Stokes equations in a half-plane with periodic tangential direction, with a boundary-layer profile attached to a no-slip wall. As a byproduct they prove the first low-Mach-number limit in the presence of a Prandtl boundary layer, at rate $O(m^2\\nu^{-5/8-})$. If correct, the result gives a Sobolev-level validation of Prandtl's boundary-layer expansion for steady subsonic compressible flows.","feed_headline":"Subsonic boundary layers hold for every Mach number below 1","feed_subtitle":"Steady compressible shear layers are stable uniformly in the Mach number, with the first low-Mach Prandtl-layer limit.","key_machinery":"The load-bearing device is the quasi-compressible approximation (3.1): an artificial viscosity $\\sqrt\\nu\\Delta_\\alpha$ acting separately on the velocity components replaces the physical viscosity, which breaks the density-velocity coupling. Writing $u=\\partial_Y\\varphi-U_s\\varrho$, $v=-i\\alpha\\varphi$ for an effective stream function $\\varphi$ reduces the approximation to the compressible Orr-Sommerfeld equation $i\\epsilon\\Lambda(\\Delta_\\alpha\\varphi)+U_s\\Lambda(\\varphi)-\\varphi\\partial_Y(A^{-1}\\partial_Y U_s)=f$, where $\\Lambda(\\varphi)=\\partial_Y(A^{-1}\\partial_Y\\varphi)-\\alpha^2\\varphi$ and $A(Y)=1-m^2U_s(Y)^2$. The factor $A$, treated as uniformly positive because $m<1$, is what makes $\\Lambda$ elliptic. Rayleigh-Airy iteration solves the equation at low and middle frequencies; a direct energy method covers high frequencies. A Stokes regularizing system and boundary-layer correctors built from homogeneous quasi-compressible solutions then recover the no-slip velocity condition, and a modified linear system controls the nonlinear iteration.","core_discovery":"This paper's central claim is Theorem 1.1: for every Mach number $m\\in(0,1)$ and any shear profile $U_s$ satisfying the structural conditions (1.3)-(1.4), there exists a torus length $L_0$ such that for $L\\in(0,L_0)$, small $\\nu$, and external perturbations of weighted norm at most $\\nu^{9/8+}$, the steady compressible Navier-Stokes system (1.6) has a unique solution $(\\rho,u,v)$ near the boundary layer $(1,U_s(y/\\sqrt\\nu),0)$, obeying $\\|(\\rho,u,v)\\|_X\\le C\\|(F_{\\mathrm{ext},1},F_{\\mathrm{ext},2})\\|_w$ together with the zero-mass condition. Theorem 1.2 then sends $m\\to0$: the same family of solutions converges to the incompressible Navier-Stokes solution with Prandtl boundary layer, at rate $O(m^2\\nu^{-5/8-})$, which the authors identify as the first such low-Mach limit in the presence of a Prandtl layer. The argument splits the linearized problem into zero and non-zero Fourier modes, solves a compressible Orr-Sommerfeld equation by Rayleigh-Airy iteration in low and middle frequencies and by an energy method in high frequencies, and uses boundary-layer corrections to restore the no-slip condition.","pith_inferences":["Beyond the paper, the same proof structure implies the theorem should be stated with the extra structural condition $\\sup_Y U_s(Y)<1/m$; without it the claimed uniform positivity of $A=1-m^2U_s^2$ is not guaranteed by $m\\in(0,1)$ alone.","A natural testable extension is the non-isentropic case, where the Mach factor becomes $1-m^2U_s^2$ times a thermodynamic function of the background profile; the quasi-compressible iteration should still close at all subsonic Mach numbers.","The authors' remark that the torus length can be large when the shear amplitude is small points to a long-torus variant: fix any period and make $\\|\\partial_Y U_s\\|_{L^\\infty}$ small, connecting structural stability to spectral conditions rather than to shortness of the spatial period.","The reported rate $O(m^2\\nu^{-5/8-})$ likely reflects the norm's worst derivative weights; a refined cancellation between the pressure and divergence terms could lift the negative power of $\\nu$."],"forward_implications":["For every $m\\in(0,1)$, a sufficiently short torus supports a unique steady solution near the shear boundary layer, with error controlled linearly by the weighted external force.","The uniform-in-$m$ estimates make the entire subsonic regime a single parameter range: no separate treatment is needed at any $m<1$.","Letting $m\\to0$ recovers the incompressible Navier-Stokes solution with a Prandtl layer, and the boundary layer survives the low-Mach limit at an explicit rate.","The zero-mass condition and the extra boundary condition $\\operatorname{div}_{x,y}(u,v)|_{y=0}=0$ are preserved by the nonlinear solution, which underpins the higher-order estimates.","The stability holds in Sobolev regularity rather than analytic regularity, so the result applies to nonzero modes across all frequency ranges."],"supporting_citations":[{"why":"supplies the quasi-compressible-Stokes iteration and the linear instability framework on which Sections 3-6 are built.","marker":"[38]"},{"why":"provides the steady Rayleigh-Airy iteration and slow/fast mode construction that the low-frequency analysis adapts to compressible flow.","marker":"[14]"},{"why":"gives the shear-flow stability framework and the high-frequency boundary-layer profile used for the middle- and high-frequency correctors.","marker":"[5]"},{"why":"motivates the modified linear system (8.1) that absorbs the derivative loss in the nonlinear continuity equation.","marker":"[24]"},{"why":"supplies the penalty-term trick that enforces the zero-mass condition for the zero Fourier mode.","marker":"[11]"},{"why":"identifies Tollmien-Schlichting destabilization across the subsonic regime, the instability that the present stability theorem must overcome.","marker":"[28]"}],"fun_headline_variants":["Entire subsonic regime: boundary-layer stability proven","First low-Mach Prandtl-layer limit comes with full stability","All Mach numbers below 1: shear layers structurally stable"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"For every estimate in the proof, $A(Y)=1-m^2U_s(Y)^2$ must stay uniformly positive, but conditions (1.3)-(1.4) do not force $U_s\\le 1$; a profile that rises above $1/m$ makes $A$ change sign and destroys the ellipticity on which the Orr-Sommerfeld analysis rests.","fun_headline_variants_meta":{"raw":{"variants":["Entire subsonic regime: boundary-layer stability proven","First low-Mach Prandtl-layer limit comes with full stability","All Mach numbers below 1: shear layers structurally stable"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000644,"raw_usage":{"total_tokens":2984,"prompt_tokens":994,"completion_tokens":1990,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":610,"completion_tokens_details":{"reasoning_tokens":1935}},"tokens_in":610,"tokens_out":1990,"duration_ms":14334,"temperature":1.0,"reasoning_tokens":1935,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:35:04.975736+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Let $m=0.95$ and $U_s(Y)=1-e^{-Y}+\\frac12 Y^2e^{-Y}\\sin(10Y)$, which satisfies $U_s(0)=0$, $U_s'(0)=1$, positivity on $\\mathbb R_+$, and the algebraic decay (1.4). Near $Y=2$ this profile exceeds $1/m\\approx1.053$, so $A(Y)<0$ on an interval; the Rayleigh and Airy estimates in Section 3, which all use $A^{-1}$ as an elliptic weight, then fail. Evaluating $A$ on any admissible profile with $\\sup U_s>1/m$ is therefore a direct check of whether Theorem 1.1 needs an extra hypothesis.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the quasi-compressible-Stokes iteration and the linear instability framework on which Sections 3-6 are built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the steady Rayleigh-Airy iteration and slow/fast mode construction that the low-frequency analysis adapts to compressible flow."},{"cited_title":"F.: On the stability of shear flow s of Prandtl type for the steady Navier-Stokes equations, Sc i","cited_arxiv_id":null,"evidence_quote":"gives the shear-flow stability framework and the high-frequency boundary-layer profile used for the middle- and high-frequency correctors."},{"cited_title":"Rational Mech","cited_arxiv_id":null,"evidence_quote":"motivates the modified linear system (8.1) that absorbs the derivative loss in the nonlinear continuity equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the penalty-term trick that enforces the zero-mass condition for the zero Fourier mode."},{"cited_title":"F.: Tollmien-S chlichting waves in the subsonic regime, Proc","cited_arxiv_id":null,"evidence_quote":"identifies Tollmien-Schlichting destabilization across the subsonic regime, the instability that the present stability theorem must overcome."}],"review_version":1}