An official website of the United States government
Abstract Binary black holes are the most abundant source of gravitational-wave observations. Gravitational-wave observatories in the next decade will require tremendous increases in the accuracy of numerical waveforms modeling binary black holes, compared to today’s state of the art. One approach to achieving the required accuracy is using spectral-type methods that scale to many processors. Using theSpECTREnumerical-relativity (NR) code, we present the first simulations of a binary black hole inspiral, merger, and ringdown using discontinuous Galerkin (DG) methods. The efficiency of DG methods allows us to evolve the binary through ∼ 18 orbits at reasonable computational cost. We then useSpECTRE’s Cauchy Characteristic Evolution (CCE) code to extract the gravitational waves at future null infinity. The open-source nature ofSpECTREmeans this is the first time a spectral-type method for simulating binary black hole evolutions is available to the entire NR community.
more »
« less
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.
more »
« less
- PAR ID:
- 10554591
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Classical and Quantum Gravity
- Volume:
- 41
- Issue:
- 24
- ISSN:
- 0264-9381
- Format(s):
- Medium: X Size: Article No. 245002
- Size(s):
- Article No. 245002
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present a discontinuous Galerkin (DG)–finite difference (FD) hybrid scheme that allows high-order shock capturing with the DG method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple. An unlimited DG candidate solution is computed for the next time step. If the candidate solution is inadmissible, the time step is retaken using robust FD methods. Because of its a posteriori nature, the hybrid scheme inherits the best properties of both methods. It is high-order with exponential convergence in smooth regions, while robustly handling discontinuities. We give a detailed description of how we transfer the solution between the DG and FD solvers, and the troubled-cell indicators necessary to robustly handle slow-moving discontinuities and simulate magnetized neutron stars. We demonstrate the efficacy of the proposed method using a suite of standard and very challenging 1D, 2D, and 3D relativistic magnetohydrodynamics test problems. The hybrid scheme is designed from the ground up to efficiently simulate astrophysical problems such as the inspiral, coalescence, and merger of two neutron stars.more » « less
-
Abstract Interpreting gravitational wave observations and understanding the physics of astrophysical compact objects such as black holes or neutron stars requires accurate theoretical models. Here, we present a new numerical relativity computer program, called Nmesh , that has the design goal to become a next generation program for the simulation of challenging relativistic astrophysics problems such as binary black hole or neutron star mergers. In order to efficiently run on large supercomputers, Nmesh uses a discontinuous Galerkin method together with a domain decomposition and mesh refinement that parallelizes and scales well. In this work, we discuss the various numerical methods we use. We also present results of test problems such as the evolution of scalar waves, single black holes and neutron stars, as well as shock tubes. In addition, we introduce a new positivity limiter that allows us to stably evolve single neutron stars without an additional artificial atmosphere, or other more traditional limiters.more » « less
-
null (Ed.)Binary neutron star mergers provide a unique probe of the dense-matter equation of state (EoS) across a wide range of parameter space, from the zero-temperature EoS during the inspiral to the high-temperature EoS following the merger. In this paper, we implement a new model for calculating parametrized finite-temperature EoS effects into numerical relativity simulations. This "M* model" is based on a two-parameter approximation of the particle effective mass and includes the leading-order effects of degeneracy in the thermal pressure and energy. We test our numerical implementation by performing evolutions of rotating single stars with zero- and non-zero temperature gradients, as well as evolutions of binary neutron star mergers. We find that our new finite-temperature EoS implementation can support stable stars over many dynamical timescales. We also perform a first parameter study to explore the role of the M* parameters in binary neutron star merger simulations. All simulations start from identical initial data with identical cold EoSs, and differ only in the thermal part of the EoS. We find that both the thermal profile of the remnant and the post-merger gravitational wave signal depend on the choice of M* parameters, but that the total merger ejecta depends only weakly on the finite-temperature part of the EoS across a wide range of parameters. Our simulations provide a first step toward understanding how the finite-temperature properties of dense matter may affect future observations of binary neutron star mergers.more » « less
-
ABSTRACT Precursors have been observed seconds to minutes before some short gamma-ray bursts. While the precursor origins remain unknown, one explanation relies on the resonance of neutron star pulsational modes with the tidal forces during the inspiral phase of a compact binary merger. In this paper, we present a model for short gamma-ray burst precursors that relies on tidally resonant neutron star oceans. In this scenario, the onset of tidal resonance in the crust–ocean interface mode ignites the precursor flare, possibly through the interaction between the excited neutron star ocean and the surface magnetic fields. From just the precursor total energy, the time before the main event, and a detected quasi-periodic oscillation frequency, we may constrain the binary parameters and neutron star ocean properties. Our model can immediately distinguish neutron star–black hole mergers from binary neutron star mergers without gravitational wave detection. We apply our model to GRB 211211A, the recently detected long duration short gamma-ray burst with a quasi-periodic precursor, and explore the parameters of this system. The precursor of GRB 211211A is consistent with a tidally resonant neutron star ocean explanation that requires an extreme mass ratio neutron star–black hole merger and a high-mass neutron star. While difficult to reconcile with the main gamma-ray burst and associated kilonova, our results constrain the possible precursor mechanisms in this system. A systematic study of short gamma-ray burst precursors with the model presented here can test precursor origin and probe the possible connection between gamma-ray bursts and neutron star–black hole mergers.more » « less
