{"id":"9b692eec-2b1e-448a-b3f1-5d4dc7e1eb70","arxiv_id":"1908.01694","paper_version":2,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors prove existence, uniqueness, and structural stability of transonic shocks for the 3D axisymmetric Euler system with swirl under small perturbations of the incoming supersonic flow, nozzle boundary, and exit pressure.","lead":"This paper proves that transonic shocks in a 3D axisymmetric nozzle remain stable under small perturbations of the incoming flow, the nozzle wall, and the exit pressure, even when the flow has swirl. The result matters because it extends the mathematical foundation of supersonic-to-subsonic shock stability in realistic engine nozzles to three dimensions with rotation.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Proposition 3 is the sole gate to the C^{2,α} fixed point; its proof defers the coercivity/Fredholm and corner-regularity checks to [19,22] without verifying them, so the main theorem is only as solid as that unshown estimate.","rationale":"I found no explicit mathematical error. The Lagrangian coordinate choice is natural and the reduction to (80) is detailed; the smoothness of the coefficients near the axis after the ζ-transformation is credible. However, the proof of Proposition 3 is a one-paragraph sketch, and Proposition 4 explicitly defers again to [19, Lemma B.3], while Theorem 2 is stated without proof. These are exactly the spots where a hidden compatibility assumption could invalidate the uniqueness and the contraction argument. The reader's weakest-assumption analysis already targets Proposition 3; my stress test agrees and sharpens it to the coercivity/Fredholm and corner-regularity checks. This does not warrant rejection: the prior literature makes the result plausible, and the burden is completion of a deferred proof, not correction of a demonstrated flaw. The verdict should remain CONDITIONAL, hence unchanged relative to the reader's verdict, with the condition being a complete proof of Proposition 3 and, if Theorem 2 is part of the claim, its proof.","tokens_in":26738,"tokens_out":10919,"duration_ms":118860,"concrete_test":"Independently write out the H1(E1) bilinear form for (80), including the trace term a3(z1)Υ*(0,z') and the Robin condition at ζ1=0, and verify coercivity and the Fredholm alternative with explicit signs and bounds on a3 and a4. Then check the hypotheses of [22, Thms. 5.36, 5.45, 4.6] at the axis ρ=0 and at the two corner circles (ζ1=0,ρ=M) and (ζ1=N,ρ=M), stating the compatibility conditions on G1,G2,G3. If coercivity or the corner conditions fail, or require assumptions not present in Theorem 1, the main theorem is unsupported; if the proof goes through unchanged, the conditional concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central fixed-point iteration closes only if the linear elliptic problem (80) has a unique solution satisfying the weighted Hölder bound (81). Proposition 3 is the only place this is established, and its proof is a sketch: it asserts a variational structure 'similar' to [19, Lemma 4.3], invokes Lax-Milgram/Fredholm, and then cites [22, Thms. 5.36, 5.45, 4.6] for global regularity and the final estimate. What is not checked is (i) coercivity of the bilinear form with the nonlocal term a3(z1)Υ*(0,z2,z3) and the Robin coefficient a4 in the ζ1=0 boundary condition; (ii) the compactness/Fredholm step after the trace term is moved to the right side; (iii) the corner compatibility conditions at ζ1=0, ζ1=N, and ρ=M needed for C^{2,α} up to the boundary, since the weighted norm only weights Γw,ζ; and (iv) the hypotheses of Lieberman's theorems in this singular cylindrical geometry. Proposition 4 then bootstraps W2,W4 to C^{2,α} using a second 'similar' argument and [19, Lemma B.3], and Theorem 2's proof is omitted. Every contraction estimate in Section 3.4 passes through Proposition 4, so a failure or a need for stronger compatibility in any of these steps would prevent the fixed point from closing. The paper gives no evidence of such a failure; the concern is that the central estimate is not actually demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the steady three-dimensional axisymmetric compressible Euler system with nonzero swirl in a divergent nozzle and claims a structural stability result for the spherical symmetric transonic shock. Theorem 1 asserts that, under small perturbations of the incoming supersonic flow, the nozzle wall, and the exit pressure, there is a unique transonic shock solution whose shock front and subsonic flow satisfy weighted Hölder estimates of size O(ε). The proof introduces an invertible Lagrangian transformation that straightens streamlines, fixes the unknown shock as a free boundary, decomposes the system into an elliptic system for the flow angle and pressure plus transport equations, and then runs a contraction mapping in weighted Hölder spaces. Theorem 2 states a C^{2,α} regularity improvement in the case of a straight nozzle under additional compatibility conditions. The central structure of the iteration is coherent and detailed, but the key elliptic estimate (Proposition 3) is only sketched, and the proof of Theorem 2 is omitted.","tokens_in":27058,"tokens_out":3782,"duration_ms":35538,"significance":"If the main theorem is correct, it would be a meaningful extension of the transonic shock stability literature: the axisymmetric setting with swirl has both the artificial singularity at the axis and corner singularities at the shock-wall intersection, and the proposed Lagrangian transformation appears to be a genuine new ingredient that overcomes the earlier non-invertibility near the axis. The fixed-point scheme is standard in spirit but is adapted with care to the singular coefficients, and the paper gives credit to the prior framework of Li-Xin-Yin. The main weakness is that the load-bearing elliptic regularity result is not actually proved in the manuscript: Proposition 3 defers the existence, uniqueness, and weighted Hölder estimate to earlier papers without checking their hypotheses, and Proposition 4 and Theorem 2 contain further deferred arguments. Thus the paper is not yet self-contained for its central claims, though no direct error in the fixed-point estimates is apparent.","major_comments":[{"comment":"The proof of Proposition 3, which is the only route to the weighted C^{2,α} estimate needed for the fixed-point iteration, is a sketch. It asserts a variational structure 'similar' to [19, Lemma 4.3] and obtains existence and uniqueness by Lax-Milgram and Fredholm, but it does not identify the bilinear form, verify its coercivity in the presence of the nonlocal term a3(ζ1)Υ*(0,ζ2,ζ3) and the Robin coefficient a4 at ζ1=0, or justify the compactness/Fredholm step after moving the trace term. It then invokes [22, Theorems 5.36 and 5.45] and [22, Theorem 4.6] for global regularity and the final estimate without checking the hypotheses of those theorems in the present cylindrical, mixed-boundary, corner domain. Because estimate (81) feeds directly into Proposition 4 and hence into every contraction estimate in Section 3.4, the main theorem is not fully demonstrated until this proposition is proved in detail.","section":"Section 3.2, Proposition 3 (Eqs. (80)-(81))"},{"comment":"The bootstrap from C^{1,α} regularity to the stated C^{2,α} regularity for W2 and W4 is not self-contained. The estimate (86) for W4 is asserted by saying that the argument is 'similar to the proof of Proposition 3,' and the final regularity for W2 is attributed to an argument 'similar to [19, Lemma B.3].' The passage from the first-order normal form (85) to the weighted C^{2,α} bound (82) requires an explicit treatment of compatibility conditions at the corners (0,0), (N,0), (0,M), and (N,M), since the norm (20) weights only a portion of the boundary. Without this, the iteration map is not shown to take the space Ξδ into itself.","section":"Section 3.2, Proposition 4 (Eqs. (82)-(86))"},{"comment":"Theorem 2 is stated as a main result, but its proof is omitted with the comment that it is 'very similar to the proof for [19, Theorem 1.1].' The compatibility lemmas (Lemmas 6 and 7) provide useful information, but the final assertions (24)-(25) require global C^{2,α} regularity up to the intersection point of the shock with the nozzle wall and up to the axis, and the reference to a 'standard even extension' is not enough to establish that the transformed equations and boundary conditions preserve that regularity. Since Theorem 2 is advertised as a new result, it should either be proved in the paper or explicitly stated as a corollary whose full proof appears elsewhere.","section":"Section 4, Proof of Theorem 2"}],"minor_comments":[{"comment":"The word 'probelm' in the statement of Proposition 4 is a typo for 'problem.'","section":"Proposition 4 statement"},{"comment":"Reference [11] is listed as 'Supsonic flow and Shock Waves'; the correct spelling is 'Supersonic.'","section":"References"},{"comment":"The weighted Hölder spaces are defined in the introduction, but the notation Γw,s is introduced only in (20), while the shock front estimate (18) uses the weight {θ*}; it would be clearer to define both sets precisely in the same paragraph.","section":"Section 1, Theorem 1"},{"comment":"The remainder terms Ri in (42)-(44) depend on Φ− evaluated at the shock position, but this dependence is encoded only in the notation; writing the arguments explicitly, as done for the listing of Φ± in (43), would improve readability.","section":"Section 2.2, Eq. (44)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious contribution to a well-established line of work on transonic shocks, and the proposed Lagrangian transformation is a credible new idea. The main obstacle to acceptance is that the central elliptic estimate (Proposition 3) and the higher regularity result (Theorem 2) are deferred rather than proved; these are load-bearing, not cosmetic, so the revision should provide the missing proofs or a precise reduction to existing theorems with hypotheses verified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a real paper with a real new idea, and if you are in the transonic shock business you should read it. The authors prove structural stability of the spherical transonic shock in a 3D axisymmetric divergent nozzle under small perturbations of the incoming flow, the nozzle wall, and the exit pressure, for Euler with swirl. That extends the prior two-dimensional and swirl-free axisymmetric results [19,20]. The genuinely new ingredient is an invertible Lagrangian transformation that avoids the sinθ degeneracy at the axis, so the streamlines can be straightened even with swirl. That is not cosmetic; the earlier transformation simply failed there. The reformulation is careful, and the paper is written in the standard style of this area, reducing the free boundary problem to a fixed point for a first-order elliptic system plus transport equations.\n\nWhat it does well: the nonlinear estimates are presented in detail. The quadratic error terms, the singular terms near the axis, and the boundary conditions are tracked. The fixed-point iteration in Theorem 1 is coherent; if the linear elliptic estimate holds, the contraction argument closes. The paper does not hide behind fitted parameters or numerics; it is a proof-driven paper.\n\nWhere the soft spots are: Proposition 3 is the load-bearing gate. It asserts existence, uniqueness, and the weighted C^{2,α} bound for the linearized problem (80), and its proof is a sketch: a variational structure 'similar' to [19, Lemma 4.3], an appeal to Lax-Milgram and Fredholm, then Lieberman's book for global regularity. What is not checked is coercivity with the nonlocal term a3(z1)Υ*(0,z2,z3) and the Robin coefficient, the compactness step after moving the trace term, and the corner compatibility conditions needed for C^{2,α} up to the boundary in this cylindrical domain. Proposition 4 then uses a second 'similar' argument plus [19, Lemma B.3] to bootstrap to W2,W4, and every contraction estimate in Section 3.4 passes through Proposition 4. So the main theorem is only as solid as that unshown estimate. I want to be fair: nothing in the paper suggests the estimate is wrong, and the regularity theory of Lieberman likely does apply, but the proof is not actually demonstrated here. Theorem 2's proof is omitted entirely; that is a secondary result, but it is still an omission.\n\nWho should read this: PDE analysts working on free boundary problems for steady Euler, transonic shocks, or corner/axis singularities. This deserves a serious referee. The referee should be asked specifically to verify Proposition 3, and the authors should be asked to expand that proof or supply an appendix. It is a significant result, conditional on that estimate being true.","headline":"A genuine new idea—an axis-regular Lagrangian transformation—powers a credible transonic shock stability result for 3D axisymmetric Euler with swirl, but the main theorem leans on a key elliptic estimate whose proof is sketched rather than shown.","tokens_in":27572,"tokens_out":3073,"would_cite":true,"duration_ms":28914,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35L65","35L67","76N15"],"pacs":[],"model":"deepseek-v4-flash","headline":"Small axisymmetric perturbations of the incoming flow, the nozzle wall, and the exit pressure do not destroy the spherical transonic shock: a unique nearby solution with swirl exists.","keywords":["transonic shock","structural stability","axisymmetric Euler system","swirl","Lagrangian transformation","divergent nozzle","free boundary problem","weighted Hölder spaces"],"falsifier":"Solve the homogeneous version of the linear elliptic problem (80) with the background coefficients in the weighted Hölder spaces used in the paper: a nonzero solution, or a sequence of solutions whose weighted norms blow up as the axis is approached, would contradict Proposition 3 and the fixed-point argument for Theorem 1 would not close.","tokens_in":26544,"feed_emoji":"🌀","tokens_out":13963,"duration_ms":124078,"temperature":0.7,"pith_summary":"The paper establishes that the spherical symmetric transonic shock in a divergent axisymmetric nozzle survives small axisymmetric perturbations of the incoming supersonic flow, the nozzle wall, and the exit pressure, even when the flow carries swirl. Previous stability results for multidimensional transonic shocks typically excluded swirl or worked in geometries where the corner singularity at the shock-wall intersection and the artificial singularity at the symmetry axis do not appear together. The authors construct a new invertible Lagrangian coordinate system that straightens the streamlines through the whole nozzle, including near the axis, and reformulate the free-boundary problem as an elliptic system for the flow angle and pressure coupled to an ordinary differential equation for the shock front. A fixed-point contraction then yields a unique perturbed solution whose shock location and subsonic flow lie within a constant multiple of the perturbation size of the spherical background. This gives a rigorous foundation for the standard picture of divergent-nozzle operation: small changes in upstream conditions, wall shape, or back pressure move the shock slightly without destroying it.","feed_headline":"Transonic shocks hold with swirl in 3D nozzles","feed_subtitle":"Small changes in inflow, wall shape, or back pressure leave a unique shock front close to the spherical baseline.","key_machinery":"The central object is a new invertible Lagrangian coordinate map $(r,\\theta)\\mapsto(y_1,y_2)$ with $y_1=r$ and $y_2=\\tilde y_2^{1/2}$, where $\\tilde y_2$ is the axisymmetric stream function determined by $d\\tilde y_2=r^2\\rho U_1\\sin\\theta\\,d\\theta-r\\rho U_2\\sin\\theta\\,dr$. Because $\\sin\\theta=O(\\theta)$ near the axis, the Jacobian $r^2\\rho U_1\\sin\\theta/(2y_2)$ stays bounded away from zero, so the map is invertible even at the axis, unlike the usual Lagrangian coordinates used for two-dimensional nozzles. This map straightens the streamlines into rectangles, lets the entropy, swirl, and Bernoulli equations be integrated explicitly, and turns the free-boundary problem into a first-order elliptic system for the flow angle $\\varpi=U_2/U_1$ and pressure $P$, coupled to an ODE for the shock front. The elliptic system is rewritten through a potential function and solved in a cylindrical domain in weighted Hölder spaces, with the nonlocal terms controlled by the fixed-point contraction.","core_discovery":"The paper proves structural stability of the spherical symmetric transonic shock for the three-dimensional axisymmetric Euler system with swirl: for every sufficiently small axisymmetric perturbation of the incoming supersonic flow, the nozzle wall, and the exit pressure, there exists a unique solution consisting of a supersonic part, a transonic shock surface $r=\\xi(\\theta)$, and a subsonic part, with the shock front and the subsonic flow satisfying $\\|\\xi-r_b\\|\\le C_0\\epsilon$ and $\\|\\Psi^+-\\Psi_b^+\\|\\le C_0\\epsilon$ in the appropriate weighted Hölder norms. The physical entropy condition across the shock is preserved, and the shock location varies continuously with the data. When the wall is straight and the data satisfy extra compatibility conditions, the solution gains higher regularity up to the wall and the axis.","pith_inferences":["The same square-root Lagrangian coordinate should apply to other axisymmetric free-boundary problems in gas dynamics wherever streamlines need straightening without losing the axis, since the only structural input is $\\sin\\theta=O(\\theta)$ near the axis.","Because the contraction constant in the proof is proportional to $\\epsilon$, the argument yields an explicit, if conservative, radius of stability; a numerical continuation in $\\epsilon$ could test how far the linear bound $C_0\\epsilon$ remains valid beyond the regime covered by the proof.","Theorem 2's compatibility conditions suggest that, for straight walls, regularity is limited mainly by corner compatibility; relaxing them would require resolving the corner singularities directly, which is the part of the elliptic estimate the paper sketches rather than fully details."],"forward_implications":["For every sufficiently small axisymmetric perturbation of the incoming supersonic flow, the nozzle wall, and the exit pressure, a unique transonic shock solution exists, with the shock front $\\xi(\\theta)$ satisfying $\\|\\xi-r_b\\|_{3,\\alpha}^{(-1-\\alpha;\\{\\theta^*\\})}\\le C_0\\epsilon$.","The subsonic flow behind the shock stays within $C_0\\epsilon$ of the spherical background in the weighted Hölder norm, and the entropy condition $P^+>P^-$ across the shock is preserved.","The shock front and the subsonic flow depend continuously on the three perturbation data, with constants depending only on the background solution and the boundary data, not on the particular perturbation.","If the nozzle wall is unperturbed and the data satisfy additional compatibility conditions, the shock front gains $C^{3,\\alpha}$ regularity and the subsonic flow gains $C^{2,\\alpha}$ regularity up to the wall and the axis.","Because the transport equations for entropy, swirl, and the Bernoulli function are solved explicitly and the shock front is governed by an ODE, the free-boundary problem reduces to a fixed-boundary elliptic problem whose uniqueness follows from the contraction map."],"supporting_citations":[{"why":"Supplies the overall strategy: straighten streamlines with Lagrangian coordinates and reduce the shock problem to a second-order elliptic equation with a nonlocal term plus an ODE for the shock front.","marker":"[20]"},{"why":"Provides the elliptic boundary-value framework for transonic shocks in conic nozzles and the lemma invoked for existence and uniqueness of the analogous elliptic problem.","marker":"[19]"},{"why":"Gives the weighted Hölder regularity theorems used to control the solution of the elliptic problem near the axis and the wall.","marker":"[22]"},{"why":"Establishes the optimal boundary regularity for subsonic nozzle flows and the even-extension lemma used in the higher-regularity part.","marker":"[26]"},{"why":"Supplies the characteristic-method existence theory for the perturbed supersonic part from the inlet data.","marker":"[15]"},{"why":"Describes the basic transonic shock problem in a divergent nozzle and the background shock solution whose stability is at issue.","marker":"[11]"},{"why":"Constructs the spherical symmetric transonic shock solution used as the unperturbed background for the stability estimates.","marker":"[29]"}],"fun_headline_variants":["Transonic shocks stay stable with swirl in 3D nozzles","Swirling transonic shocks persist under small perturbations","3D axisymmetric nozzle shocks stable to flow and wall changes","Shock front stability proven for swirling nozzle flows","Unique transonic shock survives swirl and nozzle twists"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the linear elliptic boundary-value problem obtained after fixing the shock front has a unique solution with the stated weighted Hölder bounds; the proof invokes standard elliptic theory for this, so the singular corner and axis behavior must be covered by those estimates for the fixed-point iteration to close.","fun_headline_variants_meta":{"raw":{"variants":["Transonic shocks stay stable with swirl in 3D nozzles","Swirling transonic shocks persist under small perturbations","3D axisymmetric nozzle shocks stable to flow and wall changes","Shock front stability proven for swirling nozzle flows","Unique transonic shock survives swirl and nozzle twists"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000654,"raw_usage":{"total_tokens":2938,"prompt_tokens":829,"completion_tokens":2109,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":445,"completion_tokens_details":{"reasoning_tokens":2031}},"tokens_in":445,"tokens_out":2109,"duration_ms":15114,"temperature":1.0,"reasoning_tokens":2031,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:05:26.178359+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the homogeneous version of the linear elliptic problem (80) with the background coefficients in the weighted Hölder spaces used in the paper: a nonzero solution, or a sequence of solutions whose weighted norms blow up as the axis is approached, would contradict Proposition 3 and the fixed-point argument for Theorem 1 would not close.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the overall strategy: straighten streamlines with Lagrangian coordinates and reduce the shock problem to a second-order elliptic equation with a nonlocal term plus an ODE for the shock front."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the elliptic boundary-value framework for transonic shocks in conic nozzles and the lemma invoked for existence and uniqueness of the analogous elliptic problem."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the weighted Hölder regularity theorems used to control the solution of the elliptic problem near the axis and the wall."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the optimal boundary regularity for subsonic nozzle flows and the even-extension lemma used in the higher-regularity part."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the characteristic-method existence theory for the perturbed supersonic part from the inlet data."},{"cited_title":"Courant and K","cited_arxiv_id":null,"evidence_quote":"Describes the basic transonic shock problem in a divergent nozzle and the background shock solution whose stability is at issue."},{"cited_title":"Xin and H","cited_arxiv_id":null,"evidence_quote":"Constructs the spherical symmetric transonic shock solution used as the unperturbed background for the stability estimates."}],"review_version":1}