Abstract

This paper describes a combination of automatic $hp$-adaptive finite elements and domain decomposition. The combination is based on the Bank-Holst parallel adaptive meshing paradigm. The $hp$-adaptivity is based on a derivative recovery technique while the domain decomposition method is formulated using a mortar-like formulation with Dirac $\delta$ functions for the mortar element space. Numerical results show that the approach can scale to hundreds of processors and the convergence of the domain decomposition method is independent of the number of processors as well as the distribution of element sizes and element degrees among the local meshes.