2010/03/23 by Lennaert van Veen, Genta KAWAHARA, van Veen, Lennaert +4 · 1 citation
Engineering · Mathematics · Physics and Astronomy · #34K19 #37N10 #65L10 #65P99 #Chaos control and synchronization #Computational Physics (physics.comp-ph) #Dynamical Systems (math.DS) #FOS: Mathematics #FOS: Physical sciences #Fluid Dynamics (physics.flu-dyn) #Fluid Dynamics and Turbulent Flows #Quantum chaos and dynamical systems #math.DS #msc:34K19 #msc:37N10 #msc:65L10 #msc:65P99 #physics.comp-ph #physics.flu-dyn
paper · pdf · doi:10.48550/arxiv.1003.4463
arxiv created 2010/03/23 · openalex publication_date 2010/03/23 · arxiv updated 2010/03/24 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
Recently, a flexible and stable algorithm was introduced for the computation of 2D unstable manifolds of periodic solutions to systems of ordinary differential equations. The main idea of this approach is to represent orbits in this manifold as the solutions of an appropriate boundary value problem. The boundary value problem is under determined and a one parameter family of solutions can be found by means of arclength continuation. This family of orbits covers a piece of the manifold. The quality of this covering depends on the way the boundary value problem is discretised, as do the tractability and accuracy of the computation. In this paper, we describe an implementation of the orbit continuation algorithm which relies on multiple shooting and Newton-Krylov continuation. We show that the number of time integrations necessary for each continuation step scales only with the number of shooting intervals but not with the number of degrees of freedom of the dynamical system. The number of shooting intervals is chosen based on linear stability analysis to keep the conditioning of the boundary value problem in check. We demonstrate our algorithm with two test systems: a low-order model of shear flow and a well-resolved simulation of turbulent plane Couette flow.