{"id":"c505935c-7703-4a0d-ae1e-d686a9ae7bec","arxiv_id":"2411.15575","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The two-parameter relaxation model (1.3) converges to the incompressible Navier-Stokes equations; velocity converges for all scalings and pressure converges when delta = o(sqrt(epsilon)), with explicit rates.","lead":"The paper introduces a two-parameter hyperbolic approximation to the incompressible Navier-Stokes equations, combining relaxation and artificial compressibility, and proves convergence to the Navier-Stokes solution as both parameters vanish. The new part is a rigorous error estimate for the pressure variable, achieved by comparing with an auxiliary linear system.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified; central proof closes under the stated δ=O(√ε) well-prepared assumptions. (Minor fixable sign typo in (4.11a) and δ/√ε typo in Cor. 2.5.)","rationale":"The manuscript's central claim is the two-parameter convergence of (1.3) to incompressible Navier-Stokes, with pressure convergence in the regime δ=o(√ε) and well-prepared data. I followed the proof structure: Section 3 gives velocity convergence via symmetric hyperbolic energy estimates and a blow-up criterion; Section 4 introduces the intermediate linear system (4.2), proves the pressure estimate (4.7) through the f,g formulation, proves the velocity estimate (4.4) via divergence/curl splitting, and then applies Proposition 4.1 to transfer the errors to the original system. I checked the main algebraic steps: the derivation of Proposition 4.2, the ODE (4.13), the energy estimate for (4.11), the initial-energy bound (4.16), the vorticity system (4.22), and the H¹ closure via Lemma 2.9. The only definite issue I found is a sign error in the displayed equation (4.11a): with f=δ∂tp'+p'−p^NS and g=√(εδ)∂t∇p', Proposition 4.2 yields LHS = −ε∂t²p^NS, not +ε∂t²p^NS. This is not load-bearing because the subsequent estimate bounds the source term by its absolute value and the crucial cancellation of the √(ε/δ)∇·g terms is unaffected. Corollary 2.5 also contains a typo ('δ√ε→0' should be 'δ/√ε→0'), but the preceding statement δ=o(√ε) is unambiguous. The restrictive δ=O(√ε) assumption is genuinely needed: in (4.16), the term (δ²/ε)||∇div u0||² only fits into C(ε+δ)² when δ≤C√ε, and the application of Theorem 2.2 to obtain the uniform L∞ bound (4.1) requires the pressure initial error to stay bounded, again using δ≤C√ε. These are stated conditions, so they do not undermine the theorem. No circularity appears: (4.4) is proved in Section 4.2 using (4.15) and (4.17), while Proposition 4.1 uses (4.4) only after it has been established. I therefore concur with the reader's ACCEPT verdict and would not adjust it.","tokens_in":20428,"tokens_out":57279,"duration_ms":464444,"concrete_test":"Independently re-derive (4.11a) from Proposition 4.2 with f,g as in (4.10), then re-run the dE/dt computation with the corrected sign and confirm the same bound (4.15) follows. As a second check, trace the (δ²/ε)||∇div u0||² term in (4.16) and verify that δ≤C√ε is exactly what makes it O((ε+δ)²); if either check fails, the pressure estimate would need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I find no load-bearing flaw in the central claim. The velocity convergence follows from the energy estimates of Section 3, and the pressure convergence rests on the auxiliary linear system (4.2), with Proposition 4.1 plus (4.4) and (4.7) closing the triangle without circularity. The sensitive point is exactly what the reader identified: the δ=O(√ε) scaling is used in (4.16) to control (δ²/ε)||∇div u0||², and the well-prepared initial-data conditions are stated hypotheses rather than hidden assumptions. Two minor typos do not affect the estimates: independently re-deriving (4.11a) from Proposition 4.2 and the definitions (4.10) gives (ε∂t²−Δ)f − √(ε/δ)∇·g = −ε∂t²p^NS, not +ε∂t²p^NS, but the sign is irrelevant to the subsequent absolute-value bound and energy estimate; and Corollary 2.5's parenthetical should read δ/√ε→0, not δ√ε→0. Neither issue changes the theorem or its proof.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a two-parameter hyperbolic relaxation approximation (1.3) of the two-dimensional incompressible Navier-Stokes equations, combining a first-order relaxation term with the artificial compressibility method. The main results are: (i) Theorem 2.1 and Theorem 2.2, which prove that smooth solutions exist on a time interval that grows to infinity as (ε, δ) → 0 and that the velocity component converges to the Navier-Stokes velocity in L2, with rates O(ε + δ), under suitable initial-data assumptions; and (ii) Theorem 2.4 and Corollary 2.5, which, under the scaling δ = O(√ε) and well-prepared initial data, prove convergence of the velocity in H1 and of the pressure in H1 with rate √ε + δ/√ε. The key technical device is an auxiliary linear system (4.2). The paper estimates the difference between the original system and the auxiliary system, and separately the difference between the auxiliary system and Navier-Stokes, using energy estimates, an auxiliary pressure equation (4.11), and a vorticity initial-layer correction.","tokens_in":20601,"tokens_out":22209,"duration_ms":178946,"significance":"If the main theorem is correct, this is a valuable contribution: it provides a rigorous two-parameter singular limit for a hyperbolic approximation of the incompressible Navier-Stokes equations, and it gives pressure convergence, which the earlier one-parameter relaxation results did not provide. The proof is self-contained and follows standard energy methods; all assumptions are stated explicitly, and no parameter is fitted to make the convergence work. The auxiliary- system argument is coherent and the term-by-term estimates are plausible. The main limitation is the regime δ = O(√ε) and the well-prepared initial-data hypotheses; the paper explicitly acknowledges that more general parameter relations are left to future work. The manuscript does not provide machine-checked proofs, but the estimates are detailed enough for a human referee to trace. I judge the central claim sound.","major_comments":[],"minor_comments":[{"comment":"The displayed convergence rate and the parenthetical equivalence are incorrect as printed: δ = o(√ε) means δ/√ε → 0, not δ√ε → 0, and the proof of Theorem 2.4 gives ||p^{ε,δ} − p^{NS}||_{H1} ≤ C_T(√ε + δ/√ε), not C_T(√ε + δ√ε). The final conclusion that the pressure converges is still true, but the displayed rate should be corrected.","section":"§2, Corollary 2.5"},{"comment":"There is a sign error in (4.11a): independently re-deriving it from Proposition 4.2 and the definitions (4.10) gives (ε∂_t^2 − Δ)f − √(ε/δ) ∇·g = −ε∂_t^2 p^{NS}, not +ε∂_t^2 p^{NS}. The sign is irrelevant for the subsequent absolute-value and energy estimates, but the displayed equation should be corrected.","section":"§4, Eq. (4.11a)"},{"comment":"In the energy estimate for the vorticity system, the expression 'δtΩ^{NS}' should read '∂_tΩ^{NS}'.","section":"§4.2, text after Eq. (4.22)"},{"comment":"Remark 4.4 contains a LaTeX artifact: 'δ /greaterorsimilar√ε' should read 'δ \\gtrsim √ε' or 'with δ not smaller than √ε'.","section":"§4, Remark 4.4"}],"recommendation":"minor_revision","confidential_remarks":"The paper is within the journal's scope and the central proof closes under the stated assumptions. The only issues I found are local: a sign typo in (4.11a) and a wrong displayed rate in Corollary 2.5. Both should be fixed before publication, but neither undermines the main theorem. I recommend minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a strong analytical paper. It introduces the two-parameter relaxation–artificial compressibility system (1.3), which is new, and proves the first joint convergence result for both parameters, including the first pressure error estimate. The main theorems (2.4 and Corollary 2.5) give H1 convergence of velocity and pressure to the incompressible Navier–Stokes solution under δ = O(√ε) and well-prepared initial data. The proof strategy is sensible: energy estimates for the velocity, then an auxiliary linear system to control the pressure. I went through the key steps, particularly Section 4, and the argument closes. The assumptions are stated explicitly, and the paper does not overclaim.\n\nThe genuinely new content is the pressure estimate. Prior work on the Brenier–Natalini–Puel relaxation did not characterize the pressure, and the classical artificial compressibility estimates for pressure did not have rates. The auxiliary system (4.2) and the decomposition into divergence and curl in Section 4.2 are the right tools, and the proof is careful.\n\nSoft spots are real but not damning. First, the pressure result only works for δ = O(√ε), and the convergent rate for pressure is √ε after fixing δ = o(√ε). Other parameter regimes are not analyzed—the paper says so. Second, the well-prepared initial data conditions are fairly strong, e.g. ||∇div u0|| ≤ C(ε+δ) and δ(||Δp0||+||∇div U0||) ≤ C(ε+δ). They are stated clearly, but they limit the practical applicability. Third, the lower bound on the existence time is logarithmic in the parameters, which means for a fixed time horizon the parameters have to be exponentially small; that is a practical caveat for any numerical use. Fourth, there is no numerical verification, and the proofs are not machine-checked; however, the estimates are detailed enough that I found no hidden assumption.\n\nThere are two typos that should be fixed: in (4.11a) the sign in front of ε∂t²p^NS appears to be wrong, and in Corollary 2.5 the parenthetical should read δ/√ε → 0, not δ√ε → 0. Neither affects the results.\n\nWho is this for? Researchers working on hyperbolic relaxation approximations to incompressible flows, singular limits, and numerical methods based on such systems. It deserves a serious referee; I expect it to be accepted after the minor typos are corrected. I would engage with it.","headline":"Solid, new two-parameter convergence result with the first pressure error estimate; the central proof closes under the stated well-prepared assumptions, with only minor typos to fix.","tokens_in":21177,"tokens_out":3488,"would_cite":true,"duration_ms":29719,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q30","76D05"],"pacs":[],"model":"deepseek-v4-flash","headline":"A new two-parameter hyperbolic approximation is proved to converge to the incompressible Navier-Stokes equations, with explicit pressure error bounds when the relaxation parameter is much smaller than the square root of the…","keywords":["incompressible Navier-Stokes equations","hyperbolic relaxation approximation","artificial compressibility","two-parameter singular limit","pressure convergence","energy estimates","initial-layer correction"],"falsifier":"Solve the approximate system numerically on the two-dimensional periodic square for a smooth test flow with $\\delta=\\sqrt{\\varepsilon}$ (so $\\delta$ is not $o(\\sqrt{\\varepsilon})$) and well-prepared initial data, and measure $\\sup_{0<t<T}\\|p^{\\varepsilon,\\delta}-p^{NS}\\|_{H^1}$ as $\\varepsilon\\to0$. The theorem only guarantees $\\sqrt{\\varepsilon}\\|p^{\\varepsilon,\\delta}-p^{NS}\\|_{H^1}=O(\\varepsilon+\\delta)$, so the unweighted error could diverge like $1/\\sqrt{\\varepsilon}$; observing it instead vanish would show the $\\delta=o(\\sqrt{\\varepsilon})$ threshold is not necessary. Conversely, keeping $\\delta=o(\\sqrt{\\varepsilon})$ but taking data with $\\|\\nabla\\,\\mathrm{div}\\,u_0^{\\varepsilon,\\delta}\\|_{L^2}$ of order one rather than $O(\\varepsilon+\\delta)$ would test whether the well-prepared condition is needed for pressure convergence.","tokens_in":20185,"feed_emoji":"🌊","tokens_out":8421,"duration_ms":66003,"temperature":0.7,"pith_summary":"The paper introduces a two-parameter hyperbolic approximation to the two-dimensional incompressible Navier-Stokes equations, obtained by combining the artificial compressibility method (parameter $\\varepsilon$) with a first-order relaxation of the velocity gradient (parameter $\\delta$). The central claim is that, as $\\varepsilon$ and $\\delta$ both go to zero, the smooth solutions of this new system converge to the Navier-Stokes solution, with explicit error bounds for both velocity and pressure, provided $\\delta = O(\\sqrt{\\varepsilon})$ and the initial data are well prepared. Velocity convergence is obtained by direct energy estimates for a symmetrized residual system; pressure convergence, which the one-parameter methods in the literature did not deliver, is obtained by inserting a linear auxiliary system and estimating the two error legs separately. A sympathetic reader should care because the approximating system is genuinely hyperbolic, hence has finite signal speeds and access to hyperbolic numerical methods, while still provably reproducing incompressible flow.","feed_headline":"Hyperbolic two-parameter model provably converges to Navier-Stokes","feed_subtitle":"Rigorous error bounds cover velocity and pressure on any finite time window when the relaxation parameter shrinks faster than √ε.","key_machinery":"The load-bearing object is the linear auxiliary system (4.2): it is obtained from the two-parameter model by replacing the quadratic term $u\\otimes u$ with the known $u^{NS}\\otimes u^{NS}$, keeping the same initial data. The argument first bounds the difference between the auxiliary system and the incompressible Navier-Stokes solution, then bounds the difference between the original approximation and the auxiliary system, and combines the two. The auxiliary pressure is further controlled through the derived scalar PDE $\\varepsilon\\delta\\,\\partial_t^3 p' - (\\varepsilon+\\delta)\\partial_t\\Delta p' + \\varepsilon\\,\\partial_t^2 p' - \\Delta(p'-p^{NS})=0$, reformulated as a first-order system in $f=\\delta\\partial_t p'+p'-p^{NS}$ and $g=\\sqrt{\\varepsilon\\delta}\\,\\partial_t\\nabla p'$ with an explicit energy; the velocity difference is split into divergence and curl parts, with the vorticity leg requiring an initial-layer correction for the auxiliary vorticity variables.","core_discovery":"On the unit periodic square, for smooth Navier-Stokes solutions and smooth solutions of the new system, the paper proves (Theorem 2.4) that if $\\delta \\leq C\\sqrt{\\varepsilon}$ and the initial data satisfy the stated $H^1$ and divergence-preparation conditions, then for any fixed time $T$ before the approximation's lifespan, $\\sup_{0<t<T}\\big(\\|u^{\\varepsilon,\\delta}(t)-u^{NS}(t)\\|_{H^1}+\\sqrt{\\varepsilon}\\,\\|p^{\\varepsilon,\\delta}(t)-p^{NS}(t)\\|_{H^1}\\big) \\leq C_T(\\varepsilon+\\delta)$. Corollary 2.5 then yields unweighted pressure convergence in $H^1$ with rate $O(\\sqrt{\\varepsilon}+\\delta/\\sqrt{\\varepsilon})$ whenever $\\delta=o(\\sqrt{\\varepsilon})$. Earlier theorems (2.1 and 2.2) establish $L^2$ velocity convergence at rate $O(\\varepsilon+\\delta)$ and show the lifespan tends to infinity as $\\varepsilon+\\delta\\to0$. The proof of pressure convergence relies on a linear intermediate system (4.2) and energy estimates for its differences to both the original approximation and Navier-Stokes; a key vorticity estimate uses an initial-layer correction for the auxiliary vorticity variables.","pith_inferences":["The restriction $\\delta=O(\\sqrt{\\varepsilon})$ is likely essential to the pressure argument, because the $\\sqrt{\\varepsilon}$ in front of the pressure error forces the auxiliary system to track $p^{NS}$ at a faster rate; testing $\\delta\\approx\\sqrt{\\varepsilon}$ would reveal whether the unweighted pressure error genuinely fails to vanish in that regime, or whether another intermediate system could","The techniques rely on the 2D periodic torus through Ladyzhenskaya-type interpolation, Poincar\\'e inequalities, Helmholtz decomposition, and constant mean-zero properties; extending the result to 3D or to bounded domains with physical boundary conditions will require new boundary-layer and interpolation arguments.","The explicit rates suggest a testable numerical prediction: for well-prepared initial data and $\\delta\\ll\\sqrt{\\varepsilon}$, the $H^1$ pressure error over a fixed time interval should scale like $\\sqrt{\\varepsilon}$; a modest numerical experiment on a periodic shear flow could check whether the predicted rate appears.","For non-well-prepared initial data the paper proves nothing about pressure; a natural extension would be to introduce initial-layer corrections in the original system itself, not only in the auxiliary vorticity equations, to remove the preparation conditions."],"forward_implications":["When $\\delta=o(\\sqrt{\\varepsilon})$, the approximate pressure $p^{\\varepsilon,\\delta}$ converges to $p^{NS}$ in $H^1$ at rate $O(\\sqrt{\\varepsilon}+\\delta/\\sqrt{\\varepsilon})$, so the model can be used as a provably consistent pressure proxy, not just a velocity proxy.","For any fixed finite time horizon, smooth solutions of the approximate system exist for all sufficiently small parameters, so the asymptotic statement is uniform in time up to $T$.","Since the approximation is hyperbolic with finite propagation speed, it can be discretized with established hyperbolic balance-law schemes, and the proven rates give a concrete parameter guideline: pick $\\delta$ much smaller than $\\sqrt{\\varepsilon}$ to keep pressure errors small.","Adding a linear friction term to the Navier-Stokes target does not destroy any of the convergence results (Proposition 2.6).","The proof gives a template for two-parameter singular limits in hyperbolic relaxation systems: use an intermediate linear system to transfer convergence to a variable, here pressure, that direct energy estimates cannot control."],"supporting_citations":[{"why":"Introduces the first-order relaxation system (1.2) whose single-parameter limit is the incompressible Navier-Stokes equations, supplying the baseline velocity convergence rate that the new two-parameter result extends.","marker":"[4]"},{"why":"Gives the classical artificial compressibility convergence result for velocity and pressure gradient, motivating the AC component of the new system but lacking the error estimates needed for pressure convergence.","marker":"[30]"},{"why":"Provides the well-established local existence theory for quasilinear symmetric hyperbolic systems, used to guarantee existence and uniqueness for (1.3) and to justify the blow-up criterion in Lemma 2.8.","marker":"[20]"},{"why":"Supplies the Ladyzhenskaya inequality and the two-dimensional interpolation framework used throughout the energy estimates.","marker":"[22]"},{"why":"Presents the basic aspects of hyperbolic relaxation systems and a convergence-stability principle for singularly parametrized hyperbolic systems, organizing the two-parameter limit analysis.","marker":"[33]"},{"why":"Advances the analysis of the one-parameter relaxation model (1.2), providing a comparison point for the velocity convergence result.","marker":"[26]"}],"fun_headline_variants":["Hyperbolic relaxation provably hits Navier-Stokes limit","Two-parameter relaxation achieves Navier-Stokes convergence","Pressure and velocity converge in hyperbolic relaxation model","Hyperbolic two-parameter system proven to reach Navier-Stokes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The weakest link is the combination of the parameter relation $\\delta \\leq C\\sqrt{\\varepsilon}$ with the well-prepared initial data conditions (in particular $\\|\\nabla\\,\\mathrm{div}\\,u_0^{\\varepsilon,\\delta}\\|_{L^2} \\leq C(\\varepsilon+\\delta)$ and $\\delta(\\|\\Delta p_0^{\\varepsilon,\\delta}\\|_{L^2}+\\|\\nabla\\,\\mathrm{div}\\,U_0^{\\varepsilon,\\delta}\\|_{L^2}) \\leq C(\\varepsilon+\\delta)$); if either fails, the proof gives no convergence of the pressure.","fun_headline_variants_meta":{"raw":{"variants":["Hyperbolic relaxation provably hits Navier-Stokes limit","Two-parameter relaxation achieves Navier-Stokes convergence","Pressure and velocity converge in hyperbolic relaxation model","Hyperbolic two-parameter system proven to reach Navier-Stokes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001191,"raw_usage":{"total_tokens":4894,"prompt_tokens":906,"completion_tokens":3988,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":522,"completion_tokens_details":{"reasoning_tokens":3923}},"tokens_in":522,"tokens_out":3988,"duration_ms":26228,"temperature":1.0,"reasoning_tokens":3923,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:10:23.195390+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the approximate system numerically on the two-dimensional periodic square for a smooth test flow with $\\delta=\\sqrt{\\varepsilon}$ (so $\\delta$ is not $o(\\sqrt{\\varepsilon})$) and well-prepared initial data, and measure $\\sup_{0<t<T}\\|p^{\\varepsilon,\\delta}-p^{NS}\\|_{H^1}$ as $\\varepsilon\\to0$. The theorem only guarantees $\\sqrt{\\varepsilon}\\|p^{\\varepsilon,\\delta}-p^{NS}\\|_{H^1}=O(\\varepsilon+\\delta)$, so the unweighted error could diverge like $1/\\sqrt{\\varepsilon}$; observing it instead vanish would show the $\\delta=o(\\sqrt{\\varepsilon})$ threshold is not necessary. Conversely, keeping $\\delta=o(\\sqrt{\\varepsilon})$ but taking data with $\\|\\nabla\\,\\mathrm{div}\\,u_0^{\\varepsilon,\\delta}\\|_{L^2}$ of order one rather than $O(\\varepsilon+\\delta)$ would test whether the well-prepared condition is needed for pressure convergence.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the first-order relaxation system (1.2) whose single-parameter limit is the incompressible Navier-Stokes equations, supplying the baseline velocity convergence rate that the new two-parameter result extends."},{"cited_title":"343, American Mathematical Soc., 2001","cited_arxiv_id":null,"evidence_quote":"Gives the classical artificial compressibility convergence result for velocity and pressure gradient, motivating the AC component of the new system but lacking the error estimates needed for pressure convergence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the well-established local existence theory for quasilinear symmetric hyperbolic systems, used to guarantee existence and uniqueness for (1.3) and to justify the blow-up criterion in Lemma 2.8."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Ladyzhenskaya inequality and the two-dimensional interpolation framework used throughout the energy estimates."},{"cited_title":"259–305, Birkh¨ auser Boston, Boston, MA, 2001","cited_arxiv_id":null,"evidence_quote":"Presents the basic aspects of hyperbolic relaxation systems and a convergence-stability principle for singularly parametrized hyperbolic systems, organizing the two-parameter limit analysis."},{"cited_title":"21, EDP Sciences, 2007, pp","cited_arxiv_id":null,"evidence_quote":"Advances the analysis of the one-parameter relaxation model (1.2), providing a comparison point for the velocity convergence result."}],"review_version":1}