Related papers: GRTresna: An open-source code to solve the initial…
Novel applications of Numerical Relativity demand for more flexible algorithms and tools. In this paper, I develop and test a multigrid solver, based on the infrastructure provided by the Einstein Toolkit, for elliptic partial differential…
GRChombo is an open-source code for performing Numerical Relativity time evolutions, built on top of the publicly available Chombo software for the solution of PDEs. Whilst GRChombo uses standard techniques in NR, it focusses on…
Numerical simulations are becoming a more effective tool for conducting detailed investigations into the evolution of our universe. In this article, we show how the framework of numerical relativity can be used for studying cosmological…
GRFolres is an open-source code for performing simulations in modified theories of gravity, based on the publicly available 3+1D numerical relativity code GRChombo. Note: Submitted for review in the Journal of Open Source Software; Comments…
GRDzhadzha is an open-source code for relativistic simulations of matter fields on curved spacetimes that admit an analytic description (e.g. stationary black holes). It is based on the publicly available 3+1D numerical relativity code…
In the extreme violence of merger and mass accretion, compact objects like black holes and neutron stars are thought to launch some of the most luminous outbursts of electromagnetic and gravitational wave energy in the Universe. Modeling…
Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudo-spectral code Elliptica, an…
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…
Numerical relativity became a powerful tool to investigate the dynamics of binary problems with black holes or neutron stars as well as the very structure of General Relativity. Although public numerical relativity codes are available to…
Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and…
The production of numerical relativity waveforms that describe quasicircular binary black hole mergers requires high-quality initial data, and an algorithm to iteratively reduce residual eccentricity. To date, these tools remain closed…
Discontinuous Galerkin (DG) methods for solving elliptic equations are gaining popularity in the computational physics community for their high-order spectral convergence and their potential for parallelization on computing clusters.…
While numerical simulations offer unparalleled precision and robustness in studying complex physical systems, their execution is often hindered by complexity, costliness, and time consumption due to the intricate equations involved. This…
The astrophysics of compact objects, which requires Einstein's theory of general relativity for understanding phenomena such as black holes and neutron stars, is attracting increasing attention. In general relativity, gravity is governed by…
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.…
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…
Many problems in computational science and engineering involve partial differential equations and thus require the numerical solution of large, sparse (non)linear systems of equations. Multigrid is known to be one of the most efficient…
Though the main applications of computer simulations in relativity are to astrophysical systems such as black holes and neutron stars, nonetheless there are important applications of numerical methods to the investigation of general…
Two new methods have been proposed for solving the gravitational constraints without using elliptic solvers by formulating them as either an algebraic-hyperbolic or parabolic-hyperbolic system. Here, we compare these two methods and present…
Numerical relativity is finally coming of age with the development of massively parallel computers. 3D problems, which were completely intractable several years ago due to limited computer power, can now be performed with grid sizes of…