Strong stability preserving (SSP) Runge-Kutta methods are often desired when evolving in time problems that have two components that have very different time scales. Where the SSP property is needed, it has been shown that implicit and implicit-explicit methods have very restrictive time-steps and are therefore not efficient. For this reason, SSP integrating factor methods may offer an attractive alternative to traditional time-stepping methods for problems with a linear component that is stiff and a nonlinear component that is not. However, the strong stability properties of integrating factor Runge-Kutta methods have not been established. In this work we show that it is possible to define explicit integrating factor Runge-Kutta methods that preserve the desired strong stability properties satisfied by each of the two components when coupled with forward Euler time-stepping, or even given weaker conditions. We define sufficient conditions for an explicit integrating factor Runge-Kutta method to be SSP, namely that they are based on explicit SSP Runge-Kutta methods with non-decreasing abscissas. We find such methods of up to fourth order and up to ten stages, analyze their SSP coefficients, and prove their optimality in a few cases. We test these methods to demonstrate their convergence and to show that the SSP time-step predicted by the theory is generally sharp, and that the non-decreasing abscissa condition is needed in our test cases. Finally, we show that on typical total variation diminishing linear and nonlinear test-cases our new explicit SSP integrating factor Runge-Kutta methods out-perform the corresponding explicit SSP Runge-Kutta methods, implicit-explicit SSP Runge-Kutta methods, and some well-known exponential time differencing methods.
Problems with components that feature significantly different time scales, where the stiff time-step restriction comes from a linear component, implicit-explicit (IMEX) methods alleviate this restriction if the concern is linear stability. However, when nonlinear non-inner-product stability properties are of interest, such as in the evolution of hyperbolic partial differential equations with shocks or sharp gradients, linear inner-product stability is no longer sufficient for convergence, and so strong stability preserving (SSP) methods are often needed. Where the SSP property is needed, IMEX SSP Runge-Kutta (SSP-IMEX) methods have very restrictive time-steps. An alternative to SSP-IMEX schemes is to adopt an integrating factor approach to handle the linear component exactly and step the transformed problem forward using some time-evolution method. The strong stability properties of integrating factor Runge-Kutta methods were established in [15], where it was shown that it is possible to define explicit integrating factor Runge-Kutta methods that preserve strong stability properties satisfied by each of the two components when coupled with forward Euler time-stepping. It was proved that the solution will be SSP if the transformed problem is stepped forward with an explicit SSP Runge-Kutta method that has non-decreasing abscissas. However, explicit SSP Runge-Kutta methods have an order barrier of p = 4, and sometimes higher order is desired. In this work we consider explicit SSP two-step Runge-Kutta integrating factor methods to raise the order. We show that strong stability is ensured if the two-step Runge-Kutta method used to evolve the transformed problem is SSP and has non-decreasing abscissas. We find such methods up to eighth order and present their SSP coefficients. Adding a step allows us to break the fourth order barrier on explicit SSP Runge-Kutta methods; furthermore, our explicit SSP two-step Runge-Kutta methods with non-decreasing abscissas typically have larger SSP coefficients than the corresponding one-step methods. A selection of our methods are tested for convergence and demonstrate the design order. We also show, for selected methods, that the SSP time-step predicted by the theory is a lower bound of the allowable time-step for linear and nonlinear problems that satisfy the total variation diminishing (TVD) condition. We compare some of the non-decreasing abscissa SSP two-step Runge-Kutta methods to previously found methods that do not satisfy this criterion on linear and nonlinear TVD test cases to show that this non-decreasing abscissa condition is indeed necessary in practice as well as theory. We also compare these results to our SSP integrating factor Runge-Kutta methods designed in [15]. This paper is dedicated to the memory of Saul Abarbanel. His wisdom, humor, and kindness were a gift to all who knew him.
Strong stability preserving (SSP) Runge-Kutta methods are desirable when evolving in time problems that have discontinuities or sharp gradients and require nonlinear non-inner-product stability properties to be satisfied. Unlike the case for L 2 linear stability, implicit methods do not significantly alleviate the time-step restriction when the SSP property is needed. For this reason, when handling problems with a linear component that is stiff and a nonlinear component that is not, SSP integrating factor Runge-Kutta methods may offer an attractive alternative to traditional time-stepping methods. The strong stability properties of integrating factor Runge-Kutta methods where the transformed problem is evolved with an explicit SSP Runge-Kutta method with non-decreasing abscissas was recently established. However, these methods typically have smaller SSP coefficients (and therefore a smaller allowable time-step) than the optimal SSP Runge-Kutta methods, which often have some decreasing abscissas. In this work, we consider the use of downwinded spatial operators to preserve the strong stability properties of integrating factor Runge-Kutta methods where the Runge-Kutta method has some decreasing abscissas. We present the SSP theory for this approach and present numerical evidence to show that such an approach is feasible and performs as expected. However, we also show that in some cases the integrating factor approach with explicit SSP Runge-Kutta methods with non-decreasing abscissas performs nearly as well, if not better, than with explicit SSP Runge-Kutta methods with downwinding. In conclusion, while the downwinding approach can be rigorously shown to guarantee the SSP property for a larger time-step, in practice using the integrating factor approach by including downwinding as needed with optimal explicit SSP Runge-Kutta methods does not necessarily provide significant benefit over using explicit SSP Runge-Kutta methods with non-decreasing abscissas.
We develop and use a novel mixed-precision weighted essentially nonoscillatory (WENO) method for solving the Teukolsky equation, which arises when modeling perturbations of Kerr black holes. We show that WENO methods outperform higher-order finite-difference methods, standard in the discretization of the Teukolsky equation, due to the need to add dissipation for stability purposes in the latter. In particular, as the WENO scheme uses no additional dissipation it is well-suited for scenarios requiring long-time evolution such as the study of Price tails and gravitational wave emission from extreme mass ratio binaries. In the mixed-precision approach, the expensive computation of the WENO weights is performed in reduced floating-point precision that results in a significant speedup factor of ≈ 3.3. In addition, we use state-of-the-art Nvidia general-purpose graphics processing units and cluster parallelism to further accelerate the WENO computations. Our optimized WENO solver can be used to quickly generate accurate results of significance in the field of black hole and gravitational wave physics. We apply our solver to study the behavior of the Aretakis charge -a conserved quantity, that if detected by a gravitational wave observatory like LIGO/Virgo would prove the existence of extremal black holes.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.