Computation of slow manifolds using the isostable coordinate system

Koopman analysis can be used to understand the dynamics of a nonlinear dynamical system in terms of a linear, but generally infinite-dimensional operator. The isostable coordinate system focuses on the slowest decaying principal Koopman eigenmodes. This work leverages the isostable coordinate framework in the identification of slow manifolds for dynamical systems with fixed point attractors, defined as surfaces for which the fastest decaying isostable coordinates are zero. Numerical challenges associated with separation between fast and slow timescales necessitate the development of new computational approaches to identify these slow manifolds. Two such strategies are developed which approximate backward-time solutions on the slow manifold starting near the fixed point and extending far beyond the linear regime. Application to a variety of examples illustrates the utility of these methods and their potential use for model order reduction purposes.