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

We discuss an adaptive local refinement finite element method of lines for solving vector systems of parabolic partial differential equations on two-dimensional rectangular points. 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.