Showing 1–24 of 24 results for author: Saarelma, S

The JET hybrid H-mode scenario from a pedestal turbulence perspective

Authors: L. A. Leppin, T. Görler, L. Frassinetti, S. Saarelma, J. Hobirk, F. Jenko, JET contributors

Abstract: Turbulent transport is a decisive factor in determining the pedestal structure of H-modes. Here, we present the first comprehensive characterization of gyrokinetic turbulent transport in a JET hybrid H-mode pedestal. Local, linear simulations are performed to identify instabilities and global, nonlinear electromagnetic simulations reveal the turbulent heat and particle flux structure of the pedest… ▽ More Turbulent transport is a decisive factor in determining the pedestal structure of H-modes. Here, we present the first comprehensive characterization of gyrokinetic turbulent transport in a JET hybrid H-mode pedestal. Local, linear simulations are performed to identify instabilities and global, nonlinear electromagnetic simulations reveal the turbulent heat and particle flux structure of the pedestal. Our analysis focuses on the Deuterium reference discharge \#97781 performed in the scenario development for the Deuterium-Tritium campaign. We find the pedestal top transport to be dominated by ion temperature gradient (ITG) modes. In the pedestal center turbulent ion transport is suppressed and electron transport is driven by multi-faceted electron temperature gradient (ETG) modes, which extend down to ion-gyroradius scales. A strong impact of $E\times B$ shear on the absolute turbulence level is confirmed by the global, nonlinear simulations. Furthermore, impurities are shown to reduce the main ion transport. Dedicated density and ion temperature profile variations test the sensitivity of the results and do not find strong differences in the turbulent transport in more reactor-like conditions. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Flat-top plasma operational space of the STEP power plant

Authors: E. Tholerus, F. J. Casson, S. P. Marsden, T. Wilson, D. Brunetti, P. Fox, S. J. Freethy, T. C. Hender, S. S. Henderson, A. Hudoba, K. K. Kirov, F. Koechl, H. Meyer, S. I. Muldrew, C. Olde, B. S. Patel, C. M. Roach, S. Saarelma, G. Xia

Abstract: STEP is a spherical tokamak prototype power plant that is being designed to demonstrate net electric power. The design phase involves the exploitation of plasma models to optimise fusion performance subject to satisfying various physics and engineering constraints. A modelling workflow, including integrated core plasma modelling, MHD stability analysis, SOL and pedestal modelling, coil set and fre… ▽ More STEP is a spherical tokamak prototype power plant that is being designed to demonstrate net electric power. The design phase involves the exploitation of plasma models to optimise fusion performance subject to satisfying various physics and engineering constraints. A modelling workflow, including integrated core plasma modelling, MHD stability analysis, SOL and pedestal modelling, coil set and free boundary equilibrium solvers, and whole plant design, has been developed to specify the design parameters and to develop viable scenarios. The integrated core plasma model JETTO is used to develop individual flat-top operating points that satisfy imposed criteria for fusion power performance within operational constraints. Key plasma parameters such as normalised beta, Greenwald density fraction, auxiliary power and radiated power have been scanned to scope the operational space and to derive a collection of candidate non-inductive flat-top points. The assumed auxiliary heating and current drive is either from electron cyclotron systems only or a combination of electron cyclotron and electron Bernstein waves. At present stages of transport modelling, there is a large uncertainty in overall confinement for relevant parameter regimes. For each of the two auxiliary heating and current drive systems scenarios, two candidate flat-top points have been developed based on different confinement assumptions, totalling to four operating points. A lower confinement assumption generally suggests operating points in high-density, high auxiliary power regimes, whereas higher confinement would allow access to a broader parameter regime in density and power while maintaining target fusion power performance. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 44 pages, 15 figures. Submitted to Nucl. Fusion (IOP)

EuroPED-NN: Uncertainty aware surrogate model

Authors: A. Panera Alvarez, A. Ho, A. Jarvinen, S. Saarelma, S. Wiesen, JET Contributors

