Nullspaces yield new explicit Runge–Kutta pairs
J. H. Verner
DOI: https://doi.org/10.1007/s11075-023-01731-6
IF: 2.37
2024-02-17
Numerical Algorithms
Abstract:Sixty years ago, Butcher (Butcher Math. Soc. 3, 185–201 1963) characterized a natural tabulation of the order conditions for Runge–Kutta methods of order p as an isomorphism from the set of rooted trees having up to p nodes and provided examples of explicit and implicit methods of several orders. Within a few years, Fehlberg (Fehlberg 1968) derived pairs of explicit methods of successive orders that could be implemented efficiently by using the difference of each pair of estimates to control the local error. Unfortunately, Fehlberg's pairs were deficient for quadrature problems. Subsequently, this author (Verner SIAM J. Numer. Anal. 15, 772–790 1978, Verner SIAM J. Numer. Anal. 16, 857–875 1979) derived parametric families of explicit Runge–Kutta pairs of increasing orders 6 to 9 that avoided this problem altogether. These, and most known explicit methods, have been derived by exploiting certain "simplifying conditions" suggested by Butcher (Butcher Math. Soc. 3, 185–201 1963, Butcher 1987 [pp. 194 f.]) that imposed constraints on subsets of the coefficients and thereby simplified the solution of the order conditions for moderate to high-order methods. "RK-Test-21", a MAPLE program developed recently by Butcher (Butcher 2021), was applied to derive known 13-stage pairs of orders 7 and 8. Unexpectedly, results of this application revealed the existence of some previously unknown pairs—i.e., some with the correct orders that satisfied most, but not all, of the previously known simplifying conditions. This present study develops formulas for directly computing exact coefficients of these new pairs together with others lying within this new parametric family of (13,7-8) pairs. While the best of these new pairs falls short of the best of pairs already known, the properties discovered might be utilized to precisely characterize recently reported higher-order methods found using other approaches by Khashin (Khashin 2009) and Zhang (Zhang 2019) and possibly lead to finding other Runge–Kutta and related yet unknown methods.
mathematics, applied