Design of High-Order Decoupled Multirate GARK Schemes

Multirate time integration methods apply different step sizes to resolve different components of a system based on local activity levels. This local selection of step sizes allows increased computational efficiency while achieving the desired solution accuracy. While the multirate idea is elegant and has been around for decades, multirate methods are not yet widely used in applications. This is due, in part, to the difficulties created by the construction of high-order multirate schemes. Seeking to overcome these challenges, this work focuses on the design of practical high-order multirate methods using the theoretical framework of multirate general-structure additive Runge--Kutta (MrGARK) methods [M. Gunther and A. Sandu, Numer. Math., 133 (2016), pp. 497--524], which provides generic order conditions and linear and nonlinear stability analyses. A set of design criteria for practical multirate methods is defined herein: method coefficients should be generic in the step size ratio but should not depend st...