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

DOI10.1088/0264-9381/22/24/006zbMATH Open1088.83002arXivgr-qc/0507004OpenAlexW3105039587MaRDI QIDQ3367261FDOQ3367261

Authors: Luis Lehner, Manuel Tiglio, O. Reula 

Publication date: 24 January 2006

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

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.

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

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

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)

Cited In (18)

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

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)