Numerical computation and continuation of invariant manifolds connecting fixed points

A numerical method for the computation of an invariant manifold that connects two fixed points of a vector field in $\mathbb{R}^n $ is given, extending the results of an earlier paper [Comput. Appl. Math., 26 (1989), pp. 159’170] by the authors. Basically, a boundary value problem on the real line is truncated to a finite interval. The method applies, in particular, to the computation of heteroclinic orbits. The emphasis is on the systematic computation of such orbits by continuation. Using the fact that the linearized operator of our problem is Fredholm in appropriate Banach spaces, the general theory of approximation of nonlinear problems is employed to show that the errors in the approximate solution decay exponentially with the length of the approximating interval. Several applications are considered, including the computation of traveling wave solutions to reaction diffusion problems. Computations were done using the software package AUTO.

[1]  M. Cross,et al.  One dimensional spatial structure near a Hopf bifurcation at finite wavenumber , 1992 .

[2]  Alan C. Lazer,et al.  Large-Amplitude Periodic Oscillations in Suspension Bridges: Some New Connections with Nonlinear Analysis , 1990, SIAM Rev..

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

[4]  J. Descloux,et al.  Approximation of Solution Branches of Non-Linear Equations , 1982 .

[5]  S. Wiggins.  Introduction to Applied Nonlinear Dynamical Systems and Chaos , 1989 .

[6]  S. Hogan.  On the dynamics of rigid-block motion under harmonic forcing , 1989, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[7]  Kenneth J. Palmer,et al.  Exponential dichotomies and transversal homoclinic points , 1984 .

[8]  N. Rashevsky,et al.  Mathematical biology , 1961, Connecticut medicine.

[9]  W. Beyn.  Numerical methods for dynamical systems , 1991 .

[10]  F. V. Vleck,et al.  Stability and Asymptotic Behavior of Differential Equations , 1965 .

[11]  Yu. A. Kuznetsov,et al.  NUMERICAL DETECTION AND CONTINUATION OF CODIMENSION-TWO HOMOCLINIC BIFURCATIONS , 1994 .

[12]  Mark J. Friedman,et al.  Numerical computation of heteroclinic orbits , 1989 .

[13]  Sam D. Howison,et al.  Macroscopic Models for Superconductivity , 1992, SIAM Rev..

[14]  Lee A. Segel,et al.  Modeling Dynamic Phenomena in Molecular and Cellular Biology , 1984 .

[15]  Swinney,et al.  Convective versus absolute instability in flow between counterrotating cylinders. , 1990, Physical review. A, Atomic, molecular, and optical physics.

[16]  David H. Sattinger,et al.  On the stability of waves of nonlinear parabolic systems , 1976 .

[17]  D. P. Stoten,et al.  Robustness of a minimal controller synthesis algorithm , 1990 .

[18]  J. Chomaz,et al.  Bifurcations to local and global modes in spatially developing flows. , 1988, Physical review letters.

[19]  Brian D. Hassard.  Computation of invariant manifolds , 1980 .

[20]  Eugene L. Allgower,et al.  Continuation and path following , 1993, Acta Numerica.

[21]  H. B. Keller,et al.  Boundary Value Problems on Semi-Infinite Intervals and Their Numerical Solution , 1980 .

[22]  R. Seydel.  From Equilibrium to Chaos: Practical Bifurcation and Stability Analysis , 1988 .

[23]  Michael Tinkham,et al.  Introduction to Superconductivity , 1975 .