Multiscale turbulent dynamics of tokamak plasmas

Lead Research Organisation: University of York
Department Name: Physics

Abstract

Plasma turbulence underpins a wide range of phenomena, including the formation of stars and galaxies; the properties of the solar wind, and - the focus of this programme - the confinement of plasmas in tokamaks. It is complicated by feedback mechanisms that couple space and time scales spanning several orders of magnitude. The full problem is extremely challenging, and so to make progress for real world applications we must develop reduced models that capture the essential physics. The goal of our proposed programme is to address this by advancing our understanding of these multi-scale interactions at a fundamental science level. This will be achieved by coupling analytic theory, advanced computation and experimental capabilities, including the newly upgraded MAST-U tokamak.

Plasma turbulence is complicated by the fact that there are at least two types of interacting "fluids" - electrons and ions - and these are charged. Fluctuations in density therefore drive charge separation and hence fluctuations in the electrostatic field, while fluctuations in velocity drive currents and hence fluctuations in the magnetic field. These fields then couple the relative motions of the electron and ion "fluids". The situation is further complicated by the rich variety of waves that a magnetised plasma supports, and the resonances that exist when the phase velocity of a wave matches the particle velocity. To properly treat these resonances requires knowledge of the particle velocity distribution; this, in turn, requires either a kinetic or an advanced fluid approach - a daunting task.

Turbulence, typically at the millimetre-centimetre scale in tokamaks, interacts in a complex way with the global equilibrium profiles (density, temperature and flow gradients, for example), which are on the metre-scale. To quantify the complex, multi-scale feedback mechanisms between tokamak plasma turbulence and profiles, and so provide a predictive capability for the quasi-steady final states, we will address and integrate a number of topics. We will first learn how mean flows interact with electrostatic turbulence (ie neglecting fluctuations in the magnetic field), requiring coupling between fluctuations with characteristic scales ranging from the electron Larmor radius (sub-mm) through to the ion Larmor radius (few mm) and beyond (cm), to the system length scale of the profiles (m). Our new theory and simulations will inform experiments on MAST-U, exploiting two diagnostic instruments already planned for the device (beam emission spectroscopy and doppler back-scattering). It is likely there will be gaps in the wavelength range that these instruments can measure, so we anticipate a need to develop and install a new microwave imaging system. This will be designed using knowledge gained from the early phase of the programme, and deployed for further experiments towards the end.

Understanding of electromagnetic turbulence is less developed and new theoretical models will be required. Building on the knowledge gained from the electrostatic turbulence, we will seek to again understand the multi-scale interactions and feedbacks, including flows. However, now the situation is more complicated as electromagnetic turbulence can drive large scale currents, modifying the magnetic field which confines the plasma, and coupling into large scale electromagnetic modes.

A key motivation is to optimise tokamak plasmas for fusion performance, and this requires us to understand the impact of fast particles. These can drive turbulence directly through the instabilities they excite, or influence the turbulence driven by the thermal particles. Our simulations will assess the impact of the fast particles created by the neutral beam heating systems on MAST-U, and also the impact of energetic alpha particles from fusion reactions on future devices like ITER, as well as experiments planned on JET with the deuterium-tritium mix fusion fuel.

Planned Impact

The largest impact will be on the fusion programme. This impact will occur at several levels, as follows.
(1) A greater understanding of turbulence and confinement on MAST-U will enable us to identify a wider range of plasma scenarios on MAST-U. These, in turn, will broaden the range of plasma parameters achievable on MAST-U, further extending its flexibility for a wide variety of plasma science experiments. This will benefit many UK researchers, but will also make the device even more attractive for international participation which will increase opportunities for international funding for MAST-U operations (EUROfusion, and the wider international community).
(2) Our improved understanding will feed into the development of operational scenarios for the DT operation of JET and ITER. This will, in turn, help maintain UK influence in ITER, and contribute towards ITER meeting its fusion objectives in a timely way. This is important, as the design of the first demonstration power plants, and hence the commercialisation of fusion energy, cannot be finalised until the first results are obtained from burning fusion plasmas in ITER.
(3) The size and capital cost of a tokamak fusion power plant depend to a large extent on the effectiveness of the confinement which, in turn, depends on turbulence. We will exploit our improved understanding of turbulence to seek techniques to further enhance tokamak confinement, and feed these results into the integrated design studies of fusion reactors which are presently under way (e.g. at UKAEA). If we are successful, our research results will underpin designs for more compact fusion reactors, that could form the basis for small modular fusion power plants, or for a fusion components test facility that would be valuable for mitigating risks of delays to commercial fusion power.

