ATTENZIONE! IL REPOSITORY E’ APERTO ESCLUSIVAMENTE ALL’IMMISSIONE DEI RAPPORTI TECNICI (COLLEZIONE 4.01). NON E’ POSSIBILE INSERIRE PRODOTTI IN TUTTE LE ALTRE COLLEZIONI, A CAUSA DI PROBLEMI TECNICI. THE REPOSITORY IS OPEN EXCLUSIVELY FOR THE SUBMISSION OF TECHNICAL REPORTS (COLLECTION 4.01). SUBMISSIONS TO ALL OTHER COLLECTIONS ARE TEMPORARILY UNAVAILABLE, DUE TO TECHNICAL ISSUES.
 

A Particle Module for the PLUTO Code. I. An Implementation of the MHD-PIC Equations

Journal
Date Issued
2018
Author(s)
Mignone, A.
•
•
Vaidya, B.
•
Mattia, G.
DOI
Abstract
We describe an implementation of a particle physics module available for the PLUTO code appropriate for the dynamical evolution of a plasma consisting of a thermal fluid and a nonthermal component represented by relativistic charged particles or cosmic rays (CRs). While the fluid is approached using standard numerical schemes for magnetohydrodynamics, CR particles are treated kinetically using conventional Particle-In-Cell (PIC) techniques. The module can be used either to describe test-particle motion in the fluid electromagnetic field or to solve the fully coupled magnetohydrodynamics (MHD)-PIC system of equations with particle backreaction on the fluid as originally introduced by Bai et al. Particle backreaction on the fluid is included in the form of momentum-energy feedback and by introducing the CR-induced Hall term in Ohm’s law. The hybrid MHD-PIC module can be employed to study CR kinetic effects on scales larger than the (ion) skin depth provided that the Larmor gyration scale is properly resolved. When applicable, this formulation avoids resolving microscopic scales, offering substantial computational savings with respect to PIC simulations. We present a fully conservative formulation that is second-order accurate in time and space, and extends to either the Runge-Kutta (RK) or the corner transport upwind time-stepping schemes (for the fluid), while a standard Boris integrator is employed for the particles. For highly energetic relativistic CRs and in order to overcome the time-step restriction, a novel subcycling strategy that retains second-order accuracy in time is presented. Numerical benchmarks and applications including Bell instability, diffusive shock acceleration, and test-particle acceleration in reconnecting layers are discussed.
Volume
859
Issue
1
Start page
13
Uri
Url
Issn Identifier
0004-637X
Ads BibCode
2018ApJ...859...13M
Rights
open.access
File(s)
Loading...
Thumbnail Image
Name

Mignone_2018_ApJ_859_13.pdf

Description
Pdf editoriale
Size

2.64 MB

Format

Adobe PDF

Checksum (MD5)

373c6a95fe99bb3008e9fdfc039ae5fd