TY - JOUR
T1 - Issues with positivity-preserving Patankar-type schemes
AU - Torlo, Davide
AU - Öffner, Philipp
AU - Ranocha, Hendrik
N1 - KAUST Repository Item: Exported on 2022-09-14
Acknowledgements: D. T. was funded by Team CARDAMOM in Inria–Bordeaux Sud–Ouest, France and by a SISSA Mathematical Fellowship, Italy. P.Ö. gratefully acknowledge support of the Gutenberg Research College, JGU Mainz and the UZH Postdoc Scholarship (Number FK-19-104). H. R. was supported by the German Research Foundation (DFG, Deutsche Forschungsgemeinschaft) under Germany's Excellence Strategy EXC 2044-390685587, Mathematics Münster: Dynamics-Geometry-Structure. We would like to thank Stefan Kopecz and David Ketcheson for fruitful discussion at the beginning of this project. This project has started with the visit by H.R. in Zurich in 2019 which was supported by the SNF project (Number 175784) and the King Abdullah University of Science and Technology (KAUST).
This publication acknowledges KAUST support, but has no KAUST affiliated authors.
PY - 2022/8/17
Y1 - 2022/8/17
N2 - 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.
AB - 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.
UR - http://hdl.handle.net/10754/680560
UR - https://linkinghub.elsevier.com/retrieve/pii/S016892742200188X
UR - http://www.scopus.com/inward/record.url?scp=85135957395&partnerID=8YFLogxK
U2 - 10.1016/j.apnum.2022.07.014
DO - 10.1016/j.apnum.2022.07.014
M3 - Article
SN - 0168-9274
VL - 182
SP - 117
EP - 147
JO - Applied Numerical Mathematics
JF - Applied Numerical Mathematics
ER -