The main pathway to realising the above impacts will be through our close links with the UK national fusion programme led by UKAEA. However, we also have close collaborative links with international programmes, including those in Europe, USA, S Korea and China. These collaborative links will provide us with a broader range of fusion facilities and expertise to help meet our research objectives, but also provide a pathway to international impact. We will, of course, seek to strengthen our impact via high profile presentations at international conferences and publications in high impact journals. While our research has a clear application, it is also developing fundamental plasma science. As such, there is potential for publications in journals with a wider readership, beyond the fusion community, such as in Nature and Science.

Turbulence is a feature of most magnetised plasma systems, and plays a role in many solar and astrophysical phenomena, as well as laser-produced plasmas. We have contacted key individuals from astrophysics groups from 14 UK universities, asking about their interest in a 'National Centre for Turbulence in Magnetised Plasmas' that would cut across astrophysical and laboratory plasmas. The fundamental science has similarities across the disciplines, so the different systems, when brought together, could provide opportunities to advance the field beyond the capabilities of the individual communities. We will launch this new Centre as a key part of our Programme, organising three national meetings which, we anticipate, will act as a platform on which to build a strong interdisciplinary plasma turbulence community that will last well beyond the term of the Programme Grant.

A sustainable, cost-effective source of energy is fundamental to our standard of living and the national economy. We will ensure that the general public, industrialists and policy makers are kept aware of the significant contribution that fusion energy can make as a key part of a portfolio of possible solutions through public lectures and targeted events at our universities.

Publications

10 25 50
 
Description Turbulence in tokamak plasmas is driven by small-scale micro-instabilities. We have delivered a number of key findings::
(1) We have developed a novel code, global STELLA, to investigate the rotation generated by gradients of turbulence intensity and other effects; this has now been extended to include electromagnetic fluctuations.
(2) We have studied the pedestal of JET plasmas, finding a role for so-called electron temperature gradient modes using the new STELLA code, and also kinetic ballooning modes using the global ORB5 code.
(3) We have developed "TDoTP baseline equilibria" representative of STEP plasmas and characterised the different classes of micro-instabilities in those, identifying the most dangerous for turbulence; developing non-linear simulations for the electromagnetic turbulence is proving challenging, but we have made significant progress in understanding the challenge and improving convergence.
(4) We have identified extremely extended micro-instability structures in ballooning space (a kind of radial wavenumber space), and developed the analytic theory to simplify their simulation.
(5) We have justified the need for a short wavelength scattering system to resolve turbulence approaching the electron Larmor radius on MAST-U, and prepared an initial specification and design for the instrument. A proposal for funding to build the instrument and implement on MAST-U has been submitted to EPSRC.
(6) We have advanced the SAMI-2 microwave scattering diagnostic, and have installed it ready for experiments on MAST-U.
(7) Commissioning of the beam-emission spectroscopy imaging system on MAST-U is complete, with data taken in the 2022/23 MAST-U campaign.
(8) Simulations of negative triangularity STEP-like plasmas indicate robust kinetic ballooning modes, which may preclude this geometry as a mode of operation for STEP; paper accepted for publication.
(9) We have identified a novel mechanism for turbulence suppression via parallel-to-the-magnetic-field shear in turbulent flows and shown that this mechanism leads to stabilisation of short-wavelength instabilities by long-wavelength turbulence.
(10) We have developed a new model to interpret the data provided by one of the most widely used turbulence diagnostics, Doppler BackScattering (DBS). Previous models ignored effects that become important in spherical tokamaks such as MAST-U.
(11) We have proposed a new explanation for the observation that turbulence is suppressed for sufficiently small gradients, a phenomenon known as the Dimits shift because it can be interpreted as a shift of the gradient at which the plasma becomes unstable.
(12) We have discovered via numerical simulations the coexistence of multiple saturated turbulent states for the same sets of plasma parameters. These states correspond to levels of turbulent transport that differ by an order of magnitude and appear to be needed to match experimental flux levels in the studied JET discharge.
Exploitation Route We work in close collaboration with the STEP team, where the outcomes will be a direct benefit for quantifying the confinement of that demonstration power plant design. We are also working with EUROfusion to find applications in JET and in preparation for ITER and EU DEMO. We also anticipate applications in the stellarator community. We have collaborations with the company Tokamak Energy to support their programme in designing future fusion devices.
Sectors Energy

