Einstein's Equations with Asymptotically Stable Constraint Propagation
  We introduce a proposal to modify Einstein's equations by embedding them in a larger symmetric hyperbolic system. The additional dynamical variables of the modified system are essentially first integrals of the original constraints. The extended system of equations reproduces the usual dynamics on the constraint surface of general relativity, and therefore naturally includes the solutions to Einstein gravity. The main feature of this extended system is that, at least for a linearized version of it, the constraint surface is an attractor of the time evolution. This feature suggests that this system may be a useful alternative to Einstein's equations when obtaining numerical solutions to full, non-linear gravity. 
Othmar Brodbeck, Simonetta Frittelli, Peter Huebner, Oscar A. Reula