pith. sign in

arxiv: 2404.02768 · v2 · submitted 2024-04-03 · 🧮 math.NA · cs.NA

Locking-free hybrid high-order method for linear elasticity

Carsten Carstensen , Ngoc Tien Tran This is my paper

Pith reviewed 2026-05-24 02:33 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords hybrid high-order methodlinear elasticitylocking-freea priori error analysisa posteriori error estimatessimplicial meshesincompressible limit
0
0 comments X

The pith

A single reconstruction operator for the Green strain yields locking-free hybrid high-order approximations for linear elasticity.

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

The paper presents a hybrid high-order scheme for linear elasticity that applies one reconstruction operator directly to the linear Green strain without any deviatoric-spherical split. This construction produces a priori error bounds whose equivalence constants remain independent of the Lamé parameter lambda, including in the nearly incompressible regime. The accompanying a posteriori error estimator is reliable and (up to data oscillations) efficient, requires no stabilization term, and stays robust as lambda grows. All analysis is performed on simplicial meshes so that the stabilization kernel contains conforming piecewise-polynomial elements. Numerical experiments on adaptive meshes confirm optimal convergence rates down to the incompressible limit and supply corresponding statements for the Stokes problem.

Core claim

The hybrid-high order scheme employs a single reconstruction operator for the linear Green strain and provides quasi-best approximation with lambda-independent equivalence constants. The reliable and efficient a posteriori error estimates are stabilization-free and lambda-robust on simplicial meshes.

What carries the argument

Single reconstruction operator for the linear Green strain, which replaces the classical split into deviatoric and spherical contributions.

If this is right

  • The method achieves optimal convergence rates independent of lambda.
  • A posteriori estimates guide adaptive mesh refinement effectively even in the incompressible limit.
  • The scheme applies directly to the Stokes problem with corresponding assertions.
  • Quasi-best approximation property holds with constants independent of lambda.

Where Pith is reading between the lines

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

  • Similar single-operator approaches might extend locking-free properties to other mixed problems or nonlinear elasticity.
  • The restriction to simplicial meshes suggests potential for generalization to polytopal meshes if kernel conditions can be met differently.
  • Empirical evidence from benchmarks could motivate theoretical extensions to three-dimensional cases or other boundary conditions.

Load-bearing premise

The error analysis requires simplicial meshes so that conforming piecewise polynomial finite elements lie in the kernel of the stabilization terms.

What would settle it

Numerical computation showing that the equivalence constants in the a priori estimates grow with increasing lambda would contradict the lambda-independent claim.

Figures

Figures reproduced from arXiv: 2404.02768 by Carsten Carstensen, Ngoc Tien Tran.

