{"id":"17ddc94d-e2c9-4670-9f62-bb2a855a82b8","arxiv_id":"1908.04945","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Axisymmetric transonic shocks with nonzero angular momentum in divergent nozzles are shown to be stable under small perturbations of the incoming supersonic flow and exit pressure.","lead":"This paper proves that three-dimensional axisymmetric supersonic flow through a widening nozzle can form a stable shock when the exit pressure is slightly perturbed, including flows that swirl. The result extends an earlier zero-swirl proof and gives a new stream-function method to locate the shock.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 2.17 is proved only with B fixed at B0 on the entrance; Remark 2.13's general-perturbation claim is unsupported, so the advertised stability is narrower than the proof.","rationale":"The main construction has substantial independent support: after the reduction to B = B0, the paper supplies a full fixed-point scheme, an elliptic-system treatment of the singular stream-function equation, and a Fredholm/invertibility argument for the shock-shape map. I found no algebraic inconsistency that would invalidate Theorem 2.16 or Theorem 2.17 for data satisfying (2.3.2). The soft spot is the scope: Problem 1 fixes B = B0 on Γen, and §2.5 makes that restriction structural by reducing away the B-equation and the [B] Rankine-Hugoniot condition. All subsequent steps operate on the reduced system. Remark 2.13 is the only claim of generality, and it is a bare assertion. Since the abstract advertises small perturbations of the incoming supersonic flow without the B caveat, this is a load-bearing gap for the advertised claim, though not for the theorem as stated. The suggested check isolates whether the restriction is cosmetic or essential. The reader's weakest assumption identified this same spot, and my recommendation does not move the CONDITIONAL verdict.","tokens_in":91426,"tokens_out":12257,"duration_ms":143609,"concrete_test":"Perform the following analytical check: write the full stream-function system (2.4.15) with B = B0 + εb, keep the Rankine-Hugoniot condition [B]Γ = 0 from (2.4.20), and linearize at the radial shock for a representative b(θ) = sin(πθ/θ1). Compare the resulting equations with the reduced linearization (3.1.7)–(3.1.14). If the εb-terms enter F1 or the boundary formula (3.1.11) at first order, then the existing quadratic estimate (3.1.16) cannot absorb them and the fixed-point/implicit-function argument in §3.4–§4.1 must be reworked. If the terms are identically zero or lower order, Remark 2.13 can be turned into a proof by adding the corresponding estimates.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Problem 1 fixes B = B0 on Γen, equation (2.3.2), and §2.5 uses this to reduce the Euler system and the Rankine-Hugoniot conditions: the full R-H set (2.4.20) contains [B]Γ = 0, but the reduced Problem 2 instead imposes B = B0 throughout N+f. Every later ingredient—Problem 3, the linearized system (3.1.7)–(3.1.14), the Fréchet derivative in Lemma 4.2, and the implicit-function argument in §4.1—depends on that reduction. Remark 2.13 asserts that the result is unchanged for a general perturbation of (ρ−0, u−0er, p−0), but it supplies no transport equation for B − B0, no treatment of [B]Γ = 0, and no fixed-point reformulation. This is not a consistency error in Theorem 2.17 as written, since the theorem does contain (2.3.2). It does, however, make the abstract's 'small perturbations of an incoming radial supersonic flow' and the advertised generality narrower than what is proved. If Remark 2.13 is false, Theorem 2.17 covers only a codimension-one subclass of small incoming data; if it is true, the proof is missing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies 3-D axisymmetric transonic shock solutions of the steady full Euler system in divergent nozzles. The main result, Theorem 2.17 (via Theorem 2.16 for a stream-function formulation), asserts that a given radial transonic shock solution is structurally stable under sufficiently small axisymmetric perturbations of the incoming supersonic radial flow and of the exit pressure, including perturbations with nonzero angular momentum, with no restriction on the nozzle tip angle and without an S-condition. The proof recasts the downstream subsonic problem in a vector-potential stream-function formulation, splits the free boundary into an initial shock position and a shock shape, solves a Pseudo Free Boundary Problem by a Schauder fixed point argument after linearization, and then determines the shock shape by a weak implicit function theorem. The paper also resolves the axis singularity of the stream-function formulation by solving the singular elliptic equation as an elliptic system and proves an orthogonal completeness statement for an associated Legendre problem.","tokens_in":91712,"tokens_out":3069,"duration_ms":35879,"significance":"If the main theorem is correct, this is a significant contribution to the transonic-shock stability literature: it extends prior 3-D axisymmetric results to nonzero swirl, removes small-tip-angle restrictions, and avoids the S-condition used in some earlier general 3-D work. The manuscript contains substantial original technical machinery: the vector-potential stream-function reformulation, the monotonicity argument (Lemma 2.10) that fixes the initial shock position through the exit-pressure solvability condition, and an implicit-function scheme for the shock shape that avoids the nonlocal elliptic equation appearing in earlier iteration schemes. These are real strengths and make the paper worth serious consideration. However, the advertised range of admissible incoming perturbations is narrower than what is proved unless Remark 2.13 is supplied with a proof, and several regularity statements that are load-bearing for the fixed-point argument are either omitted or delegated to previous work.","major_comments":[{"comment":"Problem 1 fixes the Bernoulli invariant exactly, B = B0 on Γen, and §2.5 uses this assumption to impose B = B0 throughout the downstream region and to reduce the Rankine-Hugoniot conditions, replacing [B]Γ = 0 in (2.4.20) by B = B0 in N+f. Every subsequent object—Problem 2, Problem 3, the linearized system (3.1.7)–(3.1.14), the Fréchet derivative in Lemma 4.2, and the implicit-function argument in §4—depends on that reduction. Remark 2.13 asserts that the result is unchanged for a general perturbation of (ρ−0, u−0er, p−0), but no transport equation for B−B0, no treatment of the Rankine-Hugoniot condition [B]Γ = 0, and no modified fixed-point formulation are provided. As written, Theorem 2.17 covers only the codimension-one class of entrance data satisfying (2.3.2), which is weaker than the abstract's claim of stability under small perturbations of an incoming radial supersonic flow. This is a load-bearing gap: please either prove Remark 2.13 or revise the abstract and introduction to state the theorem with the restriction B = B0.","section":"§2.3–§2.5, Problem 1, Eq. (2.3.2) and Remark 2.13"},{"comment":"Corollary 3.12 is stated with the sentence 'The following Corollary is obtained from Lemma 3.11 in the same way that Corollary 3.11 is obtained from Lemma 3.10 in [17]. We omit the proof.' This corollary is then used directly in the proofs of Lemma 3.9 and Lemma 3.10 to obtain the Cβ and C1,α regularity up to the corners Γf ∩ Γ+w and Γ+w ∩ Γex. Lemma 3.13, which upgrades this to C^{2,α}_{(-1-α,Γ+w)}(N+f), is also stated without proof, with only a reference to a scaling argument from [2] and Theorem 5.21 of [16]. These regularity results are essential for the Schauder fixed-point argument in Proposition 3.1 and for the Fréchet differentiability in Lemma 4.2. The manuscript should either provide the full proofs, or give precise theorem statements in the cited works with a verification that all hypotheses (especially the corner-reflection condition and the structure of the boundary data) are met in this setting.","section":"§3.2, Corollary 3.12 and Lemma 3.13"},{"comment":"The existence proof of Proposition 3.1 constructs a Schauder map on the compact convex set P(M1) and then states that continuity of the map J follows by 'the standard argument.' Given that J involves solving the transport equations (B′), the linear elliptic problem (3.4.25)–(3.4.26), and the intermediate-value step defining f(0), continuity is not entirely immediate and depends on the uniqueness estimates for Problems 3.1 and on the regularity of the elliptic solver. Since the entire fixed-point step rests on this continuity assertion, a short proof or a precise reference for the continuity of this composition should be supplied.","section":"§3.4, proof of Proposition 3.1, existence step"}],"minor_comments":[{"comment":"There are numerous typographical errors, including 'incomming' for 'incoming', 'Prolem' for 'Problem', 'Lemam' for 'Lemma', and 'funtion' for 'function'. These should be corrected in a revision.","section":"Throughout"},{"comment":"The abstract and the introduction state the result for 'small perturbations of an incoming radial supersonic flow' without mentioning the condition B = B0 on the entrance. Until Remark 2.13 is proved, the statement of the main theorem should explicitly include this restriction.","section":"Abstract and Introduction"},{"comment":"In the uniqueness part of Proposition 3.1, the estimates labelled (3.4.59) and (3.4.60) are cited before they are stated; renumbering or moving these displayed estimates would improve readability.","section":"§3.4, after (3.4.59)–(3.4.60)"},{"comment":"The uniqueness proof of Theorem 2.16 refers to 'the arguments in Step 1 in the proof of Lemma 4.2' for the unique solvability of (4.1.8)–(4.1.12) at low regularity; since that system is derived by a formal limiting process, a direct statement of the low-regularity well-posedness would help the reader verify the contraction argument.","section":"§4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper has already appeared in ARMA according to the acknowledgements, so the editor may wish to treat this report as an assessment of the arXiv version. The main concern is not the internal consistency of Theorem 2.17 as stated with (2.3.2), but the mismatch between the advertised generality and the proved theorem. This is fixable by either proving Remark 2.13 or by restating the abstract and introduction. I would not recommend rejection, as the central construction appears sound and the technical contributions are substantial."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a substantial paper and the central theorem is proved for the class it is actually stated on. The advertised generality is another matter: Remark 2.13 claims the result is unchanged for arbitrary small incoming perturbations, and that claim has no proof behind it.\n\nThe genuinely new pieces are the vector-potential stream-function formulation, which sidesteps the axis singularity by solving a singular elliptic equation as an elliptic system, and the two-step iteration that determines the initial shock position and the shock shape separately. The entropy monotonicity in Lemma 2.10 feeding an intermediate-value argument for the shock location is clean, and the invertibility argument for the linearized shock-shape operator (Lemma 4.4, via eigenfunction expansion and maximum principle) is serious. No free parameters are adjusted; the shock position comes out of the equations. That is exactly what a stability proof in this area should look like.\n\nThe soft spots are real, but they are mostly about scope rather than soundness of the core argument. The proof does not quite do what the abstract promises. Problem 1 fixes B = B0 on the entrance, and from §2.5 onward the entire reduction — the simplified Rankine-Hugoniot conditions, Problems 2 and 3, the linearized system, the Fréchet derivative, the implicit function argument — depends on B remaining equal to B0 throughout the downstream region. The abstract's \"small perturbations of an incoming radial supersonic flow\" therefore overstates the result: the theorem covers only incoming perturbations that preserve B on the entrance, which is a genuine infinite-dimensional restriction, not a byproduct of smallness. Remark 2.13 asserts the general case is unchanged, but supplies no transport equation, no treatment of [B]_Γ = 0, and no fixed-point reformulation. If the assertion is true, the proof is missing; if it is false, the theorem is narrower than advertised. The theorem itself is honest — (2.3.2) is there — but the packaging is not.\n\nTwo smaller items: Corollary 3.12 and Lemma 3.13 are deferred to prior work inside the regularity machinery, which is a noticeable hand-wave in a proof this long, though not unusual for this literature.\n\nFor whom: researchers in transonic shock stability, free-boundary problems for the full Euler system, or stream-function/axisymmetric methods. They will get a serious proof of a real missing case — nonzero swirl, no small tip-angle restriction — and they should know to read Theorem 2.17 as restricted to B = B0 on the entrance. The subsequent ARMA publication suggests referees found the structure acceptable. I would send this to peer review; the right referee will ask for Remark 2.13 to be proved or cut to match the proof.","headline":"A substantial stability theorem for 3-D axisymmetric transonic shocks with swirl, honestly proved for the restricted data class B = B0 on the entrance, but Remark 2.13 advertises a generality the proof does not deliver.","tokens_in":92171,"tokens_out":4300,"would_cite":true,"duration_ms":44781,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35A01","35A02","35J57","35J62","35M10","35Q31","35R35","76H05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Radial transonic shock solutions in divergent nozzles are stable under small axisymmetric perturbations of the incoming supersonic flow and the exit pressure, including perturbations with nonzero angular momentum.","keywords":["transonic shock","full Euler system","stream function","free boundary problem","axisymmetric flow","divergent nozzle","exit pressure","nonzero angular momentum"],"falsifier":"Evaluate the derivative in Lemma 2.10, $\\frac{d}{dt}\\bar S|_{D_t^+}(t;t)$, for a strongly divergent nozzle with $\\gamma$ close to $1$; a nonpositive value at some $t\\in(r_0,r_1)$ would break the shock-position adjustment mechanism. Alternatively, take a small entrance perturbation with $B\\ne B_0$ and try to close the fixed-point scheme; if no solution satisfying estimate (2.5.14) exists, the unproved claim in Remark 2.13 fails.","tokens_in":91169,"feed_emoji":"🌀","tokens_out":10930,"duration_ms":100696,"temperature":0.7,"pith_summary":"This paper establishes that radial transonic shocks in divergent nozzles are structurally stable: if the incoming supersonic flow and the exit pressure are changed by a small axisymmetric perturbation, there is a nearby shock surface and downstream subsonic solution of the steady full Euler equations, unique in a natural weighted Hölder class, and the perturbation may carry nonzero angular momentum. The proof works in a stream-function formulation that avoids the singularity at the symmetry axis, and it treats the shock's initial radius and its shape as two separate unknowns determined by two different mechanisms. This matters because previous full-Euler results either excluded angular momentum, required small nozzle opening angle, or imposed additional conditions to pin down the shock location.","feed_headline":"Transonic shocks in divergent nozzles are stable in 3-D","feed_subtitle":"Small perturbations of the incoming flow and exit pressure still give a unique shock surface, now with swirl allowed.","key_machinery":"The argument runs through a stream-function formulation of the axisymmetric full Euler system written with the vector-potential stream function $\\Phi e_\\phi$, the swirl variable $L=2\\pi r\\sin\\theta\\,u_\\phi$, the entropy $S$, and the Bernoulli invariant $B$. The divergence equation is built into the identity $\\mathrm{div}(\\nabla\\times(\\Phi e_\\phi))=0$, and the remaining equations become one singular elliptic equation for $\\Phi$ coupled to two transport equations for $L$ and $S$. The singularity at the axis is resolved by converting the singular scalar equation into an elliptic system of the form $\\mathrm{div}(A D(\\Psi e_\\phi))-d\\,\\Psi e_\\phi=-F$, whose coefficients are positive-definite because the downstream reference flow is subsonic; existence and regularity follow from Lax-Milgram, freezing of coefficients, and a reflection argument at the two perpendicular corners. The transport equations are solved by following level surfaces of the modified stream function $V=2\\pi r\\sin\\theta(\\Phi_0+\\Psi)$, giving solutions of the form $Q=Q_{\\mathrm{en}}(L)$. The shock location is split into two unknowns: the initial shock radius is fixed by a compatibility equation whose solvability rests on monotonicity of the downstream entropy in the shock position, and the shock shape is fixed by inverting the Fréchet derivative of the tangential Rankine-Hugoniot condition, using eigenfunctions of the associated Legendre problem of type $m=1$ on the interval $(0,\\theta_1)$.","core_discovery":"The central claim is Theorem 2.17: for any Hölder exponent $\\alpha\\in(2/3,1)$ there is a $\\sigma_2>0$ such that whenever the axisymmetric incoming supersonic flow differs from the fixed radial background by $\\sigma$ in $C^{2,\\alpha}$, with the Bernoulli invariant fixed to $B_0$ on the entrance, and the exit pressure differs from $p_c$ by $\\sigma$ in the appropriate weighted norm, Problem 1 has a unique axisymmetric transonic shock solution with shock front $r=f(\\theta)$ and downstream density, velocity and pressure satisfying estimate (2.5.14), so all deviations from the radial shock are controlled by $C\\sigma$. The proof first establishes the equivalent stream-function Problem 2, Theorem 2.16, and then recovers the physical variables through an implicit relation that determines the density from the stream-function unknowns. In the author's terms, this is structural stability of radial transonic shock solutions in divergent nozzles, with no restriction on the nozzle tip angle and no assumption that the incoming supersonic solution itself be radial.","pith_inferences":["Beyond the paper: if the assertion in Remark 2.13 can be upgraded to a proof, the theorem becomes full structural stability under arbitrary small axisymmetric entrance perturbations; currently the proven statement is conditional on the entrance Bernoulli invariant being fixed to $B_0$.","Beyond the paper: the two-step scheme, mass-flux compatibility for the initial shock radius and tangential Rankine-Hugoniot for the shock shape, suggests a numerical shock-fitting algorithm that alternates these two updates instead of solving the full free-boundary problem monolithically.","Beyond the paper: the completeness proof for the associated Legendre problem of type $m=1$ on a general interval is a technical tool that likely transfers to other axisymmetric free-boundary elliptic systems with conical corners.","Beyond the paper: the lower bound $\\alpha>2/3$ appears tied to the corner regularity obtained by the reflection argument; testing whether the Hölder exponent can be lowered with a different corner treatment would clarify the optimal regularity threshold."],"forward_implications":["For every sufficiently small $\\sigma$, Problem 1 has a solution and it is unique in the class satisfying (2.5.14): a nearby shock front and nearby subsonic flow exist, with all deviations controlled by $C\\sigma$.","The shock location is determined by the exit pressure alone, with no assumption that the shock passes through a prescribed wall point; the initial shock radius and the shock shape are fixed by two different conditions.","Perturbations with nonzero angular momentum are admitted: the swirl variable $L$ is transported along the streamlines, and the resulting velocity field remains regular across the shock.","The result covers nozzles of arbitrary opening angle $\\theta_1<\\pi$ and flows with $C^{1,\\alpha}$ interior and $C^\\alpha$ boundary regularity, so general axisymmetric nozzle walls are within reach."],"supporting_citations":[{"why":"Supplies the unique radial transonic shock solutions in divergent nozzles, the background states this paper perturbs.","marker":"[31]"},{"why":"Provides the non-isentropic potential model stability result and the weak implicit function theorem used in Section 4.","marker":"[3]"},{"why":"Introduces the two-mechanism view that initial shock position and shock shape are fixed by different conditions.","marker":"[24]"},{"why":"Gives the 2-D full Euler stability result with C^{1,alpha} interior regularity and general nozzle perturbations that the 3-D result extends.","marker":"[21]"},{"why":"Treats 3-D axisymmetric flows with zero angular momentum; this paper removes that restriction.","marker":"[20]"},{"why":"Supplies the regularity characterization of axisymmetric solenoidal vector fields used in Lemma 2.8.","marker":"[22]"},{"why":"Shows how to turn a singular axisymmetric elliptic equation into a solvable elliptic system; the same device is central here.","marker":"[4]"},{"why":"Provides the freezing-coefficients elliptic regularity theory used in Lemmas 3.9 and 3.10.","marker":"[15]"}],"fun_headline_variants":["3-D transonic shocks stable in divergent nozzles","Axisymmetric shock stability in divergent nozzles proven","Unique transonic shocks in divergent nozzles under perturbation","Divergent nozzle transonic shocks resist small perturbations","3-D axisymmetric shocks remain stable in divergent nozzles"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the incoming and downstream flows share the exact Bernoulli constant $B_0$ of the reference radial flow; the paper asserts in Remark 2.13 that this restriction can be dropped without changing the result, but supplies no proof, so if that assertion fails the theorem covers only entrance perturbations with $B=B_0$.","fun_headline_variants_meta":{"raw":{"variants":["3-D transonic shocks stable in divergent nozzles","Axisymmetric shock stability in divergent nozzles proven","Unique transonic shocks in divergent nozzles under perturbation","Divergent nozzle transonic shocks resist small perturbations","3-D axisymmetric shocks remain stable in divergent nozzles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000198,"raw_usage":{"total_tokens":1337,"prompt_tokens":884,"completion_tokens":453,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":500,"completion_tokens_details":{"reasoning_tokens":378}},"tokens_in":500,"tokens_out":453,"duration_ms":4493,"temperature":1.0,"reasoning_tokens":378,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:27:37.852663+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the derivative in Lemma 2.10, $\\frac{d}{dt}\\bar S|_{D_t^+}(t;t)$, for a strongly divergent nozzle with $\\gamma$ close to $1$; a nonpositive value at some $t\\in(r_0,r_1)$ would break the shock-position adjustment mechanism. Alternatively, take a small entrance perturbation with $B\\ne B_0$ and try to close the fixed-point scheme; if no solution satisfying estimate (2.5.14) exists, the unproved claim in Remark 2.13 fails.","supporting_citations":[{"cited_title":"Real World Appl","cited_arxiv_id":null,"evidence_quote":"Supplies the unique radial transonic shock solutions in divergent nozzles, the background states this paper perturbs."},{"cited_title":"Ra- tion","cited_arxiv_id":null,"evidence_quote":"Provides the non-isentropic potential model stability result and the weak implicit function theorem used in Section 4."},{"cited_title":"Hyperbolic Differ","cited_arxiv_id":null,"evidence_quote":"Introduces the two-mechanism view that initial shock position and shock shape are fixed by different conditions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the 2-D full Euler stability result with C^{1,alpha} interior regularity and general nozzle perturbations that the 3-D result extends."},{"cited_title":"Differential Equations 248 (2010), no","cited_arxiv_id":null,"evidence_quote":"Treats 3-D axisymmetric flows with zero angular momentum; this paper removes that restriction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the regularity characterization of axisymmetric solenoidal vector fields used in Lemma 2.8."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows how to turn a singular axisymmetric elliptic equation into a solvable elliptic system; the same device is central here."},{"cited_title":"105, Princeton University Press , Princeton, NJ, 1983","cited_arxiv_id":null,"evidence_quote":"Provides the freezing-coefficients elliptic regularity theory used in Lemmas 3.9 and 3.10."}],"review_version":1}