{"id":"112c5470-edc4-4eaa-9a15-c9d21239859c","arxiv_id":"1908.02463","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For almost flat nozzles, the paper proves that a scalar solvability condition determines admissible transonic shock locations, with multiple locations possible when the nozzle wall has both expanding and contracting parts.","lead":"In a nearly flat nozzle where gas enters faster than sound and exits at a fixed pressure, this paper derives a condition on the wall shape that decides where the internal shock wave can stand. The same exit pressure can allow several shock positions when the wall is not monotone, so the flow pattern is not unique.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Contractiveness of the iteration mapping in Lemma 5.3 is asserted rather than proved; the existence theorem depends on this missing estimate.","rationale":"The reader's weakest-assumption analysis correctly identifies the unproved contraction estimate in Lemma 5.3 as the load-bearing gap. The paper's main novelty is the construction of the initial approximation via the solvability condition (2.72) and the nonlinear iteration around it, both of which are substantial. The linear well-posedness theory in the appendix and the expansion of the solvability functional I in Lemma 5.1 are plausible and largely carried out. The missing step is the verification that the nonlinear iteration is contractive, which is the point where the fixed-point argument would actually produce a solution. Since this step is explicitly deferred to 'analogous computations', the argument as written does not establish the main theorem. However, there is no internal inconsistency or obvious counterexample; the missing estimates appear likely to hold because most differences scale with a positive power of σ. I therefore do not recommend rejection, but the paper should be conditional on the authors supplying the full proof of Lemma 5.3. The reader's CONDITIONAL verdict is appropriate and unchanged.","tokens_in":38680,"tokens_out":11170,"duration_ms":110490,"concrete_test":"Complete the omitted computation in Lemma 5.3. Starting from (5.21), apply Theorem 4.3 to bound the norm of (δU*2−δU*1, δψ*2′−δψ*1′) by the right-hand side of (4.32)–(4.34) with data (F2−F1, G2−G1, δP3(δU2)−δP3(δU1), δΘ4(ξ;δξ2*)−δΘ4(ξ;δξ1*)). Then use (5.22) and the definitions of f_j, g_j, δP3, δΘ4 to check that each term is bounded by C σ^α times (||δU2−δU1|| + ||δψ'2−δψ'1||) with α>0, uniformly over Kσ. If the optimal α is positive and Cσ^α < 1/2 on the claimed σ-range, the contraction is verified and the gap is fillable. If any essential term yields α=0, the iteration is not contractive and the proof needs a new mechanism.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, Theorem 2.4, is established by a fixed-point argument whose convergence rests entirely on Lemma 5.3, which asserts that the iteration mapping Js is contractive on Kσ. The proof of Lemma 5.3, however, is not carried out. After deriving the bound (5.22) on |δξ2*−δξ1*|, the text states: 'by analogous computations as in Lemma 5.2, with the help of the estimate (5.22), we can show that the inequality (5.20) holds'. Inequality (5.20) is exactly the contraction estimate: the distance between the two iterates is bounded by half the distance between the inputs. Without this inequality, the iteration scheme may fail to converge and the existence of the transonic shock solution with the stated estimates is not proved. The missing estimates are not a formality: (5.21) involves differences of the nonlinear source terms f_j, g_j, δP3, and δΘ4, each a nonlinear function of (δU, δψ', δξ*). The mapping depends on the updated shock location through ψ and on δξ*, which is itself determined implicitly by (4.23). The passage from (5.22) to (5.20) requires Lipschitz estimates in the W^{1−1/β,β} trace norms for compositions with ψ and for products of O(σ) quantities; these are delicate and are not shown. Thus the manuscript leaves a load-bearing step of the proof unverified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies two-dimensional steady compressible Euler flows in a finite nozzle that is a small perturbation of a flat nozzle, with a prescribed receiver pressure at the exit. Because a flat nozzle admits normal shocks at arbitrary locations, the authors propose a linearized free boundary problem whose solvability condition determines an initial approximation ξ* of the shock location via equation (2.72), R(ξ*) = P*. They prove that for monotone nozzles this condition has a unique root, while for nozzles with both expanding and contracting portions multiple roots can occur. The main result, Theorem 2.4, asserts that for any root ξ* satisfying Θ(ξ*)≠0 and sufficiently small wall perturbation amplitude σ, there exists a transonic shock solution to the nonlinear free boundary problem JNSL K, with the shock front close to ξ* and satisfying the estimates (2.75)–(2.80). The proof is based on a nonlinear iteration scheme in which the linearized problem (4.11)–(4.19) is solved first, the shock-position correction δξ* is determined from the solvability condition (4.23), and the iteration mapping Js is shown to be well-defined and contractive on a small ball Kσ(˙U+;˙ψ′).","tokens_in":38989,"tokens_out":4146,"duration_ms":43037,"significance":"If Theorem 2.4 is correct, it provides a rigorous first-principles criterion for admissible transonic shock locations in almost flat nozzles under the physical Courant–Friedrichs exit-pressure condition, and it predicts genuine non-uniqueness of shock solutions for non-monotone nozzle walls. The linearized analysis is a genuine strength: the solvability condition (2.72) is derived from the elliptic system rather than imposed, the coefficients in Lemma 2.1 are computed explicitly, and the elliptic boundary-value theory in Appendix A is developed in a self-contained manner with explicit solvability conditions. The paper also clearly identifies the dependence of the smallness constant on 1/|Θ(ξ*)| and states the range condition for the receiver pressure. However, the central existence proof relies on a contraction estimate in Lemma 5.3 that is asserted rather than proved; until that estimate is supplied, the main theorem must be regarded as conditional.","major_comments":[{"comment":"The proof of the contraction estimate (5.20) is not completed. After establishing (5.22), the text states that 'by analogous computations as in Lemma 5.2, with the help of the estimate (5.22), we can show that the inequality (5.20) holds.' Inequality (5.20) is exactly the contraction bound needed for the fixed-point argument: it requires the difference of two iterates to be bounded by half the distance between the input states. This is the load-bearing step for Theorem 2.4, since without contractiveness the iteration scheme does not produce the existence of a solution to JNSL K. The passage from (5.22) to (5.20) involves estimating differences of the nonlinear source terms f_j, g_j, δP3, and δΘ4 in the norms appearing in Theorem 4.3; these include compositions such as δP3(η;δU) through Y(L,η;δU) and δΘ4(ξ;δξ*) through the map Πψ, as well as products of O(σ) quantities in the trace norms W^{1−1/β,β}. These estimates are not shown, and they are not immediate consequences of (5.22). The authors should write out the 'analogous computations' in full; as it stands, the manuscript does not prove the existence theorem it announces.","section":"Section 5, Lemma 5.3"}],"minor_comments":[{"comment":"In the proof of Lemma 5.2, the estimate for the term involving g4 is delegated by the sentence 'Analogous computations show also that' before (5.18). While this estimate is less central than the contraction estimate, providing the explicit bound for g4 would make the proof more complete and would illustrate the pattern that is later claimed in Lemma 5.3.","section":"Section 5, Lemma 5.2"},{"comment":"There is a typo in the first paragraph: 'constant efficients' should read 'constant coefficients.'","section":"Section 3.1"},{"comment":"The phrase 'For a expanding nozzle' should be 'For an expanding nozzle.'","section":"Remark 2.6"},{"comment":"The abstract states that for strictly expanding or contracting nozzles 'there exists only one solution,' which could be misread as global uniqueness of the nonlinear shock problem. Remark 2.6 clarifies that only one solution is established by the present argument and that global uniqueness has not been proved; the abstract should be rephrased to match this qualification.","section":"Abstract and Remark 2.6"},{"comment":"The function R(ξ) is defined twice, in (2.70) and in (3.21), with identical formulas. The duplication is harmless but could be streamlined by defining it once in Section 2 and referring back to it in Section 3.","section":"Section 2.4 and Section 3.3"}],"recommendation":"major_revision","confidential_remarks":"The main technical gap is the missing proof of the contraction estimate in Lemma 5.3. The framework suggests that the estimate can be supplied by writing out the 'analogous computations' from Lemma 5.2, but as submitted the central existence theorem is not proven. I recommend major revision, with the expectation that the authors provide a complete proof of Lemma 5.3 and, ideally, of the delegated g4 estimate in Lemma 5.2."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper for one reason: it attacks the flat-nozzle degeneracy head-on. Instead of pinning the shock to a wall point as earlier work did, the authors design a linearized free-boundary problem whose solvability condition (2.72) selects the initial shock location from the wall geometry and exit pressure. That is a real idea, and it yields a clean prediction: monotone nozzles give one admissible location, non-monotone nozzles can give several. The explicit R(ξ) = P* equation and the sin^2 example are easy to grasp and honestly presented. The elliptic solvability in Section 3 and Appendix A is careful and appears sound; the derivation of (2.72) from the linearized Rankine-Hugoniot conditions is first-principles and not fitted.\n\nWhere I agree with your reader is the soft spot: Lemma 5.3, the contractiveness of the iteration map Js, is not proved. The text jumps from the estimate (5.22) on |δξ2*−δξ1*| to 'by analogous computations as in Lemma 5.2' for the key inequality (5.20). That is not a mere cosmetic omission. The right-hand side of (5.21) involves differences of nonlinear compositions with ψ, products of O(σ) quantities, and the implicitly defined δξ*, and the trace-norm Lipschitz estimates are delicate. Without (5.20), the existence theorem does not follow from the written argument. That said, I do not see a fundamental obstruction: the structure of the iteration is standard for this type of free-boundary problem, and the authors have done the hard work of establishing the linear theory and the derivative of the solvability condition. The gap looks fixable with a substantial but routine computation, not a new idea.\n\nThe paper deserves a serious referee. It is novel, the main mechanism is transparent, and the missing contraction estimate is exactly the sort of thing a referee can ask to be written out. I would not desk-reject it; I would send it out with a request to complete Lemma 5.3. For a reader working on transonic shock stability or nozzle flows, this is worth citing for the selection mechanism even before the gap is closed, though I would be careful to distinguish the conditional existence theorem from the unconditioned linear analysis.\n\nBring it to reading group if your group tolerates a long paper with a known gap; the first three sections stand on their own. My recommendation: engage with it, and push the authors to supply the missing contraction estimates.","headline":"A genuinely new approach to shock-location selection in almost flat nozzles, with a real but plausibly fixable gap in the contraction proof of Lemma 5.3.","tokens_in":39489,"tokens_out":883,"would_cite":true,"duration_ms":12408,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35A01","35A02","35B20","35B35","35B65","35J56","35L65","35L67"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that in a nearly flat nozzle, an explicit solvability condition selects admissible transonic shock locations, and wavy walls can yield several shocks for one exit pressure.","keywords":["2-D steady Euler system","transonic shocks","nozzle flow","receiver pressure","free boundary problem","shock front location","non-uniqueness","almost flat nozzle"],"falsifier":"Evaluate the missing estimate behind (5.20): for two nearby states in $K_\\sigma$, compute the right-hand side of (5.21) and check whether the output distance is at most half the input distance; a contraction factor larger than $1/2$, or a smallness condition on $\\sigma$ stronger than the one allowed in Lemma 5.1, would invalidate the proof. As an independent check, solve the 2D steady Euler equations numerically for the wall angle $\\Theta(\\xi)=\\sin^2(k\\pi\\xi/L)$ with receiver pressure in the admissible range: the theorem predicts a shock near each of the $2k$ algebraic roots, so failure to find a shock near any one of them would refute the claim.","tokens_in":38478,"feed_emoji":"💨","tokens_out":13982,"duration_ms":135513,"temperature":0.7,"pith_summary":"The paper aims to fix the location of a transonic shock front in a two-dimensional steady compressible Euler flow through a nearly flat finite nozzle when the receiver pressure at the exit is prescribed. Because a flat nozzle allows the normal shock to sit anywhere, the authors replace the flat background by a linearized free boundary problem and show that its solvability condition reduces to the algebraic equation $R(\\xi_*)=\\dot P_*$, in which $R$ records the integrated effect of the wall-angle perturbation and $\\dot P_*$ records the averaged exit-pressure perturbation. When the wall is strictly expanding or contracting, that equation has one root and the paper proves a transonic shock solution exists near it; when the wall has both expanding and contracting portions, several roots can occur, giving several different shock locations for the same exit pressure. This provides a concrete mechanism for non-uniqueness and for the instability of the unperturbed normal shock under generic small wall perturbations.","feed_headline":"Wall shape selects where transonic shocks stand in a nozzle","feed_subtitle":"One algebraic condition fixes shock locations; wavy walls admit several shocks for one exit pressure.","key_machinery":"The load-bearing object is the linearized free boundary problem formed after the Lagrange transformation straightens the streamlines into a rectangle: the subsonic part of the linearized Euler system becomes a first-order elliptic system for the linearized pressure $\\dot p$ and angle $\\dot\\theta$, with boundary data on the two walls, the exit pressure, and the linearized Rankine-Hugoniot conditions. For such a system on a rectangle, solvability forces an integral compatibility condition; here that condition collapses to $R(\\xi_*)=\\dot P_*$, where $R(\\xi)=\\int_0^L \\Theta(\\tau)\\,d\\tau - \\dot K\\int_0^\\xi \\Theta(\\tau)\\,d\\tau$ and $\\dot P_*$ is the averaged exit-pressure perturbation. This equation determines the initial approximating shock position $\\xi_*$. The nonlinear argument then defines an iteration $\\mathcal{J}_s$ on a small ball around the linearized solution and shows that the iteration stays in that ball (Lemma 5.2) and asserts contractiveness in Lemma 5.3, which yields the actual shock solution via a fixed point.","core_discovery":"The central claim is Theorem 2.4: if $\\xi_*$ solves (2.72) and $\\Theta(\\xi_*)\\neq 0$, then for any sufficiently small wall deviation $\\sigma$ there is a genuine transonic shock solution whose front lies within order $\\sigma$ of the vertical line $\\xi=\\xi_*$, with the subsonic state and shock slope close to the linearized approximation. The discovery is that the shock position is selected by the solvability condition of the first-order elliptic system satisfied by the linearized pressure and flow angle behind the shock. In a strictly monotone nozzle this selection is unique, while in a nozzle with alternating expanding and contracting parts, for instance $\\Theta(\\xi)=\\sin^2(k\\pi \\xi/L)$, the equation can have $2k$ roots, and each root with nonzero wall slope generates an admissible shock solution. Thus the same prescribed exit pressure can support multiple transonic shocks, and the flat nozzle's one-parameter family of normal shock locations breaks into finitely many admissible positions determined by the wall shape.","pith_inferences":["Because $R$ depends on $\\Theta$ only through its integrals, two nozzle shapes with the same cumulative wall-angle profile should have identical leading-order shock positions; this is a testable prediction the paper does not state.","The sign of $\\Theta(\\xi_*)$ makes $R$ locally monotone increasing or decreasing, which suggests that the selected shock locations might carry different dynamical stability properties in an unsteady setting, mirroring the quasi-one-dimensional picture; the paper does not address time dependence.","Varying the exit-pressure perturbation $\\dot P_*$ across the range of $R$ should make admissible shock positions appear and disappear in pairs at extrema of $R$, so the algebraic condition supplies a one-dimensional bifurcation diagram for shock location versus receiver pressure."],"forward_implications":["For a strictly expanding or contracting almost flat nozzle with receiver pressure in the admissible range, the theorem yields existence of a transonic shock solution and the leading-order shock location is unique.","For a non-monotone wall, the same receiver pressure can produce several transonic shock solutions, one near each algebraic root of (2.72) at which the wall slope is nonzero.","The flat nozzle's continuum of normal shock positions is not structurally stable: generic small wall perturbations reduce the possible shock locations to a discrete set selected by the integrated wall angle and exit pressure.","The shock front and subsonic flow are, up to order $\\sigma$, exactly what the linearized free boundary problem predicts, with error of order $\\sigma^{3/2}$ between the nonlinear solution and the linearized approximation."],"supporting_citations":[{"why":"provides the physical boundary condition of prescribed receiver pressure at the exit and the flat-nozzle normal shock background.","marker":"[14]"},{"why":"introduces the Lagrange transformation and the elliptic-hyperbolic decomposition used to separate the subsonic problem.","marker":"[8]"},{"why":"supplies the structural-stability framework for transonic shocks in diverging nozzles that the nonlinear iteration adapts.","marker":"[22]"},{"why":"extends that framework to general two-dimensional de Laval nozzles and provides the regularity and iteration estimates.","marker":"[27]"},{"why":"gives the classical quasilinear hyperbolic boundary-value theory used for the supersonic region ahead of the shock.","marker":"[28]"},{"why":"provides the elliptic estimates and trace theorems for nonsmooth domains used to handle corner singularities.","marker":"[21]"},{"why":"establishes regularity limitations at the shock-wall intersection, justifying the chosen function spaces.","marker":"[38]"},{"why":"demonstrates multiple steady transonic shock states in quasi-one-dimensional nozzle flow, the phenomenology the paper extends to two dimensions.","marker":"[17]"}],"fun_headline_variants":["Wall geometry decides transonic shock stand","Algebraic rule pins shock front in nozzles","Nozzle shape controls shock placement","Multiple transonic shocks from one exit pressure","One exit pressure, many shock fronts"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the nonlinear iteration mapping $\\mathcal{J}_s$ is contractive on its small ball (Lemma 5.3); the text does not prove that contractiveness directly, saying only that it follows by analogous computations to Lemma 5.2, so the existence theorem stands or falls on that omitted estimate.","fun_headline_variants_meta":{"raw":{"variants":["Wall geometry decides transonic shock stand","Algebraic rule pins shock front in nozzles","Nozzle shape controls shock placement","Multiple transonic shocks from one exit pressure","One exit pressure, many shock fronts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000761,"raw_usage":{"total_tokens":3408,"prompt_tokens":1006,"completion_tokens":2402,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":622,"completion_tokens_details":{"reasoning_tokens":2338}},"tokens_in":622,"tokens_out":2402,"duration_ms":18203,"temperature":1.0,"reasoning_tokens":2338,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:43:53.705500+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the missing estimate behind (5.20): for two nearby states in $K_\\sigma$, compute the right-hand side of (5.21) and check whether the output distance is at most half the input distance; a contraction factor larger than $1/2$, or a smallness condition on $\\sigma$ stronger than the one allowed in Lemma 5.1, would invalidate the proof. As an independent check, solve the 2D steady Euler equations numerically for the wall angle $\\Theta(\\xi)=\\sin^2(k\\pi\\xi/L)$ with receiver pressure in the admissible range: the theorem predicts a shock near each of the $2k$ algebraic roots, so failure to find a shock near any one of them would refute the claim.","supporting_citations":[{"cited_title":"Courant, K.O","cited_arxiv_id":null,"evidence_quote":"provides the physical boundary condition of prescribed receiver pressure at the exit and the flat-nozzle normal shock background."},{"cited_title":"Chen; Stability of transonic shock fronts in two-dimensional Euler systems.Trans","cited_arxiv_id":null,"evidence_quote":"introduces the Lagrange transformation and the elliptic-hyperbolic decomposition used to separate the subsonic problem."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the structural-stability framework for transonic shocks in diverging nozzles that the nonlinear iteration adapts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"extends that framework to general two-dimensional de Laval nozzles and provides the regularity and iteration estimates."},{"cited_title":"Li, W.-C","cited_arxiv_id":null,"evidence_quote":"gives the classical quasilinear hyperbolic boundary-value theory used for the supersonic region ahead of the shock."},{"cited_title":"Grisvard,Elliptic Problems in Nonsmooth Domains, Monographs and Studies in Mathe- matics, 24, Pitman: Boston, 1985","cited_arxiv_id":null,"evidence_quote":"provides the elliptic estimates and trace theorems for nonsmooth domains used to handle corner singularities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"establishes regularity limitations at the shock-wall intersection, justifying the chosen function spaces."},{"cited_title":"Embid, J","cited_arxiv_id":null,"evidence_quote":"demonstrates multiple steady transonic shock states in quasi-one-dimensional nozzle flow, the phenomenology the paper extends to two dimensions."}],"review_version":1}