Pith. sign in

REVIEW 2 major objections 5 minor 56 references

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.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-12 08:40 UTC pith:IPX7LVKV

load-bearing objection 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. 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

classification math.NA cs.NA MSC 65N3065N9965Z05
keywords stream function formulationStokes–Darcy interfacevirtual element methodC1-conforming VEMBeavers–Joseph–Saffman conditionpolygonal meshesa priori error analysis
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 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.

Core claim

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.

What carries the argument

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.

Load-bearing premise

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.

What would settle it

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^δ).

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

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

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.

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 (2)
  1. 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.
  2. 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.
minor comments (5)
  1. 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.
  2. 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.
  3. 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.
  4. 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.
  5. 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.”

Circularity Check

0 steps flagged

No circularity: standard a-priori VEM analysis of a new stream-function–pressure Stokes–Darcy scheme; self-citations supply background spaces only.

full rationale

The derivation chain is self-contained and non-circular. Continuous well-posedness (Theorem 1) follows from operator bounds (Lemma 1), kernel characterisation (Lemma 2), and a global inf-sup obtained by the elementary test (ξ,ψ)=(χ,−φ) (Lemma 3); none of these steps is definitional or imported as a uniqueness theorem. The discrete scheme (18)–(19) inherits the same structure with standard VEM consistency/stability of the projected bilinear forms plus positive-semidefinite stabilisations (17); unique solvability is proved afresh via the discrete inf-sup (Lemma 6 / Theorem 2). The energy-error estimate (Lemma 8) is a classical consistency-plus-interpolation argument that absorbs interface terms by discrete trace/inverse inequalities under the explicit mesh hypotheses (M1)–(M2); rates (Theorem 3) then follow from the polynomial and interpolation estimates of Lemmas 4 and 7. Manufactured solutions used in the numerical section are independent of the scheme and no parameters are fitted to the error tables. Self-citations ([1,3,45,48] for uncoupled stream-function VEM spaces, [8,9,10] for basic VEM projections/stabilisations, [43] for a related but different Stokes–Darcy VEM) supply background constructions that are not load-bearing for the coupled analysis or the claimed rates. Mesh regularity is an assumption controlling constants, not a circular reduction. Consequently the central claims do not reduce to their inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 5 axioms · 2 invented entities

The central claim rests on classical continuum interface conditions, standard VEM mesh and projection theory, and the stream-function representation in simply-connected 2D domains. No free parameters are fitted to data; physical coefficients are inputs. Invented entities are limited to the discrete spaces and the mesh-filter class, which are constructive definitions rather than postulated physical objects.

axioms (5)
  • domain assumption Beavers–Joseph–Saffman tangential slip condition with coefficient α at the interface Σ
    Enters the continuous weak form (5c) and the bilinear form a_Σ; taken from classical interface modeling [7,50].
  • domain assumption Stokes domain Ω_S is simply connected so that u = curl χ with χ unique up to constants under the given boundary conditions
    Used to introduce the stream-function space H²_ΓS (Section 2.2); multi-connected case is deferred to Remark 2.2.
  • standard math Mesh regularity (M1)–(M2): star-shaped polygons and edge-length lower bound with constant η > 0
    Standard VEM assumption (Section 3); required for polynomial approximation and inverse/trace estimates.
  • standard math Stabilization forms S^K_S and S^K_D are positive semi-definite and spectrally equivalent to the consistency forms on the kernels of the projections
    Assumed in (17a)–(17b); any such stabilizations work; concrete recipes used only in numerics.
  • standard math Zero-mean constraint on Darcy pressure (space H¹_⋆(Ω_D)) to restore uniqueness
    Follows from kernel characterization Lemma 2; imposed by Lagrange multiplier in the code.
invented entities (2)
  • Enhanced local C¹ virtual space Ξ_h(K) with DoFs (vertex values + scaled gradients) no independent evidence
    purpose: Conformingly approximate the stream function on polygons while making H² and L² projections computable from DoFs
    Constructive discrete space built from standard VEM enhancement; independent_evidence false because it is a numerical construct, not a physical object.
  • vem::vemMesh2dDofFilter class no independent evidence
    purpose: Assign different VEM spaces to subdomains and keep only interface-shared DoFs
    Software device enabling the mixed C¹–C⁰ discretization; no external physical prediction.

pith-pipeline@v1.1.0-grok45 · 29064 in / 2929 out tokens · 30113 ms · 2026-07-12T08:40:45.906372+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.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

