56th Annual Meeting of the APS Division of Plasma Physics
Volume 59, Number 15
Monday–Friday, October 27–31, 2014;
New Orleans, Louisiana
Session BI1: Edge Pedestal Physics
9:30 AM–12:30 PM,
Monday, October 27, 2014
Room: Acadia
Chair: Anne White, Massachusetts Institute of Technology
Abstract ID: BAPS.2014.DPP.BI1.5
Abstract: BI1.00005 : 4D Fokker-Planck calculations of neoclassical effects in tokamak pedestals and stellarators*
11:30 AM–12:00 PM
Preview Abstract
Abstract
Author:
Matt Landreman
(University of Maryland)
In this work, techniques are discussed for efficiently solving time-independent 4D (2 space + 2 velocity) linear drift-kinetic equations with Fokker-Planck collisions, illustrating several applications in stellarators and tokamaks. Conventional calculations of neoclassical flows, current, and fluxes can become inadequate in the tokamak pedestal if strong gradients violate the thin-orbit-width ordering, yet accurate calculation of these quantities is important for understanding edge stability and confinement. We have therefore developed a new radially global continuum neoclassical code PERFECT [1] which allows some radial scale lengths to be as small as the poloidal ion gyroradius. Using this tool, we demonstrate the first precise verification of recent analytic theory that accounts for finite orbit width effects in the limit of large aspect ratio [2]. The main-ion and impurity flows predicted by the finite-orbit-width calculations can have significantly different magnitude and poloidal variation compared to conventional (local) calculations. PERFECT solves 4D drift-kinetic equations for each species, using a preconditioned Krylov-space algorithm to accelerate solution of the time-independent problem. A related new code SFINCS implements these same algorithms but with a toroidal angle in place of the radial coordinate, enabling radially local calculations for stellarators and rippled tokamaks [3]. For W7-X and LHD scenarios, we use SFINCS to compare the common incompressible-$E\times B$ drift-kinetic equation to more accurate drift-kinetic equations. For radial electric fields below roughly 1/3 of the resonant value, the different kinetic equations lead to similar predictions for the fluxes and flows, whereas the results can differ significantly for larger electric fields. In SFINCS calculations for expected W7-X conditions, the various kinetic models predict similar levels of bootstrap current, improving confidence in past predictions for this quantity that strongly impacts the W7-X divertor.\\[4pt]
[1] Landreman et al, PPCF 56, 045005 (2014).\\[0pt]
[2] Catto et al, PPCF 55, 045009 (2013).\\[0pt]
[3] Landreman et al, Phys Plasmas 21, 042503 (2014).
*Supported by US DOE.
To cite this abstract, use the following reference: http://meetings.aps.org/link/BAPS.2014.DPP.BI1.5