Numerical Analysis and Scientific Computing Seminar

Lyapunov Inverse Iteration for Identifying Hopf Bifurcations in Models of Incompressible Flow

Speaker: Howard Elman, University of Maryland

Location: Warren Weaver Hall 1302

Date: March 2, 2012, 10 a.m.


Identification of instability caused by Hopf bifurcation requires computing the rightmost pair of complex eigenvalues of large eigenvalue problems. A new method replaces this computation with inverse iteration for a certain set of real eigenvalues. We enhance the robustness and efficiency of this approach and demonstrate its utility on examples from fluid dynamics. We explore inexact inner iterations, the impact of solution strategy for a subsidiary Lyapunov system, and effects of eigenvalue distribution on performance.

Joint work with

Karl Meerbergen, Katholieke Universiteit Leuven, Belgium Alastair Spence, University of Bath, United Kingdom Minghao Wu, University of Maryland