Blending low-order stabilised finite element methods : a positivity-preserving local projection method for the convection-diffusion equation

Barrenechea, Gabriel and Burman, Erik and Karakatsani, Fotini (2017) Blending low-order stabilised finite element methods : a positivity-preserving local projection method for the convection-diffusion equation. Computer Methods in Applied Mechanics and Engineering, 317. pp. 1169-1193. ISSN 0045-7825 (https://doi.org/10.1016/j.cma.2017.01.016)

[thumbnail of Barrenechea-etal-CMAME-2017-Blending-low-order-stabilised-finite-element-methods]
Preview
Text. Filename: Barrenechea_etal_CMAME_2017_Blending_low_order_stabilised_finite_element_methods.pdf
Accepted Author Manuscript
License: Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 logo

Download (1MB)| Preview

Abstract

In this work we propose a nonlinear blending of two low-order stabilisation mechanisms for the convection-diffusion equation. The motivation for this approach is to preserve monotonicity without sacrificing accuracy for smooth solutions. The approach is to blend a first-order artificial diffusion method, which will be active only in the vicinity of layers and extrema, with an optimal order local projection stabilisation method that will be active on the smooth regions of the solution. We prove existence of discrete solutions, as well as convergence, under appropriate assumptions on the nonlinear terms, and on the exact solution. Numerical examples show that the discrete solution produced by this method remains within the bounds given by the continuous maximum principle, while the layers are not smeared significantly.