Statistics – Computation
Scientific paper
Jan 1993
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1993phdt........21t&link_type=abstract
PhD Dissertation, British Columbia Univ. Vancouver, British Columbia Canada
Statistics
Computation
19
Black Holes (Astronomy), Finite Difference Theory, Relativity, Space-Time Functions, Computational Astrophysics, Event Horizon, Algorithms, Hyperbolic Differential Equations, Singularity (Mathematics), Tensors, Wave Equations, Partial Differential Equations
Scientific paper
This thesis is concerned with the development of better techniques for the 3 + 1 numerical relativity study of black hole spacetimes. The main result of this thesis is the development of a new technique for avoiding singularities in such spacetimes. In this technique the slices are allowed to penetrate the black hole, but only the region of spacetime outside the apparent horizon is numerically evolved. This allows the slicing to be chosen to avoid the 'grid stretching' problems commonly encountered when freezing slicings are used. To implement this scheme, we have developed a robust and efficient apparent-horizon-finding algorithm. We use this at each time step during a numerical evolution to monitor the apparent horizon's position; we then dynamically adjust the region of spacetime excluded from the numerical evolution so this region tracks the apparent horizon's motion. In this thesis we use coordinates in which all components of the metric and other 3 + 1 field tensors are (generally) nonzero. This makes the 3 + 1 equations very complex, so we have developed a prototype 'PDE Compiler' to automatically finite difference them and generate the required code. This automation of the finite differencing process allows us to work with and think about the 3 + 1 equations almost entirely at the tensor-differential-operator level. We have developed a new initial data solver, which numerically solves the full four-vector York equations on slices which are generally not maximal and not three-conformally-flat. Our numerical methods are based on fourth-order finite differencing, using the method of lines for hyperbolic partial differential equations (PDE's). To study these and our black hole exclusion technique in a simple setting, we have made a series of model problem studies using a one-dimensional flat-space scalar wave equation. These have been very successful, yielding a stable and highly accurate finite differencing scheme. To test these techniques in a more realistic setting, we have written a prototype numerical relativity code to simulate the time evolution of axisymmetric (single) black hole spacetimes. At present, this code suffers from severe finite differencing instabilities. We have identified potential causes for some of these instabilities, but we have not yet resolved them.
No associations
LandOfFree
Numerical relativity in black hole spacetimes does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Numerical relativity in black hole spacetimes, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Numerical relativity in black hole spacetimes will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1264245