-
Runaway electron beam formation, vertical motion, termination and wall loads in EU-DEMO
Authors:
F. Vannini,
V. Bandaru,
H. Bergstroem,
N. Schwarz,
F. J. Artola,
M. Hoelzl,
G. Pautasso,
E. Nardon,
F. Maviglia,
M. L. Richiusa,
E. Emanuelli,
the JOREK team
Abstract:
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…
▽ More
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.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
3D MHD modelling of plasmoid drift following massive material injection in a tokamak
Authors:
M. Kong,
E. Nardon,
D. Bonfiglio,
M. Hoelzl,
D. Hu,
the JOREK team,
the EUROfusion Tokamak Exploitation Team
Abstract:
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 on the low field side midplane of a JET H-mode plasma to clarify basic processes and compare with existing theory. The simulations confirm the important role of the propagation of shear Alfvén wave (SAW) packets from both ends…
▽ More
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 on the low field side midplane of a JET H-mode plasma to clarify basic processes and compare with existing theory. The simulations confirm the important role of the propagation of shear Alfvén wave (SAW) packets from both ends of the plasmoid ("SAW braking") and the development of external resistive currents along magnetic field lines ("Pégourié 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égourié braking is shown to set in earlier with larger toroidal extent of the source, both in good agreement with existing theory. 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 also discussed.
△ Less
Submitted 18 July, 2024; v1 submitted 1 July, 2024;
originally announced July 2024.
-
Plasmoid drift and first wall heat deposition during ITER H-mode dual-SPIs in JOREK simulations
Authors:
D. Hu,
F. J. Artola,
E. Nardon,
M. Lehnen,
M. Kong,
D. Bonfiglio,
M. Hoelzl,
G. T. A. Huijsmans,
JOREK Team
Abstract:
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 "degr…
▽ More
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.
△ Less
Submitted 7 June, 2024;
originally announced June 2024.
-
Runaway Electron Dynamics in ITER Disruptions with Shattered Pellet Injections
Authors:
Oskar Vallhagen,
Lise Hanebring,
Javier Artola,
Michael Lehnen,
Eric Nardon,
Tünde Fülöp,
Mathias Hoppe,
Sarah Newton,
Istvan Pusztai
Abstract:
This study systematically explores the parameter space of disruption mitigation through shattered pellet injection in ITER with a focus on runaway electron dynamics, using the disruption modelling tool DREAM. The physics fidelity is considerably increased compared to previous studies, by e.g., using realistic magnetic geometry, resistive wall configuration, thermal quench onset criteria, as well a…
▽ More
This study systematically explores the parameter space of disruption mitigation through shattered pellet injection in ITER with a focus on runaway electron dynamics, using the disruption modelling tool DREAM. The physics fidelity is considerably increased compared to previous studies, by e.g., using realistic magnetic geometry, resistive wall configuration, thermal quench onset criteria, as well as including additional effects, such as ion transport and enhanced runaway electron transport during the thermal quench. The work aims to provide a fairly comprehensive coverage of experimentally feasible scenarios, considering plasmas representative of both non-activated and high-performance DT operation, different thermal quench onset criteria and transport levels, a wide range of hydrogen and neon quantities injected in one or two stages, and pellets with various characteristic shard sizes. Using a staggered injection scheme, with a pure hydrogen injection preceding a mixed hydrogen-neon injection, we find injection parameters leading to acceptable runaway electron currents in all investigated discharges without activated runaway sources. Dividing the injection into two stages is found to significantly enhance the assimilation and minimize runaway electron generation due to the hot-tail mechanism. However, while a staggered injection outperforms a single stage injection also in cases with radioactive runaway electron sources, no cases with acceptable runaway electron currents are found for a DT-plasma with a 15 MA plasma current.
△ Less
Submitted 25 January, 2024;
originally announced January 2024.
-
Lagrangian particle simulation of hydrogen pellets and SPI into runaway electron beam in ITER
Authors:
Shaohua Yuan,
Nizar Naitlho,
Roman Samulyak,
Bernard Pégourié,
Eric Nardon,
Eric Hollmann,
Paul Parks,
Michael Lehnen
Abstract:
Numerical studies of the ablation of pellets and shattered pellet injection (SPI) fragments into a runaway electron beam in ITER have been performed using a time-dependent pellet ablation code [R. Samulyak at el., Nucl Fusion, 61 (4), 046007 (2021)]. The code resolves detailed ablation physics near pellet fragments and large-scale expansion of ablated clouds. The study of a single fragment ablatio…
▽ More
Numerical studies of the ablation of pellets and shattered pellet injection (SPI) fragments into a runaway electron beam in ITER have been performed using a time-dependent pellet ablation code [R. Samulyak at el., Nucl Fusion, 61 (4), 046007 (2021)]. The code resolves detailed ablation physics near pellet fragments and large-scale expansion of ablated clouds. The study of a single fragment ablation quantifies the influence of various factors, in particular the impact ionization by runaway electrons and cross-field transport models, on the dynamics of ablated plasma and its penetration into the runaway beam. Simulations of SPI performed using different numbers of pellet fragments study the formation and evolution of ablation clouds and their large-scale dynamics in ITER. The penetration depth of ablation clouds is found to be of the order of 50 cm.
△ Less
Submitted 20 July, 2022;
originally announced July 2022.
-
Collisional-radiative non-equilibrium impurity treatment for JOREK simulations
Authors:
D. Hu,
G. T. A. Huijsmans,
E. Nardon,
M. Hoelzl,
M. Lehnen,
D. Bonfiglio,
JOREK Team
Abstract:
A collisional-radiative non-equilibrium impurity treatment for JOREK 3D nonlinear magneto-hydrodynamic (MHD) simulations has been developed. The impurities are represented by super-particles flowing along the fluid velocity field lines, while ionizing and recombining independently according to ADAS data and local fluid density and temperature. The non-equilibrium impurity contributions are then pr…
▽ More
A collisional-radiative non-equilibrium impurity treatment for JOREK 3D nonlinear magneto-hydrodynamic (MHD) simulations has been developed. The impurities are represented by super-particles flowing along the fluid velocity field lines, while ionizing and recombining independently according to ADAS data and local fluid density and temperature. The non-equilibrium impurity contributions are then projected back to the fluid field for self-consistent time evolution. A 2D test case is used to compare the new non-equilibrium impurity model against previous Coronal Equilibrium (CE) impurity treatment, as well as to compare the non-equilibrium impurity behavior between the single and the two temperature model. Further, we conduct benchmark with previously published coronal non-equilibrium results by other 3D nonlinear MHD codes such as M3D-C1 and NIMROD. The new non-equilibrium treatment is shown to successfully capture the early phase cooling by weakly ionized impurities which the CE model missed. The benchmarks with M3D-C1 and NIMROD show general agreement in both the integrated quantities and the 2D profile evolution, despite the difference in the atomic model used. The above comparison and benchmark cases demonstrate the capability of the non-equilibrium impurity model for JOREK, paving the way for more sophisticated 3D non-linear Massive Material Injection (MMI) simulations which have important applications in disruption mitigation studies.
△ Less
Submitted 24 June, 2021;
originally announced June 2021.
-
Latest results on quiescent and post-disruption runaway electron mitigation experiments at Frascati Tokamak Upgrade
Authors:
D. Carnevale,
P. Buratti,
M. Baruzzo,
W. Bin,
F. Bombarda,
L. Boncagni,
C. Paz-Soldan,
L. Calacci,
M. Cappelli,
C. Castaldo,
S. Ceccuzzi,
C. Centioli,
C. Cianfarani,
S. Coda,
F. Cordella,
O. D Arcangelo,
J. Decker,
B. Duval,
B. Esposito,
L. Gabellieri,
S. Galeani,
S. Garavaglia,
C. Galperti,
G. Ghillardi,
G. Granucci
, et al. (16 additional authors not shown)
Abstract:
Results from the last FTU campaigns on the deuterium large (wrt FTU volume) pellet REs suppression capability, mainly due to the induced burst MHD activity expelling REs seed are presented for discharges with 0.5 MA and 5.3T. Clear indications of avalanche multiplication of REs following single pellet injection on 0.36 MA flat-top discharges is shown together with quantitative indications of dissi…
▽ More
Results from the last FTU campaigns on the deuterium large (wrt FTU volume) pellet REs suppression capability, mainly due to the induced burst MHD activity expelling REs seed are presented for discharges with 0.5 MA and 5.3T. Clear indications of avalanche multiplication of REs following single pellet injection on 0.36 MA flat-top discharges is shown together with quantitative indications of dissipative effects in terms of critical electrical field increase due to fan-like instabilities. Analysis of large fan-like instabilities on post-disruption RE beams, that seem to be correlated with low electrical field and background density drops, reveal their strong RE energy suppression capability suggesting a new strategy for RE energy suppression controlling large fan instabilities. We demonstrate how such density drops can be induced using modulated ECRH power on post-disruption beams.
△ Less
Submitted 25 May, 2021; v1 submitted 10 May, 2021;
originally announced May 2021.
-
The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas
Authors:
M Hoelzl,
GTA Huijsmans,
SJP Pamela,
M Becoulet,
E Nardon,
FJ Artola,
B Nkonga,
CV Atanasiu,
V Bandaru,
A Bhole,
D Bonfiglio,
A Cathey,
O Czarny,
A Dvornova,
T Feher,
A Fil,
E Franck,
S Futatani,
M Gruca,
H Guillard,
JW Haverkort,
I Holod,
D Hu,
SK Kim,
SQ Korving
, et al. (28 additional authors not shown)
Abstract:
JOREK is a massively parallel fully implicit non-linear extended MHD code for realistic tokamak X-point plasmas. It has become a widely used versatile code for studying large-scale plasma instabilities and their control developed in an international community. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and phys…
▽ More
JOREK is a massively parallel fully implicit non-linear extended MHD code for realistic tokamak X-point plasmas. It has become a widely used versatile code for studying large-scale plasma instabilities and their control developed in an international community. 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 and current quench triggered by massive gas injection (MGI) and shattered pellet injection (SPI), runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events (VDEs). Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring thermal quenches triggered by locked modes, and radiative collapses are being studied.
△ Less
Submitted 21 April, 2021; v1 submitted 18 November, 2020;
originally announced November 2020.
-
Radiation asymmetry and MHD destabilization during the thermal quench after impurity Shattered Pellet Injection
Authors:
D. Hu,
E. Nardon,
M. Hoelzl,
F. Wieschollek,
M. Lehnen,
G. T. A. Huijsmans,
D. C. van Vugt,
S. -H. Kim,
JET contributors,
JOREK Team
Abstract:
The radiation response and the MHD destabilization during the thermal quench after a mixed species Shattered Pellet Injection (SPI) with impurity species neon and argon are investigated via 3D non-linear MHD simulation using the JOREK code. Both the $n=0$ global current profile contraction and the local helical cooling at each rational surface caused by the pellet fragments are found to be respons…
▽ More
The radiation response and the MHD destabilization during the thermal quench after a mixed species Shattered Pellet Injection (SPI) with impurity species neon and argon are investigated via 3D non-linear MHD simulation using the JOREK code. Both the $n=0$ global current profile contraction and the local helical cooling at each rational surface caused by the pellet fragments are found to be responsible for MHD destabilization after the injection. Significant current driven mode growth is observed as the fragments cross low order rational surfaces, resulting in rapidly inward propagating stochastic magnetic field, ultimately causing the core temperature collapse. The Thermal Quench (TQ) is triggered as the fragments arrive on the $q=1$ or $q=2$ surface depending on the exact $q$ profile and thus mode structure. When injecting from a single toroidal location, strong radiation asymmetry is found before and during the TQ as a result of the unrelaxed impurity density profile along the field line and asymmetric outward heat flux. Such asymmetry gradually relaxes over the course of the TQ, and is entirely eliminated by the end of it. Simulation results indicate that the aforementioned asymmetric radiation behavior could be significantly mitigated by injection from toroidally opposite locations, provided that the time delay between the two injectors is shorter than $1ms$. It is also found that the MHD response are sensitive to the relative timing and injection configuration in these multiple injection cases.
△ Less
Submitted 25 November, 2020; v1 submitted 6 September, 2020;
originally announced September 2020.
-
Post-Thermal Quench Shattered Pellet Injection for small Runaway Electron seed depletion in ITER
Authors:
E. Nardon,
A. Matsuyama,
D. Hu,
F. Wieschollek
Abstract:
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 Cur…
▽ More
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.
△ Less
Submitted 3 July, 2020;
originally announced July 2020.
-
Fast plasma dilution in ITER with pure Deuterium Shattered Pellet Injection
Authors:
E. Nardon,
D. Hu,
M. Hoelzl,
D. Bonfiglio
Abstract:
JOREK 3D non-linear MagnetoHydroDynamic (MHD) simulations of pure Deuterium Shattered Pellet Injection in ITER are presented. It is shown that such a scheme could allow diluting the plasma by more than a factor 10 without immediately triggering large MHD activity, provided the background impurity density is low enough. This appears as a promising strategy to reduce the risk of hot tail Runaway Ele…
▽ More
JOREK 3D non-linear MagnetoHydroDynamic (MHD) simulations of pure Deuterium Shattered Pellet Injection in ITER are presented. It is shown that such a scheme could allow diluting the plasma by more than a factor 10 without immediately triggering large MHD activity, provided the background impurity density is low enough. This appears as a promising strategy to reduce the risk of hot tail Runaway Electron (RE) generation and to avoid RE beams altogether in ITER, motivating further studies in this direction.
△ Less
Submitted 14 September, 2020; v1 submitted 29 June, 2020;
originally announced June 2020.
-
Assessing energy dependence of the transport of relativistic electrons in perturbed magnetic fields with orbit-following simulations
Authors:
Konsta Särkimäki,
Ola Embreus,
Eric Nardon,
Tünde Fülöp,
JET Contributors
Abstract:
Experimental observations, as well as theoretical predictions, indicate that the transport of energetic electrons decreases with energy. This reduction in transport is attributed to finite orbit width (FOW) effects. Using orbit-following simulations in perturbed tokamak magnetic fields that have an ideal homogeneous stochastic layer at the edge, we quantify the energy dependence of energetic elect…
▽ More
Experimental observations, as well as theoretical predictions, indicate that the transport of energetic electrons decreases with energy. This reduction in transport is attributed to finite orbit width (FOW) effects. Using orbit-following simulations in perturbed tokamak magnetic fields that have an ideal homogeneous stochastic layer at the edge, we quantify the energy dependence of energetic electrons transport and confirm previous theoretical estimates. However, using magnetic configurations characteristic of JET disruptions, we find no reduction in RE transport at higher energies, which we attribute to the mode widths being comparable to the minor radius, making the FOW effects negligible. Instead, the presence of islands and nonuniform magnetic perturbations are found to be more important. The diffusive-advective transport coefficients calculated in this work, based on simulations for electron energies 10 keV -- 100 MeV, can be used in reduced kinetic models to account for the transport due to the magnetic field perturbations.
△ Less
Submitted 5 June, 2020;
originally announced June 2020.
-
Understanding the reduction of the edge safety factor during hot VDEs and fast edge cooling events
Authors:
F. J. Artola,
K. Lackner,
G. T. A. Huijsmans,
M. Hoelzl,
E. Nardon,
A. Loarte
Abstract:
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, whic…
▽ More
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.
△ Less
Submitted 9 March, 2020;
originally announced March 2020.
-
First predictive simulations for deuterium shattered pellet injection in ASDEX Upgrade
Authors:
M. Hoelzl,
D. Hu,
E. Nardon,
G. T. A. Huijsmans
Abstract:
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 influen…
▽ More
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.
△ Less
Submitted 29 December, 2019; v1 submitted 14 October, 2019;
originally announced October 2019.
-
Electron acceleration in a JET disruption simulation
Authors:
Cristian Sommariva,
Eric Nardon,
Peter Beyer,
Matthias Hoelzl,
Guido Huijsman
Abstract:
Runaways are suprathermal electrons having sufficiently high energy to be continuously accelerated up to tens of MeV by a driving electric field [1]. Highly energetic runaway electron (RE) beams capable of damaging the tokamak first wall can be observed after a plasma disruption [2]. Therefore, it is of primary importance to fully understand their generation mechanisms in order to design mitigatio…
▽ More
Runaways are suprathermal electrons having sufficiently high energy to be continuously accelerated up to tens of MeV by a driving electric field [1]. Highly energetic runaway electron (RE) beams capable of damaging the tokamak first wall can be observed after a plasma disruption [2]. Therefore, it is of primary importance to fully understand their generation mechanisms in order to design mitigation systems able to guarantee safe tokamak operations. In a previous work, [3], a test particle tracker was introduced in the JOREK 3D non-linear MHD code and used for studying the electron confinement during a simulated JET-like disruption. It was found in [3] that relativistic electrons are not completely deconfined by the stochastic magnetic field taking place during the disruption thermal quench (TQ). This is due to the reformation of closed magnetic surfaces at the beginning of the current quench (CQ). This result was obtained neglecting the inductive electric field in order to avoid the unrealistic particle acceleration which otherwise would have happened due to the absence of collision effects. The present paper extends [3] analysing test electron dynamics in the same simulated JET-like disruption using the complete electric field. For doing so, a simplified collision model is introduced in the particle tracker guiding center equations. We show that electrons at thermal energies can become RE during or promptly after the TQ due to a combination of three phenomena: a first REs acceleration during the TQ due to the presence of a complex MHD-induced electric field, particle reconfinement caused by the fast reformation of closed magnetic surfaces after the TQ and a secondary acceleration induced by the CQ electric field.
△ Less
Submitted 13 August, 2018; v1 submitted 15 May, 2018;
originally announced May 2018.
-
Test particles dynamics in the JOREK 3D non-linear MHD code and application to electron transport in a disruption simulation
Authors:
Cristian Sommariva,
Eric Nardon,
Matthias Hoelzl,
Guido Huijsmans,
Daan van Vugt
Abstract:
In order to contribute to the understanding of runaway electron generation mechanisms during tokamak disruptions, a test particle tracker is introduced in the JOREK 3D non-linear MHD code, able to compute both full and guiding center relativistic orbits. Tests of the module show good conservation of the invariants of motion and consistency between full orbit and guiding center solutions. A first a…
▽ More
In order to contribute to the understanding of runaway electron generation mechanisms during tokamak disruptions, a test particle tracker is introduced in the JOREK 3D non-linear MHD code, able to compute both full and guiding center relativistic orbits. Tests of the module show good conservation of the invariants of motion and consistency between full orbit and guiding center solutions. A first application is presented where test electron confinement properties in a JET massive gas injection-triggered disruption simulation are investigated. It is found that electron populations initialised before the thermal quench (TQ) are typically not fully deconfined in spite of the global stochasticity of the magnetic field during the TQ. The fraction of "survivors" decreases from few tens down to a few tenths of percent as the electron energy varies from 1keV to 10MeV. The underlying mechanism for electron "survival" is the prompt reformation of closed magnetic surfaces at the plasma core and, to a smaller extent, the subsequent reappearance of a magnetic surface at the edge. It is also found that electrons are less deconfined at 10MeV than at 1MeV which appears consistent with a phase averaging effect due to orbit shifts at high energy.
△ Less
Submitted 12 September, 2017; v1 submitted 28 April, 2017;
originally announced April 2017.
-
Non-linear Simulations of MHD Instabilities in Tokamaks Including Eddy Current Effects and Perspectives for the Extension to Halo Currents
Authors:
M. Hoelzl,
G. T. A. Huijsmans,
P. Merkel,
C. Atanasiu,
K. Lackner,
E. Nardon,
K. Aleynikova,
F. Liu,
E. Strumberger,
R. McAdams,
I. Chapman,
A. Fil
Abstract:
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)…
▽ More
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.
△ Less
Submitted 23 September, 2014; v1 submitted 27 August, 2014;
originally announced August 2014.
-
Self-consistent simulation of plasma scenarios for ITER using a combination of 1.5D transport codes and free-boundary equilibrium codes
Authors:
V Parail,
R Albanese,
R Ambrosino,
J-F Artaud,
K Besseghir,
M Cavinato,
G Corrigan,
J Garcia,
L Garzotti,
Y Gribov,
F Imbeaux,
F Koechl,
C V Labate,
J Lister,
X Litaudon,
A Loarte,
P Maget,
M Mattei,
D McDonald,
E Nardon,
G Saibene,
R Sartori,
J Urban
Abstract:
Self-consistent transport simulation of ITER scenarios is a very important tool for the exploration of the operational space and for scenario optimisation. It also provides an assessment of the compatibility of developed scenarios (which include fast transient events) with machine constraints, in particular with the poloidal field (PF) coil system, heating and current drive (H&CD), fuelling and pa…
▽ More
Self-consistent transport simulation of ITER scenarios is a very important tool for the exploration of the operational space and for scenario optimisation. It also provides an assessment of the compatibility of developed scenarios (which include fast transient events) with machine constraints, in particular with the poloidal field (PF) coil system, heating and current drive (H&CD), fuelling and particle and energy exhaust systems. This paper discusses results of predictive modelling of all reference ITER scenarios and variants using two suite of linked transport and equilibrium codes. The first suite consisting of the 1.5D core/2D SOL code JINTRAC [1] and the free boundary equilibrium evolution code CREATE-NL [2,3], was mainly used to simulate the inductive D-T reference Scenario-2 with fusion gain Q=10 and its variants in H, D and He (including ITER scenarios with reduced current and toroidal field). The second suite of codes was used mainly for the modelling of hybrid and steady state ITER scenarios. It combines the 1.5D core transport code CRONOS [4] and the free boundary equilibrium evolution code DINA-CH [5].
△ Less
Submitted 31 October, 2013;
originally announced October 2013.
-
Characteristics of X-point Lobe Structures in Single-Null Discharges on MAST
Authors:
J. R. Harrison,
A. Kirk,
I. T. Chapman,
P. Cahyna,
Yueqiang Liu,
E. Nardon,
A. J. Thornton
Abstract:
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…
▽ More
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.
△ Less
Submitted 5 July, 2013;
originally announced July 2013.
-
Understanding ELM mitigation by resonant magnetic perturbations on MAST
Authors:
A. Kirk,
I. T. Chapman,
Yueqiang Liu,
P. Cahyna,
P. Denner,
G. Fishpool,
C. J. Ham,
J. R. Harrison,
Yunfeng Liang,
E. Nardon,
S. Saarelma,
R. Scannell,
A. J. Thornton,
the MAST team
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 June, 2013;
originally announced June 2013.
-
Effect of resonant magnetic perturbations with toroidal mode numbers of 4 and 6 on ELMs in single null H-mode plasmas in MAST
Authors:
A. Kirk,
I. T. Chapman,
J. Harrison,
Yueqiang Liu,
E. Nardon,
S. Saarelma,
R. Scannell,
A. J. Thornton,
the MAST team
Abstract:
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 line…
▽ More
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.
△ Less
Submitted 14 May, 2013;
originally announced May 2013.
-
Effect of resonant magnetic perturbations on ELMs in connected double null plasmas in MAST
Authors:
A. Kirk,
Yueqiang Liu,
I. T. Chapman,
J. Harrison,
E. Nardon,
R. Scannell,
A. J. Thornton,
the MAST team
Abstract:
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 linearl…
▽ More
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.
△ Less
Submitted 1 March, 2013;
originally announced March 2013.
-
Coupling JOREK and STARWALL for Non-linear Resistive-wall Simulations
Authors:
M. Hoelzl,
P. Merkel,
G. T. A. Huysmans,
E. Nardon,
E. Strumberger,
R. McAdams,
I. Chapman,
S. Guenter,
K. Lackner
Abstract:
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…
▽ More
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.
△ Less
Submitted 4 October, 2012; v1 submitted 13 June, 2012;
originally announced June 2012.
-
Observation of lobes near the X-point in resonant magnetic perturbation experiments on MAST
Authors:
A. Kirk,
J. Harrison,
Yueqiang Liu,
E. Nardon,
I. T. Chapman,
P. Denner,
the MAST team
Abstract:
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…
▽ More
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.
△ Less
Submitted 1 February, 2012;
originally announced February 2012.
-
Model for screening of resonant magnetic perturbations by plasma in a realistic tokamak geometry and its impact on divertor strike points
Authors:
Pavel Cahyna,
Eric Nardon
Abstract:
This work addresses the question of the relation between strike-point splitting and magnetic stochasticity at the edge of a poloidally diverted tokamak in the presence of externally imposed magnetic perturbations. More specifically, ad-hoc helical current sheets are introduced in order to mimic a hypothetical screening of the external resonant magnetic perturbations by the plasma. These current sh…
▽ More
This work addresses the question of the relation between strike-point splitting and magnetic stochasticity at the edge of a poloidally diverted tokamak in the presence of externally imposed magnetic perturbations. More specifically, ad-hoc helical current sheets are introduced in order to mimic a hypothetical screening of the external resonant magnetic perturbations by the plasma. These current sheets, which suppress magnetic islands, are found to reduce the amount of splitting expected at the target, which suggests that screening effects should be observable experimentally. Multiple screening current sheets reinforce each other, i.e. less current relative to the case of only one current sheet is required to screen the perturbation.
△ Less
Submitted 27 January, 2011; v1 submitted 17 December, 2010;
originally announced December 2010.
-
Resonant magnetic perturbations and divertor footprints in poloidally diverted tokamaks
Authors:
Pavel Cahyna,
Eric Nardon
Abstract:
General formula describing both the divertor strike point splitting and width of magnetic islands created by resonant magnetic perturbations (RMPs) in a poloidally diverted tokamak equilibrium is derived. Under the assumption that the RMP is produced by coils at the low-field side such as those used to control edge localized modes (ELMs) it is demonstrated that the width of islands on different ma…
▽ More
General formula describing both the divertor strike point splitting and width of magnetic islands created by resonant magnetic perturbations (RMPs) in a poloidally diverted tokamak equilibrium is derived. Under the assumption that the RMP is produced by coils at the low-field side such as those used to control edge localized modes (ELMs) it is demonstrated that the width of islands on different magnetic surfaces at the edge and the amount of divertor splitting are related to each other. Explanation is provided of aligned maxima of the perturbation spectra with the safety factor profile - an effect empirically observed in models of many perturbation coil designs.
△ Less
Submitted 20 May, 2010;
originally announced May 2010.