Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Univ. of Science and Technology of China, Hefei (China)
The hybrid method combining particle-in-cell and magnetohydrodynamics can be used to study the interaction between energetic particles and global plasma modes. In this paper we introduce the M3D-C1-K code, which is developed based on the M3D-C1 finite element code solving the magnetohydrodynamics equations, with a newly developed kinetic module simulating energetic particles. The particle pushing is done using a new algorithm by applying the Boris pusher to the classical Pauli particles to simulate the slow-manifold of particle orbits, with long-term accuracy and fidelity. The particle pushing can be accelerated using GPUs with a significant speedup. The moments of the particles are calculated using the δƒ method, and are coupled into the magnetohydrodynamics simulation through pressure or current coupling schemes. Several linear simulations of magnetohydrodynamics modes driven by energetic particles have been conducted using M3D-C1-K with the δƒ method, including fishbone, toroidal Alfvén eigenmodes and reversed shear Alfvén eigenmodes. Good agreement with previous results from other eigenvalue, kinetic and hybrid codes have been achieved.
- Research Organization:
- Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- USDOE
- Grant/Contract Number:
- AC02-09CH11466; AC05-00OR22725
- OSTI ID:
- 1869943
- Journal Information:
- Computer Physics Communications, Journal Name: Computer Physics Communications Vol. 275; ISSN 0010-4655
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Hybrid simulations of fishbone instabilities and Alfvén eigenmodes in DIII-D tokamak
Three-dimensional hybrid gyrokinetic-magnetohydrodynamics simulation