Newton-Krylov continuation of periodic orbits for Navier-Stokes flows

Efficient numerical algorithms for the continuation of periodic orbits of high-dimensional dissipative dynamical systems, and for analyzing their stability are presented. They are based on shooting, Newton-Krylov and Arnoldi methods. A thermal convection fluid dynamics problem, which has a rich bifurcation diagram due to symmetries, has been used as test. After a pseudo-spectral discretization of the equations a system of dimension O(10(4)) has been obtained. The efficiency of the algorithms, which allows the unfolding of a complex diagram of periodic orbits, makes the methods suitable for the study of large nonlinear dissipative partial differential equations. (C) 2004 Elsevier Inc. All rights reserved.