Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

Andres E. Rubiano; David Mora; Franco Dassi; Rekha Khot; Ricardo Ruiz-Baier

Reviewed by Pith at T0; open to challenge.

T0 means a machine referee read the full paper against a public rubric. The mark states how deep the mechanical check went, never who wrote it. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.3

A stream function virtual element method solves the Stokes-Darcy interface problem on general polygonal meshes while enforcing mass conservation and stress balance.

2026-07-03 08:59 UTC pith:IPX7LVKV

load-bearing objection This VEM paper couples a stream-function Stokes discretization to Darcy pressure across interfaces on general polygons, but offers no analysis and thin numerical evidence. the 2 major comments →

arxiv 2607.01622 v2 pith:IPX7LVKV submitted 2026-07-02 math.NA cs.NA

Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

Franco Dassi , Rekha Khot , David Mora , Andres E. Rubiano , Ricardo Ruiz-Baier This is my paper

classification math.NA cs.NA

keywords virtual element methodStokes-Darcy interfacestream function formulationBeavers-Joseph-Saffman conditionpolygonal meshesmass conservationprimal-primal formulation

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a virtual element method for the coupled Stokes-Darcy system in primal-primal form. In the Stokes domain a stream function formulation is used that automatically satisfies the incompressibility constraint. The biharmonic stream function equation is coupled to the Darcy pressure equation by enforcing mass conservation, normal stress balance, and the Beavers-Joseph-Saffman condition at the interface. The virtual element approach accommodates irregular interface geometries on polygonal meshes without remeshing. Numerical simulations on dead-end filtration and bioartificial organ flows confirm the method's accuracy.

Core claim

The proposed virtual element method for the Stokes-Darcy interface problem in primal-primal form with stream function formulation satisfies mass conservation, normal stress balance, and Beavers-Joseph-Saffman conditions and is accurate on general polygonal meshes as shown by numerical simulations.

What carries the argument

Stream function virtual element discretization of the biharmonic equation in the Stokes domain coupled to a pressure equation in the Darcy domain, with interface conditions imposed directly.

Load-bearing premise

Virtual element spaces can discretize the biharmonic stream function equation while stably coupling to the Darcy pressure equation across irregular interfaces without added stabilization.

What would settle it

A simulation on a polygonal mesh with a non-smooth interface that shows loss of mass conservation or failure to match normal stress would disprove the central claim.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

Mass is conserved across the interface by construction.
Normal stress balance and the Beavers-Joseph-Saffman slip condition hold at the discrete level.
The method works on arbitrary polygonal meshes without requiring interface-aligned refinement.
Computational cost is reduced because the stream function formulation removes the velocity divergence constraint.
The same framework applies to dead-end filtration and network flow in bioartificial organs.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same stream-function VEM coupling could be tested on time-dependent or nonlinear interface problems.
Convergence rates on successively refined polygonal meshes would provide quantitative error bounds.
Extension to three-dimensional polyhedral meshes would require only replacement of the two-dimensional virtual spaces.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

This VEM paper couples a stream-function Stokes discretization to Darcy pressure across interfaces on general polygons, but offers no analysis and thin numerical evidence.

read the letter

The core contribution is a primal-primal VEM that replaces the usual velocity-pressure Stokes system with a stream-function formulation, then weakly imposes the three interface conditions against the Darcy pressure space. The method is presented as working directly on arbitrary polygonal meshes without remeshing or extra interface fitting.

What stands out is the choice to stay in primal variables and exploit VEM's flexibility on non-matching interface edges. The applications to dead-end filtration and bioartificial organ networks are reasonable test cases for the target community.

The main weaknesses are the absence of any a priori error estimates or stability proof for the coupled system, and the lack of reported convergence rates, mesh details, or quantitative error tables. The stress-test concern about discrete inf-sup or coercivity for the biharmonic VEM traces is not obviously resolved by the abstract; if the paper simply transplants the continuous interface integrals without additional stabilization or mortar terms, stability on distorted interface polygons remains an open question. Numerical validation is asserted but not documented enough to check mass conservation or condition-number behavior.

This is for specialists in VEM or Stokes-Darcy numerics who want to see one more discretization option. A reader needing rigorous guarantees or reproducible benchmarks will not get much. It is worth sending to referees because the formulation is concrete and the application examples are relevant, even though the theoretical and experimental gaps are large enough that heavy revision would be expected.

Referee Report

2 major / 1 minor

Summary. The paper proposes a novel virtual element discretization for the Stokes-Darcy interface problem in primal-primal form. A stream-function formulation is used for the biharmonic Stokes problem to enforce incompressibility exactly, while the Darcy region is treated via a pressure equation; the three interface conditions (normal-velocity continuity, normal-stress balance, and Beavers-Joseph-Saffman slip) are imposed weakly. The method is asserted to handle arbitrary polygonal meshes without remeshing or stabilization and is validated by numerical simulations on applications including dead-end filtration and bioartificial-organ network flow.

Significance. If the discrete spaces satisfy the necessary stability and convergence properties on general meshes, the approach would offer a practical tool for coupled free/porous-media flows on complex geometries, reducing degrees of freedom via the stream-function reduction and exploiting VEM's polygonal flexibility. No machine-checked proofs or parameter-free derivations are claimed.

major comments (2)

