Binary neutron star mergers using a discontinuous Galerkin-finite difference hybrid method
Abstract
We present a discontinuous Galerkin-finite difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics in dynamical spacetimes. We present several optimizations and stability improvements to our algorithm that allow the hybrid method to successfully simulate single, rotating, and binary neutron stars. The hybrid method achieves the efficiency of discontinuous Galerkin methods throughout almost the entire spacetime during the inspiral phase, while being able to robustly capture shocks and resolve the stellar surfaces. We also use Cauchy-Characteristic evolution to compute the first gravitational waveforms at future null infinity from binary neutron star mergers. The simulations presented here are the first successful binary neutron star inspiral and merger simulations using discontinuous Galerkin methods.
Keywords
Cite
@article{arxiv.2406.19038,
title = {Binary neutron star mergers using a discontinuous Galerkin-finite difference hybrid method},
author = {Nils Deppe and Francois Foucart and Marceline S. Bonilla and Michael Boyle and Nicholas J. Corso and Matthew D. Duez and Matthew Giesler and François Hébert and Lawrence E. Kidder and Yoonsoo Kim and Prayush Kumar and Isaac Legred and Geoffrey Lovelace and Elias R. Most and Jordan Moxon and Kyle C. Nelli and Harald P. Pfeiffer and Mark A. Scheel and Saul A. Teukolsky and William Throwe and Nils L. Vu},
journal= {arXiv preprint arXiv:2406.19038},
year = {2024}
}
Comments
31 pages, 8 figures, minor wording updates