Generation Of Runaway Electrons Research Articles

Fluid and kinetic studies of tokamak disruptions using Bayesian optimization

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 (Hoppe et al., Comput. Phys. Commun., vol. 268, 2021, 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 runaway dynamics and the seed.

Experimental research of ECW pre-ionization and assisted startup in EAST

Experimental research on the electron cyclotron wave (ECW) pre-ionization and assisted start-up was carried out systematically for the first time in EAST tokamak, which is a superconducting device with ITER-like full metal wall. Breakdown and plasma initiation at low toroidal electric fields (<0.3 V m−1) with ECW pre-ionization and startup assistance has been demonstrated. Also, the parameter domain of breakdown is significantly extended towards higher prefill gas pressure. The effect of ECW injection timing, power, toroidal injection angle on breakdown were also investigated. Injecting ECW earlier leads to an earlier breakdown and a higher plasma current ramp rate. The electron cyclotron heating (ECH) power threshold for breakdown in EAST is approximately 0.4 MW. In the range of ECH power tested in this work, higher ECH power is advantageous for achieving earlier and faster breakdown. Furthermore, the breakdown with radial ECW injection occurs earlier compared with oblique injections (co-current and counter-current). During the ECW-assisted startup, the process of burn-through is prolonged by the higher pre-filled gas pressure even though it enhances the ease of breakdown. In addition, compared to the low hybrid wave assistance, the ECW assistance has an effect in averting the generation of runaway electrons and improving the safety of device during startup. Moreover, the ECW assistance exhibits a high tolerance to the impurity and thus ensures a high ramp rate of plasma current even with a high impurity level.

Optimal Conditions for the Generation of Runaway Electrons in High-Pressure Gases

Runaway electron (RAE) generation in high-pressure gases is an important physical phenomenon that significantly influences discharge shapes and properties of initiated plasma. The diffuse discharges formed due to RAEs in the air and other gases at atmospheric pressure find wide applications. In the present review, theoretical and experimental results that explain the reason for RAE occurrence at high pressures are analyzed, and recommendations are given for the implementation of conditions under which the runaway electron beam (RAEB) with the highest current can be obtained at atmospheric pressure. The experimental results were obtained using subnanosecond, nanosecond, and submicrosecond generators, including those specially developed for runaway electron generation. The RAEBs were recorded using oscilloscopes and collectors with picosecond time resolution. To theoretically describe the phenomenon of continuous electron acceleration, the method of physical kinetics was used based on the Boltzmann kinetic equation that takes into account the minimum but sufficient number of elementary processes, including shock gas ionization and elastic electron scattering. The results of modeling allowed the main factors to be established that control the RAE appearance, the most important of which is electron scattering on neutral atoms and/or molecules. Theoretical modeling has allowed the influence of various parameters (including the voltage, pressure, gas type, and geometrical characteristics of the discharge gap) to be taken into account. The results of the research presented here allow RAE accelerators with desirable parameters to be developed and the possibility of obtaining diffuse discharges to be accessed under various conditions. The review consists of the Introduction, five sections, the Conclusion, and the References.

Generation of Runaway Electrons near Micro-Inhomogeneities on the Cathode Surface in Subnanosecond Self-Sustained Discharges in a Wide Range of High Pressures

Generation of Runaway Electrons near Micro-Inhomogeneities on the Cathode Surface in Subnanosecond Self-Sustained Discharges in a Wide Range of High Pressures

Generation of Runaway Electrons near Micro-Inhomogeneities on the Cathode Surface in Subnanosecond Self-Sustained Discharges in a Wide Range of High Pressures

The results of numerical 3D modeling of the development of an electron avalanche initiated by a field emission electron in a small-sized region of an amplified electric field near the microinhomogeneities at the cathode have been presented. The simulation has been carried out for the discharge gaps with an initially homogeneous distribution of the electric field with a reduced intensity significantly lower than that required by the electron runaway criterion. The possibility of the transition of the field emission electrons initiating avalanches and the electrons in these avalanches into runaway regime has been investigated. The microinhomogeneities in the form of a cone, metal droplets, and boundaries between pores or microcraters have been considered. The calculations were carried out for nitrogen in the pressure range from atmospheric to 40 atm. It has been shown that the initial energy obtained near the microinhomogeneity can significantly facilitate the transition of the electron into the runaway mode. And the electron will continue to run away in a discharge gap electric field weak according to the runaway criterion. It has been shown that this effect is especially noticeable at gas pressures above 10 atm. A comparative analysis of the simulation results with the experimental data obtained by us on the switching characteristics of a discharge gap filled with nitrogen when exposed to voltage pulses with subnanosecond fronts of different steepness has been carried out. This made it possible to divide the ranges of experimental conditions into those when only the amplification of the electric field near the microinhomogeneities is sufficient for the runaway of electrons and when the electric field of an avalanche of critical or close to critical size is additionally necessary for the runaway.

Formation of the hot tail seeds for runaway electron generation during disruptions with lower hybrid waves in the HL-2A tokamak

The hot tail generation is expected to be the dominant mechanism for the runaway electron (RE) seed formation during disruptions, especially in large devices with high electron temperature such as international thermonuclear experimental reactor. This issue has been studied in the HL-2A tokamak by using the superthermal electrons produced by lower hybrid waves (LHWs), which can adjust the hot tail distribution. It was observed that RE generation was significantly enhanced during disruptions with LHWs. The measurements show that the multitudinous superthermal electrons with energy of 40–60 keV created by LHWs greatly transform the landscape of hot tail distribution. The tail electrons can be directly converted into REs under the acceleration of the high toroidal electric field during disruptions. Runaway current plateaus are more likely to be formed than in normal disruptions without LHWs. However, some abnormal phenomena have also been observed, that is, RE generation was not enhanced and no runaway current plateau was formed during some disruptions with LHWs. It is found that this is attributed to the complete loss of RE seeds caused by strong magnetic fluctuations, which prevents the generation of REs during disruptions. This may provide a way to avoid the generation of REs during disruptions by actively exciting magnetic fluctuations.

A numerical investigation on electron runaway threshold at the initial stage of atmospheric streamer development

Pre-ionization caused by runaway electrons is an important mechanism for negative streamer development. The aim of this paper is to investigate the runaway criteria and overvoltage threshold of electrons at the initial stage of streamer development in air, with the self-developed 3D particle-in-cell with Monte Carlo Collision code. First, numerical simulations are performed with fixed number of electrons to study the runaway criteria in nonrelativistic cases. This method takes the stochastic fluctuations of collisions into account and solves the major shortcomings of theoretical approach. The simulated critical electric field is less than that of the theoretical approach, and the amplitude of the difference increases with electron energy, due to the “tunneling effect” caused by the stochastic fluctuations of collisions. Then, simulations of negative streamers at various applied voltages are performed to investigate the overvoltage threshold. A more intuitive method, searching energetic electrons in front of the negative streamer head, which corresponds to the nature of runaway electrons, is applied to determine the generation of runaway electrons. Electrons that escaped a certain distance ahead of the streamer can be observed at 30 kV. Thus, the overvoltage threshold for runaway electrons can be roughly estimated as 3.3 in our simulations, which is about three times less than the previously published one. At last, with the redefined overvoltage threshold, the figure of regions of breakdown development for various mechanisms depending upon the overvoltage in air is updated.

Kinetics of the fast ionization waves with runaway electrons

In nanosecond gas breakdown, runaway electrons are responsible for the interelectrode space pre-ionization, which supports the propagation of fast ionization waves or streamers. However, a detailed understanding of mechanisms of how this happens is still elusive because of the wide range of energies of runaway electrons, from the hundreds of eV to the hundreds of keV. While the low-energy part of the runaway electron population can only ionize the background gas, the high-energy part can produce x-rays that can also ionize the gas. The objective of this paper is to understand the influence of different electron groups on the fast ionization wave dynamics for the conditions typical for runaway electrons generation. We find that there are two main groups, one having energies below 100 eV and another having energies above 1 keV. The first group defines the establishment of the fast ionization waves and their plasma parameters, while the second group defines the speed of the ionization front propagation.

Cross-machine comparison of runaway electron generation during tokamak start-up for extrapolation to ITER

A cross-machine comparison of global parameters that determine the runaway electron (RE) generation and loss process during tokamak start-up was carried out with the aim to extrapolate these to ITER. The study found that all considered discharges, also those that do not show signs of RE, are non-thermal at the start, i.e. have a streaming parameter larger than 0.1. During the current ramp-up the electric field, E, remains above the critical value, E c, that allows RE in the plasma. The distinction to be made is not if RE can form but, if sufficient RE can form fast enough such that they are detected or start to dominate the dynamics of the tokamak discharge. The dynamics of the value of E, density and temperature during tokamak are key to the formation of RE. It was found that larger devices operate with E closer to E c, due to their higher temperatures, hence the RE generation is relatively slower. The slower time scales for the formation of RE, estimated to be of the order of 100s of ms in ITER simplifies the development of avoidance schemes. The RE confinement time is also an important determinant of the entire process and is found to increase with the device size. The study also revealed that drift orbit losses, a mechanism often attributed as the main RE loss mechanism during the early tokamak discharge, are actually more difficult to achieve. RE losses might be more likely attributed to RE diffusion due to magnetic turbulence.

The effects of pulse voltage rise time on the nanosecond pulsed breakdown of nitrogen spark switch at atmospheric pressure with 3D PIC-MCC model

In this paper, the effects of pulse voltage rise time on the nanosecond pulsed breakdown of the nitrogen spark switch at atmospheric pressure are analyzed. Based on the assumption of initial electrons generation due to the field emission, the simulations are performed using a three-dimensional particle-in-cell, Monte Carlo-collision model for the pulse voltage with a rise time of 40, 60, and 100 ns, respectively. The breakdown experiments of the nitrogen spark switch are carried out for three different rise times. The results obtained are as follows. First, the nanosecond pulsed breakdown of the switch includes the formation and fast propagation of the streamer, which depend on the multiplication of the electron avalanche, and the intense ionization due to photoelectrons and energetic electrons, respectively. Second, with the rise time of pulse voltage increasing, the generation of runaway electrons becomes more difficult and the streamer branches, which are mainly caused by photoionization and captured energetic electrons, become more obvious. Finally, the breakdown time delay of the switch becomes shorter and the breakdown voltage becomes higher at the same pressure for the decreasing rise time of pulse voltage, which is consistent with the measurement results.

Conditions for Inception of Relativistic Runaway Discharges in Air

AbstractTerrestrial gamma‐ray flashes are linked to growth of long bidirectional lightning leader system consisting of positive and stepping negative leaders. The spatial extent of streamer zones of a typical lightning leader with tip potential exceeding several tens of megavolts is on the order of 10–100 m. The photoelectric absorption of bremsstrahlung radiation generated by avalanching relativistic runaway electrons occurs efficiently on the same spatial scales. The intense multiplication of these electrons is triggered when the size of the negative leader streamer zone crosses a threshold of approximately 100 m (for sea‐level air pressure conditions) allowing self‐replication of these avalanches due to the upstream relativistic electron seeds generated by the photoelectric absorption. The model results also highlight importance of electrode effects in interpretation of X‐ray emissions from centimeter to meter long laboratory discharges, in particular, a similar feedback effect produced by generation of runaway electrons from the cathode material.

Development of a Shattered Pellet Injector test bench for the ITER DMS support laboratory

In the ITER Shattered Pellet Injection system large, Ø28.5x57mm cryogenic pellets have to be created and preserved in a barrel, until the control system detects a forthcoming disruption event and triggers the propellant gas valve to launch the pellet. The high-pressure propellant gas accelerates the pellet, which shatters on an angled plate near the edge of the plasma. The ablation of the resulting pellet fragments is intended to mitigate the consequences of the disruption and to prevent the generation of runaway electrons. To support the Shattered Pellet Injector technology for the ITER Disruption Mitigation System, a new support laboratory has been built at the Centre for Energy Research. The aim of the laboratory is to develop a test bench with ITER relevant dimensions, produce d=28.5mm Hydrogen, Neon and mixture pellets, observe the pellet acceleration and shattering with automated control. This contribution describes the cryogenic pellet injector of the test bench.

On the minimum transport required to passively suppress runaway electrons in SPARC disruptions

In Izzo et al (2022 Nucl. Fusion 62 096029), state-of-the-art modeling of thermal and current quench (CQ) magnetohydrodynamics (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 compared to of pre-disruption plasma current. Further investigation finds that core-localized electrons, within and with kinetic energies –, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient , is needed in the core to fully mitigate these REs. Properly accounting for (a) the CQ electric field’s effect on RE transport in islands and (b) the contribution of significant RE currents to disruption MHD may help achieve this.

Mechanism of runaway electron generation in nanosecond pulsed plate-plate discharge at atmospheric-pressure air

Classical discharge theory (Townsend theory and streamer theory) has limitations in explaining nanosecond pulsed gas discharge. In recent years, the research on nanosecond pulsed gas discharge theory based on the high-energy runaway electrons has attracted extensive attention. But so far, there have been few studies of the generation mechanism of runaway electrons in atmospheric-pressure-air nanosecond pulsed plate-to-plate discharge, which seriously hinders the application and development of nanosecond pulse discharge plasma. In this paper, a one-dimensional implicit particle-in-cell/Monte Carlo collision (PIC/MCC) model is developed to investigate the mechanism of runaway electron generation and breakdown in a 1 mm-long atmospheric-pressure-air gap between the plate electrode and plate electrode driven by a negative nanosecond pulse voltage with an amplitude of 20 kV. The results show that under the influence of space charge dynamic behavior, the electric field enhancement region appears between the plate electrode and plate electrode, so that electrons can satisfy the electron runaway criteria and behaves in the runaway mode. In addition, it is also observed that the pre-ionization effect of the runaway electrons in front of the discharge channel can cause the secondary electron avalanches. As the secondary electrons avalanche and the discharge channel continues to converge, the discharge is guided and accelerated, eventually leading to the breakdown of the air gap. This study further reveals the mechanism of nanosecond pulsed plate-plate discharge, expands the basic theory of nanosecond pulsed gas discharge, and opens up new opportunities for the application and development of nanosecond pulsed discharge plasma.

Positive and negative streamers in air and nitrogen in a sharply inhomogeneous electric field under conditions of runaway electron generation

AbstractThe development of positive and negative streamers in a point‐to‐plane gap filled with air and nitrogen at various pressures (50–200 kPa) and voltages (8–25 kV) was studied. A four‐channel intensified charge‐coupled device (ICCD) and a streak‐camera were used. Electrical parameters were measured with high resolution (10 GHz). An original method of measuring a displacement current caused by a streamer was applied. As was expected, positive streamer branches in nitrogen at low voltages and/or elevated pressures while a large‐diameter streamer is formed in air. However, at high voltages or negative polarity, the large‐diameter streamer is formed both in nitrogen and air in the entire pressure range. It was found that runaway electrons (REs) are generated in the very first picoseconds of gas ionisation near the pointed cathode. It was assumed that REs can be generated near the pointed anode and produce bremsstrahlung radiation due to a sharply inhomogeneous distribution of electrical potential. It was found that in the final stage of negative streamer development in air and nitrogen, the gas between the streamer front and the opposite electrode is ionised almost simultaneously in the entire volume when the streamer diameter is almost equal to the interelectrode gap.

Runaway electron dynamics in Experimental Advanced Superconducting Tokamak helium plasmas

The generation of runaway electrons (REs) is observed during the low-density helium ohmic plasma discharge in the Experimental Advanced Superconducting Tokamak (EAST). The growth rate of hard x-ray (HXR) is inversely proportional to the line-average density. Besides, the RE generation in helium plasma is higher than that in deuterium plasma at the same density, which is obtained by comparing the growth rate of HXR with the same discharge conditions. The potential reason is the higher electron temperature of helium plasma in the same current and electron density plateau. Furthermore, two Alfvén eigenmodes driven by REs have been observed. The frequency evolution of the mode is not fully satisfied with the Alfvén scaling and when extension of the Alfvén frequency is towards 0, the high frequency branch is ∼50 kHz. The different spatial position of the two modes and the evolution of the helium concentration could be used to understand deviation between theoretical and experimental observation.

Runaway dynamics in reactor-scale spherical tokamak disruptions

Understanding generation and mitigation of runaway electrons in disruptions is important for the safe operation of future tokamaks. In this paper we investigate the runaway dynamics in reactor-scale spherical tokamaks, focusing on a compact nominal design with a plasma current of 21 megaamperes (MA), 1.8 T magnetic field on axis and major radius of approximately 3 m. We study both the severity of runaway generation during unmitigated disruptions, and 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-MA 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. However, this does not account for the runaway source due to wall activation, which has been found to severely limit successful mitigation at conventional aspect ratios, but whose definition requires a more complete wall specification. Furthermore, the majority of the thermal energy loss is found to happen through radial transport rather than radiation, which poses a risk of unacceptable localised heat loads.

Confinement of passing and trapped runaway electrons in the simulation of an ITER current quench

Runaway electrons (REs) present a high-priority R&D issue for ITER but little is known about the extent to which RE generation is affected by the stochastic field intrinsic to disrupting plasmas. RE generation can be modelled with reduced kinetic models and there has been recent progress in involving losses due to field stochasticity, either via a loss-time parameter or radial transport coefficients which can be estimated by tracing test electrons in 3D fields. We evaluate these terms in ITER using a recent JOREK 3D MHD simulation of plasma disruption to provide the stochastic magnetic fields where RE markers are traced with the built-in particle tracing module. While the MHD simulation modelled only the current quench phase, the case is MHD unstable and exhibits similar relaxation as would be expected during the thermal quench. Therefore, the RE simulations can be considered beginning right after the thermal quench but before the MHD relaxation is complete. The plasma is found to become fully stochastic for 8 ms and the resulting transport is sufficient to overcome RE avalanche before flux surfaces are reformed. We also study transport mechanisms for trapped REs and find those to be deconfined as well during this phase. While the results presented here are not sufficient to assess the magnitude of the formed RE beam, we show that significant RE losses could be expected to arise due to field stochasticity.

Numerical Study on the Effect of Plasma Density on Runaway Electron Suppression in the ADITYA-U Tokamak

Runaway Electrons (REs) generated during plasma disruptions in fusion grade tokamaks have the potential to severely damage the plasma-facing components. Designing optimal plasma discharge scenarios for RE suppression in future experiments requires interpretative modeling of current experiments. Multiple experiments have been carried out on ADITYA-U tokamak to design optimal plasma discharge scenarios for RE avoidance and suppression. In this article, we have numerically studied two representative pairs of ADITYA-U plasma discharges. In the first pair of discharges, plasma density was increased by gas puffs during the flap-top phase which is shown to suppress RE generation. In the second pair of the representative discharges, the effect of a lower ratio of peak electric field during breakdown to the pre-fill pressure is shown on RE generation. We simulate these plasma discharges using the PREDICT code to study the dynamics of Res. The results are consistent with the experimental Hard X-Ray diagnostic observations for both the pair of representative discharges. Additionally, the lowest RE-current is obtained in the discharges with a low ratio of the peak electric field during breakdown to the pre-fill pressure and high plasma density during the flap-top phase due to gas puffs. The suppression of REs is demonstrated by showing the effect of increasing plasma density on the separatrix in momentum space between thermal electrons and REs.

Runaway electron generation during tokamak start-up

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.