Aliasing instabilities in the numerical evolution of the Einstein field equations

The Einstein field equations of gravitation are characterized by cross-scale, high-order nonlinear terms, representing a challenge for numerical modeling. In an exact spectral decomposition, high-order nonlinearities correspond to a convolution that numerically might lead to aliasing instabilities....

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:General relativity and gravitation Jg. 53; H. 10
Hauptverfasser: Meringolo, C., Servidio, S.
Format: Journal Article
Sprache:Englisch
Veröffentlicht: New York Springer US 01.10.2021
Springer Nature B.V
Schlagworte:
Algorithms
Aliasing
Astronomy
Astrophysics and Cosmology
Boundary conditions
Cartesian coordinates
Classical and Quantum Gravitation
Decomposition
Differential Geometry
Einstein equations
Gravitation
Gravity
Mathematical and Computational Physics
Nonlinearity
Physics
Physics and Astronomy
Quantum Physics
Relativity Theory
Research Article
Spectra
Theoretical
Wavelengths
Numerical relativity
Black hole dynamics
BSSN decomposition
Gravitational testbeds
ISSN:0001-7701, 1572-9532
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:The Einstein field equations of gravitation are characterized by cross-scale, high-order nonlinear terms, representing a challenge for numerical modeling. In an exact spectral decomposition, high-order nonlinearities correspond to a convolution that numerically might lead to aliasing instabilities. We present a study of this problem, in vacuum conditions, based on the 3 + 1 Baumgarte–Shibata–Shapiro–Nakamura (BSSN) formalism. We inspect the emergence of numerical artifacts, in a variety of conditions, by using the Spectral-FIltered Numerical Gravity codE (SFINGE)—a pseudo-spectral algorithm, based on a classical (Cartesian) Fourier decomposition. By monitoring the highest k - modes of the dynamical fields, we identify the culprits of the aliasing and propose procedures that cure such instabilities, based on the suppression of the aliased wavelengths. This simple algorithm, together with appropriate treatment of the boundary conditions, can be applied to a variety of gravitational problems, including those related to massive objects dynamics.
Bibliographie:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0001-7701
1572-9532
DOI:10.1007/s10714-021-02865-5