Search

Runaway electron loads onto material structures are a major concern for future large tokamaks due to the efficient avalanching at high plasma currents. Here, we perform predictive studies using the JOREK code for a plausible plasma configuration in the European DEMO fusion power plant with focus on a pessimistic scenario in which a multi mega-ampere runaway electron beam is formed. The work first comprises axisymmetric predictions of runaway electron beam formation in a mitigated scenario and of the simultaneous vertical motion of the beam due to loss of position control. The subsequent runaway electron beam termination triggered by a burst of MHD activity during the course of the vertical motion is then simulated in 3D with the runaway electron fluid self-consistently coupled to the MHD modes. Finally, the resulting deposition pattern of the runaway electrons onto wall structures is calculated with a relativistic test particle approach. This way, the suitability of a possible sacrificial limiter concept for the protection of first wall components is assessed.

Mechanisms of plasmoid drift following massive material injection are studied via 3D non-linear MHD modelling with the JOREK code, using a transient neutral source deposited at the low field side midplane of a JET H-mode plasma to clarify basic processes and compare with existing theories. The simulations confirm the important role of the propagation of shear Alfv\'en wave (SAW) packets from both ends of the plasmoid (``SAW braking'') and the development of external resistive currents along magnetic field lines (``P\'egouri\'e braking'') in limiting charge separation and thus the $\mathbf{E}\times \mathbf{B}$ plasmoid drift, where $\mathbf{E}$ and $\mathbf{B}$ are the electric and magnetic fields, respectively. The drift velocity is found to be limited by the SAW braking on the few microseconds timescale for cases with relatively small source amplitude while the P\'egouri\'e braking acting on a longer timescale is shown to set in earlier with larger toroidal extent of the source, both in good agreement with existing theories. The simulations also identify the key role of the size of the $\mathbf{E}\times \mathbf{B}$ flow region on plasmoid drift and show that the saturated velocity caused by dominant SAW braking agrees well with theory when considering an effective pressure within the $\mathbf{E}\times \mathbf{B}$ flow region. The existence of SAWs in the simulations is demonstrated and the 3D picture of plasmoid drift is discussed.

The heat flux mitigation during the Thermal Quench (TQ) by the Shattered Pellet Injection (SPI) is one of the major elements of disruption mitigation strategy for ITER. It's efficiency greatly depends on the SPI and the target plasma, and is ultimately checked by the heat deposition on to the PFCs. To investigate this, JOREK simulations of neon-mixed dual-SPIs into ITER baseline H-mode and a "degraded H-mode" with and without good injector synchronization are performed with focus on the first wall heat flux and its energy impact. It is found that low neon fraction SPIs into the baseline H-mode plasmas exhibit strong major radial plasmoid drift as the fragments arrive at the pedestal, accompanied by edge stochasticity. Significant density expulsion and outgoing heat flux occurs as a result, reducing the mitigation efficiency. Such drift motion could be mitigated by injecting higher neon fraction pellets', or by considering the pre-disruption confinement degradation, thus improving the radiation fraction. The radiation heat flux is found to peak in the vicinity of the fragment injection location in the early injection phase, while it relaxes later on due to parallel impurity transport. The overall radiation asymmetry could be significantly mitigated by good synchronization. Time integration of the local heat flux is carried out to provide its energy impact for wall heat damage assessment. For the baseline H-mode case with full pellet injection, melting of the stainless steel of the diagnostic port could occur near the injection port, which is acceptable, without any melting of the first wall tungsten tiles. For the degraded H-mode cases with quarter-pellet SPIs, which have 1/4 total volume of a full pellet, the maximum energy impact approaches the tolerable limit of the stainless steel with un-synchronized SPIs, and stays well below such limit for the perfectly synchronized ones.