The Sparse-Grid-Based Adaptive Spectral Koopman Method

SIAM Journal on Scientific Computing, Volume 46, Issue 5, Page A2925-A2950, October 2024.
Abstract. The adaptive spectral Koopman (ASK) method was introduced to numerically solve autonomous dynamical systems that laid the foundation for numerous applications across different fields in science and engineering. Although ASK achieves high accuracy, it is computationally more expensive for multidimensional systems compared with conventional time integration schemes like Runge–Kutta. In this work, we combine the sparse grid and ASK to accelerate the computation for multidimensional systems. This sparse-grid-based ASK (SASK) method uses the Smolyak structure to construct multidimensional collocation points as well as associated polynomials that are used to approximate eigenfunctions of the Koopman operator of the system. In this way, the number of collocation points is reduced compared with using the tensor product rule. We demonstrate that SASK can be used to solve ordinary differential equations (ODEs) and partial differential equations (PDEs) based on their semidiscrete forms. Numerical experiments are illustrated to compare the performance of SASK and state-of-the-art ODE solvers.