High-order parameterization of (un)stable manifolds for hybrid maps: Implementation and applications

In this work we study, from a numerical point of view, the (un)stable manifolds of a certain class of dynamical systems called hybrid maps. The dynamics of these systems are generated by a two stage procedure: the first stage is continuous time advection under a given vector field, the second stage is discrete time advection under a given diffeomorphism. Such hybrid systems model physical processes where a differential equation is occasionally kicked by a strong disturbance. We propose a numerical method for computing local (un)stable manifolds, which leads to high order polynomial parameterization of the embedding. The parameterization of the invariant manifold is not the graph of a function and can follow folds in the embedding. Moreover we obtain a representation of the dynamics on the manifold in terms of a simple conjugacy relation. We illustrate the utility of the method by studying a planar example system.