The implementation of hp-adaptivity is challenging as hanging nodes, edges, and faces have to be constrained to ensure compatibility of the shape functions. For this reason, most hp-code frameworks restrict themselves to 1-irregular meshes to ease the implementational effort. This work alleviates these difficulties by introducing a new formulation for high-order mesh adaptivity that provides full local hp-refinement capabilities at a comparably small implementational effort. Its main idea is the extension of the hp-d-method such that it allows for high-order overlay meshes yielding a hierarchical, multi-level hp-formulation of the Finite Element Method. This concept enables intuitive refinement and coarsening procedures, while linear independence and compatibility of the shape functions are guaranteed by construction. The proposed method is demonstrated to achieve exponential rates of convergence—both in terms of degrees of freedom and in run-time—for problems with non-smooth solutions. Furthermore, the scheme is used alongside the Finite Cell Method to simulate the heat flow around moving objects on a non-conforming background mesh and is combined with an energy-based refinement indicator for automatic hp-adaptivity.
Bibliographical noteFunding Information:
The first and the last author gratefully acknowledge the financial support of the German Research Foundation (DFG) under Grants RA 624/19-2 and RA 624/22-1.
© 2015, Springer-Verlag Berlin Heidelberg.
- Arbitrary hanging nodes
- Automatic hp-adaptivity
- Finite Cell Method
- High-order FEM