Geometrically motivated hyperbolic coordinate conditions for numerical relativity: Analysis, issues and implementations
Abstract
We study the implications of adopting hyperbolicdriver coordinate conditions motivated by geometrical considerations. In particular, conditions that minimize the rate of change of the metric variables. We analyze the properties of the resulting system of equations and their effect when implementing excision techniques. We find that commonly used coordinate conditions lead to a characteristic structure at the excision surface where some modes are not of outflow type with respect to any excision boundary chosen inside the horizon. Thus, boundary conditions are required for these modes. Unfortunately, the specification of these conditions is a delicate issue as the outflow modes involve both gauge and main variables. As an alternative to these driver equations, we examine conditions derived from extremizing a scalar constructed from Killing's equation and present specific numerical examples.
 Authors:
 Departament de Fisica, Universitat de les Illes Balears, Palma de Mallorca (Spain)
 Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 708034001 (United States)
 Publication Date:
 OSTI Identifier:
 20711562
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Physical Review. D, Particles Fields; Journal Volume: 72; Journal Issue: 10; Other Information: DOI: 10.1103/PhysRevD.72.104009; (c) 2005 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; BOUNDARY CONDITIONS; COORDINATES; COSMOLOGY; GENERAL RELATIVITY THEORY; GEOMETRY; SCALARS
Citation Formats
Bona, Carles, Lehner, Luis, and PalenzuelaLuque, Carlos. Geometrically motivated hyperbolic coordinate conditions for numerical relativity: Analysis, issues and implementations. United States: N. p., 2005.
Web. doi:10.1103/PhysRevD.72.104009.
Bona, Carles, Lehner, Luis, & PalenzuelaLuque, Carlos. Geometrically motivated hyperbolic coordinate conditions for numerical relativity: Analysis, issues and implementations. United States. doi:10.1103/PhysRevD.72.104009.
Bona, Carles, Lehner, Luis, and PalenzuelaLuque, Carlos. Tue .
"Geometrically motivated hyperbolic coordinate conditions for numerical relativity: Analysis, issues and implementations". United States.
doi:10.1103/PhysRevD.72.104009.
@article{osti_20711562,
title = {Geometrically motivated hyperbolic coordinate conditions for numerical relativity: Analysis, issues and implementations},
author = {Bona, Carles and Lehner, Luis and PalenzuelaLuque, Carlos},
abstractNote = {We study the implications of adopting hyperbolicdriver coordinate conditions motivated by geometrical considerations. In particular, conditions that minimize the rate of change of the metric variables. We analyze the properties of the resulting system of equations and their effect when implementing excision techniques. We find that commonly used coordinate conditions lead to a characteristic structure at the excision surface where some modes are not of outflow type with respect to any excision boundary chosen inside the horizon. Thus, boundary conditions are required for these modes. Unfortunately, the specification of these conditions is a delicate issue as the outflow modes involve both gauge and main variables. As an alternative to these driver equations, we examine conditions derived from extremizing a scalar constructed from Killing's equation and present specific numerical examples.},
doi = {10.1103/PhysRevD.72.104009},
journal = {Physical Review. D, Particles Fields},
number = 10,
volume = 72,
place = {United States},
year = {Tue Nov 15 00:00:00 EST 2005},
month = {Tue Nov 15 00:00:00 EST 2005}
}

Einstein evolution equations are written as a hyperbolic system of balance laws. A harmonic time coordinate is used with zero shift vector (harmonic slicing). The principal part of the evolution system reduces to a set of uncoupled wave equations in first order form. The relevance for threedimensional numerical relativity of both the harmonic slicing and the resulting evolution system is stressed.

Iterated CrankNicolson method for hyperbolic and parabolic equations in numerical relativity
The iterated CrankNicolson is a predictorcorrector algorithm commonly used in numerical relativity for the solution of both hyperbolic and parabolic partial differential equations. We here extend the recent work on the stability of this scheme for hyperbolic equations by investigating the properties when the average between the predicted and corrected values is made with unequal weights and when the scheme is applied to a parabolic equation. We also propose a variant of the scheme in which the coefficients in the averages are swapped between two corrections leading to systematically larger amplification factors and to a smaller numerical dispersion. 
Comparative studies on nonlinear hyperbolic and parabolic heat conduction for various boundary conditions: Analytic and numerical solutions
With the advent of lasers with very short pulse durations and their use in materials processing, the effect of thermal wave propagation velocity becomes important. Also, localized heating in laseraided materials processing causes significant variations in the material properties. To account for these two effects, hyperbolic heat conduction is studied in this paper by considering all the thermophysical properties, except the thermal diffusivity, to be temperature dependent. The resulting nonlinear hyperbolic equations are linearized by using Kirchhoff transformation. Both analytical and numerical solutions are obtained for finite domains. Results are presented and compared with parabolic conduction results. 
Boundary Conditions for Constrained Systems of Evolution Equations in Numerical Relativity
In numerical relativity, to study astrophysical events such as blackhole collisions, the Einstein equations for geometry of spacetime are solved as system of partial differential equations. Current simulations are usually based on socalled BSSN system of 3+1 constrained hyperbolic evolution equations for tensorial fields of various ranks. The boundary conditions for evolved fields are given by the fact that the simulated events happen in an empty space and that far from the center the waves should propagate outwards. We analyze the usual outgoing radiation boundary conditions for the BSSN variables using potentials. Introduction of these potentials simplifies the linearized BSSNmore »