Continuum electromagnetic gyrokinetic simulations of turbulence in the tokamak scrape-off layer and laboratory devices
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Princeton Univ., NJ (United States)
- General Atomics, San Diego, CA (United States); Univ. of Texas, Austin, TX (United States)
- MIT Plasma Science and Fusion Center, Cambridge, MA (United States); Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Here, we present algorithms and results from Gkeyll, a full-f continuum, electromagnetic gyrokinetic code, designed to study turbulence in the edge region of fusion devices. The edge is computationally very challenging, requiring robust algorithms that can handle large-amplitude fluctuations and stable interactions with plasma sheaths. We present an energy-conserving high-order discontinuous Galerkin scheme that solves gyrokinetic equations in Hamiltonian form. Efficiency is improved by a careful choice of basis functions and automatically generated computation kernels. Previous verification tests were performed in the straight-field-line large plasma device [Shi et al., J. Plasma Phys. 83, 905830304 (2017)] and the Texas Helimak, a simple magnetized torus [Bernard et al., Phys. Plasmas 26, 042301 (2019)], including the effect of end-plate biasing on turbulence. Results for the scrape-off layer for NSTX parameters with a model helical magnetic geometry with bad curvature have been obtained [Shi et al., Phys. Plasmas 26, 012307 (2019)]. In this paper, we present algorithms for the two formulations of electromagnetic gyrokinetics: the Hamiltonian and the symplectic. We describe each formulation and show results of benchmark tests. Although our scheme works for the Hamiltonian formulation, the presence of spurious numerical modes for high-β and large $$k^2_⊥ρ^2_s$$ regimes shows that the symplectic formulation is more robust. We then review our recent algorithm for the symplectic formulation [Mandell et al., J. Plasma Phys. 86, 905860109 (2020)], along with example application of this new capability. Maintaining positivity of the distribution function can be challenging, and we describe a new and novel exponential recovery based algorithm to address this.
- Research Organization:
- Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- USDOE
- Grant/Contract Number:
- FG02-97ER25308; FG02-95ER54309; FG02-04ER-54742; AC02-09CH11466; FC02-08ER54966
- OSTI ID:
- 1668050
- Alternate ID(s):
- OSTI ID: 1608504
- Journal Information:
- Physics of Plasmas, Vol. 27, Issue 4; ISSN 1070-664X
- Publisher:
- American Institute of Physics (AIP)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
Investigating shear flow through continuum gyrokinetic simulations of limiter biasing in the Texas Helimak
Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems