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

Title: Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

Abstract

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited duemore » to the stability properties of the linear least squares projective integrator.« less

Authors:
ORCiD logo [1];  [2]; ORCiD logo [1]; ORCiD logo [1]
  1. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
  2. Univ. of Colorado, Boulder, CO (United States). Dept. of Physics
Publication Date:
Research Org.:
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Sponsoring Org.:
USDOE
Contributing Org.:
17-SC-20-SC - exascale computing project
OSTI Identifier:
1606354
Alternate Identifier(s):
OSTI ID: 1604235
Grant/Contract Number:  
AC02-09CH11466; 17-SC-20-SC
Resource Type:
Accepted Manuscript
Journal Name:
Physics of Plasmas
Additional Journal Information:
Journal Volume: 27; Journal Issue: 3; Journal ID: ISSN 1070-664X
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Gyrokinetic simulations; Particle-in-cell method; Tokamaks; Multiscale methods

Citation Formats

Sturdevant, Benjamin J., Parker, Scott E., Chang, C. S., and Hager, Robert. Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry. United States: N. p., 2020. Web. doi:10.1063/1.5141337.
Sturdevant, Benjamin J., Parker, Scott E., Chang, C. S., & Hager, Robert. Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry. United States. doi:https://doi.org/10.1063/1.5141337
Sturdevant, Benjamin J., Parker, Scott E., Chang, C. S., and Hager, Robert. Sun . "Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry". United States. doi:https://doi.org/10.1063/1.5141337. https://www.osti.gov/servlets/purl/1606354.
@article{osti_1606354,
title = {Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry},
author = {Sturdevant, Benjamin J. and Parker, Scott E. and Chang, C. S. and Hager, Robert},
abstractNote = {A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.},
doi = {10.1063/1.5141337},
journal = {Physics of Plasmas},
number = 3,
volume = 27,
place = {United States},
year = {2020},
month = {3}
}

Works referenced in this record:

Plateau regime dynamics of the relaxation of poloidal rotation in tokamak plasmas
journal, August 1996

  • Lebedev, V. B.; Yushmanov, P. N.; Diamond, P. H.
  • Physics of Plasmas, Vol. 3, Issue 8
  • DOI: 10.1063/1.871638

Equation Free Projective Integration: A multiscale method applied to a plasma ion acoustic wave
journal, September 2007


A fast low-to-high confinement mode bifurcation dynamics in the boundary-plasma gyrokinetic code XGC1
journal, May 2018

  • Ku, S.; Chang, C. S.; Hager, R.
  • Physics of Plasmas, Vol. 25, Issue 5
  • DOI: 10.1063/1.5020792

Variational principles of guiding centre motion
journal, February 1983


Direct multiscale coupling of a transport code to gyrokinetic turbulence codes
journal, May 2010

  • Barnes, M.; Abel, I. G.; Dorland, W.
  • Physics of Plasmas, Vol. 17, Issue 5
  • DOI: 10.1063/1.3323082

Tokamak profile prediction using direct gyrokinetic and neoclassical simulation
journal, June 2009

  • Candy, J.; Holland, C.; Waltz, R. E.
  • Physics of Plasmas, Vol. 16, Issue 6
  • DOI: 10.1063/1.3167820

Coarse projective kMC integration: forward/reverse initial and boundary value problems
journal, May 2004

  • Rico-Martı́nez, R.; Gear, C. W.; Kevrekidis, Ioannis G.
  • Journal of Computational Physics, Vol. 196, Issue 2
  • DOI: 10.1016/j.jcp.2003.11.005

ComPat framework for multiscale simulations applied to fusion plasmas
journal, June 2019


Information Theory and Statistical Mechanics. II
journal, October 1957


Information Theory and Statistical Mechanics
journal, May 1957


Projective Methods for Stiff Differential Equations: Problems with Gaps in Their Eigenvalue Spectrum
journal, January 2003


Using the maximum entropy distribution to describe electrons in reconnecting current sheets
journal, August 2018

  • Ng, Jonathan; Hakim, Ammar; Bhattacharjee, A.
  • Physics of Plasmas, Vol. 25, Issue 8
  • DOI: 10.1063/1.5041758

Elimination of fast variables
journal, July 1985


The moment-guided Monte Carlo method
journal, May 2010

  • Degond, Pierre; Dimarco, Giacomo; Pareschi, Lorenzo
  • International Journal for Numerical Methods in Fluids, Vol. 67, Issue 2
  • DOI: 10.1002/fld.2345

Bringing global gyrokinetic turbulence simulations to the transport timescale using a multiscale approach
journal, March 2018


Numerical study of tokamak equilibria with arbitrary flow
journal, February 2004

  • Guazzotto, L.; Betti, R.; Manickam, J.
  • Physics of Plasmas, Vol. 11, Issue 2
  • DOI: 10.1063/1.1637918

Kinetic theory of low-frequency Alfvén modes in tokamaks
journal, November 1996

  • Zonca, Fulvio; Chen, Liu; Santoro, Robert A.
  • Plasma Physics and Controlled Fusion, Vol. 38, Issue 11
  • DOI: 10.1088/0741-3335/38/11/011

Gyrokinetic approach in particle simulation
journal, January 1983


Stability of general plasma equilibria - I formal theory
journal, January 1968


Nonlinear gyrokinetic equations for low-frequency electromagnetic waves in general plasma equilibria
journal, January 1982


A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma
journal, June 2016


Gyrokinetic particle simulation model
journal, September 1987


Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles
journal, January 2012

  • Crestetto, Anaïs; Crouseilles, Nicolas; Lemou, Mohammed
  • Kinetic & Related Models, Vol. 5, Issue 4
  • DOI: 10.3934/krm.2012.5.787

Self-consistent modeling of turbulence and transport
journal, March 2003


Equation-Free Multiscale Computation: Algorithms and Applications
journal, May 2009


Equation-Free, Coarse-Grained Multiscale Computation: Enabling Mocroscopic Simulators to Perform System-Level Analysis
journal, January 2003

  • Gear, C. William; Hyman, James M.; Kevrekidid, Panagiotis G.
  • Communications in Mathematical Sciences, Vol. 1, Issue 4
  • DOI: 10.4310/CMS.2003.v1.n4.a5

Numerical Linear Algebra
book, January 1997


Drift Instabilities in General Magnetic Field Configurations
journal, January 1968


Transport of momentum in full f gyrokinetics
journal, May 2010

  • Parra, Felix I.; Catto, Peter J.
  • Physics of Plasmas, Vol. 17, Issue 5
  • DOI: 10.1063/1.3327127