A backward Monte-Carlo method for time-dependent runaway electron simulations

Research output: Contribution to journalArticlepeer-review

19 Scopus citations

Abstract

Kinetic descriptions of runaway electrons (REs) are usually based on Fokker-Planck models that determine the probability distribution function of REs in 2-dimensional momentum space. Despite the simplification involved, the Fokker-Planck equation can rarely be solved analytically and direct numerical approaches [e.g., continuum and particle-based Monte Carlo (MC)] can be time consuming, especially in the computation of asymptotic-type observables including the runaway probability, the slowing-down and runaway mean times, and the energy limit probability. Here, we present a novel backward MC approach to these problems based on backward stochastic differential equations that describe the dynamics of the runaway probability by means of the Feynman-Kac theory. The key ingredient of the backward MC algorithm is to place all the particles in a runaway state and simulate them backward from the terminal time to the initial time. As such, our approach can provide much faster convergence than direct MC methods (by significantly reducing the number of particles required to achieve a prescribed accuracy) while at the same time maintaining the advantages of particle-based methods (compared to continuum approaches). The proposed algorithm is unconditionally stable and can be parallelized as easy as the direct MC method, and its extension to dimensions higher than two is straightforward, thus paving the way for conducting large-scale RE simulation.

Original languageEnglish
Article number092511
JournalPhysics of Plasmas
Volume24
Issue number9
DOIs
StatePublished - Sep 1 2017

Funding

We thank E. Hirvijoki for valuable comments and suggestions. This material is based on work supported in part by the U.S. Department of Energy, Office of Science, Offices of Advanced Scientific Computing Research and Fusion Energy Science; and by the Laboratory Directed Research and Development program at the Oak Ridge National Laboratory, which is operated by UT-Battelle, LLC., for the U.S. Department of Energy under Contract No. DE-AC05-00OR22725.

Fingerprint

Dive into the research topics of 'A backward Monte-Carlo method for time-dependent runaway electron simulations'. Together they form a unique fingerprint.

Cite this