Gradient sub-grid-scale model for relativistic MHD Large Eddy Simulations
arXiv:1908.01419 · doi:10.1103/PhysRevD.101.063003
Abstract
MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the computational resources are used to evolve the system only up to a finite grid size. The resolution is not fine enough to capture all the relevant small-scale physics at play, which is instead effectively modeled by a set of additional terms in the evolution equations, dubbed as sub-grid-scale model. Here we extend such approach, commonly used in non-relativistic/non-magnetic/incompressible fluid dynamics, applying the so-called gradient model to a general set of balance-law equations, that includes the relevant case in which a non-trivial inversion of conserved to primitive fields is needed. In particular, we focus on the relativistic compressible ideal MHD scenario, providing for the first time (and for any equation of state) all the additional sub-grid-scale terms. As an application, we consider box simulations of the relativistic Kelvin-Helmholtz instability, which is also the first mechanism responsible for the magnetic field amplification in binary neutron star mergers and cannot yet be fully captured by the finest-grid and longest simulations available. The performance of our model is numerically assessed by comparing it to the residuals arising from the filtering of high-resolution simulations. We find that the model can fit very well those residuals from resolutions a few times higher. Although the application shown here explicitly considers the Minkowski metric, it can be directly extended to general relativity, thus settling the basis to implement the gradient sub-grid model in a GRMHD binary merger. Our results suggest that this approach will be potentially able to unveil much better the small-scale dynamics achievable in full GRMHD simulations.
15 pages, 6 figures. Final version. arXiv admin note: text overlap with arXiv:1904.04099
References in corpus (18)
- GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
- Multi-messenger Observations of a Binary Neutron Star Merger
- Constraining the Maximum Mass of Neutron Stars From Multi-Messenger Observations of GW170817
- GW170817: Joint Constraint on the Neutron Star Equation of State from Multimessenger Observations
- Neutron-star radius constraints from GW170817 and future detections
- Producing ultra-strong magnetic fields in neutron star mergers
- AT2017gfo: an anisotropic and three-component kilonova counterpart of GW170817
- Estimating the Contribution of Dynamical Ejecta in the Kilonova Associated with GW170817
- Magnetized Neutron Star Mergers and Gravitational Wave Signals
- General-Relativistic Large-Eddy Simulations of Binary Neutron Star Mergers
- Binary Neutron Star Mergers and Short Gamma-Ray Bursts: Effects of Magnetic Field Orientation, Equation of State, and Mass Ratio
- Gravitational waves from remnant massive neutron stars of binary neutron star merger: Viscous hydrodynamics effects
- Energy transfer in compressible magnetohydrodynamic turbulence
- Universality and intermittency in relativistic turbulent flows of a hot plasma
- Comparative statistics of selected subgrid-scale models in large eddy simulations of decaying, supersonic MHD turbulence
- Simflowny 3: An upgraded platform for scientific modelling and simulation
- Nonlinear closures for scale separation in supersonic magnetohydrodynamic turbulence
- Massively parallel implementation in Python of a pseudo-spectral DNS code for turbulent flows
Cited by in corpus (24)
- Turbulent magnetic field amplification in binary neutron star mergers
- The key role of magnetic fields in binary neutron star mergers
- Turbulent magnetic-field amplification in the first 10 milliseconds after a binary neutron star merger: comparing high-resolution and large eddy simulations
- Universality of the turbulent magnetic field in hypermassive neutron stars produced by binary mergers
- GRMHD large eddy simulations with gradient subgrid-scale model
- Formulating bulk viscosity for neutron star simulations
- Large Eddy Simulations of Magnetized Mergers of Neutron Stars with Neutrinos
- Impact of a mean field dynamo on neutron star mergers leading to magnetar remnants
- Artificial neural network subgrid models of 2-D compressible magnetohydrodynamic turbulence
- Magnetohydrodynamics with Physics Informed Neural Operators
- A covariant approach to relativistic large-eddy simulations: The fibration picture
- Electromagnetic Counterparts of Compact Binary Mergers
- A comparison of momentum transport models for numerical relativity
- Turbulence modelling in neutron star merger simulations
- A hybrid approach to long-term binary neutron-star simulations
- Assessment of a new sub-grid model for magneto-hydrodynamical turbulence. I. Magnetorotational instability
- Large Eddy Simulations of Magnetized Mergers of Black Holes and Neutron Stars
- Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations
- Assessment of a new sub-grid model for magnetohydrodynamical turbulence. II. Kelvin-Helmholtz instability
- Covariant approach to relativistic large-eddy simulations: Lagrangian filtering
- Higher-level large-eddy filtering strategy for general relativistic fluid simulations
- Local magneto-shear instability in Newtonian gravity
- Impact of magnetic field gradients on the development of the MRI: Applications to binary neutron star mergers and proto-planetary disks
- General relativistic magnetohydrodynamics simulations for binary neutron star mergers