12 research outputs found
Gravitational perturbations of Schwarzschild spacetime at null infinity and the hyperboloidal initial value problem
We study gravitational perturbations of Schwarzschild spacetime by solving a
hyperboloidal initial value problem for the Bardeen-Press equation.
Compactification along hyperboloidal surfaces in a scri-fixing gauge allows us
to have access to the gravitational waveform at null infinity in a general
setup. We argue that this hyperboloidal approach leads to a more accurate and
efficient calculation of the radiation signal than the common approach where a
timelike outer boundary is introduced. The method can be generalized to study
perturbations of Kerr spacetime using the Teukolsky equation.Comment: 14 pages, 9 figure
Implementation of standard testbeds for numerical relativity
We discuss results that have been obtained from the implementation of the
initial round of testbeds for numerical relativity which was proposed in the
first paper of the Apples with Apples Alliance. We present benchmark results
for various codes which provide templates for analyzing the testbeds and to
draw conclusions about various features of the codes. This allows us to sharpen
the initial test specifications, design a new test and add theoretical insight.Comment: Corrected versio
Numerical Relativity: A review
Computer simulations are enabling researchers to investigate systems which
are extremely difficult to handle analytically. In the particular case of
General Relativity, numerical models have proved extremely valuable for
investigations of strong field scenarios and been crucial to reveal unexpected
phenomena. Considerable efforts are being spent to simulate astrophysically
relevant simulations, understand different aspects of the theory and even
provide insights in the search for a quantum theory of gravity. In the present
article I review the present status of the field of Numerical Relativity,
describe the techniques most commonly used and discuss open problems and (some)
future prospects.Comment: 2 References added; 1 corrected. 67 pages. To appear in Classical and
Quantum Gravity. (uses iopart.cls
Stability Analysis of Difference Methods for Parabolic Initial Value Problems
Abstract A decomposition of the numerical solution can be defined by the normal mode representation, that generalizes further the spatial eigenmode decomposition of the von Neumann analysis by taking into account the boundary conditions which are not periodic. In this paper we present some new theoretical results on normal mode analysis for a linear and parabolic initial value problem. Furthermore we suggest an algorithm for the calculation of stability regions based on the normal mode theory
Speech development in children aged 3-4 years
We present a high-order difference method for problems in elastodynamics involving the interaction of waves with highly nonlinear frictional interfaces. We restrict our attention to two-dimensional antiplane problems involving deformation in only one direction. Jump conditions that relate tractions on the interface, or fault, to the relative sliding velocity across it are of a form closely related to those used in earthquake rupture models and other frictional sliding problems. By using summation-by-parts (SBP) finite difference operators and weak enforcement of boundary and interface conditions, a strictly stable method is developed. Furthermore, it is shown that unless the nonlinear interface conditions are formulated in terms of characteristic variables, as opposed to the physical variables in terms of which they are more naturally stated, the semi-discretized system of equations can become extremely stiff, preventing efficient solution using explicit time integrators. The use of SBP operators also provides a rigorously defined energy balance for the discretized problem that, as the mesh is refined, approaches the exact energy balance in the continuous problem. This enables one to investigate earthquake energetics, for example the efficiency with which elastic strain energy released during rupture is converted to radiated energy carried by seismic waves, rather than dissipated by frictional sliding of the fault. These theoretical results are confirmed by several numerical tests in both one and two dimensions demonstrating the computational efficiency, the high-order convergence rate of the method, the benefits of using strictly stable numerical methods for long time integration, and the accuracy of the energy balance.funding agencies|NSF|EAR-0910574EAR-0529922|Southern California Earthquake Center (SCEC)||USGS|07HQAG00081417