Related papers: ExaGRyPE: Numerical General Relativity Solvers Bas…
In this work numerical methods for solving Einstein's equations are developed and applied to the study of inhomogeneous cosmological models. A two-dimensional computer code is described which implements two advanced numerical methods:…
The Einstein and Maxwell equations are both systems of hyperbolic equations which need to satisfy a set of elliptic constraints throughout evolution. However, while electrodynamics (EM) and magnetohydrodynamics (MHD) have benefited from a…
Many problems in geophysical and atmospheric modelling require the fast solution of elliptic partial differential equations (PDEs) in "flat" three dimensional geometries. In particular, an anisotropic elliptic PDE for the pressure…
Numerical solvers of Partial Differential Equations (PDEs) are of fundamental significance to science and engineering. To date, the historical reliance on legacy techniques has circumscribed possible integration of big data knowledge and…
We present GRaM-X (General Relativistic accelerated Magnetohydrodynamics on AMReX), a new GPU-accelerated dynamical-spacetime general relativistic magnetohydrodynamics (GRMHD) code which extends the GRMHD capability of Einstein Toolkit to…
Numerical relativity is an essential tool for solving Einstein's equations of general relativity for dynamical systems characterized by high velocities and strong gravitational fields. The implementation of new algorithms that can solve…
It is well known that multigrid methods are optimally efficient for solution of elliptic equations (O(N)), which means that effort is proportional to the number of points at which the solution is evaluated). Thus this is an ideal method to…
The demand for substantial increases in the spatial resolution of global weather- and climate- prediction models makes it necessary to use numerically efficient and highly scalable algorithms to solve the equations of large scale…
The computational complexity of classical numerical methods for solving Partial Differential Equations (PDE) scales significantly as the resolution increases. As an important example, climate predictions require fine spatio-temporal…
We present a novel spectral solver for general relativistic magnetohydrodynamics on dynamical spacetimes. By combining a high order discontinuous spectral method on mapped Chebyshev Fourier grids, our scheme attains exponential convergence.…
We use rigorous techniques from numerical analysis of hyperbolic equations in bounded domains to construct stable finite-difference schemes for Numerical Relativity, in particular for their use in black hole excision. As an application, we…
This paper presents an overview of the functionalities and applications of Exasim, an open-source code for generating high-order discontinuous Galerkin codes to numerically solve parametrized partial differential equations (PDEs). The…
Apparent horizon (AH) finders are essential for characterizing black holes and excising their interiors in numerical relativity (NR) simulations. However, open-source AH finders to date are tightly coupled to individual NR codes. We…
A new numerical scheme to solve the Einstein field equations based upon the generalized harmonic decomposition of the Ricci tensor is introduced. The source functions driving the wave equations that define generalized harmonic coordinates…
Detailed observations of phenomena involving black holes, be it via gravitational waves or more traditional electromagnetic means, can probe the strong field regime of the gravitational interaction. The prediction of features in such…
Solving partial differential equations (PDEs) is an important research means in the fields of physics, biology, and chemistry. As an approximate alternative to numerical methods, PINN has received extensive attention and played an important…
A new algorithm for solving the general relativistic MHD equations is described in this paper. We design our scheme to incorporate black hole excision with smooth boundaries, and to simplify solving the combined Einstein and MHD equations…
We describe a numerical code that solves Einstein's equations for a Schwarzschild black hole in spherical symmetry, using a hyperbolic formulation introduced by Choquet-Bruhat and York. This is the first time this formulation has been used…
We construct new static, spherically symmetric non-extremal black hole solutions of four-dimensional ${\cal N}=2$ supergravity, using a systematic technique based on dimensional reduction over time (the c-map) and the real formulation of…
Associated to every stationary extremal black hole is a unique near-horizon geometry, itself a solution of the field equations. These latter spacetimes are more tractable to analyze and most importantly, retain properties of the original…