Convergence and round-off errors in a two-dimensional eigenvalue problem using spectral methods and Arnoldi-Chebyshev algorithm
From MaRDI portal
Publication:2372937
Abstract: An efficient way of solving 2D stability problems in fluid mechanics is to use, after discretization of the equations that cast the problem in the form of a generalized eigenvalue problem, the incomplete Arnoldi-Chebyshev method. This method preserves the banded structure sparsity of matrices of the algebraic eigenvalue problem and thus decreases memory use and CPU-time consumption. The errors that affect computed eigenvalues and eigenvectors are due to the truncation in the discretization and to finite precision in the computation of the discretized problem. In this paper we analyze those two errors and the interplay between them. We use as a test case the two-dimensional eigenvalue problem yielded by the computation of inertial modes in a spherical shell. This problem contains many difficulties that make it a very good test case. It turns out that that single modes (especially most-damped modes i.e. with high spatial frequency) can be very sensitive to round-off errors, even when apparently good spectral convergence is achieved. The influence of round-off errors is analyzed using the spectral portrait technique and by comparison of double precision and extended precision computations. Through the analysis we give practical recipes to control the truncation and round-off errors on eigenvalues and eigenvectors.
Recommendations
- Removal of spurious modes encountered in solving stability problems by spectral methods
- Numerical solution of eigenvalue problems using spectral techniques
- Spectral Galerkin discretization for hydrodynamic stability problems
- Efficient solution of a generalized eigenvalue problem arising in a thermoconvective instability
- Hydrodynamic stability, the Chebyshev tau method and spurious eigenvalues
Cites work
- scientific article; zbMATH DE number 432500 (Why is no real title available?)
- scientific article; zbMATH DE number 1260326 (Why is no real title available?)
- scientific article; zbMATH DE number 770270 (Why is no real title available?)
- scientific article; zbMATH DE number 1405268 (Why is no real title available?)
- Highly nonnormal eigenproblems in the aeronautical industry
- Hydrodynamic stability without eigenvalues
- Inertial waves in a rotating spherical shell
- Inertial waves in a rotating spherical shell: Attractors and asymptotic spectrum
- Jacobi--Davidson Style QR and QZ Algorithms for the Reduction of Matrix Pencils
- Large-scale computation of pseudospectra using ARPACK and eigs
- Lectures on Finite Precision Computations
- Linear theory of rotating fluids using spherical harmonics part I: Steady flows
- Spectral methods in general relativistic astrophysics
Cited in
(11)- Fingering convection in a spherical shell
- The impact of eigenvalue locality on the convergence behavior of the PSD method for two-cyclic matrices
- Viscous dissipation by tidally forced inertial modes in a rotating spherical shell
- A note on high-precision approximation of asymptotically decaying solution and orthogonal decomposition
- Axisymmetric inertial modes in a spherical shell at low Ekman numbers
- Gauss-Lobatto to Bernstein polynomials transformation
- Efficient solution of a generalized eigenvalue problem arising in a thermoconvective instability
- Inertial waves in a differentially rotating spherical shell
- Gravito-inertial waves in a differentially rotating spherical shell
- High-precision numerical simulations on a CUDA GPU: Kerr black hole tails
- Compressible fluid modes in rigid ellipsoids: towards modal acoustic velocimetry
This page was built for publication: Convergence and round-off errors in a two-dimensional eigenvalue problem using spectral methods and Arnoldi-Chebyshev algorithm
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2372937)
