DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

Abstract

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

Authors:
ORCiD logo [1]; ORCiD logo [1]; ORCiD logo [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC), Fusion Energy Sciences (FES)
OSTI Identifier:
1827585
Alternate Identifier(s):
OSTI ID: 1829947
Report Number(s):
LA-UR-21-22936
Journal ID: ISSN 0021-9991; TRN: US2216042
Grant/Contract Number:  
89233218CNA000001
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 449; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Mathematics; Magnetic Fusion Energy

Citation Formats

Chacon, Luis, Daniel, Don, and Taitano, William T. An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks. United States: N. p., 2021. Web. doi:10.1016/j.jcp.2021.110772.
Chacon, Luis, Daniel, Don, & Taitano, William T. An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks. United States. https://doi.org/10.1016/j.jcp.2021.110772
Chacon, Luis, Daniel, Don, and Taitano, William T. Wed . "An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks". United States. https://doi.org/10.1016/j.jcp.2021.110772. https://www.osti.gov/servlets/purl/1827585.
@article{osti_1827585,
title = {An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks},
author = {Chacon, Luis and Daniel, Don and Taitano, William T.},
abstractNote = {We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.},
doi = {10.1016/j.jcp.2021.110772},
journal = {Journal of Computational Physics},
number = ,
volume = 449,
place = {United States},
year = {Wed Oct 13 00:00:00 EDT 2021},
month = {Wed Oct 13 00:00:00 EDT 2021}
}

Works referenced in this record:

Spatial transport of runaway electrons in axisymmetric tokamak plasmas
journal, January 2019

  • McDevitt, Christopher J.; Guo, Zehua; Tang, Xian-Zhu
  • Plasma Physics and Controlled Fusion, Vol. 61, Issue 2
  • DOI: 10.1088/1361-6587/aaf4d1

An asymptotic-preserving semi-Lagrangian algorithm for the time-dependent anisotropic heat transport equation
journal, September 2014

  • Chacón, L.; del-Castillo-Negrete, D.; Hauck, C. D.
  • Journal of Computational Physics, Vol. 272
  • DOI: 10.1016/j.jcp.2014.04.049

Asymptotic Preserving schemes for highly oscillatory Vlasov–Poisson equations
journal, September 2013

  • Crouseilles, Nicolas; Lemou, Mohammed; Méhats, Florian
  • Journal of Computational Physics, Vol. 248
  • DOI: 10.1016/j.jcp.2013.04.022

NORSE: A solver for the relativistic non-linear Fokker–Planck equation for electrons in a homogeneous plasma
journal, March 2017


Differential form of the collision integral for a relativistic plasma
journal, October 1987


Phase-space dynamics of runaway electrons in magnetic fields
journal, February 2017

  • Guo, Zehua; McDevitt, Christopher J.; Tang, Xian-Zhu
  • Plasma Physics and Controlled Fusion, Vol. 59, Issue 4
  • DOI: 10.1088/1361-6587/aa5952

Numerical characterization of bump formation in the runaway electron tail
journal, January 2016


Kinetic modelling of runaway electron avalanches in tokamak plasmas
journal, July 2015


Runaway electron production in DIII-D killer pellet experiments, calculated with the CQL3D/KPRAD model
journal, November 2000

  • Harvey, R. W.; Chan, V. S.; Chiu, S. C.
  • Physics of Plasmas, Vol. 7, Issue 11
  • DOI: 10.1063/1.1312816

Conservative magnetic moment of runaway electrons and collisionless pitch-angle scattering
journal, August 2018


Relativistic limitations on runaway electrons
journal, June 1975


A fully implicit, scalable, conservative nonlinear relativistic Fokker–Planck 0D-2P solver for runaway electrons
journal, September 2020


Nonlinear relativistic gyrokinetic Vlasov-Maxwell equations
journal, December 1999

  • Brizard, Alain J.; Chan, Anthony A.
  • Physics of Plasmas, Vol. 6, Issue 12
  • DOI: 10.1063/1.873742

Runaway electron generation in axisymmetric tokamak geometry
journal, August 2019


Numerical calculation of the runaway electron distribution function and associated synchrotron emission
journal, March 2014

  • Landreman, Matt; Stahl, Adam; Fülöp, Tünde
  • Computer Physics Communications, Vol. 185, Issue 3
  • DOI: 10.1016/j.cpc.2013.12.004

Effect of partially ionized impurities and radiation on the effective critical electric field for runaway generation
journal, June 2018

  • Hesslow, L.; Embréus, O.; Wilkie, G. J.
  • Plasma Physics and Controlled Fusion, Vol. 60, Issue 7
  • DOI: 10.1088/1361-6587/aac33e

A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems
journal, May 1954


Runaway Electrons in a Plasma
journal, September 1973

  • Kulsrud, Russell M.; Sun, Yung-Chiun; Winsor, Niels K.
  • Physical Review Letters, Vol. 31, Issue 11
  • DOI: 10.1103/PhysRevLett.31.690

Electron and Ion Runaway in a Fully Ionized Gas. I
journal, July 1959


Asymptotic-Preserving Scheme for the Resolution of Evolution Equations with Stiff Transport Terms
journal, January 2019

  • Fedele, Baptiste; Negulescu, Claudia; Possanner, Stefan
  • Multiscale Modeling & Simulation, Vol. 17, Issue 1
  • DOI: 10.1137/18M1175732