Source author record

E. Nardon

E. Nardon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2025arXiv

Runaway electron avalanche and macroscopic beam formation: simulations of the DTT full power scenario

The transition of the Divertor Tokamak Test (DTT) facility from its initial commissioning phase (Day-0, plasma current $I_{p}=2$ MA) to the full power scenario ($I_{p}=5.5$ MA) introduces a critical shift in the dynamics of runaway electrons (REs) generation. While previous predictive studies of the low-current scenario indicated a robust safety margin against RE beam formation, this work reveals that the exponential scaling of the RE avalanche gain with plasma current severely narrows the safe operational window in the full power scenario. Using the non-linear magnetohydrodynamic code JOREK, we perform comprehensive 2D simulations of the current quench (CQ) phase of several disruption scenarios, systematically scanning initial RE seed currents and injected impurity levels. The results demonstrate that in the full power scenario, the avalanche multiplication factor is sufficiently high ($G_\text{av} \approx 1.3 \cdot 10^5$) to convert a mere 5.5 A seed current into macroscopic RE beams of $\approx 0.7$ MA when large amounts of impurities are present. For even higher RE seeds, the RE current can peak at $ \approx 3.2$ MA, constituting up to $\approx$ 80% of the total plasma current during the CQ. These findings suggest that, unlike the Day-0 phase, the disruption mitigation strategy for the full power scenario involves a careful balance between thermal load mitigation and RE avoidance, necessitating a well-chosen quantity of injected impurities. This work provides the baseline needed for future estimations of RE loads on the plasma-facing components of DTT, which will be essential for designing and positioning mitigation components like sacrificial limiters.

preprint2020arXiv

Post-Thermal Quench Shattered Pellet Injection for small Runaway Electron seed depletion in ITER

The possibility of using Shattered Pellet Injection(s) after the Thermal Quench phase of an ITER disruption in order to deplete Runaway Electron (RE) seeds before they can substantially avalanche is studied. Analytical and numerical estimates of the required injection rate for shards to penetrate into the forming RE beam and stop REs are given. How much material could be assimilated before the Current Quench (CQ) becomes too short is also estimated. It appears that, if Hydrogen pellets were used, the required number of pellets to be injected during the CQ would be prohibitive, at least considering the present design of the ITER Disruption Mitigation System (DMS). For Neon or Argon, the required number of pellets, although large, might be within reach of the ITER DMS, but the assimilated fraction would have to be very small. Other materials may be better suited but would require a modification of the ITER DMS.

preprint2020arXiv

Understanding the reduction of the edge safety factor during hot VDEs and fast edge cooling events

In the present work a simple analytical approach is presented in order to clarify the physics behind the edge current density behaviour of a hot plasma entering in contact with a resistive conductor. When a plasma enters in contact with a highly resistive wall, large current densities appear at the edge of the plasma. The model shows that this edge current originates from the plasma response, which attempts to conserve the poloidal magnetic flux ($Ψ$) when the outer current is being lost. The loss of outer current is caused by the high resistance of the outer current path compared to the plasma core resistance. The resistance of the outer path may be given by plasma contact with a very resistive structure or by a sudden decrease of the outer plasma temperature (e.g. due to a partial thermal quench or due to a cold front penetration caused by massive gas injection). For general plasma geometries and current density profiles the model shows that given a small change of minor radius ($δa$) the plasma current is conserved to first order ($δI_p = 0 + \mathcal{O}(δa^2)$). This conservation comes from the fact that total inductance remains constant ($δL = 0$) due to an exact compensation of the change of external inductance with the change of internal inductance ($δL_\text{ext}+δL_\text{int} = 0$). As the total current is conserved and the plasma volume is reduced, the edge safety factor drops according to $q_a \propto a^2/I_p$. Finally the consistency of the resulting analytical predictions is checked with the help of free-boundary MHD simulations.

preprint2019arXiv

First predictive simulations for deuterium shattered pellet injection in ASDEX Upgrade