56 extracted references · 1 linked inside Pith

  1. [1]

    Adak, D., Mora, D., Silgado, A.: The Morley-type virtual element method for the Navier–Stokes equations in stream- function form. Comput. Methods Appl. Mech. Engrg.419, 116573 (2024)

  2. [2]

    Anaya, V., Mora, D., Reales, C., Ruiz-Baier, R.: Vorticity-pressure formulations for the Brinkman-Darcy coupled problem. Numer. Methods Partial Differ. Equ.35(2), 528–544 (2019)

  3. [3]

    Antonietti, P.F., Beir˜ ao da Veiga, L., Mora, D., Verani, M.: A stream virtual element formulation of the Stokes problem on polygonal meshes. SIAM J. Numer. Anal.52(1), 386–404 (2014)

  4. [4]

    Antonietti, P.F., Beir˜ ao da Veiga, L., Scacchi, S., Verani, M.: A C 1 virtual element method for the Cahn–Hilliard equation with polygonal meshes. SIAM J. Numer. Anal.54(1), 34–56 (2016)

  5. [5]

    arXiv preprint2606.11935(2026)

    Antonietti, P.F., Botti, M., Parolini, N., Pederzoli, V., Verani, M.: Polytopal discontinuous Galerkin discretizations of coupled non-Newtonian Stokes-Darcy systems. arXiv preprint2606.11935(2026). Available at:https://arxiv. org/abs/2606.11935

  6. [6]

    Badea, L., Discacciati, M., Quarteroni, A.: Numerical analysis of the Navier–Stokes/Darcy coupling. Numer. Math. 115(2), 195–227 (2010)

  7. [7]

    Beavers, G.S., Joseph, D.D.: Boundary conditions at a naturally permeable wall. J. Fluid Mech.30(1), 197–207 (1967)

  8. [8]

    Beir˜ ao da Veiga, L., Brezzi, F., Cangiani, A., Manzini, G., Marini, L.D., Russo, A.: Basic principles of virtual element methods. Math. Models Methods Appl. Sci.23(01), 199–214 (2013)

  9. [9]

    Beir˜ ao da Veiga, L., Brezzi, F., Marini, L.D., Russo, A.: The Hitchhiker’s guide to the virtual element method. Math. Models Methods Appl. Sci.24(08), 1541–1573 (2014)

  10. [10]

    Beir˜ ao da Veiga, L., Dassi, F., Russo, A.: High-order virtual element method on polyhedral meshes. Comput. Math. Appl.74(5), 1110–1122 (2017)

  11. [11]

    Beir˜ ao da Veiga, L., Lovadina, C., Russo, A.: Stability analysis for the virtual element method. Math. Models Methods Appl. Sci.27(13), 2557–2594 (2017)

  12. [12]

    ESAIM: Math

    Bernardi, C., Hecht, F., Pironneau, O.: Coupling Darcy and Stokes equations for porous media with cracks. ESAIM: Math. Model. Numer. Anal.39(1), 7–35 (2005)

  13. [13]

    Boon, W.M., Koch, T., Kuchta, M., Mardal, K.A.: Robust monolithic solvers for the Stokes–Darcy problem with the Darcy equation in primal form. SIAM J. Sci. Comput.44(4), B1148–B1174 (2022)

  14. [14]

    Springer (2008)

    Brenner, S., Scott, R.: The mathematical theory of finite element methods. Springer (2008)

  15. [15]

    PLOS Comput

    Bukaˇ c, M.,ˇCani´ c, S., Muha, B., Wang, Y.: A computational algorithm for optimal design of a bioartificial organ scaffold architecture. PLOS Comput. Biol.20(11) (2024)

  16. [16]

    ESAIM: Math

    Bukaˇ c, M., Muha, B., Salgado, A.J.: Analysis of a diffuse interface method for the Stokes–Darcy coupled problem. ESAIM: Math. Model. Numer. Anal.57(5), 2623–2658 (2023)

  17. [17]

    Burman, E., Delay, G., Ern, A.: An unfitted hybrid high-order method for the Stokes interface problem. IMA J. Numer. Anal.41(4), 2362–2387 (2021)

  18. [18]

    Preprint 2026-13, Centro de Investigaci´ on en Ingenier´ ıa Matem´ atica (CI²MA), Universidad de Concepci´ on (2026)

    Cama˜ no, J., Oyarz´ ua, R., Rojo, K., Villa-Fuentes, S.: A mixed finite element method based on pseudostress and stream-function for the Navier–Stokes problem in 2D. Preprint 2026-13, Centro de Investigaci´ on en Ingenier´ ıa Matem´ atica (CI²MA), Universidad de Concepci´ on (2026). Available at:http://www.ci2ma.udec.cl/publications/ preprints/

  19. [19]

    Cama˜ no, J., Gatica, G.N., Oyarz´ ua, R., Ruiz-Baier, R., Venegas, P.: New fully-mixed finite element methods for the Stokes–Darcy coupling. Comput. Methods Appl. Mech. Eng.295, 362–395 (2015)

  20. [20]

    Cangiani, A., Georgoulis, E.H., Pryer, T., Sutton, O.J.: A posteriori error estimates for the virtual element method. Numer. Math.137, 857–893 (2017)

  21. [21]

    Carstensen, C., Mallik, G., Nataraj, N.: Nonconforming finite element discretization for semilinear problems with trilinear nonlinearity. IMA J. Numer. Anal.41(1), 164–205 (2021)

  22. [22]

    Cayco, M.E., Nicolaides, R.A.: Finite element technique for optimal pressure recovery from stream function formu- lation of viscous flows. Math. Comp.46(174), 371–377 (1986)

  23. [23]

    Calcolo55(1), 5 (2018)

    Chen, L., Huang, J.: Some error analysis on virtual element methods. Calcolo55(1), 5 (2018)

  24. [24]

    Cimolin, F., Discacciati, M.: Navier–Stokes/Forchheimer models for filtration through porous media. Appl. Numer. Math.72, 205–224 (2013)

  25. [25]

    ESAIM: Math

    D’Angelo, C., Zunino, P.: Robust numerical approximation of coupled Stokes’ and Darcy’s flows applied to vascular hemodynamics and biochemical transport. ESAIM: Math. Model. Numer. Anal.45(3), 447–476 (2011)

  26. [26]

    Dassi, F.: VEM++, a C++ library to handle and play with the virtual element method. Numer. Algor.101, 1633–1675 (2026)

  27. [27]

    Dassi, F., Khot, R., Rubiano, A.E., Ruiz-Baier, R.: Analysis and virtual element discretisation of a Stokes/Biot– Kirchhoff bulk–surface model. Comput. Methods Appl. Mech. Eng.449, 118545 (2026)

  28. [28]

    Discacciati, M., Miglio, E., Quarteroni, A.: Mathematical and numerical models for coupling surface and groundwater flows. Appl. Numer. Math.43(1-2), 57–74 (2002)

  29. [29]

    ESAIM: Math

    Eggenweiler, E., Discacciati, M., Rybak, I.: Analysis of the Stokes–Darcy problem with generalised interface condi- tions. ESAIM: Math. Model. Numer. Anal.56(2), 727–742 (2022)

  30. [30]

    Ervin, V., Jenkins, E.W., Sun, S.: Coupled generalized nonlinear Stokes flow with flow through a porous medium. SIAM J. Numer. Anal.47(2), 929–952 (2009)

  31. [31]

    SIAM Journal on Scientific Computing24(6), 1919–1929 (2003)

    Fairag, F.: Numerical computations of viscous, incompressible flow problems using a two-level finite element method. SIAM Journal on Scientific Computing24(6), 1919–1929 (2003)

  32. [32]

    Gatica, G., Oyarz´ ua, R., Sayas, F.J.: Analysis of fully-mixed finite element methods for the Stokes–Darcy coupled problem. Math. Comp.80, 1911–1948 (2011)

  33. [33]

    Gatica, G.N., Oyarz´ ua, R., Sayas, F.J.: A twofold saddle point approach for the coupling of fluid flow with nonlinear porous media flow. IMA J. Numer. Anal.32(3), 845–887 (2012)

  34. [34]

    Springer Berlin, Heidelberg, Berlin (1986)

    Girault, V., Raviart, P.A.: Finite Element Methods for Navier–Stokes Equations. Springer Berlin, Heidelberg, Berlin (1986)

  35. [35]

    Hanspal, N.S., Waghode, A.N., Nassehi, V., Wakeman, R.J.: Numerical analysis of coupled Stokes/Darcy flows in industrial filtrations. Trans. Porous Media64(1), 73–101 (2006)

  36. [36]

    Huang, J., Yu, Y.: A medius error analysis for nonconforming virtual element methods for Poisson and biharmonic equations. J. Comput. Appl. Math.386, 113229 (2021) Stream function – pressure virtual element methods for the Stokes–Darcy interface problem 23

  37. [37]

    Kanschat, G., Rivi` ere, B.: A strongly conservative finite element method for the coupling of Stokes and Darcy flow. J. Comput. Phys.229(17), 5933–5943 (2010)

  38. [38]

    Khot, R., Mora, D., Ruiz-Baier, R.: Virtual element methods for Biot–Kirchhoff poroelasticity. Math. Comp.94(353), 1101–1146 (2025)

  39. [39]

    Com- put

    Kim, D., Pani, A.K., Park, E.J.: Morley finite element methods for the stationary quasi-geostrophic equation. Com- put. Methods Appl. Mech. Engrg.375, 113639 (2021)

  40. [40]

    Gels9(2) (2023)

    Krishani, M., Shin, W.Y., Suhaimi, H., Sambudi, N.S.: Development of scaffolds from bio-based natural materials for tissue regeneration applications: A review. Gels9(2) (2023)

  41. [41]

    Lippke, A., Wagner, H.: Numerical solution of the Navier-Stokes equations in multiply connected domains. Comput. Fluids20(1), 19–28 (1991)

  42. [42]

    Liu, J.G., Shu, C.W.: A high-order discontinuous Galerkin method for 2D incompressible flows. J. Comput. Phys. 160(2), 577–596 (2000)

  43. [43]

    Calcolo56(4), 48 (2019)

    Liu, X., Li, R., Chen, Z.: A virtual element method for the coupled Stokes–Darcy problem with the Beaver–Joseph– Saffman interface condition. Calcolo56(4), 48 (2019)

  44. [44]

    Mishra, S., Natarajan, E.: A unified local projection-based stabilized virtual element method for the coupled Stokes– Darcy problem. Adv. Comput. Math.50(6), 106 (2024)

  45. [45]

    Mora, D., Reales, C., Silgado, A.: AC 1-virtual element method of high order for the Brinkman equations in stream function formulation with pressure recovery. IMA J. Numer. Anal.42, 3632–3674 (2022)

  46. [46]

    Mora, D., Rivera, G., Rodr´ ıguez, R.: A virtual element method for the Steklov eigenvalue problem. Math. Models Methods Appl. Sci.25(08), 1421–1445 (2015)

  47. [47]

    ESAIM: Math

    Mora, D., Rivera, G., Vel´ asquez, I.: A virtual element method for the vibration problem of Kirchhoff plates. ESAIM: Math. Model. Numer. Anal.52(4), 1437–1456 (2018)

  48. [48]

    Mora, D., Silgado, A.: Stream virtual elements for the Navier-Stokes system: Nonstandard boundary conditions and variable recovery algorithms. SIAM J. Sci. Comput.47, 207–237 (2025)

  49. [49]

    Qi, F., Liu, Y., Chen, R., Hu, P., Cai, X.C.: Computational modeling of the kidney hemodynamics with a coupled unsteady Stokes–Darcy model. Int. J. Numer. Methods Biomed. Eng.41(11), e70108 (2025)

  50. [50]

    Saffman, P.G.: On the boundary condition at the surface of a porous medium. Stud. Appl. Math.50(2), 93–101 (1971)

  51. [51]

    Song, S., Faleo, G., Yeung, R., Kant, R., Posselt, A.M., Desai, T.A., Tang, Q., Roy, S.: Silicon Nanopore Membrane (SNM) for islet encapsulation and immunoisolation under convective transport. Sci. Rep.6, 23679 (2016)

  52. [52]

    ESAIM: Math

    Vexler, B., Leykekhman, D., Wagner, J.: Fully discrete error analysis of finite element discretizations of time- dependent Stokes equations in a stream-function formulation. ESAIM: Math. Model. Numer. Anal. (2026). In press

  53. [53]

    Wang, G., Wang, F., Chen, L., He, Y.: A divergence free weak virtual element method for the Stokes–Darcy problem on general meshes. Comput. Methods Appl. Mech. Eng.344, 998–1020 (2019)

  54. [54]

    Cell Trans.28(1), 5–17 (2019)

    Wang, X.: Bioartificial organ manufacturing technologies. Cell Trans.28(1), 5–17 (2019)

  55. [55]

    Wen, J., Su, J., He, Y., Chen, H.: A discontinuous Galerkin method for the coupled Stokes and Darcy problem. J. Sci. Comput.85(2), 26 (2020)

  56. [56]

    Xie, X., Xu, J., Xue, G.: Uniformly-stable finite element methods for Darcy-Stokes-Brinkman models. J. Comput. Math. pp. 437–455 (2008)