URL https://www.tdotp.ac.uk
 
Description This programme provides underpinning research for the UKAEA STEP programme, which seeks to build a fusion reactor by 2040 that will deliver net electricity. The effectiveness of the confinement system is a key driver for the design of the STEP fusion reactor, and this confinement is determined by plasma turbulence. Team members have been funded by STEP to provide expertise on the micro-instabilities responsible for driving turbulence in STEP plasmas. We have advised that negative triangularity is not a viable option for STEP because of the enhanced KBM drive. In addition, our expertise in the interaction of microwaves in plasmas has been employed to understand the issues related to current driven by microwaves in STEP. We are also collaborating with the private fusion company Tokamak Energy to inform their future strategy
First Year Of Impact 2020
Sector Energy,Environment
Impact Types Societal,Policy & public services

 
Description EPSRC studentship
Amount £50,000 (GBP)
Funding ID 1991712 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 10/2019 
End 03/2023
 
Description Embedding FAIRness in Plasma Science
Amount £992,754 (GBP)
Funding ID EP/V051822/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 07/2021 
End 06/2026
 
Description Kinetic Theory of MHD Instabilities in Tokamak Plasmas: Neoclassical Tearing Modes and Kink-Ballooning Modes
Amount € 191,368 (EUR)
Funding ID AWP22-ERG-UKAEA/Dudkovskaia 
Organisation EUROfusion 
Sector Public
Country European Union (EU)
Start 03/2022 
End 02/2024
 
Description SAMI-2: two-dimensional Doppler imaging of tokamak plasmas
Amount £213,056 (GBP)
Funding ID EP/S018867/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 01/2019 
End 06/2020
 
Description STEP Plasma Modelling Lot 1: Gyrokinetic plasma modelling for STEP
Amount £443,253 (GBP)
Funding ID CMF/11441 
Organisation UK Atomic Energy Authority 
Sector Public
Country United Kingdom
Start 07/2019 
End 03/2024
 
Description STEP Plasma Modelling Lot 3: Simulations of non-inductive plasma start-up and non-inductive current drive using waves in the electron cyclotron range of frequencies
Amount £253,402 (GBP)
Funding ID CMF/11441 (Lot 3) 
Organisation UK Atomic Energy Authority 
Sector Public
Country United Kingdom
Start 07/2019 
End 03/2021
 
Description STEP Plasma Physics Research
Amount £168,171 (GBP)
Organisation UK Atomic Energy Authority 
Sector Public
Country United Kingdom
Start 04/2022 
End 03/2024
 
Description General Atomics - extensions to TGLF to tight aspect ratio, high beta equilibria 
Organisation General Atomics
Department Fusion Theory and Computational Sciences
Country United States 
Sector Private 
PI Contribution Howard Wilson and Colin Roach of TDoTP, together with a PhD student linked to the programme (Bhavin Patel) have been extending the TGLF tokamak plasma transport model to high beta equilibria characteristic of tight aspect ratio. Bhavin has visited General Atomics (funded from elsewhere) during 2019 to develop the extended model in collaboration.
Collaborator Contribution Expertise of Gary Staebler on TGLF transport model, and also Jeff Candy on Gyro code. Gary Staebler has visited the team in the UK two times, funded from outside the programme)
Impact An abstract was submitted to 2020 EPS plasma physics conference on the work. Bhavin Patel has submitted his PhD thesis, which includes research arising from this collaboration. Bhavin Patel has led a journal publication that included aspects of this research
Start Year 2019
 
