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

Title: A conservative scheme of drift kinetic electrons for gyrokinetic simulation of kinetic-MHD processes in toroidal plasmas

Authors:
 [1]; ORCiD logo [2]; ORCiD logo [3]
  1. Department of Physics and Astronomy, University of California, Irvine, California 92697, USA, Fusion Simulation Center, Peking University, Beijing 100871, China
  2. College of Physical Science and Technology, Sichuan University, Chengdu 610064, China
  3. Department of Physics and Astronomy, University of California, Irvine, California 92697, USA
Publication Date:
Sponsoring Org.:
USDOE
OSTI Identifier:
1402100
Grant/Contract Number:
AC02-05CH11231; AC05-00OR22725
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Physics of Plasmas
Additional Journal Information:
Journal Volume: 24; Journal Issue: 10; Related Information: CHORUS Timestamp: 2018-02-14 13:13:55; Journal ID: ISSN 1070-664X
Publisher:
American Institute of Physics
Country of Publication:
United States
Language:
English

Citation Formats

Bao, J., Liu, D., and Lin, Z. A conservative scheme of drift kinetic electrons for gyrokinetic simulation of kinetic-MHD processes in toroidal plasmas. United States: N. p., 2017. Web. doi:10.1063/1.4995455.
Bao, J., Liu, D., & Lin, Z. A conservative scheme of drift kinetic electrons for gyrokinetic simulation of kinetic-MHD processes in toroidal plasmas. United States. doi:10.1063/1.4995455.
Bao, J., Liu, D., and Lin, Z. 2017. "A conservative scheme of drift kinetic electrons for gyrokinetic simulation of kinetic-MHD processes in toroidal plasmas". United States. doi:10.1063/1.4995455.
@article{osti_1402100,
title = {A conservative scheme of drift kinetic electrons for gyrokinetic simulation of kinetic-MHD processes in toroidal plasmas},
author = {Bao, J. and Liu, D. and Lin, Z.},
abstractNote = {},
doi = {10.1063/1.4995455},
journal = {Physics of Plasmas},
number = 10,
volume = 24,
place = {United States},
year = 2017,
month =
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on October 23, 2018
Publisher's Accepted Manuscript

Save / Share:
  • A gyrokinetic-magnetohydrodynamic (MHD) hybrid simulation code has been developed in order to study high-{ital n} (where {ital n} is the toroidal mode number) MHD instabilities driven by energetic particles in finite-{beta} tokamaks. Here {beta} is the ratio between plasma and magnetic pressures. Specifically, it is observed that as the core plasma {beta} value increases, there is a corresponding transition from the toroidal Alfv{acute e}n eigenmode (TAE) to the kinetic ballooning mode (KBM). The energetic particle mode (EPM) branches of both the toroidal Alfv{acute e}n mode (TAM) and KBM are shown to be important in this transition. KBM are preferentially excitedmore » when the energetic particle velocity is small compared to the Alfv{acute e}n velocity. {copyright} {ital 1996 American Institute of Physics.}« less
  • The eigenmode stability properties of three-dimensional lower-hybrid-drift-instabilities (LHDI) in a Harris current sheet with a small but finite guide magnetic field have been systematically studied by employing the gyrokinetic electron and fully kinetic ion (GeFi) particle-in-cell (PIC) simulation model with a realistic ion-to-electron mass ratio m i/m e. In contrast to the fully kinetic PIC simulation scheme, the fast electron cyclotron motion and plasma oscillations are systematically removed in the GeFi model, and hence one can employ the realistic m i/m e. The GeFi simulations are benchmarked against and show excellent agreement with both the fully kinetic PIC simulation and the analytical eigenmode theory. Our studies indicate that, for small wavenumbers, ky, along the current direction, the most unstable eigenmodes are peaked at the location wheremore » $$\vec{k}$$• $$\vec{B}$$ =0, consistent with previous analytical and simulation studies. Here, $$\vec{B}$$ is the equilibrium magnetic field and $$\vec{k}$$ is the wavevector perpendicular to the nonuniformity direction. As ky increases, however, the most unstable eigenmodes are found to be peaked at $$\vec{k}$$ •$$\vec{B}$$ ≠0. Additionally, the simulation results indicate that varying m i/m e, the current sheet width, and the guide magnetic field can affect the stability of LHDI. Simulations with the varying mass ratio confirm the lower hybrid frequency and wave number scalings.« less
  • Transport processes and resultant entropy production in magnetically confined plasmas are studied in detail for toroidal systems with gyrokinetic electromagnetic turbulence. The kinetic equation including the turbulent fluctuations are double averaged over the ensemble and the gyrophase. The entropy balance equation is derived from the double-averaged kinetic equation with the nonlinear gyrokinetic equation for the fluctuating distribution function. The result clarifies the spatial transport and local production of the entropy due to the classical, neoclassical and anomalous transport processes, respectively. For the anomalous transport process due to the electromagnetic turbulence as well as the classical and neoclassical processes, the kineticmore » form of the entropy production is rewritten as the thermodynamic form, from which the conjugate pairs of the thermodynamic forces and the transport fluxes are identified. The Onsager symmetry for the anomalous transport equations is shown to be valid within the quasilinear framework. The complete energy balance equation, which takes account of the anomalous transport and exchange of energy due to the fluctuations, is derived from the ensemble-averaged kinetic equation. The intrinsic ambipolarity of the anomalous particle fluxes is shown to hold for the self-consistent turbulent electromagnetic fields satisfying Poisson{close_quote}s equation and Amp{grave e}re{close_quote}s law. {copyright} {ital 1996 American Institute of Physics.}« less