-
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.
-
Neural-Parareal: Dynamically Training Neural Operators as Coarse Solvers for Time-Parallelisation of Fusion MHD Simulations
Authors:
S. J. P. Pamela,
N. Carey,
J. Brandstetter,
R. Akers,
L. Zanisi,
J. Buchanan,
V. Gopakumar,
M. Hoelzl,
G. Huijsmans,
K. Pentland,
T. James,
G. Antonucci,
the JOREK Team
Abstract:
The fusion research facility ITER is currently being assembled to demonstrate that fusion can be used for industrial energy production, while several other programmes across the world are also moving forward, such as EU-DEMO, CFETR, SPARC and STEP. The high engineering complexity of a tokamak makes it an extremely challenging device to optimise, and test-based optimisation would be too slow and to…
▽ More
The fusion research facility ITER is currently being assembled to demonstrate that fusion can be used for industrial energy production, while several other programmes across the world are also moving forward, such as EU-DEMO, CFETR, SPARC and STEP. The high engineering complexity of a tokamak makes it an extremely challenging device to optimise, and test-based optimisation would be too slow and too costly. Instead, digital design and optimisation must be favored, which requires strongly-coupled suites of High-Performance Computing calculations. In this context, having surrogate models to provide quick estimates with uncertainty quantification is essential to explore and optimise new design options. Furthermore, these surrogates can in turn be used to accelerate simulations in the first place. This is the case of Parareal, a time-parallelisation method that can speed-up large HPC simulations, where the coarse-solver can be replaced by a surrogate. A novel framework, Neural-Parareal, is developed to integrate the training of neural operators dynamically as more data becomes available. For a given input-parameter domain, as more simulations are being run with Parareal, the large amount of data generated by the algorithm is used to train new surrogate models to be used as coarse-solvers for future Parareal simulations, leading to progressively more accurate coarse-solvers, and thus higher speed-up. It is found that such neural network surrogates can be much more effective than traditional coarse-solver in providing a speed-up with Parareal. This study is a demonstration of the convergence of HPC and AI which simply has to become common practice in the world of digital engineering design.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
Probing non-linear MHD stability of the EDA H-mode in ASDEX Upgrade
Authors:
A Cathey,
M Hoelzl,
L Gil,
MG Dunne,
GF Harrer,
GTA Huijsmans,
J Kalis,
K Lackner,
SJP Pamela,
E Wolfrum,
S Günter,
the JOREK team,
the ASDEX Upgrade Team,
the EUROfusion MST1 Team
Abstract:
Regimes of operation in tokamaks that are devoid of large ELMs have to be better understood to extrapolate their applicability to reactor-relevant devices. This paper describes non-linear extended MHD simulations that use an experimental equilibrium from an EDA H-mode in ASDEX Upgrade. Linear ideal MHD analysis indicates that the operational point lies slightly inside of the stable region. The non…
▽ More
Regimes of operation in tokamaks that are devoid of large ELMs have to be better understood to extrapolate their applicability to reactor-relevant devices. This paper describes non-linear extended MHD simulations that use an experimental equilibrium from an EDA H-mode in ASDEX Upgrade. Linear ideal MHD analysis indicates that the operational point lies slightly inside of the stable region. The non-linear simulations with the visco-resistive extended MHD code, JOREK, sustain non-axisymmetric perturbations that are linearly most unstable with toroidal mode numbers of n = \{6 \dots 9\}, but non-linearly higher and lower n become driven and the low-n become dominant. The poloidal mode velocity during the linear phase is found to correspond to the expected velocity for resistive ballooning modes. The perturbations that exist in the simulations have somewhat smaller poloidal wavenumbers (k_θ \sim 0.1 to 0.5 cm^{-1} ) than the experimental expectations for the quasi-coherent mode in EDA, and cause non-negligible transport in both the heat and particle channels. In the transition from linear to non-linear phase, the mode frequency chirps down from approximately 35 kHz to 13 kHz, which corresponds approximately to the lower end of frequencies that are typically observed in EDA H-modes in ASDEX Upgrade.
△ Less
Submitted 22 January, 2023;
originally announced January 2023.
-
MHD simulations of formation, sustainment and loss of Quiescent H-mode in the all-tungsten ASDEX Upgrade
Authors:
Lorenz Meier,
Matthias Hoelzl,
Andres Cathey,
Guido Huijsmans,
Eleonora Viezzer,
Mike Dunne,
Jan van Dijk,
Diego José Cruz Zabala,
Karl Lackner,
Sibylle Günter,
ASDEX Upgrade Team,
EUROfusion MST1 Team,
JOREK Team
Abstract:
Periodic edge localized modes (ELMs) are the non-linear consequences of pressure-gradient-driven ballooning modes and current-driven peeling modes becoming unstable in the pedestal region of high confinement fusion plasmas. In future tokamaks like ITER, large ELMs are foreseen to severely affect the lifetime of wall components as they transiently deposit large amounts of heat onto a narrow region…
▽ More
Periodic edge localized modes (ELMs) are the non-linear consequences of pressure-gradient-driven ballooning modes and current-driven peeling modes becoming unstable in the pedestal region of high confinement fusion plasmas. In future tokamaks like ITER, large ELMs are foreseen to severely affect the lifetime of wall components as they transiently deposit large amounts of heat onto a narrow region at the divertor targets. Several strategies exist for avoidance, suppression, or mitigation of these instabilities, such as the naturally ELM-free quiescent H-mode (QH-mode). In the present article, an ASDEX Upgrade equilibrium that features a QH-mode is investigated through non-linear extended MHD simulations covering the dynamics over tens of milliseconds. The equilibrium is close to the ideal peeling limit and non-linearly develops saturated modes at the edge of the plasma. A dominant toroidal mode number of $n=1$ is found, for which the characteristic features of the edge harmonic oscillation are recovered. The saturated modes contribute to heat and particle transport preventing pedestal build-up to the ELM triggering threshold. The non-linear dynamics of the mode, in particular its interaction with the evolution of the edge safety factor is studied, which suggest a possible new saturation mechanism for the QH-mode. The simulations show good qualitative and quantitative agreement to experiments in AUG. In particular, the processes leading to the termination of QH-mode above a density threshold is studied, which results in the transition into an ELM regime. In the vicinity of this threshold, limit cycle oscillations are observed.
△ Less
Submitted 12 July, 2023; v1 submitted 18 January, 2023;
originally announced January 2023.
-
Development and application of a hybrid MHD-kinetic model in JOREK
Authors:
T. J. Bogaarts,
M. Hoelzl,
G. T. A. Huijsmans,
X. Wang
Abstract:
Energetic particle (EP) driven instabilities will be of strongly increased relevance in future burning plasmas as the EP pressure will be very large compared to the thermal plasma. Understanding the interaction of EPs and bulk plasma is crucial for developing next-generation fusion devices. In this work, the JOREK MHD code is extended to allow for the simulation of EP instabilities at high EP pres…
▽ More
Energetic particle (EP) driven instabilities will be of strongly increased relevance in future burning plasmas as the EP pressure will be very large compared to the thermal plasma. Understanding the interaction of EPs and bulk plasma is crucial for developing next-generation fusion devices. In this work, the JOREK MHD code is extended to allow for the simulation of EP instabilities at high EP pressures using realistic plasma and EP parameter in a full-f formulation with anisotropic pressure coupling to the bluid background. The code is first benchmarked linearly for the ITPA-TAE as well as the experiment based AUG-NLED cases, obtaining good agreement to other codes. Then, it is applied to a high energetic particle pressure discharge in the ASDEX Upgrade tokamak using a realistic non-Maxwellian distribution of EPs, reproducing aspects of the experimentally observed instabilities. Non-linear applications are possible based on the implentation, but will require dedicated verification and validation left for future work.
△ Less
Submitted 10 November, 2022; v1 submitted 3 August, 2022;
originally announced August 2022.
-
MHD simulations of small ELMs at low triangularity in ASDEX Upgrade
Authors:
A. Cathey,
M. Hoelzl,
G. Harrer,
M. G. Dunne,
G. T. A. Huijsmans,
K. Lackner,
S. J. P. Pamela,
E. Wolfrum,
S. Günter,
the JOREK team,
the ASDEX Upgrade Team,
the EUROfusion MST1 Team
Abstract:
The development of small- and no-ELM regimes for ITER is a high priority topic due to the risks associated to type-I ELMs. By considering non-linear extended MHD simulations of the ASDEX Upgrade tokamak with the JOREK code, we probe a regime that avoids type-I ELMs completely provided that the separatrix density is high enough. The dynamics of the pedestal in this regime are observed to be qualita…
▽ More
The development of small- and no-ELM regimes for ITER is a high priority topic due to the risks associated to type-I ELMs. By considering non-linear extended MHD simulations of the ASDEX Upgrade tokamak with the JOREK code, we probe a regime that avoids type-I ELMs completely provided that the separatrix density is high enough. The dynamics of the pedestal in this regime are observed to be qualitatively similar to the so-called quasi-continuous exhaust (QCE) regime in several ways. Repetitive type-I ELMs are substituted by roughly constant levels of outwards transport caused by peeling-ballooning modes (with dominant ballooning characteristics) which are localised in the last 5\% of the confined region (in normalised poloidal flux). The simulated low triangularity plasma transitions to a type-I ELMy H-mode if the separatrix density is sufficiently reduced or if the input heating power is sufficiently increased. The stabilising factors that play a role in the suppression of the small ELMs are also investigated by analysing the simulations, and the importance of including diamagnetic effects in the simulations is highlighted. By considering a scan in the pedestal resistivity and by measuring the poloidal velocity of the modes (and comparing to theoretical estimates for ideal and resistive modes), we identify the underlying instabilities as resistive peeling-ballooning modes. Decreasing the resistivity below experimentally-relevant conditions (i.e., going towards ideal MHD), the peeling-ballooning modes that constrain the pedestal below the type-I ELM stability boundary display sharply decreasing growth rates.
△ Less
Submitted 15 October, 2021;
originally announced October 2021.
-
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.
-
New plasma regimes with small ELMs and high confinement at the Joint European Torus
Authors:
J. Garcia,
E. de la Luna,
M. Sertoli,
F. Casson,
S. Mazzi,
Z. Stancar,
G. Szepesi,
D. Frigione,
L. Garzotti,
F. Rimini,
D. van Eester,
P. Lomas,
C. Sozzi,
R. N. Aiba,
R. Coelho,
L. Frasinetti,
G. Huijsmans,
F. Liu
Abstract:
New plasma regimes with high confinement, low core impurity accumulation and small Edge localized mode (ELMs) perturbations have been obtained close to ITER conditions in magnetically confined plasmas from the Joint European torus (JET) tokamak. Such regimes are achieved by means of optimized particle fuelling conditions which trigger a self-organize state with a strong increase in rotation and io…
▽ More
New plasma regimes with high confinement, low core impurity accumulation and small Edge localized mode (ELMs) perturbations have been obtained close to ITER conditions in magnetically confined plasmas from the Joint European torus (JET) tokamak. Such regimes are achieved by means of optimized particle fuelling conditions which trigger a self-organize state with a strong increase in rotation and ion temperature and a decrease of the edge density. An interplay between core and edge plasma regions leads to reduced turbulence levels and outward impurity convection. These results pave the way to an attractive alternative to the standard plasmas considered for fusion energy generation in a tokamak with metallic wall environment such as the ones expected in ITER
△ Less
Submitted 8 March, 2021; v1 submitted 3 March, 2021;
originally announced March 2021.
-
Comparing spontaneous and pellet-triggered ELMs via non-linear extended MHD simulations
Authors:
A. Cathey,
M. Hoelzl,
S. Futatani,
P. T. Lang,
K. Lackner,
G. T. A. Huijsmans,
S. J. P. Pamela,
S. Günter,
the JOREK team,
the ASDEX Upgrade Team,
the EUROfusion MST1 Team
Abstract:
Injecting frozen deuterium pellets into an ELMy H-mode plasma is a well established scheme for triggering edge localized modes (ELMs) before they naturally occur. Based on an ASDEX Upgrade H-mode plasma, this article presents a comparison of extended MHD simulations of spontaneous type-I ELMs and pellet-triggered ELMs allowing to study their non-linear dynamics in detail. In particular, pellet-tri…
▽ More
Injecting frozen deuterium pellets into an ELMy H-mode plasma is a well established scheme for triggering edge localized modes (ELMs) before they naturally occur. Based on an ASDEX Upgrade H-mode plasma, this article presents a comparison of extended MHD simulations of spontaneous type-I ELMs and pellet-triggered ELMs allowing to study their non-linear dynamics in detail. In particular, pellet-triggered ELMs are simulated by injecting deuterium pellets into different time points during the pedestal build-up described in [A. Cathey et al. Nuclear Fusion 60, 124007 (2020)]. Realistic ExB and diamagnetic background plasma flows as well as the time dependent bootstrap current evolution are included during the build-up to capture the balance between stabilising and destabilising terms for the edge instabilities accurately. Dependencies on the pellet size and injection times are studied. The spatio-temporal structures of the modes and the resulting divertor heat fluxes are compared in detail between spontaneous and triggered ELMs. We observe that the premature excitation of ELMs by means of pellet injection is caused by a helical perturbation described by a toroidal mode number of n = 1. In accordance with experimental observations, the pellet-triggered ELMs show reduced thermal energy losses and narrower divertor wetted area with respect to spontaneous ELMs. The peak divertor energy fluency is seen to decrease when ELMs are triggered by pellets injected earlier during the pedestal build-up.
△ Less
Submitted 11 February, 2021;
originally announced February 2021.
-
Enhanced Preconditioner for JOREK MHD Solver
Authors:
I Holod,
M Hoelzl,
P S Verma,
GTA Huijsmans,
R Nies,
JOREK Team
Abstract:
The JOREK extended magneto-hydrodynamic (MHD) code is a widely used simulation code for studying the non-linear dynamics of large-scale instabilities in divertor tokamak plasmas. Due to the large scale-separation intrinsic to these phenomena both in space and time, the computational costs for simulations in realistic geometry and with realistic parameters can be very high, motivating the investmen…
▽ More
The JOREK extended magneto-hydrodynamic (MHD) code is a widely used simulation code for studying the non-linear dynamics of large-scale instabilities in divertor tokamak plasmas. Due to the large scale-separation intrinsic to these phenomena both in space and time, the computational costs for simulations in realistic geometry and with realistic parameters can be very high, motivating the investment of considerable effort for optimization. In this article, a set of developments regarding the JOREK solver and preconditioner is described, which lead to overall significant benefits for large production simulations. This comprises in particular enhanced convergence in highly non-linear scenarios and a general reduction of memory consumption and computational costs. The developments include faster construction of preconditioner matrices, a domain decomposition of preconditioning matrices for solver libraries that can handle distributed matrices, interfaces for additional solver libraries, an option to use matrix compression methods, and the implementation of a complex solver interface for the preconditioner. The most significant development presented consists in a generalization of the physics based preconditioner to "mode groups", which allows to account for the dominant interactions between toroidal Fourier modes in highly non-linear simulations. At the cost of a moderate increase of memory consumption, the technique can strongly enhance convergence in suitable cases allowing to use significantly larger time steps. For all developments, benchmarks based on typical simulation cases demonstrate the resulting improvements.
△ Less
Submitted 31 March, 2021; v1 submitted 20 January, 2021;
originally announced January 2021.
-
Simulations of COMPASS Vertical Displacement Events with a self-consistent model for halo currents including neutrals and sheath boundary conditions
Authors:
F. J. Artola,
A. Loarte,
E. Matveeva,
J. Havlicek,
T. Markovic,
J. Adamek,
J. Cavalier,
L. Kripner,
G. T. A. Huijsmans,
M. Lehnen,
M. Hoelzl,
R. Panek,
COMPASS team,
JOREK team
Abstract:
The understanding of the halo current properties during disruptions is key to design and operate large scale tokamaks in view of the large thermal and electromagnetic loads that they entail. For the first time, we present a fully self-consistent model for halo current simulations including neutral particles and sheath boundary conditions. The model is used to simulate Vertical Displacement Events…
▽ More
The understanding of the halo current properties during disruptions is key to design and operate large scale tokamaks in view of the large thermal and electromagnetic loads that they entail. For the first time, we present a fully self-consistent model for halo current simulations including neutral particles and sheath boundary conditions. The model is used to simulate Vertical Displacement Events (VDEs) occurring in the COMPASS tokamak. Recent COMPASS experiments have shown that the parallel halo current density at the plasma-wall interface is limited by the ion saturation current during VDE-induced disruptions. We show that usual MHD boundary conditions can lead to the violation of this physical limit and we implement this current density limitation through a boundary condition for the electrostatic potential. Sheath boundary conditions for the density, the heat flux, the parallel velocity and a realistic parameter choice (e.g. Spitzer $η$ and Spitzer-Härm $χ_\parallel$ values) extend present VDE simulations beyond the state of the art. Experimental measurements of the current density, temperature and heat flux profiles at the COMPASS divertor are compared with the results obtained from axisymmetric simulations. Since the ion saturation current density ($J_{sat}$) is shown to be essential to determine the halo current profile, parametric scans are performed to study its dependence on different quantities such as the plasma resistivity and the particle and heat diffusion coefficients. In this respect, the plasma resistivity in the halo region broadens significantly the $J_{sat}$ profile, increasing the halo width at a similar total halo current.
△ Less
Submitted 5 January, 2021;
originally announced January 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.
-
Transition from no-ELM response to pellet ELM triggering during pedestal build-up -- insights from extended MHD simulations
Authors:
S Futatani,
A Cathey,
M Hoelzl,
P T Lang,
G T A Huijsmans,
M. Dunne,
JOREK Team,
ASDEX Upgrade Team,
EUROfusion MST1 Team
Abstract:
Pellet ELM triggering is a well established scheme for decreasing the time between two successive ELM crashes below its natural value. Reliable ELM pacing has been demonstrated experimentally in several devices increasing the ELM frequency considerably. However, it was also shown that the frequency cannot be increased arbitrarily due to a so-called lag-time. During this time after a preceding natu…
▽ More
Pellet ELM triggering is a well established scheme for decreasing the time between two successive ELM crashes below its natural value. Reliable ELM pacing has been demonstrated experimentally in several devices increasing the ELM frequency considerably. However, it was also shown that the frequency cannot be increased arbitrarily due to a so-called lag-time. During this time after a preceding natural or triggered ELM crash, neither a natural ELM crash occurs nor the triggering of an ELM crash by pellet injection is possible. For this article, pellet ELM triggering simulations are advanced beyond previous studies in two ways. Firstly, realistic ExB and diamagnetic background flows are included. And secondly, the pellet is injected at different stages of the pedestal build-up. This allows to recover the lag-time for the first time in simulations and investigate it in detail. A series of non-linear extended MHD simulations is performed to investigate the plasma dynamics resulting from an injection at different time points during the pedestal build-up. The experimentally observed lag-time is qualitatively reproduced well. In particular, a sharp transition is observed between the regime where no ELMs can be triggered and the regime where pellet injection causes an ELM crash. Via variations of pellet parameters and injection time, the two regimes are studied and compared in detail revealing pronounced differences in the non-linear dynamics. The toroidal mode spectrum is significantly broader when an ELM crash is triggered enhancing the stochasticity and therefore also the losses of thermal energy along magnetic field lines. In the heat fluxes to the divertor targets, pronounced toroidal asymmetries are observed. In case of high injection velocities leading to deep penetration, also the excitation of core modes like the $2/1$ neoclassical tearing mode is observed.
△ Less
Submitted 24 December, 2020; v1 submitted 16 September, 2020;
originally announced September 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.
-
Non-linear extended MHD simulations of type-I edge localised mode cycles in ASDEX Upgrade and their underlying triggering mechanism
Authors:
Andres Cathey,
M. Hoelzl,
K. Lackner,
G. T. A. Huijsmans,
M. G. Dunne,
E. Wolfrum,
S. J. P. Pamela,
F. Orain,
S. Günter,
the JOREK team,
the ASDEX Upgrade Team,
the EUROfusion MST1 Team
Abstract:
A triggering mechanism responsible for the explosive onset of edge localised modes (ELMs) in fusion plasmas is identified by performing, for the first time, non-linear magnetohydrodynamic simulations of repetitive type-I ELMs. Briefly prior to the ELM crash, destabilising and stabilising terms are affected at different timescales by an increasingly ergodic magnetic field caused by non-linear inter…
▽ More
A triggering mechanism responsible for the explosive onset of edge localised modes (ELMs) in fusion plasmas is identified by performing, for the first time, non-linear magnetohydrodynamic simulations of repetitive type-I ELMs. Briefly prior to the ELM crash, destabilising and stabilising terms are affected at different timescales by an increasingly ergodic magnetic field caused by non-linear interactions between the axisymmetric background plasma and growing non-axisymmetric perturbations. The separation of timescales prompts the explosive, i.e. faster than exponential, growth of an ELM crash which lasts ${\sim}$ 500 $μ$s. The duration and size of the simulated ELM crashes compare qualitatively well with type-I ELMs in ASDEX Upgrade. As expected for type-I ELMs, a direct proportionality between the heating power in the simulations and the ELM repetition frequency is obtained. The simulations presented here are a major step forward towards predictive modelling of ELMs and of the assessment of mitigation techniques in ITER and other future tokamaks.
△ Less
Submitted 26 October, 2020; v1 submitted 20 July, 2020;
originally announced July 2020.
-
The initial step towards JOREK integration in IMAS
Authors:
Dejan Penko,
Leon Kos,
Guido Huijsmans,
Simon D. Pinches
Abstract:
JOREK is being adapted to work with the Integrated Modelling & Analysis Suite (IMAS) which is being actively developed and used by the ITER Organization, the EUROfusion community and other ITER Members. The list of codes adapted to use the IMAS Data Model is gradually increasing with examples including SOLPS-ITER and JINTRAC. The main goal of the integration of JOREK with IMAS is to enable interac…
▽ More
JOREK is being adapted to work with the Integrated Modelling & Analysis Suite (IMAS) which is being actively developed and used by the ITER Organization, the EUROfusion community and other ITER Members. The list of codes adapted to use the IMAS Data Model is gradually increasing with examples including SOLPS-ITER and JINTRAC. The main goal of the integration of JOREK with IMAS is to enable interaction with the plasma scenarios stored in the IMAS databases in the form of Interface Data Structures (IDSs): input conditions can be read from the databases and nonlinear plasma states determined by JOREK stored. IDSs provide a uniform way of representing data within the IMAS framework and allow to transfer data between codes and to storage within larger integrated modelling workflows. In order to integrate JOREK within IMAS it is therefore necessary that transformation tools are developed to facilitate the reading and writing of the relevant IDSs, including the MHD IDS, with its underlying Generalized Grid Description (GGD). For this purpose, utilities have been developed that extract JOREK simulation plasma state, namely the grid geometry and computed physical quantities for each time slice, and then transform them to the appropriate output IDSs. In this article, these initial steps towards full JOREK integration into IMAS is presented.
△ Less
Submitted 16 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.
-
Simulating the nonlinear interaction of relativistic electrons and tokamak plasma instabilities: Implementation and validation of a fluid model
Authors:
V. Bandaru,
M. Hoelzl,
F. J. Artola,
G. Papp,
G. T. A. Huijsmans
Abstract:
For the simulation of disruptions in tokamak fusion plasmas, a fluid model describing the evolution of relativistic runaway electrons and their interaction with the background plasma is presented. The overall aim of the model is to self-consistently describe the nonlinear coupled evolution of runaway electrons (REs) and plasma instabilities during disruptions. In this model, the runaway electrons…
▽ More
For the simulation of disruptions in tokamak fusion plasmas, a fluid model describing the evolution of relativistic runaway electrons and their interaction with the background plasma is presented. The overall aim of the model is to self-consistently describe the nonlinear coupled evolution of runaway electrons (REs) and plasma instabilities during disruptions. In this model, the runaway electrons are considered as a separate fluid species in which the initial seed is generated through the Dreicer source, which eventually grows by the avalanche mechanism (further relevant source mechanisms can easily be added). Advection of the runaway electrons is considered primarily along field lines, but also taking into account the ExB drift. The model is implemented in the nonlinear magnetohydrodynamic (MHD) code JOREK based on Bezier finite elements, with current coupling to the thermal plasma. Benchmarking of the code with the one-dimensional runaway electron code GO is done using an artificial thermal quench on a circular plasma. As a first demonstration, the code is applied to the problem of an axisymmetric cold vertical displacement event in an ITER plasma, revealing significantly different dynamics between cases computed with and without runaway electrons. Though it is not yet feasible to achieve fully realistic runaway electron velocities close to the speed of light in complete simulations of slowly evolving plasma instabilities, the code is demonstrated to be suitable to study various kinds of MHD-RE interactions in MHD-active and disruption relevant plasmas.
△ Less
Submitted 28 June, 2019;
originally announced June 2019.
-
Numerical study of tearing mode seeding in tokamak X-point plasma
Authors:
Dmytro Meshcheriakov,
Matthias Hoelzl,
Valentin Igochine,
Sina Fietz,
Francois Orain,
Guido T. A. Huijsmans,
Marc Maraschek,
Mike Dunne,
Rachael McDermott,
Hartmut Zohm,
Karl Lackner,
Sibylle Guenter,
ASDEX Upgrade Team,
EUROfusion MST1 Team
Abstract:
A detailed understanding of island seeding is crucial to avoid (N)TMs and their negative consequences like confinement degradation and disruptions. In the present work, we investigate the growth of 2/1 islands in response to magnetic perturbations. Although we use externally applied perturbations produced by resonant magnetic perturbation (RMP) coils for this study, results are directly transferab…
▽ More
A detailed understanding of island seeding is crucial to avoid (N)TMs and their negative consequences like confinement degradation and disruptions. In the present work, we investigate the growth of 2/1 islands in response to magnetic perturbations. Although we use externally applied perturbations produced by resonant magnetic perturbation (RMP) coils for this study, results are directly transferable to island seeding by other MHD instabilities creating a resonant magnetic field component at the rational surface. Experimental results for 2/1 island penetration from ASDEX Upgrade are presented extending previous studies. Simulations are based on an ASDEX Upgrade L-mode discharge with low collisionality and active RMP coils. Our numerical studies are performed with the 3D, two fluid, non-linear MHD code JOREK. All three phases of mode seeding observed in the experiment are also seen in the simulations: first a weak response phase characterized by large perpendicular electron flow velocities followed by a fast growth of the magnetic island size accompanied by a reduction of the perpendicular electron velocity, and finally the saturation to a fully formed island state with perpendicular electron velocity close to zero. Thresholds for mode penetration are observed in the plasma rotation as well as in the RMP coil current. A hysteresis of the island size and electron perpendicular velocity is observed between the ramping up and down of the RMP amplitude consistent with an analytically predicted bifurcation. The transition from dominant kink/bending to tearing parity during the penetration is investigated.
△ Less
Submitted 16 April, 2019;
originally announced April 2019.
-
Non-linear modeling of the threshold between ELM mitigation and ELM suppression by Resonant Magnetic Perturbations in ASDEX Upgrade
Authors:
François Orain,
M. Hoelzl,
F. Mink,
M. Willensdorfer,
M. Bécoulet,
M. Dunne,
S. Günter,
G. T. A. Huijsmans,
K. Lackner,
S. Pamela,
W. Suttrop,
E. Viezzer
Abstract:
The interaction between Edge Localized Modes (ELMs) and Resonant Magnetic Perturbations (RMPs) is modeled with the magnetohydrodynamic code JOREK using experimental parameters from ASDEX Upgrade discharges. The ELM mitigation or suppression is optimal when the amplification of both tearing and peeling-kink responses result in a better RMP penetration. The ELM mitigation or suppression is not only…
▽ More
The interaction between Edge Localized Modes (ELMs) and Resonant Magnetic Perturbations (RMPs) is modeled with the magnetohydrodynamic code JOREK using experimental parameters from ASDEX Upgrade discharges. The ELM mitigation or suppression is optimal when the amplification of both tearing and peeling-kink responses result in a better RMP penetration. The ELM mitigation or suppression is not only due to the reduction of the pressure gradient, but predominantly arises from the toroidal coupling between the ELMs and the RMP-induced mode at the plasma edge, forcing the edge modes to saturate at a low level. The bifurcation from ELM mitigation to ELM suppression is observed when the RMP amplitude is increased. ELM mitigation is characterized by rotating modes at the edge, while the mode locking to RMPs is induced by the resonant braking of the electron perpendicular flow in the ELM suppression regime.
△ Less
Submitted 3 June, 2019; v1 submitted 1 February, 2019;
originally announced February 2019.
-
Performance analysis and optimization of the JOREK code for many-core CPUs
Authors:
T. B. Fehér,
M. Hölzl,
G. Latu,
G. T. A. Huijsmans
Abstract:
This report investigates the performance of the JOREK code on the Intel Knights Landing and Skylake processor architectures. The OpenMP scaling of the matrix construction part of the code was analyzed and improved synchronization methods were implemented. A new switch was implemented to control the number of threads used for the linear equation solver independently from other parts of the code. Th…
▽ More
This report investigates the performance of the JOREK code on the Intel Knights Landing and Skylake processor architectures. The OpenMP scaling of the matrix construction part of the code was analyzed and improved synchronization methods were implemented. A new switch was implemented to control the number of threads used for the linear equation solver independently from other parts of the code. The matrix construction subroutine was vectorized, and the data locality was also improved. These steps led to a factor of two speedup for the matrix construction.
△ Less
Submitted 10 October, 2018;
originally announced October 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.
-
Early evolution of electron cyclotron driven current during suppression of tearing modes in a circular tokamak
Authors:
J. Pratt,
G. T. A. Huijsmans,
E. Westerhof
Abstract:
When electron cyclotron (EC) driven current is first applied to the inside of a magnetic island, the current spreads throughout the island and after a short period achieves a steady level. Using a two equation fluid model for the EC current that allows us to examine this early evolution in detail, we analyze high-resolution simulations of a 2/1 classical tearing mode in a low-beta large aspect-rat…
▽ More
When electron cyclotron (EC) driven current is first applied to the inside of a magnetic island, the current spreads throughout the island and after a short period achieves a steady level. Using a two equation fluid model for the EC current that allows us to examine this early evolution in detail, we analyze high-resolution simulations of a 2/1 classical tearing mode in a low-beta large aspect-ratio circular tokamak. These simulations use a nonlinear 3D reduced-MHD fluid model and the JOREK code. During the initial period where the EC driven current grows and spreads throughout the magnetic island, it is not a function of the magnetic flux. However, once it has reached a steady-state, it should be a flux function. We demonstrate numerically that if sufficiently resolved toroidally, the steady-state EC driven current becomes approximately a flux function. We discuss the physics of this early period of EC evolution and its impact on the size of the magnetic island.
△ Less
Submitted 28 September, 2016;
originally announced September 2016.
-
Non-linear modeling of the plasma response to RMPs in ASDEX Upgrade
Authors:
F. Orain,
M. Hoelzl,
E. Viezzer,
M. Dunne,
M. Becoulet,
P. Cahyna,
G. T. A. Huijsmans,
J. Morales,
M. Willensdorfer,
W. Suttrop,
A. Kirk,
S. Pamela,
E. Strumberger,
S. Guenter,
A. Lessig,
the ASDEX Upgrade Team,
the EUROfusion MST1 Team
Abstract:
The plasma response to Resonant Magnetic Perturbations (RMPs) in ASDEX Upgrade is modeled with the non-linear resistive MHD code JOREK, using input profiles that match those of the experiments as closely as possible. The RMP configuration for which Edge Localized Modes are best mitigated in experiments is related to the largest edge kink response observed near the X-point in modeling. On the edge…
▽ More
The plasma response to Resonant Magnetic Perturbations (RMPs) in ASDEX Upgrade is modeled with the non-linear resistive MHD code JOREK, using input profiles that match those of the experiments as closely as possible. The RMP configuration for which Edge Localized Modes are best mitigated in experiments is related to the largest edge kink response observed near the X-point in modeling. On the edge resonant surfaces $q = m/n$, the coupling between the m + 2 kink component and the m resonant component is found to induce the amplification of the resonant magnetic perturbation. The ergodicity and the 3D-displacement near the X-point induced by the resonant amplification can only partly explain the density pumpout observed in experiments.
△ Less
Submitted 24 February, 2016;
originally announced February 2016.
-
Modelling of spatial structure of divertor footprints caused by edge-localized modes mitigated by magnetic perturbations
Authors:
Pavel Cahyna,
Marina Becoulet,
Guido T. A. Huijsmans,
Francois Orain,
Jorge Morales,
Andrew Kirk,
Andrew J. Thornton,
Stanislas Pamela,
Radomir Panek,
Matthias Hoelzl
Abstract:
Resonant magnetic perturbations (RMPs) can mitigate the edge-localized modes (ELMs), i.e. cause a change of the ELM character towards smaller energy loss and higher frequency. During mitigation a change of the spatial structure of ELM loads on divertor was observed on DIII-D and MAST: the power is deposited predominantly in the footprint structures formed by the magnetic perturbation. In the prese…
▽ More
Resonant magnetic perturbations (RMPs) can mitigate the edge-localized modes (ELMs), i.e. cause a change of the ELM character towards smaller energy loss and higher frequency. During mitigation a change of the spatial structure of ELM loads on divertor was observed on DIII-D and MAST: the power is deposited predominantly in the footprint structures formed by the magnetic perturbation. In the present contribution we develop a theory explaining this effect, based on the idea that part of the ELM loss is caused by parallel transport in the homoclinic tangle formed by the magnetic perturbation of the ELM. The modified tangle resulting from the combination of the ELM perturbation and the applied RMP has the expected property of bringing open field lines in the same areas as the tangle from the RMP alone. We show that this explanation is consistent with features of the mitigated ELMs on MAST.
We in addition validated our theory by an analysis of simulations of mitigated ELMs using the code JOREK. We produced detailed laminar plots of field lines on the divertor in the JOREK runs with an ELM, an applied RMP, and an ELM mitigated by the presence of the RMP. The results for an ELM clearly show a high-n rotating footprint structure appearing during the nonlinear stage of the ELM, which is not present in the early stage of the ELM. The results for a n=2 RMP from the ELM control coils show the expected n=2 footprint structure. The results for the mitigated ELM show a similar structure, modulated by a higher n perturbation of the ELM, consistent with our theory.
△ Less
Submitted 12 January, 2016;
originally announced January 2016.
-
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.
-
Current understanding of the processes underlying the triggering of and energy loss associated with type I ELMs
Authors:
A. Kirk,
D. Dunai,
M. Dunne,
G. Huijsmans,
S. Pamela,
M. Becoulet,
J. R. Harrison,
J. Hillesheim,
C. Roach,
S. Saarelma
Abstract:
The type I ELMy H-mode is the baseline operating scenario for ITER. While it is known that the type I ELM ultimately results from the peeling-ballooning instability, there is growing experimental evidence that a mode grows up before the ELM crash that may modify the edge plasma, which then leads to the ELM event due to the peeling-ballooning mode. The triggered mode results in the release of a lar…
▽ More
The type I ELMy H-mode is the baseline operating scenario for ITER. While it is known that the type I ELM ultimately results from the peeling-ballooning instability, there is growing experimental evidence that a mode grows up before the ELM crash that may modify the edge plasma, which then leads to the ELM event due to the peeling-ballooning mode. The triggered mode results in the release of a large number of particles and energy from the core plasma but the precise mechanism by which these losses occur is still not fully understood and hence makes predictions for future devices uncertain. Our current understanding of the processes that trigger type I ELMs and the size of the resultant energy loss are reviewed and compared to experimental data and ideas for further development are discussed.
△ Less
Submitted 16 December, 2013;
originally announced December 2013.
-
Understanding the effect resonant magnetic perturbations have on ELMs
Authors:
A. Kirk,
I. T. Chapman,
T. E. Evans,
C. Ham,
J. R. Harrison,
G. Huijsmans,
Y. Liang,
Y. Q. Liu,
A. Loarte,
W. Suttrop,
A. J. Thornton
Abstract:
All current estimations of the energy released by type I ELMs indicate that, in order to ensure an adequate lifetime of the divertor targets on ITER, a mechanism is required to decrease the amount of energy released by an ELM, or to eliminate ELMs altogether. One such amelioration mechanism relies on perturbing the magnetic field in the edge plasma region, either leading to more frequent, smaller…
▽ More
All current estimations of the energy released by type I ELMs indicate that, in order to ensure an adequate lifetime of the divertor targets on ITER, a mechanism is required to decrease the amount of energy released by an ELM, or to eliminate ELMs altogether. One such amelioration mechanism relies on perturbing the magnetic field in the edge plasma region, either leading to more frequent, smaller ELMs (ELM mitigation) or ELM suppression. This technique of Resonant Magnetic Perturbations (RMPs) has been employed to suppress type I ELMs at high collisionality/density on DIII-D, ASDEX Upgrade, KSTAR and JET and at low collisionality on DIII-D. At ITER-like collisionality the RMPs enhance the transport of particles or energy and keep the edge pressure gradient below the 2D linear ideal MHD critical value that would trigger an ELM, whereas at high collisionality/density the type I ELMs are replaced by small type II ELMs. Although ELM suppression only occurs within limitied operational ranges, ELM mitigation is much more easily achieved. The exact parameters that determine the onset of ELM suppression are unknown but in all cases the magnetic perturbations produce 3D distortions to the plasma and enhanced particle transport. The incorporation of these 3D effects in codes will be essential in order to make quantitative predictions for future devices.
△ Less
Submitted 28 June, 2013;
originally announced June 2013.