Dynamics of Gravitational Waves in 3D: Formulations, Methods, and Tests

The dynamics of gravitational waves is investigated in full 3+1 dimensional numerical relativity, emphasizing the difficulties that one might encounter in numerical evolutions, particularly those arising from non-linearities and gauge degrees of freedom. Using gravitational waves with amplitudes low enough that one has a good understanding of the physics involved, but large enough to enable non-linear effects to emerge, we study the coupling between numerical errors, coordinate effects, and the nonlinearities of the theory. We discuss the various strategies used in identifying specific features of the evolution. We show the importance of the flexibility of being able to use different numerical schemes, different slicing conditions, different formulations of the Einstein equations (standard ADM vs. first order hyperbolic), and different sets of equations (linearized vs. full Einstein equations). A non-linear scalar field equation is presented which captures some properties of the full Einstein equations, and has been useful in our understanding of the coupling between finite differencing errors and non-linearites. We present a set of monitoring devices which have been crucial in our studying of the waves, including Riemann invariants, pseudo-energy momentum tensor, hamiltonian constraint violation, and fourier spectrum analysis.