30.1NAMay 13
$\ell$FEM: An efficient loop-free Matlab implementation of isoparametric bulk and surface finite elementsBalázs Kovács, Michael Lantelme
The $\ell$FEM MATLAB package provides a simple, efficient, and flexible implementation of isoparametric finite elements in bulk domains and on surfaces. The finite element matrix assemblies are based on MATLAB's paged operators and therefore completely loop-free. We give a short and conscious description of high-order isoparametric surface finite elements, which is then used to describe the assembly process and the implementation. We report on relevant numerical experiments (runtime comparisons, modifications for non-linear problems, etc.), and on additional functions, examples, and a testing unit which are all part of the $\ell$FEM package.
20.6NAMar 30
A posteriori error estimates for parabolic partial differential equations on stationary surfacesBalázs Kovács, Michael Lantelme
This paper develops and discusses a residual-based a posteriori error estimator for parabolic surface partial differential equations on closed stationary surfaces. The full discretization uses the surface finite element method in space and the backward Euler method in time. The proposed error indicator bounds the error quantities globally in space from above and below, and globally in time from above and locally from below. Based on the derived error indicator, a space-time adaptive algorithm is proposed. Numerical experiments illustrate and complement the theory.