BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Darin R. Ernst appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species interactions are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, satisfies the H-Theorem unlike analogous Bhatnagar-Gross-Krook operators.
Numerical techniques for discretization of velocity space in continuum kinetic calculations are described. An efficient spectral collocation method is developed for the speed coordinate - the radius in velocity space - employing a novel set of non-classical orthogonal polynomials. For problems in which Fokker-Planck collisions are included, a common situation in plasma physics, a procedure is detailed to accurately and efficiently treat the field term in the collision operator (in the absence of gyrokinetic corrections). When species with disparate masses are included simultaneously, a careful extrapolation of the Rosenbluth potentials is performed. The techniques are demonstrated in neoclassical calculations of the bootstrap current and plasma flows in a tokamak.
In a tokamak pedestal, radial scale lengths can become comparable to the ion orbit width, invalidating conventional neoclassical calculations of flow and bootstrap current. In this work we illustrate a non-local approach that allows strong radial density variation while maintaining small departures from a Maxwellian distribution. Non-local effects alter the magnitude and poloidal variation of the flow and current. The approach is implemented in a new global delta-f continuum code using the full linearized Fokker-Planck collision operator. Arbitrary collisionality and aspect ratio are allowed as long as the poloidal magnetic field is small compared to the total magnetic field. Strong radial electric fields, sufficient to electrostatically confine the ions, are also included. These effects may be important to consider in any comparison between experimental pedestal flow measurements and theory.
In transport barriers, particularly H-mode edge pedestals, radial scale lengths can become comparable to the ion orbit width, causing neoclassical physics to become radially nonlocal. In this work, the resulting changes to neoclassical flow and current are examined both analytically and numerically. Steep density gradients are considered, with scale lengths comparable to the poloidal ion gyroradius, together with strong radial electric fields sufficient to electrostatically confine the ions. Attention is restricted to relatively weak ion temperature gradients (but permitting arbitrary electron temperature gradients), since in this limit a delta-f (small departures from a Maxwellian distribution) rather than full-f approach is justified. This assumption is in fact consistent with measured inter-ELM H-Mode edge pedestal density and ion temperature profiles in many present experiments, and is expected to be increasingly valid in future lower collisionality experiments. In the numerical analysis, the distribution function and Rosenbluth potentials are solved for simultaneously, allowing use of the exact field term in the linearized Fokker-Planck collision operator. In the pedestal, the parallel and poloidal flows are found to deviate strongly from the best available conventional neoclassical prediction, with large poloidal variation of a different form than in the local theory. These predicted effects may be observable experimentally. In the local limit, the Sauter bootstrap current formulae appear accurate at low collisionality, but they can overestimate the bootstrap current near the plateau regime. In the pedestal ordering, ion contributions to the bootstrap and Pfirsch-Schluter currents are also modified.