Abstract:One of the main challenges in the numerical modeling of binary neutron-star mergers are long-term simulations of the postmerger remnant over timescales of the order of seconds. When this modeling includes all the aspects of the complex physics accompanying the remnant, the computational costs can easily become enormous. To address this challenge in part, we have developed a novel hybrid approach in which the solution from a general-relativistic magnetohydrodynamics (GRMHD) code solving the full set of the Einstein equations in Cartesian coordinates is coupled with another GRMHD code in which the Einstein equations are solved under the conformally flat condition (CFC). The latter approximation has a long history and has been shown to provide an accurate description of compact objects in nonvacuum spacetimes. An important aspect of the CFC approximation is that the elliptic equations need to be solved only for a fraction of the steps needed for the underlying hydrodynamical/magnetohydrodynamical evolution, thus allowing for a gain in computational efficiency that can be up to a factor of ∼6(230) in three-dimensional (two-dimensional) simulations. We present here the basic features of the new code, the strategies necessary to interface it when importing both two- and three-dimensional data, and a novel and robust approach to the recovery of the primitive variables. To validate our new framework, we have carried out a number of tests with various coordinates systems and different numbers of spatial dimensions, involving a variety of astrophysical scenarios, including the evolution of the postmerger remnant of a binary neutron-star merger over a timescale of one second. Overall, our results show that the new code, bhac +, is able to accurately reproduce the evolution of compact objects in nonvacuum spacetimes and that, when compared with the evolution in full general relativity, the CFC approximation reproduces accurately both the gravitational fields and the matter variables at a fraction of the computational costs. This opens the way for the systematic study of the secular matter and electromagnetic emission from binary-merger remnants. https://doi.org/10.1103/PhysRevD.109.064061 © 2024 American Physical Society

Numerical relativity simulations of compact binaries: comparison of cell- and vertex-centered adaptive meshes

GR-Athena++: Puncture Evolutions on Vertex-centered Oct-tree Adaptive Mesh Refinement

Lessons for adaptive mesh refinement in numerical relativity

Adaptive mesh refinement in binary black holes simulations

A hybrid approach to long-term binary neutron-star simulations

Hybrid approach to long-term binary neutron-star simulations

Performance-Portable Binary Neutron Star Mergers with AthenaK

Improved Moving-Puncture Techniques for Compact Binary Simulations

GR-Athena++: General-relativistic magnetohydrodynamics simulations of neutron star spacetimes

Thinking Outside the Box: Numerical Relativity with Particles

Binary black hole simulation with an adaptive finite element method II: Application of local discontinuous Galerkin method to Einstein equations

Simulating binary black hole mergers using discontinuous Galerkin methods

Insights into Binary Neutron Star Merger Simulations: A Multi-Code Comparison

Inspiraling black-hole binary spacetimes: Challenges in transitioning from analytical to numerical techniques

Binary neutron star mergers using a discontinuous Galerkin-finite difference hybrid method

Initial Data and Eccentricity Reduction Toolkit for Binary Black Hole Numerical Relativity Waveforms

New public code for initial data of unequal-mass, spinning compact-object binaries

Recent Advances in the Numerical Simulations of Binary Black Holes

General Relativistic Magneto-Hydrodynamic Simulations with BAM: Implementation and Code Comparison