Uniformly $hp$-stable elements for the elasticity complex

Francis R. A. Aznaran,Kaibo Hu,Charles Parker

2024-09-26

Abstract:For the discretization of symmetric, divergence-conforming stress tensors in continuum mechanics, we prove inf-sup stability bounds which are uniform in polynomial degree and mesh size for the Hu--Zhang finite element in two dimensions. This is achieved via an explicit construction of a bounded right inverse of the divergence operator, with the crucial component being the construction of bounded Poincaré operators for the stress elasticity complex which are polynomial-preserving, in the Bernstein--Gelfand--Gelfand framework of the finite element exterior calculus. We also construct $hp$-bounded projection operators satisfying a commuting diagram property and $hp$-stable Hodge decompositions. Numerical examples are provided.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in continuum mechanics, for the discretization of symmetric and divergence - conforming stress tensors, how to prove that the Hu - Zhang finite element has uniform inf - sup stability with respect to polynomial order and mesh size in two - dimensional cases. Specifically, the author achieves this by constructing a bounded right inverse of the divergence operator, where the key part is to construct a bounded Poincaré operator that preserves polynomials in the Bernstein - Gelfand - Gelfand (BGG) framework. In addition, they also construct hp - bounded projection operators that satisfy the commutative diagram property and an hp - stable Hodge decomposition.

### Core problems of the paper
1. **Discretization of symmetric stress tensors**:
- In linear elasticity, while directly approximating deformation - description variables (such as displacement or velocity vector fields), some form of stress tensor (usually a matrix - valued dual variable) is also approximated. This follows the Hellinger - Reissner variational principle.
- For example, consider a linear isotropic elastic body \(\Omega\subset\mathbb{R}^d\) that deforms under the action of a body force \(f\), and the boundary \(\Gamma\) is divided into disjoint subsets \(\Gamma_N\) and \(\Gamma_D\). This problem can be described by the displacement \(u\) and the Cauchy stress tensor \(\sigma\) using the Hellinger - Reissner principle.

2. **Bounded right inverse of the divergence operator**:
- Construct a bounded right inverse of the divergence operator, ensuring that it is independent of the polynomial order \(p\) and the mesh size \(h\). This is achieved by constructing a bounded Poincaré operator that preserves polynomials.
- Specifically, the author proves the inf - sup stability of the Hu - Zhang finite element in two - dimensional cases and shows the superiority of this stability in high - order finite - element methods.

3. **Numerical experiment verification**:
- Provide numerical examples to verify the proposed theoretical results.

### Main contributions
1. **Uniform hp - stability**:
- Prove that the Hu - Zhang finite element has uniform inf - sup stability with respect to polynomial order and mesh size in two - dimensional cases.

2. **Construction of bounded Poincaré operator**:
- Construct a bounded Poincaré operator that preserves polynomials, which is achieved by the finite - element exterior calculus in the BGG framework.

3. **Projection operators and Hodge decomposition**:
- Construct hp - bounded projection operators that satisfy the commutative diagram property and show the existence of these projection operators and an hp - stable Hodge decomposition.

4. **Applications and extensions**:
- These results are applicable not only to linear elasticity problems but also to other problems that require the same inf - sup conditions, such as the linear Reissner - Mindlin plate problem.

Through these works, the author provides a theoretical basis and numerical tools for efficiently and stably solving complex mechanical problems involving symmetric stress tensors.