Abstract:The high-fidelity BOUT++ two-fluid code suite has demonstrated significant recent progress toward integrated multi-scale simulations of tokamak pedestal, including Edge-Localized-Mode (ELM) dynamics, evolution of ELM cycles, and continuous fluctuations, as observed in experiments. Nonlinear ELM simulations show three stages of an ELM event: (1) a linear growing phase; (2) a fast crash phase; and (3) a slow inward turbulence spreading phase lasting until the core heating flux balances the ELM energy loss and the ELM is terminated. A new coupling/splitting model has been developed to perform simulations of multi-scale ELM dynamics. Simulation tracks five ELM cycles for 10 000 Alfvén times for small ELMs. The temporal evolution of the pedestal pressure is similar to that of experimental measurements for the pedestal pressure profile collapses and recovers to a steep gradient during ELM cycles. To validate BOUT++ simulations against experimental data and develop physics understanding of the fluctuation characteristics for different tokamak operation regimes, both quasi-coherent fluctuations (QCFs) in ELMy H-modes and Weakly Coherent Modes in I-modes have been simulated using three dimensional 6-field 2-fluid electromagnetic model. The H-mode simulation results show that (1) QCFs are localized in the pedestal region having a predominant frequency at f≃300−400 kHz and poloidal wavenumber at kθ≃0.7 cm−1, and propagate in the electron diamagnetic direction in the laboratory frame. The overall signatures of simulation results for QCFs show good agreement with C-Mod and DIII-D measurements. (2) The pedestal profiles giving rise to QCFs are near the marginal instability threshold for ideal peeling-ballooning modes for both C-Mod and DIII-D, while the collisional electromagnetic drift-Alfvén wave appears to be dominant for DIII-D. (3) Particle diffusivity is either smaller than the heat diffusivity for DIII-D or similar to the heat diffusivity for C-Mod. Key I-mode simulation results are that (1) a strong instability exists at n≥20 for resistive ballooning mode and drift-Alfvén wave; (2) the frequency spectrum of nonlinear BOUT++ simulation features a peak around 300 kHz for the mode number n = 20, consistent with a reflectometer measurement at nearby position; (3) the calculated particle diffusivity is larger than the calculated heat diffusivity, which is consistent with a key feature of the I-mode pedestal with no particle barrier.

Kinetic-Ballooning-Bifurcation in Tokamak Pedestals Across Shaping and Aspect-Ratio

Kinetic-Ballooning-Limited Pedestals in Spherical Tokamak Plasmas

Microturbulence in Diii-D Tokamak Pedestal. Iii. Effects of Collisions

Stability and transport of gyrokinetic critical pedestals

Sensitivity of Kinetic Ballooning Mode Instability to Tokamak Equilibrium Implementations

Microturbulence in Diii-D Tokamak Pedestal. I. Electrostatic Instabilities

Global Gyrokinetic Simulation of Tokamak Edge Pedestal Instabilities

TH/P2-32 Gyrokinetic Simulation of Tokamak Edge Plasmas

Parallel Ion Compressibility Effects on Kinetic Ballooning Mode for Different Magnetic Shears

Electromagnetic gyrokinetic instabilities in STEP

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

Toward Integrated Multi-Scale Pedestal Simulations Including Edge-Localized-Mode Dynamics, Evolution of Edge-Localized-Mode Cycles, and Continuous Fluctuations

Broadening of the Divertor Heat Flux Profile in High Confinement Tokamak Fusion Plasmas with Edge Pedestals Limited by Turbulence in DIII-D

A fast low-to-high mode bifurcation dynamics in a tokamak edge plasma gyrokinetic simulation

The Effect of Triangularity on Predicted Pedestals for the CFETR

The Shaping Effect of Cross Section on Ballooning Modes in Tokamak Plasmas

Linear gyrokinetic stability of a high $β$ non-inductive spherical tokamak

Unconventional Ballooning Structures for Toroidal Drift Waves

Impact of triangularity on edge peeling–ballooning modes in H-mode plasmas