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

      Barrenechea, Gabriel; Burman, Erik; Karakatsani, Fotini; University of Strathclyde; UCL; University of Chester (Elsevier, 2017-01-20)
      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.
    • Edge-based nonlinear diffusion for finite element approximations of convection–diffusion equations and its relation to algebraic flux-correction schemes

      Barrenechea, Gabriel; Burman, Erik; Karakatsani, Fotini; University of Strathclyde; University College London; University of Chester (Springer, 2016-05-07)
      For the case of approximation of convection–diffusion equations using piecewise affine continuous finite elements a new edge-based nonlinear diffusion operator is proposed that makes the scheme satisfy a discrete maximum principle. The diffusion operator is shown to be Lipschitz continuous and linearity preserving. Using these properties we provide a full stability and error analysis, which, in the diffusion dominated regime, shows existence, uniqueness and optimal convergence. Then the algebraic flux correction method is recalled and we show that the present method can be interpreted as an algebraic flux correction method for a particular definition of the flux limiters. The performance of the method is illustrated on some numerical test cases in two space dimensions.