Abstract

This paper presents for the first time the derivation of an hp a posteriori error estimator for the symmetric interior penalty discontinuous Galerkin finite element method for linear elastic analysis. Any combination of Neumann and Dirichlet boundary conditions are admissible in the formulation, including applying Neumann and Dirichlet on different components on the same region of the boundary. Therefore, the error estimator is applicable to a variety of physical problems. The error estimator is incorporated into an hp-adaptive finite element solver and verified against smooth and non-smooth problems with closed-form analytical solutions, as well as, being demonstrated on a non-smooth problem with complex boundary conditions. The hp-adaptive finite element analyses achieve exponential rates of convergence. The performances of the hp-adaptive scheme are contrasted against uniform and adaptive h refinement. This paper provides a complete framework for adaptivity in the symmetric interior penalty discontinuous Galerkin finite element method for linear elastic analysis.

Highlights

  • The discretization of partial differential equations for numerical computation facilitates the solution of physical systems, it introduces approximation errors

  • We introduce the continuous weak formulation of problem (1), this formulation is only used in the analysis for the error estimator and it is not implemented in the code

  • We introduce the polynomial degrees for the approximation in our Discontinuous Galerkin (DG) method

Read more

Summary

Introduction

The discretization of partial differential equations for numerical computation facilitates the solution of physical systems, it introduces approximation errors. We present a new error estimator for linear elasticity based on a Discontinuous Galerkin (DG) finite element method. DG methods were introduced in the early 1970s as a way to numerically solve first-order hyperbolic problems [2]. A variety of DG methods have been developed in the decades [6], among them is the subclass of interior penalty DG methods which are stable and, in the authors’ opinion, easy to implement. In [7] the first residual-based energy-norm error estimator for the symmetric interior penalty discontinuous Galerkin finite element method is presented. The symmetric interior penalty discontinuous Galerkin finite element method is introduced, and Section 5 provides the a priori convergence results of the method.

Model problem
Symmetric interior penalty discontinuous Galerkin method
Reliability of the error estimator
Efficiency of the error estimator
Numerical examples
Smooth solution problem
L-shaped non-smooth problem
Crack in a plate problem
Conclusions
Findings
Methods
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call