On the Discrete Maximum Principle for the Local Projection Scheme with Shock Capturing

Abstract

It is a well known fact that finite element solutions of convection dominated problems can exhibit spurious oscillations in the vicinity of boundary layers. One way to overcome this numerical instability is to use schemes that satisfy the discrete maximum principle. There are monotone methods for piecewise linear elements on simplices based on the upwind techniques or artificial diffusion. In order to satisfy the discrete maximum principle for the local projection scheme, we add an edge oriented shock capturing term to the bilinear form. The analysis of the proposed stabilisation method is complemented with numerical examples in 2D.