New, efficient, and accurate high order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions
arXiv:gr-qc/0512001 · doi:10.1007/s10915-006-9123-7
Abstract
We construct new, efficient, and accurate high-order finite differencing operators which satisfy summation by parts. Since these operators are not uniquely defined, we consider several optimization criteria: minimizing the bandwidth, the truncation error on the boundary points, the spectral radius, or a combination of these. We examine in detail a set of operators that are up to tenth order accurate in the interior, and we surprisingly find that a combination of these optimizations can improve the operators' spectral radius and accuracy by orders of magnitude in certain cases. We also construct high-order dissipation operators that are compatible with these new finite difference operators and which are semi-definite with respect to the appropriate summation by parts scalar product. We test the stability and accuracy of these new difference and dissipation operators by evolving a three-dimensional scalar wave equation on a spherical domain consisting of seven blocks, each discretized with a structured grid, and connected through penalty boundary conditions.
16 pages, 9 figures. The files with the coefficients for the derivative and dissipation operators can be accessed by downloading the source code for the document. The files are located in the "coeffs" subdirectory
References in corpus (5)
- Evolutions in 3D numerical relativity using fixed mesh refinement
- Gravitational wave astronomy
- A multi-block infrastructure for three-dimensional time-dependent numerical relativity
- Multi-block simulations in general relativity: high order discretizations, numerical stability, and applications
- Cauchy-perturbative matching revisited: tests in spherical symmetry
Cited by in corpus (46)
- Review of Summation-by-parts schemes for initial-boundary-value problems
- Conformal and covariant formulation of the Z4 system with constraint-violation damping
- Compact binary evolutions with the Z4c formulation
- Eccentric, nonspinning, inspiral, Gaussian-process merger approximant for the detection and characterization of eccentric binary black hole mergers
- High accuracy binary black hole simulations with an extended wave zone
- Continuum and Discrete Initial-Boundary-Value Problems and Einstein's Field Equations
- Summation-By-Parts Operators and High-Order Quadrature
- Stability of general-relativistic accretion disks
- Characteristic Evolution and Matching
- Three-Dimensional General-Relativistic Hydrodynamic Simulations of Binary Neutron Star Coalescence and Stellar Collapse with Multipatch Grids
- A numerical study of the quasinormal mode excitation of Kerr black holes
- Constraint preserving boundary conditions for the Z4c formulation of general relativity
- Triggering magnetar outbursts in 3D force-free simulations
- Nonlinear/non-iterative treatment of EFT-motivated gravity
- Physics of eccentric binary black hole mergers: A numerical relativity perspective
- Multi-patch methods in general relativistic astrophysics - I. Hydrodynamical flows on fixed backgrounds
- Self-force with (3+1) codes: a primer for numerical relativists
- Numerical relativity with characteristic evolution, using six angular patches
- Striped Blandford/Znajek jets from advection of small scale magnetic field
- Fusing numerical relativity and deep learning to detect higher-order multipole waveforms from eccentric binary black hole mergers
- Magnetosphere of an orbiting neutron star
- Orbiting binary black hole evolutions with a multipatch high order finite-difference approach
- Pulsar magnetospheres in General Relativity
- Magnetospheres of black hole-neutron star binaries
- General relativistic null-cone evolutions with a high-order scheme
- A novel scheme for simulating the force-free equations: boundary conditions and the evolution of solutions towards stationarity
- The second order spin-2 system in flat space near space-like and null-infinity
- Toward a Consistent Framework for High Order Mesh Refinement Schemes in Numerical Relativity
- A spectral solver for evolution problems with spatial S3-topology
- Numerical space-times near space-like and null infinity. The spin-2 system on Minkowski space
- Astrophysical jets from boosted compact objects
- Solving the Einstein constraint equations on multi-block triangulations using finite element methods
- Scalar self-force for eccentric orbits around a Schwarzschild black hole
- Effects of Onset of Phase Transition on Binary Neutron Star Mergers
- Evolution of scalar fields surrounding black holes on compactified constant mean curvature hypersurfaces
- Supporting High-Performance and High-Throughput Computing for Experimental Science
- Asymptotics and singularities in cosmological models with positive cosmological constant
- COFFEE -- An MPI-parallelized Python package for the numerical evolution of differential equations
- Summation by parts methods for the spherical harmonic decomposition of the wave equation in arbitrary dimensions
- An algorithmic exploration of the existence of high-order summation by parts operators with diagonal norm
- Reflecting boundary conditions in numerical relativity as a model for black hole echoes
- Computational General Relativistic Force-Free Electrodynamics: I. Multi-Coordinate Implementation and Testing
- Gravitational wave signatures of phase transition from hadronic to quark matter in isolated neutron stars and binaries
- Stable and Accurate Interpolation Operators for High-Order Multi-Block Finite-Difference Methods
- Numerical treatment of interfaces for second-order wave equations
- A new class of high-order summation by parts finite-difference schemes