Description General Atomics - pedestal modelling for STs 
Organisation General Atomics
Department Fusion Theory and Computational Sciences
Country United States 
Sector Private 
PI Contribution The collaborator funds Dr Koki Imada, employed by University of York, to provide research on spherical tokamak pedestal stability and, later, transport. Prof Howard Wilson provides line management and technical guidance for Dr Imada, as well as joining collaboration meetings with colleagues from General Atomics. The research is relevant to the TDoTP objective to develop a pedestal model valid for the spherical tokamak.
Collaborator Contribution There is some in-kind contribution from their staff, leading the project, but quantifying the contribution to TDoTP is difficult, and the biggest contribution is the cash value. This includes guidance on running the ELITE stability code for real tokamak equilibria, and will include designing experiments on the MAST-U tokamak. The collaborator provides funding to cover the salary costs of Dr Koki Imada. Present funding is due to end in June 2021, but the intention is to continue with a contract extension beyond that,
Impact No outputs yet - too early.
Start Year 2020
 
Description Gyrokinetic analysis of JET ILW shots 
Organisation EUROfusion
Country European Union (EU) 
Sector Public 
PI Contribution We have analysed the stability of JET ILW pedestals. Pedestals are regions of reduced turbulence with large temperature and density gradients. Pedestals are considered to be crucial to achieving fusion energy, but are poorly understood. We have conducted a careful stability analysis of the experimental profiles, discovering a new type of mode. The results will be submitted for publication.
Collaborator Contribution EUROfusion funded 5 researcher-days and provided the data that was used for the gyrokinetic analysis.
Impact The outcome of the collaboration is expected to be reported in a publication
Start Year 2019
 
Description Modelling of the Doppler-Back-Scattering Diagnostic in MAST-Upgrade 
Organisation Culham Centre for Fusion Energy
Country United Kingdom 
Sector Academic/University 
PI Contribution We have developed a model to interpret the data obtained with a Doppler-Back-Scattering (DBS) system. This diagnostic is a non-perturbative method to measure turbulence in tokamaks. Previous analyses of DBS made assumptions about the magnetic field that are not valid for spherical tokamaks such as MAST-Upgrade, and hence needed upgrading.
Collaborator Contribution The Culham Centre for Fusion Energy (CCFE) provided time of one of his researchers, Dr. Jon Hillesheim, and promised run time in MAST-Upgrade in 2020 to check the model. The estimated value of CCFE's contribution does not include the MAST-Upgrade run time yet.
Impact A publication on the model is expected
Start Year 2019
 
Description Partnership on high frequency circuits and components 
Organisation Rutherford Appleton Laboratory
Country United Kingdom 
Sector Academic/University 
PI Contribution Strathclyde have been involved in designing novel high frequency waveguides for many years, and in the measuring of these waveguides in both cold and hot tests. In the context of TDoTP Strathclyde (with York and others) are actively developing a new microwave diagnostic. We believe RAL Space expertise will be important in delivering this diagnostic.
Collaborator Contribution RAL Space's precision machine shop have made critical contributions to the fabrication of the novel waverguides. They have expertise in high performance sources and detectors in the sub-mm wave range that will be vital for future microwave diagnostic developement.
Impact There have been many successful outcomes from this partnership between Strathclyde and RAL Space in the past, including a world leading result in fast wave microwave amplifiers at W band. In the context though of the TDoTP project the outcomes are envisioned for the future.
Start Year 2019
 
