Superconvergent and Divergence-Free Mixed Finite Element Methods for The Stokes Equation

Chao Zhang; Long Chen; Xinyue Zhao; Xuehai Huang

arxiv: 2510.14192 · v3 · submitted 2025-10-16 · 🧮 math.NA · cs.NA

Superconvergent and Divergence-Free Mixed Finite Element Methods for The Stokes Equation

Long Chen , Xuehai Huang , Chao Zhang , Xinyue Zhao This is my paper

Pith reviewed 2026-05-18 06:58 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords mixed finite element methodsStokes equationsdivergence-free methodssuperconvergenceinf-sup conditionH(div) spacesvector Laplacian

0 comments

The pith

A mixed finite element method for the Stokes equation transfers the distributional divergence-free property from discrete to continuous velocity spaces to achieve superconvergence and stability without stabilization.

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

This paper develops a divergence-free mixed finite element method for the Stokes equation by pairing H(div)-conforming velocities with discontinuous pressures to ensure pointwise divergence-free velocities. It introduces tangential-normal continuous traceless tensor elements to discretize the stress and proves an inf-sup condition for the weak divergence operator between these spaces. The setup provides stable discretization of the vector Laplacian without extra stabilization terms. The key innovation is the property that a stress field distributionally divergence-free against the discrete divergence-free velocity space remains so against the continuous space, decoupling errors. This yields optimal error estimates for the stress and convergence rates for velocity and pressure that exceed the approximation orders of the spaces.

Core claim

The authors establish that for their mixed finite element scheme, a stress field that is distributionally divergence-free against the discrete divergence-free velocity space is also distributionally divergence-free against the continuous divergence-free velocity space. This transfer property, together with the inf-sup stability of the weak div operator, decouples the stress and velocity errors and produces optimal-order error estimates for the stress while the velocity and pressure converge at rates higher than the approximation orders of the chosen spaces.

What carries the argument

The tangential-normal continuous traceless tensor stress space together with the weak divergence operator whose inf-sup condition is proved and the distributional divergence-free transfer property between discrete and continuous spaces.

If this is right

Optimal-order error estimates are obtained for the stress.
Velocity and pressure converge at rates higher than the approximation orders of the chosen spaces.
The discretization of the vector Laplacian remains stable without additional stabilization terms.
Velocities are exactly pointwise divergence-free.

Where Pith is reading between the lines

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

The error decoupling could extend the approach to related incompressible flow models where similar distributional orthogonality arguments apply.
Absence of stabilization parameters may simplify code development and parameter tuning for practical simulations of Stokes problems.
The superconvergence property suggests potential use in a posteriori error estimators that exploit the transfer between discrete and continuous spaces.

Load-bearing premise

The inf-sup condition for the weak div operator between the tangential-normal continuous traceless tensor stress space and the H(div)-conforming velocity space holds along with the distributional divergence-free transfer property.

What would settle it

A direct numerical check computing the distributional action of the discrete stress divergence on a continuous divergence-free velocity test function outside the discrete space; nonzero values when the discrete pairing is zero would falsify the transfer property and the resulting superconvergence.

read the original abstract

This paper develops divergence-free mixed finite element methods for the Stokes equation. Using H(div)-conforming velocities and discontinuous pressures ensures the inf-sup condition for the velocity--pressure pair and yields pointwise divergence-free velocities. However, this choice makes the vector Laplacian difficult to discretize. Inspired by mass-conserving mixed formulations with stresses, tangential--normal continuous traceless tensor elements are introduced to discretize the vector Laplacian. An inf-sup condition for the weak div operator between the stress and velocity spaces is then proved. Two key properties characterize the scheme. First, the stress--velocity inf-sup stability gives a stable discretization of the vector Laplacian without additional stabilization, unlike discontinuous Galerkin or virtual element methods. Second, the scheme has the property that if a stress field is distributionally divergence-free against the discrete divergence-free velocity space, then it is also distributionally divergence-free against the continuous divergence-free velocity space. This property decouples the stress and velocity errors and leads to superconvergence. As a result, optimal-order error estimates are obtained for the stress, while the velocity and pressure converge at rates higher than the approximation orders of the chosen spaces. Numerical experiments confirm the theoretical results.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a mixed FEM for Stokes with exact div-free velocities and superconvergence via new traceless tensor elements for the stress.

read the letter

