Particle pinch in global tokamak edge simulations
Abstract
The inward particle flux, or particle pinch, is routinely observed in magnetically confined fusion experiments, yet its mechanism is not fully understood. We study the particle pinch in the tokamak edge with the flux-driven global turbulence model GDB. Starting from a flat density profile fueled only near the last closed flux surface, the simulation develops a strong inward particle flux that builds a centrally peaked density profile over O(10) ms. The pinch coexists with the usual outward turbulent heat transport; two distinct mechanisms carry it. In the early stage, when density and temperature gradients oppose each other (ηα=Ln/LTα<0), drift-wave turbulence drives the inward flux through electron thermal diffusion. In the late stage, once the density profile has flattened, a persistent inward equilibrium E× B flux, carried by the poloidally asymmetric density and electrostatic potential, drives the central peaking. While the Pfirsch--Schlüter neoclassical transport sets the amplitude of the up-down asymmetric density, classical theory predicts no net radial flux from this asymmetry at leading order. The observed flux flows in a separate non-ideal channel, opened by parallel resistivity, electron inertia, and electromagnetic induction in the equilibrium electron force balance. This channel shifts the equilibrium potential poloidally against the density by δs≃−0.08π and supplies the late-phase density build-up. These results show that a flux-driven global edge simulation can self-consistently produce a centrally peaked density profile without ad hoc assumptions, and bear on longstanding edge questions such as density pedestal formation.
Related
- Spontaneous Enhancement of Inward Impurity Pinch during Massive Impurity Injection in Tokamak Plasmas
- Two-fluid boundary turbulence simulations in reversed field pinch plasmas
- Modeling of plasma transport during edge-localized mode in tokamak using a kinetic Vlasov-Poisson code
- Implicit-Explicit time integration scheme with Physics-based preconditioning for two-fluid tokamak boundary simulations
- Impact of magnetic drift configuration on the edge radial electric field in the TCV tokamak
BibTeX
@misc{zhu2026particle,
title = {{Particle pinch in global tokamak edge simulations}},
author = {Ben Zhu},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.06744},
eprint = {2609.06744},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.06744},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Particle pinch in global tokamak edge simulations
AU - Ben Zhu
PY - 2026
DA - 2026/09/06
PB - arXiv
DO - 10.48550/arxiv.2609.06744
UR - https://arxiv.org/abs/2609.06744
AB - The inward particle flux, or particle pinch, is routinely observed in magnetically confined fusion experiments, yet its mechanism is not fully understood. We study the particle pinch in the tokamak edge with the flux-driven global turbulence model GDB. Starting from a flat density profile fueled only near the last closed flux surface, the simulation develops a strong inward particle flux that builds a centrally peaked density profile over O(10) ms. The pinch coexists with the usual outward turbulent heat transport; two distinct mechanisms carry it. In the early stage, when density and temperature gradients oppose each other (η_α=L_n/L_T_α<0), drift-wave turbulence drives the inward flux through electron thermal diffusion. In the late stage, once the density profile has flattened, a persistent inward equilibrium E× B flux, carried by the poloidally asymmetric density and electrostatic potential, drives the central peaking. While the Pfirsch--Schlüter neoclassical transport sets the amplitude of the up-down asymmetric density, classical theory predicts no net radial flux from this asymmetry at leading order. The observed flux flows in a separate non-ideal channel, opened by parallel resistivity, electron inertia, and electromagnetic induction in the equilibrium electron force balance. This channel shifts the equilibrium potential poloidally against the density by δ_s≃−0.08π and supplies the late-phase density build-up. These results show that a flux-driven global edge simulation can self-consistently produce a centrally peaked density profile without ad hoc assumptions, and bear on longstanding edge questions such as density pedestal formation.
ER -