Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Journal Article · Thu Nov 20 00:00:00 EST 2025 · Physics of Plasmas

DOI:https://doi.org/10.1063/5.0294040· OSTI ID:3006374

Liu, D. ^[1]; Juno, J. ^[2]; Hammett, G. W. ^[2]; Hakim, A. ^[2]; Shukla, A. ^[3]; Francisquez, M. ^[2]

---

1. Princeton Univ., NJ (United States)
2. Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
3. Univ. of Texas, Austin, TX (United States)

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

• Research Organization: Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
• Sponsoring Organization: National Energy Research Scientific Computing Center (NERSC); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR). Scientific Discovery through Advanced Computing (SciDAC); USDOE Office of Science (SC), Basic Energy Sciences (BES). Scientific User Facilities (SUF)
• Grant/Contract Number: AC02-09CH11466
• OSTI ID: 3006374
• Journal Information: Physics of Plasmas, Journal Name: Physics of Plasmas Journal Issue: 11 Vol. 32; ISSN 1070-664X; ISSN 1089-7674
• Publisher: American Institute of Physics (AIP)
• Country of Publication: United States
• Language: English

Related Subjects

Gyrokinetic simulations
Plasma collisions
Plasma confinement
Plasma transport properties
Tokamaks