-
Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations
Authors:
Manuel R. Izquierdo,
Lorenzo Pareschi,
Borja Miñano,
Joan Massó,
Carlos Palenzuela
Abstract:
Modeling correctly the transport of neutrinos is crucial in some astrophysical scenarios such as core-collapse supernovae and binary neutron star mergers. In this paper, we focus on the truncated-moment formalism, considering only the first two moments (M1 scheme) within the grey approximation, which reduces Boltzmann seven-dimensional equation to a system of $3+1$ equations closely resembling the…
▽ More
Modeling correctly the transport of neutrinos is crucial in some astrophysical scenarios such as core-collapse supernovae and binary neutron star mergers. In this paper, we focus on the truncated-moment formalism, considering only the first two moments (M1 scheme) within the grey approximation, which reduces Boltzmann seven-dimensional equation to a system of $3+1$ equations closely resembling the hydrodynamic ones. Solving the M1 scheme is still mathematically challenging, since it is necessary to model the radiation-matter interaction in regimes where the evolution equations become stiff and behave as an advection-diffusion problem. Here, we present different global, high-order time integration schemes based on Implicit-Explicit Runge-Kutta (IMEX) methods designed to overcome the time-step restriction caused by such behavior while allowing us to use the explicit RK commonly employed for the MHD and Einstein equations. Finally, we analyze their performance in several numerical tests.
△ Less
Submitted 29 May, 2023; v1 submitted 31 October, 2022;
originally announced November 2022.
-
Large Eddy Simulations of Magnetized Mergers of Neutron Stars with Neutrinos
Authors:
C. Palenzuela,
S. L. Liebling,
B. Miñano
Abstract:
Neutron star mergers are very violent events involving extreme physical processes: dynamic, strong-field gravity, large magnetic field, very hot, dense matter, and the copious production of neutrinos. Accurate modeling of such a system and its associated multi-messenger signals, such as gravitational waves, short gamma ray burst, and kilonova, requires the inclusion of all these processes, and is…
▽ More
Neutron star mergers are very violent events involving extreme physical processes: dynamic, strong-field gravity, large magnetic field, very hot, dense matter, and the copious production of neutrinos. Accurate modeling of such a system and its associated multi-messenger signals, such as gravitational waves, short gamma ray burst, and kilonova, requires the inclusion of all these processes, and is increasingly important in light of advancements in multi-messenger astronomy generally, and in gravitational wave astronomy in particular (such as the development of third-generation detectors). Several general relativistic codes have been incorporating some of these elements with different levels of realism. Here, we extend our code MHDuet, which can perform large eddy simulations of magnetohydrodynamics to help capture the magnetic field amplification during the merger, and to allow for realistic equations of state and neutrino cooling via a leakage scheme. We perform several tests involving isolated and binary neutron stars demonstrating the accuracy of the code.
△ Less
Submitted 6 May, 2022; v1 submitted 6 April, 2022;
originally announced April 2022.
-
Turbulent magnetic field amplification in binary neutron star mergers
Authors:
C. Palenzuela,
R. Aguilera-Miret,
F. Carrasco,
R. Ciolfi,
J. V. Kalinani,
W. Kastaun,
B. Miñano,
D. Viganò
Abstract:
Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insu…
▽ More
Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insufficient to fully capture the small-scale dynamo, initially triggered by the Kelvin-Helmholtz instability, and later sourced by several MHD processes involving differential rotation. Here, we alleviate this limitation by using explicit large-eddy simulations, a technique where the unresolved dynamics occurring at the sub-grid scales (SGS) is modeled by extra terms, which are functions of the resolved fields and their derivatives. The combination of high-order numerical schemes, high resolutions, and the gradient SGS model allow us to capture the small-scale dynamos produced during the binary neutron star mergers. Here we follow the first 50 milliseconds after the merger and, for the first time, we find numerical convergence on the magnetic field amplification, in terms of integrated energy and spectral distribution over spatial scales. We also find that the average intensity of the magnetic field in the remnant saturates at $\sim 10^{16}$~G around $5$~ms after the merger. After $20-30$~ms, both toroidal and poloidal magnetic field components grow continuously, fed by the winding mechanism that provides a slow inverse cascade. We find no clear hints for magneto-rotational instabilities, and no significant impact of the magnetic field on the redistribution of angular momentum in the remnant in our simulations, probably due to the very turbulent and dynamical topology of the magnetic field at all stages, with small-scale components largely dominating over the large-scale ones.
△ Less
Submitted 3 July, 2022; v1 submitted 15 December, 2021;
originally announced December 2021.
-
Turbulent magnetic-field amplification in the first 10 milliseconds after a binary neutron star merger: comparing high-resolution and large eddy simulations
Authors:
Ricard Aguilera-Miret,
Daniele Viganò,
Federico Carrasco,
Borja Miñano,
Carlos Palenzuela
Abstract:
The detection of binary neutron star mergers represents one of the most important and complex astrophysical discoveries of the recent years. One of the unclear aspects of the problem is the turbulent magnetic field amplification, initially triggered by the Kelvin-Helmholtz instability at much smaller scales than any reachable numerical resolution nowadays. Here we present numerical simulations of…
▽ More
The detection of binary neutron star mergers represents one of the most important and complex astrophysical discoveries of the recent years. One of the unclear aspects of the problem is the turbulent magnetic field amplification, initially triggered by the Kelvin-Helmholtz instability at much smaller scales than any reachable numerical resolution nowadays. Here we present numerical simulations of the first ten milliseconds of a binary neutron star merger. First, we confirm in detail how the simulated amplification depends on the numerical resolution and is distributed on a broad range of scales, as expected from turbulent MHD theory. We find that an initial large-scale magnetic field of $10^{11}\,$G inside each star is amplified in the remnant to root-mean-square values above $10^{16}\,$G within the first $5$ milliseconds for our highest-resolution run. Then, we run large eddy simulations, exploring the performance of the subgrid-scale gradient model, already tested successfully in previous turbulent box simulations. We show that the addition of this model is especially important in the induction equation, since it leads to an amplification of the magnetic field comparable to a higher-resolution run, but with a greatly reduced computational cost. In the first 10 milliseconds, there is no clear hint for an ordered, large-scale magnetic field, which should indeed occur in longer timescales through magnetic winding and the magneto-rotational instability.
△ Less
Submitted 14 September, 2020;
originally announced September 2020.
-
GRMHD large eddy simulations with gradient subgrid-scale model
Authors:
Daniele Viganò,
Ricard Aguilera-Miret,
Federico Carrasco,
Borja Miñano,
Carlos Palenzuela
Abstract:
The detection of binary neutron star mergers represents one of the most important astrophysical discoveries of the recent years. Due to the extreme matter and gravity conditions and the rich dynamics developed, it becomes a tremendous challenge to accurately simulate numerically all the scales present during the collision. Here we present how to study such systems by using large eddy simulations w…
▽ More
The detection of binary neutron star mergers represents one of the most important astrophysical discoveries of the recent years. Due to the extreme matter and gravity conditions and the rich dynamics developed, it becomes a tremendous challenge to accurately simulate numerically all the scales present during the collision. Here we present how to study such systems by using large eddy simulations with a self-consistent subgrid-scale gradient model, that we generalized to the special relativistic case in a previous work and now extend to the general relativistic case. Adapted from nonrelativistic scenarios, the so-called gradient model allows to capture part of the effects of the hidden dynamics on the resolved scales, by means of a physically-agnostic, mathematically-based Taylor expansion of the nonlinear terms in the conservative evolution equations' fluxes. We assess the validity of this approach in bounding-box simulations of the magnetic Kelvin-Helmholtz instability. Several resolutions and a broad range of scenarios are considered in order to carefully test the performance of the model under three crucial aspects: (i) highly curved backgrounds, (ii) jumps on the fluid density profiles and (iii) strong shocks. The results suggest our extension of the gradient subgrid-scale model to general relativistic magnetohydrodynamics is a promising approach for studying binary neutron stars mergers, and potentially to other relevant astrophysical scenarios.
△ Less
Submitted 2 April, 2020;
originally announced April 2020.
-
A Simflowny-based high-performance 3D code for the generalized induction equation
Authors:
Daniele Viganò,
David Martínez-Gómez,
José A. Pons,
Carlos Palenzuela,
Federico Carrasco,
Borja Miñano,
Antoni Arbona,
Carles Bona,
Joan Massó
Abstract:
In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage…
▽ More
In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage and a consensus on the general dynamics at least for some simple initial data, only few attempts have been performed in 3D, due to the computational costs and the need for a proper numerical treatment of the intrinsic non-linearity of the equation. Here, we carefully analyze the general induction equation, studying its characteristic structure, and we present a new Cartesian 3D code, generated by the user-friendly, publicly available {\em Simflowny} platform. The code uses high-order numerical schemes for the time and spatial discretization, and relies on the highly-scalable {\em SAMRAI} architecture for the adaptive mesh refinement. We present the application of the code to several benchmark tests, showing the high order of convergence and accuracy achieved and the capabilities in terms of magnetic shock resolution and three-dimensionality. This paper paves the way for the applications to a realistic, 3D long-term evolution of neutron stars interior and, possibly, of other astrophysical sources.
△ Less
Submitted 11 December, 2018; v1 submitted 20 November, 2018;
originally announced November 2018.