Related papers: The computation of seismic normal modes with rotat…
Properties of inertial modes of Jupiter are investigated for an n=1 polytropic description of the planet interior. We use the anelastic approximation to overcome the usual handicap of a severe spherical harmonics truncation. A powerful…
The quasinormal modes (QNMs) of a rotating quantum corrected black hole (RQCBH) are studied by employing the hyperboloidal framework for the scalar perturbation. This framework is used to cast the QNMs spectra problem into a two-dimensional…
Quantum algorithms for estimating the eigenvalues of matrices, including the phase estimation algorithm, serve as core subroutines in a wide range of quantum algorithms, including those in quantum chemistry and quantum machine learning. The…
We analyse and compare several algorithms to compute numerically periodic solutions of high-dimensional dynamical systems and investigate their Floquet stability without building the monodromy matrix. The solution and its perturbation are…
The generalized pseudospectral method is employed for the accurate calculation of eigenvalues, densities and expectation values for the spiked harmonic oscillators. This allows \emph{nonuniform} and \emph{optimal} spatial discretization of…
In this article, we investigate inertial modes of rigidly rotating neutron stars, i.e. modes for which the Coriolis force is dominant. This is done using the assumption of a fixed spacetime (Cowling approximation). We present frequencies…
We demonstrate an application of the spectral method as a numerical approximation for solving Hyperbolic PDEs. In this method a finite basis is used for approximating the solutions. In particular, we demonstrate a set of such solutions for…
We present a semi-analytical approach to compute quasi-guided elastic wave modes in horizontally layered structures radiating into unbounded fluid or solid media. This problem is of relevance, e.g., for the simulation of guided ultrasound…
Quartic eigenvalue problem $(\lambda^4 A + \lambda^3 B + \lambda^2C + \lambda D + E)x = \mathbf{0}$ naturally arises e.g. when solving the Orr-Sommerfeld equation in the analysis of the stability of the {Poiseuille} flow, in theoretical…
For large-scale eigenvalue problems requiring many mutually orthogonal eigenvectors, traditional numerical methods suffer substantial computational and communication costs with limited parallel scalability, primarily due to explicit…
A new solution strategy for quadratic eigenvalue problems, and the derivatives of the eigenvalues, is proposed, by combining the generalized reduction method with dual numbers. To demonstrate the method, we use the quadratic eigenvalue…
In this letter, a matrix method is employed to study the scalar quasinormal modes of Kerr as well as Kerr-Sen black holes. Discretization is applied to transfer the scalar perturbation equation into a matrix form eigenvalue problem, where…
A solid object's geometry, density, and elastic moduli completely determine its spectrum of normal modes. Solving the inverse problem - determining a material's elastic moduli given a set of resonance frequencies and sample geometry -…
We compute the inertial modes of a freely rotating two-layer planetary model with an ellipsoidal inviscid fluid core and a perfectly rigid mantle. We present a method to derive analytical formulae for the frequencies of the Free Core…
We describe the Reversibility Error Method (REM) and its applications to planetary dynamics. REM is based on the time-reversibility analysis of the phase-space trajectories of conservative Hamiltonian systems. The round-off errors break the…
In a previous paper (J. Comp. Phys. 230 (2011), 3668--3694), the authors proposed a new practical method for computing expected values of functionals of solutions for certain classes of elliptic partial differential equations with random…
A quasimodal expansion method (QMEM) is developed to model and understand the scattering properties of arbitrary shaped two-dimensional (2-D) open structures. In contrast with the bounded case which have only discrete spectrum (real in the…
Self-sustained convective dynamos in planetary systems operate in an asymptotic regime of rapid rotation, where a balance is thought to hold between the Coriolis, pressure, buoyancy and Lorentz forces (the MAC balance). Classical numerical…
We consider the bi-Laplacian eigenvalue problem for the modes of vibration of a thin elastic plate with a discrete set of clamped points. A high-order boundary integral equation method is developed for efficient numerical determination of…
Concerning the anelastic nature of Earth, the quasi-static extended Burgers model (abbreviated by q-EBM), an integro-differential system, is used to study the free oscillation of Earth (abbreviated by FOE). In this paper, we first provide a…