Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published August 1993 | Published
Journal Article Open

Stabilization of Unstable Procedures: The Recursive Projection Method


Fixed-point iterative procedures for solving nonlinear parameter dependent problems can converge for some interval of parameter values and diverge as the parameter changes. The Recursive Projection Method (RPM), which stabilizes such procedures by computing a projection onto the unstable subspace is presented. On this subspace a Newton or special Newton iteration is performed, and the fixed-point iteration is used on the complement. As continuation in the parameter proceeds, the projection is efficiently updated, possibly increasing or decreasing the dimension of the unstable subspace. The method is extremely effective when the dimension of the unstable subspace is small compared to the dimension of the system. Convergence proofs are given and pseudo-arclength continuation on the unstable subspace is introduced to allow continuation past folds. Examples are presented for an important application of the RPM in which a "black-box" time integration scheme is stabilized, enabling it to compute unstable steady states. The RPM can also be used to accelerate iterative procedures when slow convergence is due to a few slowly decaying modes.

Additional Information

© 1993 Society for Industrial and Applied Mathematics. Received by the editors June 3, 1991; accepted for publication (in revised form) September 11, 1992. This work was supported in part by the National Science Foundation Cooperative Agreement CCR-9120008 and the Department of Energy Project Agreement No. DE-FG03-89ER25073. The government has certain rights on this material.

Attached Files

Published - SHROsiamjna93.pdf


Files (2.4 MB)
Name Size Download all
2.4 MB Preview Download

Additional details

August 22, 2023
October 24, 2023