First simulations of deuterium shattered pellet injection (SPI) into an ASDEX Upgrade H-Mode plasma with the JOREK MHD code are presented. Resistivity is increased by one order of magnitude in most simulations to reduce computational costs and allow for extensive parameter scans. The effect of various physical parameters onto MHD activity and thermal quench (TQ) dynamics is studied and the influence of MHD onto ablation is shown. TQs are obtained quickly after injection in most simulations with a typical duration of 100 microseconds, which slows down at lower resistivity. Although the n=1 magnetic perturbation dominates in the simulations, toroidal harmonics up to n=10 contribute to stochastization and stochastic transport in the plasma core. The post-TQ density profile remains hollow for a few hundred microseconds. However, when flux surfaces re-form around the magnetic axis, the density has become monotonic again suggesting a beneficial behaviour for runaway electron avoidance/mitigation. With $10^{21}$ atoms injected, the TQ is typically incomplete and triggered when the shards reach the q=2 rational surface. At a larger number of injected atoms, the TQ can set in even before the shards reach this surface. For low field side injection considered here, repeated formation of outward convection cells is observed in the ablation region reducing material assimilation. Responsible is a sudden rise of pressure in the high density cloud when the stochastic region expands further releasing heat from the hot core. After the TQ, strong sheared poloidal rotation is created by Maxwell stress, which contributes to re-formation of flux surfaces.

preprint2014arXiv

Non-linear Simulations of MHD Instabilities in Tokamaks Including Eddy Current Effects and Perspectives for the Extension to Halo Currents

The dynamics of large scale plasma instabilities can strongly be influenced by the mutual interaction with currents flowing in conducting vessel structures. Especially eddy currents caused by time-varying magnetic perturbations and halo currents flowing directly from the plasma into the walls are important. The relevance of a resistive wall model is directly evident for Resistive Wall Modes (RWMs) or Vertical Displacement Events (VDEs). However, also the linear and non-linear properties of most other large-scale instabilities may be influenced significantly by the interaction with currents in conducting structures near the plasma. The understanding of halo currents arising during disruptions and VDEs, which are a serious concern for ITER as they may lead to strong asymmetric forces on vessel structures, could also benefit strongly from these non-linear modeling capabilities. Modeling the plasma dynamics and its interaction with wall currents requires solving the magneto-hydrodynamic (MHD) equations in realistic toroidal X-point geometry consistently coupled with a model for the vacuum region and the resistive conducting structures. With this in mind, the non-linear finite element MHD code JOREK has been coupled with the resistive wall code STARWALL, which allows to include the effects of eddy currents in 3D conducting structures in non-linear MHD simulations. This article summarizes the capabilities of the coupled JOREK-STARWALL system and presents benchmark results as well as first applications to non-linear simulations of RWMs, VDEs, disruptions triggered by massive gas injection, and Quiescent H-Mode. As an outlook, the perspectives for extending the model to halo currents are described.

preprint2013arXiv

Characteristics of X-point Lobe Structures in Single-Null Discharges on MAST

Lobe structures due to the application of resonant magnetic perturbations (RMPs) have been observed using wide-angle imaging of light from He1+ ions in the vicinity of the lower X-point in MAST. The data presented are from lower single-null discharges where RMPs of toroidal mode number, n, of 4 and 6 were applied. It has been found that, above a threshold value, the lobe structures extend radially, linearly with the coil current, both in L-mode and H-mode. It is observed that after the application of the RMP, as the toroidal rotation in the confined plasma decreases, the lobes extend radially, suggesting the plasma is less effectively screening the RMP field. Comparing the imaging data with results from vacuum modelling shows that this technique can accurately predict the number and poloidal location of the lobes, but over-estimates their radial extent. More accurate estimates of the extent of the lobes can be made by accounting for plasma screening of the RMP field. Qualitative agreement between simulation and experiment is found if it is assumed that the RMP penetrates 2% in normalised radius from the last closed flux surface.

preprint2013arXiv

Effect of resonant magnetic perturbations on ELMs in connected double null plasmas in MAST