The main takeaway is a mixed finite element scheme for the Stokes equations that keeps velocities pointwise divergence-free and delivers superconvergent rates for velocity and pressure. It uses H(div)-conforming velocities paired with tangential-normal continuous traceless tensor elements to handle the vector Laplacian in a mixed stress formulation. The key is proving an inf-sup condition for the weak div operator between these spaces, plus a distributional divergence-free transfer property that maps discrete kernel membership to the continuous case. This decouples the errors and pushes convergence orders higher than the basic approximation degrees of the spaces. Optimal rates hold for the stress, with the velocity and pressure doing better than expected. Numerical tests are reported to match the theory. The construction avoids extra stabilization terms that show up in some DG or VEM approaches, which is a practical plus. The analysis rests on the new element spaces and the two proved properties rather than on fitted parameters or self-referential definitions. That keeps the circularity burden low. Soft spots are modest. Everything turns on the specific choice of traceless tensor elements and the transfer property holding in general. The abstract does not spell out mesh restrictions or implementation details, so one would want to check those in the full proofs for any hidden limitations on degree or regularity. Still, the claims line up without obvious internal contradictions. This work is aimed at numerical analysts working on stable discretizations for incompressible flows. Readers who need exact mass conservation and higher accuracy without heavy stabilization will find it relevant. It has enough new analysis and verification to merit a serious referee.

Referee Report

0 major / 3 minor

Summary. The manuscript develops a divergence-free mixed finite element scheme for the Stokes equations. It employs H(div)-conforming velocity spaces and discontinuous pressures to enforce pointwise divergence-free velocities, while introducing tangential-normal continuous traceless tensor elements to discretize the stress and handle the vector Laplacian. An inf-sup condition is proved for the weak divergence operator between the new stress space and the velocity space. The scheme is shown to satisfy a distributional divergence-free transfer property: if a stress is distributionally divergence-free against the discrete kernel, it is also against the continuous kernel. This property decouples the stress and velocity errors, yielding stability without stabilization and superconvergent rates (optimal order for stress, higher than the approximation order for velocity and pressure). Numerical experiments are presented to confirm the theoretical estimates.

Significance. If the inf-sup proof and the distributional transfer property hold as stated, the work offers a notable contribution to structure-preserving discretizations for incompressible flows. The avoidance of stabilization terms while maintaining stability and achieving superconvergence distinguishes the approach from many DG and VEM alternatives. The error decoupling mechanism could prove useful in related mixed formulations for other saddle-point problems.

minor comments (3)

The introduction would benefit from a short paragraph explicitly comparing the new tangential-normal continuous traceless tensor space to existing stress-based mixed methods for Stokes or elasticity, to better situate the novelty.
In the numerical section, the convergence plots would be clearer if the observed rates were labeled directly on the figures rather than only stated in the caption.
A brief remark on the dependence of the inf-sup constant on the polynomial degree would help readers assess the practical range of applicability.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading and positive assessment of our manuscript. The recommendation for minor revision is appreciated, and we will incorporate improvements to enhance clarity and presentation in the revised version.

Circularity Check

0 steps flagged

No significant circularity; derivation relies on proved properties for new spaces

full rationale

The paper constructs new tangential-normal continuous traceless tensor stress spaces and H(div)-conforming velocity spaces, then proves an inf-sup condition for the weak divergence operator and a distributional divergence-free transfer property mapping discrete kernel membership to continuous kernel membership. These proofs directly enable stability without stabilization and error decoupling for superconvergence. No load-bearing step reduces by definition, fitted parameter, or self-citation chain to the target result itself; the claimed optimal stress estimates and superconvergent velocity/pressure rates follow from the established properties rather than circular redefinition.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The method introduces one new finite element space and relies on standard Sobolev-space assumptions plus two proved stability properties; no free parameters or invented physical entities appear.

axioms (2)

domain assumption The inf-sup condition holds for the weak divergence operator between the tangential-normal continuous traceless tensor stress space and the H(div)-conforming velocity space.
Invoked to guarantee stability of the vector Laplacian discretization without additional terms.
domain assumption If a stress is distributionally divergence-free against the discrete divergence-free velocity space then it is distributionally divergence-free against the continuous divergence-free velocity space.
Central to decoupling stress and velocity errors and obtaining superconvergence.

invented entities (1)

tangential-normal continuous traceless tensor elements no independent evidence
purpose: Discretize the vector Laplacian while preserving the required continuity and traceless properties for the mixed formulation.
Newly introduced finite element space whose definition and approximation properties are part of the contribution.

pith-pipeline@v0.9.0 · 5746 in / 1678 out tokens · 39331 ms · 2026-05-18T06:58:06.393394+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

57 extracted references · 57 canonical work pages

[1]

