Higher-order additive Runge–Kutta schemes for ordinary differential equations

By Christopher A. Kennedy, Mark H. Carpenter. Applied Numerical Mathematics, 2019

73 citations recorded by the source corpus.

Abstract Two new implicit–explicit, additive Runge–Kutta ( ARK 2 ) methods are given with fourth- and fifth-order formal accuracies, respectively. Both combine explicit Runge–Kutta (ERK) methods with explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) methods and include an embedded method for error control. The two methods have ESDIRKs which are internally L-stable on stages three and higher, have only modestly negative eigenvalues to the stage and step algebraic-stability matrices and have stage-order two. To improve computational efficiency, the fourth-order method has a diagonal coefficient of 0.1235. This is done to offset much of the extra computational cost of an extra stage by facilitating iterative convergence at each stage. Linear stability domains for both ERK methods have been made quite large and the dominant coupling stability term between the stability of the ESDIRK and ERK for very stiff modes has been removed. Though the fourth-order method is one of the best all-around fourth-order IMEX ARK 2 of which we are aware, the fifth-order method is likely best suited to mildly stiff problems with tight error tolerances. Methods are tested using the Van der Pol and Kaps' singular-perturbation problems. Results suggest that these new methods represent an improvement over existing methods of the same class.

View Paper: https://www.semanticscholar.org/paper/5acbdf194d662ef054c08beb1addbdb8740a7609

How researchers see this paper

The researchers whose trusted networks rank it highest, strongest advocate first. Every researcher gets a different order.

Citations in the corpus

Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers (2020, methodology): “For stif, nonstif, and mixed stif/nonstif problems, where f1 = f E may contain the nonstif components, which will be integrated using an explicit method, and f2 = f I the stif components, which will be integrated implicitly, ARKODE provides explicit, implicit, and implicit-explicit (IMEX) additive R”