Constraint Relaxation

Full relativistic simulations in three dimensions are known to develop runaway modes that grow exponentially and are accompanied by violations of the Hamiltonian and momentum constraints. We present here a method that controls the violation of these constraints and is tested with simulations of binary neutron stars in circular orbits. We show that this technique improves the overall quality of the simulations.