D. N. Arnold.Finite element exterior calculus, volume 93 ofCBMS-NSF Regional Conference Series in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2018. 14

work page 2018
[2]

D. N. Arnold and F. Brezzi. Mixed and nonconforming finite element methods: Implementation, postpor- cessing and error estimates.RAIRO Model Math. Anal. Numer., 19:7–32, 1985. 5 21 TABLE4. Errors∥u−u ∗ h∥and∥grad h(u−u ∗ h)∥for Example 5.2 with ℓ= 0andk= 0,1in three dimensions. h (k, ℓ) ∥u−u ∗ h∥order ∥gradh(u−u ∗ h)∥order 2−1 (0,0) 4.167e-04 – 4.576e-03 ...

work page 1985
[3]

D. N. Arnold, F. Brezzi, and M. Fortin. A stable finite element for the Stokes equations.Calcolo, 21(4):337– 344, 1984. 1

work page 1984
[4]

D. N. Arnold, R. S. Falk, and J. Gopalakrishnan. Mixed finite element approximation of the vector Laplacian with Dirichlet boundary conditions.Math. Models Methods Appl. Sci., 22(9):1250024, 26, 2012. 2

work page 2012
[5]

D. N. Arnold, R. S. Falk, and R. Winther. Finite element exterior calculus, homological techniques, and applications.Acta Numer., 15:1–155, 2006. 9, 14

work page 2006
[6]

D. N. Arnold and J. Qin. Quadratic velocity/linear pressure Stokes elements. InAdvances in Computer Methods for Partial Differential Equations-VII, IMACS, pages 28—34, 1992. 1

work page 1992
[7]

Ayuso de Dios, K

B. Ayuso de Dios, K. Lipnikov, and G. Manzini. The nonconforming virtual element method.ESAIM Math. Model. Numer. Anal., 50(3):879–904, 2016. 16

work page 2016
[8]

G. A. Baker, W. N. Jureidini, and O. A. Karakashian. Piecewise solenoidal vector fields and the Stokes problem.SIAM J. Numer. Anal., 27(6):1466–1485, 1990. 1

work page 1990
[9]

Beir ˜ao da Veiga, C

L. Beir ˜ao da Veiga, C. Lovadina, and G. Vacca. Divergence free virtual elements for the Stokes problem on polygonal meshes.ESAIM Math. Model. Numer. Anal., 51(2):509–535, 2017. 1

work page 2017
[10]

Boffi, F

D. Boffi, F. Brezzi, and M. Fortin.Mixed finite element methods and applications, volume 44 ofSpringer Series in Computational Mathematics. Springer, Heidelberg, 2013. 1

work page 2013
[11]

Braess and J

D. Braess and J. Sch ¨oberl. Equilibrated residual error estimator for edge elements.Math. Comp., 77(262):651–672, 2008. 4

work page 2008
[12]

Brezzi, J

F. Brezzi, J. Douglas, Jr., and L. D. Marini. Two families of mixed finite elements for second order elliptic problems.Numer. Math., 47(2):217–235, 1985. 1, 6

work page 1985
[13]

C. Chen, X. Huang, and H. Wei. Virtual element methods without extrinsic stabilization.SIAM J. Numer. Anal., 62(1):567–591, 2024. 14, 17

work page 2024
[14]

Chen.ifem: an integrated finite element methods package in MATLAB

L. Chen.ifem: an integrated finite element methods package in MATLAB. Technical report, University of California at Irvine, 2009. 19

work page 2009
[15]

Chen and X

L. Chen and X. Huang. Nonconforming virtual element method for2mth order partial differential equations inR n.Math. Comp., 89(324):1711–1744, 2020. 16, 17

work page 2020
[16]

Chen and X

L. Chen and X. Huang. Finite elements for div- and divdiv-conforming symmetric tensors in arbitrary di- mension.SIAM J. Numer. Anal., 60(4):1932–1961, 2022. 6

work page 1932
[17]

Chen and X

L. Chen and X. Huang. Finite element de Rham and Stokes complexes in three dimensions.Math. Comp., 93(345):55–110, 2024. 1

work page 2024
[18]

Chen and X

L. Chen and X. Huang. Finite element complexes in two dimensions (in Chinese).Sci. Sin. Math., 55(8):1593–1626, 2025. 1

work page 2025
[19]

L. Chen, X. Huang, and C. Zhang. Distributional finite element curl div complexes and application to quad curl problems.SIAM J. Numer. Anal., 63(3):1078–1104, 2025. 6

work page 2025
[20]

Chen and F

