Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical simulations of high Lundquist number relativistic magnetic reconnection
210
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
We use extensive 3D resistive MHD simulations to study how large-scale current sheets will undergo fast reconnection in the high Lundquist number $S$ limit (above $sim 10^4$), when the system is subject to different externally driven turbulence levels and the self-generated turbulence produced by 3D reconnection dynamics. We find that the normalized global reconnection rate $sim 0.01 - 0.13$, weakly dependent on $S$. Global reconnection with the classic inflow/outflow configurations is observed, and 3D flux ropes are hierarchically formed and ejected from reconnection regions. A statistical separation of the reconnected magnetic field lines follows a super-diffusive behavior, from which the rate is measured to be very similar to that obtained from the mixing of tracer populations. We find that the reconnection rate scales roughly linearly with the turbulence level during the peak of reconnection. This scaling is consistent with the turbulence properties produced by both the externally driven and self-generation processes. These results imply that large-scale thin current sheets tend to undergo rigorous reconnection.
A numerical study of magnetic reconnection in the large-Lundquist-number ($S$), plasmoid-dominated regime is carried out for $S$ up to $10^7$. The theoretical model of Uzdensky {it et al.} [Phys. Rev. Lett. {bf 105}, 235002 (2010)] is confirmed and partially amended. The normalized reconnection rate is $ ormEeffsim 0.02$ independently of $S$ for $Sgg10^4$. The plasmoid flux ($Psi$) and half-width ($w_x$) distribution functions scale as $f(Psi)sim Psi^{-2}$ and $f(w_x)sim w_x^{-2}$. The joint distribution of $Psi$ and $w_x$ shows that plasmoids populate a triangular region $w_xgtrsimPsi/B_0$, where $B_0$ is the reconnecting field. It is argued that this feature is due to plasmoid coalescence. Macroscopic monster plasmoids with $w_xsim 10$% of the system size are shown to emerge in just a few Alfven times, independently of $S$, suggesting that large disruptive events are an inevitable feature of large-$S$ reconnection.
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.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
