arrow
Volume 3, Issue 2
Efficient Solution of a Generalized Eigenvalue Problem Arising in a Thermoconvective Instability

M. C. Navarro, H. Herrero, A. M. Mancho & A. Wathen

Commun. Comput. Phys., 3 (2008), pp. 308-329.

Published online: 2008-03

Export citation
  • Abstract

The aim of this paper is to develop an efficient numerical method to compute the eigenvalues of the stability analysis of a problem describing the motion of a fluid within a cylindrical container heated non-homogeneously from below. An axisymmetric stationary motion settles in, at certain values of the external parameters appearing in the set of partial differential equations modeling the problem. This basic solution is computed by discretizing the equations with a Chebyshev collocation method. Its linear stability is formulated with a generalized eigenvalue problem. The numerical approach (generalized Arnoldi method) uses the idea of preconditioning the eigenvalue problem with a modified Cayley transformation before applying the Arnoldi method. Previous works have dealt with transformations requiring regularity to one of the submatrices. In this article we extend those results to the case in which that submatrix is singular. This method allows a fast computation of the critical eigenvalues which determine whether the steady flow is stable or unstable. The algorithm based on this method is compared to the QZ method and is found to be computationally more efficient. The reliability of the computed eigenvalues in terms of stability is confirmed via pseudospectra calculations.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-3-308, author = {}, title = {Efficient Solution of a Generalized Eigenvalue Problem Arising in a Thermoconvective Instability}, journal = {Communications in Computational Physics}, year = {2008}, volume = {3}, number = {2}, pages = {308--329}, abstract = {

The aim of this paper is to develop an efficient numerical method to compute the eigenvalues of the stability analysis of a problem describing the motion of a fluid within a cylindrical container heated non-homogeneously from below. An axisymmetric stationary motion settles in, at certain values of the external parameters appearing in the set of partial differential equations modeling the problem. This basic solution is computed by discretizing the equations with a Chebyshev collocation method. Its linear stability is formulated with a generalized eigenvalue problem. The numerical approach (generalized Arnoldi method) uses the idea of preconditioning the eigenvalue problem with a modified Cayley transformation before applying the Arnoldi method. Previous works have dealt with transformations requiring regularity to one of the submatrices. In this article we extend those results to the case in which that submatrix is singular. This method allows a fast computation of the critical eigenvalues which determine whether the steady flow is stable or unstable. The algorithm based on this method is compared to the QZ method and is found to be computationally more efficient. The reliability of the computed eigenvalues in terms of stability is confirmed via pseudospectra calculations.

}, issn = {1991-7120}, doi = {https://doi.org/}, url = {http://global-sci.org/intro/article_detail/cicp/7855.html} }
TY - JOUR T1 - Efficient Solution of a Generalized Eigenvalue Problem Arising in a Thermoconvective Instability JO - Communications in Computational Physics VL - 2 SP - 308 EP - 329 PY - 2008 DA - 2008/03 SN - 3 DO - http://doi.org/ UR - https://global-sci.org/intro/article_detail/cicp/7855.html KW - AB -

The aim of this paper is to develop an efficient numerical method to compute the eigenvalues of the stability analysis of a problem describing the motion of a fluid within a cylindrical container heated non-homogeneously from below. An axisymmetric stationary motion settles in, at certain values of the external parameters appearing in the set of partial differential equations modeling the problem. This basic solution is computed by discretizing the equations with a Chebyshev collocation method. Its linear stability is formulated with a generalized eigenvalue problem. The numerical approach (generalized Arnoldi method) uses the idea of preconditioning the eigenvalue problem with a modified Cayley transformation before applying the Arnoldi method. Previous works have dealt with transformations requiring regularity to one of the submatrices. In this article we extend those results to the case in which that submatrix is singular. This method allows a fast computation of the critical eigenvalues which determine whether the steady flow is stable or unstable. The algorithm based on this method is compared to the QZ method and is found to be computationally more efficient. The reliability of the computed eigenvalues in terms of stability is confirmed via pseudospectra calculations.

M. C. Navarro, H. Herrero, A. M. Mancho & A. Wathen. (2020). Efficient Solution of a Generalized Eigenvalue Problem Arising in a Thermoconvective Instability. Communications in Computational Physics. 3 (2). 308-329. doi:
Copy to clipboard
The citation has been copied to your clipboard