Related papers: The computation of seismic normal modes with rotat…
Three new Arnoldi-type methods are presented to accelerate the modal analysis and critical speed analysis of the damped rotor dynamics finite element (FE) model. They are the linearized quadratic eigenvalue problem (QEP) Arnoldi method, the…
Non-Hermitian generalized eigenvalue problems (GEPs) play a significant role in many practical applications, such as mechanical engineering. Based on the generalized Schur decomposition, we propose a variational quantum algorithm for…
Analyzing the stability of quasinormal modes (QNM) is essential for understanding black hole dynamics, particularly in the context of gravitational wave emissions and black hole spectroscopy. In this study, we employ the hyperboloidal…
We extend the conforming virtual element method to the numerical resolution of eigenvalue problems with potential terms on a polytopal mesh. An important application is that of the Schrodinger equation with a pseudopotential term. This…
In this paper the relationship between the problem of constructing the ground state energy for the quantum quartic oscillator and the problem of computing mean eigenvalue of large positively definite random hermitean matrices is…
We solve Poisson's equation by combining a spectral-element method with a mapped infinite-element method. We focus on problems in geostatics and geodynamics, where Earth's gravitational field is determined by Poisson's equation inside the…
This paper proposes a finite element method for solving the periodic steady-state problem for the scalar-valued and vector-valued Poisson equations, a simple reduction model of the Maxwell equations under the Coulomb gauge. Introducing a…
Spectral methods are now common in the solution of ordinary differential eigenvalue problems in a wide variety of fields, such as in the computation of black hole quasinormal modes. Most of these spectral codes are based on standard…
Estimating the eigenvalues of non-normal matrices is a foundational problem with far-reaching implications, from modeling non-Hermitian quantum systems to analyzing complex fluid dynamics. Yet, this task remains beyond the reach of standard…
In this work, we describe, analyze, and implement a pseudospectral quadrature method for a global computer modeling of the incompressible surface Navier-Stokes equations on the rotating unit sphere. Our spectrally accurate numerical error…
The rectangular multiparameter eigenvalue problem (RMEP) involves rectangular coefficient matrices (usually with more rows than columns) and may potentially have no solution in its original form. A minimal perturbation framework is proposed…
Rapidly rotating spherical kinematic dynamos are computed using the combination of a quasi geostrophic (QG) model for the velocity field and a classical spectral 3D code for the magnetic field. On one hand, the QG flow is computed in the…
Many of exoplanetary systems consist of more than one planet and the study of planetary orbits with respect to their long-term stability is very interesting. Furthermore, many exoplanets seem to be locked in a mean-motion resonance (MMR),…
This paper aims to illustrate the applications of resonant Hamiltonian normal forms to some problems of galactic dynamics. We detail the construction of the 1:1 resonant normal form corresponding to a wide class of potentials with…
The eigenvalue density of a quantum-mechanical system exhibits oscillations, determined by the closed orbits of the corresponding classical system; this relationship is simple and strong for waves in billiards or on manifolds, but becomes…
This paper deals with the spectral element modeling of seismic wave propagation at the global scale. Two aspects relevant to low-frequency studies are particularly emphasized. First, the method is generalized beyond the Cowling…
High-precision calculations of hadron spectroscopy are a crucial task for Lattice QCD. State-of-the-art techniques are needed to disentangle the contributions from different energy states, such as solving the generalized eigenvalue problem…
Elastic wave manipulation using large arrays of resonators is driving the need for advanced simulation and optimization methods. To address this we introduce and explore a robust framework for wave control: Quasi-normal modes (QNMs).…
We present a new algorithm that computes eigenvalues and eigenvectors of a Hermitian positive definite matrix while solving a linear system of equations with Conjugate Gradient (CG). Traditionally, all the CG iteration vectors could be…
In this paper, we propose a unified approach for solving structure-preserving eigenvalue embedding problem (SEEP) for quadratic regular matrix polynomials with symmetry structures. First, we determine perturbations of a quadratic matrix…