Arbitrary precision computation of hydrodynamic stability eigenvalues
Fuente:
arXiv
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866912354773499904 |
|---|---|
| author | Dondl, Patrick Striet, Ludwig Straughan, Brian |
| author_facet | Dondl, Patrick Striet, Ludwig Straughan, Brian |
| contents | We show that by using higher order precision arithmetic, i.e., using floating point types with more significant bits than standard double precision numbers, one may accurately compute eigenvalues for non-normal matrices arising in hydrodynamic stability problems. The basic principle is illustrated by a classical example of two $7\times 7$ matrices for which it is well known that eigenvalue computations fail when using standard double precision arithmetic. We then present an implementation of the Chebyshev tau-QZ method allowing the use of a large number of Chebyshev polynomials together with arbitrary precision arithmetic. This is used to compute the behavior of the spectra for Couette and Poiseuille flow at high Reynolds number. An experimental convergence analysis finally makes it evident that high order precision is required to obtain accurate results. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2504_21511 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Arbitrary precision computation of hydrodynamic stability eigenvalues Dondl, Patrick Striet, Ludwig Straughan, Brian Numerical Analysis 35Q35, 76E06, 76E05, 76D99 We show that by using higher order precision arithmetic, i.e., using floating point types with more significant bits than standard double precision numbers, one may accurately compute eigenvalues for non-normal matrices arising in hydrodynamic stability problems. The basic principle is illustrated by a classical example of two $7\times 7$ matrices for which it is well known that eigenvalue computations fail when using standard double precision arithmetic. We then present an implementation of the Chebyshev tau-QZ method allowing the use of a large number of Chebyshev polynomials together with arbitrary precision arithmetic. This is used to compute the behavior of the spectra for Couette and Poiseuille flow at high Reynolds number. An experimental convergence analysis finally makes it evident that high order precision is required to obtain accurate results. |
| title | Arbitrary precision computation of hydrodynamic stability eigenvalues |
| topic | Numerical Analysis 35Q35, 76E06, 76E05, 76D99 |
| url | https://arxiv.org/abs/2504.21511 |