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

Title: Global linear gyrokinetic simulation of energetic particle-driven instabilities in the LHD stellarator

Abstract

Energetic particles are inherent to toroidal fusion systems and can drive instabilities in the Alfvén frequency range, leading to decreased heating efficiency, high heat fluxes on plasma-facing components, and decreased ignition margin. The applicability of global gyrokinetic simulation methods to macroscopic instabilities has now been demonstrated and it is natural to extend these methods to 3D configurations such as stellarators, tokamaks with 3D coils and reversed field pinch helical states. This has been achieved by coupling the GTC global gyrokinetic PIC model to the VMEC equilibrium model, including 3D effects in the field solvers and particle push. Here, this paper demonstrates the application of this new capability to the linearized analysis of Alfvénic instabilities in the LHD stellarator. For normal shear iota profiles, toroidal Alfvén instabilities in the n = 1 and 2 toroidal mode families are unstable with frequencies in the 75 to 110 kHz range. Also, an LHD case with non-monotonic shear is considered, indicating reductions in growth rate for the same energetic particle drive. Finally, since 3D magnetic fields will be present to some extent in all fusion devices, the extension of gyrokinetic models to 3D configurations is an important step for the simulation of future fusionmore » systems.« less

Authors:
ORCiD logo [1];  [2];  [3];  [3]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  2. Univ. of California, Irvine, CA (United States). Dept. of Physics and Astronomy; Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  3. National Inst. for Fusion Science, Toki (Japan)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1366391
Grant/Contract Number:
AC05-00OR22725; AC02-05CH11231
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Nuclear Fusion
Additional Journal Information:
Journal Volume: 57; Journal Issue: 8; Journal ID: ISSN 0029-5515
Publisher:
IOP Science
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; stellarator; LHD; gyrokinetics; Alfven wave; energetic particles; fast ion instability linear

Citation Formats

Spong, Donald A., Holod, Ihor, Todo, Y., and Osakabe, M. Global linear gyrokinetic simulation of energetic particle-driven instabilities in the LHD stellarator. United States: N. p., 2017. Web. doi:10.1088/1741-4326/aa7601.
Spong, Donald A., Holod, Ihor, Todo, Y., & Osakabe, M. Global linear gyrokinetic simulation of energetic particle-driven instabilities in the LHD stellarator. United States. doi:10.1088/1741-4326/aa7601.
Spong, Donald A., Holod, Ihor, Todo, Y., and Osakabe, M. 2017. "Global linear gyrokinetic simulation of energetic particle-driven instabilities in the LHD stellarator". United States. doi:10.1088/1741-4326/aa7601.
@article{osti_1366391,
title = {Global linear gyrokinetic simulation of energetic particle-driven instabilities in the LHD stellarator},
author = {Spong, Donald A. and Holod, Ihor and Todo, Y. and Osakabe, M.},
abstractNote = {Energetic particles are inherent to toroidal fusion systems and can drive instabilities in the Alfvén frequency range, leading to decreased heating efficiency, high heat fluxes on plasma-facing components, and decreased ignition margin. The applicability of global gyrokinetic simulation methods to macroscopic instabilities has now been demonstrated and it is natural to extend these methods to 3D configurations such as stellarators, tokamaks with 3D coils and reversed field pinch helical states. This has been achieved by coupling the GTC global gyrokinetic PIC model to the VMEC equilibrium model, including 3D effects in the field solvers and particle push. Here, this paper demonstrates the application of this new capability to the linearized analysis of Alfvénic instabilities in the LHD stellarator. For normal shear iota profiles, toroidal Alfvén instabilities in the n = 1 and 2 toroidal mode families are unstable with frequencies in the 75 to 110 kHz range. Also, an LHD case with non-monotonic shear is considered, indicating reductions in growth rate for the same energetic particle drive. Finally, since 3D magnetic fields will be present to some extent in all fusion devices, the extension of gyrokinetic models to 3D configurations is an important step for the simulation of future fusion systems.},
doi = {10.1088/1741-4326/aa7601},
journal = {Nuclear Fusion},
number = 8,
volume = 57,
place = {United States},
year = 2017,
month = 6
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on June 23, 2018
Publisher's Version of Record

Citation Metrics:
Cited by: 1work
Citation information provided by
Web of Science

Save / Share:
  • Comprehensive understanding of energetic-ion-driven global instabilities such as Alfven eigenmodes (AEs) and their impact on energetic ions and bulk plasma is crucially important for tokamak and stellarator/helical plasmas and in the future for deuterium-tritium (DT) burning plasma experiments. Various types of global modes and their associated enhanced energetic ion transport are commonly observed in toroidal plasmas. Toroidicity-induced AEs and ellipticity-induced AEs, whose gaps are generated through poloidal mode coupling, are observed in both tokamak and stellarator/helical plasmas. Global AEs and reversed shear AEs, where toroidal couplings are not as dominant were also observed in those plasmas. Helicity induced AEs thatmore » exist only in 3D plasmas are observed in the large helical device (LHD) and Wendelstein 7 Advanced Stellarator plasmas. In addition, the geodesic acoustic mode that comes from plasma compressibility is destabilized by energetic ions in both tokamak and LHD plasmas. Nonlinear interaction of these modes and their influence on the confinement of the bulk plasma as well as energetic ions are observed in both plasmas. In this paper, the similarities and differences in these instabilities and their consequences for tokamak and stellarator/helical plasmas are summarized through comparison with the data sets obtained in LHD. In particular, this paper focuses on the differences caused by the rotational transform profile and the 2D or 3D geometrical structure of the plasma equilibrium. Important issues left for future study are listed.« less
  • The unstable spectrum of Alfven eigenmodes (AEs) driven by neutral beam-sourced energetic particles (EPs) in a benchmark DIII-D discharge (142111) is calculated in a fully gyrokinetic model using the GYRO code's massively parallel linear eigenvalue solver. One cycle of the slow (equilibrium scale) frequency sweep of the reverse shear Alfven eigenmode (RSAE) at toroidal mode number n=3 is mapped. The RSAE second harmonic and an unstable beta-induced Alfven eigenmode (BAE) are simultaneously tracked alongside the primary RSAE. An observed twist in the eigenmode pattern, caused mostly by shear in the driving EP profile, is shown through artificially varying the Emore » Multiplication-Sign B rotational velocity shear to depend generally on shear in the local wave phase velocity. Coupling to the BAE and to the toroidal Alfven eigenmode limit the RSAE frequency sweeps at the lower and upper end, respectively. While the present fully gyrokinetic model (including thermal ions and electrons) constitutes the best treatment of compressibility physics available, the BAE frequency is overpredicted by about 20% against experiment here and is found to be sensitive to energetic beam ion pressure. The RSAE frequency is more accurately matched except when it is limited by the BAE. Simulations suggest that the experiment is very close to marginal AE stability at points of RSAE-BAE coupling. A recipe for comparing the radial profile of quasilinear transport flux from local modes to that from global modes paves the way for the development of a stiff (critical gradient) local AE transport model based on local mode stability thresholds.« less
  • No abstract prepared.
  • No abstract prepared.