A multidimensional numerical scheme for two-fluid relativistic magnetohydrodynamics

Barkov, Maxim V.; Komissarov, S. S.; Korolev, Vitaliy; Zankovich, Andrey M.

doi:10.1093/mnras/stt2247

Cited by 23 publications

(35 citation statements)

References 22 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The friction term R μ used in this paper differs from that used in earlier studies adopting the same set of equations (e.g., Zenitani et al 2009aZenitani et al , 2009bBarkov et al 2014). However, it is readily seen that our model reduces to the earlier work for a perfectly symmetric case (i.e., a pair plasma with ρ p =ρ e and…”

Section: Basic Equationsmentioning

confidence: 88%

“…Although application of the RTFED model to astrophysical problems so far has been very limited (Zenitani et al 2009a(Zenitani et al , 2009bAmano & Kirk 2013;Barkov et al 2014;Barkov & Komissarov 2016), we believe that it has the potential for more widespread use in the astrophysical community. Part of the reason is that numerical methods that can be used for the present system of equations have not adequately been explored.…”

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

A Second-Order Divergence-Constrained Multidimensional Numerical Scheme for Relativistic Two-Fluid Electrodynamics

Amano

2016

ApJ

View full text Add to dashboard Cite

A new multidimensional simulation code for relativistic two-fluid electrodynamics (RTFED) is described. The basic equations consist of the full set of Maxwell's equations coupled with relativistic hydrodynamic equations for separate two charged fluids, representing the dynamics of either an electron-positron or an electron-proton plasma. It can be recognized as an extension of conventional relativistic magnetohydrodynamics (RMHD). Finite resistivity may be introduced as a friction between the two species, which reduces to resistive RMHD in the long wavelength limit without suffering from a singularity at infinite conductivity. A numerical scheme based on HLL (HartenLax-Van Leer) Riemann solver is proposed that exactly preserves the two divergence constraints for Maxwell's equations simultaneously. Several benchmark problems demonstrate that it is capable of describing RMHD shocks/discontinuities at long wavelength limit, as well as dispersive characteristics due to the two-fluid effect appearing at small scales. This shows that the RTFED model is a promising tool for high energy astrophysics application.

show abstract

Section: Basic Equationsmentioning

confidence: 88%

Section: Introductionmentioning

confidence: 99%

A Second-Order Divergence-Constrained Multidimensional Numerical Scheme for Relativistic Two-Fluid Electrodynamics

Amano

2016

ApJ

View full text Add to dashboard Cite

show abstract

“…The conservation of particle number and of momentum for each species, and of total energy, then read as (see e.g., Mihalas & Mihalas 1984;Barkov et al 2014 These equations are, however, not well suited to analysis of particle simulations. First, because accessible quantities are those defined in the simulation (or lab) frame, while those in the comoving frame of the plasma must be obtained by a boost at the local mean speedū s .…”

Section: B1 Fluid Equationsmentioning

confidence: 99%

Relativistic magnetic reconnection in collisionless ion-electron plasmas explored with particle-in-cell simulations

Melzani

Walder

Folini

et al. 2014

A&A

View full text Add to dashboard Cite

Magnetic reconnection is a leading mechanism for magnetic energy conversion and high-energy non-thermal particle production in a variety of high-energy astrophysical objects, including ones with relativistic ion-electron plasmas (e.g., microquasars or AGNs), a regime where first principle studies are scarce. We present 2D particle-in-cell (PIC) simulations of low β ion-electron plasmas under relativistic conditions, i.e., with inflow magnetic energy exceeding the plasma restmass energy. We identify outstanding properties: (i) For relativistic inflow magnetizations (here 10 ≤ σ e ≤ 360), the reconnection outflows are dominated by thermal agitation instead of bulk kinetic energy. (ii) At high inflow electron magnetization (σ e ≥ 80), the reconnection electric field is sustained more by bulk inertia than by thermal inertia. It challenges the thermal-inertia paradigm and its implications. (iii) The inflows feature sharp transitions at the entrance of the diffusion zones. These are not shocks but results from particle ballistic motions, all bouncing at the same location, provided that the thermal velocity in the inflow is far lower than the inflow E × B bulk velocity. (iv) Island centers are magnetically isolated from the rest of the flow and can present a density depletion at their center. (v) The reconnection rates are slightly higher than in non-relativistic studies. They are best normalized by the inflow relativistic Alfvén speed projected in the outflow direction, which then leads to rates in a close range (0.14-0.25), thus allowing for an easy estimation of the reconnection electric field.

show abstract

“…An additional direction of research has been represented by the inclusion of dissipative effects, namely non-ideal resistive magnetohydrodynamics, with encouraging results (Komissarov 2007;Palenzuela et al 2009;Dumbser & Zanotti 2009;Zenitani et al 2010;Takamoto & Inoue 2011;Bucciantini & Del Zanna 2013). Moreover, high order numerical schemes have also been pursued (Del Zanna et al 2003;Anderson et al 2006), while simulations of multi-fluids in RMHD are emerging as a new frontier (Zenitani et al 2009;Barkov et al 2014). Finally, Adaptive Mesh Refinement (AMR) within RMHD codes has been also considered (Balsara 2001b;Neilsen et al 2006;Etienne et al 2010;Mignone et al 2012;Keppens et al 2012; and it is an active field of research.…”

Section: Introductionmentioning

confidence: 99%

Solving the relativistic magnetohydrodynamics equations with ADER discontinuous Galerkin methods, a posteriori subcell limiting and adaptive mesh refinement

Zanotti

Fambri

Dumbser

2015

Mon. Not. R. Astron. Soc.

114

View full text Add to dashboard Cite

We present a new numerical tool for solving the special relativistic ideal MHD equations that is based on the combination of the following three key features: (i) a one-step ADER discontinuous Galerkin (DG) scheme that allows for an arbitrary order of accuracy in both space and time, (ii) an a posteriori subcell finite volume limiter that is activated to avoid spurious oscillations at discontinuities without destroying the natural subcell resolution capabilities of the DG finite element framework and finally (iii) a space-time adaptive mesh refinement (AMR) framework with time-accurate local time-stepping.The divergence-free character of the magnetic field is instead taken into account through the so-called "divergence-cleaning" approach. The convergence of the new scheme is verified up to 5 th order in space and time and the results for a set of significant numerical tests including shock tube problems, the RMHD rotor and blast wave problems, as well as the Orszag-Tang vortex system are shown. We also consider a simple case of the relativistic Kelvin-Helmholtz instability with a magnetic field, emphasizing the potential of the new method for studying turbulent RMHD flows. We discuss the advantages of our new approach when the equations of relativistic MHD need to be solved with high accuracy within various astrophysical systems.

show abstract

A multidimensional numerical scheme for two-fluid relativistic magnetohydrodynamics

Cited by 23 publications

References 22 publications

A Second-Order Divergence-Constrained Multidimensional Numerical Scheme for Relativistic Two-Fluid Electrodynamics

A Second-Order Divergence-Constrained Multidimensional Numerical Scheme for Relativistic Two-Fluid Electrodynamics

Relativistic magnetic reconnection in collisionless ion-electron plasmas explored with particle-in-cell simulations

Solving the relativistic magnetohydrodynamics equations with ADER discontinuous Galerkin methods, a posteriori subcell limiting and adaptive mesh refinement

Contact Info

Product

Resources

About