×

A Newton-Picard shooting method for computing periodic solutions of large-scale dynamical systems. (English) Zbl 1080.65542

Summary: A numerical method is presented for the efficient computation and continuation of periodic solutions of large systems of ordinary differential equations. The method is particularly useful when a shooting approach based on full Newton or quasi-Newton iteration is prohibitively expensive. Both stable and unstable solutions can be computed. The basic idea is to split the eigenspace of the monodromy matrix around the periodic orbit into two orthogonal components which contain, respectively, all eigenvalues with norm less than or greater than some threshold value less than unity. In the former ‘stable’ subspace one performs time integration, which is equivalent to a Picard iteration; in the much smaller, ‘unstable’ subspace a Newton step is carried out. A strategy for computing the stable and unstable subspaces without forming the full monodromy matrix is briefly discussed. Numerical results are presented for the one-dimensional Brusselator model. The Newton-Picard method with an appropriately chosen threshold value proves to be accurate and much more efficient than shooting based on full Newton iteration.

MSC:

65L10 Numerical solution of boundary value problems involving ordinary differential equations
34C25 Periodic solutions to ordinary differential equations

Software:

LOCBIF; AUTO; AUTO-86
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Doedel, E., AUTO: Software for continuation and bifurcation problems in ordinary differential equations, (Report Applied Mathematics (1986), California Institute of Technology: California Institute of Technology Pasadena, USA)
[2] Holodniok, M.; Knedlik, P.; Kubicek, M., Continuation of periodic solutions in parabolic differential equations, (Küpper, T.; Seydel, R.; Troger, H., Bifurcation: Analysis, Algorithms, Applications. Bifurcation: Analysis, Algorithms, Applications, ISNM 79 (1987), Birkhäuser: Birkhäuser Basel), 122-130 · Zbl 0637.65094
[3] Holodniok, M.; Kubicek, M., DERPER — an algorithm for continuation of periodic solutions in ordinary differential equations, J. Comput. Phys., 55, 254-267 (1984) · Zbl 0575.65076
[4] Jarausch, H.; Mackens, W., Numerical treatment of bifurcation problems by adaptive condensation, (Küpper, T.; Mittelman, H. D.; Weber, H., Numerical Methods in Bifurcation Problems. Numerical Methods in Bifurcation Problems, ISNM 70 (1984), Birkhäuser: Birkhäuser Basel), 296-309 · Zbl 0576.65043
[5] Jarausch, H.; Mackens, W., Solving large nonlinear systems of equations by an adaptive condensation process, Numer. Math., 50, 633-653 (1987) · Zbl 0647.65036
[6] Jennings, A.; Stewart, W. J., Simultaneous iteration for partial eigensolutions of real matrices, J. Inst. Math. Applic., 15, 351-361 (1975) · Zbl 0307.65042
[7] Khibnik, A. I.; Kuznetsov, Yu. A.; Levitin, V. V.; Nikolaev, E. N., Continuation techniques and interactive software for bifurcation analysis of ODEs and iterated maps, Physica D, 62, 360-371 (1993) · Zbl 0784.34030
[8] Lust, K.; Roose, D.; Champneys, A.; Spence, A., An adaptive Newton-Picard algorithm with subspace iteration for computing periodic solutions, (Report TW218 (1994), Dept. of Computer Science, K. U. Leuven) · Zbl 1080.65542
[9] Parker, T. S.; Chua, L. O., (Practical Numerical Algorithms for Chaotic Systems (1989), Springer: Springer New-York) · Zbl 0692.58001
[10] Roose, D.; Vandewalle, S., Efficient parallel computation of periodic solutions of partial differential equations, (Seydel, R.; Schneider, F. W.; Küpper, T.; Troger, H., Bifurcations and Chaos: Analysis, Algorithms, Applications. Bifurcations and Chaos: Analysis, Algorithms, Applications, ISNM 97 (1991), Birkhäuser: Birkhäuser Basel), 278-288
[11] Saad, Y., (Numerical Methods for Large Eigenvalue Problems (1992), Manchester University Press: Manchester University Press Manchester) · Zbl 0991.65039
[12] Seydel, R., (From Equilibrium to Chaos. Practical Bifurcation and Stability Analysis (1988), Elsevier: Elsevier New York) · Zbl 0652.34059
[13] Shroff, G. M.; Keller, H. B., Stabilization of unstable procedures: the recursive projection method, SIAM J. Numer. Anal., 30, 4, 1099-1120 (1993) · Zbl 0789.65037
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.