Abstract: This work successfully generates uncertainty aware surrogate models, via the Bayesian neural network with noise contrastive prior (BNN-NCP) technique, of the EuroPED plasma pedestal model using data from the JET-ILW pedestal database and subsequent model evaluations. All this conform EuroPED-NN. The BNN-NCP technique is proven to be a good fit for uncertainty aware surrogate models, matching the o… ▽ More This work successfully generates uncertainty aware surrogate models, via the Bayesian neural network with noise contrastive prior (BNN-NCP) technique, of the EuroPED plasma pedestal model using data from the JET-ILW pedestal database and subsequent model evaluations. All this conform EuroPED-NN. The BNN-NCP technique is proven to be a good fit for uncertainty aware surrogate models, matching the output results as a regular neural network, providing prediction's confidence as uncertainties, and highlighting the out of distribution (OOD) regions using surrogate model uncertainties. This provides critical insights into model robustness and reliability. EuroPED-NN has been physically validated, first, analyzing electron density $n_e\!\left(ψぷさい_{\text{pol}}=0.94\right)$ with respect to increasing plasma current, $I_p$, and second, validating the $Δでるた-βべーた_{p,ped}$ relation associated with the EuroPED model. Affirming the robustness of the underlying physics learned by the surrogate model. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Flexible, integrated modeling of tokamak stability, transport, equilibrium, and pedestal physics

Authors: B. C. Lyons, J. McClenaghan, T. Slendebroek, O. Meneghini, T. F. Neiser, S. P. Smith, D. B. Weisberg, E. A. Belli, J. Candy, J. M. Hanson, L. L. Lao, N. C. Logan, S. Saarelma, O. Sauter, P. B. Snyder, G. M. Staebler, K. E. Thome, A. D. Turnbull

Abstract: The STEP (Stability, Transport, Equilibrium, and Pedestal) integrated-modeling tool has been developed in OMFIT to predict stable, tokamak equilibria self-consistently with core-transport and pedestal calculations. STEP couples theory-based codes to integrate a variety of physics, including MHD stability, transport, equilibrium, pedestal formation, and current-drive, heating, and fueling. The inpu… ▽ More The STEP (Stability, Transport, Equilibrium, and Pedestal) integrated-modeling tool has been developed in OMFIT to predict stable, tokamak equilibria self-consistently with core-transport and pedestal calculations. STEP couples theory-based codes to integrate a variety of physics, including MHD stability, transport, equilibrium, pedestal formation, and current-drive, heating, and fueling. The input/output of each code is interfaced with a centralized ITER-IMAS data structure, allowing codes to be run in any order and enabling open-loop, feedback, and optimization workflows. This paradigm simplifies the integration of new codes, making STEP highly extensible. STEP has been verified against a published benchmark of six different integrated models. Core-pedestal calculations with STEP have been successfully validated against individual DIII-D H-mode discharges and across more than 500 discharges of the $H_{98,y2}$ database, with a mean error in confinement time from experiment less than 19%. STEP has also reproduced results in less conventional DIII-D scenarios, including negative-central-shear and negative-triangularity plasmas. Predictive STEP modeling has been used to assess performance in several tokamak reactors. Simulations of a high-field, large-aspect-ratio reactor show significantly lower fusion power than predicted by a zero-dimensional study, demonstrating the limitations of scaling-law extrapolations. STEP predictions have found promising EXCITE scenarios, including a high-pressure, 80%-bootstrap-fraction plasma. ITER modeling with STEP has shown that pellet fueling enhances fusion gain in both the baseline and advanced-inductive scenarios. Finally, STEP predictions for the SPARC baseline scenario are in good agreement with published results from the physics basis. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 15 pages, 11 figures Associated with invited talk at 63nd Annual Meeting of the APS Division of Plasma Physics: https://meetings.aps.org/Meeting/DPP21/Session/NI02.1 . The following article has been submitted to Physics of Plasmas. After it is published, it will be found at https://publishing.aip.org/resources/librarians/products/journals/

H-mode inhibition in negative triangularity tokamak reactor plasmas

Authors: A. O. Nelson, C. Paz-Soldan, S. Saarelma

