170 likes | 325 Vues
Numerical Relativity is still Relativity. ERE Salamanca 2008 Palma Group Alic, Dana · Bona, Carles · Bona-Casas, Carles. Most recent successful stories in BH simulations. Long term evolutions: Harmonic (4D spacetime, excision, harmonic gauge source functions)
E N D
Numerical Relativity is still Relativity ERE Salamanca 2008 Palma Group Alic, Dana · Bona, Carles · Bona-Casas, Carles
Most recent successful stories in BH simulations • Long term evolutions: • Harmonic (4D spacetime, excision, harmonic gauge source functions) • BSSN (3+1 decomposition, punctures/excision, 1+log and gamma freezing) • Isn’t the gauge choice too limited? Shouldn’t numerical relativity be relativity?
Do we have any choice? • Reported experiences: • No long term simulations with normal coordinates (zero shift). • Generalised harmonic slicing but strictly harmonic shift. • BSSN normal coordinates (zero shift) and 1+log slicing crashes at 30-40M (gr-qc/0206072). • Gaugewave test: gauge imposed is harmonic, so harmonic code succeeds, but BSSN crashes.
Looking for a gauge polyvalent code • Z4 formalism • MoL with 3rd order SSP Runge-Kutta. • Powerful 3rd order FD algorithm (submitted to JCP). See a variant in http://arxiv.org/abs/0711.4685 (ERE 2007) • Scalar field stuffing. • Cactus. Single grid calculation. Logarithmic grid for long runs.
Gaugewave Test • Minkowski spacetime: • Harmonic coordinates x,y,z,t.
BSSN Comparison t=30 t=1000
Single BH Test • Singularity avoidant conditions (Bona-Massó) Q = f (trK-2) • 1+log (f=2/) slicing with normal coordinates (zero shift) up to 1000M and more! Never done before (BSSN reported to crash at 30-40M without shift). • Unigrid simulation. Logcoords =1.5.
R/M=20 r/M=463000
More gauges (zero shift) • Isotropic coords. Boundaries at 20M. • Logcoords f=1/ 150M.
Shift • 1st order conditions. • Vectorial. • Harmonic? xi = 0. 1st order version
Advection terms • Lie derivative “advection/damping” • Covariant advection term
1st order vector ingredients • Time-independent coordinate transformations.