81 research outputs found
An advection-diffusion model for cross-field runaway electron transport in perturbed magnetic fields
Disruption-generated runaway electrons (RE) present an outstanding issue for
ITER. The predictive computational studies of RE generation rely on
orbit-averaged computations and, as such, they lack the effects from the
magnetic field stochasticity. Since stochasiticity is naturally present in
post-disruption plasma, and externally induced stochastization offers a
prominent mechanism to mitigate RE avalanche, we present an advection-diffusion
model that can be used to couple an orbit-following code to an orbit-averaged
tool in order to capture the cross-field transport and to overcome the latter's
limitation. The transport coefficients are evaluated via a Monte Carlo method.
We show that the diffusion coefficient differs significantly from the
well-known Rechester-Rosenbluth result. We also demonstrate the importance of
including the advection: it has a two-fold role both in modelling transport
barriers created by magnetic islands and in amplifying losses in regions where
the islands are not present
Adjoint Monte Carlo Simulation of Fusion Product Activation Probe Experiment in ASDEX Upgrade tokamak
The activation probe is a robust tool to measure flux of fusion products from
a magnetically confined plasma. A carefully chosen solid sample is exposed to
the flux, and the impinging ions transmute the material making it radioactive.
Ultra-low level gamma-ray spectroscopy is used post mortem to measure the
activity and, thus, the number of fusion products.
This contribution presents the numerical analysis of the first measurement in
the ASDEX Upgrade tokamak, which was also the first experiment to measure a
single discharge. The ASCOT suite of codes was used to perform adjoint/reverse
Monte Carlo calculations of the fusion products. The analysis facilitates, for
the first time, a comparison of numerical and experimental values for
absolutely calibrated flux. The results agree to within a factor of about two,
which can be considered a quite good result considering the fact that all
features of the plasma cannot be accounted in the simulations.
Also an alternative to the present probe orientation was studied. The results
suggest that a better optimized orientation could measure the flux from a
significantly larger part of the plasma.Comment: Contribution in 1st EPS Conference on Plasma Diagnostics. First two
versions are for PoS(ECPD 2015)055. This 3rd version was accepted for
publishing in Journal of Instrumentatio
ASCOT: solving the kinetic equation of minority particle species in tokamak plasmas
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\'en 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.Comment: 13 pages, 9 figures, submitted to Computer Physics Communication
Calculating the 3D magnetic field of ITER for European TBM studies
The magnetic perturbation due to the ferromagnetic test blanket modules
(TBMs) may deteriorate fast ion confinement in ITER. This effect must be
quantified by numerical studies in 3D. We have implemented a combined finite
element method (FEM) -- Biot-Savart law integrator method (BSLIM) to calculate
the ITER 3D magnetic field and vector potential in detail. Unavoidable geometry
simplifications changed the mass of the TBMs and ferritic inserts (FIs) up to
26%. This has been compensated for by modifying the nonlinear ferromagnetic
material properties accordingly. Despite the simplifications, the computation
geometry and the calculated fields are highly detailed. The combination of
careful FEM mesh design and using BSLIM enables the use of the fields
unsmoothed for particle orbit-following simulations. The magnetic field was
found to agree with earlier calculations and revealed finer details. The vector
potential is intended to serve as input for plasma shielding calculations.Comment: In proceedings of the 28th Symposium on Fusion Technolog
Modelling of 3D fields due to ferritic inserts and test blanket modules in toroidal geometry at ITER
Computations in toroidal geometry are systematically performed for the plasma response to 3D magnetic perturbations produced by ferritic inserts (FIs) and test blanket modules (TBMs) for four ITER plasma scenarios: the 15 MA baseline, the 12.5 MA hybrid, the 9 MA steady state, and the 7.5 MA half-field helium plasma. Due to the broad toroidal spectrum of the FI and TBM fields, the plasma response for all the n = 1-6 field components are computed and compared. The plasma response is found to be weak for the high-n (n > 4) components. The response is not globally sensitive to the toroidal plasma flow speed, as long as the latter is not reduced by an order of magnitude. This is essentially due to the strong screening effect occurring at a finite flow, as predicted for ITER plasmas. The ITER error field correction coils (EFCC) are used to compensate the n = 1 field errors produced by FIs and TBMs for the baseline scenario for the purpose of avoiding mode locking. It is found that the middle row of the EFCC, with a suitable toroidal phase for the coil current, can provide the best correction of these field errors, according to various optimisation criteria. On the other hand, even without correction, it is predicted that these n = 1 field errors will not cause substantial flow damping for the 15 MA baseline scenario
- …