Search NASA⌕ Search

Engineering topics

Sheikh, U. (ORCID:0000000162072489)

Publications and source records attributed to Sheikh, U. (ORCID:0000000162072489).

MHD, disruptions and control physics: Chapter 4 of the special issue: on the path to tokamak burning plasma operation

In this chapter, we review the progress in MHD stability, disruptions and control in magnetic fusion research that has occurred over the past (more than) one and a half decades since the publication by Hender et al in 2007 on the same topic as part of the update of ITER Physics Basis. During this period, remarkable progress has been achieved in the understanding of the basic physics and overall control of MHD instabilities through a wide spectrum of dedicated experiments, theory and modeling. The sawtooth activities are probably today one of the best understood of MHD events and very robust control schemes have been developed for reliable operation of tokamaks through core heating. Similarly, significant improvements have been achieved in understanding and control of neoclassical tearing modes, resistive wall modes or locked modes and their control through ECCD or error field control. The field of disruption prediction through application of artificial intelligence, machine learning or deep learning methods, which had already started at the time of the 2007 review, has progressed significantly due to general progress in these fields and application of newer, more sophisticated algorithms. However, although remarkable progress has been achieved in the field of Disruptions, their understanding, prediction, possible avoidance and mitigation still remain probably the most active fields of R&D globally in this field. This is especially because reactor grade machines like ITER and DEMO will be much less tolerant in respect of disruptions and runaway currents, and their occurrences must be either avoided altogether or minimized to an acceptable value without causing any significant hindrance to robust machine operations. This review is intended to present a broad spectrum of the R&D that has occurred in this field in support of ITER, which will also be of immense significance for all future machines, especially reactors like DEMO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An upper pressure limit for low- Z benign termination of runaway electron beams in TCV

Abstract We present a model for the particle balance in the post-disruption runaway electron plateau phase of a tokamak discharge. The model is constructed with the help of, and applied to, experimental data from TCV discharges investigating the so-called ‘low- Z benign termination’ runaway electron mitigation scheme. In the benign termination scheme, the free electron density is first reduced in order for a subsequently induced MHD instability to grow rapidly and spread the runaway electrons widely across the wall. We show that the observed non-monotonic dependence of the free electron density with the measured neutral pressure is due to plasma re-ionization induced by runaway electron impact ionization. At higher neutral pressures, more target particles are present in the plasma for runaway electrons to collide with and ionize. Parameter scans are conducted to clarify the role of the runaway electron density and energy on the free electron density, and it is found that only the runaway electron density has a noticeable impact. While the free electron density is shown to be related to the spread of heat fluxes at termination, the exact cause for the upper neutral pressure limit remains undetermined and an object for further study.

Hoppe, M. (ORCID:0000000339948977)↗

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

Abstract The pre-thermal quench (pre-TQ) dynamics of a pure deuterium ( D 2 ) shattered pellet injection (SPI) into a 3 MA / 7 MJ JET H-mode plasma is studied via 3D non-linear MHD modelling with the JOREK code. The interpretative modelling captures the overall evolution of the measured density and radiated power. The simulations also identify the importance of the drifts of ablation plasmoids towards the tokamak low field side (LFS) and the impurities in the background plasma in fragment penetration, assimilation, radiative cooling and MHD activity in D 2 SPI experiments. It is found that plasmoid drifts lead to an about 70% reduction of the central line-integrated density (compared to a simulation without drifts) in the JET D 2 SPI discharge considered. Impurities that pre-exist before the SPI as well as those from possible impurity influxes related to the SPI are shown to dominate the radiation in the considered discharge. With inputs from JOREK simulations, modelling with the Lagrangian particle-based pellet code PELOTON reproduces the deviation of the SPI fragments in the direction of the major radius as observed by the fast camera. This confirms the role of rocket effects and plasmoid drifts in the considered discharge and reinforces the validity of the JOREK modelling. The limited core density rise due to plasmoid drifts and the strong radiative cooling and MHD activity with impurities (depending on their species and concentration) could limit the effectiveness of LFS D 2 SPI in runaway electron avoidance and are worth considering in the design of the ITER disruption mitigation system.

Physics↗

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 ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Benign termination of runaway electron beams on ASDEX Upgrade and TCV

