Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical simulations of relativistic magnetic reconnection with Galerkin methods
208
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present the results of two-dimensional magnetohydrodynamical numerical simulations of relativistic magnetic reconnection, with particular emphasis on the dynamics of Petschek-type configurations with high Lundquist numbers, S ~ 10^5-10^8. The numerical scheme adopted, allowing for unprecedented accuracy for this type of calculations, is based on high order finite volume and discontinuous Galerkin methods as recently proposed by Dumbser & Zanotti (2009). The possibility of producing high Lorentz factors is discussed, by studying the effects produced on the dynamics by different magnetization and resistivity regimes. We show that Lorentz factors close to ~4 can be produced for a plasma magnetization parameter sigma=20. Moreover, we find that the Sweet-Parker layers are unstable, generating secondary magnetic islands, but only for S>S_c~10^8, much larger than what is reported in the Newtonian regime.
rate research
Read More
We present the results of two-dimensional and three-dimensional magnetohydrodynamical numerical simulations of relativistic magnetic reconnection, with particular emphasis on the dynamics of the plasma in a Petschek-type configuration with high Lundquist numbers, Ssim 10^5-10^8. The numerical scheme adopted, allowing for unprecedented accuracy for this type of calculations, is based on high order finite volume and discontinuous Galerkin methods as recently proposed by citet{Dumbser2009}. The possibility of producing high Lorentz factors is discussed, showing that Lorentz factors close to sim 4 can be produced for a plasma parameter sigma_m=20. Moreover, we find that the Sweet-Parker layers are unstable, generating secondary magnetic islands, but only for S > S_c = 10^8, much larger than what is reported in the Newtonian regime. Finally, the effects of a mildly anisotropic Ohm law are considered in a configuration with a guide magnetic field. Such effects produce only slightly faster reconnection rates and Lorentz factors of about 1% larger with respect to the perfectly isotropic Ohm law.
First results are presented from kinetic numerical simulations of relativistic collisionless magnetic reconnection in pair plasma that include radiation reaction from both synchrotron and inverse Compton (IC) processes, motivated by non-thermal high-energy astrophysical sources, including in particular blazars. These simulations are initiated from a configuration known as ABC fields that evolves due to coalescence instability and generates thin current layers in its linear phase. Global radiative efficiencies, instability growth rates, time-dependent radiation spectra, lightcurves, variability statistics and the structure of current layers are investigated for a broad range of initial parameters. We find that the IC radiative signatures are generally similar to the synchrotron signatures. The luminosity ratio of IC to synchrotron spectral components, the Compton dominance, can be modified by more than one order of magnitude with respect to its nominal value. For very short cooling lengths, we find evidence for modification of the temperature profile across the current layers, no systematic compression of plasma density, and very consistent profiles of E.B. We decompose the profiles of E.B with the use of the Vlasov momentum equation, demonstrating a contribution from radiation reaction at the thickness scale consistent with the temperature profile.
We present the results of 2D particle-in-cell (PIC) simulations of relativistic magnetic reconnection (RMR) in electron-positron plasma, including the dynamical influence of the synchrotron radiation process, and integrating the observable emission signatures. The simulations are initiated with a single Harris current layer with a central gap that triggers the RMR process. We achieve a steady-state reconnection with unrestricted outflows by means of open boundary conditions. The radiative cooling efficiency is regulated by the choice of initial plasma temperature Theta. We explore different values of Theta and of the background magnetisation sigma_0. Throughout the simulations, plasmoids are generated in the central region of the layer, and they evolve at different rates, achieving a wide range of sizes. The gaps between plasmoids are filled by smooth relativistic outflows called minijets, whose contribution to the observed radiation is very limited due to their low particle densities. Small-sized plasmoids are rapidly accelerated, however, they have lower contributions to the observed emission, despite stronger relativistic beaming. Large-sized plasmoids are slow, but produce most of the observed synchrotron emission, with major part of their radiation produced within the central cores, the density of which is enhanced by radiative cooling. Synchrotron lightcurves show rapid bright flares that can be identified as originating from mergers between small/fast plasmoids and large/slow targets moving in the same direction. In the high-magnetisation case, the accelerated particles form a broken power-law energy distribution with a soft tail produced by particles accelerated in the minijets.
We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27-solar-mass progenitor was studied in 2D by B. Mueller et al. (ApJ 761:72, 2012), who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.
Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection from the non-relativistic to ultra-relativistic limit. In the anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter $sigma$ and approaches the light speed when $sigma > O(100)$, leading to an enhanced reconnection rate. In all regimes, the divergence of pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region remains $sim 0.1$ in both the non-relativistic and relativistic limits.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