Figure 1
Figure 1. Figure 1: (a) Cook’s membrane and (b) convergence history plot of η for k = 1, . . . , 5 in Subsection 5.2 102 103 104 105 10−2 10−1 100 101 102 O(ndof−3/2 ) ndof η(ν = 0.3) η(ν = 0.4999) [PITH_FULL_IMAGE:figures/full_fig_p021_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (a) Adaptive triangulation into 592 triangles generated with k = 2 and (b) comparison between different ν with k = 2 in Subsection 5.2 5.2. Cook’s membrane. The tapered panel Ω := conv{(0, 0),(48, 44),(48, 60), (0, 44)} of [PITH_FULL_IMAGE:figures/full_fig_p021_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) Rotated L-shaped domain and (b) convergence history plot of ηe for k = 1, . . . , 5 in Subsection 5.3 102 103 104 105 10−7 10−5 10−3 10−1 101 O(ndof−1/4 ) O(ndof−3 ) O(ndof−1 ) ndof ∥σ − σh∥(k = 1) ∥σ − σh∥(k = 2) ∥σ − σh∥(k = 3) ∥σ − σh∥(k = 4) ∥σ − σh∥(k = 5) [PITH_FULL_IMAGE:figures/full_fig_p022_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (a) Convergence history plot of ∥σ − σh∥ for k = 1, . . . , 5 and (b) adaptive triangulation into 545 triangles generated with k = 2 in Subsection 5.3 102 103 104 105 0 1 2 ndof k = 1 k = 2 k = 3 k = 4 k = 5 102 103 104 105 10−4 10−3 10−2 10−1 100 101 O(ndof−3/2 ) ndof η˜(ν = 0.3) η˜(ν = 0.4999) ∥σ − σh∥(ν = 0.3) ∥σ − σh∥(ν = 0.4999) [PITH_FULL_IMAGE:figures/full_fig_p022_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (a) Efficiency index η/e ∥σ − σh∥ for k = 1, . . . , 5 and (b) comparison between different ν with k = 2 in Subsection 5.3 force f = 0 and g = 0 vanish but inhomogeneous Dirichlet boundary conditions apply with ∥σ − σh∥ ≲ ηe from (4.17). (In particular, the local contributions of osc(uD, FD) 2 are added to (5.1).) [PITH_FULL_IMAGE:figures/full_fig_p022_5.png] view at source ↗
read the original abstract

The hybrid-high order (HHO) scheme has many successful applications including linear elasticity as the first step towards computational solid mechanics. The striking advantage is the simplicity among other higher-order nonconforming schemes and its geometric flexibility as a polytopal method on the expanse of a parameter-free refined stabilization. This paper utilizes just one reconstruction operator for the linear Green strain and therefore does not rely on a split in deviatoric and spherical behaviour as in the classical HHO discretization. The a priori error analysis provides quasi-best approximation with $\lambda$-independent equivalence constants. The reliable and (up to data oscillations) efficient a posteriori error estimates are stabilization-free and $\lambda$-robust. The error analysis is carried out on simplicial meshes to allow conforming piecewise polynomials finite elements in the kernel of the stabilization terms. Numerical benchmarks provide empirical evidence for optimal convergence rates of the a posteriori error estimator in some associated adaptive mesh-refining algorithm also in the incompressible limit, where this paper provides corresponding assertions for the Stokes problem.

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 develops a hybrid high-order (HHO) discretization for linear elasticity that employs a single reconstruction operator applied directly to the symmetric gradient (linear Green strain), without the conventional deviatoric-volumetric split. On simplicial meshes the a priori analysis establishes a quasi-best-approximation result whose equivalence constants are independent of the Lamé parameter λ. The a posteriori analysis yields reliable and (up to data oscillations) efficient estimators that are free of stabilization parameters and remain robust as λ → ∞. Numerical experiments on adaptively refined meshes confirm optimal convergence rates, including in the incompressible limit, and the paper supplies corresponding statements for the Stokes problem.

Significance. If the central claims hold, the work supplies a technically simplified, parameter-free HHO scheme whose a priori constants and a posteriori estimators are explicitly λ-independent. The restriction to simplicial meshes is stated up front and is used to place conforming piecewise polynomials in the kernel of the stabilization, thereby securing the robustness; this is a clear, self-consistent scope choice rather than an internal inconsistency. The combination of a single reconstruction, stabilization-free estimators, and numerical evidence for adaptive performance in the incompressible regime constitutes a concrete advance for higher-order methods in computational solid mechanics.

minor comments (2)
  1. [Abstract] Abstract, final paragraph: the statement that 'corresponding assertions for the Stokes problem' are provided should be expanded in the main text to indicate whether these are complete proofs or reduced arguments obtained by specialization of the elasticity estimates.
  2. The notation for the single reconstruction operator (and its action on the symmetric gradient) should be introduced with an explicit equation number in the introduction so that later references to 'the reconstruction' are unambiguous.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary, significance assessment, and recommendation of minor revision. The report accurately captures the key contributions of the manuscript regarding the locking-free HHO scheme, λ-independent analysis, and stabilization-free a posteriori estimates on simplicial meshes. No major comments are provided in the report.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper derives quasi-best approximation and λ-robust a posteriori estimates from a single reconstruction operator on simplicial meshes, with the kernel property for conforming polynomials stated explicitly as the reason for the mesh restriction. No equation reduces a claimed prediction or constant to a fitted input by construction, and no load-bearing premise rests on a self-citation chain or ansatz smuggled from prior work by the same authors. The central claims remain independent of the target results.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The scheme rests on the standard HHO reconstruction framework and the restriction to simplicial meshes; no new free parameters are introduced because the stabilization is stated to be parameter-free.

axioms (2)
  • domain assumption Standard HHO reconstruction operators exist and satisfy the usual approximation properties
    Invoked implicitly for the error analysis on simplicial meshes.
  • domain assumption Meshes are simplicial
    Explicitly required to place conforming piecewise polynomials in the kernel of the stabilization.

pith-pipeline@v0.9.0 · 5697 in / 1215 out tokens · 25786 ms · 2026-05-24T02:33:12.544281+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Lower eigenvalue bounds with hybrid high-order methods

    math.NA 2024-06 unverdicted novelty 4.0

    Hybrid high-order eigensolvers compute guaranteed lower eigenvalue bounds with high-order convergence for linear elasticity and Steklov problems.

Reference graph

Works this paper leans on

44 extracted references · 44 canonical work pages · cited by 1 Pith paper

  1. [1]

    Alberty, C

    J. Alberty, C. Carstensen and S. A. Funken. Remarks around 50 lines of Matlab: short finite element implementation. Numer. Algorithms 20.2-3 (1999), 117–137

    work page 1999

  2. [2]

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

    work page 2002

  3. [3]

    Babuˇ ska and M

    I. Babuˇ ska and M. Suri. Locking effects in the finite element approximation of elasticity problems. Numer. Math. 62.4 (1992), 439–463

    work page 1992

  4. [4]

    Bertrand, C

    F. Bertrand, C. Carstensen, B. Gr¨ aßle and N. T. Tran. Stabilization-free HHO a posteriori error control. Numer. Math. 154.3-4 (2023), 369–408

    work page 2023

  5. [5]

    R. E. Bird, W. M. Coombs and S. Giani. A posteriori discontinuous Galerkin error estimator for linear elasticity. Appl. Math. Comput. 344/345 (2019), 78–96

    work page 2019

  6. [6]

    Boffi, F

    D. Boffi, F. Brezzi and M. Fortin. Mixed finite element methods and applications. Vol. 44. Springer, Heidelberg, 2013

    work page 2013

  7. [7]

    Botti, D

    M. Botti, D. A. Di Pietro and A. Guglielmana. A low-order nonconforming method for linear elasticity on general meshes. Comput. Methods Appl. Mech. Engrg. 354 (2019), 96– 118

    work page 2019

  8. [8]

    Botti, D

    M. Botti, D. A. Di Pietro and P. Sochala. A hybrid high-order method for nonlinear elasti- city. SIAM J. Numer. Anal. 55.6 (2017), 2687–2717

    work page 2017

  9. [9]

    D. Braess. Finite elements. Third. Theory, fast solvers, and applications in elasticity the- ory, Translated from the German by Larry L. Schumaker. Cambridge University Press, Cambridge, 2007, xviii+365

    work page 2007

  10. [10]

    S. C. Brenner. Korn’s inequalities for piecewise H 1 vector fields. Math. Comp. 73.247 (2004), 1067–1087

    work page 2004

  11. [11]

    S. C. Brenner and L. R. Scott. The mathematical theory of finite element methods. Third. Vol. 15. Springer, New York, 2008

    work page 2008

  12. [12]

    S. C. Brenner and L.-Y. Sung. Linear finite element methods for planar linear elasticity. Math. Comp. 59.200 (1992), 321–338

    work page 1992

  13. [13]

    Carstensen, M

    C. Carstensen, M. Eigel and J. Gedicke. Computational competition of symmetric mixed FEM in linear elasticity. Comput. Methods Appl. Mech. Engrg. 200.41-44 (2011), 2903– 2915

    work page 2011

  14. [14]

    Carstensen, D

    C. Carstensen, D. Gallistl and M. Schedensack. L2 best approximation of the elastic stress in the Arnold-Winther FEM. IMA J. Numer. Anal. 36.3 (2016), 1096–1119

    work page 2016

  15. [15]

    Carstensen and N

    C. Carstensen and N. Nataraj. A priori and a posteriori error analysis of the Crouzeix– Raviart and Morley FEM with original and modified right-hand sides. Comput. Methods Appl. Math. 21.2 (2021), 289–315

    work page 2021

  16. [16]

    Carstensen and M

    C. Carstensen and M. Schedensack. Medius analysis and comparison results for first-order finite element methods in linear elasticity. IMA J. Numer. Anal. 35.4 (2015), 1591–1621

    work page 2015

  17. [17]

    Carstensen and G

    C. Carstensen and G. Dolzmann. A posteriori error estimates for mixed FEM in elasticity. Numer. Math. 81.2 (1998), 187–209

    work page 1998

  18. [18]

    Carstensen, D

    C. Carstensen, D. Gallistl and J. Gedicke. Residual-based a posteriori error analysis for symmetric mixed Arnold-Winther FEM. Numer. Math. 142.2 (2019), 205–234

    work page 2019

  19. [19]

    Carstensen, B

    C. Carstensen, B. Gr¨ aß le and N. Nataraj. Unifying a posteriori error analysis of five piecewise quadratic discretisations for the biharmonic equation. J. Numer. Math. 32.1 (2024), 77–109

    work page 2024

  20. [20]

    Carstensen and N

    C. Carstensen and N. Heuer. A fractional-order trace-dev-div inequality. arXiv:2403.01291 (2024). 24 REFERENCES

    work page arXiv 2024

  21. [21]

    Carstensen, Q

    C. Carstensen, Q. Zhai and R. Zhang. A skeletal finite element method can compute lower eigenvalue bounds. SIAM J. Numer. Anal. 58.1 (2020), 109–124

    work page 2020

  22. [22]

    Cicuttin, A

    M. Cicuttin, A. Ern and N. Pignet. Hybrid high-order methods—a primer with applications to solid mechanics. Springer Briefs in Mathematics. Springer, Cham, 2021, viii+136

    work page 2021

  23. [23]

    D. A. Di Pietro and A. Ern. A hybrid high-order locking-free method for linear elasticity on general meshes. Comput. Methods Appl. Mech. Engrg. 283 (2015), 1–21

    work page 2015

  24. [24]

    D. A. Di Pietro, A. Ern and S. Lemaire. An arbitrary-order and compact-stencil discret- ization of diffusion on general meshes based on local reconstruction operators. Comput. Methods Appl. Math. 14.4 (2014), 461–472

    work page 2014

  25. [25]

    D. A. Di Pietro, A. Ern, A. Linke and F. Schieweck. A discontinuous skeletal method for the viscosity-dependent Stokes problem. Comput. Methods Appl. Mech. Engrg. 306 (2016), 175–195

    work page 2016

  26. [26]

    D. A. Di Pietro and J. Droniou. The hybrid high-order method for polytopal meshes. Vol. 19. MS&A. Modeling, Simulation and Applications. Design, analysis, and applications. Springer, 2020, xxxi+525

    work page 2020

  27. [27]

    D¨ orfler

    W. D¨ orfler. A convergent adaptive algorithm for Poisson’s equation.SIAM J. Numer. Anal. 33.3 (1996), 1106–1124

    work page 1996

  28. [28]

    Ern and J.-L

    A. Ern and J.-L. Guermond. Finite elements I—Approximation and interpolation. Vol. 72. Springer, 2021

    work page 2021

  29. [29]

    Ern and P

    A. Ern and P. Zanotti. A quasi-optimal variant of the hybrid high-order method for elliptic partial differential equations with H −1 loads. IMA J. Numer. Anal. 40.4 (2020), 2163– 2188

    work page 2020

  30. [30]

    Girault and P.-A

    V. Girault and P.-A. Raviart. Finite element methods for Navier-Stokes equations. Vol. 5. Theory and algorithms. Springer-Verlag, Berlin, 1986

    work page 1986

  31. [31]

    Kouhia and R

    R. Kouhia and R. Stenberg. A linear nonconforming finite element method for nearly in- compressible elasticity and Stokes flow.Comput. Methods Appl. Mech. Engrg. 124.3 (1995), 195–212

    work page 1995

  32. [32]

    Kreuzer and A

    C. Kreuzer and A. Veeser. Oscillation in a posteriori error estimation. Numer. Math. 148.1 (2021), 43–78

    work page 2021

  33. [33]

    P. L. Lederer and R. Stenberg. Analysis of weakly symmetric mixed finite elements for elasticity. Math. Comp. 93.346 (2024), 523–550

    work page 2024

  34. [34]

    Hybrid Discontinuous Galerkin methods for solving incompressible flow problems

    C. Lehrenfeld. “Hybrid Discontinuous Galerkin methods for solving incompressible flow problems”. PhD thesis. Rheinisch-Westf¨ alischen Technischen Hochschule Aachen, 2010

    work page 2010

  35. [35]

    J. M. Maubach. Local bisection refinement for n-simplicial grids generated by reflection. SIAM J. Sci. Comput. 16.1 (1995), 210–227

    work page 1995

  36. [36]

    Mora and G

    D. Mora and G. Rivera. A priori and a posteriori error estimates for a virtual element spectral analysis for the elasticity equations. IMA J. Numer. Anal. 40.1 (2020), 322–357

    work page 2020

  37. [37]

    R¨ ossle

    A. R¨ ossle. Corner singularities and regularity of weak solutions for the two-dimensional Lam´ e equations on domains with angular corners.J. Elasticity 60 (2000), 57–75

    work page 2000

  38. [38]

    L. R. Scott and S. Zhang. Finite element interpolation of nonsmooth functions satisfying boundary conditions. Math. Comp. 54.190 (1990), 483–493

    work page 1990

  39. [39]

    Stevenson

    R. Stevenson. The completion of locally refined simplicial partitions created by bisection. Math. Comput. 77.261 (2008), 227–241

    work page 2008

  40. [40]

    Veeser and P

    A. Veeser and P. Zanotti. Quasi-optimal nonconforming methods for symmetric elliptic problems. II—Overconsistency and classical nonconforming elements. SIAM J. Numer. Anal. 57.1 (2019), 266–292

    work page 2019

  41. [41]

    Veeser and P

    A. Veeser and P. Zanotti. Quasi-optimal nonconforming methods for symmetric elliptic problems. III—Discontinuous Galerkin and other interior penalty methods. SIAM J. Nu- mer. Anal. 56.5 (2018), 2871–2894

    work page 2018

  42. [42]

    Beir˜ ao da Veiga, F

    L. Beir˜ ao da Veiga, F. Brezzi, L. D. Marini and A. Russo. The virtual element method. Acta Numer. 32 (2023), 123–202

    work page 2023

  43. [43]

    Beir˜ ao da Veiga, C

    L. Beir˜ ao da Veiga, C. Canuto, R. H. Nochetto, G. Vacca and M. Verani. Adaptive VEM: stabilization-free a posteriori error analysis and contraction property. SIAM J. Numer. Anal. 61.2 (2023), 457–494

    work page 2023

  44. [44]

    Verf¨ urth

    R. Verf¨ urth. A posteriori error estimation techniques for finite element methods. Oxford University Press, Oxford, 2013. (C. Carstensen) Humboldt-Universit¨at zu Berlin, 10117 Berlin, Germany Email address: cc@math.hu-berlin.de (N. T. Tran) Universit¨at Augsburg, 86159 Augsburg, Germany Email address: ngoc1.tran@uni-a.de

    work page 2013