Toroidal Angular Momentum Research Articles

Calculation of tearing mode stability in an inverse aspect-ratio expanded tokamak plasma equilibrium

The tearing mode stability of an inverse aspect-ratio expanded tokamak plasma equilibrium of general shape is investigated using asymptotic matching techniques. Particular emphasis is placed on the conservation of toroidal electromagnetic angular momentum. The TJ code, which is a specific implementation of the results of the investigation, is described.

Full-f gyrokinetic simulations of hydrogen isotope mixing in tokamak plasmas

Hydrogen isotope mixing phenomena in tokamak plasmas are analyzed using global full-f gyrokinetic simulations. Model plasma parameters are chosen based on hydrogen isotope pellet experiments on JET, in which fast hydrogen isotope mixing in the timescale of the energy confinement time τE occurred after injecting deuterium (D) pellets into hydrogen (H) plasmas. Two numerical experiments are conducted using plasma profiles before and after the pellet injection. In both cases, turbulent fluctuations in the plasma core are characterized by ion temperature gradient driven turbulence. In the former case, the density profile of bulk H ions is kept in a quasi-steady state, and the particle confinement time of bulk H ions τH is an order of magnitude longer than τE. In the latter case, the density profiles of bulk H ions and pellet D ions show transient relaxation in the timescale of τE, and both τH and the timescale of particle pinch τD become comparable to τE, indicating the fast hydrogen isotope mixing. In the toroidal angular momentum balance, it is found that the enhanced ion particle transport leading to the fast hydrogen isotope mixing is driven by the toroidal field stress.

Orbit tomography in constants-of-motion phase-space

Tomographic reconstructions of a 3D fast-ion constants-of-motion phase-space distribution function are computed by inverting synthetic signals based on projected velocities of the fast ions along the diagnostic lines of sight. A spectrum of projected velocities is a key element of the spectrum formation in fast-ion D-alpha spectroscopy, collective Thomson scattering, and gamma-ray and neutron emission spectroscopy, and it can hence serve as a proxy for any of these. The fast-ion distribution functions are parameterised by three constants of motion, the kinetic energy, the magnetic moment and the toroidal canonical angular momentum. The reconstructions are computed using both zeroth-order and first-order Tikhonov regularisation expressed in terms of Bayesian inference to allow uncertainty quantification. In addition to this, a discontinuity appears to be present in the solution across the trapped-passing boundary surface in the three-dimensional phase space due to a singularity in the Jacobian of the transformation from position and velocity space to phase space. A method to allow for this apparent discontinuity while simultaneously penalising large gradients in the solution is demonstrated. Finally, we use our new methods to optimise the diagnostic performance of a set of six fans of sightlines by finding where the detectors contribute most complementary diagnostic information for the future COMPASS-Upgrade tokamak.

Measurement of toroidal Alfvén eigenmode-driven fast-ion flows using an imaging neutral particle analyzer at ASDEX Upgrade

Alfvén Eigenmodes-driven fast-ion flows have been measured for the first time at the ASDEX Upgrade tokamak using an imaging neutral particle analyzer. The flow is such that particles are expelled from the mode location, losing energy as they move outwards. This flow aligns well with the projection of the lines of constant magnetic moment and constant E′ ( E′=E−ωPϕ/n , with E being the particle energy, ω,n the mode frequency and toroidal number and P φ the toroidal canonical angular momentum). The observed redistributions are consistent with full-orbit simulations performed with the ASCOT5 code, using the mode structure predicted by the linear gyro-kinetic code LIGKA.

Diagnostic weight functions in constants-of-motion phase-space

The fast-ion phase-space distribution function in axisymmetric tokamak plasmas is completely described by the three constants of motion: energy, magnetic moment and toroidal canonical angular momentum. In this work, the observable regions of constants-of-motion phase-space, given a diagnostic setup, are identified and explained using projected velocities of the fast ions along the diagnostic lines-of-sight as a proxy for several fast-ion diagnostics, such as fast-ion Dα spectroscopy, collective Thomson scattering, neutron emission spectroscopy and gamma-ray spectroscopy. The observable region in constants-of-motion space is given by a position condition and a velocity condition, and the diagnostic sensitivity is given by a gyro-orbit and a drift-orbit weighting. As a practical example, 3D orbit weight functions quantifying the diagnostic sensitivity to each point in phase-space are computed and investigated for the future COMPASS-Upgrade and MAST-Upgrade tokamaks.

