Source author record

Gioel Calabrese

Gioel Calabrese appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

gr-qc

Catalog footprint

What is connected

3works

1topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2005arXiv

Constraint damping in the Z4 formulation and harmonic gauge

We show that by adding suitable lower-order terms to the Z4 formulation of the Einstein equations, all constraint violations except constant modes are damped. This makes the Z4 formulation a particularly simple example of a lambda-system as suggested by Brodbeck et al. We also show that the Einstein equations in harmonic coordinates can be obtained from the Z4 formulation by a change of variables that leaves the implied constraint evolution system unchanged. Therefore the same method can be used to damp all constraints in the Einstein equations in harmonic gauge.

preprint2004arXiv

Finite differencing second order systems describing black hole spacetimes

Keeping Einstein's equations in second order form can be appealing for computational efficiency, because of the reduced number of variables and constraints. Stability issues emerge, however, which are not present in first order formulations. We show that a standard discretization of the second order ``shifted'' wave equation leads to an unstable semi-discrete scheme if the shift parameter is too large. This implies that discretizations obtained using integrators such as Runge-Kutta, Crank-Nicholson, leap-frog are unstable for any fixed value of the Courant factor. We argue that this situation arises in numerical relativity, particularly in simulations of spacetimes containing black holes, and discuss several ways of circumventing this problem. We find that the first order reduction in time based on ``ADM'' type variables is very effective.

preprint2002arXiv

Convergence and stability in numerical relativity

It is often the case in numerical relativity that schemes that are known to be convergent for well posed systems are used in evolutions of weakly hyperbolic (WH) formulations of Einstein's equations. Here we explicitly show that with several of the discretizations that have been used through out the years, this procedure leads to non-convergent schemes. That is, arbitrarily small initial errors are amplified without bound when resolution is increased, independently of the amount of numerical dissipation introduced. The lack of convergence introduced by this instability can be particularly subtle, in the sense that it can be missed by several convergence tests, especially in 3+1 dimensional codes. We propose tests and methods to analyze convergence that may help detect these situations.