Abstract: Instability to high toroidal mode number ($n$) ballooning modes has been proposed as the primary gradient-limiting mechanism for tokamak equilibria with negative triangularity ($δでるた$) shaping, preventing access to strong H-mode regimes when $δでるた\ll0$. To understand how this mechanism extrapolates to reactor conditions, we model the infinite-$n$ ballooning stability as a function of internal profiles a… ▽ More Instability to high toroidal mode number ($n$) ballooning modes has been proposed as the primary gradient-limiting mechanism for tokamak equilibria with negative triangularity ($δでるた$) shaping, preventing access to strong H-mode regimes when $δでるた\ll0$. To understand how this mechanism extrapolates to reactor conditions, we model the infinite-$n$ ballooning stability as a function of internal profiles and equilibrium shape using a combination of the CHEASE and BALOO codes. While the critical $δでるた$ required for avoiding $2^\mathrm{nd}$ stability to high-$n$ modes is observed to depend in a complicated way on various shaping parameters, including the equilibrium aspect ratio, elongation and squareness, equilibria with negative triangularity are robustly prohibited from accessing the $2^\mathrm{nd}$ stability region, offering the prediction that that negative triangularity reactors should maintain L-mode-like operation. In order to access high-$n$ $2^\mathrm{nd}$ stability, the local shear over the entire bad curvature region must be sufficiently negative to overcome curvature destabilization on the low field side. Scalings of the ballooning-limited pedestal height are provided as a function of plasma parameters to aid future scenario design. △ Less

Submitted 26 April, 2022; originally announced April 2022.

Three-Dimensional Inhomogeneity of Electron-Temperature-Gradient Turbulence in the Edge of Tokamak Plasmas

Authors: J. F. Parisi, F. I. Parra, C. M. Roach, M. R. Hardman, A. A. Schekochihin, I. G. Abel, N. Aiba, J. Ball, M. Barnes, B. Chapman-Oplopoiou, D. Dickinson, W. Dorland, C. Giroud, D. R. Hatch, J. C. Hillesheim, J. Ruiz Ruiz, S. Saarelma, D. St-Onge

Abstract: Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a… ▽ More Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport -- in our case by $\sim$40 \% -- through multiscale interactions. △ Less

Submitted 2 July, 2022; v1 submitted 1 March, 2022; originally announced March 2022.

Comments: 17 pages, 14 figures

Microtearing modes as the source of magnetic fluctuations in the JET pedestal

Authors: D. R. Hatch, M. Kotschenreuther, S. M. Mahajan, M. J. Pueschel, C. Michoski, G. Merlo, E. Hassan, A. R. Field, L. Frassinetti, C. Giroud, J. C. Hillesheim, C. F. Maggi, C. Perez von Thun, C. M. Roach, S. Saarelma, D. Jarema, F. Jenko, JET contributors

Abstract: We report on a detailed study of magnetic fluctuations in the JET pedestal, employing basic theoretical considerations, gyrokinetic simulations, and experimental fluctuation data, to establish the physical basis for their origin, role, and distinctive characteristics. We demonstrate quantitative agreement between gyrokinetic simulations of microtearing modes (MTMs) and two magnetic frequency bands… ▽ More We report on a detailed study of magnetic fluctuations in the JET pedestal, employing basic theoretical considerations, gyrokinetic simulations, and experimental fluctuation data, to establish the physical basis for their origin, role, and distinctive characteristics. We demonstrate quantitative agreement between gyrokinetic simulations of microtearing modes (MTMs) and two magnetic frequency bands with corresponding toroidal mode numbers n=4 and 8. Such disparate fluctuation scales, with substantial gaps between toroidal mode numbers, are commonly observed in pedestal fluctuations. Here we provide a clear explanation, namely the alignment of the relevant rational surfaces (and not others) with the peak in the omega star profile, which is localized in the steep gradient region of the pedestal. We demonstrate that a global treatment is required to capture this effect. Nonlinear simulations suggest that the MTM fluctuations produce experimentally-relevant transport levels and saturate by relaxing the background electron temperature gradient, slightly downshifting the fluctuation frequencies from the linear predictions. Scans in collisionality are compared with simple MTM dispersion relations. At the experimental points considered, MTM growth rates can either increase or decrease with collision frequency depending on the parameters thus defying any simple characterization of collisionality dependence. △ Less

Submitted 14 July, 2020; originally announced July 2020.

Gyrokinetic analysis and simulation of pedestals, to identify the culprits for energy losses using fingerprints

Authors: M. Kotschenreuther, X. Liu, D. R. Hatch, S. Mahajan, L. Zheng, A. Diallo, R. Groebner, the DIII-D TEAM, J. C. Hillesheim, C. F. Maggi, C. Giroud, F. Koechl, V. Parail, S. Saarelma, E. Solano, JET Contributors, A. Chankin

