Numerical studies of constraints and gravitational wave extraction in general relativity

dc.contributor.advisorMisner, Charles Wen_US
dc.contributor.authorFiske, David Roberten_US
dc.contributor.publisherDigital Repository at the University of Marylanden_US
dc.contributor.publisherUniversity of Maryland (College Park, Md.)en_US
dc.description.abstractWithin classical physics, general relativity is the theory of gravity. Its equations are non-linear partial differential equations for which relatively few closed form solutions are known. Because of the growing observational need for solutions representing gravitational waves from astrophysically plausible sources, a subfield of general relativity, numerical relativity, has a emerged with the goal of generating numerical solutions to the Einstein equations. This dissertation focuses on two fundamental problems in modern numerical relativity: (1) Creating a theoretical treatment of the constraints in the presence of constraint-violating numerical errors, and (2) Designing and implementing an algorithm to compute the spherical harmonic decomposition of radiation quantities for comparison with observation. On the issue of the constraints, I present a novel and generic procedure for incorporating the constraints into the equations of motion of the theory in a way designed to make the constraint hypersurface an attractor of the evolution. In principle, the prescription generates non-linear corrections for the Einstein equations. The dissertation presents numerical evidence that the correction terms do work in the case of two formulations of the Maxwell equations and two formulations of the linearized Einstein equations. On the issue of radiation extraction, I provide the first in-depth analysis of a novel algorithm, due originally to Misner, for computing spherical harmonic components on a cubic grid. I compute explicitly how the truncation error in the algorithm depends on its various parameters, and I also provide a detailed analysis showing how to implement the method on grids in which explicit symmetries are enforced via boundary conditions. Finally, I verify these error estimates and symmetry arguments with a numerical study using a solution of the linearized Einstein equations known as a Teukolsky wave. The algorithm performs well and the estimates prove true both in simulations run on a uniform grid and in simulations that make use of fixed mesh refinement techniques.en_US
dc.format.extent905503 bytes
dc.subject.pqcontrolledPhysics, Astronomy and Astrophysicsen_US
dc.subject.pqcontrolledComputer Scienceen_US
dc.subject.pquncontrolledgeneral relativityen_US
dc.subject.pquncontrolledgravitational waveen_US
dc.subject.pquncontrolledspherical harmonicsen_US
dc.subject.pquncontrolledpartial differential equationsen_US
dc.subject.pquncontrolledcomputer simulationen_US
dc.titleNumerical studies of constraints and gravitational wave extraction in general relativityen_US


Original bundle
Now showing 1 - 1 of 1
Thumbnail Image
884.28 KB
Adobe Portable Document Format