An Adaptive Newton-Picard Algorithm with Subspace Iteration for Computing Periodic Solutions

This paper is concerned with the efficient computation of periodic orbits in large-scale dynamical systems that arise after spatial discretization of partial differential equations (PDEs). A hybrid Newton--Picard scheme based on the shooting method is derived, which in its simplest form is the recursive projection method (RPM) of Shroff and Keller [SIAM J. Numer. Anal., 30 (1993), pp. 1099--1120] and is used to compute and determine the stability of both stable and unstable periodic orbits. The number of time integrations needed to obtain a solution is shown to be determined only by the system's dynamics. This contrasts with traditional approaches based on Newton's method, for which the number of time integrations grows with the order of the spatial discretization. Two test examples are given to show the performance of the methods and to illustrate various theoretical points.

[1]  William L. Kath,et al.  Bifurcation with memory , 1986 .

[2]  W. Rheinboldt.  Numerical analysis of parametrized nonlinear equations , 1986 .

[3]  W. Mackens,et al.  Solving large nonlinear systems of equations by an adaptive condensation process , 1986 .

[4]  H. Jarausch,et al.  Computing Bifurcation Diagrams for Large Nonlinear Variational Problems , 1987 .

[5]  K. Meerbergen,et al.  Matrix transformations for computing rightmost eigenvalues of large sparse non-symmetric eigenvalue problems , 1996 .

[6]  Wolfgang Mackens,et al.  Numerical Treatment of Bifurcation Branches by Adaptive Condensation , 1984 .

[7]  N. G. Parke,et al.  Ordinary Differential Equations. , 1958 .

[8]  G. Stewart.  Simultaneous iteration for computing invariant subspaces of non-Hermitian matrices , 1976 .

[9]  Alan R. Champneys,et al.  A Newton-Picard shooting method for computing periodic solutions of large-scale dynamical systems , 1995 .

[10]  Milan Kubíček,et al.  Continuation of Periodic Solutions in Parabolic Partial Differential Equations , 1987 .

[11]  Allan D. Jepson,et al.  The computation and use of Floquet multipliers for bifurcation analysis , 1994 .

[12]  R. Seydel.  Practical Bifurcation and Stability Analysis , 1994 .

[13]  Allan D. Jepson,et al.  O. K. Floquet multipliers , 1991 .

[14]  Y. Kuznetsov,et al.  Continuation techniques and interactive software for bifurcation analysis of ODEs and iterated maps: physics , 1993 .

[15]  Karl Meerbergen,et al.  Matrix transformations for computing rightmost eigenvalues of large sparse nonsymmetric matrices , 1994 .

[16]  H. Jarausch,et al.  Analyzing stationary and periodic solutions of systems of parabolic partial differential equations by using singular subspaces as reduced basis , 1994 .

[17]  W. Rheinboldt.  Numerical analysis of continuation methods for nonlinear structural problems , 1981 .

[18]  Gautam M. Shroff,et al.  You have printed the following article : Stabilization of Unstable Procedures : The Recursive Projection Method , 2007 .