pith. sign in

arxiv: 2607.01933 · v1 · pith:XFL3EKGVnew · submitted 2026-07-02 · 🧮 math.NA · cs.NA

Low-regularity finite element elasticity complexes with hybridizable stresses on tetrahedral Alfeld splits

Johnny Guzm\'an , Xuehai Huang This is my paper

Pith reviewed 2026-07-03 08:07 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords finite elementselasticity complexAlfeld splithybridizable stresssymmetric tensorde Rham complexBernstein-Gelfand-Gelfandlow regularity
0
0 comments X

The pith

Finite element elasticity complexes of low regularity are constructed on tetrahedral Alfeld splits, ending in hybridizable H(div;S) symmetric stress spaces with no vertex degrees of freedom.

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

The paper constructs finite element elasticity complexes on tetrahedral Alfeld splits that require lower Sobolev regularity and lower polynomial degrees than prior versions on the same meshes. These complexes terminate in a symmetric stress space that is conforming in H(div), hybridizable, and free of vertex degrees of freedom. The construction proceeds by applying local Bernstein-Gelfand-Gelfand arguments to polynomial de Rham complexes on the Alfeld split, yielding two exact local sequences whose bubble subcomplexes and dimensions are computed explicitly. If the sequences are exact, the resulting global subcomplexes supply unisolvent elements for displacement and incompatibility fields together with commuting interpolation operators. This matters for discretizing three-dimensional linear elasticity because the reduced regularity and absence of vertex degrees of freedom can simplify implementation while preserving the structure of the continuous complex.

Core claim

Two local polynomial elasticity complexes are proved on the Alfeld split: an H²-H¹(inc) complex and a lower-regularity H¹(curl)-H(inc⁺) complex. Their bubble subcomplexes and dimension formulas are derived. These local exact sequences produce unisolvent finite elements for the displacement and incompatibility spaces and global finite element subcomplexes of the corresponding elasticity sequences. In the lowest-order H¹(curl)-H(inc⁺) case the terminal stress-displacement pair recovers the Johnson-Mercier-Křížek element, while the construction supplies hybridizable symmetric stresses of all orders k≥1. A second family yields a low-regularity H¹-H(inc) complex for the standard elasticity sequen

What carries the argument

Local Bernstein-Gelfand-Gelfand arguments applied to polynomial de Rham complexes on the Alfeld split, which generate exact elasticity sequences together with their bubble subcomplexes and dimension counts.

If this is right

  • Unisolvent finite elements are obtained for the displacement and incompatibility spaces.
  • Global finite element subcomplexes of the elasticity sequences are obtained on Alfeld-split tetrahedral meshes.
  • The lowest-order terminal pair recovers the Johnson-Mercier-Křížek element.
  • Commuting interpolation diagrams hold for both global complexes.
  • Hybridizable symmetric stress elements are available for every polynomial degree k≥1.

Where Pith is reading between the lines

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

  • The hybridizable stress space without vertex degrees of freedom may simplify the assembly of hybridizable discontinuous Galerkin schemes for elasticity.
  • The same local-sequence technique could be tested on other tetrahedral splits or on polyhedral meshes to produce analogous low-regularity complexes.
  • Because the stress space is piecewise cubic at lowest order, the method may reduce the number of global degrees of freedom relative to higher-regularity complexes of comparable accuracy.

Load-bearing premise

The local Bernstein-Gelfand-Gelfand arguments applied to polynomial de Rham complexes on the Alfeld split produce exact sequences with the claimed bubble subcomplexes and dimension formulas.

What would settle it

An explicit low-degree polynomial on a single Alfeld-split tetrahedron that lies in the kernel of one operator in the claimed local sequence but is not in the range of the preceding operator, or a global interpolator that fails to commute with the differential operators at the lowest order.

read the original abstract

Finite element elasticity complexes of low regularity are constructed on tetrahedral Alfeld splits. In comparison with existing three-dimensional elasticity complexes on such splits, the complexes constructed here lower both the Sobolev regularity and the polynomial degrees, while ending in a hybridizable $H({\rm div};\mathbb S)$-conforming symmetric stress space with no vertex degrees of freedom. The construction is obtained from local Bernstein-Gelfand-Gelfand arguments applied to polynomial de Rham complexes on the Alfeld split. Two local polynomial elasticity complexes are proved: an $H^2$-$H^1({\rm inc})$ complex and a lower-regularity $H^1({\rm curl})$-$H({\rm inc}^+)$ complex. Their bubble subcomplexes and dimension formulas are derived. These local exact sequences lead to unisolvent finite elements for the displacement and incompatibility spaces and to global finite element subcomplexes of the corresponding elasticity sequences. In the lowest-order $H^1({\rm curl})$-$H({\rm inc}^+)$ finite element complex, the $H({\rm inc}^+;\mathbb S)$-conforming tensor space is piecewise cubic. At the same order, the terminal stress-displacement pair recovers the Johnson-Mercier-K\v{r}\'{i}\v{z}ek element, while the construction covers higher-order hybridizable symmetric stresses for all $k\ge1$. A second family gives a low-regularity $H^1$-$H({\rm inc})$ finite element complex for the standard elasticity sequence for all $k\ge2$. Commuting interpolation diagrams are established for both global complexes.

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.

