skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Geometrically motivated hyperbolic coordinate conditions for numerical relativity: Analysis, issues and implementations

Abstract

We study the implications of adopting hyperbolic-driver 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:
 [1]; ;  [2]
  1. Departament de Fisica, Universitat de les Illes Balears, Palma de Mallorca (Spain)
  2. Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803-4001 (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 Palenzuela-Luque, 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, & Palenzuela-Luque, 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 Palenzuela-Luque, 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 Palenzuela-Luque, Carlos},
abstractNote = {We study the implications of adopting hyperbolic-driver 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 three-dimensional numerical relativity of both the harmonic slicing and the resulting evolution system is stressed.
  • No abstract prepared.
  • The iterated Crank-Nicolson is a predictor-corrector 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.
  • 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 laser-aided 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.
  • In numerical relativity, to study astrophysical events such as black-hole collisions, the Einstein equations for geometry of spacetime are solved as system of partial differential equations. Current simulations are usually based on so-called 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 » system into a set of wave equations. We also devise modifications of boundary conditions which prevent creation of constraint violations at boundary due to the conversion of outgoing constrained waves into ingoing unconstrained ones.« less