A nodally bound-preserving finite element method for hyperbolic convection-reaction problems

In this article, we present a numerical approach to ensure the preservation of physical bounds on the solutions to linear and nonlinear hyperbolic convection-reaction problems at the discrete level. We provide a rigorous framework for error analysis, formulating the discrete problem as a variational inequality and demonstrate optimal convergence rates in a natural norm. We summarise extensive numerical experiments validating the effectiveness of the proposed methods in preserving physical bounds and preventing unphysical oscillations, even in challenging scenarios involving highly nonlinear reaction terms.