An operator-customized wavelet-finite element approach for the adaptive solution of second-order partial differential equations on unstructured meshes

Abstract

The Finite Element Method (FEM) is a widely popular method for the numerical solution of Partial Differential Equations (PDE), on multi-dimensional unstructured meshes. Lagrangian finite elements, which preserve C⁰ continuity with interpolating piecewise-polynomial shape functions, are a common choice for second-order PDEs. Conventional single-scale methods often have difficulty in efficiently capturing fine-scale behavior (e.g. singularities or transients), without resorting to a prohibitively large number of variables. This can be done more effectively with a multi-scale method, such as the Hierarchical Basis (HB) method. However, the HB FEM generally yields a multi-resolution stiffness matrix that is coupled across scales. We propose a powerful generalization of the Hierarchical Basis: a second-generation wavelet basis, spanning a Lagrangian finite element space of any given polynomial order.