Singularity Avoidance in Numerical Black Hole Spacetimes
  Spacetime singularities in numerical relativity can be avoided by excising a region of the computational domain from inside the apparent horizon. We report on results of such a scheme that is based on using ({\it i}) a horizon locking coordinate which locks the coordinate system to the geometry, and ({\it ii}) a finite differencing scheme which respects the causal structure of the spacetime. With this technique a black hole can be evolved accurately well beyond $t=1000M$, where $M$ is the black hole mass. 
Peter Anninos, Greg Daues, Joan Masso, Edward Seidel, Wai-Mo Suen