On the convergence of higher-order finite element methods to weak solutions

Abstract

The ability to handle discontinuities appropriately is essential when solving nonlinear hyperbolic partial differential equations (PDEs). Discrete solutions to the PDE must converge to weak solutions in order for the discontinuity propagation speed to be correct. As shown by the Lax-Wendroff theorem, one method to guarantee that convergence, if it occurs, will be to a weak solution is to use a discretely conservative scheme. However, discrete conservation is not a strict requirement for convergence to a weak solution. This suggests a hierarchy of discretizations, where discretely conservative schemes are a subset of the larger class of methods that converge to the weak solution. We show here that a range of finite element methods converge to the weak solution without using discrete conservation arguments. The effect of using quadrature rules to approximate integrals is also considered. In addition, we show that solutions using non-conservation working variables also converge to weak solutions.