{"id":"ea0a9de3-db02-4660-beee-c24a7e53f44f","arxiv_id":"2411.17095","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"FV-PINN solves steady incompressible Navier-Stokes with a finite-volume, divergence-theorem loss that lowers the automatic-differentiation order from third to second derivatives of the stream function.","lead":"A neural network solver called FV-PINN solves steady fluid flows by checking fluxes across cell boundaries instead of pointwise equation residuals. On two 2D benchmarks it reports better accuracy and 40 percent faster training than a standard PINN, but gives no quantitative error metrics.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 16 defines the physical loss as a raw sum of two identical x-momentum flux brackets, with no y-momentum term and no squaring, contradicting the MSE definition in Eq. 13. The actual optimized objective is therefore unspecified.","rationale":"The reader's weakest assumption is on target but can be sharpened. The ambiguity is not just whether residuals are squared; as printed, Eq. 16 omits the y-momentum equation entirely because the second bracket repeats the x-momentum expression. This makes the published definition of L_G internally inconsistent with the finite-volume formulation in Eqs. 9-10. For the central claim to be true, the code must have used a corrected, squared objective; the paper neither states this nor releases code. I do not regard this as evidence of fraud—it is a documentation failure, likely a typo—but it is load-bearing because the equivalence between the trained network and the weak-form NS equations is the entire basis for the method. A concrete reimplementation with the literal and corrected losses would settle it. Other weaknesses (no error metrics, sketchy baseline PINN) are real but secondary: they affect quantification of the improvement, not the definition of the method. Hence the verdict remains CONDITIONAL pending correction and a clarifying reproducibility check.","tokens_in":7483,"tokens_out":8293,"duration_ms":78623,"concrete_test":"Re-run the pipe-bend benchmark with three variants using identical architecture and optimizer: (A) the literal Eq. 16 (raw sum, duplicated x-bracket), (B) the corrected squared loss with distinct x- and y-momentum flux residuals, and (C) the corrected raw-sum loss without squares. Compare all three against the COMSOL reference with relative L2 errors in velocity and pressure. If variant A—or C—does not reproduce the reported figure and accuracy, then the paper's stated loss is not the one that produced the results; the central claim depends on an undocumented correction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that FV-PINN enforces the steady incompressible Navier-Stokes equations in weak form with reduced derivative order—requires that the loss actually minimized is a sum of squared finite-volume flux residuals. Eq. 13 defines an MSE, but Eq. 16 (the stated L_G) is a raw sum over cells of the boundary integral of (I1+I2), with no square. Moreover, the two bracketed integrands printed in Eq. 16 are identical: both write u1(u1 n1+u2 n2) - v(∂u1/∂x1 n1 + ∂u1/∂x2 n2) + p n1. The y-momentum term from Eq. 10, u2(u1 n1+u2 n2) - v(∂u2/∂x1 n1 + ∂u2/∂x2 n2) + p n2, never appears. If a reader implements Eq. 16 literally, the network is trained with two copies of x-momentum and no y-momentum; with raw summation, residual cancellation across cells can also hide large local errors. The reported success therefore depends on an unstated corrected loss, and no code or data are provided to resolve the ambiguity. This is load-bearing because the entire equivalence between the finite-volume residual and the enforced PDE is asserted through Eq. 16.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes FV-PINN, a physics-informed neural network for steady incompressible laminar flow. Instead of evaluating strong-form Navier-Stokes residuals at collocation points, the method applies the divergence theorem to convert momentum equation residuals over finite-volume cells into boundary flux integrals, which are then evaluated at Gaussian quadrature points on cell boundaries. A stream function output enforces mass conservation, and the total loss combines a physical loss with a boundary-condition loss. The method is tested on a pipe-bend problem and a double-pipe problem and compared qualitatively with COMSOL, with one comparison against a traditional strong-form PINN. The paper claims improved prediction accuracy, faster convergence, and a 40% reduction in training time.","tokens_in":7741,"tokens_out":6740,"duration_ms":59481,"significance":"If the formulation is corrected and the reported gains are reproducible, the paper makes a modest but useful contribution: it demonstrates a finite-volume-style weak-form loss that reduces the automatic-differentiation derivative order for incompressible Navier-Stokes problems and can be trained with fewer sampling points. Validation against COMSOL and FEM is a positive feature, and the use of Gauss's theorem is standard. However, the contribution is incremental relative to existing weak/variational PINN literature (Refs. 15-19), and the current manuscript does not provide quantitative accuracy metrics, code, or a precise statement of the optimized loss, so the central claims are not yet substantiated.","major_comments":[{"comment":"The written physical loss is not the MSE defined in Eq. (13). Eq. (13) defines the loss as 1/N multiplied by a sum of squared residuals, but Eq. (16) defines L_G as a sum of per-cell flux integrals weighted by quadrature weights, with no square and no 1/N factor. Minimizing a raw sum can let residuals of opposite signs cancel across cells, so large local violations of the momentum equation could coexist with a small L_G. In addition, the second integrand in the quadrature sum in Eq. (16) repeats the x-momentum expression (with u1, du1/dx1, du1/dx2, p n1) instead of the y-momentum expression from Eq. (10), and the pressure flux terms omit the 1/rho factor carried through in Eq. (6). Because the training objective is the central object of the method, this must be corrected, and the exact loss actually minimized should be stated unambiguously, preferably with code or pseudocode.","section":"§2.4, Eq. (16) and Eq. (13)"},{"comment":"The claim that FV-PINN 'significantly improves the prediction accuracy' is not supported by quantitative evidence. Figures 4 and 6 show contour plots, and the text states that results 'closely resemble' COMSOL, but no error norms (for example, relative L2 error) are reported for velocity or pressure, and no error is reported for the traditional PINN baseline in Figure 4(c). Without quantitative comparisons, the relative accuracy claim cannot be assessed.","section":"§3, Figures 4 and 6"},{"comment":"The training-time comparison is incomplete. Table 2 gives training time and number of sampling points for FV-PINN and traditional PINN but does not report the error level achieved by each method, the final loss values, or whether both runs used the same stopping criteria. The 40% speedup is therefore not an apples-to-apples comparison; it should be presented as an accuracy-versus-training-time curve or as error at matched epochs.","section":"§3.1 and Table 2"},{"comment":"The explanation that pressure differences arise from the arbitrary constant is not consistent with the problem setup, which specifies a zero-pressure outlet (p=0) at the bottom boundary. This Dirichlet condition fixes the pressure level, so the observed pressure mismatch should be quantified and explained rather than attributed to indeterminacy. A likely contributor is the very small BC weight beta=0.0001 in Eq. (15), which makes the pressure boundary condition only weakly enforced.","section":"§3.1, pressure-field discussion"}],"minor_comments":[{"comment":"The method is called both FV-PINN and FVI-PINN in different parts of the manuscript; please use one name consistently throughout.","section":"Abstract, §2.3, Table 1"},{"comment":"The kinematic viscosity is written as mu/rho in Eqs. (9)-(10) but as v in Eq. (16); please define nu=mu/rho and use a single symbol to avoid confusion with the velocity component u1.","section":"Eqs. (9), (10), (16)"},{"comment":"The stream-function relation should be written with explicit components, u1=phi_y and u2=-phi_x, and the sign convention should be checked against the vorticity definition or the momentum equation to ensure consistency.","section":"Eq. (3)"},{"comment":"The text refers to 'Chapter 2' and 'Chapter 3'; these should be 'Section 2' and 'Section 3' in a journal article.","section":"Section headings"},{"comment":"The loss curves in Figure 5 are unlabeled and the axes are not described; please identify which curve corresponds to which method and state whether the plotted quantity is L_PINN or one of its components.","section":"Figure 5"},{"comment":"The manuscript does not include a data or code availability statement; for a methods paper, providing the code or pseudocode for the loss computation is strongly recommended.","section":"Reproducibility"}],"recommendation":"major_revision","confidential_remarks":"The core idea is sound and the application is relevant, but the manuscript as written does not specify the loss function actually minimized, and the accuracy claims are qualitative. If the authors can correct Eq. (16), state the exact MSE form, provide quantitative error tables, and address the pressure-boundary inconsistency, the paper could be acceptable. The absence of code is a particular concern given the typographical issues in Eq. (16)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core idea is sound. Replacing collocation-point residuals with boundary flux integrals via Gauss's theorem, and evaluating those with Gaussian quadrature, is a legitimate finite-volume-style weak form. It cuts the derivative order for velocity from second to first (third to second in the stream function), which is a real practical benefit for AD-based training. The comparison against COMSOL is independent and the two benchmark problems are standard. The 40% training-time reduction is plausible, though not rigorously established.\n\nThe soft spots are real but addressable. The big one is Eq. (16). As printed, it sums raw per-cell fluxes rather than squares, it repeats the x-momentum integrand twice, and it drops the y-momentum term entirely. The pressure flux also lacks the 1/rho factor that appears in Eq. (6). That makes the actual optimized objective undefined. The reader's stress-test claim holds up on reading the paper. This is load-bearing because the whole method is the loss, and without a corrected equation the reported success cannot be reproduced. A serious referee should flag this.\n\nThe paper also gives no quantitative error metrics—only visual comparisons with COMSOL. The baseline PINN is under-described, and beta is hand-picked without sensitivity analysis. No code or data are released, which makes the ambiguity in Eq. (16) harder to resolve.\n\nWhat is genuinely new: using finite-volume boundary quadrature points for PINN training, with the stream-function output, is a close relative of hp-VPINN and wPINNs but not explicitly in the cited literature. The derivative-order reduction is correctly derived. The limitations are mostly presentation and documentation, not a wrong central idea.\n\nThis paper deserves serious peer review after revision. The loss definition must be corrected, quantitative errors reported, and ideally code released. I would not cite it yet in my own work without the corrected loss and numbers.\n\nFor a reading group, it is worth one session to discuss weak-form PINN variants and the perils of writing losses carelessly. Send it back for major revision.","headline":"A useful finite-volume weak-form PINN with a real derivative-order reduction, but the written loss function is ambiguous and the reported gains need quantitative support.","tokens_in":8309,"tokens_out":1276,"would_cite":false,"duration_ms":14566,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proposes FV-PINN, which applies Gauss's theorem to convert the Navier-Stokes momentum equation into boundary flux integrals evaluated at Gaussian quadrature points, reducing the required derivative order of velocity, and…","keywords":["Physics-informed neural networks","Finite volume method","Navier-Stokes equations","Incompressible flow","Steady-state problems","Gaussian quadrature","Weak form","Stream function"],"falsifier":"Train FV-PINN on a case with a known analytical solution (e.g., channel flow) and record the per-cell flux residual from Eq. (16) at convergence; if the signed sum across cells is near zero but the individual cell residuals are not, the loss as written is not equivalent to enforcing the equations cell by cell. Alternatively, re-run the paper's pipe-bend case with the loss defined as the squared sum of per-cell residuals and compare the predicted pressure field; a noticeable change would confirm that the sign convention matters.","tokens_in":7259,"feed_emoji":"🌊","tokens_out":9278,"duration_ms":74277,"temperature":0.7,"pith_summary":"The paper proposes a finite-volume physics-informed neural network (FV-PINN) for steady, incompressible, laminar flow. The key move is to rewrite the Navier-Stokes momentum equation, using Gauss's theorem, as surface flux integrals over control-volume boundaries, and to evaluate these fluxes at Gaussian quadrature points rather than computing PDE residuals at scattered collocation points. This reduces the highest derivative of the network's stream-function output from third order to second order, because the viscous term now requires only first derivatives of velocity. The paper applies FV-PINN to the pipe-bend and double-pipe benchmark problems and reports that it predicts velocity and pressure fields that agree with a commercial finite-volume solver, while outperforming a traditional residual-based PINN in accuracy, sampling efficiency, and training speed. If the method holds up, it gives a simple recipe for lowering the automatic-differentiation cost of PINNs for incompressible flow.","feed_headline":"Finite-volume rewrite cuts derivative order for fluid-flow PINNs","feed_subtitle":"Converting Navier-Stokes residuals to boundary flux integrals lowers automatic-differentiation cost and lifts accuracy.","key_machinery":"The central object is the divergence theorem, applied term-by-term to the steady incompressible Navier-Stokes equations to convert each volume integral into a surface flux integral over the boundary of a control volume. The network is trained by assigning a loss to the sum of convective, pressure, and viscous flux integrals at Gaussian quadrature points on each cell boundary. This is the mechanism that reduces the derivative order: the viscous term becomes an integral of $\\nabla\\mathbf{u}$, so the network's second derivatives (of the stream function) are the highest needed, eliminating third-derivative computations. The stream-function ansatz, which automatically enforces mass conservation, works together with the flux-form loss to define the optimization problem.","core_discovery":"The paper's central claim is that converting the differential Navier-Stokes equations into boundary-flux integrals and using those integrals as the PINN loss yields more accurate velocity and pressure predictions for steady incompressible flows than using pointwise residuals, while needing fewer sampling points and less training time. In the FV-PINN, the network outputs the stream function and pressure; velocity is obtained by first derivatives of the stream function, which automatically satisfies the continuity equation. Gauss's theorem turns the inertial, pressure, and viscous terms of the momentum equation into surface integrals of fluxes across each control volume, and Gaussian quadrature on the cell boundaries discretizes those integrals. The viscous flux contains only first derivatives of velocity, so the highest derivative appearing in the loss is a second derivative of the stream function, in contrast to the third derivative needed when the residual is written in strong form. On the pipe-bend and double-pipe problems, the resulting fields are reported to match the commercial solver closely and to be markedly better than those of the traditional PINN.","pith_inferences":["The order reduction is a consequence of the weak form itself, not of the finite-volume discretization, so the same integration-by-parts step should work for other second-order PDEs solved with PINNs, replacing third derivatives with second derivatives.","The loss definition carries a hidden choice: if the implementation follows Eq. (16) literally and sums raw per-cell flux residuals, positive and negative residuals from neighboring cells could cancel and leave large local errors hidden in a small total. The paper does not state whether the loss is a sum of residuals or a sum of their squares, so code inspection would be needed to settle this.","The reported advantage over the traditional PINN is shown only for two relatively simple steady benchmarks; the claim that FV-PINN generalizes to unsteady flows or evolving fluid-solid interfaces would need separate tests, since the time derivative would reset the derivative-order accounting."],"forward_implications":["Training a PINN for steady incompressible flow no longer requires third derivatives of the network output, which lowers the cost and error accumulation of automatic differentiation.","The finite-volume flux form means the loss is assembled from surface integrals over cells, so a converged solution approximately satisfies momentum conservation in integral form on each cell, not only at scattered points.","On the two benchmark problems, FV-PINN used about 30% fewer sampling points and roughly 40% less training time than the traditional PINN while giving better agreement with the commercial solver.","The same divergence-theorem reduction could be applied to the momentum equations used in density-based fluid topology optimization, where accurate pressure fields are needed at many design iterations."],"supporting_citations":[{"why":"Introduces the PINN framework and the stream-function formulation for the Navier-Stokes equations that FV-PINN builds on.","marker":"[1]"},{"why":"The mixed-variable PINN for incompressible laminar flows is the direct rival formulation that FV-PINN compares against.","marker":"[8]"},{"why":"Documents automatic differentiation in the deep learning library used for training, grounding the claim that high-order derivatives are costly.","marker":"[13]"},{"why":"A PINN library reference that discusses the overhead of high-order automatic differentiation, supporting the motivation for order reduction.","marker":"[14]"},{"why":"Establishes variational PINNs with weak-form residuals, the line of work from which FV-PINN's integral loss descends.","marker":"[16]"},{"why":"The divergence-theorem reference that supplies the mathematical identity converting volume integrals into surface flux integrals.","marker":"[20]"},{"why":"The numerical-analysis reference for Gaussian quadrature, which FV-PINN uses to sample and weight cell-boundary points.","marker":"[23]"},{"why":"Introduces the pipe-bend and double-pipe problems used as benchmark cases for validation.","marker":"[24]"}],"fun_headline_variants":["Flux integrals in PINN loss cut derivative order, lift accuracy","Finite-volume PINN reduces derivative order, speeds training","Gauss's theorem makes PINN loss cheaper and more accurate","Stream-function loss cuts derivative order for flow PINNs","Boundary integrals in PINN loss boost speed and accuracy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The main unstated assumption is that the loss actually sums the squares of the per-cell flux residuals, so that a small total loss means each cell's momentum flux is small; if instead the implementation sums raw signed fluxes, the reported accuracy could hide large local errors through cancellation.","fun_headline_variants_meta":{"raw":{"variants":["Flux integrals in PINN loss cut derivative order, lift accuracy","Finite-volume PINN reduces derivative order, speeds training","Gauss's theorem makes PINN loss cheaper and more accurate","Stream-function loss cuts derivative order for flow PINNs","Boundary integrals in PINN loss boost speed and accuracy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000611,"raw_usage":{"total_tokens":2866,"prompt_tokens":993,"completion_tokens":1873,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":609,"completion_tokens_details":{"reasoning_tokens":1790}},"tokens_in":609,"tokens_out":1873,"duration_ms":11898,"temperature":1.0,"reasoning_tokens":1790,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:30:57.389883+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Train FV-PINN on a case with a known analytical solution (e.g., channel flow) and record the per-cell flux residual from Eq. (16) at convergence; if the signed sum across cells is near zero but the individual cell residuals are not, the loss as written is not equivalent to enforcing the equations cell by cell. Alternatively, re-run the paper's pipe-bend case with the loss defined as the squared sum of per-cell residuals and compare the predicted pressure field; a noticeable change would confirm that the sign convention matters.","supporting_citations":[{"cited_title":"This integration enhances model interpretability and ensure s that predictions are more consistent with physical laws","cited_arxiv_id":null,"evidence_quote":"Introduces the PINN framework and the stream-function formulation for the Navier-Stokes equations that FV-PINN builds on."},{"cited_title":"Tackling the curse of dimensionality with physics -informed neural networks","cited_arxiv_id":null,"evidence_quote":"The mixed-variable PINN for incompressible laminar flows is the direct rival formulation that FV-PINN compares against."},{"cited_title":"On physics-informed neural networks for quantum computers","cited_arxiv_id":null,"evidence_quote":"Documents automatic differentiation in the deep learning library used for training, grounding the claim that high-order derivatives are costly."},{"cited_title":"Physics -informed deep learning for incompressible laminar flows","cited_arxiv_id":null,"evidence_quote":"A PINN library reference that discusses the overhead of high-order automatic differentiation, supporting the motivation for order reduction."},{"cited_title":"Flow over an espresso cup: inferring 3 -D velocity and pressure fields from tomographic background oriented Schlieren via physics -informed neural networks","cited_arxiv_id":null,"evidence_quote":"Establishes variational PINNs with weak-form residuals, the line of work from which FV-PINN's integral loss descends."},{"cited_title":"DeepXDE: A deep learning library for solving differential equations","cited_arxiv_id":null,"evidence_quote":"The divergence-theorem reference that supplies the mathematical identity converting volume integrals into surface flux integrals."},{"cited_title":"MIM: A deep mixed residual method for solving high -order partial differential equations","cited_arxiv_id":null,"evidence_quote":"The numerical-analysis reference for Gaussian quadrature, which FV-PINN uses to sample and weight cell-boundary points."},{"cited_title":"wPINNs: Weak physics informed neural networks for approximating entropy solutions of hyperbolic conservation laws","cited_arxiv_id":null,"evidence_quote":"Introduces the pipe-bend and double-pipe problems used as benchmark cases for validation."}],"review_version":1}