{"id":"16874b8c-49e2-49ff-b580-39dc0eabbe07","arxiv_id":"2603.27714","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Divergence-free BDM fields on triangulated surfaces split L2-orthogonally into streamfunction rotations plus a Betti-number-dimensional harmonic space, enabling pressure-free high-order surface Navier-Stokes.","lead":"A discrete Helmholtz-Hodge split for high-order BDM elements on surfaces of any topology turns incompressible surface flow into a streamfunction plus a few harmonic coefficients, removing pressure entirely. This keeps exact tangentiality, pointwise divergence-free velocity and pressure-robustness while cutting the saddle-point system.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The reader correctly isolates the dimension count as the only potential soft spot, yet that count is elementary, standard in FEEC, and already verified by the authors for both closed and open surfaces. The Piola-mapped polynomial spaces preserve the necessary inclusion rot(S^{k+1}_0)\\subset BDM^k_0 and the divergence relation, so geometric mappings do not alter the combinatorial dimensions used in the argument. Consequently the L2-orthogonal splitting holds, the pressure-free streamfunction-harmonic formulation inherits exact tangentiality, pointwise divergence-freeness and pressure-robustness, and the numerical examples on the trefoil and the pierced sculpture are faithful illustrations rather than mere demonstrations. No adjustment to the ACCEPT verdict is warranted.","tokens_in":22856,"tokens_out":442,"duration_ms":5086,"concrete_test":"Independently recompute the two DOF formulae in the proof of Theorem Appendix B.1 for a concrete closed genus-1 mesh (e.g., a regularly refined torus triangulation) and for a genus-0 surface with four boundary components; verify that the resulting dim(H^k_BDM) equals 2 and 3 respectively for several polynomial degrees k. If the counts match, the topological exactness of the discrete complex is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Theorem 3.2) rests on a dimension-counting argument (Appendix B.1) that dim(H^k_BDM)=b1(M). The DOF tallies for mapped BDM and continuous Lagrange spaces under Piola transforms, together with the Euler-Poincaré formula, are standard and carefully checked; the paper also notes that the lowest-order RT0 component already carries the full harmonic structure (Remark 3.3). No internal inconsistency or missing hypothesis that would invalidate the splitting or the subsequent pressure-free reformulation is apparent. The continuous background, the randomized basis construction, hybridization/Schur treatment, and the numerical illustrations on non-trivial topology are all consistent with the claim.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper establishes a discrete L2-orthogonal Helmholtz–Hodge decomposition for the divergence-free subspace of H(div)-conforming BDM elements of degree k on triangulated surfaces of arbitrary topology: J^k_BDM = rot(S^{k+1}_0) ⊕_L2 H^k_BDM, with dim(H^k_BDM) equal to the first Betti number b1(M) (Theorem 3.2). The proof is by elementary dimension counting that recovers the Euler–Poincaré formula (Appendix B.1). Consequently any incompressible surface flow discretized in this subspace can be rewritten with a continuous streamfunction and finitely many harmonic coefficients as the only unknowns, eliminating the pressure and the saddle-point structure while retaining exact tangentiality, pointwise divergence-freeness and pressure-robustness. A randomized algorithm constructs an L2-orthonormal harmonic basis; hybridization, a Schur-complement treatment of the few harmonic unknowns, and a post-processed pressure reconstruction are described. Numerical experiments for unsteady surface Navier–Stokes on a trefoil knot and a multiply-connected sculpture surface illustrate the method and the physical role of the harmonic component.","tokens_in":23007,"tokens_out":826,"duration_ms":8150,"significance":"If the discrete splitting holds, the work supplies a structure-preserving, pressure-free high-order method for surface Stokes/Navier–Stokes that inherits all geometric exactness properties of the BDM framework while removing the velocity-pressure saddle point. The dimension count is elementary and topology-exact; the lowest-order RT0 component already carries the full harmonic structure (Remark 3.3). The randomized basis construction, hybridization and Schur treatment of the O(b1) harmonic unknowns are practical, and the numerical illustrations on non-trivial topology make the physical content of the harmonic fields transparent. The approach therefore extends classical streamfunction methods from simply-connected flat domains to surfaces of arbitrary genus in a way that is both theoretically clean and implementable.","major_comments":[],"minor_comments":[{"comment":"Section 3.3 and Remark 3.5: the discussion of incomplete decompositions and the obstruction for curved triangulations is interesting but somewhat peripheral; a short forward pointer that the main applications only need the complete splitting of J^k_BDM would help the reader.","section":"Section 3.3"},{"comment":"Table 1: the hybrid-unknown counts appear in gray parentheses; a one-sentence clarification in the caption that these are the condensed facet unknowns would improve readability.","section":"Table 1"},{"comment":"Figures 5 and 7: the color scale is stated in the caption but the absolute magnitude of the harmonic component is hard to judge visually; a brief remark on relative L2 norms of u_rot versus u_H at selected times would strengthen the physical interpretation.","section":"Figures 5, 7"},{"comment":"A few typographical inconsistencies remain (e.g., “hol(e)y” in the caption of Figure 1, occasional missing spaces around operators). A final copy-edit pass would be beneficial.","section":null}],"recommendation":"accept","confidential_remarks":"The manuscript is a clean, self-contained contribution that sits comfortably in the journal’s numerical-analysis scope. The dimension-counting argument is standard and carefully checked; I see no load-bearing technical gap. The only possible novelty question is the precise relation to the fluid-cohomology literature [50,52], but the authors already cite those works and clearly distinguish the high-order BDM setting. I would accept as is."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This paper gives a practical, high-order way to drop the pressure from surface Navier-Stokes while keeping every geometric exactness property of BDM. The core is Theorem 3.2: the divergence-free BDM subspace splits L2-orthogonally into rot of continuous Lagrange streamfunctions plus a harmonic space whose dimension is exactly the first Betti number. That is the new piece. Streamfunction formulations already exist for simply-connected surfaces and for lowest-order DEC; the high-order BDM realization that correctly captures cohomology on arbitrary topology, together with the randomized orthonormal basis construction and the Schur/hybridization treatment of the few harmonic unknowns, is not in the earlier literature.\n\nThe math is elementary and clean. They prove the dimension count by DOF tallies plus Euler-Poincaré (Appendix B.1); the continuous background is taken from their own prior work but used only as scaffolding. Remark 3.3 notes that the lowest-order RT0 component already carries the full harmonic structure, so higher-order moments are unnecessary for topology. Implementation details (hybridization, Schur complement for the tiny harmonic block, pressure reconstruction by residual projection) are concrete and the code is released. The two numerical examples (trefoil knot, pierced sculpture) are chosen precisely to show the physical role of the harmonic component; they are illustrative rather than exhaustive, which is fine given that the velocity solution is equivalent to the earlier velocity-pressure scheme.\n\nSoft spots are minor. The dimension argument assumes standard DOF counts under Piola maps and shape-regular triangulations; that is standard and carefully checked, not a hidden flaw. There is no new convergence theory, only inheritance from the parent BDM method. Incomplete discrete Helmholtz-Hodge remarks are honest about what cannot be characterized as pure gradients. Citation pattern is appropriate; self-citations are background, not circular.\n\nThis is for people who actually write structure-preserving surface FEM codes for thin films or biomembranes. It is a solid incremental method paper inside an active program. I would send it to referees without hesitation; the central claim holds and the engineering is usable.","headline":"Clean high-order BDM Helmholtz-Hodge on surfaces of any topology that removes the pressure saddle-point while keeping exact tangentiality and pointwise div-free structure; solid and ready for referees.","tokens_in":23630,"tokens_out":538,"would_cite":true,"duration_ms":6291,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N30","65N12","76D05","58A14"],"pacs":[],"model":"grok-4.5","headline":"Any incompressible surface flow on BDM elements can be rewritten with only a streamfunction and a few harmonic coefficients, eliminating pressure.","keywords":["surface Navier-Stokes","Helmholtz-Hodge decomposition","streamfunction","harmonic fields","BDM finite elements","divergence-free discretization","pressure-robustness"],"falsifier":"Compute dim(H^k_BDM) by linear algebra on a closed genus-g surface for several polynomial degrees k and mesh sizes; if the dimension is not constantly equal to 2g, the central claim is false.","tokens_in":23750,"feed_emoji":"∿️","tokens_out":707,"duration_ms":7610,"temperature":0.7,"pith_summary":"On a triangulated surface of any topology, the divergence-free part of the high-order BDM finite-element space splits exactly into the rotated gradient of a continuous streamfunction plus a finite-dimensional space of discrete harmonic fields whose dimension equals the first Betti number of the surface. Because of that splitting, every incompressible flow problem that lives in the BDM space can be rewritten with only the streamfunction and a handful of harmonic coefficients as unknowns. The resulting formulation stays exactly tangential, pointwise divergence-free and pressure-robust, yet never solves a velocity-pressure saddle-point system. The paper supplies a randomized algorithm that builds an orthonormal harmonic basis in a few sparse solves, together with hybridization and Schur-complement techniques that keep the extra harmonic unknowns cheap. Numerical experiments on a thick trefoil knot and a pierced sculpture surface show that the harmonic component carries the non-contractible circulation that a pure streamfunction cannot represent.","feed_headline":"Surface flows drop pressure, keep only streamfunction and harmonics","feed_subtitle":"BDM elements on any topology rewrite incompressible flow as a scalar plus b1 coefficients","key_machinery":"Discrete Helmholtz-Hodge decomposition of the BDM complex (Theorem 3.2): the L2-orthogonal splitting of the divergence-free subspace into rotated continuous streamfunctions of one degree higher and a harmonic complement of dimension b1(M).","core_discovery":"The divergence-free BDM subspace of degree k on a triangulated surface admits the L2-orthogonal splitting J^k_BDM = rot(S^{k+1}_0) ⊕ H^k_BDM, where dim(H^k_BDM) equals the first Betti number of the surface. Consequently every incompressible flow discretized in that subspace can be reformulated with a scalar streamfunction and finitely many harmonic coefficients as the only unknowns, eliminating pressure while retaining exact tangentiality, pointwise divergence-freeness and pressure-robustness.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Surface BDM flows drop pressure for streamfunction plus harmonics","Incompressible surface flows reduce to streamfunction and b1 harmonics","Discrete Helmholtz-Hodge splits BDM subspace into streamfunction and harmonics","Pressure-free surface Navier-Stokes via streamfunction and harmonic fields","BDM on surfaces: streamfunction and Betti harmonics replace pressure"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The dimension count that proves the harmonic space has exactly the right size rests on the Euler-Poincaré formula together with the precise degrees-of-freedom counts of the mapped polynomial spaces; if the geometric mapping or boundary conditions change those counts, the topological exactness fails.","fun_headline_variants_meta":{"raw":{"variants":["Surface BDM flows drop pressure for streamfunction plus harmonics","Incompressible surface flows reduce to streamfunction and b1 harmonics","Discrete Helmholtz-Hodge splits BDM subspace into streamfunction and harmonics","Pressure-free surface Navier-Stokes via streamfunction and harmonic fields","BDM on surfaces: streamfunction and Betti harmonics replace pressure"]},"model":"grok-4.5","effort":"low","cost_usd":0.005068,"raw_usage":{"total_tokens":1424,"prompt_tokens":773,"num_sources_used":0,"completion_tokens":73,"cost_in_usd_ticks":50680000,"prompt_tokens_details":{"text_tokens":773,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":578,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":773,"tokens_out":73,"duration_ms":4599,"temperature":1.0,"reasoning_tokens":578,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-13T16:46:47.953253+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Compute dim(H^k_BDM) by linear algebra on a closed genus-g surface for several polynomial degrees k and mesh sizes; if the dimension is not constantly equal to 2g, the central claim is false.","supporting_citations":[],"review_version":1}