skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Differential formulation of the gyrokinetic Landau operator

Abstract

Subsequent to the recent rigorous derivation of an energetically consistent gyrokinetic collision operator in the so-called Landau representation, this work investigates the possibility of finding a differential formulation of the gyrokinetic Landau collision operator. It is observed that, while a differential formulation is possible in the gyrokinetic phase space, reduction of the resulting system of partial differential equations to five dimensions via gyroaveraging poses a challenge. Finally, based on the present work, it is likely that the gyrocentre analogues of the Rosenbluth–MacDonald–Judd potential functions must be kept gyroangle dependent.

Authors:
ORCiD logo [1];  [2]; ORCiD logo [1]
  1. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
  2. Saint Michael's College, Colchester, VT (United States). Dept. of Physics
Publication Date:
Research Org.:
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Saint Michael's College, Colchester, VT (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Fusion Energy Sciences (FES) (SC-24)
OSTI Identifier:
1351535
Report Number(s):
PPPL-5340
Journal ID: ISSN 0022-3778; PII: S0022377816001203
Grant/Contract Number:
AC02-09CH11466; SC0014032
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Plasma Physics
Additional Journal Information:
Journal Volume: 83; Journal Issue: 1; Journal ID: ISSN 0022-3778
Publisher:
Cambridge University Press
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; fusion plasma; plasma nonlinear phenomena; plasma simulation

Citation Formats

Hirvijoki, Eero, Brizard, Alain J., and Pfefferlé, David. Differential formulation of the gyrokinetic Landau operator. United States: N. p., 2017. Web. doi:10.1017/s0022377816001203.
Hirvijoki, Eero, Brizard, Alain J., & Pfefferlé, David. Differential formulation of the gyrokinetic Landau operator. United States. doi:10.1017/s0022377816001203.
Hirvijoki, Eero, Brizard, Alain J., and Pfefferlé, David. Thu . "Differential formulation of the gyrokinetic Landau operator". United States. doi:10.1017/s0022377816001203. https://www.osti.gov/servlets/purl/1351535.
@article{osti_1351535,
title = {Differential formulation of the gyrokinetic Landau operator},
author = {Hirvijoki, Eero and Brizard, Alain J. and Pfefferlé, David},
abstractNote = {Subsequent to the recent rigorous derivation of an energetically consistent gyrokinetic collision operator in the so-called Landau representation, this work investigates the possibility of finding a differential formulation of the gyrokinetic Landau collision operator. It is observed that, while a differential formulation is possible in the gyrokinetic phase space, reduction of the resulting system of partial differential equations to five dimensions via gyroaveraging poses a challenge. Finally, based on the present work, it is likely that the gyrocentre analogues of the Rosenbluth–MacDonald–Judd potential functions must be kept gyroangle dependent.},
doi = {10.1017/s0022377816001203},
journal = {Journal of Plasma Physics},
number = 1,
volume = 83,
place = {United States},
year = {Thu Jan 05 00:00:00 EST 2017},
month = {Thu Jan 05 00:00:00 EST 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:
  • Achieving plasmas with good stability and confinement properties is a key research goal for magnetic fusion devices. The underlying equations are the Vlasov–Poisson and Vlasov–Maxwell (VPM) equations in three space variables, three velocity variables, and one time variable. Even in those somewhat academic cases where global equilibrium solutions are known, studying their stability requires the analysis of the spectral properties of the linearized operator, a daunting task. We have identified a model, for which not only equilibrium solutions can be constructed, but many of their stability properties are amenable to rigorous analysis. It uses a class of solution to themore » VPM equations (or to their gyrokinetic approximations) known as waterbag solutions which, in particular, are piecewise constant in phase-space. It also uses, not only the gyrokinetic approximation of fast cyclotronic motion around magnetic field lines, but also an asymptotic approximation regarding the magnetic-field-induced anisotropy: the spatial variation along the field lines is taken much slower than across them. Together, these assumptions result in a drastic reduction in the dimensionality of the linearized problem, which becomes a set of two nested one-dimensional problems: an integral equation in the poloidal variable, followed by a one-dimensional complex Schrödinger equation in the radial variable. We show here that the operator associated to the poloidal variable is meromorphic in the eigenparameter, the pulsation frequency. We also prove that, for all but a countable set of real pulsation frequencies, the operator is compact and thus behaves mostly as a finite-dimensional one. The numerical algorithms based on such ideas have been implemented in a companion paper [D. Coulette and N. Besse, “Numerical resolution of the global eigenvalue problem for gyrokinetic-waterbag model in toroidal geometry” (submitted)] and were found to be surprisingly close to those for the original gyrokinetic-Vlasov equations. The purpose of the present paper is to make these new ideas accessible to two readerships: applied mathematicians and plasma physicists.« less
  • Fluidlike models have long been used to develop qualitative understanding of the drift-wave class of instabilities (such as the ion temperature gradient mode and various trapped-particle modes) which are prime candidates for explaining anomalous transport in plasmas. Here, the fluid approach is improved by developing fairly realistic models of kinetic effects, such as Landau damping and gyroradius orbit averaging, which strongly affect both the linear mode properties and the resulting nonlinear turbulence. Central to this work is a simple but effective fluid model (Phys. Rev. Lett. {bold 64}, 3019 (1990)) of the collisionless phase mixing responsible for Landau damping (andmore » inverse Landau damping). This model is based on a nonlocal damping term with a damping rate {similar to} {ital v}{sub {ital t}}{vert bar}{ital k}{sub {parallel}}{vert bar} in the closure approximation for the {ital n}th velocity space moment of the distribution function {ital f}, resulting in an {ital n}-pole approximation of the plasma dispersion function {ital Z}. Alternatively, this closure approximation is linearly exact (and therefore physically realizable) for a particular {ital f}{sub 0} which is close to Maxwellian. Gyrofluid'' equations (conservation laws for the {ital guiding}-{ital center} density {ital n}, momentum {ital mnu}{sub {parallel}}, and parallel and perpendicular pressures {ital p}{sub {parallel}} and {ital p}{sub {perpendicular}}) are derived by taking moments of the gyrokinetic equation in guiding-center coordinates rather than particle coordinates. This naturally yields nonlinear gyroradius terms and an important gyroaveraging of the shear. The gyroradius effects in the Bessel functions are modeled with robust Pade-like approximations.« less
  • This paper presents a class of collision operators, which reproduce neoclassical transport and comply with the constraints of a full-f global gyrokinetic code. The assessment of these operators is based on a variational entropy method, which allows a fast calculation of the neoclassical diffusivity and poloidal velocity.
  • Conventional neoclassical predictions are successfully recovered within a gyrokinetic framework using a minimal Fokker-Planck collision operator. This operator is shown to accurately describe some essential features of neoclassical theory, namely the neoclassical transport, the poloidal rotation and the linear damping of axisymmetric flows while interestingly preserving a high numerical efficiency. Its form makes it especially adapted to Eulerian or Semi-Lagrangian schemes.