Solution of the eigenvalue problems resulting from global non-parallel flow stability analysis

The eigensolution method, more flexible and efficient than previously used for global non-parallel flow stability analysis, is presented. The differential eigenvalue problem, resulting from linear stability equations of non-parallel Row is discretized with the penalty FEM. The large-dimensional algebraic eigenvalue problem is solved with the use of Subspace Iteration Method. Separation of the eigenvalues interesting for the Bow stability analysis is performed with the inverse Cayley transformation. It is shown that certain barriers, resulting from very large dimensions of the eigenvalue problem and limiting the non-parallel Row stability method can be overcome with this approach. To demonstrate the algorithm, the stability of the flow around circular cylinder is analyzed. (C) 1999 Elsevier Science S.A. All rights reserved.