Description Research Collaboration on Microwave emission due to kinetic instabilities in an overdense mirror-confined plasma (supported by Royal Society in UK and RFBR in Russia) 
Organisation Russian Academy of Sciences
Department Institute of Applied Physics
Country Russian Federation 
Sector Academic/University 
PI Contribution Strathclyde are contributing the simulations of complex microwave emissions from overdense mirror confined and microwave driven plasma
Collaborator Contribution IAP are undertaking experiments driving plasma in a mirror confined plasma volume using microwaves, studying the micorwave emissions
Impact A joint conference paper has arisen from bilateral visits.
Start Year 2018
 
Title 2D Dynamic Time Warping (DTW) algorithm for python 
Description DTW is a dynamic time warping library designed for tokamak velocimetry measurements. The algorithm works by distorting image 2 into image 1, then using the distorted coordinates to compute the displacements required to distort image 1 into image 2, given at the positions in image 1. 
Type Of Technology Software 
Year Produced 2023 
Open Source License? Yes  
Impact This implementation of an orthogonal dynamic programming approach to optical flow measurement based on the dynamic time warping algorithm is being used by students and researchers at the University of York, as well as UKAEA. Benchmarking and detailed testing and uncertainty quantification of the software on realistic simulated data has recently been reported. It has been found that in optical diagnostics with spatial resolution and size (number of channels) of at least that of the MAST-U beam emission spectroscopy (BES) system, DTW can increase the time-resolution of velocimetry over other methods, often by an order of magnitude, yielding up to sampling frequency speed in the best situations. The method is currently in active use for research on turbulence-scale flow variations including zonal flows and a search for ExB staircases. 
URL https://zenodo.org/record/7649022
 
Title GS2 
Description GS2 is a physics application, developed to study low-frequency turbulence in magnetized plasma. It is typically used to assess the microstability of plasmas produced in the laboratory and to calculate key properties of the turbulence which results from instabilities. It is also used to simulate turbulence in plasmas which occur in nature, such as in astrophysical and magnetospheric systems. 
Type Of Technology Software 
Year Produced 2020 
Open Source License? Yes  
Impact GS2 is a vital tool for a large fraction of the simulations performed under this project, as well as work at labs across the world, including the USA (PPPL, University of Maryland), Japan (NIFS, SOKENDAI), and the UK (CCFE, University of Oxford, University of York). 
URL https://gyrokinetics.gitlab.io/gs2/
 
Title GS2 v8.1.2 
Description GS2 is a physics application, developed to study low-frequency turbulence in magnetized plasma. It is typically used to assess the microstability of plasmas produced in the laboratory and to calculate key properties of the turbulence which results from instabilities. It is also used to simulate turbulence in plasmas which occur in nature, such as in astrophysical and magnetospheric systems. 
Type Of Technology Software 
Year Produced 2022 
Impact This release fixed several important bugs 
URL https://zenodo.org/record/2551066
 
Title stella gyrokinetic code 
Description stella is a piece of software that solves the gyrokinetic-Poisson system of equations to simulate the turbulent dynamics of magnetised plasma inside small 'flux tubes' encompassing magnetic field lines in magnetic confinement fusion plasmas. 
Type Of Technology Software 
Year Produced 2019 
Open Source License? Yes  
Impact This software is still in its infancy but is being used by students and researchers at Oxford, as well as researchers at the national fusion laboratories in the UK (CCFE), Spain (CIEMAT) and Germany (Max Planck Institute for Plasma Physics in Greifswald). This software has facilitated the first simulations of impurity transport in the W7X stellarator and the first nonlinear simulations of a novel turbulence regime in the edge of tokamaks (dominated by so-called 'toroidal ETG' modes). In both stellarators and tokamaks it has been found to reduce the computation time for micro stability calculations and turbulent transport by greater than an order or magnitude for certain physically-relevant parameters. 
URL https://github.com/mabarnes/stella
 
Description Engagement with Strathclyde ANRC/UKAEA senior officer fusion liaison group 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Awareness briefing of senior University Officers and UKAEA staff. These are staff consulted by policy makers and politicians. Awareness of senior UKAEA staff, and awareness of other University practicioners - potential to develop new collaborations
Year(s) Of Engagement Activity 2023