{"id":"877b7e35-bedc-420b-987f-5ca903c81b93","arxiv_id":"2607.01622","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A C1–C0 lowest-order virtual element method for the Stokes–Darcy interface problem in stream-function–pressure form is well-posed, optimally convergent, and mesh-flexible.","lead":"A new virtual-element scheme solves free-flow/porous-media coupling by replacing Stokes velocity-pressure with a stream function, cutting degrees of freedom and removing inf-sup constraints. It works on general polygons and irregular interfaces, with proofs and tests on filtration and bioartificial-organ flow.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"Mesh regularity (M1)–(M2) is the load-bearing hypothesis controlling all constants and rates in the central energy-error claim of Theorem 3.","rationale":"The reader correctly isolated the mesh-regularity hypotheses as the single softest point supporting the strongest claim (unique solvability plus the energy-error bound of Theorem 3). All other ingredients—continuous and discrete well-posedness via the coercive testing pair (χ,−φ), consistency of the projected forms a_h and c_h, and the exact interface coupling b—are standard and free of hidden gaps. The numerical tables confirm the predicted rates on meshes that satisfy (M1)–(M2), and the application experiments are consistent with the theory. No deeper inconsistency or missing estimate was found that would justify lowering the verdict from ACCEPT.","tokens_in":25061,"tokens_out":511,"duration_ms":35201,"concrete_test":"Generate a family of Voronoi meshes for the unit-square geometry of Experiment 1 in which a controlled fraction of interface-adjacent cells are deliberately made non-star-shaped (or contain edges with h_e/h_K=10^{-3}); recompute the relative energy errors of Table 1; if the observed order on the finest levels drops below 0.8, the O(h^δ) claim of Theorem 3 is not robust outside (M1)–(M2).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The energy-error argument of Lemma 8 repeatedly invokes discrete trace inequalities (||∇ξ_h·n||_{0,e}≲h_K^{-1/2}|ξ_h|_{1,K}) and the local trace inequality of [11, Lemma 6.4] to absorb interface consistency terms into the H^{2}(Ω_S) and H^{1}(Ω_D) norms. Those inequalities, together with the interpolation estimates of Lemma 7 and the polynomial projection bounds of Lemma 4, hold only when every polygon is star-shaped with respect to a ball of radius ≥η h_K and every edge satisfies h_e≥η h_K. If either condition fails, the hidden constants in the bound of Theorem 3 become uncontrolled, even though the discrete scheme (18) itself remains well-defined and uniquely solvable by the discrete inf-sup of Lemma 6.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper develops a continuous stream-function–pressure formulation of the Stokes–Darcy interface problem and a lowest-order C^{1}–C^{0} virtual element discretisation on general polygonal meshes. Incompressibility is built into the stream function, so the Stokes pressure is eliminated; mass conservation, normal-stress balance and the Beavers–Joseph–Saffman condition couple the biharmonic stream-function equation to the Darcy pressure equation. Continuous well-posedness is obtained via a kernel characterisation and a global inf-sup condition (Theorem 1). The discrete scheme (18)–(19) is shown to be uniquely solvable (Theorem 2) and to satisfy the energy-error bound of Theorem 3 under standard mesh-regularity assumptions (M1)–(M2). Numerical tests on manufactured solutions, a quarter-annulus dead-end filter and a bioartificial-organ flow network confirm the predicted rates and the geometric flexibility of the method.","tokens_in":25256,"tokens_out":1009,"duration_ms":9364,"significance":"The work fills a genuine gap: stream-function VEMs exist for pure Stokes/Navier–Stokes, and velocity–pressure VEMs exist for Stokes–Darcy, but a primal-primal stream-function–pressure VEM for the coupled interface problem has not previously been analysed. The formulation is completely free of discrete inf-sup conditions, uses only 3N_S + N_D + 4N_Σ degrees of freedom, and inherits VEM’s ability to treat non-matching polygonal interfaces without remeshing. The a-priori analysis is standard but carefully executed, and the open-source extension of vem++ (vem::vemMesh2dDofFilter) makes the scheme immediately usable. The two application-oriented experiments demonstrate practical relevance beyond manufactured solutions.","major_comments":[{"comment":"Theorem 3 and the energy-error argument of Lemma 8 rest on the global mesh-regularity hypotheses (M1)–(M2) of Section 3 (star-shapedness with radius ≥η h_K and edge-length lower bound h_e ≥ η h_K). These enter every projection, interpolation and discrete-trace estimate used to absorb the interface consistency terms. The manuscript never states that the constants in Theorem 3 become uncontrolled if small edges or non-star-shaped cells appear, even though the discrete scheme itself remains well-defined. A short remark after Theorem 3 (or in the conclusions) acknowledging this limitation would make the claim precise.","section":null},{"comment":"In the continuous analysis the Beavers–Joseph–Saffman term a_Σ is absorbed into the coercivity of A only under the standing assumption α>0 (see (10b) and Remark 2.2). When α=0 the interface contribution vanishes and the argument for coercivity on H^{2}_⋆(Ω_S) needs a different Poincaré-type control. The paper never discusses the limiting case α=0, which is occasionally of interest. A one-sentence clarification of the range of α would remove the ambiguity.","section":null}],"minor_comments":[{"comment":"Section 3.2: the enhancement condition that makes Π^∇^{2},K_2 and Π^{0},K_2 coincide is stated without a reference; a pointer to [47] or [4] would help the reader.","section":null},{"comment":"Figure 3 caption and surrounding text: the schematic of vem::vemMesh2dDofFilter is useful, but the colour coding of the two DoF sets is not explained in the caption itself.","section":null},{"comment":"Tables 1–2: the column headers use over-bars (ē_h, etc.) that are never defined; a short sentence after the definition of e_h would clarify that the tabulated quantities are the relative energy errors.","section":null},{"comment":"Experiment 4: the conversion κ=Kμ and α=γ√κ/μ is given, but the numerical values of γ and μ are stated only in the text; listing them once in a parameter table would improve reproducibility.","section":null},{"comment":"A few typographical slips: “a a priori” (p. 2), “Poincar´e” with inconsistent accents, and the arXiv identifier in the header still carries the temporary “Noname manuscript No.”","section":null}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid, carefully written contribution that fits the journal’s scope. The two major points I raise are genuinely minor clarifications rather than structural flaws; once they are addressed the paper can be accepted. The open-source filter class is a welcome practical addition."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is the first continuous analysis and lowest-order C^{1}–C^{0} VEM that couples a biharmonic stream function on the Stokes side to Darcy pressure with the full Beavers–Joseph–Saffman interface conditions. That is the real novelty: prior VEM work either did stream-function Stokes alone or velocity-pressure Stokes–Darcy. They get a uniquely solvable scheme with only 3N_S + N_D + 4N_Σ DoFs, prove the energy-error bound of Theorem 3, and ship a clean DoF-filter class in vem++ that makes the interface implementation straightforward.\n\nThe math is standard but carefully done. Continuous well-posedness follows from kernel characterization plus the global inf-sup obtained by testing with (χ,−φ). Discrete stability and the a-priori estimate reuse the usual VEM projections, diagonal/DOFI-DOFI stabilizations, and interpolation estimates. Manufactured-solution tables on quads, non-convex, Voronoi and perturbed Voronoi meshes all show the expected first-order rates; the two application runs (dead-end filter, bioartificial-organ network) recover the expected physics for a range of permeabilities. No free parameters, no circular fitting.\n\nThe soft spots are the ones the authors themselves flag. Global mesh regularity (M1)–(M2) controls every constant and rate that appears in Lemma 8 and Theorem 3; if small edges or non-star-shaped cells appear, the analysis no longer guarantees the rates even though the discrete scheme remains well-defined. The simply-connected assumption is also load-bearing for the stream-function representation; they sketch the multiply-connected and immersed-domain cases but leave them open. Those are ordinary limitations for this literature, not fatal ones. Code is referenced via the open library rather than shipped with commit-level artifacts, which is a minor reproducibility gap.\n\nThis paper is for people who already work on polytopal methods for interface flows or who need a cheap, inf-sup-free Stokes–Darcy solver on irregular meshes. It deserves a serious referee. I would accept it for peer review and would cite the formulation and the DoF count when I next need a stream-function interface scheme.","headline":"Solid first stream-function–pressure VEM for Stokes–Darcy: complete analysis, low DoF count, and working experiments; mesh regularity is the usual load-bearing hypothesis, not a hidden crack.","tokens_in":25946,"tokens_out":551,"would_cite":true,"duration_ms":5565,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N30","65N99","65Z05"],"pacs":[],"model":"grok-4.5","headline":"A stream-function virtual element method solves Stokes–Darcy coupling without velocity–pressure inf-sup pairs and with few degrees of freedom on general polygons.","keywords":["stream function formulation","Stokes–Darcy interface","virtual element method","C1-conforming VEM","Beavers–Joseph–Saffman condition","polygonal meshes","a priori error analysis"],"falsifier":"Compute the energy error on a sequence of polygonal meshes that deliberately violate the star-shapedness or edge-length condition (for example by inserting arbitrarily small edges or non-star-shaped cells) and check whether the observed convergence rate still matches the predicted O(h^δ).","tokens_in":25922,"feed_emoji":"📐","tokens_out":746,"duration_ms":6231,"temperature":0.7,"pith_summary":"The paper builds a continuous and discrete stream-function–pressure formulation of the classical Stokes–Darcy interface problem. In the free-flow region the velocity is recovered as the curl of a scalar stream function, so incompressibility holds identically and the Stokes pressure disappears; the resulting biharmonic equation is coupled to a Darcy pressure equation through mass conservation, normal-stress balance and the Beavers–Joseph–Saffman slip condition. The discrete scheme is a lowest-order C^{1}–C^{0} virtual-element method on arbitrary polygonal meshes, needing only three degrees of freedom per interior Stokes vertex, one per interior Darcy vertex and four per interface vertex. The analysis proves unique solvability without any discrete inf-sup condition and supplies an optimal energy-error bound of order h^δ. Numerical tests on manufactured solutions, a dead-end filter and a bioartificial-organ blood-flow network confirm the rates and show that irregular interfaces can be treated without remeshing.","feed_headline":"Stream-function VEM couples Stokes and Darcy on polygons","feed_subtitle":"Incompressibility is automatic, degrees of freedom drop, and irregular interfaces need no remeshing.","key_machinery":"The global discrete bilinear form A_h((χ_h,φ_h),(ξ_h,ψ_h)) = a_h(χ_h,ξ_h) + b(ξ_h,φ_h) + b(χ_h,ψ_h) – c_h(φ_h,ψ_h), whose coercivity on the product space of C^{1} stream-function and C^{0} pressure virtual elements yields a discrete inf-sup condition free of any velocity–pressure pair.","core_discovery":"The lowest-order C^{1}–C^{0} virtual-element discretisation of the stream-function–pressure Stokes–Darcy system is well-posed and converges optimally in the energy norm on general polygonal meshes that meet a uniform star-shapedness and edge-length condition, while automatically satisfying incompressibility and using far fewer degrees of freedom than a classical velocity–pressure scheme.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Stream-function VEM couples Stokes-Darcy on polygons with auto-incompressibility","C1-C0 virtual elements discretise Stokes-Darcy optimally on star-shaped meshes","Polygonal VEM stream function meets Darcy pressure at irregular interfaces","Lowest-order VEM for Stokes-Darcy reduces DOFs while preserving mass balance","Biharmonic stream function VEM enforces Beavers-Joseph-Saffman on polygons"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"Every mesh polygon must be star-shaped with respect to a ball of radius at least a fixed fraction of its diameter, and every edge must be at least that same fraction long; if those shape-regularity constants fail, the projection, interpolation and trace estimates that underwrite the error bound no longer hold.","fun_headline_variants_meta":{"raw":{"variants":["Stream-function VEM couples Stokes-Darcy on polygons with auto-incompressibility","C1-C0 virtual elements discretise Stokes-Darcy optimally on star-shaped meshes","Polygonal VEM stream function meets Darcy pressure at irregular interfaces","Lowest-order VEM for Stokes-Darcy reduces DOFs while preserving mass balance","Biharmonic stream function VEM enforces Beavers-Joseph-Saffman on polygons"]},"model":"grok-4.5","effort":"low","cost_usd":0.004462,"raw_usage":{"total_tokens":1259,"prompt_tokens":677,"num_sources_used":0,"completion_tokens":111,"cost_in_usd_ticks":44620000,"prompt_tokens_details":{"text_tokens":677,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":471,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":677,"tokens_out":111,"duration_ms":4910,"temperature":1.0,"reasoning_tokens":471,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-12T08:40:45.906372+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Compute the energy error on a sequence of polygonal meshes that deliberately violate the star-shapedness or edge-length condition (for example by inserting arbitrarily small edges or non-star-shaped cells) and check whether the observed convergence rate still matches the predicted O(h^δ).","supporting_citations":[],"review_version":2}