GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems
SIAM Journal on Scientific and Statistical Computing
Krylov methods for the incompressible Navier-Stokes equations
Journal of Computational Physics
Iterative Methods for Sparse Linear Systems
Iterative Methods for Sparse Linear Systems
Newton-Krylov continuation of periodic orbits for Navier-Stokes flows
Journal of Computational Physics
An overview of the Trilinos project
ACM Transactions on Mathematical Software (TOMS) - Special issue on the Advanced CompuTational Software (ACTS) Collection
Hi-index | 0.01 |
Recently, a flexible and stable algorithm was introduced for the computation of two-dimensional (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 (BVP). The BVP is underdetermined, 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 BVP is discretized, 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 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 BVP 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.