Impact of velocity space distribution on hybrid kinetic-magnetohydrodynamic simulation of the (1,1) mode

Abstract

Numeric studies of the impact of the velocity space distribution on the stabilization of (1,1) internal kink mode and excitation of the fishbone mode are performed with a hybrid kinetic-magnetohydrodynamic model. These simulations demonstrate an extension of the physics capabilities of NIMROD [C. R. Sovinec et al., J. Comput. Phys. 195, 355 (2004)], a three-dimensional extended magnetohydrodynamic (MHD) code, to include the kinetic effects of an energetic minority ion species. Kinetic effects are captured by a modification of the usual MHD momentum equation to include a pressure tensor calculated from the δf particle-in-cell method [S. E. Parker and W. W. Lee, Phys. Fluids B 5, 77 (1993)]. The particles are advanced in the self-consistent NIMROD fields. We outline the implementation and present simulation results of energetic minority ion stabilization of the (1,1) internal kink mode and excitation of the fishbone mode. A benchmark of the linear growth rate and real frequency is shown to agree well with another code. The impact of the details of the velocity space distribution is examined; particularly extending the velocity space cutoff of the simulation particles. Modestly increasing the cutoff strongly impacts the (1,1) mode. Numeric experiments are performed to study the impact of passing versus trapped particles. Observations of these numeric experiments suggest that assumptions of energetic particle effects should be re-examined.

A low-frequency variational model for energetic particle effects in the pressure-coupling scheme

“…However, many variants have been proposed and a direct comparison is not always clear. For example, in the works [31,32,47] the authors use the definition of pressure involving the relative velocity, that is…”

Time-domain simulation of nonlinear radiofrequency phenomena

“…Such methods find use in the study of plasma wakefield accelerators, 25 astrophysical problems, 26 and propulsion, 27 to name only a few applications. In magnetic fusion research, they have been used to investigate the behavior of ion temperature gradient-driven turbulence 28 and the kinetic stabilization of internal kink modes 29 in tokamaks, as well as to model wave-particle interactions in the plasma core. 30 In the most rudimentary form (the "total-f" method) of PIC, one writes the collisionless kinetic equation for the distribution function f s…”

A model of energetic ion effects on pressure driven tearing modes in tokamaks

“…where (27) and i ∈ [1,2]. The perturbed magnetic flux solutionψ has the standard cylindrical r ±m form, and can be written as a linear combination of basis functions (in radial order)…”