Lyapunov inverse iteration for identifying Hopf bifurcations in models of incompressible flow

Howard C. Elman, Karl Meerbergen, Alastair Spence, Minghao Wu

Research output: Contribution to journalArticle

138 Downloads (Pure)

Abstract

The identification of instability in large-scale dynamical systems caused by Hopf bifurcation is difficult because of the problem of identifying the rightmost pair of complex eigenvalues of large sparse generalized eigenvalue problems. A new method developed in [K. Meerbergen and A. Spence, SIAM J. Matrix Anal. Appl., 31 (2010), pp. 1982--1999] avoids this computation, instead performing an inverse iteration for a certain set of real eigenvalues that requires the solution of a large-scale Lyapunov equation at each iteration. In this study, we refine the Lyapunov inverse iteration method to make it more robust and efficient, and we examine its performance on challenging test problems arising from fluid dynamics. Various implementation issues are discussed, including the use of inexact inner iterations and the impact of the choice of iterative solution for the Lyapunov equations, and the effect of eigenvalue distribution on performance. Numerical experiments demonstrate the robustness of the algorithm.


Read More: http://epubs.siam.org/doi/abs/10.1137/110827600
Original languageEnglish
Pages (from-to)A1584-A1606
Number of pages23
JournalSIAM Journal on Scientific Computing
Volume34
Issue number3
Early online date11 Jun 2012
DOIs
Publication statusPublished - 2012

Fingerprint Dive into the research topics of 'Lyapunov inverse iteration for identifying Hopf bifurcations in models of incompressible flow'. Together they form a unique fingerprint.

  • Cite this