Many-Stage Optimal Stabilized Runge–Kutta Methods for Hyperbolic Partial Differential Equations

Journal of Scientific Computing, vol. 99

Abstract

Abstract A novel optimization procedure for the generation of stability polynomials of stabilized explicit Runge–Kutta methods is devised. Intended for semidiscretizations of hyperbolic partial differential equations, the herein developed approach allows the optimization of stability polynomials with more than hundred stages. A potential application of these high degree stability polynomials are problems with locally varying characteristic speeds as found for non-uniformly refined meshes and spatially varying wave speeds. To demonstrate the applicability of the stability polynomials we construct 2N-storage many-stage Runge–Kutta methods that match their designed second order of accuracy when applied to a range of linear and nonlinear hyperbolic PDEs with smooth solutions. These methods are constructed to reduce the amplification of round off errors which becomes a significant concern for these many-stage methods.

Authors 3

1. Daniel Doehring corresponding Aachen Applied and Computational Mathematics

   RWTH Aachen University

   Affiliation as printed

   Applied and Computational Mathematics, RWTH Aachen University, Schinkelstrasse 2, 52062, Aachen, North Rhine-Westphalia, Germany

2. Gregor J. Gassner

   University of Cologne

   Affiliation as printed

   Department of Mathematics and Computer Science, Center for Data and Simulation Science, University of Cologne, Weyertal 86-90, 50931, Cologne, North Rhine-Westphalia, Germany

3. Manuel Torrilhon Aachen Applied and Computational Mathematics

   RWTH Aachen University

   Affiliation as printed

   Applied and Computational Mathematics, RWTH Aachen University, Schinkelstrasse 2, 52062, Aachen, North Rhine-Westphalia, Germany