Related papers: The BSSN formulation is a partially constrained ev…
Full relativistic simulations in three dimensions invariably develop runaway modes that grow exponentially and are accompanied by violations of the Hamiltonian and momentum constraints. Recently, we introduced a numerical method…
I present a new, simple method to dynamically control the growth of the discretized constraints during a free evolution of Einstein's equations. During an evolution, any given family of formulations is adjusted off the constraints surface…
We analyze the effect of pressure on the evolution of perturbations of an Einstein-de Sitter Universe in the matter dominated epoch assuming an ideal gas equation of state. For the sake of simplicity the temperature is considered uniform.…
We present a novel approach, based entirely on the gravitational potential, for studying the evolution of non-linear cosmological matter perturbations. Starting from the perturbed Einstein equations, we integrate out the non-relativistic…
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…
Numerical relativity is finally approaching a state where the evolution of rather general (3+1)-dimensional data sets can be computed in order to solve the Einstein equations. After a general introduction, three topics of current interest…
We present a set of well-posed constraint-preserving boundary conditions for a first-order in time, second-order in space, harmonic formulation of the Einstein equations. The boundary conditions are tested using robust stability, linear and…
We compare numerical evolutions performed with the BSSNOK formulation and a conformal decomposition of a Z4-like formulation of General Relativity. The important difference between the two formulations is that the Z4 formulation has a…
The evolution equations of Einstein's theory and of Maxwell's theory---the latter used as a simple model to illustrate the former--- are written in gauge covariant first order symmetric hyperbolic form with only physically natural…
We present a new formulation of the Einstein equations that casts them in an explicitly first order, flux-conservative, hyperbolic form. We show that this now can be done for a wide class of time slicing conditions, including maximal…
The puncture method for dealing with black holes in the numerical simulation of vacuum spacetimes is remarkably successful when combined with the BSSN formulation of the Einstein equations. We examine a generalized class of formulations…
Using techniques from singular perturbation theory, we explicitly calculate the cosmological evolution in a class of modified gravity models. By considering the (m)CDTT model, which aims to explain the current acceleration of the universe…
We review some old and new results about strict and non strict hyperbolic formulations of the Einstein equations.
We investigate how the accuracy and stability of numerical relativity simulations of 1D colliding plane waves depends on choices of equation formulations, gauge conditions, boundary conditions, and numerical methods, all in the context of a…
This paper presents both a numerical method for general relativity and an application of that method. The method involves the use of harmonic coordinates in a 3+1 code to evolve the Einstein equations with scalar field matter. In such…
Brown has recently introduced a covariant formulation of the BSSN equations which is well suited for curvilinear coordinate systems. This is particularly desirable as many astrophysical phenomena are symmetric with respect to the rotation…
We solve Einstein's field equations coupled to relativistic hydrodynamics in full 3+1 general relativity to evolve astrophysical systems characterized by strong gravitational fields. We model rotating, collapsing and binary stars by…
We prove that when the equations are restricted to the principal part the standard version of the BSSN formulation of the Einstein equations is equivalent to the NOR formulation for any gauge, and that the KST formulation is equivalent to…
We show that well-posed, conformally-decomposed formulations of the 3+1 Einstein equations can be obtained by densitizing the lapse and by combining the constraints with the evolution equations. We compute the characteristics structure and…
It is often the case in numerical relativity that schemes that are known to be convergent for well posed systems are used in evolutions of weakly hyperbolic (WH) formulations of Einstein's equations. Here we explicitly show that with…