-
Valid Error Bars for Neural Weather Models using Conformal Prediction
Authors:
Vignesh Gopakumar,
Joel Oskarrson,
Ander Gray,
Lorenzo Zanisi,
Stanislas Pamela,
Daniel Giles,
Matt Kusner,
Marc Deisenroth
Abstract:
Neural weather models have shown immense potential as inexpensive and accurate alternatives to physics-based models. However, most models trained to perform weather forecasting do not quantify the uncertainty associated with their forecasts. This limits the trust in the model and the usefulness of the forecasts. In this work we construct and formalise a conformal prediction framework as a post-pro…
▽ More
Neural weather models have shown immense potential as inexpensive and accurate alternatives to physics-based models. However, most models trained to perform weather forecasting do not quantify the uncertainty associated with their forecasts. This limits the trust in the model and the usefulness of the forecasts. In this work we construct and formalise a conformal prediction framework as a post-processing method for estimating this uncertainty. The method is model-agnostic and gives calibrated error bounds for all variables, lead times and spatial locations. No modifications are required to the model and the computational cost is negligible compared to model training. We demonstrate the usefulness of the conformal prediction framework on a limited area neural weather model for the Nordic region. We further explore the advantages of the framework for deterministic and probabilistic models.
△ Less
Submitted 20 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.
-
Emulation Techniques for Scenario and Classical Control Design of Tokamak Plasmas
Authors:
A. Agnello,
N. C. Amorisco,
A. Keats,
G. K. Holt,
J. Buchanan,
S. Pamela,
C. Vincent,
G. McArdle
Abstract:
The optimisation of scenarios and design of real-time-control in tokamaks, especially for machines still in design phase, requires a comprehensive exploration of solutions to the Grad-Shafranov (GS) equation over a high-dimensional space of plasma and coil parameters. Emulators can bypass the numerical issues in the GS equation, if a large enough library of equilibria is available. We train an ens…
▽ More
The optimisation of scenarios and design of real-time-control in tokamaks, especially for machines still in design phase, requires a comprehensive exploration of solutions to the Grad-Shafranov (GS) equation over a high-dimensional space of plasma and coil parameters. Emulators can bypass the numerical issues in the GS equation, if a large enough library of equilibria is available. We train an ensemble of neural networks to emulate the typical shape-control targets (separatrix at midplane, X-points, divertor strike point, flux expansion, poloidal beta) as a function of plasma parameters and active coil currents for the range of plasma configurations relevant to spherical tokamaks with a super-X divertor, with percent-level accuracy. This allows a quick calculation of the classical-control shape matrices, potentially allowing real-time calculation at any point in a shot with sub-ms latency. We devise a hyperparameter sampler to select the optimal network architectures and quantify uncertainties on the model predictions. To generate the relevant training set, we devise a Markov-Chain Monte Carlo algorithm to produce large libraries of forward Grad-Shafranov solutions without the need for user intervention. The algorithm promotes equilibria with desirable properties, while avoiding parameter combinations resulting in problematic profiles or numerical issues in the integration of the GS equation.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
Data efficiency and long term prediction capabilities for neural operator surrogate models of core and edge plasma codes
Authors:
N. Carey,
L. Zanisi,
S. Pamela,
V. Gopakumar,
J. Omotani,
J. Buchanan,
J. Brandstetter
Abstract:
Simulation-based plasma scenario development, optimization and control are crucial elements towards the successful deployment of next-generation experimental tokamaks and Fusion power plants. Current simulation codes require extremely intensive use of HPC resources that make them unsuitable for iterative or real time applications. Neural network based surrogate models of expensive simulators have…
▽ More
Simulation-based plasma scenario development, optimization and control are crucial elements towards the successful deployment of next-generation experimental tokamaks and Fusion power plants. Current simulation codes require extremely intensive use of HPC resources that make them unsuitable for iterative or real time applications. Neural network based surrogate models of expensive simulators have been proposed to speed up such costly workflows. Current efforts in this direction in the Fusion community are mostly limited to point estimates of quantities of interest or simple 1D PDE models, with a few notable exceptions. While the AI literature on methods for neural PDE surrogate models is rich, performance benchmarks for Fusion-relevant 2D fields has so far remained flimited. In this work neural PDE surrogates are trained for the JOREK MHD code and the STORM scrape-off layer code using the PDEArena library (https://github.com/microsoft/pdearena). The performance of these surrogate models is investigated as a function of training set size as well as for long-term predictions. The performance of surrogate models that are trained on either one variable or multiple variables at once is also considered. It is found that surrogates that are trained on more data perform best for both long- and short-term predictions. Additionally, surrogate models trained on multiple variables achieve higher accuracy and more stable performance. Downsampling the training set in time may provide stability in the long term at the expense of the short term predictive capability, but visual inspection of the resulting fields suggests that multiple metrics should be used to evaluate performance.
△ Less
Submitted 13 February, 2024;
originally announced February 2024.
-
Plasma Surrogate Modelling using Fourier Neural Operators
Authors:
Vignesh Gopakumar,
Stanislas Pamela,
Lorenzo Zanisi,
Zongyi Li,
Ander Gray,
Daniel Brennand,
Nitesh Bhatia,
Gregory Stathopoulos,
Matt Kusner,
Marc Peter Deisenroth,
Anima Anandkumar,
JOREK Team,
MAST Team
Abstract:
Predicting plasma evolution within a Tokamak reactor is crucial to realizing the goal of sustainable fusion. Capabilities in forecasting the spatio-temporal evolution of plasma rapidly and accurately allow us to quickly iterate over design and control strategies on current Tokamak devices and future reactors. Modelling plasma evolution using numerical solvers is often expensive, consuming many hou…
▽ More
Predicting plasma evolution within a Tokamak reactor is crucial to realizing the goal of sustainable fusion. Capabilities in forecasting the spatio-temporal evolution of plasma rapidly and accurately allow us to quickly iterate over design and control strategies on current Tokamak devices and future reactors. Modelling plasma evolution using numerical solvers is often expensive, consuming many hours on supercomputers, and hence, we need alternative inexpensive surrogate models. We demonstrate accurate predictions of plasma evolution both in simulation and experimental domains using deep learning-based surrogate modelling tools, viz., Fourier Neural Operators (FNO). We show that FNO has a speedup of six orders of magnitude over traditional solvers in predicting the plasma dynamics simulated from magnetohydrodynamic models, while maintaining a high accuracy (MSE in the normalised domain $\approx$ $10^{-5}$). Our modified version of the FNO is capable of solving multi-variable Partial Differential Equations (PDE), and can capture the dependence among the different variables in a single model. FNOs can also predict plasma evolution on real-world experimental data observed by the cameras positioned within the MAST Tokamak, i.e., cameras looking across the central solenoid and the divertor in the Tokamak. We show that FNOs are able to accurately forecast the evolution of plasma and have the potential to be deployed for real-time monitoring. We also illustrate their capability in forecasting the plasma shape, the locations of interactions of the plasma with the central solenoid and the divertor for the full (available) duration of the plasma shot within MAST. The FNO offers a viable alternative for surrogate modelling as it is quick to train and infer, and requires fewer data points, while being able to do zero-shot super-resolution and getting high-fidelity solutions.
△ Less
Submitted 18 June, 2024; v1 submitted 10 November, 2023;
originally announced November 2023.
-
Efficient training sets for surrogate models of tokamak turbulence with Active Deep Ensembles
Authors:
L. Zanisi,
A. Ho,
T. Madula,
J. Barr,
J. Citrin,
S. Pamela,
J. Buchanan,
F. Casson,
V. Gopakumar,
JET contributors
Abstract:
Model-based plasma scenario development lies at the heart of the design and operation of future fusion powerplants. Including turbulent transport in integrated models is essential for delivering a successful roadmap towards operation of ITER and the design of DEMO-class devices. Given the highly iterative nature of integrated models, fast machine-learning-based surrogates of turbulent transport ar…
▽ More
Model-based plasma scenario development lies at the heart of the design and operation of future fusion powerplants. Including turbulent transport in integrated models is essential for delivering a successful roadmap towards operation of ITER and the design of DEMO-class devices. Given the highly iterative nature of integrated models, fast machine-learning-based surrogates of turbulent transport are fundamental to fulfil the pressing need for faster simulations opening up pulse design, optimization, and flight simulator applications. A significant bottleneck is the generation of suitably large training datasets covering a large volume in parameter space, which can be prohibitively expensive to obtain for higher fidelity codes.
In this work, we propose ADEPT (Active Deep Ensembles for Plasma Turbulence), a physics-informed, two-stage Active Learning strategy to ease this challenge. Active Learning queries a given model by means of an acquisition function that identifies regions where additional data would improve the surrogate model. We provide a benchmark study using available data from the literature for the QuaLiKiz quasilinear transport model. We demonstrate quantitatively that the physics-informed nature of the proposed workflow reduces the need to perform simulations in stable regions of the parameter space, resulting in significantly improved data efficiency. We show an up to a factor of 20 reduction in training dataset size needed to achieve the same performance as random sampling. We then validate the surrogates on multichannel integrated modelling of ITG-dominated JET scenarios and demonstrate that they recover the performance of QuaLiKiz to better than 10\%. This matches the performance obtained in previous work, but with two orders of magnitude fewer training data points.
△ Less
Submitted 13 October, 2023;
originally announced October 2023.
-
Fourier Neural Operator for Plasma Modelling
Authors:
Vignesh Gopakumar,
Stanislas Pamela,
Lorenzo Zanisi,
Zongyi Li,
Anima Anandkumar,
MAST Team
Abstract:
Predicting plasma evolution within a Tokamak is crucial to building a sustainable fusion reactor. Whether in the simulation space or within the experimental domain, the capability to forecast the spatio-temporal evolution of plasma field variables rapidly and accurately could improve active control methods on current tokamak devices and future fusion reactors. In this work, we demonstrate the util…
▽ More
Predicting plasma evolution within a Tokamak is crucial to building a sustainable fusion reactor. Whether in the simulation space or within the experimental domain, the capability to forecast the spatio-temporal evolution of plasma field variables rapidly and accurately could improve active control methods on current tokamak devices and future fusion reactors. In this work, we demonstrate the utility of using Fourier Neural Operator (FNO) to model the plasma evolution in simulations and experiments. Our work shows that the FNO is capable of predicting magnetohydrodynamic models governing the plasma dynamics, 6 orders of magnitude faster than the traditional numerical solver, while maintaining considerable accuracy (NMSE $\sim 10^{-5})$. Our work also benchmarks the performance of the FNO against other standard surrogate models such as Conv-LSTM and U-Net and demonstrate that the FNO takes significantly less time to train, requires less parameters and outperforms other models. We extend the FNO approach to model the plasma evolution observed by the cameras positioned within the MAST spherical tokamak. We illustrate its capability in forecasting the formation of filaments within the plasma as well as the heat deposits. The FNO deployed to model the camera is capable of forecasting the full length of the plasma shot within half the time of the shot duration.
△ Less
Submitted 13 February, 2023;
originally announced February 2023.
-
Fourier-RNNs for Modelling Noisy Physics Data
Authors:
Vignesh Gopakumar,
Stanislas Pamela,
Lorenzo Zanisi
Abstract:
Classical sequential models employed in time-series prediction rely on learning the mappings from the past to the future instances by way of a hidden state. The Hidden states characterise the historical information and encode the required temporal dependencies. However, most existing sequential models operate within finite-dimensional Euclidean spaces which offer limited functionality when employe…
▽ More
Classical sequential models employed in time-series prediction rely on learning the mappings from the past to the future instances by way of a hidden state. The Hidden states characterise the historical information and encode the required temporal dependencies. However, most existing sequential models operate within finite-dimensional Euclidean spaces which offer limited functionality when employed in modelling physics relevant data. Alternatively recent work with neural operator learning within the Fourier space has shown efficient strategies for parameterising Partial Differential Equations (PDE). In this work, we propose a novel sequential model, built to handle Physics relevant data by way of amalgamating the conventional RNN architecture with that of the Fourier Neural Operators (FNO). The Fourier-RNN allows for learning the mappings from the input to the output as well as to the hidden state within the Fourier space associated with the temporal data. While the Fourier-RNN performs identical to the FNO when handling PDE data, it outperforms the FNO and the conventional RNN when deployed in modelling noisy, non-Markovian data.
△ Less
Submitted 13 February, 2023;
originally announced February 2023.
-
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.
-
Filamentary plasma eruptions and the heating and acceleration of electrons
Authors:
Heinz Isliker,
Andres Cathey,
Matthias Hoelzl,
Stanislas Pamela,
Loukas Vlahos
Abstract:
We present test-particle simulations of electrons during a nonlinear MHD simulation of a type-I edge localized mode (ELM) to explore the effect of an eruptive plasma filament on the kinetic level. The electrons are moderately heated and accelerated during the filamentary eruption on a fast time scale of the order of 0.5 ms. A clearly non-thermal tail is formed in the distribution of the kinetic en…
▽ More
We present test-particle simulations of electrons during a nonlinear MHD simulation of a type-I edge localized mode (ELM) to explore the effect of an eruptive plasma filament on the kinetic level. The electrons are moderately heated and accelerated during the filamentary eruption on a fast time scale of the order of 0.5 ms. A clearly non-thermal tail is formed in the distribution of the kinetic energy that is of power-law shape and reaches 90 keV for some particles. The acceleration is exclusively observed in the direction parallel to the magnetic field, i.e. with a clear preference in counter-current direction, and we show that the parallel electric field is the cause of the observed acceleration. Most particles that escape from the system leave at one distinct strike-line in the outer divertor leg at some time during their energization. The escaping high energy electrons in the tail of the energy distribution are not affected by collisions, they thus show characteristics of runaway electrons. The mean square displacement indicates that transport in energy space clearly is super-diffusive, and interpreting the acceleration process as a random walk, we find that the distributions of energy-increments exhibit exponential tails, and transport in energy space is equally important of convective (systematic) and diffusive (stochastic) nature. By analyzing the MHD simulations per se, it turns out that the histograms of the parallel electric field in the edge region exhibit power-law shapes, and this clearly non-Gaussian statistics is ultimately one of the reasons for the moderately anomalous phenomena of particle transport that we find in energy space.
△ Less
Submitted 21 October, 2022; v1 submitted 5 August, 2022;
originally announced August 2022.
-
Loss Landscape Engineering via Data Regulation on PINNs
Authors:
Vignesh Gopakumar,
Stanislas Pamela,
Debasmita Samaddar
Abstract:
Physics-Informed Neural Networks have shown unique utility in parameterising the solution of a well-defined partial differential equation using automatic differentiation and residual losses. Though they provide theoretical guarantees of convergence, in practice the required training regimes tend to be exacting and demanding. Through the course of this paper, we take a deep dive into understanding…
▽ More
Physics-Informed Neural Networks have shown unique utility in parameterising the solution of a well-defined partial differential equation using automatic differentiation and residual losses. Though they provide theoretical guarantees of convergence, in practice the required training regimes tend to be exacting and demanding. Through the course of this paper, we take a deep dive into understanding the loss landscapes associated with a PINN and how that offers some insight as to why PINNs are fundamentally hard to optimise for. We demonstrate how PINNs can be forced to converge better towards the solution, by way of feeding in sparse or coarse data as a regulator. The data regulates and morphs the topology of the loss landscape associated with the PINN to make it easily traversable for the minimiser. Data regulation of PINNs helps ease the optimisation required for convergence by invoking a hybrid unsupervised-supervised training approach, where the labelled data pushes the network towards the vicinity of the solution, and the unlabelled regime fine-tunes it to the solution.
△ Less
Submitted 16 May, 2022;
originally announced May 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.
-
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.
-
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.
-
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.
-
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.
-
Overview of recent physics results from MAST
Authors:
A Kirk,
J Adamek,
RJ Akers,
S Allan,
L Appel,
F Arese Lucini,
M Barnes,
T Barrett,
N Ben Ayed,
W Boeglin,
J Bradley,
P K Browning,
J Brunner,
P Cahyna,
M Carr,
F Casson,
M Cecconello,
C Challis,
IT Chapman,
S Chapman,
S Conroy,
N Conway,
WA Cooper,
M Cox,
N Crocker
, et al. (138 additional authors not shown)
Abstract:
New results from MAST are presented that focus on validating models in order to extrapolate to future devices. Measurements during start-up experiments have shown how the bulk ion temperature rise scales with the square of the reconnecting field. During the current ramp up models are not able to correctly predict the current diffusion. Experiments have been performed looking at edge and core turbu…
▽ More
New results from MAST are presented that focus on validating models in order to extrapolate to future devices. Measurements during start-up experiments have shown how the bulk ion temperature rise scales with the square of the reconnecting field. During the current ramp up models are not able to correctly predict the current diffusion. Experiments have been performed looking at edge and core turbulence. At the edge detailed studies have revealed how filament characteristic are responsible for determining the near and far SOL density profiles. In the core the intrinsic rotation and electron scale turbulence have been measured. The role that the fast ion gradient has on redistributing fast ions through fishbone modes has led to a redesign of the neutral beam injector on MAST Upgrade. In H-mode the turbulence at the pedestal top has been shown to be consistent with being due to electron temperature gradient modes. A reconnection process appears to occur during ELMs and the number of filaments released determines the power profile at the divertor. Resonant magnetic perturbations can mitigate ELMs provided the edge peeling response is maximised and the core kink response minimised. The mitigation of intrinsic error fields with toroidal mode number n>1 has been shown to be important for plasma performance.
△ Less
Submitted 18 November, 2016;
originally announced November 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.
-
Effect of resonant magnetic perturbations on low collisionality discharges in MAST and a comparison with ASDEX Upgrade
Authors:
A. Kirk,
W. Suttrop,
Yueqiang Liu,
I. T. Chapman,
P. Cahyna,
T. Eich,
C. Fuchs,
C. Ham,
J. R. Harrison,
M W. Jakubowski,
S. Pamela,
M. Peterka,
D. Ryan,
S. Saarelma,
R. Scannell,
A. J. Thornton,
M. Valovic,
B. Sieglin,
L. Barrera Orte,
M. Willensdorfer,
B. Kurzan,
R. Fischer
Abstract:
Sustained ELM mitigation has been achieved on MAST and AUG using RMPs with a range of toroidal mode numbers over a wide region of low to medium collisionality discharges. The ELM energy loss and peak heat loads at the divertor targets have been reduced. The ELM mitigation phase is typically associated with a drop in plasma density and overall stored energy. In one particular scenario on MAST, by c…
▽ More
Sustained ELM mitigation has been achieved on MAST and AUG using RMPs with a range of toroidal mode numbers over a wide region of low to medium collisionality discharges. The ELM energy loss and peak heat loads at the divertor targets have been reduced. The ELM mitigation phase is typically associated with a drop in plasma density and overall stored energy. In one particular scenario on MAST, by carefully adjusting the fuelling it has been possible to counteract the drop in density and to produce plasmas with mitigated ELMs, reduced peak divertor heat flux and with minimal degradation in pedestal height and confined energy. While the applied resonant magnetic perturbation field can be a good indicator for the onset of ELM mitigation on MAST and AUG there are some cases where this is not the case and which clearly emphasise the need to take into account the plasma response to the applied perturbations. The plasma response calculations show that the increase in ELM frequency is correlated with the size of the edge peeling-tearing like response of the plasma and the distortions of the plasma boundary in the X-point region.
△ Less
Submitted 24 December, 2014;
originally announced December 2014.
-
Electron kinetics inferred from observations of microwave bursts during edge localised modes in the Mega-Amp Spherical Tokamak
Authors:
S. J. Freethy,
K. G. McClements,
S. C. Chapman,
R. O. Dendy,
W. N. Lai,
S. J. P. Pamela,
V. F. Shevchenko,
R. G. L. Vann
Abstract:
Recent measurements of microwave and X-ray emission during edge localised mode (ELM) activity in tokamak plasmas provide a fresh perspective on ELM physics. It is evident that electron kinetics, which are not incorporated in standard (fluid) models for the instability that drives ELMs, play a key role in the new observations. These effects should be included in future models for ELMs and the ELM c…
▽ More
Recent measurements of microwave and X-ray emission during edge localised mode (ELM) activity in tokamak plasmas provide a fresh perspective on ELM physics. It is evident that electron kinetics, which are not incorporated in standard (fluid) models for the instability that drives ELMs, play a key role in the new observations. These effects should be included in future models for ELMs and the ELM cycle. The observed radiative effects paradoxically imply acceleration of electrons parallel to the magnetic field combined with rapid acquisition of perpendicular momentum. It is shown that this paradox can be resolved by the action of the anomalous Doppler instability which enables fast collective radiative relaxation, in the perpendicular direction, of electrons accelerated in the parallel direction by inductive electric fields generated by the initial ELM instability.
△ Less
Submitted 4 November, 2014;
originally announced November 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.
-
Towards understanding edge localised mode mitigation by resonant magnetic perturbations in MAST
Authors:
I. T. Chapman,
A. Kirk,
C. J. Ham,
J. R. Harrison,
Y. Q. Liu,
S. Saarelma,
R. Scannell,
A. J. Thornton,
M. Becoulet,
F. Orain,
W. A. Cooper,
S. Pamela,
MAST Team
Abstract:
Type-I Edge Localised Modes (ELMs) have been mitigated in MAST through the application of n = 3, 4 and 6 resonant magnetic perturbations (RMPs). For each toroidal mode number of the non-axisymmetric applied fields, the frequency of the ELMs has been increased significantly, and the peak heat flux on the divertor plates reduced commensurately. This increase in ELM frequency occurs despite a signifi…
▽ More
Type-I Edge Localised Modes (ELMs) have been mitigated in MAST through the application of n = 3, 4 and 6 resonant magnetic perturbations (RMPs). For each toroidal mode number of the non-axisymmetric applied fields, the frequency of the ELMs has been increased significantly, and the peak heat flux on the divertor plates reduced commensurately. This increase in ELM frequency occurs despite a significant drop in the edge pressure gradient, which would be expected to stabilise the peeling-ballooning modes thought to be responsible for type-I ELMs. Various mechanisms which could cause a destabilisation of the peeling-ballooning modes are presented, including pedestal widening, plasma rotation braking, three dimensional corrugation of the plasma boundary and the existence of radially extended lobe structures near to the X-point. This leads to a model aimed at resolving the apparent dichotomy of ELM control, that is to say ELM suppression occurring due to the pedestal pressure reduction below the peeling-ballooning stability boundary, whilst the reduction in pressure can also lead to ELM mitigation, which is ostensibly a destabilisation of peeling-ballooning modes. In the case of ELM mitigation, the pedestal broadening, 3d corrugation or lobes near the X-point degrade ballooning stability so much that the pedestal recovers rapidly to cross the new stability boundary at lower pressure more frequently, whilst in the case of suppression, the plasma parameters are such that the particle transport reduces the edge pressure below the stability boundary which is only mildly affected by negligible rotation braking, small edge corrugation or short, broad lobe structures.
△ Less
Submitted 16 May, 2013;
originally announced May 2013.