Impact of magnetic ripple on neoclassical equilibrium in gyrokinetic simulations

The effect of magnetic field ripple on tokamak plasma without turbulence is studied numerically and augmented with a reduced analytical model that includes neoclassical processes in the presence of non-axisymmetric perturbation and stochastic transport. For this study, a magnetic field ripple perturbation has been implemented in the GYSELA gyrokinetic code. This implementation has been verified thanks to a test of toroidal angular momentum conservation. The GYSELA code was then successfully benchmarked against the NEO code, which solves the drift kinetic equation, and against the reduced model in the collisionality range [0.05–0.5] for several amplitudes of the magnetic ripple. An observation, shared by the model, the NEO code and GYSELA simulations is that the thermal drive of the mean poloidal velocity—measured by the coefficient—decreases sharply for large yet experimentally relevant magnetic ripple amplitudes, and may even change sign.

Metriplectic foundations of gyrokinetic Vlasov–Maxwell–Landau theory

This Letter reports on a metriplectic formulation of a collisional, nonlinear full-f electromagnetic gyrokinetic theory compliant with energy conservation and monotonic entropy production. In an axisymmetric background magnetic field, the toroidal angular momentum is also conserved. Notably, a new collisional current, contributing to the gyrokinetic Maxwell–Ampère equation and the gyrokinetic charge conservation law, is discovered.

Simulation study of fast ion losses associated with the rotating n = 1 resonant magnetic perturbations in KSTAR

This paper describes a simulation framework for testing the fast ion loss mechanism associated with the experimentally observed beam ion losses when an externally applied toroidally rotating perturbed magnetic field is used to control edge localized modes in the KSTAR tokamak. The simulations reproduce the key qualitative features of neutral beam injection (NBI) ion detection by a fast ion loss detector. The NBI ion losses in the simulation mainly occur for passing particles due to orbit stochastization, which is caused by orbit resonance with the 3D field perturbations. The relative toroidal angle of NBI ion deposition to the orbit island determines the radial path of the particles between confinement or loss. The fast ion loss quantity depends on the toroidal angle rotation of the 3D field with respect to the beam deposition position. The averaged transport of toroidal canonical angular momentum reveals that linear and nonlinear resonance of the NBI ions within the 3D field are the main factors determining fast ion transport and loss.

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E × B drift to cause a nondiffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should, in principle, relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear ŝ was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

Development of a gyrokinetic hyperbolic solver based on discontinuous Galerkin method in tokamak geometry

A hyperbolic solver is developed for the gyrokinetic equation in tokamak geometry. Aiming a whole device modeling of fusion plasma surrounded by first walls of tokamak device, an unstructured spatial mesh is introduced. The discontinuous Galerkin (DG) method is used to discretize the gyrokinetic equation on the mesh and test various numerical elements for the discretization. Based on the conservations of physical quantities such as mass, kinetic energy, and toroidal canonical angular momentum in an axisymmetric configuration of toroidal plasma, we investigate the effects of basis functions for the DG method on the numerical solutions. With proper choices of the basis functions and spatial grid resolutions, the conservations of the key physical quantities are shown to be satisfied nearly to machine accuracies in the simplified circular magnetic geometry. Even in realistic tokamak geometry with machine wall boundaries, it is shown that a good conservation property can be demonstrated if the flux across boundaries of a test domain is carefully accounted for. Also, the invariance of the canonical Maxwellian distribution function in time is well satisfied with the developed solver. The effect of weighting functions for the basis is investigated too. Overall, the Maxwellian weighted basis shows a similar conservation property with the polynomial basis. On the other hand, the Maxwellian weighted basis shows better performance in resolving small scale structures in velocity space, which can be utilized to set up an efficient basis set to simulate fine structures with less computational costs. The parallelization of the newly developed solver is also reported. Employing MPI for the parallelization, the solver shows good performances up to a few thousand CPU cores.

Toroidal modeling of plasma response to RMP fields for HL-2M

