Numerical relativity from a gauge theory perspective
Author(s)Farr, Will M. (Will Meierjurgen)
Massachusetts Institute of Technology. Dept. of Physics.
MetadataShow full item record
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.
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).
DepartmentMassachusetts Institute of Technology. Dept. of Physics.; Massachusetts Institute of Technology. Department of Physics
Massachusetts Institute of Technology