Referee Report

0 major / 2 minor

Summary. The manuscript constructs finite element elasticity complexes of reduced Sobolev regularity and polynomial degree on tetrahedral Alfeld splits. Local Bernstein-Gelfand-Gelfand arguments applied to polynomial de Rham complexes yield two exact local polynomial sequences—an H²-H¹(inc) complex and a lower-regularity H¹(curl)-H(inc⁺) complex—together with their bubble subcomplexes and dimension formulas. These local sequences are assembled into global hybridizable subcomplexes of the elasticity sequence, with unisolvent elements for the displacement and incompatibility spaces and commuting interpolation diagrams. The lowest-order H¹(curl)-H(inc⁺) case recovers the Johnson-Mercier-Křížek element (piecewise-cubic symmetric stress), while the construction supplies hybridizable H(div; S)-conforming stresses without vertex degrees of freedom for all k ≥ 1 and a second family for the standard elasticity sequence when k ≥ 2.

Significance. If the local exactness statements and dimension counts hold, the work supplies new, lower-regularity finite-element spaces for three-dimensional elasticity that reduce both the polynomial degree and the number of degrees of freedom relative to existing Alfeld-split complexes. The hybridizable terminal stress space and the explicit recovery of the Johnson-Mercier-Křížek pair at lowest order provide concrete computational advantages and an independent consistency check. The commuting diagrams further enable stable mixed discretizations for problems with limited regularity.

minor comments (2)
  1. The precise polynomial degrees of each space in the two families (beyond the lowest-order case) should be tabulated explicitly, e.g., in a table following the dimension formulas, to facilitate implementation and comparison with other elements.
  2. Notation for the incompatibility operator inc⁺ and the precise definition of the hybridizable trace space on faces should be recalled in the global-assembly section for readers who skip the local-construction details.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading and positive evaluation of the manuscript. We are pleased that the referee recognizes the construction of the low-regularity complexes, the recovery of the Johnson-Mercier-Křížek element, and the computational advantages of the hybridizable stress spaces.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained in homological algebra

full rationale

The paper constructs the complexes by proving two local polynomial elasticity sequences via BGG arguments applied to de Rham complexes on the Alfeld split, then assembles them into global hybridizable elements. These proofs supply explicit bubble subcomplexes and dimension counts. The lowest-order case recovers the independent Johnson-Mercier-Křížek element as a consistency check. No load-bearing step reduces to a fitted parameter, self-definition, or self-citation chain; all central claims rest on direct verification of exactness within standard finite-element exterior calculus. The work is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper relies on standard properties of polynomial de Rham complexes and BGG arguments with no free parameters fitted to data, no ad-hoc axioms beyond domain assumptions in finite element theory, and no invented entities.

axioms (1)
  • domain assumption Polynomial de Rham complexes on Alfeld splits admit the required exact sequences under the stated regularity.
    Invoked to apply local BGG arguments and derive the elasticity complexes.

pith-pipeline@v0.9.1-grok · 5825 in / 1405 out tokens · 27229 ms · 2026-07-03T08:07:42.298590+00:00 · methodology

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

