ASCOT: solving the kinetic equation of minority particle species in tokamak plasmas
arXiv:1308.1904 · doi:10.1016/j.cpc.2014.01.014
Abstract
A comprehensive description of methods, suitable for solving the kinetic equation for fast ions and impurity species in tokamak plasmas using Monte Carlo approach, is presented. The described methods include Hamiltonian orbit-following in particle and guiding center phase space, test particle or guiding center solution of the kinetic equation applying stochastic differential equations in the presence of Coulomb collisions, neoclassical tearing modes and Alfvén eigenmodes as electromagnetic perturbations relevant to fast ions, together with plasma flow and atomic reactions relevant to impurity studies. Applying the methods, a complete reimplementation of the well-established minority species code ASCOT is carried out as a response both to the increase in computing power during the last twenty years and to the weakly structured growth of the code, which has made implementation of additional models impractical. Also, a benchmark between the previous code and the reimplementation is accomplished, showing good agreement between the codes.
13 pages, 9 figures, submitted to Computer Physics Communications
References in corpus (2)
Cited by in corpus (32)
- Nonlinear Energetic Particle Transport in the Presence of Multiple Alfvenic Waves in ITER
- A quasi-isodynamic configuration with good confinement of fast ions at low plasma
- Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil
- Optimizing beam-ion confinement in ITER by adjusting the toroidal phase of the 3-D magnetic fields applied for ELM control
- Degenerate Variational Integrators for Magnetic Field Line Flow and Guiding Center Trajectories
- Mechanics of ELM control coil induced alpha particle transport
- An advection-diffusion model for cross-field runaway electron transport in perturbed magnetic fields
- Validating the ASCOT modelling of NBI fast ions in Wendelstein 7-X stellarator
- Monte Carlo method and High Performance Computing for solving Fokker-Planck equation of minority plasma particles
- Gyrokinetic modelling of anisotropic energetic particle driven instabilities in tokamak plasmas
- On the minimum transport required to passively suppress runaway electrons in SPARC disruptions
- Calculating the 3D magnetic field of ITER for European TBM studies
- Drift kinetic theory of alpha transport by tokamak perturbations
- Validating neutral-beam current drive simulations in the TJ-II stellarator
- Effects of resonant magnetic perturbations on neutral beam heating in a tokamak
- Role of the edge electric field in the resonant mode-particle interactions and the formation of transport barriers in toroidal plasmas
- Validating fast-ion wall-load IR analysis-methods against W7-X NBI empty-torus experiment
- Kinetic versus Magnetic Chaos in Toroidal Plasmas: A systematic quantitative comparison
- Simulation of neutral beam current drive on EAST tokamak
- Testing the conservative character of particle simulations: II. Spurious heating of guiding centers and full orbits subject to fluctuations expressed in terms of and
- Performance of the BGSDC integrator for computing fast ion trajectories in nuclear fusion reactors
- Adjoint Monte Carlo Simulation of Fusion Product Activation Probe Experiment in ASDEX Upgrade tokamak
- Energetic Particle Tracing in Optimized Quasisymmetric Stellarator Equilibria
- Synthetic fast ion diagnostics in tokamaks: comparing the Monte Carlo test particle code ASCOT against experiments
- Validation of neutron emission and neutron energy spectrum calculations on MAST with DRESS
- Moment-Preserving and Mesh-Adaptive Reweighting Method for Rare-Event Sampling in Monte-Carlo Algorithms
- Adjoint Monte Carlo calculation of charged plasma particle flux to wall
- FOXTAIL: Modeling the nonlinear interaction between Alfvén eigenmodes and energetic particles in tokamaks
- Quasilinear drift kinetic theory of alpha particle transport by neoclassical tearing modes
- High-order coupling of shear and sonic continua in JET plasmas
- Eliminating poor statistics in Monte-Carlo simulations of fast-ion losses to plasma-facing components and detectors
- Modeling the complete prevention of disruption-generated runaway electron beam formation with a passive 3D coil in SPARC