Convergence and round-off errors in a two-dimensional eigenvalue problem using spectral methods and Arnoldi-Chebyshev algorithm

DOI10.1016/J.CAM.2006.05.009zbMATH Open1246.76119arXivphysics/0604219OpenAlexW2117636168MaRDI QIDQ2372937FDOQ2372937

Authors: M. Rieutord, T. Braconnier, Valérie Frayssé, L. Valdettaro 

Publication date: 17 July 2007

Published in: Journal of Computational and Applied Mathematics (Search for Journal in Brave)

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.

Full work available at URL: https://arxiv.org/abs/physics/0604219

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

zbMATH Keywords

eigenvalue problemrigidly rotating fluidsperical shell

Mathematics Subject Classification ID

Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Spectral methods applied to problems in fluid mechanics (76M22) General theory of rotating fluids (76U05)

Cites Work

• Large-scale computation of pseudospectra using ARPACK and eigs
• Hydrodynamic stability without eigenvalues
• Spectral methods in general relativistic astrophysics
• Title not available (Why is that?)
• Jacobi--Davidson Style QR and QZ Algorithms for the Reduction of Matrix Pencils
• Title not available (Why is that?)
• Inertial waves in a rotating spherical shell
• Title not available (Why is that?)
• Lectures on Finite Precision Computations
• Inertial waves in a rotating spherical shell: Attractors and asymptotic spectrum
• Highly nonnormal eigenproblems in the aeronautical industry
• Linear theory of rotating fluids using spherical harmonics part I: Steady flows
• Title not available (Why is that?)

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

Uses Software

• eigs
• JDQR
• ARPACK
• Eigtool

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)