0000000000115210
AUTHOR
Tomasz Rembiasz
On the Measurements of Numerical Viscosity and Resistivity in Eulerian MHD Codes
We propose a simple ansatz for estimating the value of the numerical resistivity and the numerical viscosity of any Eulerian MHD code. We test this ansatz with the help of simulations of the propagation of (magneto)sonic waves, Alfven waves, and the tearing mode instability using the MHD code Aenus. By comparing the simu- lation results with analytical solutions of the resistive-viscous MHD equations and an empirical ansatz for the growth rate of tearing modes we measure the numerical viscosity and resistivity of Aenus. The comparison shows that the fast-magnetosonic speed and wavelength are the characteristic velocity and length, respectively, of the aforementioned (relatively simple) syst…
Termination of the magnetorotational instability via parasitic instabilities in core-collapse supernovae
The magnetorotational instability (MRI) can be a powerful mechanism amplifying the magnetic field in core-collapse supernovae. Whether initially weak magnetic fields can be amplified by this instability to dynamically relevant strengths is still a matter of debate. One of the main uncertainties concerns the process that terminates the growth of the instability. Parasitic instabilities of both Kelvin-Helmholtz and tearing-mode type have been suggested to play a crucial role in this process, disrupting MRI channel flows and quenching magnetic field amplification. We perform two-dimensional and three-dimensional sheering-disc simulations of a differentially rotating protoneutron star layer in …
Employee shirking and overworking: modelling the unintended consequences of work organisation
Underworking (i.e. shirking) and overworking of employees can have detrimental effects for the individual and the organisation. We develop a computational model to investigate how work structure, specifically the way in which managers distribute work tasks amongst employees, impacts work intensity and working time. The model draws on theories from economics, psychology and management, and on empirical observations. The simulations show that when managers correctly estimate task difficulty, but undervalue the employee’s competence, opportunities for shirking are provided due to longer deadlines. Similarly, if managers overvalue the employee’s competence, they set tighter deadlines leading to…
How to form a millisecond magnetar? Magnetic field amplification in protoneutron stars
Extremely strong magnetic fields of the order of $10^{15}\,{\rm G}$ are required to explain the properties of magnetars, the most magnetic neutron stars. Such a strong magnetic field is expected to play an important role for the dynamics of core-collapse supernovae, and in the presence of rapid rotation may power superluminous supernovae and hypernovae associated to long gamma-ray bursts. The origin of these strong magnetic fields remains, however, obscure and most likely requires an amplification over many orders of magnitude in the protoneutron star. One of the most promising agents is the magnetorotational instability (MRI), which can in principle amplify exponentially fast a weak initia…
Magnetorotational Instability in Core-Collapse Supernovae
We discuss the relevance of the magnetorotational instability (MRI) in core-collapse supernovae (CCSNe). Our recent numerical studies show that in CCSNe, the MRI is terminated by parasitic instabilities of the Kelvin-Helmholtz type. To determine whether the MRI can amplify initially weak magnetic fields to dynamically relevant strengths in CCSNe, we performed three-dimensional simulations of a region close to the surface of a differentially rotating proto-neutron star in non-ideal magnetohydrodynamics with two different numerical codes. We find that under the conditions prevailing in proto-neutron stars, the MRI can amplify the magnetic field by (only) one order of magnitude. This severely …
Heavy sterile neutrinos in stellar core-collapse
We perform spherically symmetric simulations of the core collapse of a single progenitor star of zero age main sequence mass $M_{\rm ZAMS} = 15 \, M_{\odot}$ with two models of heavy sterile neutrinos in the mass range of hundred MeV$/c^2$. According to both models, these hypothetical particles are copiously produced in the center, stream outwards a subsequently decay releasing energy into final states (including neutrinos) of the Standard Model. We find that they can lead to a successful explosion in otherwise non-exploding progenitors. Depending on their unknown parameters (e.g., mass and coupling constants with matter), we obtain either no explosion or an explosion of one of two types, i…
High-order methods for the simulation of hydromagnetic instabilities in core-collapse supernovae
AbstractWe present an assessment of the accuracy of a recently developed MHD code used to study hydromagnetic flows in supernovae and related events. The code, based on the constrained transport formulation, incorporates unprecedented ultra-high-order methods (up to 9th order) for the reconstruction and the most accurate approximate Riemann solvers. We estimate the numerical resistivity of these schemes in tearing instability simulations.
On the maximum magnetic field amplification by the magnetorotational instability in core-collapse supernovae
Whether the magnetorotational instability (MRI) can amplify initially weak magnetic fields to dynamically relevant strengths in core collapse supernovae is still a matter of active scientific debate. Recent numerical studies have shown that the first phase of MRI growth dominated by channel flows is terminated by parasitic instabilities of the Kelvin-Helmholtz type that disrupt MRI channel flows and quench further magnetic field growth. However, it remains to be prop- erly assessed by what factor the initial magnetic field can be amplified and how it depends on the initial field strength and the amplitude of the perturbations. Different termination criteria leading to different estimates of…