L. Chen and F. Wang. A divergence free weak virtual element method for the Stokes problem on polytopal meshes.J. Sci. Comput., 78(2):864–886, 2019. 1, 2

work page 2019
[21]

L. Chen, M. Wang, and L. Zhong. Convergence analysis of triangular MAC schemes for two dimensional Stokes equations.J. Sci. Comput., 63(3):716–744, 2015. 2, 3 22

work page 2015
[22]

S. H. Christiansen and K. Hu. Generalized finite element systems for smooth differential forms and Stokes’ problem.Numer. Math., 140(2):327–371, 2018. 1

work page 2018
[23]

Cockburn and F.-J

B. Cockburn and F.-J. Sayas. Divergence-conforming HDG methods for Stokes flows.Math. Comp., 83(288):1571–1598, 2014. 2

work page 2014
[24]

Costabel and A

M. Costabel and A. McIntosh. On Bogovski ˘ı and regularized Poincar´e integral operators for de Rham com- plexes on Lipschitz domains.Math. Z., 265(2):297–320, 2010. 14

work page 2010
[25]

Crouzeix and P.-A

M. Crouzeix and P.-A. Raviart. Conforming and nonconforming finite element methods for solving the stationary Stokes equations. I.Rev. Franc ¸aise Automat. Informat. Recherche Op ´erationnelle S ´er. Rouge, 7:33–75, 1973. 1, 16

work page 1973
[26]

Dubois, M

F. Dubois, M. Sala ¨un, and S. Salmon. First vorticity-velocity-pressure numerical scheme for the Stokes problem.Comput. Methods Appl. Mech. Engrg., 192(44-46):4877–4907, 2003. 2

work page 2003
[27]

Dubois, M

F. Dubois, M. Sala ¨un, and S. Salmon. V orticity-velocity-pressure and stream function-vorticity formulations for the Stokes problem.J. Math. Pures Appl. (9), 82(11):1395–1451, 2003. 2

work page 2003
[28]

G. Fu, J. Guzm ´an, and M. Neilan. Exact smooth piecewise polynomial sequences on Alfeld splits.Math. Comp., 89(323):1059–1091, 2020. 1

work page 2020
[29]

Girault and P.-A

V . Girault and P.-A. Raviart.Finite element methods for Navier-Stokes equations, volume 5 ofSpringer Series in Computational Mathematics. Springer-Verlag, Berlin, 1986. Theory and algorithms. 13

work page 1986
[30]

Gopalakrishnan, P

J. Gopalakrishnan, P. L. Lederer, and J. Sch ¨oberl. A mass conserving mixed stress formulation for Stokes flow with weakly imposed stress symmetry.SIAM J. Numer. Anal., 58(1):706–732, 2020. 2, 3, 6, 19

work page 2020
[31]

Gopalakrishnan, P

J. Gopalakrishnan, P. L. Lederer, and J. Sch¨oberl. A mass conserving mixed stress formulation for the Stokes equations.IMA J. Numer. Anal., 40(3):1838–1874, 2020. 2, 3, 18, 19

work page 2020
[32]

Guzm ´an, A

J. Guzm ´an, A. Lischke, and M. Neilan. Exact sequences on Powell-Sabin splits.Calcolo, 57(2):Paper No. 13, 25, 2020. 1

work page 2020
[33]

Guzm ´an and M

J. Guzm ´an and M. Neilan. Conforming and divergence-free Stokes elements in three dimensions.IMA J. Numer. Anal., 34(4):1489–1508, 2014. 1

work page 2014
[34]

Guzm ´an and M

J. Guzm ´an and M. Neilan. Conforming and divergence-free Stokes elements on general triangular meshes. Math. Comp., 83(285):15–36, 2014. 1

work page 2014
[35]

F. H. Harlow, J. E. Welch, et al. Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface.Physics of fluids, 8(12):2182, 1965. 3

work page 1965
[36]

J. Hu, Y . Huang, and Q. Lin. Lower bounds for eigenvalues of elliptic operators: by nonconforming finite element methods.J. Sci. Comput., 61(1):196–221, 2014. 16

work page 2014
[37]

Hu and R

J. Hu and R. Ma. The enriched Crouzeix-Raviart elements are equivalent to the Raviart-Thomas elements. J. Sci. Comput., 63(2):410–425, 2015. 16

work page 2015
[38]

K. Hu, Q. Zhang, and Z. Zhang. A family of finite element Stokes complexes in three dimensions.SIAM J. Numer. Anal., 60(1):222–243, 2022. 1

work page 2022
[39]

Huang and Z

