3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract

Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power (Prad) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an underdetermined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling ("guess and check") approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 +0.03/-0.29 and 1.01 +0.02/-0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

Radiation asymmetry in JET disruption mitigation experiments with shattered pellet injection

“…Studies on the pellet injector trajectories, MHD modes, role of plasma rotation, complimented with JOREK [24] simulations and Emis3D analyses [25,26] are ongoing to justify this physics insight and will be presented in a dedicated manuscript.…”

Interpretative 3D MHD modelling of deuterium SPI into a JET H-mode plasma

“…It can be seen that core T e starts to drop from t ≈ 2 ms (line y) along with the penetration and assimilation of the first group of fragments. P rad starts to increase rapidly from t ≈ 6 ms (line z), based on either an inversion using the vertical bolometry and assuming axisymmetric radiation (solid blue trace in (b)) or a 3D tomographic reconstruction without assuming axisymmetric radiation (dashed red) using the Emis3D tool [45,46]. W th gradually drops mainly via conduction up to t ≈ 6 ms (line z), after which radiation plays an increasingly important role, as seen from panel (d).…”

Plasmoid drift and first wall heat deposition during ITER H-mode dual-SPIs in JOREK simulations

“…There have been intensive studies, both experimentally [2][3][4][5][6] and numerically [7][8][9][10], on both the total radiated energy and the radiation power density within the plasma volume after SPIs. It has been found in general that an apparent toroidal asymmetry in the radiation power density exists for the single SPI cases, especially in the early injection phase.…”