[Abstract and §3] Abstract and §3 (discrete formulation): the central claim that the scheme remains stable and mass-conserving on arbitrary polygonal interfaces rests on the assertion that the standard VEM bilinear form for the biharmonic stream-function problem, augmented only by the continuous interface integrals, satisfies a discrete inf-sup condition with the Darcy pressure space. No proof or estimate is supplied showing that the virtual degrees of freedom and projection operators control the required traces sufficiently to prevent loss of coercivity or mass conservation when interface edges are highly distorted or non-aligned.
[Numerical results] Numerical results section: the abstract states that accuracy is validated by simulations, yet no convergence rates, error tables, mesh-refinement studies, or quantitative comparison against analytic solutions or other methods are referenced. Without these data the claim that the method is accurate on general meshes cannot be assessed and is load-bearing for the paper's main assertion.

minor comments (1)

[§2-3] Notation for the interface conditions and the precise definition of the discrete spaces (including the polynomial degree and projection operators) should be stated explicitly in the formulation section to allow reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important aspects regarding theoretical justification and quantitative validation. We address each major comment below and indicate the revisions we will make to strengthen the manuscript.

read point-by-point responses

Referee: [Abstract and §3] Abstract and §3 (discrete formulation): the central claim that the scheme remains stable and mass-conserving on arbitrary polygonal interfaces rests on the assertion that the standard VEM bilinear form for the biharmonic stream-function problem, augmented only by the continuous interface integrals, satisfies a discrete inf-sup condition with the Darcy pressure space. No proof or estimate is supplied showing that the virtual degrees of freedom and projection operators control the required traces sufficiently to prevent loss of coercivity or mass conservation when interface edges are highly distorted or non-aligned.

Authors: We acknowledge that the manuscript does not contain a self-contained proof or a priori estimate establishing the discrete inf-sup condition for the coupled system on arbitrarily distorted polygonal interfaces. The formulation relies on the known stability of the individual VEM spaces for the biharmonic stream-function problem and the Darcy pressure equation, together with the weak enforcement of the interface conditions. A full analysis of the coupled inf-sup constant under general interface distortions would require additional technical work on the trace and projection operators. In the revision we will add a dedicated remark in §3 clarifying the stability assumptions inherited from the literature and include a brief numerical check of the inf-sup constant on a sequence of increasingly distorted interface meshes. A complete rigorous proof is not feasible within the scope of the present revision. revision: partial
Referee: [Numerical results] Numerical results section: the abstract states that accuracy is validated by simulations, yet no convergence rates, error tables, mesh-refinement studies, or quantitative comparison against analytic solutions or other methods are referenced. Without these data the claim that the method is accurate on general meshes cannot be assessed and is load-bearing for the paper's main assertion.

Authors: The referee is correct that the current numerical section presents only qualitative results from application examples without tabulated errors or observed convergence rates. We will revise the numerical results section to include a new subsection containing a manufactured-solution test with known analytic solution. This subsection will report L2 and H1 error tables, computed convergence rates under successive mesh refinement, and results on both quasi-uniform and highly distorted polygonal meshes. These additions will provide the quantitative evidence needed to support the accuracy claims. revision: yes

Circularity Check

0 steps flagged

No circularity: direct VEM discretization with independent numerical validation

full rationale

The paper introduces a primal-primal VEM discretization for the Stokes-Darcy system using a stream-function formulation in the free-flow region, with interface conditions (mass conservation, normal stress balance, Beavers-Joseph-Saffman) imposed directly in the variational form. No equations reduce to fitted parameters renamed as predictions, no self-definitional loops, and no load-bearing self-citations that substitute for independent justification. The central claims rest on the construction of the discrete spaces and bilinear forms plus numerical experiments on polygonal meshes; these steps are self-contained against external benchmarks and do not collapse to the inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Review based on abstract only; no explicit free parameters, invented entities, or ad-hoc axioms are identifiable from the given text. The work relies on established VEM theory and standard physical interface conditions.

axioms (2)

standard math Virtual element methods can be applied to biharmonic problems on polygonal meshes
Invoked implicitly by using VEM for the stream function equation.
domain assumption Interface conditions (mass conservation, normal stress balance, Beavers-Joseph-Saffman slip) are the appropriate coupling for Stokes-Darcy
These conditions are stated as enforced across the interface.

pith-pipeline@v0.9.1-grok · 5667 in / 1263 out tokens · 24312 ms · 2026-07-03T08:59:24.918295+00:00 · methodology

0 comments

read the original abstract

This paper introduces a novel Virtual Element Method (VEM) for the coupled Stokes--Darcy system in primal-primal form. In the free-flow Stokes domain, we implement a stream function formulation that inherently satisfies the incompressibility constraint and reduces computational cost. Across the interface, mass conservation, normal stress balance, and the Beavers--Joseph--Saffman slip condition are enforced to couple the biharmonic stream function equation with the Darcy's pressure equation. Leveraging VEM's ability to handle general polygonal meshes, the proposed method naturally accommodates irregular interface geometries without requiring remeshing or adaptive refinement. The accuracy of the method is validated through several numerical simulations that include applications to dead-end filtration, and network flow in bioartificial organs.

Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

Core claim

What carries the argument

Load-bearing premise

What would settle it

If this is right

Where Pith is reading between the lines

discussion (0)