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

A critical issue for ITER is to evaluate the forces produced on the surrounding conducting structures during plasma disruptions. We calculate the non-axisymmetric ‘sideways’ wall force Fx, produced in disruptions. Simulations were carried out of disruptions produced by destabilization of n = 1 modes by a vertical displacement event (VDE). The force depends strongly on γτwall, where γ is the mode growth rate and τwall is the wall penetration time, and is largest for γτwall = constant, which depends on initial conditions. Simulations of disruptions caused by a model of massive gas injection were also performed. It was found that the wall force increases approximately offset linearly with the displacement from the magnetic axis produced by a VDE. These results are also obtained with an analytical model. Disruptions are accompanied by toroidal variation of the plasma current Iφ. This is caused by toroidal variation of the halo current, as verified computationally and analytically.

http://iopscience.iop.org/article/10.1088/0029-5515/53/7/073018/pdf

Sideways wall force produced during tokamak disruptions

https://fire.pppl.gov/ARIES_AT_Equilibrium_kessel.pdf

Plasma profile and shape optimization for the advanced tokamak power plant, ARIES-AT

In the absence of a conducting wall, typical spheromak plasmas are unstable to tilting and/or shifting modes. The effects of the cross-sectional shape, aspect ratio, and the location of a conducting wall on the stability of these modes are investigated. A circular cross-section (b/a ~ 1) configuration with a flux hole δ = 0.5 (aspect ratio R/a = 2) will will be stabilized by an ellipsoidal wall of mean separation of 1.3 minor radii from the plasma.

/pdf/tilting-and-shifting-modes-in-a-spheromak-c8zobbsxfn.pdf

Tilting and shifting modes in a spheromak

Calculations of axisymmetric stability of Tokamak plasmas with active and passive feedback

The multilevel physics, massively parallel plasma simulation code, M3D has been used to study spherical toris (STs) and tokamaks. The magnitude of outboard shift of density profiles relative to electron temperature profiles seen in NSTX under strong toroidal flow is explained. Internal reconnection events in ST discharges can be classified depending on the crash mechanism, just as in tokamak discharges; a sawtooth crash, disruption due to stochasticity, or high-β disruption. Toroidal shear flow can reduce linear growth of internal kink. It has a strong stabilizing effect nonlinearly and causes mode saturation if its profile is maintained, e.g. through a fast momentum source. Normally, however, the flow profile itself flattens during the reconnection process, allowing a complete reconnection to occur. In some cases, the maximum density and pressure spontaneously occur inside the island and cause mode saturation. Gyrokinetic hot particle/MHD hybrid studies of NSTX show the effects of fluid compression on a fast-ion driven n = 1 mode. MHD studies of recent tokamak experiments with a central current hole indicate that the current clamping is due to sawtooth-like crashes, but with n = 0.

/pdf/nonlinear-simulation-studies-of-tokamaks-and-sts-2ar8yiaqq1.pdf

Stephen Jardin

Papers

Sideways wall force produced during tokamak disruptions

Plasma profile and shape optimization for the advanced tokamak power plant, ARIES-AT

Tilting and shifting modes in a spheromak

Calculations of axisymmetric stability of Tokamak plasmas with active and passive feedback

Nonlinear Simulation Studies of Tokamaks and STs