TY - JOUR
T1 - Multidimensional manifold continuation for adaptive boundary-value problems
AU - Dankowicz, Harry
AU - Wang, Yuqing
AU - Schilder, Frank
AU - Henderson, Michael E.
N1 - Publisher Copyright:
Copyright © 2020 by ASME.
PY - 2020/5
Y1 - 2020/5
N2 - Parameter continuation of finitely parameterized, approximate solutions to integro-differential boundary-value problems typically involves regular adaptive updates to the number and meaning of the unknowns and/or the associated constraints. Different continuation steps produce solutions with different discretizations or to formally different sets of equations. Existing general-purpose, multidimensional continuation algorithms fail to account for such differences without significant additional coding and are therefore prone to redundant coverage of the set of solutions. We describe a new algorithm, implemented in the software package COCO, which overcomes this problem by characterizing the solution set in an invariant, finite dimensional, projected geometry rather than in the space of unknowns corresponding to any particular discretization. It is in this geometry that distances between solutions and angles between tangent spaces are quantified and used to construct possible directions of outward expansion. A pointwise lift identifies such directions in the projected geometry with directions of continuation in the full set of unknowns, used by a nonlinear predictor-corrector algorithm to expand into uncharted parts of the solution set. Several benchmark problems from the analysis of periodic orbits in autonomous dynamical systems are used to illustrate the theory.
AB - Parameter continuation of finitely parameterized, approximate solutions to integro-differential boundary-value problems typically involves regular adaptive updates to the number and meaning of the unknowns and/or the associated constraints. Different continuation steps produce solutions with different discretizations or to formally different sets of equations. Existing general-purpose, multidimensional continuation algorithms fail to account for such differences without significant additional coding and are therefore prone to redundant coverage of the set of solutions. We describe a new algorithm, implemented in the software package COCO, which overcomes this problem by characterizing the solution set in an invariant, finite dimensional, projected geometry rather than in the space of unknowns corresponding to any particular discretization. It is in this geometry that distances between solutions and angles between tangent spaces are quantified and used to construct possible directions of outward expansion. A pointwise lift identifies such directions in the projected geometry with directions of continuation in the full set of unknowns, used by a nonlinear predictor-corrector algorithm to expand into uncharted parts of the solution set. Several benchmark problems from the analysis of periodic orbits in autonomous dynamical systems are used to illustrate the theory.
UR - http://www.scopus.com/inward/record.url?scp=85088378727&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85088378727&partnerID=8YFLogxK
U2 - 10.1115/1.4046498
DO - 10.1115/1.4046498
M3 - Article
AN - SCOPUS:85088378727
SN - 1555-1415
VL - 15
JO - Journal of Computational and Nonlinear Dynamics
JF - Journal of Computational and Nonlinear Dynamics
IS - 5
M1 - 051002-1
ER -