Related papers: GRTresna: An open-source code to solve the initial…
Einstein's field equation of General Relativity (GR) has been known for over 100 years, yet it remains challenging to solve analytically in strongly relativistic regimes, particularly where there is a lack of a priori symmetry. Numerical…
We explore whether a new method to solve the constraints of Einstein's equations, which does not involve elliptic equations, can be applied to provide initial data for black holes. We show that this method can be successfully applied to a…
We solve the elliptic equations associated with the Hamiltonian and momentum constraints, corresponding to a system composed of two black holes with arbitrary linear and angular momentum. These new solutions are based on a Kerr-Schild…
We present a new scheme for constructing initial data for the Einstein field equations using the conformal thin-sandwich formulation that does not assume conformal flatness or approximate Killing vectors. This includes a method for…
Spherically symmetric (1D) black-hole spacetimes are considered as a test for numerical relativity. A finite difference code, based in the hyperbolic structure of Einstein's equations with the harmonic slicing condition is presented.…
This work consists of two distinct parts. In the first part we present a new method for solving the initial value problem of general relativity. Given any spatial metric with a surface orthogonal Killing field and two freely specified…
We present a new approach for setting initial Cauchy data for multiple black hole spacetimes. The method is based upon adopting an initially Kerr-Schild form of the metric. In the case of non-spinning holes, the constraint equations take a…
We present an overview of recent developments in the numerical solution of Horndeski gravity theories, which are the class of all scalar-tensor theories of gravity that have second order equations of motion. We review several methods that…
We extend a recently developed numerical code to obtain stationary, axisymmetric solutions that describe rotating black hole spacetimes in a wide class of modified theories of gravity. The code utilizes a relaxed Newton-Raphson method to…
Elliptic partial differential equations must be solved numerically for many problems in numerical relativity, such as initial data for every simulation of merging black holes and neutron stars. Existing elliptic solvers can take multiple…
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…
We describe the special relativistic extension of the CRONOS code, which has been used for studies of gamma-ray binaries in recent years. The code was designed to be easily adaptable, allowing the user to easily change existing…
There are many complementary approaches to the construction of solutions to the field equations of general relativity. Among these, numerical approximation offers the only possibility to compute a variety of dynamical spacetimes, and so has…
In order to generate initial data for nonlinear relativistic simulations, one needs to solve the Einstein constraints, which can be cast into a coupled set of nonlinear elliptic equations. Here we present an approach for solving these…
Given the compact binary evolution problem of numerical relativity, in the finite-difference, block-based, adaptive mesh refinement context, choices must be made on how evolved fields are to be discretized. In GR-Athena++, the space-time…
We present a new two-dimensional numerical code called Nada designed to solve the full Einstein equations coupled to the general relativistic hydrodynamics equations. The code is mainly intended for studies of self-gravitating accretion…
Numerical relativity codes now being developed will evolve initial data representing colliding black holes at a relatively late stage in the collision. The choice of initial data used for code development has been made on the basis of…
We report on a new open-source, user-friendly numerical relativity code package called SENR/NRPy+. Our code extends previous implementations of the BSSN reference-metric formulation to a much broader class of curvilinear coordinate systems,…
Next-generation gravitational wave detectors such as Cosmic Explorer, the Einstein Telescope, and LISA, demand highly accurate and extensive gravitational wave (GW) catalogs to faithfully extract physical parameters from observed signals.…
The equations of general relativistic magnetohydrodynamics (GRMHD) have become the standard mathematical framework for modeling high-energy plasmas in curved spacetimes. However, the fragility of the primitive variable reconstruction…