This paper proposes and analyzes a weak Galerkin (WG) finite element method with strong symmetric stresses for two- and three-dimensional linear elasticity problems on conforming or nonconforming polygon/polyhedral meshes. The WG method uses piecewise-polynomial approximations of degrees k ( ≥ 1 ${\geq 1}$ ) for the stress, k + 1 ${k+1}$ for the displacement, and k for the displacement trace on the inter-element boundaries. It is shown to be equivalent to a hybridizable discontinuous Galerkin (HDG) finite element scheme. We show that the WG methods are robust in the sense that the derived a priori error estimates are optimal and uniform with respect to the Lamé constant λ. Numerical experiments confirm the theoretical results.
1 December 2018 in Journal of Computational and Applied Mathematics
In this paper, we study a hybridizable discontinuous Galerkin (HDG) method for the Maxwell operator. The only global unknowns are defined on the inter-element boundaries, and the numerical solutions are obtained by using discontinuous polynomial approximations. The error analysis is based on a mixed curl-curl formulation for the Maxwell equations. Theoretical results are obtained under a more general regularity requirement. In particular for the low regularity case, special treatment is applied to approximate data on the boundary. The HDG method is shown to be stable and convergence in an optimal order for both high and low regularity cases. Numerical experiments with both smooth and singular analytical solutions are performed to verify the theoretical results.