Energy distribution and substructure formation in astrophysical MHD simulations

Abstract

During substructure formation in magnetized astrophysical plasma, dissipation of magnetic energy facilitated by magnetic reconnection affects the system dynamics by heating and accelerating the ejected plasmoids. Numerical simulations are a crucial tool for investigating such systems. In astrophysical simulations, the energy dissipation, reconnection rate and substructure formation critically depend on the onset of reconnection of numerical or physical origin. In this paper, we hope to assess the reliability of the state-of-the-art numerical codes, PLUTO and KORAL by quantifying and discussing the impact of dimensionality, resolution, and code accuracy on magnetic energy dissipation, reconnection rate, and substructure formation. We quantitatively compare results obtained with relativistic and non-relativistic, resistive and non-resistive, as well as two- and three-dimensional setups performing the Orszag-Tang test problem. We find the sufficient resolution in each model, for which numerical error is negligible and the resolution does not significantly affect the magnetic energy dissipation and reconnection rate. The non-relativistic simulations show that at sufficient resolution, magnetic and kinetic energies convert to internal energy and heat up the plasma. The results show that in the relativistic system, energy components undergo mutual conversion during the simulation time, which leads to a substantial increase in magnetic energy at 20\% and 90\% of the total simulation time of $10$ light-crossing times -- the magnetic field is amplified by a factor of five due to relativistic shocks. We also show that the reconnection rate in all our simulations is higher than $0.1$, indicating plasmoid-mediated regime. It is shown that in KORAL simulations magnetic energy is slightly larger and more substructures are captured than in PLUTO simulations.

17 pages, 20 Figures,

References in corpus (19)

• PLUTO: a Numerical Code for Computational Astrophysics
• Instability of current sheets and formation of plasmoid chains
• Fast TeV variability in blazars: jets in a jet
• Black hole flares: ejection of accreted magnetic flux through 3D plasmoid-mediated reconnection
• Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes
• Global simulations of axisymmetric radiative black hole accretion disks in general relativity with a sub-grid magnetic dynamo
• A fundamental mechanism of solar eruption initiation
• Sgr A* near-infrared flares from reconnection events in a magnetically arrested disc
• On the diversity of magnetic interactions in close-in star-planet systems
• Small scale structures in three-dimensional magnetohydrodynamic turbulence
• Orbital motion near Sagittarius A* -- Constraints from polarimetric ALMA observations
• Tearing instability in relativistic magnetically dominated plasmas
• On magnetic field amplification and particle acceleration near non-relativistic astrophysical shocks: Particles in MHD Cells simulations
• Puffy accretion disks: sub-Eddington, optically thick, and stable
• Modelling MHD accretion-ejection - from the launching area to propagation scales
• Beamed emission from a neutron-star ULX in a GRRMHD simulation
• Forced magnetic reconnection and plasmoid coalescence: I - MHD Simulations
• Formation of episodic jets and associated flares from black hole accretion systems
• Magnetic reconnection in a comparison of topology and helicities in two and three dimensional resistive MHD simulations