X. Huang and Z. Tang. Robust and optimal mixed methods for a fourth-order elliptic singular perturbation problem.arXiv preprint arXiv:2501.12137, 2025. 17

work page arXiv 2025
[40]

Huang and F

X. Huang and F. Wang. Analysis of divergence free conforming virtual elements for the Brinkman problem. Math. Models Methods Appl. Sci., 33(6):1245–1280, 2023. 1

work page 2023
[41]

V . John, A. Linke, C. Merdon, M. Neilan, and L. G. Rebholz. On the divergence constraint in mixed finite element methods for incompressible flows.SIAM Rev., 59(3):492–544, 2017. 1

work page 2017
[42]

A. Linke. On the role of the Helmholtz decomposition in mixed methods for incompressible flows and a new variational crime.Comput. Methods Appl. Mech. Engrg., 268:782–800, 2014. 18

work page 2014
[43]

K. A. Mardal, X.-C. Tai, and R. Winther. A robust finite element method for Darcy-Stokes flow.SIAM J. Numer. Anal., 40(5):1605–1631, 2002. 1

work page 2002
[44]

Maz’ya and J

V . Maz’ya and J. Rossmann.Elliptic equations in polyhedral domains, volume 162 ofMathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 2010. 13

work page 2010
[45]

N ´ed´elec

J.-C. N ´ed´elec. Mixed finite elements inR 3.Numer. Math., 35(3):315–341, 1980. 1, 6

work page 1980
[46]

N ´ed´elec

J.-C. N ´ed´elec. A new family of mixed finite elements inR 3.Numer. Math., 50(1):57–81, 1986. 1, 6

work page 1986
[47]

Raviart and J

P.-A. Raviart and J. M. Thomas. A mixed finite element method for 2nd order elliptic problems. InMathe- matical aspects of finite element methods (Proc. Conf., Consiglio Naz. delle Ricerche (C.N.R.), Rome, 1975), volume V ol. 606 ofLecture Notes in Math., pages 292–315. Springer, Berlin-New York, 1977. 1, 6

work page 1975
[48]

L. R. Scott and M. V ogelius. Norm estimates for a maximal right inverse of the divergence operator in spaces of piecewise polynomials.RAIRO Mod ´el. Math. Anal. Num´er., 19(1):111–143, 1985. 1

work page 1985
[49]

Tai and R

X.-C. Tai and R. Winther. A discrete de Rham complex with enhanced smoothness.Calcolo, 43(4):287–306,

work page
[50]

Taylor and P

C. Taylor and P. Hood. A numerical solution of the Navier-Stokes equations using the finite element tech- nique.Computers. & Fluids, 1(1):73–100, 1973. 1 23

work page 1973
[51]

Wang and S

C. Wang and S. Zhang. Auto-stabilized weak Galerkin finite element methods for Stokes equations on non- convex polytopal meshes.J. Comput. Phys., 533:Paper No. 114006, 22, 2025. 1

work page 2025
[52]

Wang and X

J. Wang and X. Ye. New finite element methods in computational fluid dynamics byH(div)elements. SIAM J. Numer. Anal., 45(3):1269–1286, 2007. 1, 2

work page 2007
[53]

H. Wei, X. Huang, and A. Li. Piecewise divergence-free nonconforming virtual elements for Stokes problem in any dimensions.SIAM J. Numer. Anal., 59(3):1835–1856, 2021. 1, 17

work page 2021
[54]

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

work page 2008
[55]

S. Zhang. A new family of stable mixed finite elements for the 3D Stokes equations.Math. Comp., 74(250):543–554, 2005. 1

work page 2005
[56]

S. Zhang. Divergence-free finite elements on tetrahedral grids fork≥6.Math. Comp., 80(274):669–695,

work page
[57]

1 DEPARTMENT OFMATHEMATICS, UNIVERSITY OFCALIFORNIA ATIRVINE, IRVINE, CA 92697, USA Email address:chenlong@math.uci.edu SCHOOL OFMATHEMATICS, SHANGHAIUNIVERSITY OFFINANCE ANDECONOMICS, SHANGHAI200433, CHINA Email address:huang.xuehai@sufe.edu.cn SCHOOL OFGENERALEDUCATION, WENZHOUBUSINESSCOLLEGE, WENZHOU325035, CHINA Email address:zhang.chao@wzbc.edu.cn SC...

work page

[1] [1]

D. N. Arnold.Finite element exterior calculus, volume 93 ofCBMS-NSF Regional Conference Series in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2018. 14

work page 2018

[2] [2]

