An efficient method is given to reconstruct the current profile parameters, the plasma shape, and a current profile consistent with the magnetohydrodynamic equilibrium constraint from external magnetic measurements, based on a Picard iteration approach which approximately conserves the measurements. Computational efforts are reduced by parametrizing the current profile linearly in terms of a number of physical parameters. Results of detailed comparative calculations and a sensitivity study are described. Illustrative calculations to reconstruct the current profiles and plasma shapes in ohmically and auxiliarily heated Doublet III plasmas are given which show many interesting features of the current profiles.

https://iopscience.iop.org/article/10.1088/0029-5515/25/11/007/pdf

Reconstruction of current profile parameters and plasma shapes in tokamaks

One of the scientific success stories of fusion research over the past decade is the development of the ExB shear stabilization model to explain the formation of transport barriers in magnetic confinement devices. This model was originally developed to explain the transport barrier formed at the plasma edge in tokamaks after the L (low) to H (high) transition. This concept has the universality needed to explain the edge transport barriers seen in limiter and divertor tokamaks, stellarators, and mirror machines. More recently, this model has been applied to explain the further confinement improvement from H (high)-mode to VH (very high)-mode seen in some tokamaks, where the edge transport barrier becomes wider. Most recently, this paradigm has been applied to the core transport barriers formed in plasmas with negative or low magnetic shear in the plasma core. These examples of confinement improvement are of considerable physical interest; it is not often that a system self-organizes to a higher energy state with reduced turbulence and transport when an additional source of free energy is applied to it. The transport decrease that is associated with ExB velocity shear effects also has significant practical consequences for fusion research. The fundamental physics involved in transport reduction is the effect of ExB shear on the growth, radial extent and phase correlation of turbulent eddies in the plasma. The same fundamental transport reduction process can be operational in various portions of the plasma because there are a number ways to change the radial electric field Er. An important theme in this area is the synergistic effect of ExB velocity shear and magnetic shear. Although the ExB velocity shear appears to have an effect on broader classes of microturbulence, magnetic shear can mitigate some potentially harmful effects of ExB velocity shear and facilitate turbulence stabilization.

/pdf/effects-of-e-b-velocity-shear-and-magnetic-shear-on-5f3wdz73f7.pdf

Effects of E×B velocity shear and magnetic shear on turbulence and transport in magnetic confinement devices

The major increase in discharge duration and plasma energy in a next step DT fusion reactor will give rise to important plasma-material effects that will critically influence its operation, safety and performance. Erosion will increase to a scale of several centimetres from being barely measurable at a micron scale in today's tokamaks. Tritium co-deposited with carbon will strongly affect the operation of machines with carbon plasma facing components. Controlling plasma-wall interactions is critical to achieving high performance in present day tokamaks, and this is likely to continue to be the case in the approach to practical fusion reactors. Recognition of the important consequences of these phenomena stimulated an internationally co-ordinated effort in the field of plasma-surface interactions supporting the Engineering Design Activities of the International Thermonuclear Experimental Reactor project (ITER), and significant progress has been made in better understanding these issues. The paper reviews the underlying physical processes and the existing experimental database of plasma-material interactions both in tokamaks and laboratory simulation facilities for conditions of direct relevance to next step fusion reactors. Two main topical groups of interaction are considered: (i) erosion/redeposition from plasma sputtering and disruptions, including dust and flake generation and (ii) tritium retention and removal. The use of modelling tools to interpret the experimental results and make projections for conditions expected in future devices is explained. Outstanding technical issues and specific recommendations on potential R&D avenues for their resolution are presented.

/pdf/plasma-material-interactions-in-current-tokamaks-and-their-1a6iwv89x6.pdf

