A new hybridLagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma
Abstract
In order to enable kinetic simulation of nonthermal edge plasmas at a reduced computational cost, a new hybridLagrangian δf scheme has been developed that utilizes the phase space grid in addition to the usual marker particles, taking advantage of the computational strengths from both sides. The new scheme splits the particle distribution function of a kinetic equation into two parts. Marker particles contain the fast spacetime varying, δf, part of the distribution function and the coarsegrained phasespace grid contains the slow spacetime varying part. The coarsegrained phasespace grid reduces the memoryrequirement and the computing cost, while the marker particles provide scalable computing ability for the finegrained physics. Weights of the marker particles are determined by a direct weight evolution equation instead of the differential form weight evolution equations that the conventional deltaf schemes use. The particle weight can be slowly transferred to the phase space grid, thereby reducing the growth of the particle weights. The nonLagrangian part of the kinetic equation – e.g., collision operation, ionization, charge exchange, heatsource, radiative cooling, and others – can be operated directly on the phase space grid. Deviation of the particle distribution function on the velocity grid from a Maxwellian distribution function – drivenmore »
 Authors:

 Princeton Plasma Physics Laboratory, Princeton University, Princeton, NJ 08543 (United States)
 National Fusion Research Institute (Korea, Republic of)
 University of Colorado Boulder (United States)
 Publication Date:
 OSTI Identifier:
 22572320
 Resource Type:
 Journal Article
 Journal Name:
 Journal of Computational Physics
 Additional Journal Information:
 Journal Volume: 315; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA); Journal ID: ISSN 00219991
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; CHARGE EXCHANGE; COLLISIONS; DISTRIBUTION FUNCTIONS; HEAT SOURCES; HEAT TRANSFER; IONIZATION; KINETIC EQUATIONS; LAGRANGIAN FUNCTION; LOSSES; PHASE SPACE; PLASMA; RADIATIVE COOLING; SPACETIME; TOKAMAK DEVICES; VELOCITY
Citation Formats
Ku, S., Email: sku@pppl.gov, Hager, R., Chang, C. S., Kwon, J. M., and Parker, S. E. A new hybridLagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma. United States: N. p., 2016.
Web. doi:10.1016/J.JCP.2016.03.062.
Ku, S., Email: sku@pppl.gov, Hager, R., Chang, C. S., Kwon, J. M., & Parker, S. E. A new hybridLagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma. United States. doi:10.1016/J.JCP.2016.03.062.
Ku, S., Email: sku@pppl.gov, Hager, R., Chang, C. S., Kwon, J. M., and Parker, S. E. Wed .
"A new hybridLagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma". United States. doi:10.1016/J.JCP.2016.03.062.
@article{osti_22572320,
title = {A new hybridLagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma},
author = {Ku, S., Email: sku@pppl.gov and Hager, R. and Chang, C. S. and Kwon, J. M. and Parker, S. E.},
abstractNote = {In order to enable kinetic simulation of nonthermal edge plasmas at a reduced computational cost, a new hybridLagrangian δf scheme has been developed that utilizes the phase space grid in addition to the usual marker particles, taking advantage of the computational strengths from both sides. The new scheme splits the particle distribution function of a kinetic equation into two parts. Marker particles contain the fast spacetime varying, δf, part of the distribution function and the coarsegrained phasespace grid contains the slow spacetime varying part. The coarsegrained phasespace grid reduces the memoryrequirement and the computing cost, while the marker particles provide scalable computing ability for the finegrained physics. Weights of the marker particles are determined by a direct weight evolution equation instead of the differential form weight evolution equations that the conventional deltaf schemes use. The particle weight can be slowly transferred to the phase space grid, thereby reducing the growth of the particle weights. The nonLagrangian part of the kinetic equation – e.g., collision operation, ionization, charge exchange, heatsource, radiative cooling, and others – can be operated directly on the phase space grid. Deviation of the particle distribution function on the velocity grid from a Maxwellian distribution function – driven by ionization, charge exchange and wall loss – is allowed to be arbitrarily large. The numerical scheme is implemented in the gyrokinetic particle code XGC1, which specializes in simulating the tokamak edge plasma that crosses the magnetic separatrix and is in contact with the material wall.},
doi = {10.1016/J.JCP.2016.03.062},
journal = {Journal of Computational Physics},
issn = {00219991},
number = ,
volume = 315,
place = {United States},
year = {2016},
month = {6}
}
Works referencing / citing this record:
Invited Review Article: Gas puff imaging diagnostics of edge plasma turbulence in magnetic fusion devices
journal, April 2017
 Zweben, S. J.; Terry, J. L.; Stotler, D. P.
 Review of Scientific Instruments, Vol. 88, Issue 4
Verification of long wavelength electromagnetic modes with a gyrokineticfluid hybrid model in the XGC code
journal, May 2017
 Hager, Robert; Lang, Jianying; Chang, C. S.
 Physics of Plasmas, Vol. 24, Issue 5
Analysis of equilibrium and turbulent fluxes across the separatrix in a gyrokinetic simulation
journal, July 2018
 Keramidas Charidakos, I.; Myra, J. R.; Parker, S.
 Physics of Plasmas, Vol. 25, Issue 7
Blobhole correlation model for edge turbulence and comparisons with NSTX gas puff imaging data
journal, June 2018
 Myra, J. R.; Zweben, S. J.; Russell, D. A.
 Plasma Physics and Controlled Fusion, Vol. 60, Issue 7