Related papers: Toward computing gravitational initial data withou…
In this paper we continue earlier investigations of evolutionary formulations of the Einstein vacuum constraint equations originally introduced by R\'{a}cz. Motivated by the strong evidence from these works that the resulting vacuum initial…
We investigate the initial-boundary value problem for linearized gravitational theory in harmonic coordinates. Rigorous techniques for hyperbolic systems are applied to establish well-posedness for various reductions of the system into a…
We consider the design and analysis of numerical methods for approximating positive solutions to nonlinear geometric elliptic partial differential equations containing critical exponents. This class of problems includes the Yamabe problem…
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…
A new development of the ``monodromy transform'' method for analysis of hyperbolic as well as elliptic integrable reductions of Einstein equations is presented. Compatibility conditions for some alternative representations of the…
We are concerned with hyperbolic systems of order-one linear PDEs originated on non-characteristic manifolds. We put forward a simple but effective method of transforming such initial conditions to standard initial conditions (i.e. when the…
We consider a series of distorted black hole initial data sets, and develop techniques to evolve them using the linearized equations of motion for the gravitational wave perturbations on a Schwarzschild background. We apply this to 2D and…
The purpose of this work is to construct asymptotically flat, time symmetric initial data with an apparent horizon of prescribed intrinsic geometry. To do this, we use the parabolic partial differential equation for prescribing scalar…
When numerically solving Einstein's equations for binary black holes (BBH), we must find initial data on a three-dimensional spatial slice by solving constraint equations. The construction of initial data is a multi-step process, in which…
Vacuum solutions to the Einstein equations can be viewed as the interplay between the geometry and the gravitational wave energy content. The constraints on initial data reflect this interaction. We assume we are looking at cosmological…
We propose and explore a "stationary 1+log" slicing condition for the construction of solutions to Einstein's constraint equations. For stationary spacetimes, these initial data will give a stationary foliation when evolved with "moving…
By applying the method of moving frames modelling one and two dimensional local anisotropies we construct new solutions of Einstein equations on pseudo-Riemannian spacetimes. The first class of solutions describes non-trivial deformations…
With the aim of deriving symmetric hyperbolic free-evolution systems for GR that possess Hamiltonian structure and allow for the popular puncture gauge condition we analyze the hyperbolicity of Hamiltonian systems. We develop helpful tools…
A method for studying the causal structure of space-time evolution systems is presented. This method, based on a generalization of the well known Riemann problem, provides intrinsic results which can be interpreted from the geometrical…
We investigate the Helmholtz equation with suitable boundary conditions and uncertainties in the wavenumber. Thus the wavenumber is modeled as a random variable or a random field. We discretize the Helmholtz equation using finite…
It is possible to solve the Einstein constraint equations as an evolutionary rather than an elliptic system. Here we consider the Gauss constraint in electrodynamics as a toy model for thist. We use a combination of the evolutionary method…
We introduce NRPyElliptic, an elliptic solver for numerical relativity (NR) built within the NRPy+ framework. As its first application, NRPyElliptic sets up conformally flat, binary black hole (BBH) puncture initial data (ID) on a single…
We present a new pseudospectral code, bamps, for numerical relativity written with the evolution of collapsing gravitational waves in mind. We employ the first order generalized harmonic gauge formulation. The relevant theory is reviewed…
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…
We address the generation of initial conditions (ICs) for GRAMSES, a code for nonlinear general relativistic (GR) $N$-body cosmological simulations recently introduced in Ref. [1]. GRAMSES adopts a constant mean curvature slicing with a…