Multistep and Runge-Kutta convolution quadrature methods for coupled dynamical systems

We consider the efficient numerical solution of coupled dynamical systems, consisting of a low dimensional nonlinear part and a high dimensional linear time invariant part, e.g., stemming from spatial discretization of an underlying partial differential equation. The linear subsystem can be eliminated in frequency domain and for the numerical solution of the resulting integro-differential algebraic equations, we propose a combination of Runge-Kutta or multistep time stepping methods with appropriate convolution quadrature to handle the integral terms. The resulting methods are shown to be algebraically equivalent to a Runge-Kutta or multistep solution of the coupled system and thus automatically inherit the corresponding stability and accuracy properties. After a computationally expensive pre-processing step, the online simulation can, however, be performed at essentially the same cost as solving only the low dimensional nonlinear subsystem. The proposed method is, therefore, particularly attractive, if repeated simulation of the coupled dynamical system is required. (C) 2019 Elsevier B.V. All rights reserved.

Automated summary

This study introduces an efficient numerical solution method for coupled dynamical systems, which eliminates the linear subsystem and handles integral terms with appropriate convolution methods to maintain stability and accuracy. Despite the high computational cost of pre-processing, online simulation can be performed at essentially the same cost, making it particularly suitable for situations requiring repeated simulations.