Stabilized Mixed Finite Element Methods for Linear Elasticity on Simplicial Grids in ℝn

Abstract In this paper, we design two classes of stabilized mixed finite element methods for linear elasticity on simplicial grids. In the first class of elements, we use 𝑯 ⁢ ( div , Ω ; 𝕊 ) ${\boldsymbol{H}(\operatorname{div},\Omega;\mathbb{S})}$ - P k ${P_{k}}$ and 𝑳 2 ⁢ ( Ω ; ℝ n ) ${\boldsymbol{L}^{2}(\Omega;\mathbb{R}^{n})}$ - P k - 1 ${P_{k-1}}$ to approximate the stress and displacement spaces, respectively, for 1 ≤ k ≤ n ${1\leq k\leq n}$ , and employ a stabilization technique in terms of the jump of the discrete displacement over the edges/faces of the triangulation under consideration; in the second class of elements, we use 𝑯 0 1 ⁢ ( Ω ; ℝ n ) ${\boldsymbol{H}_{0}^{1}(\Omega;\mathbb{R}^{n})}$ - P k ${P_{k}}$ to approximate the displacement space for 1 ≤ k ≤ n ${1\leq k\leq n}$ , and adopt the stabilization technique suggested by Brezzi, Fortin, and Marini [19]. We establish the discrete inf-sup conditions, and consequently present the a priori error analysis for them. The main ingredient for the analysis are two special interpolation operators, which can be constructed using a crucial 𝑯 ⁢ ( div ) ${\boldsymbol{H}(\operatorname{div})}$ bubble function space of polynomials on each element. The feature of these methods is the low number of global degrees of freedom in the lowest order case. We present some numerical results to demonstrate the theoretical estimates.