The plasma response to the n = 1, 2, 4 (n is the toroidal mode number) resonant magnetic perturbation (RMP) fields, and the consequences on the fast ion confinement, are numerically investigated for a reference high-pressure plasma in HL-2M, by utilizing the linear resistive magnetohydrodynamic code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681). The best coil current configurations, in terms of the coil phasing between the upper and lower rows of coils for controlling type-I edge localized modes (ELMs) in HL-2M, are identified as −130, −30, 180 degrees for the n = 1, 2, 4 fields, respectively, based on the edge peeling-tearing plasma response criterion. The plasma is found to substantially amplify the applied vacuum RMP field with the best coil phasing for the reference HL-2M equilibrium. The overall field amplification factor, defined as the peak-to-peak ratio of the poloidal spectra for the total field perturbation including the plasma response and the vacuum field alone, is about five for all n’s. The amplification, however, does not occur with the worst coil phasing for ELM control. This field amplification due to the high-pressure plasma response, together with the plasma screening of the resonant radial field components in the core region, have several consequences on the fast ion confinement in HL-2M during ELM control with RMP. (i) Three-dimensional fields including the plasma response, and with the best coil phasing, substantially enhance the distortion of fast ion orbits compared to the vacuum field approximation. With the n = 1 RMP, the plasma-response-induced enhancement of the orbit distortion reaches a factor of four when measured in terms of the canonical toroidal angular momentum. (ii) With the best coil phasing, the plasma response widens the stochastic region for the particle orbits on the Poincaré plane. (iii) The orbit islands, including the plasma response, remain as large as the vacuum counterparts in the plasma core where strong screening of the resonant field components occur. All these effects lead to enhanced fast ion transport (and loss) in the high-pressure HL-2M plasma, when the best RMP spectrum is applied to control ELMs.

Orbit-modulated transport and sources in time-dependent plasmas

The orbit-modulation framework for steady-state ion-orbit boundary fluxes is reformulated to apply to time-dependent plasmas. The new formulation can quantify the extent to which orbit fluxes are affected by the time dependence of the distribution function, and may be applied both to time-dependent physical scenarios and to computationally expensive simulations that are unable to reach a true steady state. As the system tends to a time-independent state, the orbit-loss fluxes approach their steady-state values, which are determined entirely by upstream transport (turbulent and collisional) and upstream sources and sinks. This result holds for the fluxes of any moments that are expressible in terms of the invariants that determine the orbits, including parallel toroidal angular momentum.

Drift orbit islands of energetic particles due to 3D fields in ITER

In this paper, we report a comparative study of the orbit islands formed by the guiding center drift motion of test energetic particles (EPs), including both high-energy deuterium ions and fusion-born α-particles, with respect to the magnetic islands formed by tracing the field lines, for an ITER plasma representing the 15 MA baseline scenario. The particle drift orbit is modified by the presence of static resonant magnetic perturbations (RMPs) of different toroidal mode numbers n (n = 1 and 3 in this study), forming orbit islands in the Poincaré plane by passing EPs. The key findings are as follows. (i) With an n = 1 RMP field, the size of the orbit islands combined with the total RMP field, including the plasma response, is about three times smaller than that obtained by assuming the corresponding vacuum field. (ii) The orbit island size is not sensitive to the EP energy, and is comparable to those of the magnetic islands. (iii) Passing EPs with outward radial orbital drift form orbit islands that shift inward with respect to the corresponding magnetic islands, and vice versa. This radial shift, ΔΨ, measured in the normalized equilibrium poloidal flux, is quantified as ΔΨ ≃ 0.03X for the ITER baseline plasma and assuming n = 1 RMPs, with , where M is the particle mass, M p is the proton mass, E is the particle energy, and Z is the particle charge number. (iv) Trapped EPs do not form drift orbit islands, even in the presence of the 3D fields produced by the 90 kAt ELM control coil current in ITER, and thus possess good confinement properties. Nevertheless, the deformation of the trapped particle drift orbits in the Poincaré plane shows that the canonical toroidal angular momentum of EPs is no longer conserved in the presence of the RMP fields.

Orbit-modulated plasma transport and sources

Using a general conservative gyrokinetic framework in axisymmetric magnetic geometry, the steady-state orbit-loss boundary fluxes are demonstrated to be predominantly determined by upstream transport (turbulent and collisional) and upstream sources and sinks. Within the framework, the details of the equilibrium orbits, even accounting for large but axisymmetric and time-independent electric fields and poloidal pressure gradients, are demonstrated to contribute only by modulating the upstream transport and source terms. This is true not only for gyrocenter fluxes, but also for any other moments that are expressible in terms of the adiabatic invariants that determine the orbits, such as parallel toroidal angular momentum. An explicit reformulation of the orbit-loss terms facilitates their evaluation via numerical diagnostics or reduced models. Additional contributions for realistic cases, including those of fast ions and time-evolving equilibria, are briefly discussed.