D. N. Arnold and F. Brezzi. Mixed and nonconforming finite element methods: Implementation, postpor- cessing and error estimates.RAIRO Model Math. Anal. Numer., 19:7–32, 1985. 5 21 TABLE4. Errors∥u−u ∗ h∥and∥grad h(u−u ∗ h)∥for Example 5.2 with ℓ= 0andk= 0,1in three dimensions. h (k, ℓ) ∥u−u ∗ h∥order ∥gradh(u−u ∗ h)∥order 2−1 (0,0) 4.167e-04 – 4.576e-03 ...

work page 1985

[3] [3]

D. N. Arnold, F. Brezzi, and M. Fortin. A stable finite element for the Stokes equations.Calcolo, 21(4):337– 344, 1984. 1

work page 1984

[4] [4]

D. N. Arnold, R. S. Falk, and J. Gopalakrishnan. Mixed finite element approximation of the vector Laplacian with Dirichlet boundary conditions.Math. Models Methods Appl. Sci., 22(9):1250024, 26, 2012. 2

work page 2012

[5] [5]

D. N. Arnold, R. S. Falk, and R. Winther. Finite element exterior calculus, homological techniques, and applications.Acta Numer., 15:1–155, 2006. 9, 14

work page 2006

[6] [6]

D. N. Arnold and J. Qin. Quadratic velocity/linear pressure Stokes elements. InAdvances in Computer Methods for Partial Differential Equations-VII, IMACS, pages 28—34, 1992. 1

work page 1992

[7] [7]

Ayuso de Dios, K

B. Ayuso de Dios, K. Lipnikov, and G. Manzini. The nonconforming virtual element method.ESAIM Math. Model. Numer. Anal., 50(3):879–904, 2016. 16

work page 2016

[8] [8]

G. A. Baker, W. N. Jureidini, and O. A. Karakashian. Piecewise solenoidal vector fields and the Stokes problem.SIAM J. Numer. Anal., 27(6):1466–1485, 1990. 1

work page 1990

[9] [9]

Beir ˜ao da Veiga, C

L. Beir ˜ao da Veiga, C. Lovadina, and G. Vacca. Divergence free virtual elements for the Stokes problem on polygonal meshes.ESAIM Math. Model. Numer. Anal., 51(2):509–535, 2017. 1

work page 2017

[10] [10]

Boffi, F

D. Boffi, F. Brezzi, and M. Fortin.Mixed finite element methods and applications, volume 44 ofSpringer Series in Computational Mathematics. Springer, Heidelberg, 2013. 1

work page 2013

[11] [11]

Braess and J

D. Braess and J. Sch ¨oberl. Equilibrated residual error estimator for edge elements.Math. Comp., 77(262):651–672, 2008. 4

work page 2008

[12] [12]

Brezzi, J

F. Brezzi, J. Douglas, Jr., and L. D. Marini. Two families of mixed finite elements for second order elliptic problems.Numer. Math., 47(2):217–235, 1985. 1, 6

work page 1985

[13] [13]

C. Chen, X. Huang, and H. Wei. Virtual element methods without extrinsic stabilization.SIAM J. Numer. Anal., 62(1):567–591, 2024. 14, 17

work page 2024

[14] [14]

Chen.ifem: an integrated finite element methods package in MATLAB

L. Chen.ifem: an integrated finite element methods package in MATLAB. Technical report, University of California at Irvine, 2009. 19

work page 2009

[15] [15]

Chen and X

L. Chen and X. Huang. Nonconforming virtual element method for2mth order partial differential equations inR n.Math. Comp., 89(324):1711–1744, 2020. 16, 17

work page 2020

[16] [16]

Chen and X

L. Chen and X. Huang. Finite elements for div- and divdiv-conforming symmetric tensors in arbitrary di- mension.SIAM J. Numer. Anal., 60(4):1932–1961, 2022. 6

work page 1932

[17] [17]

Chen and X

L. Chen and X. Huang. Finite element de Rham and Stokes complexes in three dimensions.Math. Comp., 93(345):55–110, 2024. 1

work page 2024

[18] [18]

Chen and X

L. Chen and X. Huang. Finite element complexes in two dimensions (in Chinese).Sci. Sin. Math., 55(8):1593–1626, 2025. 1

work page 2025

[19] [19]

L. Chen, X. Huang, and C. Zhang. Distributional finite element curl div complexes and application to quad curl problems.SIAM J. Numer. Anal., 63(3):1078–1104, 2025. 6

work page 2025

[20] [20]

Chen and F

L. Chen and F. Wang. A divergence free weak virtual element method for the Stokes problem on polytopal meshes.J. Sci. Comput., 78(2):864–886, 2019. 1, 2

