Engineering topics

Nardon, E.

Publications and source records attributed to Nardon, E..

5 records

DOE OSTI2024

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Developing high performance RF heating scenarios on the WEST tokamak

High power experiments, up to 9.2 MW with LHCD and ICRH, have been carried out in the full tungsten tokamak WEST. Here, quasi non inductive discharges have been achieved allowing to extend the plasma duration to 53 s with stationary conditions in particular with respect to tungsten contamination. Transitions to H mode are observed, and H-modes lasting up to 4 s have been obtained. The increase in stored energy is weak since the power crossing the separatrix is close to the L–H threshold. Hot L mode plasmas (central temperature exceeding 3 keV) with a confinement time following the ITER L96 scaling law are routinely obtained. The weak aspect ratio dependence of this scaling law is confirmed. Tungsten accumulation is generally not an operational issue on WEST. Difficulty of burning through tungsten can prevent the discharge from accessing to a hot core plasma in the ramp-up phase, or can lead to rapid collapse of the central temperature when radiation is enhanced by a slight decrease of the temperature. Except a few pulses post-boronization, the plasma radiation is rather high (P rad /P tot ~ 50%) and is dominated by tungsten. This fraction does not vary as the RF power is ramped up and is quite similar in ICRH and/or LHCD heated plasmas. An estimate of the contribution of the RF antennas to the plasma contamination in tungsten is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Developing high performance RF heating scenarios on the WEST tokamak

High power experiments, up to 9.2 MW with LHCD and ICRH, have been carried out in the full tungsten tokamak WEST. Quasi non inductive discharges have been achieved allowing to extend the plasma duration to 53s with stationary conditions in particular with respect to tungsten contamination. Transitions in H mode are obtained lasting up to 4s with weak energy increment at the power crossing the separatrix is close to the threshold. Hot L mode plasmas (Te(0)>3keV) with a confinement time following the ITER L96 scaling are routinely obtained. The weak aspect ratio dependence of this scaling law is confirmed. Tungsten accumulation is generally not an operational issue on WEST. Difficulty of burning through tungsten can prevent from accessing to a hot core plasma in the ramp-up phase or can lead to rapid collapse of the central temperature when radiation is enhanced by a slight decrease of the temperature. Apart few pulses post-boronization, the plasma radiation is rather high (Prad/Ptot~50%) and is dominated by tungsten. This fraction does not vary as the RF power is ramped up and is quite similar in ICRH and/or LHCD heated plasmas. An estimate of the contribution of the RF antennas to the plasma contamination in tungsten is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Disruption thermal load mitigation with shattered pellet injection on the Joint European Torus (JET)

Disruption mitigation remains a critical, unresolved challenge for ITER. To aid in addressing this challenge, a shattered pellet injection (SPI) system was installed on JET and experiments conducted at a range of thermal energy fractions and stored energies in excess of 7 MJ. The primary goals of these experiments were to investigate the efficacy of the SPI on JET and the ability of the plasma to assimilate multiple pellets. Single pellet injections produced a saturation in total radiated energy (Wrad) with increasing injected neon content, suggesting total radiation of stored thermal energy. Further increases in injected neon quantities resulted in reduced cooling times and current quench (CQ) durations, indicating higher impurity assimilation. No significant variation in CQ duration or W rad was observed when varying the deuterium content at fixed neon quantities. Additionally, higher assimilation, inferred by shorter CQ durations, was measured when a mechanical punch was used to launch the pellets and this was attributed to a lower pellet velocity leading to higher solid content in the pellet plume and larger fragments penetrating deeper into the plasma. Radiation asymmetries averaged over the cooling time were inferred from Emis3D and ranged from 1.6 to 1.9. Asymmetries averaged over the entire disruption sequence were found to increase at higher thermal energy fractions. The radiated energy fractions decreased with increasing thermal energy fractions but this trend was eliminated when toroidal asymmetries were accounted for with Emis3D. Pure deuterium pellets were able to produce cooling times of up to 75 ms with a gradual loss in thermal stored energy of up to 80%. Experiments with multiple pellet injection indicated W rad can be increased through pellet superposition and density can be increased with an additional D2 injection without a reduction in W rad . KPRAD modelling accurately reproduced the cooling times and the CQ duration at high thermal energies. Assimilation estimates from KPRAD indicated CQ rates scale strongly whilst W rad scales weakly and saturates with assimilated neon content. Comparable W rad can be achieved with lower assimilated neon quantities as longer cooling times are attained. Thus reduced neon content can be preferential in a thermal load mitigation scheme as it may reduce radiation asymmetries and prevent flash melting.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