Intrinsic rotation in axisymmetric devices

Toroidal rotation is critical for fusion in tokamaks, since it stabilizes instabilities that can otherwise cause disruptions or degrade confinement. Unlike present-day devices, ITER might not have enough neutral-beam torque to easily avoid these instabilities. We must therefore understand how the plasma rotates ‘intrinsically,’ that is, without applied torque. Experimentally, torque-free plasmas indeed rotate, with profiles that are often non-flat and even non-monotonic. The rotation depends on many plasma parameters including collisionality and plasma current, and exhibits sudden bifurcations (‘rotation reversals’) at critical parameter values.Since toroidal angular momentum is conserved in axisymmetric systems, and since experimentally inferred momentum transport is much too large to be neoclassical, theoretical work has focused on rotation drive by nondiffusive turbulent momentum fluxes. In the edge, intrinsic rotation relaxes to a steady state in which the total momentum outflux from the plasma vanishes. Ion drift orbits, scrape-off-layer flows, separatrix geometry, and turbulence intensity gradient all play a role. In the core, nondiffusive and viscous momentum fluxes balance to set the rotation gradient at each flux surface. Although many mechanisms have been proposed for the nondiffusive fluxes, most are treated in one of two distinct but related gyrokinetic formulations. In a radially local fluxtube, appropriate for , the lowest-order gyrokinetic formulations exhibit a symmetry that prohibits nondiffusive momentum flux for nonrotating plasmas in an up-down symmetric magnetic geometry with no shear. Many symmetry-breaking mechanisms have been identified, but none have yet been conclusively demonstrated to drive a strong enough flux to explain commonly observed experimental rotation profiles. Radially global gyrokinetic simulations naturally include many symmetry-breaking mechanisms, and have shown cases with experimentally relevant levels of nondiffusive flux. These promising early results motivate further work to analyze, verify, and validate.This article provides a pedagogical introduction to intrinsic rotation in axisymmetric devices. Intended for both newcomers to the topic and experienced practitioners, the article reviews a broad range of topics including experimental and theoretical results for both edge and core rotation, while maintaining a focus on the underlying concepts.

Numerical study of neutral beam injection prompt loss at KSTAR poloidal limiters

We use a newly developed simulation tool to numerically study fast-ion loading on plasma-facing components (PFCs) at the Korea Superconducting Tokamak Advanced Research facility in the high-poloidal-β plasma operation regime. The new code can calculate neutral beam ionization and follow the guiding center orbit of ionized particles. The results of the simulation indicate that fast ions ionized in the high-field side drift out and strike the PFCs as they rotate poloidally. Momentum projection onto a phase space defined by canonical toroidal angular momentum and magnetic moment leads to a simple criterion to avoid fast-ion loading on poloidal limiters (PLs), which are PFCs that are toroidally localized on the low-field side. Control of fast-ion loss is examined by varying the plasma current and plasma boundary. A larger plasma current and inner shifting of the outer plasma boundary is predicted to substantially reduce fast-ion loading on the PLs. The proposed phase-space criterion is qualitatively consistent with the results of the simulation of fast-ion control.

Fully Kinetic Simulation of Ion-Temperature-Gradient Instabilities in Tokamaks

The feasibility of using full ion kinetics, instead of gyrokinetics, in simulating low-frequency Ion-Temperature-Gradient (ITG) instabilities in tokamaks has recently been demonstrated. The present work extends the full ion kinetics to the nonlinear regime and investigates the nonlinear saturation of a single-n ITG instability due to the E × B trapping mechanism (n is the toroidal mode number). The saturation amplitude predicted by the E × B trapping theory is found to agree with the saturation level observed in the simulation. In extending to the nonlinear regime, we developed a toroidal Boris full orbit integrator, which proved to be accurate in capturing both the short-time scale cyclotron motion and long time scale drift motion, with good kinetic energy conservation and toroidal angular momentum conservation in tokamak equilibrium magnetic fields. This work also extends the previous work from analytic circular magnetic equilibria to general numerical magnetic equilibria, enabling simulation of realistic equilibria reconstructed from tokamak experiments.

