Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Univ. of Colorado, Boulder, CO (United States). Dept. of Physics
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.
- Research Organization:
- Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- USDOE
- Contributing Organization:
- 17-SC-20-SC - exascale computing project
- Grant/Contract Number:
- AC02-09CH11466
- OSTI ID:
- 1606354
- Alternate ID(s):
- OSTI ID: 1604235
- Journal Information:
- Physics of Plasmas, Journal Name: Physics of Plasmas Journal Issue: 3 Vol. 27; ISSN 1070-664X
- Publisher:
- American Institute of Physics (AIP)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Tokamak profile prediction using direct gyrokinetic and neoclassical simulation
Tokamak profile prediction using direct gyrokinetic and neoclassical simulation
Ideal internal kink stability in presence of plasma flow and neoclassical toroidal viscosity due to energetic particles
Journal Article
·
Wed Dec 31 23:00:00 EST 2008
· Physics of Plasmas
·
OSTI ID:1007844
Tokamak profile prediction using direct gyrokinetic and neoclassical simulation
Journal Article
·
Mon Jun 15 00:00:00 EDT 2009
· Physics of Plasmas
·
OSTI ID:21277207
Ideal internal kink stability in presence of plasma flow and neoclassical toroidal viscosity due to energetic particles
Journal Article
·
Thu Mar 18 00:00:00 EDT 2021
· Nuclear Fusion
·
OSTI ID:1849033