Plasmas with varying collisionalities occur in many applications, such as tokamak edge regions, where the flows are characterized by significant variations in density and temperature. While a kinetic model is necessary for weakly-collisional high-temperature plasmas, high collisionality in colder regions render the equations numerically stiff due to disparate time scales. In this paper, we propose an implicitexplicit algorithm for such cases, where the collisional term is integrated implicitly in time, while the advective term is integrated explicitly in time, thus allowing time step sizes that are comparable to the advective time scales. This partitioning results in a more efficient algorithm than those using explicit time integrators, where the -implicit additive Runge-Kutta methods in COGENT, a finite-volume gyrokinetic code for mapped, multiblock grids and test the accuracy, convergence, and computational cost of these semi-implicit methods for test cases with highly-collisional plasmas.Keywords IMEX time integration · plasma physics · gyrokinetic simulations · Vlasov-Fokker-Planck equations
IntroductionThe purpose of this paper is to describe the application and performance of a semiimplicit time integration algorithm for the solution of a system of Vlasov-Fokker-Planck equations, motivated by the goal of simulating the edge plasma region of tokamak fusion reactors. Plasma dynamics in the tokamak edge region is an unsteady multiscale phenomenon, characterized by a large range of spatial and temporal scales due to the density and temperature variations. Figure 1 shows the cross-section of a typical tokamak fusion reactor. The geometry is defined by the magnetic flux surfaces that contain the plasma, and the core and the edge regions are marked. Within the edge region, as the temperature decreases from the hot near-core region to the cold outeredge region, there are three scale regimes. The hot and dense plasma in the inner edge region adjacent to the core plasma is weakly collisional, and the mean free paths of the particles are significantly larger than the density and temperature gradient length scales [20,26,25]. Near the separatrix that separates the closed magnetic field lines from the open ones, the plasma is moderately collisional, and the mean free paths are comparable to the density and temperature gradient scales. At the outer edge (near the material surfaces), the cold plasma is strongly collisional. In this region, the particle mean free paths are significantly smaller than the density and temperature length scales.Due to the weak collisionality near the core, a kinetic description with an appropriate collision model is required to model accurately the perturbations to the velocity distribution from the Maxwellian distribution. The Vlasov-Fokker-Planck (VFP) equation governs the evolution of the distribution function of each charged particle in the position and velocity space [80]. In the presence of a strong, externally-applied magnetic field, the ionized particles gyrate around the magnetic fiel...