Abstract: Fusion performance in tokamaks hinges critically on the efficacy of the Edge Transport Barrier (ETB) at suppressing energy losses. The new concept of fingerprints is introduced to identify the instabilities that cause the transport losses in the ETB of many of today's experiments, from widely posited candidates. Analysis of the Gyrokinetic-Maxwell equations, and gyrokinetic simulations of experime… ▽ More Fusion performance in tokamaks hinges critically on the efficacy of the Edge Transport Barrier (ETB) at suppressing energy losses. The new concept of fingerprints is introduced to identify the instabilities that cause the transport losses in the ETB of many of today's experiments, from widely posited candidates. Analysis of the Gyrokinetic-Maxwell equations, and gyrokinetic simulations of experiments, find that each mode type produces characteristic ratios of transport in the various channels: density, heat and impurities. This, together with experimental observations of transport in some channel, or, of the relative size of the driving sources of channels, can identify or determine the dominant modes causing energy transport. In multiple ELMy H-mode cases that are examined, these fingerprints indicate that MHD-like modes are apparently not the dominant agent of energy transport; rather, this role is played by Micro-Tearing Modes (MTM) and Electron Temperature Gradient (ETG) modes, and in addition, possibly Ion Temperature Gradient (ITG)/Trapped Electron Modes (ITG/TEM) on JET. MHD-like modes may dominate the electron particle losses. Fluctuation frequency can also be an important means of identification, and is often closely related to the transport fingerprint. The analytical arguments unify and explain previously disparate experimental observations on multiple devices, including DIII-D, JET and ASDEX-U, and detailed simulations of two DIII-D ETBs also demonstrate and corroborate this. △ Less

Submitted 24 March, 2019; originally announced March 2019.

Direct Gyrokinetic Comparison of Pedestal Transport in JET with Carbon and ITER-Like Walls

Authors: D. R. Hatch, M. Kotschenreuther, S. M. Mahajan, G. Merlo, A. R. Field, C. Giroud, J. C. Hillesheim, C. F. Maggi, C. Perez von Thun, C. M. Roach, S. Saarelma

Abstract: This paper compares the gyrokinetic instabilities and transport in two representative JET pedestals, one (pulse 78697) from the JET configuration with a carbon wall (C) and another (pulse 92432) from after the installation of JET's ITER-like Wall (ILW). The discharges were selected for a comparison of JET-ILW and JET-C discharges with good confinement at high current (3 MA, corresponding also to l… ▽ More This paper compares the gyrokinetic instabilities and transport in two representative JET pedestals, one (pulse 78697) from the JET configuration with a carbon wall (C) and another (pulse 92432) from after the installation of JET's ITER-like Wall (ILW). The discharges were selected for a comparison of JET-ILW and JET-C discharges with good confinement at high current (3 MA, corresponding also to low $ρろー_*$) and retain the distinguishing features of JET-C and JET-ILW, notably, decreased pedestal top temperature for JET-ILW. A comparison of the profiles and heating power reveals a stark qualitative difference between the discharges: the JET-ILW pulse (92432) requires twice the heating power, at a gas rate of $1.9 \times 10^{22}e/s$, to sustain roughly half the temperature gradient of the JET-C pulse (78697), operated at zero gas rate. This points to heat transport as a central component of the dynamics limiting the JET-ILW pedestal and reinforces the following emerging JET-ILW pedestal transport paradigm, which is proposed for further examination by both theory and experiment. ILW conditions modify the density pedestal in ways that decrease the normalized pedestal density gradient $a/L_n$, often via an outward shift of the density pedestal. This is attributable to some combination of direct metal wall effects and the need for increased fueling to mitigate tungsten contamination. The modification to the density profile increases $ηいーた= L_n/L_T$ , thereby producing more robust ion temperature gradient (ITG) and electron temperature gradient driven instability. The decreased pedestal gradients for JET-ILW (92432) also result in a strongly reduced $E \times B$ shear rate, further enhancing the ion scale turbulence. Collectively, these effects limit the pedestal temperature and demand more heating power to achieve good pedestal performance. △ Less

Submitted 6 March, 2019; originally announced March 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.

Comments: 34 pages, 10 figures. This is an author-created, un-copyedited version of an article submitted for publication in Nuclear Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