work page 2019

[21] [21]

L. Chen, M. Wang, and L. Zhong. Convergence analysis of triangular MAC schemes for two dimensional Stokes equations.J. Sci. Comput., 63(3):716–744, 2015. 2, 3 22

work page 2015

[22] [22]

S. H. Christiansen and K. Hu. Generalized finite element systems for smooth differential forms and Stokes’ problem.Numer. Math., 140(2):327–371, 2018. 1

work page 2018

[23] [23]

Cockburn and F.-J

B. Cockburn and F.-J. Sayas. Divergence-conforming HDG methods for Stokes flows.Math. Comp., 83(288):1571–1598, 2014. 2

work page 2014

[24] [24]

Costabel and A

M. Costabel and A. McIntosh. On Bogovski ˘ı and regularized Poincar´e integral operators for de Rham com- plexes on Lipschitz domains.Math. Z., 265(2):297–320, 2010. 14

work page 2010

[25] [25]

Crouzeix and P.-A

M. Crouzeix and P.-A. Raviart. Conforming and nonconforming finite element methods for solving the stationary Stokes equations. I.Rev. Franc ¸aise Automat. Informat. Recherche Op ´erationnelle S ´er. Rouge, 7:33–75, 1973. 1, 16

work page 1973

[26] [26]

Dubois, M

F. Dubois, M. Sala ¨un, and S. Salmon. First vorticity-velocity-pressure numerical scheme for the Stokes problem.Comput. Methods Appl. Mech. Engrg., 192(44-46):4877–4907, 2003. 2

work page 2003

[27] [27]

Dubois, M

F. Dubois, M. Sala ¨un, and S. Salmon. V orticity-velocity-pressure and stream function-vorticity formulations for the Stokes problem.J. Math. Pures Appl. (9), 82(11):1395–1451, 2003. 2

work page 2003

[28] [28]

G. Fu, J. Guzm ´an, and M. Neilan. Exact smooth piecewise polynomial sequences on Alfeld splits.Math. Comp., 89(323):1059–1091, 2020. 1

work page 2020

[29] [29]

Girault and P.-A

V . Girault and P.-A. Raviart.Finite element methods for Navier-Stokes equations, volume 5 ofSpringer Series in Computational Mathematics. Springer-Verlag, Berlin, 1986. Theory and algorithms. 13

work page 1986

[30] [30]

Gopalakrishnan, P

J. Gopalakrishnan, P. L. Lederer, and J. Sch ¨oberl. A mass conserving mixed stress formulation for Stokes flow with weakly imposed stress symmetry.SIAM J. Numer. Anal., 58(1):706–732, 2020. 2, 3, 6, 19

work page 2020

[31] [31]

Gopalakrishnan, P

J. Gopalakrishnan, P. L. Lederer, and J. Sch¨oberl. A mass conserving mixed stress formulation for the Stokes equations.IMA J. Numer. Anal., 40(3):1838–1874, 2020. 2, 3, 18, 19

work page 2020

[32] [32]

Guzm ´an, A

J. Guzm ´an, A. Lischke, and M. Neilan. Exact sequences on Powell-Sabin splits.Calcolo, 57(2):Paper No. 13, 25, 2020. 1

work page 2020

[33] [33]

Guzm ´an and M

J. Guzm ´an and M. Neilan. Conforming and divergence-free Stokes elements in three dimensions.IMA J. Numer. Anal., 34(4):1489–1508, 2014. 1

work page 2014

[34] [34]

Guzm ´an and M

J. Guzm ´an and M. Neilan. Conforming and divergence-free Stokes elements on general triangular meshes. Math. Comp., 83(285):15–36, 2014. 1

work page 2014

[35] [35]

F. H. Harlow, J. E. Welch, et al. Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface.Physics of fluids, 8(12):2182, 1965. 3

work page 1965

[36] [36]

J. Hu, Y . Huang, and Q. Lin. Lower bounds for eigenvalues of elliptic operators: by nonconforming finite element methods.J. Sci. Comput., 61(1):196–221, 2014. 16

work page 2014

[37] [37]

Hu and R

J. Hu and R. Ma. The enriched Crouzeix-Raviart elements are equivalent to the Raviart-Thomas elements. J. Sci. Comput., 63(2):410–425, 2015. 16

work page 2015

[38] [38]

K. Hu, Q. Zhang, and Z. Zhang. A family of finite element Stokes complexes in three dimensions.SIAM J. Numer. Anal., 60(1):222–243, 2022. 1

