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

DOI10.1088/0264-9381/23/16/S14MaRDI QIDQ5486415zbMATH OpenWikidataFDO

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)

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

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)

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.

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

Cited in

(12)

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

Describes a project that uses

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)