Viriato : A Fourier–Hermite spectral code for strongly magnetized fluid–kinetic plasma dynamics
We report on the algorithms and numerical methods used in Viriato, a novel fluid–kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations (Zocco and Schekochihin, 2011) (which reduce to the standard Reduced-MHD equations in the appropr...
| Main Authors: | , , , , , , |
|---|---|
| Other Authors: | |
| Format: | Article |
| Published: |
Elsevier BV
2018
|
| Online Access: | http://hdl.handle.net/1721.1/118403 https://orcid.org/0000-0001-9755-6563 |
| _version_ | 1826202746693353472 |
|---|---|
| author | Dorland, W. Fazendeiro, L. Kanekar, A. Mallet, A. Vilelas, M.S. Zocco, A. Gomes Loureiro, Nuno F |
| author2 | Massachusetts Institute of Technology. Department of Nuclear Science and Engineering |
| author_facet | Massachusetts Institute of Technology. Department of Nuclear Science and Engineering Dorland, W. Fazendeiro, L. Kanekar, A. Mallet, A. Vilelas, M.S. Zocco, A. Gomes Loureiro, Nuno F |
| author_sort | Dorland, W. |
| collection | MIT |
| description | We report on the algorithms and numerical methods used in Viriato, a novel fluid–kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations (Zocco and Schekochihin, 2011) (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations (Schekochihin et al., 2009). Two main applications of these equations are magnetized (Alfvénic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge–Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor–corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag–Tang-type decaying turbulence, both in fluid and kinetic regimes. Keywords: PlasmaFourier–Hermite, Reduced gyrokinetics, Turbulence, Magnetic reconnection |
| first_indexed | 2024-09-23T12:17:05Z |
| format | Article |
| id | mit-1721.1/118403 |
| institution | Massachusetts Institute of Technology |
| last_indexed | 2024-09-23T12:17:05Z |
| publishDate | 2018 |
| publisher | Elsevier BV |
| record_format | dspace |
| spelling | mit-1721.1/1184032022-10-01T08:53:48Z Viriato : A Fourier–Hermite spectral code for strongly magnetized fluid–kinetic plasma dynamics Dorland, W. Fazendeiro, L. Kanekar, A. Mallet, A. Vilelas, M.S. Zocco, A. Gomes Loureiro, Nuno F Massachusetts Institute of Technology. Department of Nuclear Science and Engineering Massachusetts Institute of Technology. Plasma Science and Fusion Center Gomes Loureiro, Nuno F We report on the algorithms and numerical methods used in Viriato, a novel fluid–kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations (Zocco and Schekochihin, 2011) (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations (Schekochihin et al., 2009). Two main applications of these equations are magnetized (Alfvénic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge–Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor–corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag–Tang-type decaying turbulence, both in fluid and kinetic regimes. Keywords: PlasmaFourier–Hermite, Reduced gyrokinetics, Turbulence, Magnetic reconnection Fundação para a Ciência e a Tecnologia (Portugal) (UID/FIS/50010/2013) Fundação para a Ciência e a Tecnologia (Portugal) (PTDC/FIS/118187/2010) Fundação para a Ciência e a Tecnologia (Portugal) (IF/00530/2013) Leverhulme Trust 2018-10-09T19:15:09Z 2018-10-09T19:15:09Z 2016-09 2016-04 2018-09-25T18:01:34Z Article http://purl.org/eprint/type/JournalArticle 0010-4655 http://hdl.handle.net/1721.1/118403 Loureiro, N. F., et al. “Viriato : A Fourier–Hermite Spectral Code for Strongly Magnetized Fluid–Kinetic Plasma Dynamics.” Computer Physics Communications, vol. 206, Sept. 2016, pp. 45–63. https://orcid.org/0000-0001-9755-6563 http://dx.doi.org/10.1016/J.CPC.2016.05.004 Computer Physics Communications Creative Commons Attribution-NonCommercial-NoDerivs License http://creativecommons.org/licenses/by-nc-nd/4.0/ application/pdf Elsevier BV arXiv |
| spellingShingle | Dorland, W. Fazendeiro, L. Kanekar, A. Mallet, A. Vilelas, M.S. Zocco, A. Gomes Loureiro, Nuno F Viriato : A Fourier–Hermite spectral code for strongly magnetized fluid–kinetic plasma dynamics |
| title | Viriato : A Fourier–Hermite spectral code for strongly magnetized fluid–kinetic plasma dynamics |
| title_full | Viriato : A Fourier–Hermite spectral code for strongly magnetized fluid–kinetic plasma dynamics |
| title_fullStr | Viriato : A Fourier–Hermite spectral code for strongly magnetized fluid–kinetic plasma dynamics |
| title_full_unstemmed | Viriato : A Fourier–Hermite spectral code for strongly magnetized fluid–kinetic plasma dynamics |
| title_short | Viriato : A Fourier–Hermite spectral code for strongly magnetized fluid–kinetic plasma dynamics |
| title_sort | viriato a fourier hermite spectral code for strongly magnetized fluid kinetic plasma dynamics |
| url | http://hdl.handle.net/1721.1/118403 https://orcid.org/0000-0001-9755-6563 |
| work_keys_str_mv | AT dorlandw viriatoafourierhermitespectralcodeforstronglymagnetizedfluidkineticplasmadynamics AT fazendeirol viriatoafourierhermitespectralcodeforstronglymagnetizedfluidkineticplasmadynamics AT kanekara viriatoafourierhermitespectralcodeforstronglymagnetizedfluidkineticplasmadynamics AT malleta viriatoafourierhermitespectralcodeforstronglymagnetizedfluidkineticplasmadynamics AT vilelasms viriatoafourierhermitespectralcodeforstronglymagnetizedfluidkineticplasmadynamics AT zoccoa viriatoafourierhermitespectralcodeforstronglymagnetizedfluidkineticplasmadynamics AT gomesloureironunof viriatoafourierhermitespectralcodeforstronglymagnetizedfluidkineticplasmadynamics |