HYPIC: A fast hybrid EM PIC-MCC code for ion cyclotron resonance energization in cylindrical coordinate system

Ion cyclotron resonance energization (ICRE) such as ion cyclotron resonance heating (ICRH) is widely applied to magnetic confinement fusion and high-power electric propulsion. Since ICRE involves cyclotron resonance processes, a kinetic model is required. Both conventional particle-in-cell (PIC) simulations and solving the Boltzmann equation require enormous computation and memory. The hybrid simulation incorporating of adiabatic electrons and PIC ions provides a viable solution for both a substantial reduction in computation and the inclusion of cyclotron resonance effects. Under the adiabatic electron approximation, we have developed a two-dimensional (r, z) hybrid electromagnetic (EM) PIC-MCC (Monte-Carlo collision) simulation program, named HYPIC. The advantages of HYPIC are the inclusion of ion kinetic effects, electrostatic (ES) and EM effects, and collisional effects of ions and electrons, with a small computation. The HYPIC program is able to fast simulate the antenna-plasma interactions and the ion cyclotron resonance energization and/or ion cyclotron resonance heating processes in linear devices, such as high-power electric propulsion, magnetic mirror, and fieldreversed-configuration (FRC), etc. Program Summary: Program title: HYPIC. CPC Library link to program files: https://doi.org/10.17632/f6k2mx2tj9.1. Licensing provisions: BSD 3-clause. Programming language: Fortran. Nature of problem: The code solves the electromagnetic field distribution in wave-plasma interactions and the motion of ions in the presence of electrostatic, electromagnetic, background magnetic fields, and collisions in two-dimensional (2D) (r, z) cylindrical coordinate system. The code is able to fast simulate the antenna-plasma interactions, ion cyclotron resonance energization and/or ion cyclotron resonance heating processes in linear devices, such as high-power electric propulsion, magnetic mirror, and field-reversed-configuration, etc. Solution method: The PIC method with adiabatic electrons is used to fast simulate of the motion of ions by 4-stage Runge-Kutta time integral. The frequency-domain finite-difference method to solve Maxwell's equations allows for a fast solution of radio-frequency wave-plasma interactions. Monte Carlo collisions are used to deal with collisions of electrons and ions. Additional comments including restrictions and unusual features: The current version only supports solving axisymmetric 2D problems.