29 extracted references · 29 canonical work pages

  1. [1]

    Amstutz and N

    S. Amstutz and N. Van Goethem. The incompatibility operator: from Riemann’s intrinsic view of geometry to a new model of elasto-plasticity. InTopics in Applied Analysis and Optimisation, pages 33–70. Springer, 2019

    work page 2019

  2. [2]

    Arnold, G

    D. Arnold, G. Awanou, and R. Winther. Finite elements for symmetric tensors in three dimensions.Math. Comp., 77(263):1229–1251, 2008

    work page 2008

  3. [3]

    D. N. Arnold, J. Douglas, Jr., and C. P. Gupta. A family of higher order mixed finite element methods for plane elasticity.Numer. Math., 45(1):1–22, 1984

    work page 1984

  4. [4]

    D. N. Arnold, R. S. Falk, and R. Winther. Differential complexes and stability of finite element methods II: The elasticity complex. In D. Arnold, P. Bochev, R. Lehoucq, R. Nicolaides, and M. Shashkov, editors,Compatible Spatial Discretizations, volume 142 ofIMA Vol. Math. Appl., pages 47–68. Springer, Berlin, 2006

    work page 2006

  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

    work page 2006

  6. [6]

    D. N. Arnold and K. Hu. Complexes from complexes.Found. Comput. Math., 21(6):1739–1774, 2021

    work page 2021

  7. [7]

    D. N. Arnold and R. Winther. Mixed finite elements for elasticity.Numer. Math., 92(3):401–419, 2002

    work page 2002

  8. [8]

    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

    work page 1985

  9. [9]

    Chen and X

    L. Chen and X. Huang. Discrete Hessian complexes in three dimensions. InThe virtual element method and its applications, volume 31 ofSEMA SIMAI Springer Ser., pages 93–135. Springer, Cham, 2022

    work page 2022

  10. [10]

    Chen and X

    L. Chen and X. Huang. A finite element elasticity complex in three dimensions.Math. Comp., 91(337):2095–2127, 2022

    work page 2095

  11. [11]

    Chen and X

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

    work page 1932

  12. [12]

    Chen and X

    L. Chen and X. Huang. Finite elements for div div conforming symmetric tensors in three dimensions.Math. Comp., 91(335):1107–1142, 2022

    work page 2022

  13. [13]

    Chen and X

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

    work page 2025

  14. [14]

    Chen and X

    L. Chen and X. Huang. Hybridizable symmetric stress elements on the barycentric refinement in arbitrary dimensions. Math. Comp., https://arxiv.org/abs/2501.02691, 2025. LOW-REGULARITY ELASTICITY COMPLEXES 23

    work page arXiv 2025

  15. [15]

    Chen and X

    L. Chen and X. Huang. Complexes from complexes: finite element complexes in three dimensions.Math. Comp., 95(359):1083–1142, 2026

    work page 2026

  16. [16]

    S. H. Christiansen, J. Gopalakrishnan, J. Guzm´ an, and K. Hu. A discrete elasticity complex on three-dimensional Alfeld splits.Numer. Math., 156(1):159–204, 2024

    work page 2024

  17. [17]

    S. H. Christiansen and K. Hu. Finite element systems for vector bundles: elasticity and curvature.Found. Comput. Math., 23(2):545–596, 2023

    work page 2023

  18. [18]

    P. G. Ciarlet, L. Gratie, and C. Mardare. Intrinsic methods in elasticity: A mathematical survey.Discrete and Continuous Dynamical Systems, 2009

    work page 2009

  19. [19]

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

    work page 2020

  20. [20]

    Geymonat and F

    G. Geymonat and F. Krasucki. Some remarks on the compatibility conditions in elasticity.Accad. Naz. Sci. XL, 123:175–182, 2005

    work page 2005

  21. [21]

    S. Gong, J. Gopalakrishnan, J. Guzm´ an, and M. Neilan. Discrete elasticity exact sequences on Worsey-Farin splits. ESAIM Math. Model. Numer. Anal., 57(6):3373–3402, 2023

    work page 2023

  22. [22]

    Gopalakrishnan, J

    J. Gopalakrishnan, J. Guzm´ an, and J. J. Lee. The Johnson–Kˇ r´ ıˇ zek–Mercier elasticity element in higher dimensions. J. Numer. Math., 2025

    work page 2025

  23. [23]

    J. Hu. Finite element approximations of symmetric tensors on simplicial grids inR n: the higher order case.J. Comput. Math., 33(3):283–296, 2015

    work page 2015

  24. [24]

    J. Hu, R. Ma, and M. Zhang. A family of mixed finite elements for the biharmonic equations on triangular and tetrahedral grids.Sci. China Math., 64(12):2793–2816, 2021

    work page 2021

  25. [25]

    Hu and S

    J. Hu and S. Zhang. A family of symmetric mixed finite elements for linear elasticity on tetrahedral grids.Sci. China Math., 58(2):297–307, 2015

    work page 2015

  26. [26]

    Huang, C

    X. Huang, C. Zhang, Y. Zhou, and Y. Zhu. New low-order mixed finite element methods for linear elasticity.Adv. Comput. Math., 50(2):Paper No. 17, 31, 2024

    work page 2024

  27. [27]

    Johnson and B

    C. Johnson and B. Mercier. Some equilibrium finite element methods for two-dimensional elasticity problems.Numer. Math., 30(1):103–116, 1978

    work page 1978

  28. [28]

    Kˇ r´ ıˇ zek

    M. Kˇ r´ ıˇ zek. An equilibrium finite element method in three-dimensional elasticity.Apl. Mat., 27(1):46–75, 1982

    work page 1982

  29. [29]

    A. Seeger. Recent advances in the theory of defects in crystals.Physica Status Solidi (B), 1(7):669–698, 1961. Division of Applied Mathematics, Brown University, Box F, 182 George Street, Providence, RI 02912, USA Email address:johnny guzman@brown.edu School of Mathematics, Shanghai University of Finance and Economics, Shanghai 200433, China Email address...

    work page 1961