The application of resonant magnetic perturbations (RMPs) with a toroidal mode number of n=3 to connected double null plasmas in the MAST tokamak produces up to a factor of 9 increase in Edge Localized Mode (ELM) frequency and reduction in plasma energy loss associated with type-I ELMs. A threshold current for ELM mitigation is observed above which the ELM frequency increases approximately linearly with current in the coils. The effect of the RMPs is found to be scenario dependent. In one scenario the mitigation is only due to a large density pump out event and if the density is recovered by gas puffing a return to type I ELMs is observed. In another scenario sustained ELM mitigation can be achieved irrespective of the amount of fuelling. Despite a large scan of parameters complete ELM suppression has not been achieved. The results have been compared to modelling performed using either the vacuum approximation or including the plasma response. The requirement for a resonant condition, that is an optimum alignment of the perturbation with the plasma, has been confirmed by performing a scan in the pitch angle of the applied field.

preprint2013arXiv

Effect of resonant magnetic perturbations with toroidal mode numbers of 4 and 6 on ELMs in single null H-mode plasmas in MAST

The application of resonant magnetic perturbations (RMPs) with a toroidal mode number of n=4 or n=6 to lower single null plasmas in the MAST tokamak produces up to a factor of 5 increase in Edge Localized Mode (ELM) frequency and reduction in plasma energy loss associated with type-I ELMs. A threshold current for ELM mitigation is observed above which the ELM frequency increases approximately linearly with current in the coils. Despite a large scan of parameters, complete ELM suppression has not been achieved. The results have been compared to modelling performed using either the vacuum approximation or including the plasma response. During the ELM mitigated stage clear lobe structures are observed in visible-light imaging of the X-point region. The size of these lobes is correlated with the increase in ELM frequency observed. The characteristics of the mitigated ELMs are similar to those of the natural ELMs suggesting that they are type I ELMs which are triggered at a lower pressure gradient. The application of the RMPs in the n=4 and n=6 configurations before the L-H transition has little effect on the power required to achieve H-mode while still allowing the first ELM to be mitigated.

preprint2013arXiv

Understanding ELM mitigation by resonant magnetic perturbations on MAST

Sustained ELM mitigation has been achieved using RMPs with a toroidal mode number of n=4 and n=6 in lower single null and with n=3 in connected double null plasmas on MAST. The ELM frequency increases by up to a factor of eight with a similar reduction in ELM energy loss. A threshold current for ELM mitigation is observed above which the ELM frequency increases approximately linearly with current in the coils. A comparison of the filament structures observed during the ELMs in the natural and mitigated stages shows that the mitigated ELMs have the characteristics of type I ELMs even though their frequency is higher, their energy loss is reduced and the pedestal pressure gradient is decreased. During the ELM mitigated stage clear lobe structures are observed in visible-light imaging of the X-point region. The size of these lobes is correlated with the increase in ELM frequency observed. The RMPs produce a clear 3D distortion to the plasma and it is likely that these distortions explain why ELMs are destabilised and hence why ELM mitigation occurs.

preprint2012arXiv

Coupling JOREK and STARWALL for Non-linear Resistive-wall Simulations

The implementation of a resistive-wall extension to the non-linear MHD-code JOREK via a coupling to the vacuum-field code STARWALL is presented along with first applications and benchmark results. Also, non-linear saturation in the presence of a resistive wall is demonstrated. After completion of the ongoing verification process, this code extension will allow to perform non-linear simulations of MHD instabilities in the presence of three-dimensional resistive walls with holes for limited and X-point plasmas.

preprint2012arXiv

Observation of lobes near the X-point in resonant magnetic perturbation experiments on MAST

The application of non-axisymmetric resonant magnetic perturbations (RMPs) with a toroidal mode number n=6 in the MAST tokamak produces a significant reduction in plasma energy loss associated with type-I Edge Localized Modes (ELMs), the first such observation with n>3. During the ELM mitigated stage clear lobe structures are observed in visible-light imaging of the X-point region. These lobes or manifold structures, that were predicted previously, have been observed for the first time in a range of discharges and their appearance is correlated with the effect of RMPs on the plasma i.e. they only appear above a threshold when a density pump out is observed or when the ELM frequency is increased. They appear to be correlated with the RMPs penetrating the plasma and may be important in explaining why the ELM frequency increases. The number and location of the structures observed can be well described using vacuum modelling. Differences in radial extent and poloidal width from vacuum modelling are likely to be due to a combination of transport effects and plasma screening.