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 →
Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
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^δ).
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- 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.
- 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)
- 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.
- 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.
- 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.
- 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.
- 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
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
axioms (5)
- domain assumption Beavers–Joseph–Saffman tangential slip condition with coefficient α at the interface Σ
- domain assumption Stokes domain Ω_S is simply connected so that u = curl χ with χ unique up to constants under the given boundary conditions
- standard math Mesh regularity (M1)–(M2): star-shaped polygons and edge-length lower bound with constant η > 0
- 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
- standard math Zero-mean constraint on Darcy pressure (space H¹_⋆(Ω_D)) to restore uniqueness
invented entities (2)
-
Enhanced local C¹ virtual space Ξ_h(K) with DoFs (vertex values + scaled gradients)
no independent evidence
-
vem::vemMesh2dDofFilter class
no independent evidence
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.
Reference graph
Works this paper leans on
-
[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)
2024
-
[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)
2019
-
[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)
2014
-
[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)
2016
-
[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
Pith/arXiv arXiv 2026
-
[6]
Badea, L., Discacciati, M., Quarteroni, A.: Numerical analysis of the Navier–Stokes/Darcy coupling. Numer. Math. 115(2), 195–227 (2010)
2010
-
[7]
Beavers, G.S., Joseph, D.D.: Boundary conditions at a naturally permeable wall. J. Fluid Mech.30(1), 197–207 (1967)
1967
-
[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)
2013
-
[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)
2014
-
[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)
2017
-
[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)
2017
-
[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)
2005
-
[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)
2022
-
[14]
Springer (2008)
Brenner, S., Scott, R.: The mathematical theory of finite element methods. Springer (2008)
2008
-
[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)
2024
-
[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)
2023
-
[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)
2021
-
[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/
2026
-
[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)
2015
-
[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)
2017
-
[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)
2021
-
[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)
1986
-
[23]
Calcolo55(1), 5 (2018)
Chen, L., Huang, J.: Some error analysis on virtual element methods. Calcolo55(1), 5 (2018)
2018
-
[24]
Cimolin, F., Discacciati, M.: Navier–Stokes/Forchheimer models for filtration through porous media. Appl. Numer. Math.72, 205–224 (2013)
2013
-
[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)
2011
-
[26]
Dassi, F.: VEM++, a C++ library to handle and play with the virtual element method. Numer. Algor.101, 1633–1675 (2026)
2026
-
[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)
2026
-
[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)
2002
-
[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)
2022
-
[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)
2009
-
[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)
1919
-
[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)
1911
-
[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)
2012
-
[34]
Springer Berlin, Heidelberg, Berlin (1986)
Girault, V., Raviart, P.A.: Finite Element Methods for Navier–Stokes Equations. Springer Berlin, Heidelberg, Berlin (1986)
1986
-
[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)
2006
-
[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
2021
-
[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)
2010
-
[38]
Khot, R., Mora, D., Ruiz-Baier, R.: Virtual element methods for Biot–Kirchhoff poroelasticity. Math. Comp.94(353), 1101–1146 (2025)
2025
-
[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)
2021
-
[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)
2023
-
[41]
Lippke, A., Wagner, H.: Numerical solution of the Navier-Stokes equations in multiply connected domains. Comput. Fluids20(1), 19–28 (1991)
1991
-
[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)
2000
-
[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)
2019
-
[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)
2024
-
[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)
2022
-
[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)
2015
-
[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)
2018
-
[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)
2025
-
[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)
2025
-
[50]
Saffman, P.G.: On the boundary condition at the surface of a porous medium. Stud. Appl. Math.50(2), 93–101 (1971)
1971
-
[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)
2016
-
[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
2026
-
[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)
2019
-
[54]
Cell Trans.28(1), 5–17 (2019)
Wang, X.: Bioartificial organ manufacturing technologies. Cell Trans.28(1), 5–17 (2019)
2019
-
[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)
2020
-
[56]
Xie, X., Xu, J., Xue, G.: Uniformly-stable finite element methods for Darcy-Stokes-Brinkman models. J. Comput. Math. pp. 437–455 (2008)
2008
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.