work page 2022

[39] [39]

Huang and Z

X. Huang and Z. Tang. Robust and optimal mixed methods for a fourth-order elliptic singular perturbation problem.arXiv preprint arXiv:2501.12137, 2025. 17

work page arXiv 2025

[40] [40]

Huang and F

X. Huang and F. Wang. Analysis of divergence free conforming virtual elements for the Brinkman problem. Math. Models Methods Appl. Sci., 33(6):1245–1280, 2023. 1

work page 2023

[41] [41]

V . John, A. Linke, C. Merdon, M. Neilan, and L. G. Rebholz. On the divergence constraint in mixed finite element methods for incompressible flows.SIAM Rev., 59(3):492–544, 2017. 1

work page 2017

[42] [42]

A. Linke. On the role of the Helmholtz decomposition in mixed methods for incompressible flows and a new variational crime.Comput. Methods Appl. Mech. Engrg., 268:782–800, 2014. 18

work page 2014

[43] [43]

K. A. Mardal, X.-C. Tai, and R. Winther. A robust finite element method for Darcy-Stokes flow.SIAM J. Numer. Anal., 40(5):1605–1631, 2002. 1

work page 2002

[44] [44]

Maz’ya and J

V . Maz’ya and J. Rossmann.Elliptic equations in polyhedral domains, volume 162 ofMathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 2010. 13

work page 2010

[45] [45]

N ´ed´elec

J.-C. N ´ed´elec. Mixed finite elements inR 3.Numer. Math., 35(3):315–341, 1980. 1, 6

work page 1980

[46] [46]

N ´ed´elec

J.-C. N ´ed´elec. A new family of mixed finite elements inR 3.Numer. Math., 50(1):57–81, 1986. 1, 6

work page 1986

[47] [47]

Raviart and J

P.-A. Raviart and J. M. Thomas. A mixed finite element method for 2nd order elliptic problems. InMathe- matical aspects of finite element methods (Proc. Conf., Consiglio Naz. delle Ricerche (C.N.R.), Rome, 1975), volume V ol. 606 ofLecture Notes in Math., pages 292–315. Springer, Berlin-New York, 1977. 1, 6

work page 1975

[48] [48]

L. R. Scott and M. V ogelius. Norm estimates for a maximal right inverse of the divergence operator in spaces of piecewise polynomials.RAIRO Mod ´el. Math. Anal. Num´er., 19(1):111–143, 1985. 1

work page 1985

[49] [49]

Tai and R

X.-C. Tai and R. Winther. A discrete de Rham complex with enhanced smoothness.Calcolo, 43(4):287–306,

work page

[50] [50]

Taylor and P

C. Taylor and P. Hood. A numerical solution of the Navier-Stokes equations using the finite element tech- nique.Computers. & Fluids, 1(1):73–100, 1973. 1 23

work page 1973

[51] [51]

Wang and S

C. Wang and S. Zhang. Auto-stabilized weak Galerkin finite element methods for Stokes equations on non- convex polytopal meshes.J. Comput. Phys., 533:Paper No. 114006, 22, 2025. 1

work page 2025

[52] [52]

Wang and X

J. Wang and X. Ye. New finite element methods in computational fluid dynamics byH(div)elements. SIAM J. Numer. Anal., 45(3):1269–1286, 2007. 1, 2

work page 2007

[53] [53]

H. Wei, X. Huang, and A. Li. Piecewise divergence-free nonconforming virtual elements for Stokes problem in any dimensions.SIAM J. Numer. Anal., 59(3):1835–1856, 2021. 1, 17

work page 2021

[54] [54]

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

work page 2008

[55] [55]

S. Zhang. A new family of stable mixed finite elements for the 3D Stokes equations.Math. Comp., 74(250):543–554, 2005. 1

work page 2005

[56] [56]

S. Zhang. Divergence-free finite elements on tetrahedral grids fork≥6.Math. Comp., 80(274):669–695,

work page

[57] [57]

1 DEPARTMENT OFMATHEMATICS, UNIVERSITY OFCALIFORNIA ATIRVINE, IRVINE, CA 92697, USA Email address:chenlong@math.uci.edu SCHOOL OFMATHEMATICS, SHANGHAIUNIVERSITY OFFINANCE ANDECONOMICS, SHANGHAI200433, CHINA Email address:huang.xuehai@sufe.edu.cn SCHOOL OFGENERALEDUCATION, WENZHOUBUSINESSCOLLEGE, WENZHOU325035, CHINA Email address:zhang.chao@wzbc.edu.cn SC...

work page