A multi-block infrastructure for three-dimensional time-dependent numerical relativity

DOI10.1088/0264-9381/23/16/S14zbMATH Open1191.83006arXivgr-qc/0602104WikidataQ59830079 ScholiaQ59830079MaRDI QIDQ5486415FDOQ5486415

Authors: Erik Schnetter, Peter Diener, Ernst Nils Dorband, Manuel Tiglio 

Publication date: 6 September 2006

Published in: Classical and Quantum Gravity (Search for Journal in Brave)

Abstract: We describe a generic infrastructure for time evolution simulations in numerical relativity using multiple grid patches. After a motivation of this approach, we discuss the relative advantages of global and patch-local tensor bases. We describe both our multi-patch infrastructure and our time evolution scheme, and comment on adaptive time integrators and parallelisation. We also describe various patch system topologies that provide spherical outer and/or multiple inner boundaries. We employ penalty inter-patch boundary conditions, and we demonstrate the stability and accuracy of our three-dimensional implementation. We solve both a scalar wave equation on a stationary rotating black hole background and the full Einstein equations. For the scalar wave equation, we compare the effects of global and patch-local tensor bases, different finite differencing operators, and the effect of artificial dissipation onto stability and accuracy. We show that multi-patch systems can directly compete with the so-called fixed mesh refinement approach; however, one can also combine both. For the Einstein equations, we show that using multiple grid patches with penalty boundary conditions leads to a robustly stable system. We also show long-term stable and accurate evolutions of a one-dimensional non-linear gauge wave. Finally, we evolve weak gravitational waves in three dimensions and extract accurate waveforms, taking advantage of the spherical shape of our grid lines.

Full work available at URL: https://arxiv.org/abs/gr-qc/0602104

Recommendations

• Evolutions in 3D numerical relativity using fixed mesh refinement
• Multi-block simulations in general relativity: high-order discretizations, numerical stability and applications
• Black-hole excision with multiple grid patches
• Characteristic evolutions in numerical relativity using six angular patches
• NUMERICAL RELATIVITY

Mathematics Subject Classification ID

Multigrid methods; domain decomposition for initial value and initial-boundary value problems involving PDEs (65M55) Analogues of general relativity in lower dimensions (83C80) Computational methods for problems pertaining to relativity and gravitational theory (83-08)

Cited In (12)

• Parameter estimates in binary black hole collisions using neural networks
• Characteristic evolution and matching
• Solving partial differential equations numerically on manifolds with arbitrary spatial topologies
• Review of summation-by-parts schemes for initial-boundary-value problems
• Multi-block simulations in general relativity: high-order discretizations, numerical stability and applications
• Continuum and discrete initial-boundary value problems and Einstein's field equations
• The second order spin-2 system in flat space near space-like and null-infinity
• A new code for the numerical simulation of relativistic flows on supercomputers by means of a low-dissipation scheme
• On conservation and stability properties for summation-by-parts schemes
• Optimized high-order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions
• Evolutions in 3D numerical relativity using fixed mesh refinement
• Axisymmetric hydrodynamics in numerical relativity using a multipatch method

Uses Software

• LAPACK
• Cactus
• LAM-MPI
• MPI/MPICH
• Carpet

This page was built for publication: A multi-block infrastructure for three-dimensional time-dependent numerical relativity

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5486415)