Sensitivities for topological graph changes in finite element meshes and application to adaptive refinement

We propose a novel approach to adaptivity in FEM based on local sensitivities for topological mesh changes. To this end, we consider refinement as a continuous operation on the edge graph of the finite element discretization, for instance by splitting nodes along edges and expanding edges to elements. Thereby, we introduce the concept of a topological mesh derivative for a given objective function that depends on the discrete solution of the underlying PDE. These sensitivities may in turn be used as refinement indicators within an adaptive algorithm. For their calculation, we rely on the first-order asymptotic expansion of the Galerkin solution with respect to the topological mesh change. As a proof of concept, we consider the total potential energy of a linear symmetric second-order elliptic PDE, minimization of which is known to decrease the approximation error in the energy norm. In this case, our approach yields local sensitivities that are closely related to the reduction of the energy error upon refinement and may therefore be used as refinement indicators in an adaptive algorithm. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)