Hybrid approach to long-term binary neutron-star simulations
Harry Ho-Yin Ng,Jin-Liang Jiang,Carlo Musolino,Christian Ecker,Samuel D. Tootle,Luciano Rezzolla
DOI: https://doi.org/10.1103/physrevd.109.064061
IF: 5.407
2024-03-22
Physical Review D
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
astronomy & astrophysics,physics, particles & fields