Second-order finite element approximations and a posteriori error estimation for two-dimensional parabolic systems

Summary We discuss an adaptive local refinement finite element method of lines for solving vector systems of parabolic partial differential equations on two-dimensional rectangular regions. The partial differential system is discretized in space using a Galerkin approach with piecewise eight-node serendipity functions. An a posteriori estimate of the spatial discretization error of the finite element solution is obtained using piecewise fifth degree polynomials that vanish on the edges of the rectangular elements of a grid. Ordinary differential equations for the finite element solution and error estimate are integrated in time using software for stiff differential systems. The error estimate is used to control a local spatial mesh refinement procedure in an attempt to keep a global measure of the error within prescribed limits. Examples appraising the accuracy of the solution and error estimate and the computational efficiency of the procedure relative to one using bilinear finite elements are presented.Dedicated to Prof. Ivo Babuka on the occasion of his 60th birthdayThis research was partially supported by the U.S. Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant Number AFOSR 85-0156 and the U.S. Army Research Office under Contract Number DAAL 03-86-K-0112