Finite element methods for second order differential equations with significant first derivatives

Abstract

Abstract Galerkin finite element methods based on symmetric pyramid basis functions give poor accuracy when applied to second order elliptic equations with large coefficients of the first order terms. This is particularly so when the mesh size is such that oscillations are present in the numerical solution. In the present note asymmetric linear and quadratic basis functions are introduced and shown to overcome this difficulty in an appropriate two point boundary value problem. In particular symmetric quadratic basis functions are oscillation free and highly accurate for a working range of mesh sizes.