Multi-block simulations in general relativity: high-order discretizations, numerical stability and applications

DOI10.1088/0264-9381/22/24/006MaRDI QIDQ3367261zbMATH OpenOpenAlexFDO

Authors Luis Lehner, Manuel Tiglio, O. Reula

Publication date 24 January 2006

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

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

Mathematics Subject Classification ID

Einstein's equations (general structure, canonical formalism, Cauchy problems) (83C05) Research exposition (monographs, survey articles) pertaining to relativity and gravitational theory (83-02) Computational methods for problems pertaining to relativity and gravitational theory (83-08)

Abstract: The need to smoothly cover a computational domain of interest generically requires the adoption of several grids. To solve the problem of interest under this grid-structure one must ensure the suitable transfer of information among the different grids involved. In this work we discuss a technique that allows one to construct finite difference schemes of arbitrary high order which are guaranteed to satisfy linear numerical and strict stability. The technique relies on the use of difference operators satisfying summation by parts and {it penalty techniques} to transfer information between the grids. This allows the derivation of semidiscrete energy estimates for problems admitting such estimates at the continuum. We analyze several aspects of this technique when used in conjuction with high order schemes and illustrate its use in one, two and three dimensional numerical relativity model problems with non-trivial topologies, including truly spherical black hole excision.

Recommendations

• Novel finite-differencing techniques for numerical relativity: application to black-hole excision
• A multi-block infrastructure for three-dimensional time-dependent numerical relativity
• Summation by parts and dissipation for domains with excised regions
• Black-hole excision with multiple grid patches
• Optimized high-order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions

Cited in

(18)

• Numerical treatment of interfaces for second-order wave equations
• High-fidelity numerical solution of the time-dependent Dirac equation
• Optimal diagonal-norm SBP operators
• A multi-block infrastructure for three-dimensional time-dependent numerical relativity
• The second order spin-2 system in flat space near space-like and null-infinity
• Characteristic evolution and matching
• Diagonal-norm summation by parts operators for finite difference approximations of third and fourth derivatives
• High-fidelity numerical simulation of the dynamic beam equation
• Novel finite-differencing techniques for numerical relativity: application to black-hole excision
• Interface procedures for finite difference approximations of the advection-diffusion equation
• A spectral solver for evolution problems with spatial \(\mathbb S^3\)-topology
• Solving partial differential equations numerically on manifolds with arbitrary spatial topologies
• 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
• Summation by parts and dissipation for domains with excised regions
• Stable boundary treatment for the wave equation on second-order form
• Output error estimation for summation-by-parts finite-difference schemes
• A solution to the stability issues with block norm summation by parts operators

This page was built for publication: Multi-block simulations in general relativity: high-order discretizations, numerical stability and applications

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