Abstract This paper discusses the development of a benign termination scenario for runaway electron (RE) beams on ASDEX Upgrade and TCV. A systematic study revealed that a low electron density (n e ) companion plasma was required to achieve a large MHD instability, which expelled the confined REs over a large wetted area and allowed for the conversion of magnetic energy to radiation. Control of the companion plasma n e was achieved via neutral pressure regulation and was agnostic to material injection method. The neutral pressure required for recombination was found to be dependent on impurity species, quantity and RE current. On TCV,n e increased at neutral pressures above 1 Pa, indicating that higher collisionality between the REs and neutrals may lead to an upper pressure limit. The conversion of magnetic energy to radiated energy was measured on both machines and a decrease in efficiency was observed at high neutral pressure on TCV. The benign termination technique was able to prevent any significant increase in maximum heat flux on AUG from 200 to 600 kA of RE current, highlighting the ability of this approach to handle fully formed RE beams.

Physics↗

The effect of plasma shaping on high density H-mode SOL profiles and fluctuations in TCV

Abstract The impact of plasma shaping on the properties of high density H-mode scrape-off layer (SOL) profiles and transport at the outer midplane has been investigated on Tokamakà configuration variable. The experimental dataset has been acquired by evolving the upper triangularity while keeping the other parameters constant. The scan comprises δ u p values between 0.0 and 0.6, excluding negative triangularity scenarios. Within this study, a transition from type-I edge localised modes to the quasi-continuous exhaust regime takes place from low to high δ u p . The modification of the upstream SOL profiles has been assessed, in terms of separatrix quantities, within the α t turbulence control parameter theoretical framework (Eich et al 2020 Nucl. Fusion 60 056016). The target parallel heat load and the upstream near-SOL density profiles have been shown to broaden significantly for increasing α t . Correspondingly, in the far SOL a density shoulder formation is observed when moving from low to high δ u p . These behaviours have been correlated with an enhancement of the SOL fluctuation level, as registered by wall-mounted Langmuir probes as well as the thermal helium beam diagnostic. Specifically, both the background and the filamentary-induced fluctuating parts of the first wall ion saturation current signal are larger at higher δ u p , with filaments being ejected more frequently into the SOL. Comparison of two pulses at the extremes of the δ u p scan range, but with otherwise same input parameters, shows that the midplane neutral pressure does not change much during the H-mode phase of the discharge. This indicates that indirect effects of the change in geometry, linked to first wall recycling sources, should not play a significant role. The total core radiation increases at high δ u p , on account of a stronger plasma–wall interaction and resulting larger carbon impurity intake from the first wall. This is likely associated to the enhanced first wall fluctuations, as well as a smaller outer gap and the close-to-double-null magnetic topology at high shaping.

Stagni, A. (ORCID:0000000180841544)↗

MHD modeling of shattered pellet injection in JET

Abstract Nonlinear 3D MHD simulations of shattered-pellet injection (SPI) in JET show prototypical SPI-driven disruptions using the M3D-C1 and NIMROD extended-MHD codes. Initially, radiation-driven thermal quenches are accelerated by MHD activity as the pellet crosses rational surfaces, leading to a radiation spike, global stochasticization of the magnetic field, and a complete thermal quench. Eventually, current quenches, preceded by a current spike are seen as the Ohmic heating becomes equal to the radiative cooling. The results are qualitatively similar for both a single monolithic pellet, pencil-beam model, and a realistic shatter to represent the SPI plume. A scan in viscosity from 500 to 2000 m 2 s −1 for MHD simulations finds that reducing viscosity increases MHD activity and decreases thermal quench time slightly. A realistic cloud of fragments modeling shows that mixed-D–Ne pellet travels deeper into the plasma core before the thermal quench. At the slow pellet speeds, the pellet is found to be moving slowly enough inward that even the 5% neon in the mixed pellet is enough to effectively radiate the thermal energy available. Radiation toroidal peaking is predicted to be at levels consistent with experimental observations and reduced as the pellet travels deeper into the plasma. These simulations lay the ground work for more-sophisticated validative and predictive modeling of SPI in JET using both M3D-C1 and NIMROD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