Improved Confinement in JET High {beta} Plasmas with an ITER-Like Wall

Authors: C. D. Challis, J. Garcia, M. Beurskens, P. Buratti, E. Delabie, P. Drewelow, L. Frassinetti, C. Giroud, N. Hawkes, J. Hobirk, E. Joffrin, D. Keeling, D. B. King, C. F. Maggi, J. Mailloux, C. Marchetto, D. McDonald, I. Nunes, G. Pucella, S. Saarelma, J. Simpson

Abstract: The replacement of the JET carbon wall (C-wall) by a Be/W ITER-like wall (ILW) has affected the plasma energy confinement. To investigate this, experiments have been performed with both the C-wall and ILW to vary the heating power over a wide range for plasmas with different shapes. The replacement of the JET carbon wall (C-wall) by a Be/W ITER-like wall (ILW) has affected the plasma energy confinement. To investigate this, experiments have been performed with both the C-wall and ILW to vary the heating power over a wide range for plasmas with different shapes. △ Less

Submitted 16 January, 2015; originally announced January 2015.

Comments: 50 pages, 32 figures. This is an author-created, un-copyedited version of an article submitted for publication in Nuclear Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

Journal ref: Nuclear Fusion, Vol.55, No.5, May 2015, pp.053031

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.

Comments: 31 pages, 28 figures. This is an author-created, un-copyedited version of an article submitted for publication in Nuclear Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

Journal ref: Nuclear Fusion, Vol.55, No.4, April 2015, pp.043011

Contrasting H-mode behaviour with deuterium fuelling and nitrogen seeding in the all-carbon and metallic versions of JET

Authors: G. P. Maddison, C. Giroud, B. Alper, G. Arnoux, I. Balboa, M. N. A. Beurskens, A. Boboc, S. Brezinsek, M. Brix, M. Clever, R. Coelho, J. W. Coenen, I. Coffey, P. C. da Silva Aresta Belo, S. Devaux, P. Devynck, T. Eich, R. C. Felton, J. Flanagan, L. Frassinetti, L. Garzotti, M. Groth, S. Jachmich, A. Järvinen, E. Joffrin , et al. (26 additional authors not shown)

Abstract: The former all-carbon wall on JET has been replaced with beryllium in the main torus and tungsten in the divertor to mimic the surface materials envisaged for ITER. Comparisons are presented between Type I H-mode characteristics in each design by examining respective scans over deuterium fuelling and impurity seeding, required to ameliorate exhaust loads both in JET at full capability and in ITER. The former all-carbon wall on JET has been replaced with beryllium in the main torus and tungsten in the divertor to mimic the surface materials envisaged for ITER. Comparisons are presented between Type I H-mode characteristics in each design by examining respective scans over deuterium fuelling and impurity seeding, required to ameliorate exhaust loads both in JET at full capability and in ITER. △ Less

Submitted 11 June, 2014; originally announced June 2014.

Comments: 55 pages, 15 figures

Journal ref: Nuclear Fusion, Vol.54, No.7, July 2014, p.073016

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.

Comments: 35 pages, 9 figures

Journal ref: Nuclear Fusion, Vol.54, No.11, November 2014, pp.114012

Modelling of three dimensional equilibrium and stability of MAST plasmas with magnetic perturbations using VMEC and COBRA

Authors: C J Ham, I T Chapman, A Kirk, S Saarelma

Abstract: It is known that magnetic perturbations can mitigate edge localized modes (ELMs) in experiments, for example MAST (Kirk et al 2013 Nucl. Fusion 53 043007). One hypothesis is that the magnetic perturbations cause a three dimensional corrugation of the plasma and this corrugated plasma has different stability properties to peeling-ballooning modes compared to an axisymmetric plasma. It has been show… ▽ More It is known that magnetic perturbations can mitigate edge localized modes (ELMs) in experiments, for example MAST (Kirk et al 2013 Nucl. Fusion 53 043007). One hypothesis is that the magnetic perturbations cause a three dimensional corrugation of the plasma and this corrugated plasma has different stability properties to peeling-ballooning modes compared to an axisymmetric plasma. It has been shown in an up-down symmetric plasma that magnetic perturbations in tokamaks will break the usual axisymmetry of the plasma causing three dimensional displacements (Chapman et al 2012 Plasma Phys. Control. Fusion 54 105013). We produce a free boundary three-dimensional equilibrium of a lower single null MAST relevant plasma using VMEC (S P Hirshman and J C Whitson 1983 Phys. Fluids 26 3553). The current and pressure profiles used for the modelling are similar to those deduced from axisymmetric analysis of experimental data with ELMs. We focus on the effect of applying $n=3$ and $n=6$ magnetic perturbations using the RMP coils. A midplane displacement of over $\pm 1$ cm is seen when the full current is applied. The current in the coils is scaned and a linear relationship between coil current and midplane displacement is found. The effect of this non-axisymmetric equilibrium on infinite $n$ ballooning stability is investigated using COBRA (R Sanchez et al 2000 J. Comput. Phys. 161 576-588). △ Less