Non-isomorphic radial wavenumber dependencies of residual zonal flows in ion and electron Larmor radius scales, and effects of initial parallel flow and electromagnetic potentials in a circular tokamak

Radial wavenumber dependencies of the residual zonal potential for E × B flow in a circular, large aspect ratio tokamak is investigated by means of the collisionless gyrokinetic simulations of Rosenbluth–Hinton (RH) test and the semi-analytic approach using an analytic solution of the gyrokinetic equation Rosenbluth and Hinton (1998 Phys. Rev. Lett. 80 724). By increasing the radial wavenumber from an ion Larmor radius scale to an electron Larmor radius scale , the well-known level is retained, while the level remains O(1) when the wavenumber is decreased from the electron to the ion Larmor radius scale, if physically same adiabatic assumption is presumed for species other than the main species that is treated kinetically. The conclusion is not modified by treating both species kinetically, so that in the intermediate scale between the ion and electron Larmor radius scale it seems difficult to determine the level uniquely. The toroidal momentum conservation property in the RH test is also investigated by including an initial parallel flow in addition to the perpendicular flow. It is shown that by taking a balance between the initial parallel flow and perpendicular flows which include both E × B flow and diamagnetic flow in the initial condition, the mechanical toroidal angular momentum is approximately conserved despite the toroidal symmetry breaking due to the finite radial wavenumber zonal modes. Effect of electromagnetic potentials is also investigated. When the electromagnetic potentials are applied initially, fast oscillations which are faster than the geodesic acoustic modes are introduced in the decay phase of the zonal modes. Although the residual level in the long time limit is not modified, this can make the time required to reach the stationary zonal flows longer and may weaken the effectiveness of the turbulent transport suppression by the zonal flows.

Trapped particle precession and sub-bounce zonal flow dynamics in tokamaks

A drift-kinetic calculation in an axisymmetric tokamak, with super-diamagnetic flows, is presented to elucidate the relation between the radial electric field, $E_{r}$, zonal flows and the rapid precession of the trapped particle (TP) population. It has been shown earlier (Rosenbluth & Hinton, Phys. Rev. Lett., vol. 80(4), 1998, p. 724, hereafter RH) that an initial radial electric field results in geodesic acoustic mode oscillations which subsequently Landau damp, resulting in a much smaller final residual electric field, and accompanying parallel zonal flows. We observe an apparent paradox: the final angular momentum in the RH parallel zonal flow is much smaller than the angular momentum expected from the well-known rapid precession of the trapped particle population in the RH residual electric field. We reconcile this paradox by illuminating the presence of a population of reverse circulating particle flows that, dominantly, are equal and opposite to the rapid TP precession. Mathematically, the calculation is facilitated by transforming to an energy coordinate shifted from conventional by an amount proportional to $E_{r}$. We also discuss the well-known RH coefficient in the context of effective mass and show how the TP precession and the opposite circulating flows contribute to this mass. Finally, we show that in the long wavelength limit, the RH flows and RH coefficient arise as a natural consequence of conservation of toroidal angular momentum and the second adiabatic invariant.

Experimental challenges to stiffness as a transport paradigm

Two power scans were carried out in H-mode plasmas in DIII-D; one employed standard co-current neutral beam injection (NBI), while the other used a mixture of co-current and counter-current NBI to scan power while holding the torque to a low fixed value. Analysis of the ion and electron heat transport, ion toroidal angular momentum transport, and thermal deuterium transport from these scans is presented. Invariance of the gradients or gradient scalelengths, as might be expected from stiff transport, was not generally observed. When invariance was seen, it was not accompanied by a strong increase in transport, except in the case of the absolute deuterium ion transport. Conduction in the ion channel is the dominant energy loss mechanism. The variation of the ion heat transport with applied power is similar for the co-injection and fixed torque scans, indicating that E × B shearing is not determining the plasma response to additional power. There is however, a quantitative difference in the transport between the two scans, indicating E × B shearing does play a role in the transport. Comparison of these results with a previous experiment that directly probed stiffness at a single radius leads to the following conclusion: while local stiffness as formally defined may hold, invariance of the gradients or normalized scalelengths does not follow from stiff transport in more practical scaling experiments, such as the power scans discussed here. Possible reasons for the lack of correspondence between the local picture and the global expectations are discussed.