-
First numerical analysis of runaway electron generation in tungsten-rich plasmas towards ITER
Authors:
J. Walkowiak,
M. Hoppe,
I. Ekmark,
A. Jardin,
J. Bielecki,
K. Król,
Y. Savoye-Peysson,
D. Mazon,
D. Dworak,
M. Scholz
Abstract:
The disruption and runaway electron analysis model code was extended to include tungsten impurities in disruption simulations with the aim of studying the runaway electron (RE) generation. This study investigates RE current sensitivity on the following plasma parameters and modelling choices: tungsten concentration, magnetic perturbation strength, electron modelling, thermal quench time and tokama…
▽ More
The disruption and runaway electron analysis model code was extended to include tungsten impurities in disruption simulations with the aim of studying the runaway electron (RE) generation. This study investigates RE current sensitivity on the following plasma parameters and modelling choices: tungsten concentration, magnetic perturbation strength, electron modelling, thermal quench time and tokamak geometry: ITER-like or ASDEX-like. Our investigation shows that a tungsten concentration below 10-3 does not cause significant RE generation on its own. However, at higher concentrations it is possible to reach a very high RE current. Out of the two tested models of electrons in plasma: fluid and isotropic (kinetic), results from the fluid model are more conservative, which is useful when it comes to safety analysis. However, these results are overly pessimistic when compared to the isotropic model, which is based on a more reliable approach. Our results also show that the hot-tail RE generation mechanism is dominant as a primary source of RE in tungsten induced disruptions, usually providing orders of magnitude higher RE seed than Dreicer generation. We discuss best practices for simulations with tungsten-rich plasma, present the dependence of the safety limits on modelling choices and highlight the biggest shortcoming of the current simulation techniques. The obtained results pave the way for a wider analysis of tungsten impact on the disruption dynamics, including the mitigation techniques for ITER in the case of strong contamination of the plasma with tungsten.
△ Less
Submitted 23 April, 2024;
originally announced April 2024.
-
Expulsion of runaway electrons using ECRH in the TCV tokamak
Authors:
J. Decker,
M. Hoppe,
U. Sheikh,
B. P. Duval,
G. Papp,
L. Simons,
T. Wijkamp,
J. Cazabonne,
S. Coda,
E. Devlaminck,
O. Ficker,
R. Hellinga,
U. Kumar,
Y. Savoye-Peysson,
L. Porte,
C. Reux,
C. Sommariva,
A. Tema Biwolé,
B. Vincent,
L. Votta,
the TCV Team,
the EUROfusion Tokamak Exploitation Team
Abstract:
Runaway electrons (REs) are a concern for tokamak fusion reactors from discharge startup to termination. A sudden localized loss of a multi-megaampere RE beam can inflict severe damage to the first wall. Should a disruption occur, the existence of a RE seed may play a significant role in the formation of a RE beam and the magnitude of its current. The application of central electron cyclotron reso…
▽ More
Runaway electrons (REs) are a concern for tokamak fusion reactors from discharge startup to termination. A sudden localized loss of a multi-megaampere RE beam can inflict severe damage to the first wall. Should a disruption occur, the existence of a RE seed may play a significant role in the formation of a RE beam and the magnitude of its current. The application of central electron cyclotron resonance heating (ECRH) in the Tokamak à Configuration Variable (TCV) reduces an existing RE seed population by up to three orders of magnitude within only a few hundred milliseconds. Applying ECRH can also prevent the formation of a post-disruption RE beam where it would otherwise be expected. The RE expulsion rate and consequent RE current reduction are found to increase with applied ECRH power. Whereas central ECRH is effective in expelling REs, off-axis ECRH has a comparatively limited effect. A simple 0-D model for the evolution of the RE population is presented that explains the effective ECRH-induced RE expulsion results from the combined effects of increased electron temperature and enhanced RE transport.
△ Less
Submitted 15 April, 2024;
originally announced April 2024.
-
Fluid and kinetic studies of tokamak disruptions using Bayesian optimization
Authors:
Ida Ekmark,
Mathias Hoppe,
Tünde Fülöp,
Patrik Jansson,
Liam Antonsson,
Oskar Vallhagen,
Istvan Pusztai
Abstract:
When simulating runaway electron dynamics in tokamak disruptions, fluid models with lower numerical cost are often preferred to more accurate kinetic models. The aim of this work is to compare fluid and kinetic simulations of a large variety of different disruption scenarios in ITER. We consider both non-activated and activated scenarios; for the latter we derive and implement kinetic sources for…
▽ More
When simulating runaway electron dynamics in tokamak disruptions, fluid models with lower numerical cost are often preferred to more accurate kinetic models. The aim of this work is to compare fluid and kinetic simulations of a large variety of different disruption scenarios in ITER. We consider both non-activated and activated scenarios; for the latter we derive and implement kinetic sources for the Compton scattering and tritium beta decay runaway electron generation mechanisms in our simulation tool DREAM [M. Hoppe et al 2021 Comp. Phys. Commun. 268, 108098]. To achieve a diverse set of disruption scenarios, Bayesian optimization is used to explore a range of massive material injection densities for deuterium and neon. The cost function is designed to distinguish between successful and unsuccessful disruption mitigation based on the runaway current, current quench time and transported fraction of the heat loss. In the non-activated scenarios, we find that fluid and kinetic disruption simulations can have significantly different runaway electron dynamics, due to an overestimation of the runaway seed by the fluid model. The primary cause of this is that the fluid hot-tail generation model neglects superthermal electron transport losses during the thermal quench. In the activated scenarios, the fluid and kinetic models give similar predictions, which can be explained by the significant influence of the activated sources on the RE dynamics and the seed.
△ Less
Submitted 5 July, 2024; v1 submitted 8 February, 2024;
originally announced February 2024.
-
Runaway Electron Dynamics in ITER Disruptions with Shattered Pellet Injections
Authors:
Oskar Vallhagen,
Lise Hanebring,
Javier Artola,
Michael Lehnen,
Eric Nardon,
Tünde Fülöp,
Mathias Hoppe,
Sarah Newton,
Istvan Pusztai
Abstract:
This study systematically explores the parameter space of disruption mitigation through shattered pellet injection in ITER with a focus on runaway electron dynamics, using the disruption modelling tool DREAM. The physics fidelity is considerably increased compared to previous studies, by e.g., using realistic magnetic geometry, resistive wall configuration, thermal quench onset criteria, as well a…
▽ More
This study systematically explores the parameter space of disruption mitigation through shattered pellet injection in ITER with a focus on runaway electron dynamics, using the disruption modelling tool DREAM. The physics fidelity is considerably increased compared to previous studies, by e.g., using realistic magnetic geometry, resistive wall configuration, thermal quench onset criteria, as well as including additional effects, such as ion transport and enhanced runaway electron transport during the thermal quench. The work aims to provide a fairly comprehensive coverage of experimentally feasible scenarios, considering plasmas representative of both non-activated and high-performance DT operation, different thermal quench onset criteria and transport levels, a wide range of hydrogen and neon quantities injected in one or two stages, and pellets with various characteristic shard sizes. Using a staggered injection scheme, with a pure hydrogen injection preceding a mixed hydrogen-neon injection, we find injection parameters leading to acceptable runaway electron currents in all investigated discharges without activated runaway sources. Dividing the injection into two stages is found to significantly enhance the assimilation and minimize runaway electron generation due to the hot-tail mechanism. However, while a staggered injection outperforms a single stage injection also in cases with radioactive runaway electron sources, no cases with acceptable runaway electron currents are found for a DT-plasma with a 15 MA plasma current.
△ Less
Submitted 25 January, 2024;
originally announced January 2024.
-
Inelastic collisions facilitating runaway electron generation in weakly-ionized plasmas
Authors:
Y. Lee,
P. Aleynikov,
P. C. de Vries,
H. -T. Kim,
J. Lee,
M. Hoppe,
J. -K. Park,
G. J. Choi,
J. Gwak,
Y. -S. Na
Abstract:
Dreicer generation is one of the main mechanisms of runaway electrons generation, in particular during tokamak startup. In fully ionized plasma it is described as a diffusive flow from the Maxwellian core into high energies under the effect of the electric field. In this work we demonstrate a critical role of the non-differential nature of inelastic collisions in weakly ionized plasma during tokam…
▽ More
Dreicer generation is one of the main mechanisms of runaway electrons generation, in particular during tokamak startup. In fully ionized plasma it is described as a diffusive flow from the Maxwellian core into high energies under the effect of the electric field. In this work we demonstrate a critical role of the non-differential nature of inelastic collisions in weakly ionized plasma during tokamak startup, where some electrons experience virtually no collisions during acceleration to the critical energy. We show that using the Fokker-Planck collisional operator can underestimate the Dreicer generation rate by several orders of magnitude.
△ Less
Submitted 23 April, 2024; v1 submitted 24 October, 2023;
originally announced October 2023.
-
Bayesian optimization of massive material injection for disruption mitigation in tokamaks
Authors:
Istvan Pusztai,
Ida Ekmark,
Hannes Bergström,
Peter Halldestam,
Patrik Jansson,
Mathias Hoppe,
Oskar Vallhagen,
Tünde Fülöp
Abstract:
A Bayesian optimization framework is used to investigate scenarios for disruptions mitigated with combined deuterium and neon injection in ITER. The optimization cost function takes into account limits on the maximum runaway current, the transported fraction of the heat loss and the current quench time. The aim is to explore the dependence of the cost function on injected densities, and provide in…
▽ More
A Bayesian optimization framework is used to investigate scenarios for disruptions mitigated with combined deuterium and neon injection in ITER. The optimization cost function takes into account limits on the maximum runaway current, the transported fraction of the heat loss and the current quench time. The aim is to explore the dependence of the cost function on injected densities, and provide insights into the behaviour of the disruption dynamics for representative scenarios. The simulations are conducted using the numerical framework DREAM (Disruption Runaway Electron Analysis Model). We show that irrespective of the quantities of the material deposition, multi-megaampere runaway currents will be produced in the deuterium-tritium phase of operations, even in the optimal scenarios. However, the severity of the outcome can be influenced by tailoring the radial profile of the injected material; in particular if the injected neon is deposited at the edge region it leads to a significant reduction of both the final runaway current and the transported heat losses. The Bayesian approach allows us to map the parameter space efficiently, with more accuracy in favorable parameter regions, thereby providing us information about the robustness of the optima.
△ Less
Submitted 24 February, 2023; v1 submitted 2 February, 2023;
originally announced February 2023.
-
Feasibility of the EDICAM camera for runaway electron detection in JT-60SA disruptions
Authors:
Soma Olasz,
Mathias Hoppe,
Tamás Szepesi,
Kensaku Kamiya,
Peter Balazs,
Gergo I. Pokol
Abstract:
The visible camera system EDICAM (Event Detection Intelligent Camera), recently installed on JT-60SA, is simulated to assess whether it can be used for measuring synchrotron radiation from relativistic runaway electrons. In this simulation, the SOFT synthetic synchrotron diagnostic framework is used to compute the synthetic synchrotron camera images from a JT-60SA-like disruption modelled with the…
▽ More
The visible camera system EDICAM (Event Detection Intelligent Camera), recently installed on JT-60SA, is simulated to assess whether it can be used for measuring synchrotron radiation from relativistic runaway electrons. In this simulation, the SOFT synthetic synchrotron diagnostic framework is used to compute the synthetic synchrotron camera images from a JT-60SA-like disruption modelled with the DREAM disruption simulation code. In the studied scenario, a large amount of argon is added to the plasma, and a disruption is simulated by starting a prescribed exponential temperature drop and finishing with further cooling provided by the argon in a self-consistent simulation of the current quench. The background plasma evolution is calculated by DREAM self-consistently with the fast electron population, which is modelled kinetically. The resulting runaway electron distribution function along with the parameters of the EDICAM visible camera system are used as an input to the SOFT framework to assess the feasibility of the camera for runaway electron detection. We find that the runaway electron beam formed in the disruption can produce synchrotron radiation observable by the EDICAM system, thus enabling the use of the EDICAM for the characterization of runaway electron beams.
△ Less
Submitted 27 September, 2023; v1 submitted 10 January, 2023;
originally announced January 2023.
-
On the minimum transport required to passively suppress runaway electrons in SPARC disruptions
Authors:
R. A. Tinguely,
I. Pusztai,
V. A. Izzo,
K. S'"arkimäki,
T. Fülöp,
D. T. Garnier,
R. S. Granetz,
M. Hoppe,
C. Paz-Soldan,
A. Sundström,
R. Sweeney
Abstract:
In [V.A. Izzo et al 2022 Nucl. Fusion 62 096029], state-of-the-art modeling of thermal and current quench (CQ) MHD coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2…
▽ More
In [V.A. Izzo et al 2022 Nucl. Fusion 62 096029], state-of-the-art modeling of thermal and current quench (CQ) MHD coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1 MA compared to ~8.7 MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.3 and with kinetic energies 0.2-15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 $\mathrm{m^2/s}$, is needed in the core to fully mitigate these REs. Properly accounting for (i) the CQ electric field's effect on RE transport in islands and (ii) the contribution of significant RE currents to disruption MHD may help achieve this.
△ Less
Submitted 3 January, 2023;
originally announced January 2023.
-
Runaway dynamics in reactor-scale spherical tokamak disruptions
Authors:
E. Berger,
I. Pusztai,
S. L. Newton,
M. Hoppe,
O. Vallhagen,
A. Fil,
T. Fülöp
Abstract:
Understanding generation and mitigation of runaway electrons in disruptions is important for the safe operation of future tokamaks. In this paper we investigate runaway dynamics in reactor-scale spherical tokamaks. We study both the severity of runaway generation during unmitigated disruptions, as well as the effect that typical mitigation schemes based on massive material injection have on runawa…
▽ More
Understanding generation and mitigation of runaway electrons in disruptions is important for the safe operation of future tokamaks. In this paper we investigate runaway dynamics in reactor-scale spherical tokamaks. We study both the severity of runaway generation during unmitigated disruptions, as well as the effect that typical mitigation schemes based on massive material injection have on runaway production. The study is conducted using the numerical framework DREAM (Disruption Runaway Electron Analysis Model). We find that, in many cases, mitigation strategies are necessary to prevent the runaway current from reaching multi-megaampere levels. Our results indicate that with a suitably chosen deuterium-neon mixture for mitigation, it is possible to achieve a tolerable runaway current and ohmic current evolution. With such parameters, however, the majority of the thermal energy loss happens through radial transport rather than radiation, which poses a risk of unacceptable localised heat loads.
△ Less
Submitted 3 August, 2022;
originally announced August 2022.
-
Bayesian approach for validation of runaway electron simulations
Authors:
Aaro Järvinen,
Tünde Fülöp,
Eero Hirvijoki,
Mathias Hoppe,
Adam Kit,
Jan Åström
Abstract:
Plasma-terminating disruptions in future fusion reactors may result in conversion of the initial current to a relativistic runaway electron beam. Validated predictive tools are required to optimize the scenarios and mitigation actuators to avoid the excessive damage that can be caused by such events. Many of the simulation tools applied in fusion energy research require the user to specify several…
▽ More
Plasma-terminating disruptions in future fusion reactors may result in conversion of the initial current to a relativistic runaway electron beam. Validated predictive tools are required to optimize the scenarios and mitigation actuators to avoid the excessive damage that can be caused by such events. Many of the simulation tools applied in fusion energy research require the user to specify several input parameters that are not constrained by the available experimental information. Hence, a typical validation exercise requires multiparameter optimization to calibrate the uncertain input parameters for the best possible representation of the investigated physical system. The conventional approach, where an expert modeler conducts the parameter calibration based on domain knowledge, is prone to lead to an intractable validation challenge. For a typical simulation, conducting exhaustive multiparameter investigations manually to ensure a globally optimal solution and to rigorously quantify the uncertainties is an unattainable task, typically covered only partially and unsystematically. Bayesian inference algorithms offer a promising alternative approach that naturally includes uncertainty quantification and is less subjective to user bias in choosing the input parameters. The main challenge in using these methods is the computational cost of simulating enough samples to construct the posterior distributions for the uncertain input parameters. This challenge can be overcome by combining probabilistic surrogate modelling, such as Gaussian Process regression, with Bayesian optimization, which can reduce the number of required simulations by several orders of magnitude. Here, we implement this type of Bayesian optimization framework for a model for analysis of disruption runaway electrons, and explore for simulations of current quench in a JET plasma discharge with an argon induced disruption.
△ Less
Submitted 3 August, 2022;
originally announced August 2022.
-
Runaway dynamics in disruptions with current relaxation
Authors:
István Pusztai,
Mathias Hoppe,
Oskar Vallhagen
Abstract:
The safe operation of tokamak reactors requires a reliable modeling capability of disruptions, and in particular the spatio-temporal dynamics of associated runaway electron currents. In a disruption, instabilities can break up magnetic surfaces into chaotic field line regions, causing current profile relaxation, as well as a rapid radial transport of heat and particles. Using a mean-field helicity…
▽ More
The safe operation of tokamak reactors requires a reliable modeling capability of disruptions, and in particular the spatio-temporal dynamics of associated runaway electron currents. In a disruption, instabilities can break up magnetic surfaces into chaotic field line regions, causing current profile relaxation, as well as a rapid radial transport of heat and particles. Using a mean-field helicity transport model implemented in the disruption runaway modeling framework DREAM, we calculate the dynamics of runaway electrons in the presence of current relaxation events. In scenarios where flux surfaces remain intact in parts of the plasma, a skin current is induced at the boundary of the intact magnetic field region. This skin current region becomes an important center concerning the subsequent dynamics: It may turn into a hot ohmic current channel, or a sizable radially localized runaway beam, depending on the heat transport. If the intact region is in the plasma edge, runaway generation in the counter-current direction can occur, which may develop into a sizable reverse runaway beam. Even when the current relaxation extends to the entire plasma, the final runaway current density profile can be significantly affected, as the induced electric field is reduced in the core and increased in the edge, thereby shifting the center of runaway generation towards the edge.
△ Less
Submitted 2 June, 2022;
originally announced June 2022.
-
Runaway electron generation during tokamak start-up
Authors:
M. Hoppe,
I. Ekmark,
E. Berger,
T. Fülöp
Abstract:
Tokamak start-up is characterized by low electron densities and strong electric fields, in order to quickly raise the plasma current and temperature, allowing the plasma to fully ionize and magnetic flux surfaces to form. Such conditions are ideal for the formation of superthermal electrons, which may reduce the efficiency of ohmic heating and prevent the formation of a healthy thermal fusion plas…
▽ More
Tokamak start-up is characterized by low electron densities and strong electric fields, in order to quickly raise the plasma current and temperature, allowing the plasma to fully ionize and magnetic flux surfaces to form. Such conditions are ideal for the formation of superthermal electrons, which may reduce the efficiency of ohmic heating and prevent the formation of a healthy thermal fusion plasma. This is of particular concern in ITER where engineering limitations put restrictions on the allowable electric fields and limit the prefill densities during start-up. In this study, we present a new 0D burn-through simulation tool called STREAM (STart-up Runaway Electron Analysis Model), which self-consistently evolves the plasma density, temperature and electric field, while accounting for the generation and loss of relativistic runaway electrons. After verifying the burn-through model, we investigate conditions under which runaway electrons can form during tokamak start-up as well as their effects on the plasma initiation. We find that Dreicer generation plays a crucial role in determining whether a discharge becomes runaway-dominated or not, and that a large number of runaway electrons could limit the ohmic heating of the plasma, thus preventing successful burn-through or further ramp-up of the plasma current. The runaway generation can be suppressed by raising the density via gas fuelling, but only if done sufficiently early. Otherwise a large runaway seed may have already been built up, which can avalanche even at relatively low electric fields and high densities.
△ Less
Submitted 29 June, 2022; v1 submitted 18 March, 2022;
originally announced March 2022.
-
Effect of Two-Stage Shattered Pellet Injection on Tokamak Disruptions
Authors:
O. Vallhagen,
I. Pusztai,
M. Hoppe,
S. L. Newton,
T. Fülöp
Abstract:
An effective disruption mitigation system in a tokamak reactor should limit the exposure of the wall to localized heat losses and to the impact of high current runaway electron beams, and avoid excessive forces on the structure. We evaluate with respect to these aspects a two-stage deuterium-neon shattered pellet injection in an ITER-like plasma, using simulations with the DREAM framework [M. Hopp…
▽ More
An effective disruption mitigation system in a tokamak reactor should limit the exposure of the wall to localized heat losses and to the impact of high current runaway electron beams, and avoid excessive forces on the structure. We evaluate with respect to these aspects a two-stage deuterium-neon shattered pellet injection in an ITER-like plasma, using simulations with the DREAM framework [M. Hoppe et al (2021) Comp. Phys. Commun. 268, 108098]. To minimize the obtained runaway currents an optimal range of injected deuterium quantities is found. This range is sensitive to the opacity of the plasma to Lyman radiation, which affects the ionization degree of deuterium, and thus avalanche runaway generation. The two-stage injection scheme, where dilution cooling is produced by deuterium before a radiative thermal quench caused by neon, reduces both the hot-tail seed and the localized transported heat load on the wall. However, during nuclear operation, additional runaway seed sources from the activated wall and tritium make it difficult to reach tolerably low runaway currents.
△ Less
Submitted 26 January, 2022; v1 submitted 25 January, 2022;
originally announced January 2022.
-
Modeling the complete prevention of disruption-generated runaway electron beam formation with a passive 3D coil in SPARC
Authors:
RA Tinguely,
VA Izzo,
DT Garnier,
A Sundström,
K Särkimäki,
O Embréus,
T Fülöp,
RS Granetz,
M Hoppe,
I Pusztai,
R Sweeney
Abstract:
The potential formation of multi-mega-ampere beams of relativistic "runaway" electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetr…
▽ More
The potential formation of multi-mega-ampere beams of relativistic "runaway" electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetric in-vessel coil will passively prevent RE beam formation during disruptions in the SPARC tokamak, a compact, high-field, high-current device capable of achieving a fusion gain Q > 2 in deuterium-tritium plasmas.
△ Less
Submitted 20 October, 2021;
originally announced October 2021.
-
The hot-tail runaway seed landscape during the thermal quench in tokamaks
Authors:
Ida Svenningsson,
Ola Embreus,
Mathias Hoppe,
Sarah L Newton,
Tünde Fülöp
Abstract:
Runaway electron populations seeded from the hot-tail generated by the rapid cooling in plasma-terminating disruptions are a serious concern for next-step tokamak devices such as ITER. Here, we present a comprehensive treatment of the thermal quench, including the superthermal electron dynamics, heat and particle transport, atomic physics, and radial losses due to magnetic perturbations: processes…
▽ More
Runaway electron populations seeded from the hot-tail generated by the rapid cooling in plasma-terminating disruptions are a serious concern for next-step tokamak devices such as ITER. Here, we present a comprehensive treatment of the thermal quench, including the superthermal electron dynamics, heat and particle transport, atomic physics, and radial losses due to magnetic perturbations: processes that are strongly linked and essential for the evaluation of the runaway seed in disruptions mitigated by material injection. We identify limits on the injected impurity density and magnetic perturbation level for which the runaway seed current is acceptable without excessive thermal energy being lost to the wall via particle impact. The consistent modelling of generation and losses shows that runaway beams tend to form near the edge of the plasma, where they could be deconfined via external perturbations.
△ Less
Submitted 3 June, 2021; v1 submitted 7 April, 2021;
originally announced April 2021.
-
DREAM: a fluid-kinetic framework for tokamak disruption runaway electron simulations
Authors:
M. Hoppe,
O. Embreus,
T. Fülöp
Abstract:
Avoidance of the harmful effects of runaway electrons (REs) in plasma-terminating disruptions is pivotal in the design of safety systems for magnetic fusion devices. Here, we describe a computationally efficient numerical tool, that allows for self-consistent simulations of plasma cooling and associated RE dynamics during disruptions. It solves flux-surface averaged transport equations for the pla…
▽ More
Avoidance of the harmful effects of runaway electrons (REs) in plasma-terminating disruptions is pivotal in the design of safety systems for magnetic fusion devices. Here, we describe a computationally efficient numerical tool, that allows for self-consistent simulations of plasma cooling and associated RE dynamics during disruptions. It solves flux-surface averaged transport equations for the plasma density, temperature and poloidal flux, using a bounce-averaged kinetic equation to self-consistently provide the electron current, heat, density and RE evolution, as well as the electron distribution function. As an example, we consider disruption scenarios with material injection and compare the electron dynamics resolved with different levels of complexity, from fully kinetic to fluid modes.
△ Less
Submitted 10 August, 2021; v1 submitted 30 March, 2021;
originally announced March 2021.
-
Validity of models for Dreicer generation of runaway electrons in dynamic scenarios
Authors:
S. Olasz,
O. Embreus,
M. Hoppe,
M. Aradi,
D. Por,
T. Jonsson,
D. Yadikin,
G. I. Pokol,
EU-IM Team
Abstract:
Runaway electron modelling efforts are motivated by the risk these energetic particles pose to large fusion devices. The sophisticated kinetic models can capture most features of the runaway electron generation but have high computational costs which can be avoided by using computationally cheaper reduced kinetic codes. In this paper, we compare the reduced kinetic and kinetic models to determine…
▽ More
Runaway electron modelling efforts are motivated by the risk these energetic particles pose to large fusion devices. The sophisticated kinetic models can capture most features of the runaway electron generation but have high computational costs which can be avoided by using computationally cheaper reduced kinetic codes. In this paper, we compare the reduced kinetic and kinetic models to determine when the former solvers, based on analytical calculations assuming quasi-stationarity, can be used. The Dreicer generation rate is calculated by two different solvers in parallel in a workflow developed in the European Integrated Modelling framework, and this is complemented by calculations of a third code that is not yet integrated into the framework. Runaway Fluid, a reduced kinetic code, NORSE, a kinetic code using non-linear collision operator, and DREAM, a linearized Fokker-Planck solver, are used to investigate the effect of a dynamic change in the electric field for different plasma scenarios spanning across the whole tokamak-relevant range. We find that on time scales shorter than or comparable to the electron collision time at the critical velocity for runaway electron generation kinetic effects not captured by reduced kinetic models play an important role. This characteristic time scale is easy to calculate and can reliably be used to determine whether there is a need for kinetic modelling, or cheaper reduced kinetic codes are expected to deliver sufficiently accurate results. This criterion can be automated, and thus it can be of great benefit for the comprehensive self-consistent modelling frameworks that are attempting to simulate complex events such as tokamak start-up or disruptions.
△ Less
Submitted 27 May, 2021; v1 submitted 7 January, 2021;
originally announced January 2021.
-
Tomographic reconstruction of the runaway distribution function in TCV using multispectral synchrotron images
Authors:
T. A. Wijkamp,
A. Perek,
J. Decker,
B. Duval,
M. Hoppe,
G. Papp,
U. A. Sheikh,
I. G. J. Classen,
R. J. E. Jaspers,
the TCV team,
the EUROfusion MST1 team
Abstract:
Synchrotron radiation observed in a quiescent TCV runaway discharge is studied using filtered camera images targeting three distinct wavelength intervals. Through the tomographic SART procedure the high momentum, high pitch angle part of the spatial and momentum distribution of these relativistic particles is reconstructed. Experimental estimates of the distribution are important for verification…
▽ More
Synchrotron radiation observed in a quiescent TCV runaway discharge is studied using filtered camera images targeting three distinct wavelength intervals. Through the tomographic SART procedure the high momentum, high pitch angle part of the spatial and momentum distribution of these relativistic particles is reconstructed. Experimental estimates of the distribution are important for verification and refinement of formation-, decay- and transport-models underlying runaway avoidance and mitigation strategy design. Using a test distribution it is demonstrated that the inversion procedure provides estimates accurate to within a few tens of percent in the region of phase-space contributing most to the synchrotron image. We find that combining images filtered around different parts of the emission spectrum widens the probed part of momentum-space and reduces reconstruction errors. Next, the SART algorithm is used to obtain information on the spatiotemporal runaway momentum distribution in a selected TCV discharge. The momentum distribution is found to relax towards an avalanche-like exponentially decaying profile. Anomalously high pitch angles and a radial profile increasing towards the edge are found for the most strongly emitting particles in the distribution. Pitch angle scattering by toroidal magnetic field ripple is consistent with this picture. An alternative explanation is the presence of high frequency instabilities in combination with the formation of a runaway shell at the edge of the plasma.
△ Less
Submitted 5 February, 2021; v1 submitted 30 October, 2020;
originally announced November 2020.
-
Runaway electron modelling in the self-consistent core European Transport Simulator, ETS
Authors:
Gergo I. Pokol,
Soma Olasz,
Boglarka Erdos,
Gergely Papp,
Matyas Aradi,
Mathias Hoppe,
Thomas Johnson,
Jorge Ferreira,
David Coster,
Yves Peysson,
Joan Decker,
Par Strand,
Dimitriy Yadikin,
Denis Kalupin,
the EUROfusion-IM Team
Abstract:
Relativistic runaway electrons are a major concern in tokamaks. The European framework for Integrated Modelling (EU-IM), facilitates the integration of different plasma simulation tools by providing a standard data structure for communication that enables relatively easy integration of different physics codes. A three-level modelling approach was adopted for runaway electron simulations within the…
▽ More
Relativistic runaway electrons are a major concern in tokamaks. The European framework for Integrated Modelling (EU-IM), facilitates the integration of different plasma simulation tools by providing a standard data structure for communication that enables relatively easy integration of different physics codes. A three-level modelling approach was adopted for runaway electron simulations within the EU-IM. Recently, a number of runaway electron modelling modules have been integrated into this framework. The first level of modelling (Runaway Indicator) is limited to the indication if runaway electron generation is possible or likely. The second level (Runaway Fluid) adopts an approach similar to e.g. the GO code, using analytical formulas to estimate changes in the runaway electron current density. The third level is based on the solution of the electron kinetics. One such code is LUKE that can handle the toroidicity-induced effects by solving the bounce-averaged Fokker-Planck equation. Another approach is used in NORSE, which features a fully nonlinear collision operator that makes it capable of simulating major changes in the electron distribution, for example slide-away. Both codes handle the effect of radiation on the runaway distribution. These runaway-electron modelling codes are in different stages of integration into the EU-IM infrastructure, and into the European Transport Simulator (ETS), which is a fully capable modular 1.5D core transport simulator. ETS with Runaway Fluid was benchmarked to the GO code implementing similar physics. Coherent integration of kinetic solvers requires more effort on the coupling, especially regarding the definition of the boundary between runaway and thermal populations, and on consistent calculation of resistivity. Some of these issues are discussed.
△ Less
Submitted 30 September, 2020;
originally announced September 2020.
-
Spatiotemporal analysis of the runaway distribution function from synchrotron images in an ASDEX Upgrade disruption
Authors:
M. Hoppe,
L. Hesslow,
O. Embreus,
L. Unnerfelt,
G. Papp,
I. Pusztai,
T. Fülöp,
O. Lexell,
T. Lunt,
E. Macusova,
P. J. McCarthy,
G. Pautasso,
G. I. Pokol,
G. Por,
P. Svensson,
the ASDEX Upgrade team,
the EUROfusion MST1 team
Abstract:
Synchrotron radiation images from runaway electrons (REs) in an ASDEX Upgrade discharge disrupted by argon injection are analyzed using the synchrotron diagnostic tool SOFT and coupled fluid-kinetic simulations. We show that the evolution of the runaway distribution is well described by an initial hot-tail seed population, which is accelerated to energies between 25-50 MeV during the current quenc…
▽ More
Synchrotron radiation images from runaway electrons (REs) in an ASDEX Upgrade discharge disrupted by argon injection are analyzed using the synchrotron diagnostic tool SOFT and coupled fluid-kinetic simulations. We show that the evolution of the runaway distribution is well described by an initial hot-tail seed population, which is accelerated to energies between 25-50 MeV during the current quench, together with an avalanche runaway tail which has an exponentially decreasing energy spectrum. We find that, although the avalanche component carries the vast majority of the current, it is the high-energy seed remnant that dominates synchrotron emission. With insights from the fluid-kinetic simulations, an analytic model for the evolution of the runaway seed component is developed and used to reconstruct the radial density profile of the RE beam. The analysis shows that the observed change of the synchrotron pattern from circular to crescent shape is caused by a rapid redistribution of the radial profile of the runaway density.
△ Less
Submitted 5 February, 2021; v1 submitted 29 May, 2020;
originally announced May 2020.
-
Runaway electron synchrotron radiation in a vertically translated plasma
Authors:
M. Hoppe,
G. Papp,
T. Wikjamp,
A. Perek,
J. Decker,
B. Duval,
O. Embreus,
T. Fülöp,
U. A. Sheikh
Abstract:
Synchrotron radiation observed from runaway electrons (REs) in tokamaks depends upon the position and size of the RE beam, the RE energy and pitch distributions, as well as the location of the observer. We show that experimental synchrotron images of a vertically moving runaway electron beam sweeping past the detector in the TCV tokamak agree well with predictions from the synthetic synchrotron di…
▽ More
Synchrotron radiation observed from runaway electrons (REs) in tokamaks depends upon the position and size of the RE beam, the RE energy and pitch distributions, as well as the location of the observer. We show that experimental synchrotron images of a vertically moving runaway electron beam sweeping past the detector in the TCV tokamak agree well with predictions from the synthetic synchrotron diagnostic Soft. This experimental validation lends confidence to the theory underlying the synthetic diagnostics which are used for benchmarking theoretical models of and probing runaway dynamics. We present a comparison of synchrotron measurements in TCV with predictions of kinetic theory for runaway dynamics in uniform magnetic fields. We find that to explain the detected synchrotron emission, significant non-collisional pitch angle scattering as well as radial transport of REs would be needed. Such effects could be caused by the presence of magnetic perturbations, which should be further investigated in future TCV experiments.
△ Less
Submitted 18 August, 2020; v1 submitted 23 March, 2020;
originally announced March 2020.
-
Evaluation of the Dreicer runaway generation rate in the presence of high-Z impurities using a neural network
Authors:
L Hesslow,
L Unnerfelt,
O Vallhagen,
O Embreus,
M Hoppe,
G Papp,
T Fülöp
Abstract:
Integrated modelling of electron runaway requires computationally expensive kinetic models that are self-consistently coupled to the evolution of the background plasma parameters. The computational expense can be reduced by using parameterized runaway generation rates rather than solving the full kinetic problem. However, currently available generation rates neglect several important effects; in p…
▽ More
Integrated modelling of electron runaway requires computationally expensive kinetic models that are self-consistently coupled to the evolution of the background plasma parameters. The computational expense can be reduced by using parameterized runaway generation rates rather than solving the full kinetic problem. However, currently available generation rates neglect several important effects; in particular, they are not valid in the presence of partially ionized impurities. In this work, we construct a multilayer neural network for the Dreicer runaway generation rate which is trained on data obtained from kinetic simulations performed for a wide range of plasma parameters and impurities. The neural network accurately reproduces the Dreicer runaway generation rate obtained by the kinetic solver. By implementing it in a fluid runaway electron modelling tool, we show that the improved generation rates lead to significant differences in the self-consistent runaway dynamics as compared to the results using the previously available formulas for the runaway generation rate.
△ Less
Submitted 24 January, 2020; v1 submitted 1 October, 2019;
originally announced October 2019.
-
Experimental and synthetic measurements of polarized synchrotron emission from runaway electrons in Alcator C-Mod
Authors:
R. A. Tinguely,
M. Hoppe,
R. S. Granetz,
R. T. Mumgaard,
S. Scott
Abstract:
This paper presents the first experimental analysis of polarized synchrotron emission from relativistic runaway electrons (REs) in a tokamak plasma. Importantly, we show that the polarization information of synchrotron radiation can be used to diagnose spatially-localized RE pitch angle distributions. Synchrotron-producing REs were generated during low density, Ohmic, diverted plasma discharges in…
▽ More
This paper presents the first experimental analysis of polarized synchrotron emission from relativistic runaway electrons (REs) in a tokamak plasma. Importantly, we show that the polarization information of synchrotron radiation can be used to diagnose spatially-localized RE pitch angle distributions. Synchrotron-producing REs were generated during low density, Ohmic, diverted plasma discharges in the Alcator C-Mod tokamak. The ten-channel Motional Stark Effect diagnostic was used to measure spatial profiles of the polarization angle $θ_{\mathrm{pol}}$ and the fraction $\mathrm{f}_{\mathrm{pol}}$ of detected light that was linearly-polarized. Spatial transitions in $θ_{\mathrm{pol}}$ of 90$°$---from horizontal to vertical polarization and vice versa---are observed in experimental data and are well-explained by the gyro-motion of REs and high directionality of synchrotron radiation. Polarized synchrotron emission is modeled with the synthetic diagnostic SOFT; its output Green's (or detector response) functions reveal a critical RE pitch angle at which $θ_{\mathrm{pol}}$ flips by 90$°$ and $\mathrm{f}_{\mathrm{pol}}$ is minimal. Using SOFT, we determine the dominant RE pitch angle which reproduces measured $θ_{\mathrm{pol}}$ and $\mathrm{f}_{\mathrm{pol}}$ values. The spatiotemporal evolutions of $θ_{\mathrm{pol}}$ and $\mathrm{f}_{\mathrm{pol}}$ are explored in detail for one C-Mod discharge. For channels viewing REs near the magnetic axis and flux surfaces $q$ = 1 and 4/3, disagreements between synthetic and experimental signals suggest that the sawtooth instability may be influencing RE dynamics. Furthermore, other sources of pitch angle scattering, not considered in this analysis, could help explain discrepancies between simulation and experiment.
△ Less
Submitted 26 June, 2019;
originally announced June 2019.
-
Spatiotemporal evolution of runaway electrons from synchrotron images in Alcator C-Mod
Authors:
R. A. Tinguely,
R. S. Granetz,
M. Hoppe,
O. Embréus
Abstract:
In the Alcator C-Mod tokamak, relativistic runaway electron (RE) generation can occur during the flattop current phase of low density, diverted plasma discharges. Due to the high toroidal magnetic field (B = 5.4 T), RE synchrotron radiation is measured by a wide-view camera in the visible wavelength range (~400-900 nm). In this paper, a statistical analysis of over one thousand camera images is pe…
▽ More
In the Alcator C-Mod tokamak, relativistic runaway electron (RE) generation can occur during the flattop current phase of low density, diverted plasma discharges. Due to the high toroidal magnetic field (B = 5.4 T), RE synchrotron radiation is measured by a wide-view camera in the visible wavelength range (~400-900 nm). In this paper, a statistical analysis of over one thousand camera images is performed to investigate the plasma conditions under which synchrotron emission is observed in C-Mod. In addition, the spatiotemporal evolution of REs during one particular discharge is explored in detail via a thorough analysis of the distortion-corrected synchrotron images. To accurately predict RE energies, the kinetic solver CODE [Landreman et al 2014 Comput. Phys. Commun. 185 847-855] is used to evolve the electron momentum-space distribution at six locations throughout the plasma: the magnetic axis and flux surfaces q = 1, 4/3, 3/2, 2, and 3. These results, along with the experimentally-measured magnetic topology and camera geometry, are input into the synthetic diagnostic SOFT [Hoppe et al 2018 Nucl. Fusion 58 026032] to simulate synchrotron emission and detection. Interesting spatial structure near the surface q = 2 is found to coincide with the onset of a locked mode and increased MHD activity. Furthermore, the RE density profile evolution is fit by comparing experimental to synthetic images, providing important insight into RE spatiotemporal dynamics.
△ Less
Submitted 5 October, 2018;
originally announced October 2018.
-
Generalized collision operator for fast electrons interacting with partially ionized impurities
Authors:
L Hesslow,
O Embréus,
M Hoppe,
T C DuBois,
G Papp,
M Rahm,
T Fülöp
Abstract:
Accurate modelling of the interaction between fast electrons and partially ionized atoms is important for evaluating tokamak disruption mitigation schemes based on material injection. This requires accounting for the effect of screening of the impurity nuclei by the cloud of bound electrons. In this paper, we detail the derivation of a generalized collision operator including the effect of partial…
▽ More
Accurate modelling of the interaction between fast electrons and partially ionized atoms is important for evaluating tokamak disruption mitigation schemes based on material injection. This requires accounting for the effect of screening of the impurity nuclei by the cloud of bound electrons. In this paper, we detail the derivation of a generalized collision operator including the effect of partial screening. We calculate the effective ion length-scales, which are needed in the components of the collision operator, for a number of ion species commonly appearing in fusion experiments. We show that for high electric fields, the secondary runaway growth rate can be substantially larger than in a fully ionized plasma with the same effective charge, although the growth rate is significantly reduced at near-critical electric fields. Furthermore, by comparison with the Boltzmann collision operator, we show that the Fokker--Planck formalism is valid even for large impurity content.
△ Less
Submitted 27 November, 2018; v1 submitted 13 July, 2018;
originally announced July 2018.
-
Dynamics of positrons during relativistic electron runaway
Authors:
Ola Embréus,
Linnea Hesslow,
Mathias Hoppe,
Gergely Papp,
Katya Richards,
Tünde Fülöp
Abstract:
Sufficiently strong electric fields in plasmas can accelerate charged particles to relativistic energies. In this paper we describe the dynamics of positrons accelerated in such electric fields, and calculate the fraction of created positrons that become runaway accelerated, along with the amount of radiation that they emit. We derive an analytical formula that shows the relative importance of the…
▽ More
Sufficiently strong electric fields in plasmas can accelerate charged particles to relativistic energies. In this paper we describe the dynamics of positrons accelerated in such electric fields, and calculate the fraction of created positrons that become runaway accelerated, along with the amount of radiation that they emit. We derive an analytical formula that shows the relative importance of the different positron production processes, and show that above a certain threshold electric field the pair production by photons is lower than that by collisions. We furthermore present analytical and numerical solutions to the positron kinetic equation; these are applied to calculate the fraction of positrons that become accelerated or thermalized, which enters into rate equations that describe the evolution of the density of the slow and fast positron populations. Finally, to indicate operational parameters required for positron detection during runaway in tokamak discharges, we give expressions for the parameter dependencies of detected annihilation radiation compared to bremsstrahlung detected at an angle perpendicular to the direction of runaway acceleration. Using the full leading order pair production cross section, we demonstrate that previous related work has overestimated the collisional pair production by at least a factor of four.
△ Less
Submitted 12 July, 2018;
originally announced July 2018.
-
Measurements of runaway electron synchrotron spectra at high magnetic fields in Alcator C-Mod
Authors:
R. A. Tinguely,
R. S. Granetz,
M. Hoppe,
O. Embreus
Abstract:
In the Alcator C-Mod tokamak, runaway electron (RE) experiments have been performed during low density, flattop plasma discharges at three magnetic fields: 2.7, 5.4, and 7.8 T, the last being the highest field to-date at which REs have been generated and measured in a tokamak. Time-evolving synchrotron radiation spectra were measured in the visible wavelength range (~300-1000 nm) by two absolutely…
▽ More
In the Alcator C-Mod tokamak, runaway electron (RE) experiments have been performed during low density, flattop plasma discharges at three magnetic fields: 2.7, 5.4, and 7.8 T, the last being the highest field to-date at which REs have been generated and measured in a tokamak. Time-evolving synchrotron radiation spectra were measured in the visible wavelength range (~300-1000 nm) by two absolutely-calibrated spectrometers viewing co- and counter-plasma current directions. In this paper, a test particle model is implemented to predict momentum-space and density evolutions of REs on the magnetic axis and q = 1, 3/2, and 2 surfaces. Drift orbits and subsequent loss of confinement are also incorporated into the evolution. These spatiotemporal results are input into the new synthetic diagnostic SOFT [M. Hoppe, et al., Nucl. Fusion 58(2), 026032 (2018)] which reproduces experimentally-measured spectra. For these discharges, it is inferred that synchrotron radiation dominates collisional friction as a power loss mechanism and that RE energies decrease as magnetic field is increased. Additionally, the threshold electric field for RE generation, as determined by hard X-ray and photo-neutron measurements, is compared to current theoretical predictions.
△ Less
Submitted 14 May, 2018;
originally announced May 2018.
-
Interpretation of runaway electron synchrotron and bremsstrahlung images
Authors:
M. Hoppe,
O. Embréus,
C. Paz-Soldan,
R. A. Moyer,
T. Fülöp
Abstract:
The crescent spot shape observed in DIII-D runaway electron synchrotron radiation images is shown to result from the high degree of anisotropy in the emitted radiation, the finite spectral range of the camera and the distribution of runaways. The finite spectral camera range is found to be particularly important, as the radiation from the high-field side can be stronger by a factor $10^6$ than the…
▽ More
The crescent spot shape observed in DIII-D runaway electron synchrotron radiation images is shown to result from the high degree of anisotropy in the emitted radiation, the finite spectral range of the camera and the distribution of runaways. The finite spectral camera range is found to be particularly important, as the radiation from the high-field side can be stronger by a factor $10^6$ than the radiation from the low-field side in DIII-D. By combining a kinetic model of the runaway dynamics with a synthetic synchrotron diagnostic we see that physical processes not described by the kinetic model (such as radial transport) are likely to be limiting the energy of the runaways. We show that a population of runaways with lower dominant energies and larger pitch-angles than those predicted by the kinetic model provide a better match to the synchrotron measurements. Using a new synthetic bremsstrahlung diagnostic we also simulate the view of the Gamma Ray Imager (GRI) diagnostic used at DIII-D to resolve the spatial distribution of runaway-generated bremsstrahlung.
△ Less
Submitted 12 December, 2017;
originally announced December 2017.
-
SOFT: A synthetic synchrotron diagnostic for runaway electrons
Authors:
Mathias Hoppe,
Ola Embréus,
R. Alexander Tinguely,
Robert S. Granetz,
Adam Stahl,
Tünde Fülöp
Abstract:
Improved understanding of the dynamics of runaway electrons can be obtained by measurement and interpretation of their synchrotron radiation emission. Models for synchrotron radiation emitted by relativistic electrons are well established, but the question of how various geometric effects -- such as magnetic field inhomogeneity and camera placement -- influence the synchrotron measurements and the…
▽ More
Improved understanding of the dynamics of runaway electrons can be obtained by measurement and interpretation of their synchrotron radiation emission. Models for synchrotron radiation emitted by relativistic electrons are well established, but the question of how various geometric effects -- such as magnetic field inhomogeneity and camera placement -- influence the synchrotron measurements and their interpretation remains open. In this paper we address this issue by simulating synchrotron images and spectra using the new synthetic synchrotron diagnostic tool SOFT (Synchrotron-detecting Orbit Following Toolkit). We identify the key parameters influencing the synchrotron radiation spot and present scans in those parameters. Using a runaway electron distribution function obtained by Fokker-Planck simulations for parameters from an Alcator C-Mod discharge, we demonstrate that the corresponding synchrotron image is well-reproduced by SOFT simulations, and we explain how it can be understood in terms of the parameter scans. Geometric effects are shown to significantly influence the synchrotron spectrum, and we show that inherent inconsistencies in a simple emission model (i.e. not modeling detection) can lead to incorrect interpretation of the images.
△ Less
Submitted 16 November, 2017; v1 submitted 3 September, 2017;
originally announced September 2017.
-
Benford's Law and Distractors in Multiple Choice Exams
Authors:
Fred M. Hoppe
Abstract:
Suppose that in a multiple choice examination the leading digit of the correct options follows Benford's Law, while the the leading digit of the distractors are uniform. Consider a strategy for guessing at answers that selects the option with the lowest leading digit with ties broken at random. We provide an expression for the probability that this strategy selects the correct option.
Suppose that in a multiple choice examination the leading digit of the correct options follows Benford's Law, while the the leading digit of the distractors are uniform. Consider a strategy for guessing at answers that selects the option with the lowest leading digit with ties broken at random. We provide an expression for the probability that this strategy selects the correct option.
△ Less
Submitted 6 May, 2014; v1 submitted 26 November, 2013;
originally announced November 2013.