{"id":"e0276d50-d65f-48df-96af-7766bd2dbc52","arxiv_id":"2607.17036","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"A fully linear ETD-mr-SAV scheme for the 2D NSE is proposed with an advertised uniform-in-time stability theorem, but the theorem's key energy identity is algebraically wrong, so the main claim is unproven.","lead":"This paper presents a fully linear variable-step exponential time-differencing scheme for the 2D Navier-Stokes equations in vorticity form, with an embedded pair for adaptive stepping and a claimed unconditional long-time L2 stability bound. The stability proof, however, rests on an energy identity that fails because the coupled nonlinear terms do not cancel.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Energy identity (28) in Theorem 1 omits the 2q⟨B(ω~^{n+1/2},ω~^{n+1/2}),ω⟩ cross term arising from (4a)–(4b); since ⟨B(v,v),v⟩=0 does not apply when the test argument is ω≠ω~, the stated cancellation fails and the uniform-in-time bound is unproved.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing flaw: the purported cancellation in Eq. (28) is algebraically false. I rechecked the signs and arguments: both (4a) and (4b) contain the same extrapolated nonlinearity B(ω~,ω~) with the same positive sign when forming the q-weighted sum, so the cross terms add rather than cancel. The continuous conservation identity ⟨B(ω,ω),ω⟩=0 cannot be invoked because the second slot in the inner product is the current ω(t), not the extrapolated ω~ used in B. This is not a matter of missing technical regularity; it is an elementary algebraic error in the central theorem. Since Theorem 1's uniform-in-time bound is the paper's main advertised contribution, and the numerics, while promising, do not substitute for the proof, the reader's REJECT verdict is appropriate. I do not see a different, stronger concern that would change this assessment; the quadrature-vs-exact-scheme gap is real but secondary. No ad hominem is intended; the issue is purely in the argument.","tokens_in":22070,"tokens_out":3801,"duration_ms":35389,"concrete_test":"Independently re-derive the sum of ⟨(4a), ω⟩ and q·(4b) without invoking any cancellation. Verify that the term 2q⟨B(ω~^{n+1/2},ω~^{n+1/2}),ω⟩ persists and is generally nonzero by testing a simple periodic pair, e.g., ω=cos x+cos y and ω~=cos x−cos y, evaluating the inner product explicitly. If the cross term is nonzero, Eq. (28) is false and the proof of Theorem 1 fails; no numerical experiment can rescue the analytic claim without a new mechanism to control or cancel this term.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proof of Theorem 1 rests entirely on Eq. (28), which claims that summing the ω-inner product of (4a) with q times (4b) cancels the coupled nonlinear terms. But (4a) contributes +q⟨B(ω~^{n+1/2},ω~^{n+1/2}),ω⟩, while (4b) multiplied by q contributes +q⟨B(ω~^{n+1/2},ω~^{n+1/2}),ω⟩. The sum therefore contains 2q⟨B(ω~^{n+1/2},ω~^{n+1/2}),ω⟩, not zero. The only available cancellation, ⟨B(v,v),v⟩=0, requires the second argument in the inner product to be the same v as in the nonlinearity; here the second argument is ω(t), while the nonlinearity is evaluated at the extrapolated value ω~^{n+1/2}, and these differ generically. Thus Eq. (28) does not follow, and the subsequent differential inequality (29), the integrated bound (32), and the stated L∞(0,∞;L2) estimate (26) are unsupported. This is the load-bearing step: without (28) the energy estimate does not close, and the core claim of unconditional long-time stability 'for all Reynolds numbers and time-step sizes' is not established. A secondary gap—the theorem analyzes the exact subproblem (4) with continuous q, while the practical scheme uses Talbot inversion and Gauss–Lobatto quadrature—would also need bridging, but the algebraic failure in (28) already invalidates the proof as written.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a fully linear, variable-step embedded exponential time-differencing scheme for the 2D Navier–Stokes equations in vorticity–streamfunction form. The second-order scheme ETD-mr-SAV-MS2-L and an embedded first-order companion are designed so that each time step requires only linear subproblems: two heat-equation solves and a scalar auxiliary-variable equation treated by Laplace transform and Talbot inversion. The central theoretical claim is Theorem 1: under uniformly bounded L2 forcing and for any mean-reverting parameter γ>0, the discrete vorticity satisfies a uniform-in-time L∞(0,∞;L2) bound for all Reynolds numbers and all time-step sizes. Numerical experiments report second-order accuracy under fixed and perturbed variable steps, long-time boundedness, statistical consistency, and adaptive/VSVO performance.","tokens_in":22590,"tokens_out":4732,"duration_ms":47020,"significance":"If valid, the main result would be a meaningful advance: a fully linear, uniquely solvable, variable-step ETD scheme with unconditional long-time stability for the forced 2D NSE would improve on earlier mr-SAV/ETD schemes that require nonlinear scalar solves. The numerical study is broad and carefully presented, including long-time statistics and a VSVO variant, and the computational claims are plausible. However, the proof of the headline stability theorem rests on an algebraic cancellation that does not hold for the actual scheme. Since that theorem is the paper's principal theoretical contribution, the present version does not establish the claimed stability result.","major_comments":[{"comment":"The proof of Theorem 1 is invalid. Taking the H-inner product of (4a) with ω and adding q times (4b) gives, on the left-hand side, (1/2)d/dt(||ω||^2+|q|^2) + ν||L^{1/2}ω||^2 + γ|q|^2 plus the two coupled nonlinear terms: +q⟨B(ω~^{n+1/2},ω~^{n+1/2}),ω⟩ from (4a) and +q⟨B(ω~^{n+1/2},ω~^{n+1/2}),ω⟩ from (4b) multiplied by q. The sum contains 2q⟨B(ω~^{n+1/2},ω~^{n+1/2}),ω⟩, not zero. The identity ⟨B(v,v),v⟩=0 would require the second argument in the inner product to equal the argument of B, but here the second argument is ω(t), while B is evaluated at the extrapolated value ω~^{n+1/2}; these are generically different. Thus Eq. (28) does not follow, and the subsequent differential inequality (29), the integrated bound (32), the uniform bound (26), and the cumulative dissipation estimate (27) are unsupported. This is the load-bearing step of the paper.","section":"Section 4.1, Eq. (28)"}],"minor_comments":[],"recommendation":"reject","confidential_remarks":"The paper contains an extensive and interesting numerical study, but the central stability theorem is not proved: the cancellation asserted in Eq. (28) is algebraically false for the stated scheme. Because the main claim of the paper rises or falls on this theorem, I recommend rejection. A substantially revised manuscript that either proves the bound by a correct argument or modifies the scheme so that the cancellation is genuine could be reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a real step forward in scheme design, but the main theoretical claim does not hold as written. The paper's key theorem (Theorem 1) asserts a uniform-in-time L2 bound for arbitrary step sizes and Reynolds numbers. The proof depends on Eq. (28), an energy identity that claims the coupled nonlinear terms cancel. They don't. Taking the omega-inner product of (4a) gives +q<B(omega-tilde,omega-tilde),omega>; multiplying (4b) by q gives +q<B(omega-tilde,omega-tilde),omega>. Adding gives 2q times that inner product, not zero. The usual cancellation <B(v,v),v>=0 does not apply because the test argument omega is not the same as the extrapolated argument omega-tilde. So the differential inequality (29), the iteration (32), and the theorem are unsupported. This is not a cosmetic remark; without that cancellation the energy estimate does not close.\n\nWhat is genuinely good: the scheme itself is new. Replacing the earlier nonlinear mr-SAV solve [21] with a Laplace-transform/Talbot inversion for the scalar auxiliary variable gives a fully linear ETD-mr-SAV scheme, and the paper argues convincingly that only two FFTs are needed per step. The numerical work is extensive: fixed and variable step convergence, long-time Kolmogorov runs, adaptive step control, and statistical diagnostics with total variation distances. The numerics are consistent with the advertised behavior and the paper reads honestly; they even report CPU overhead.\n\nSoft spots beyond the proof gap: the theorem analyzes the exact subproblem (4) where q(t) is solved continuously, while the implementation replaces it with Talbot inversion and Gauss-Lobatto quadrature. That gap would also need closing if the stability proof were fixed. Also no convergence proof — the paper lists it as future work, which is fine, but the headline result is the stability.\n\nBottom line: the method and experiments are worth a careful look, and the Laplace/Talbot scalar solve is a useful trick. But as it stands, the central theorem is unsupported. This paper should go to peer review — a referee might help the authors repair the proof or restate a correct claim — but not be accepted on the strength of the current analysis. Cite only with caution.","headline":"The fully linear Laplace-based ETD-mr-SAV scheme is a genuine new construction with strong numerics, but the advertised uniform-in-time stability theorem is unproved: the energy identity (28) drops a sign-definite cross term.","tokens_in":23022,"tokens_out":2431,"would_cite":false,"duration_ms":22713,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M12","65M70","76D05"],"pacs":[],"model":"deepseek-v4-flash","headline":"A fully linear variable-step scheme claims unconditional long-time stability for the 2D Navier–Stokes equations: with uniformly L2-bounded forcing, the discrete vorticity remains bounded for all Reynolds numbers and all time-step sizes.","keywords":["Navier-Stokes equations","exponential time-differencing","mean-reverting scalar auxiliary variable","variable-step","adaptive time stepping","long-time stability","vorticity-streamfunction","Fourier spectral method"],"falsifier":"Compute the left side of Eq. (28) directly: test (4a) with ω(t), multiply (4b) by q(t), and add; keep the term 2q⟨B(ω̃,ω̃),ω⟩. If this term is nonzero, the claimed differential inequality (29) — and therefore the uniform bound — is not a consequence of the scheme as written.","tokens_in":21962,"feed_emoji":"🌊","tokens_out":7236,"duration_ms":58034,"temperature":0.7,"pith_summary":"The paper proposes a second-order exponential time-differencing scheme, paired with an embedded first-order companion, for the forced 2D Navier–Stokes equations in vorticity–streamfunction form on a periodic box. Its central assertion is that the scheme is unconditionally long-time stable: whenever the forcing is uniformly bounded in L2, the discrete vorticity stays uniformly bounded in L2 for every Reynolds number and every time-step size, with no CFL-type restriction. This matters because long-time simulation of fluids needs methods that inherit the equations' dissipative structure, and most higher-order linear schemes lack rigorous uniform-in-time bounds. The design keeps every step linear — two heat solves and one scalar auxiliary-variable equation solved by Laplace transform — so the stability guarantee does not come at the price of nonlinear solves.","feed_headline":"Linear adaptive scheme bounds 2D vortex flow at any step size","feed_subtitle":"Uniform-in-time stability for forced Navier–Stokes: bounded vorticity at all Reynolds numbers and step sizes.","key_machinery":"The carrying mechanism is the ETD-mr-SAV-MS2-L scheme itself: an exponential time-differencing update for the linear viscous operator, second-order midpoint extrapolation of the nonlinear advection, and a mean-reverting scalar auxiliary variable q whose equilibrium value is 1. The mr-SAV term adds a dissipative correction γ(q−1) that drives q back to equilibrium and, in the energy identity, is intended to control the non-closed discrete advection contribution. Each step requires solving two linear heat equations in Fourier space and one linear scalar convolution equation for q, evaluated by Laplace transform and Talbot's numerical inverse transform, so the whole update is uniquely solvable a","core_discovery":"The core claim is Theorem 1: for the ETD-mr-SAV-MS2-L scheme, if the forcing f lies in L∞(0,∞;H), then the discrete enstrophy obeys ||ω^{n+1}||^2 + |q^{n+1}|^2 ≤ e^{-θΣτ}(||ω^1||^2+|q^1|^2) + (1/θ)(||f||^2/(λ1ν)+γ) with θ = min{νλ1,γ} > 0, giving a bound in L∞(0,∞;L2) independent of Reynolds number and step size. The argument tests the vorticity equation against ω and the auxiliary-variable equation against q; the author asserts that the coupling terms cancel, leaving a differential inequality that is then iterated over arbitrarily many steps. The mean-reverting parameter γ is what supplies the extra dissipation when the discrete advection term fails to vanish exactly, and a cumulative dissi","pith_inferences":["Beyond the paper: if the claimed cancellation in the enstrophy identity in fact fails, a corrected scheme with an altered coupling sign or a modified auxiliary equation might still yield a similar uniform bound — but that recovery is not what this paper establishes.","Beyond the paper: the Laplace-transform/Talbot treatment of the auxiliary equation is specific to periodic settings with diagonalizable viscous operators; extending the approach to bounded domains would require an alternative solver for the scalar equation.","Beyond the paper: the mean-reverting parameter γ acts like a tunable numerical friction; the experiments suggest a trade-off between accuracy and stabilization, since small γ drifts from the reference trajectory while large γ stays close, implying γ could be selected adaptively.","Beyond the paper: the same construction should transfer to other dissipative systems with an exactly conserved nonlinearity, such as Boussinesq-type models, since only the identity ⟨B(u,u),u⟩=0 is used."],"forward_implications":["A practitioner can run the scheme to arbitrarily long times with a guaranteed uniform L2 bound on vorticity for any Reynolds number and any step size, provided the forcing is L2-bounded.","The embedded first-order companion gives a posteriori error control at negligible cost, so adaptive step selection can be used without sacrificing the long-time stability guarantee.","Because all subproblems are linear, the method sidesteps nonlinear algebraic solves that other mr-SAV–ETD schemes require, while retaining second-order accuracy.","The uniform-in-time estimate supports reliable computation of long-time statistics such as enstrophy distributions, as demonstrated by small total-variation distances between PDFs from fixed and adaptive runs."],"fun_headline_variants":["Linear adaptive ETD scheme with Reynolds-independent stability","2D vortex flow bounded for any step size and Reynolds number","Fully linear method with uniform-in-time stability for Navier-Stokes","Adaptive ETD keeps vorticity bounded at all step sizes","New linear scheme guarantees long-time stability for 2D flows"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The proof of the uniform bound depends on the two coupling terms between the vorticity and auxiliary equations cancelling exactly when the enstrophy identity is formed; if they do not cancel, the dissipative inequality and the bound do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Linear adaptive ETD scheme with Reynolds-independent stability","2D vortex flow bounded for any step size and Reynolds number","Fully linear method with uniform-in-time stability for Navier-Stokes","Adaptive ETD keeps vorticity bounded at all step sizes","New linear scheme guarantees long-time stability for 2D flows"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000225,"raw_usage":{"total_tokens":1331,"prompt_tokens":804,"completion_tokens":527,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":548,"completion_tokens_details":{"reasoning_tokens":456}},"tokens_in":548,"tokens_out":527,"duration_ms":5047,"temperature":1.0,"reasoning_tokens":456,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T19:13:50.661515+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the left side of Eq. (28) directly: test (4a) with ω(t), multiply (4b) by q(t), and add; keep the term 2q⟨B(ω̃,ω̃),ω⟩. If this term is nonzero, the claimed differential inequality (29) — and therefore the uniform bound — is not a consequence of the scheme as written.","supporting_citations":[],"review_version":1}