Predicting Negative Triangularity Edge Confinement with a Neutral-Coupled Full-f Gyrokinetic Model
Abstract
Global full-f gyrokinetic simulations coupled to a reduced neutral-fluid model reproduce key DIII-D experimental trends and identify mechanisms underlying improved confinement in negative triangularity (NT) configurations relative to positive triangularity (PT). Using Gkeyll, we simulate DIII-D inner-wall-limited discharges across the confined edge, scrape-off layer (SOL), and sheath boundaries, allowing the edge electric field to evolve self-consistently. Compared with an adaptive recycling-source model, the reduced neutral fluid model improves particle and energy balances and brings the simulated density and temperature profiles into closer agreement with measurements, with minimal additional computational cost. In NT, suppressed turbulent fluxes coincide with a broad, stationary zonal flow inside the last closed flux surface, whereas PT develops transient, fragmented shear layers. Spectral analysis indicates that modified fluctuation phase relationships and enhanced zonal shear both contribute to transport suppression despite broadly comparable fluctuation amplitudes. The larger co-current departure from a locally matched bi-Maxwellian in NT suggests a velocity-selective kinetic process that may help establish the edge radial electric field by modifying flows or current balance. These results identify kinetic and turbulent mechanisms associated with improved NT confinement within an experimentally established ELM-free regime and demonstrate the importance of coupling gyrokinetic plasma dynamics with neutral interactions across the confined plasma and SOL.
Related
- Validation of geometric proxies using negative and positive triangularity plasmas with varying elongation on TCV and multi-machine gyrokinetic simulations
- Options for heating, fueling, and current drive in MANTA, a high-field negative triangularity fusion pilot plant
- Compact Experimental Negative TriAngUlarity Reactor (CENTAUR): A design study for a compact, affordable breakeven tokamak
- Plasma-neutral collision operators for hydrogen isotopes: kinetic and gyrokinetic descriptions, comparison to fluid
- Plasma-neutral collision operators for hydrogen isotopes: kinetic and gyrokinetic descriptions, comparison to fluid
BibTeX
@misc{bernard2026predicting,
title = {{Predicting Negative Triangularity Edge Confinement with a Neutral-Coupled Full-f Gyrokinetic Model}},
author = {T. N. Bernard and A. C. D Hoffmann and M. Francisquez and F. D. Halpern and G. W. Hammett and A. Marinoni and C. Chrystal and G. McKee},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.32100},
eprint = {2610.32100},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.32100},
}
RIS
TY - UNPB
TI - Predicting Negative Triangularity Edge Confinement with a Neutral-Coupled Full-f Gyrokinetic Model
AU - T. N. Bernard
AU - A. C. D Hoffmann
AU - M. Francisquez
AU - F. D. Halpern
AU - G. W. Hammett
AU - A. Marinoni
AU - C. Chrystal
AU - G. McKee
PY - 2026
DA - 2026/10/23
PB - arXiv
DO - 10.48550/arxiv.2610.32100
UR - https://arxiv.org/abs/2610.32100
AB - Global full-f gyrokinetic simulations coupled to a reduced neutral-fluid model reproduce key DIII-D experimental trends and identify mechanisms underlying improved confinement in negative triangularity (NT) configurations relative to positive triangularity (PT). Using Gkeyll, we simulate DIII-D inner-wall-limited discharges across the confined edge, scrape-off layer (SOL), and sheath boundaries, allowing the edge electric field to evolve self-consistently. Compared with an adaptive recycling-source model, the reduced neutral fluid model improves particle and energy balances and brings the simulated density and temperature profiles into closer agreement with measurements, with minimal additional computational cost. In NT, suppressed turbulent fluxes coincide with a broad, stationary zonal flow inside the last closed flux surface, whereas PT develops transient, fragmented shear layers. Spectral analysis indicates that modified fluctuation phase relationships and enhanced zonal shear both contribute to transport suppression despite broadly comparable fluctuation amplitudes. The larger co-current departure from a locally matched bi-Maxwellian in NT suggests a velocity-selective kinetic process that may help establish the edge radial electric field by modifying flows or current balance. These results identify kinetic and turbulent mechanisms associated with improved NT confinement within an experimentally established ELM-free regime and demonstrate the importance of coupling gyrokinetic plasma dynamics with neutral interactions across the confined plasma and SOL.
ER -