Issues with Positivity-Preserving Patankar-type Schemes

Abstract

Patankar-type schemes are linearly implicit time integration methods designed to be unconditionally positivity-preserving. However, there are only little results on their stability or robustness. We suggest two approaches to analyze the performance and robustness of these methods. In particular, we demonstrate problematic behaviors of these methods that, even on very simple linear problems, can lead to undesired oscillations and order reduction for vanishing initial condition. Finally, we demonstrate in numerical simulations that our theoretical results for linear problems apply analogously to nonlinear stiff problems.

Cited by in corpus (7)

• An Arbitrary High Order and Positivity Preserving Method for the Shallow Water Equations
• Modified Patankar Linear Multistep methods for production-destruction systems
• A high-order, fully well-balanced, unconditionally positivity-preserving finite volume framework for flood simulations
• Structure-Preserving Numerical Methods for Fokker-Planck Equations
• A Boot-Strapping Technique to Design Dense Output Formulae for Modified Patankar-Runge-Kutta Methods
• Modified Patankar Semi-Lagrangian Scheme for the Optimal Control of Production-Destruction systems
• The Lax-Wendroff theorem for Patankar-type methods applied to hyperbolic conservation laws