Search

Errors due to imperfect boundary conditions in numerical relativity simulations of binary black holes can produce unphysical reflections of gravitational waves which compromise the accuracy of waveform predictions, especially for subdominant modes. A system of higher order absorbing boundary conditions which greatly reduces this problem was introduced in earlier work [arXiv:gr-qc/0608051]. In this paper, we devise two new implementations of this boundary condition system in the Spectral Einstein Code (SpEC), and test them in both linear multipolar gravitational wave and inspiralling mass ratio 7:1 binary black hole simulations. One of our implementations in particular is shown to be extremely robust and to produce accuracy superior to the standard freezing-Psi_0 boundary condition usually used by SpEC., Comment: 24 pages, 11 figures (Figs. 9 and 11 showing CCE and extrapolation results are new), accepted by Classical and Quantum Gravity

Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to formulate on complex meshes. We perform a detailed comparison of various limiting strategies presented in the literature applied to the equations of general relativistic magnetohydrodynamics. We compare the standard minmod/$\Lambda\Pi^N$ limiter, the hierarchical limiter of Krivodonova, the simple WENO limiter, the HWENO limiter, and a discontinuous Galerkin-finite-difference hybrid method. The ultimate goal is to understand what limiting strategies are able to robustly simulate magnetized TOV stars without any fine-tuning of parameters. Among the limiters explored here, the only limiting strategy we can endorse is a discontinuous Galerkin-finite-difference hybrid method., Comment: matches published version. 21 pages, 12 figures. Added KH instability, TOV runs with classical limiters. arXiv admin note: text overlap with arXiv:2109.11645