Submitted 11 December, 2013; originally announced December 2013.

Comments: 13 pages, 10 figures

Journal ref: Physics of Plasmas, Vol.21, No.10, October 2014, pp.102501

Comparison of BES measurements of ion-scale turbulence with direct, gyrokinetic simulations of MAST L-mode plasmas

Authors: A R Field, D Dunai, Y-c Ghim, P Hill, B McMillan, C M Roach, S Saarelma, A A Schekochihin, S Zoletnik, the MAST team

Abstract: Observations of ion-scale (k_y*rho_i <= 1) density turbulence of relative amplitude dn_e/n_e <= 0.2% are available on the Mega Amp Spherical Tokamak (MAST) using a 2D (8 radial x 4 poloidal channel) imaging Beam Emission Spectroscopy (BES) diagnostic. Spatial and temporal characteristics of this turbulence, i.e., amplitudes, correlation times, radial and perpendicular correlation lengths and appar… ▽ More Observations of ion-scale (k_y*rho_i <= 1) density turbulence of relative amplitude dn_e/n_e <= 0.2% are available on the Mega Amp Spherical Tokamak (MAST) using a 2D (8 radial x 4 poloidal channel) imaging Beam Emission Spectroscopy (BES) diagnostic. Spatial and temporal characteristics of this turbulence, i.e., amplitudes, correlation times, radial and perpendicular correlation lengths and apparent phase velocities of the density contours, are determined by means of correlation analysis. For a low-density, L-mode discharge with strong equilibrium flow shear exhibiting an internal transport barrier (ITB) in the ion channel, the observed turbulence characteristics are compared with synthetic density turbulence data generated from global, non-linear, gyro-kinetic simulations using the particle-in-cell (PIC) code NEMORB. This validation exercise highlights the need to include increasingly sophisticated physics, e.g., kinetic treatment of trapped electrons, equilibrium flow shear and collisions, to reproduce most of the characteristics of the observed turbulence. Even so, significant discrepancies remain: an underprediction by the simulations of the turbulence amplituide and heat flux at plasma periphery and the finding that the correlation times of the numerically simulated turbulence are typically two orders of magnitude longer than those measured in MAST. Comparison of these correlation times with various linear timescales suggests that, while the measured turbulence is strong and may be `critically balanced', the simulated turbulence is weak. △ Less

Submitted 7 January, 2014; v1 submitted 29 June, 2013; originally announced July 2013.

Comments: 27 pages, 11 figures

Understanding ELM mitigation by resonant magnetic perturbations on MAST

Authors: A. Kirk, I. T. Chapman, Yueqiang Liu, P. Cahyna, P. Denner, G. Fishpool, C. J. Ham, J. R. Harrison, Yunfeng Liang, E. Nardon, S. Saarelma, R. Scannell, A. J. Thornton, the MAST team

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

Submitted 20 June, 2013; originally announced June 2013.

Comments: 41 pages, 19 figures. arXiv admin note: text overlap with arXiv:1305.3061

Journal ref: Nuclear Fusion, 53 043007 (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.

Comments: 23 pages, 12 figures. Copyright (2013) United Kingdom Atomic Energy Authority. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics

Journal ref: Phys. Plasmas 20, 056101 (2013)

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

Authors: A. Kirk, I. T. Chapman, J. Harrison, Yueqiang Liu, E. Nardon, S. Saarelma, R. Scannell, A. J. Thornton, the MAST team

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

Submitted 14 May, 2013; originally announced May 2013.

Comments: 48 pages, 23 figures. This is an author-created, un-copyedited version of an article accepted for publication in Plasma Phys. Control. Fusion. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

Journal ref: Plasma Phys. Control. Fusion 55 015006 (2013)

MHD and Gyro-kinetic Stability of JET Pedestals

Authors: S. Saarelma, M. N. A. Beurskens, D. Dickinson, L. Frassinetti, M. J. Leyland, C. M. Roach, EFDA-JET contributors

Abstract: The pedestal profile measurements in high triangularity JET plasmas show that with low fuelling the pedestal width decreases during the ELM cycle and with high fuelling it stays constant. In the low fuelling case the pedestal pressure gradient keeps increasing until the ELM crash and in the high fuelling case it initially increases then saturates during the ELM cycle. Stability analysis reveals… ▽ More The pedestal profile measurements in high triangularity JET plasmas show that with low fuelling the pedestal width decreases during the ELM cycle and with high fuelling it stays constant. In the low fuelling case the pedestal pressure gradient keeps increasing until the ELM crash and in the high fuelling case it initially increases then saturates during the ELM cycle. Stability analysis reveals that both JET plasmas become unstable to finite-n ideal MHD peeling-ballooning modes at the end of the ELM cycle. During the ELM cycle, infinite-n ideal MHD ballooning modes and kinetic ballooning modes are found to be locally stable in most of the steep pressure gradient region of the pedestal owing to the large bootstrap current, but to be locally unstable in a narrow region of plasma at the extreme edge. Unstable micro-tearing modes are found at the JET pedestal top, but they are sub-dominant to ion temperature gradient modes. They are insensitive to collisionality and stabilised by increasing density gradient. △ Less

Submitted 4 December, 2013; v1 submitted 14 January, 2013; originally announced January 2013.

Comments: submitted to Nuclear Fusion, 23 pages, for EFDA-JET contributors see the Appendix of F. Romanelli et al., Proceedings of the 24th IAEA Fusion Energy Conference 2012, San Diego, USA

Journal ref: 2013 Nucl. Fusion 53 123012

Microtearing modes at the top of the pedestal

Authors: D. Dickinson, C. M. Roach, S. Saarelma, R. Scannell, A. Kirk, H. R. Wilson

Abstract: Microtearing modes (MTMs) are unstable in the shallow gradient region just inside the top of the pedestal in the spherical tokamak experiment MAST, and may play an important role in the pedestal evolution. The linear properties of these instabilities are compared with MTMs deeper inside the core, and further detailed investigations expose the basic drive mechanism. MTMs near the MAST pedestal top… ▽ More Microtearing modes (MTMs) are unstable in the shallow gradient region just inside the top of the pedestal in the spherical tokamak experiment MAST, and may play an important role in the pedestal evolution. The linear properties of these instabilities are compared with MTMs deeper inside the core, and further detailed investigations expose the basic drive mechanism. MTMs near the MAST pedestal top are not well described by existing theories. In particular the growth rate of the dominant edge MTM does not peak at a finite collision frequency, as frequently reported for MTMs further into the core. Our study suggests that the edge MTM is driven by a collisionless trapped particle mechanism that is sensitive to magnetic drifts. This drive is enhanced in the outer region of MAST by a high magnetic shear and a high trapped particle fraction. Observations of similar modes on conventional aspect ratio devices suggests this drive mechanism may be somewhat ubiquitous towards the edge of current day and future hot tokamaks. △ Less

Submitted 9 January, 2013; v1 submitted 17 September, 2012; originally announced September 2012.

Comments: 14 pages, 10 figures (21 unique graphics)

Kinetic instabilities that limit βべーた in the edge of a tokamak plasma: a picture of an H-mode pedestal

Authors: D. Dickinson, C. M. Roach, S. Saarelma, R. Scannell, A. Kirk, H. R. Wilson

Abstract: Plasma equilibria reconstructed from the Mega-Amp Spherical Tokamak (MAST) have sufficient resolution to capture plasma evolution during the short period between edge-localized modes (ELMs). Immediately after the ELM steep gradients in pressure, P, and density, ne, form pedestals close to the separatrix, and they then expand into the core. Local gyrokinetic analysis over the ELM cycle reveals the… ▽ More Plasma equilibria reconstructed from the Mega-Amp Spherical Tokamak (MAST) have sufficient resolution to capture plasma evolution during the short period between edge-localized modes (ELMs). Immediately after the ELM steep gradients in pressure, P, and density, ne, form pedestals close to the separatrix, and they then expand into the core. Local gyrokinetic analysis over the ELM cycle reveals the dominant microinstabilities at perpendicular wavelengths of the order of the ion Larmor radius. These are kinetic ballooning modes (KBMs) in the pedestal and microtearing modes (MTMs) in the core close to the pedestal top. The evolving growth rate spectra, supported by gyrokinetic analysis using artificial local equilibrium scans, suggest a new physical picture for the formation and arrest of this pedestal. △ Less

Submitted 14 May, 2012; v1 submitted 4 October, 2011; originally announced October 2011.

Comments: Final version as it appeared in PRL (March 2012). Minor improvements include: shortened abstract, and better colour table for figures. 4 pages, 6 figures

Towards the construction of a model to describe the inter-ELM evolution of the pedestal on MAST

Authors: D. Dickinson, S. Saarelma, R. Scannell, A. Kirk, C. M. Roach, H. R. Wilson

Abstract: Pedestal profiles that span the ELM cycle have been obtained and used to test the idea that the pedestal pressure gradient in MAST is limited by the onset of Kinetic Ballooning Modes (KBMs). During the inter-ELM period of a regularly type I ELM-ing discharge on MAST, the pressure pedestal height and width increase together while the pressure gradient increases by only 15 % during the ELM cycle. St… ▽ More Pedestal profiles that span the ELM cycle have been obtained and used to test the idea that the pedestal pressure gradient in MAST is limited by the onset of Kinetic Ballooning Modes (KBMs). During the inter-ELM period of a regularly type I ELM-ing discharge on MAST, the pressure pedestal height and width increase together while the pressure gradient increases by only 15 % during the ELM cycle. Stability analyses show that the pedestal region over which infinite-n ballooning modes are unstable also broadens during the ELM cycle. To test the relationship between the width of the region that is unstable to n = \infty ideal magnetohydrodynamic ballooning modes and KBMs the gyrokinetic code, GS2, has been used for microstability analysis of the edge plasma region in MAST. The gyrokinetic simulations find that KBM modes with twisting parity are the dominant microinstabilities in the steep pedestal region, with a transition to more slowly growing tearing parity modes in the shallower pressure gradient core region immediately inside the pedestal top. The region over which KBMs are unstable increases during the ELM cycle, and a good correlation is found between the regions unstable to KBMs and infinite-n ideal ballooning modes. △ Less

Submitted 30 September, 2011; v1 submitted 15 July, 2011; originally announced July 2011.

Comments: Version accepted and to appear in PPCF October 2011

Journal ref: Plasma Physics and Controlled Fusion 53 115010 (2011)

Observation of confined current ribbon in JET plasmas

Authors: E. R. Solano, P. J. Lomas, B. Alper, G. S. Xu, Y. Andrew, G. Arnoux, A. Boboc, L. Barrera, P. Belo, M. N. A. Beurskens, M. Brix, K. Crombe, E. de la Luna, S. Devaux, T. Eich, S. Gerasimov, C. Giroud, D. Harting, D. Howell, A. Huber, G. Kocsis, A. Korotkov, A. Lopez-Fraguas, M. F. F. Nave, E. Rachlew , et al. (7 additional authors not shown)

Abstract: we report the identification of a localised current structure inside the JET plasma. It is a field aligned closed helical ribbon, carrying current in the same direction as the background current profile (co-current), rotating toroidally with the ion velocity (co-rotating). It appears to be located at a flat spot in the plasma pressure profile, at the top of the pedestal. The structure appears sp… ▽ More we report the identification of a localised current structure inside the JET plasma. It is a field aligned closed helical ribbon, carrying current in the same direction as the background current profile (co-current), rotating toroidally with the ion velocity (co-rotating). It appears to be located at a flat spot in the plasma pressure profile, at the top of the pedestal. The structure appears spontaneously in low density, high rotation plasmas, and can last up to 1.4 s, a time comparable to a local resistive time. It considerably delays the appearance of the first ELM. △ Less

Submitted 30 October, 2009; originally announced October 2009.

Comments: 10 pages, 6 figures

Journal ref: Phys. Rev. Lett. 104, 185003 (2010)