Numerical relativity from a gauge theory perspective

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Physics, 2010.This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections.Cataloged from student submitted PDF version of thesis.Includes bibliographical references (p. 121-123).I present a new method for numerical simulations of general relativistic systems that eliminates constraint violating modes without the need for constraint damping or the introduction of extra dynamical fields. The method is a type of variational integrator. It is based on a discretization of an action for gravity (the Plebański action) on an unstructured mesh that preserves the local Lorentz transformation and diffeomorphism symmetries of the continuous action. Applying Hamilton's principle of stationary action gives discrete field equations on the mesh. For each gauge degree of freedom there is a corresponding discrete constraint; the remaining discrete evolution equations exactly preserve these constraints under time-evolution. I validate the method using simulations of several analytically solvable spacetimes: a weak gravitational wave spacetime, the Schwarzschild spacetime, and the Kerr spacetime.by Will M. Farr.Ph.D

dfm/emcee: emcee v3.0.1

Added support for long double dtypes Prepared manuscript to submit to JOSS Improved packaging and release infrastructure Fixed bug in initial linear dependence tes