Plasma{material interactions in current tokamaks and their implications for next step fusion reactors

Progress in the area of MHD stability and disruptions, since the publication of the 1999 ITER Physics Basis document (1999 Nucl. Fusion 39 2137-2664), is reviewed. Recent theoretical and experimental research has made important advances in both understanding and control of MHD stability in tokamak plasmas. Sawteeth are anticipated in the ITER baseline ELMy H-mode scenario, but the tools exist to avoid or control them through localized current drive or fast ion generation. Active control of other MHD instabilities will most likely be also required in ITER. Extrapolation from existing experiments indicates that stabilization of neoclassical tearing modes by highly localized feedback-controlled current drive should be possible in ITER. Resistive wall modes are a key issue for advanced scenarios, but again, existing experiments indicate that these modes can be stabilized by a combination of plasma rotation and direct feedback control with non-axisymmetric coils. Reduction of error fields is a requirement for avoiding non-rotating magnetic island formation and for maintaining plasma rotation to help stabilize resistive wall modes. Recent experiments have shown the feasibility of reducing error fields to an acceptable level by means of non-axisymmetric coils, possibly controlled by feedback. The MHD stability limits associated with advanced scenarios are becoming well understood theoretically, and can be extended by tailoring of the pressure and current density profiles as well as by other techniques mentioned here. There have been significant advances also in the control of disruptions, most notably by injection of massive quantities of gas, leading to reduced halo current fractions and a larger fraction of the total thermal and magnetic energy dissipated by radiation. These advances in disruption control are supported by the development of means to predict impending disruption, most notably using neural networks. In addition to these advances in means to control or ameliorate the consequences of MHD instabilities, there has been significant progress in improving physics understanding and modelling. This progress has been in areas including the mechanisms governing NTM growth and seeding, in understanding the damping controlling RWM stability and in modelling RWM feedback schemes. For disruptions there has been continued progress on the instability mechanisms that underlie various classes of disruption, on the detailed modelling of halo currents and forces and in refining predictions of quench rates and disruption power loads. Overall the studies reviewed in this chapter demonstrate that MHD instabilities can be controlled, avoided or ameliorated to the extent that they should not compromise ITER operation, though they will necessarily impose a range of constraints.

https://iopscience.iop.org/article/10.1088/0029-5515/47/6/S03/pdf

Chapter 3: MHD stability, operational limits and disruptions

Nonlinear gyrokinetic equations play a fundamental role in our understanding of the long-time behavior of strongly magnetized plasmas. The foundations of modern nonlinear gyrokinetic the- ory are based on three important pillars: (1) a gyrokinetic Vlasov equation written in terms of a gyrocenter Hamiltonian with quadratic low-frequency ponderomotive-like terms; (2) a set of gyrokinetic Maxwell (Poisson-Ampere) equations written in terms of the gyrocenter Vlasov dis- tribution that contain low-frequency polarization (Poisson) and magnetization (Ampere) terms derived from the quadratic nonlinearities in the gyrocenter Hamiltonian; and (3) an exact energy conservationlaw for the gyrokineticVlasov-Maxwell equations that includes all the relevant linear and nonlinear coupling terms. The foundations of nonlinear gyrokinetic theory are reviewed with an emphasis on the rigorous applications of Lagrangian and Hamiltonian Lie-transform perturba- tion methods used in the variationalderivationof nonlineargyrokineticVlasov-Maxwell equations. The physical motivations and applications of the nonlinear gyrokinetic equations, which describe the turbulent evolution of low-frequency electromagnetic fluctuations in a nonuniform magnetized plasmas with arbitrary magnetic geometry, are also discussed.

/pdf/foundations-of-nonlinear-gyrokinetic-theory-3l78lnlpai.pdf

Foundations of nonlinear gyrokinetic theory

An extensive numerical simulation of the Ignitor-Ult plasma dynamics has been carried out using both the free boundary equilibrium code TEQ and the transport and MHD free boundary code TSC The stability behavior of the plasma - vacuum vessel - poloidal field coils - passive plate system has been analyzed for different machine and plasma configurations The response of several feedback systems has been studied taking into account the real characteristics of the power supply The converter and its controller behavior has been extensively studied by means of the simulation program CSMP that uses a lumped parameter model of the system The results show a good stability behavior over a wide operating range

Vertical stability analysis for the ignitor configuration

The ARIES reactor study group has found an economically attractive ST-based reactor configuration with: A = 1.6, {kappa} = 3.4, {delta} = 0.65, {beta} = 50%, {beta}{sub N} = 7.3, f{sub BS} = 0.95, R{sub 0} = 3.2 meters, B{sub t0} = 2.08 Tesla, and I{sub P} = 28.5 MA which yields a cost of electricity of approximately 80mils/kWh. MHD stability analysis finds that a broad pressure profile is optimal for wall-stabilizing the pressure driven kink modes typical of such configurations, and that wall stabilization is crucial to achieving the high {beta} needed for an economical power plant. The 6MW high-harmonic fast wave system presently being installed on NSTX should allow real-time control of the plasma {beta}, and in combination with NBI may permit experimental investigations of the effect of pressure profile peaking on MHD stability in the near-term. In the longer term, ejection of ions through resonant interaction with HHFW might be used to induce a controllable edge radial electric field with potentially interesting effects on edge MHD and confinement.

/pdf/ideal-mhd-stability-characteristics-of-advanced-operating-4mgn8jubyl.pdf

Ideal MHD Stability Characteristics of Advanced Operating Regimes in Spherical Torus Plasmas and the Role of High Harmonic Fast Waves

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. We demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. We generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in [Formula: see text] tokamaks for certain q-profiles and β values.

Ideal MHD induced temperature flattening in spherical tokamaks

Alfv\'enic modes in the current quench (CQ) stage of the tokamak disruption have been observed in experiments. In DIII-D the excitation of these modes is associated with the presence of high-energy runaway electrons, and a strong mode excitation is often associated with the failure of RE plateau formation. In this work we present results of self-consistent kinetic-MHD simulations of RE-driven compressional Alfv\'en eigenmodes (CAEs) in DIII-D disruption scenarios, providing an explanation of the CQ modes. Simulation results reveal that high energy trapped REs can have resonance with the Alfv\'en mode through their precession motion, and the resonance frequency is proportional to the energy of REs. The mode frequencies and their relationship with the RE energy are consistent with experimental observation. The perturbed magnetic fields from the modes can lead to spatial diffusion of runaway electrons including the nonresonant passing ones, thus providing the theoretical basis for a potential approach for runaway electron mitigation.

Stephen Jardin

Papers

Vertical stability analysis for the ignitor configuration

Ideal MHD Stability Characteristics of Advanced Operating Regimes in Spherical Torus Plasmas and the Role of High Harmonic Fast Waves

Ideal MHD induced temperature flattening in spherical tokamaks

Self-consistent simulation of compressional Alfv\'en eigenmodes excited by runaway electrons

Simulations of Lithium Lorentz Force Accelerator (LiLFA) Plasma Flows