{"id":"00077e2c-75c1-4d14-8b41-ebc2ad1a056e","arxiv_id":"1908.09751","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"The exact construction is a restatement of known irrotational-flow solutions, and the numerical method is a fitted approximation with no proven contraction.","lead":"This paper claims to solve the steady incompressible Navier-Stokes equations through a linear system when the external force is a gradient, and proposes approximate generalized method of lines formulas for two-dimensional flows. The exact part reduces to classical irrotational flow theory, while the numerical part is a curve fit with no convergence proof.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 2.1's pressure recovery is unjustified: a curl-free momentum residual on an annulus need not produce a single-valued pressure with prescribed trace P0 on Γ1.","rationale":"The reader's weakest_assumption identifies the compatibility of boundary data with an irrotational velocity and the single-valued pressure issue on a doubly connected domain. My stress test agrees with the rejection but focuses more precisely on the proof's step from equation (11) to the existence of a pressure with a prescribed trace. This is load-bearing because Theorem 2.1's conclusion is exactly that a pressure P with the stated boundary condition exists. Even granting the algebra in (9), the argument only establishes that the momentum residual is curl-free, which is insufficient on an annulus and gives no control over the boundary trace. The uniform-flow counterexample is a decisive analytic check: the velocity is irrotational, divergence-free, satisfies (8), and the momentum equations admit only a constant pressure, so the boundary condition P0(θ)=sinθ on Γ1 cannot hold. The numerical method-of-lines section is secondary; it is validated by minimizing the same residual used as the success metric, and the Banach fixed-point step is invoked without a contraction proof. Thus the central advertised exact result is not established, and the reader's REJECT verdict should stand unchanged.","tokens_in":15806,"tokens_out":7528,"duration_ms":75709,"concrete_test":"Set Ω={1<r<2}, f=0, ν=1, and take w0=x, w1=w2=0, so that u=(1,0), v=0. This satisfies (8) with u0=1, v0=0 on Γ. The momentum equations (12) then force P to be constant. Now choose the boundary datum in (6) to be P0(θ)=sinθ on Γ1={r=2}. The theorem asserts that a solution exists with P=P0 on Γ1, but the unique pressure for this velocity is constant, and sinθ is not constant. This counterexample shows that the pressure boundary condition in Theorem 2.1 is not generally satisfiable. If the author intends P0 to be subject to a compatibility condition, the theorem must state it and the proof must verify it; this test isolates exactly where the current proof fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 2.1 asserts that, for a gradient body force h=∇f, solving the linear system (8) and defining u,v through (7) yields an exact solution of the Navier-Stokes system (5) with boundary conditions (6), including P=P0 on Γ1. The proof's key step is equation (11): the momentum residual F=(ν∇²u−u∂xu−v∂yu+∂xf, ν∇²v−u∂xv−v∂yv+∂yf) satisfies ∂yF1=∂xF2. On the domain Ω, which is explicitly an annulus between Γ0 and Γ1, a curl-free field need not be a gradient: it may have nonzero circulation around the inner boundary. The paper never verifies the period condition that would be required for a single-valued pressure. More seriously, even when a gradient potential exists, it is determined only up to an additive constant; the third boundary condition in (6), P=P0 on Γ1, is an extra Dirichlet condition on a function with no remaining freedom. The proof merely states 'we may obtain P which satisfies the concerning boundary condition' after (12), without constructing P or proving compatibility. Additionally, the constraint ∇²w2+2∂xyw1=0 forces the vorticity ∂xv−∂yu=−(∇²w2+2∂xyw1)=0, so the constructed velocity is irrotational; hence the boundary data must satisfy unstated compatibility conditions such as zero circulation. These gaps invalidate the advertised claim that solving the linear system (8) produces an exact steady Navier-Stokes solution for the stated boundary data.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript has two parts. In the first part, Theorem 2.1 claims that for a gradient body force h=(∂x f, ∂y f), any solution of the linear system (8) for potentials w0, w1, w2 yields an exact solution of the steady incompressible Navier-Stokes system with boundary conditions (6), including a prescribed pressure P=P0 on the outer boundary Γ1. The proof shows that the momentum residual is curl-free and asserts that a pressure can then be recovered. In the second part, the paper proposes a semi-discrete 'generalized method of lines' for an approximate system (25), states that the iteration is solved through the Banach fixed point theorem, presents explicit line formulas for u, v, and P, and reports two numerical examples in which coefficients are obtained by minimizing the residual functional J in (29).","tokens_in":16206,"tokens_out":6975,"duration_ms":67846,"significance":"If Theorem 2.1 were correct, it would be a surprising exact reduction of a nonlinear PDE system to a linear one for gradient forces, and the explicit line formulas could be of practical interest. The algebraic identity in equations (9)-(11) is checkable, and the manuscript is transparent that the numerical coefficients come from minimizing J. However, the central theorem is not proven as stated: the pressure recovery on an annulus is unjustified, the constructed velocity is forced to be irrotational without stating compatibility conditions on the boundary data, and no existence proof is given for the potentials. The numerical evidence is also circular because the same minimized functional J is quoted as evidence of accuracy. The paper does not supply machine-checkable proofs, reproducible code, or a falsifiable prediction that would compensate for these gaps.","major_comments":[{"comment":"Equation (11) establishes only that the momentum residual F is curl-free. On the doubly connected domain Ω, a curl-free field need not be a gradient: one must also verify zero circulation around the inner boundary Γ0. Even if F = ∇P, the potential P is determined only up to an additive constant, so the additional condition P = P0 on the whole of Γ1 is not automatically satisfiable; it requires a compatibility condition. The sentence 'we may obtain P which satisfies the concerning boundary condition' after (12) is unsupported, so the central claim that solving (8) yields a solution of (5)-(6) is not established.","section":"Theorem 2.1, proof after Eq. (11)"},{"comment":"The constraint ∇²w2 + 2∂xyw1 = 0 forces the constructed velocity to be irrotational, since ∂xv − ∂yu = −(∇²w2 + 2∂xyw1) = 0. The boundary data u0, v0 in (6) are therefore not arbitrary; they must be compatible with an irrotational, divergence-free extension, and the theorem states no such compatibility condition. Moreover, the proof begins by assuming w0, w1, w2 exist and never gives an existence argument for the system (8) with general boundary data, so the theorem is incomplete even as a conditional statement.","section":"Theorem 2.1, Eqs. (7)-(8)"},{"comment":"The paper states that the equations are solved through the Banach fixed point theorem, but it does not specify the complete metric space, the norm, or the contraction estimate for the maps T̂n. The displayed iteration is merely a Picard iteration; without a contraction proof it does not establish convergence. This gap affects the method-of-lines existence claim advertised in the abstract and the conclusion.","section":"Section 3, fixed-point iteration after Eq. (24)"},{"comment":"The coefficients {ai[n]}, {bi[n]}, {ci[n]} are obtained by minimizing J, and the same minimized J is then reported as evidence that the approximations are good ('it seems we have got good first approximations'). This is circular: a small J only shows that the ansatz (26)-(28) can fit the two chosen examples, not that the generalized method of lines predicts the Navier-Stokes solution. The general line expressions are not derived from the fixed-point scheme; they are fitting templates whose parameters are chosen to minimize the residual.","section":"Section 3.1, Eq. (29)"}],"minor_comments":[{"comment":"The phrase 'P is a solution of the system indicated in the first two lines of (25)' is confusing because (25) is introduced later and contains the regularization parameter ε; the theorem should state the momentum equations explicitly.","section":"Theorem 2.1, statement"},{"comment":"There are subscript and notation typos, such as 'd̂1(v,vn−1)' in (23), which appears to be 'd̂1(vn, vn−1)', and similar omissions in the displayed line formulas.","section":"Eq. (23) and line formulas"},{"comment":"The text says that in the numerical examples there are no boundary conditions for the pressure, but it also states that P0 must be calculated numerically in the optimization process; the role of P0 in (26)-(28) needs clarification.","section":"Section 3.1, boundary conditions"},{"comment":"The figures would be more informative with labeled axes, since only the x-axis unit convention is supplied in the captions and no comparison with a known solution is shown.","section":"Figures 1-12"},{"comment":"In the v9 expression, the coefficients −0.057 f6 u0 v0′ and −0.057 f8 v0 v0′ appear inconsistent with the monotone coefficient progression in neighboring lines; this may be a typographical error.","section":"Line 9, v9 formula"}],"recommendation":"reject","confidential_remarks":"The bibliography consists largely of the author's own prior work, and the presentation would need independent grounding of the method-of-lines claims. In view of the unproven central theorem and the circular numerical validation, I do not see a path to acceptance within the scope of the submitted manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the advertised exact result isn't proven, and the numerical part doesn't validate the method. The core observation—an irrotational, incompressible field solves steady Navier-Stokes with a potential force and Bernoulli pressure—is classical. The w0,w1,w2 representation is bookkeeping, not a new reduction. The method-of-lines formulas for annular domains are concrete and might have heuristic value, but the paper doesn't earn that.\n\nThe big hole is Theorem 2.1. Equation (11) only shows the momentum residual is curl-free. On an annulus, that does not imply the existence of a single-valued pressure; you need a zero-circulation condition, which the paper never states or verifies. Even if the residual were a gradient, the pressure is determined only up to a constant, so the extra boundary condition P=P0 on Γ1 cannot be imposed arbitrarily. The proof's sentence \"we may obtain P which satisfies the concerning boundary condition\" is exactly where the argument stops. The theorem also forces zero vorticity, so the boundary data are not arbitrary; none of this is discussed.\n\nThe numerical section has a separate, load-bearing issue. The coefficients in (26)-(28) are obtained by minimizing the residual functional J in (29), and then the smallness of J is reported as evidence that the approximations are good. That is circular: a fitted residual will be small by construction. It becomes evidence only when tested on independent data. There is also no proof of the contraction needed for the Banach fixed point step, so the convergence of the line iteration is asserted, not shown.\n\nWhat's salvageable: the line expressions might be useful as a cheap approximate solver for a narrow class of irrotational, incompressible flows in annuli, and the coordinate transformation to (t,θ) is competently written. But the central claims are unsubstantiated. I don't think this deserves referee time; a desk reject is appropriate. If the author fixes the topology/pressure issues and validates the numerical scheme on a problem with a known solution, it could become a modest contribution.","headline":"The exact theorem has a genuine topological gap (curl-free doesn't mean gradient on an annulus, and the pressure trace is overdetermined), and the numerical validation is a fit to the same residual it claims to confirm.","tokens_in":16660,"tokens_out":4513,"would_cite":false,"duration_ms":46735,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q30","76D05","65N40"],"pacs":[],"model":"deepseek-v4-flash","headline":"Solving a linear system yields exact steady Navier-Stokes solutions when the body force is a gradient.","keywords":["steady incompressible Navier-Stokes","gradient body force","three-potential ansatz","generalized method of lines","Banach fixed point theorem","potential flow","vorticity-free flow"],"falsifier":"Solve (8) with boundary data $u_0=-y$, $v_0=x$ on the outer circle and no-slip on the inner circle. Any field produced by (7) has zero vorticity, while the prescribed outer trace has non-zero circulation, so either the linear system has no solution or the resulting velocity fails the momentum equation; this would settle whether Theorem 2.1 holds for arbitrary boundary data.","tokens_in":15611,"feed_emoji":"🌀","tokens_out":15065,"duration_ms":144388,"temperature":0.7,"pith_summary":"This paper aims to show that, for the steady incompressible Navier-Stokes equations in two dimensions with a gradient body force, an exact solution can be obtained by solving a linear system for three auxiliary potentials rather than the original nonlinear system. The intended payoff is that a nonlinear fluid problem becomes, in this special case, a linear elliptic problem plus an integration for the pressure. The paper also develops approximate solutions through the generalized method of lines, which discretizes the domain into radial lines and produces line-by-line formulas for velocity and pressure, with coefficients fixed by minimizing the equation residual. If the central construction is right, it gives a direct route to exact test solutions for steady Navier-Stokes flows driven by conservative forces.","feed_headline":"Three-potential trick turns steady Navier-Stokes into a linear system","feed_subtitle":"For gradient body forces, a three-potential trick makes the nonlinear terms curl-free, leaving only linear equations.","key_machinery":"The central object is the three-potential velocity ansatz $u=\\partial_x w_0+\\partial_x w_1+\\partial_y w_2$, $v=\\partial_y w_0-\\partial_y w_1-\\partial_x w_2$. The two equations in (8) serve two purposes: $\\nabla^2 w_0+\\partial_{xx}w_1-\\partial_{yy}w_1=0$ enforces incompressibility, and $\\nabla^2 w_2+2\\partial_{xy}w_1=0$ makes the vorticity of $(u,v)$ vanish. Together they make the nonlinear convective terms curl-free, so the momentum equations reduce to a gradient and can be integrated to a pressure. The generalized method of lines discretizes one coordinate, here the radial variable, into a set of lines, writes the solution on each line as a function of the boundary data and the line index, and solves the resulting coupled system by a Banach fixed-point iteration.","core_discovery":"Theorem 2.1 claims that for a gradient body force $(\\partial_x f,\\partial_y f)$, the velocity field defined by $u=\\partial_x w_0+\\partial_x w_1+\\partial_y w_2$ and $v=\\partial_y w_0-\\partial_y w_1-\\partial_x w_2$ solves the steady incompressible Navier-Stokes system whenever $w_0,w_1,w_2$ satisfy the linear system $\\nabla^2 w_2+2\\partial_{xy}w_1=0$, $\\nabla^2 w_0+\\partial_{xx}w_1-\\partial_{yy}w_1=0$, together with the boundary conditions on $u$ and $v$. The pressure $P$ is then obtained by integrating the curl-free momentum equations with $P=P_0$ on the outer boundary. The load-bearing identity is (9): when the two auxiliary equations hold, the curl of the convective terms vanishes, so the momentum operator is curl-free and a scalar pressure exists. A direct consequence is that the resulting flow has zero vorticity. The paper's second part extends this to approximate solutions, actually for the regularized system (25) with a small parameter $\\varepsilon$, by the generalized method of lines, using a Banach fixed-point iteration on each line and minimizing the $L^2$ residual (29).","pith_inferences":["The theorem's exact solutions are necessarily potential flows: the first equation of (8) forces zero vorticity, so the boundary data admitted by the construction are only those compatible with an irrotational velocity field, even though the theorem states the boundary conditions without this restriction.","On a doubly connected domain, recovering a single-valued pressure from a curl-free momentum equation requires a compatibility condition around the inner boundary; checking whether the pressure returns to its starting value after one loop would provide a direct test of the construction.","The same three-potential trick could turn other conservative-force fluid models, such as steady magnetohydrodynamics with a gradient magnetic-pressure term, into linear problems if the additional force can be written as a gradient."],"forward_implications":["If Theorem 2.1 is correct, an exact steady incompressible Navier-Stokes solution follows from solving the linear system (8) for any gradient body force, without iterating on the nonlinear terms.","The method-of-lines formulas, with coefficients obtained by minimizing the residual $J$, give explicit approximate velocity and pressure fields on each line that can serve as initial guesses or manufactured solutions for Navier-Stokes computations.","Because the exact constructed flows are irrotational, they form a potential-flow subclass of the steady Navier-Stokes solutions; their pressure is recovered from a gradient, so the boundary pressure datum determines the integration constant.","The paper announces the extension to three-dimensional, compressible, and time-dependent cases as future work."],"supporting_citations":[{"why":"introduced the generalized method of lines and the earlier two-potential velocity ansatz that Theorem 2.1 extends.","marker":"[2]"},{"why":"presented the earlier Navier-Stokes solution ansatz that the new three-potential form complements.","marker":"[3]"},{"why":"developed the generalized method of lines for related PDE systems, giving the line-discretization template used in Section 3.","marker":"[4]"},{"why":"is cited for similar exact results for Euler and Navier-Stokes equations and for the method-of-lines framework.","marker":"[5]"},{"why":"supplies the hyper-finite-difference and software approach used to truncate the line series to order d^2.","marker":"[6]"},{"why":"underpins the finite-difference discretization used to convert the PDE system into the line equations.","marker":"[7]"}],"fun_headline_variants":["Three-potential trick linearizes steady Navier-Stokes","Gradient forces yield linear Navier-Stokes system","Method of lines meets Banach for Navier-Stokes","Steady Navier-Stokes tamed by three-potential identity","Line discretization gives Navier-Stokes approximations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The boundary data must be traces of a curl-free velocity field with a single-valued pressure around the hole; the paper states the boundary conditions without verifying this compatibility.","fun_headline_variants_meta":{"raw":{"variants":["Three-potential trick linearizes steady Navier-Stokes","Gradient forces yield linear Navier-Stokes system","Method of lines meets Banach for Navier-Stokes","Steady Navier-Stokes tamed by three-potential identity","Line discretization gives Navier-Stokes approximations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000608,"raw_usage":{"total_tokens":2829,"prompt_tokens":939,"completion_tokens":1890,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":1810}},"tokens_in":555,"tokens_out":1890,"duration_ms":15475,"temperature":1.0,"reasoning_tokens":1810,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:56:51.393057+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve (8) with boundary data $u_0=-y$, $v_0=x$ on the outer circle and no-slip on the inner circle. Any field produced by (7) has zero vorticity, while the prescribed outer trace has non-zero circulation, so either the linear system has no solution or the resulting velocity fails the momentum equation; this would settle whether Theorem 2.1 holds for arbitrary boundary data.","supporting_citations":[{"cited_title":"Botelho, Topics on Functional Analysis, Calculus of V ariations and Duality, Academic Publications, Sofia, (2011)","cited_arxiv_id":null,"evidence_quote":"introduced the generalized method of lines and the earlier two-potential velocity ansatz that Theorem 2.1 extends."},{"cited_title":"Botelho, Variational Convex Analysis, Lambert Acade mic Publishing, Berlin, June 2010","cited_arxiv_id":null,"evidence_quote":"presented the earlier Navier-Stokes solution ansatz that the new three-potential form complements."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"developed the generalized method of lines for related PDE systems, giving the line-discretization template used in Section 3."},{"cited_title":"Botelho, Functional Analysis and Applied Optimizati on in Banach Spaces, Springer Switzerland, 2014","cited_arxiv_id":null,"evidence_quote":"is cited for similar exact results for Euler and Navier-Stokes equations and for the method-of-lines framework."},{"cited_title":"On the generalized method of lines and its proximal explicit and hyper-finite difference approaches","cited_arxiv_id":"1904.12379","evidence_quote":"supplies the hyper-finite-difference and software approach used to truncate the line series to order d^2."},{"cited_title":"Strikwerda, Finite Difference Schemes and Partial Differential Equation s, SIAM, sec- ond edition (2004)","cited_arxiv_id":null,"evidence_quote":"underpins the finite-difference discretization used to convert the PDE system into the line equations."}],"review_version":1}