We present a new fully first order strongly hyperbolic representation of the BSSN formulation of Einstein's equations with optional constraint damping terms. We describe the characteristic fields of the system, discuss its hyperbolicity properties, and present two numerical implementations and simulations: one using finite differences, adaptive mesh refinement and in particular binary black holes, and another one using the discontinuous Galerkin method in spherical symmetry. The results of this paper constitute a first step in an effort to combine the robustness of BSSN evolutions with very high accuracy numerical techniques, such as spectral collocation multi-domain or discontinuous Galerkin methods.
@article{arxiv.1202.1038,
title = {Numerical simulations with a first order BSSN formulation of Einstein's field equations},
author = {J. David Brown and Peter Diener and Scott E. Field and Jan S. Hesthaven and Frank Herrmann and Abdul H. Mroué and Olivier Sarbach and Erik Schnetter and Manuel Tiglio and Michael Wagman},
journal= {arXiv preprint arXiv:1202.1038},
year = {2015}
}