Search

We develop a simple model for the generation and amplification of intrinsic axial flow in a linear device, Controlled Shear Decorrelation Experiment (CSDX). This model develops a novel dynamical symmetry breaking mechanism in drift wave turbulence, which does not require complex magnetic field structure, such as shear. Thus, the model is applicable to both tokamaks and linear devices. This mechanism is, essentially, a form of negative viscosity phenomenon. Negative compressibility ITG turbulence can also induce a negative viscosity increment. However, we show that no intrinsic axial flow can be generated by pure ITG turbulence in a straight magnetic field. When the flow gradient is steepened by any drive mechanism, the flow profile saturates at a level close to the value above which parallel shear flow instability (PSFI) becomes dominant over the ITG instability. This saturated flow gradient exceeds the PSFI linear threshold, and grows with $\nabla T_{i0}$ as $|\nabla V_\parallel| / |k_\parallel c_s| \sim|\nabla T_{i0}|^{2/3} / (k_\parallel T_{i0})^{2/3}$.The coupling of azimuthal and axial flows in CSDX--in absence of magnetic shear--is investigated.In particular, we focus on the apportionment of turbulence energy between azimuthal and axial flows, and how the azimuthal flow shear affects axial flow generation and saturation by drift wave turbulence.Detailed measurements of intrinsic axial flow parallel to the magnetic field are performed on CSDX, with no axial momentum input.The results present a direct demonstration that the broken spectral symmetry of drift wave turbulence causes the development of axial mean flows in cylindrical magnetized plasmas.Measurements suggest the axial flow is parasitic to the drift wave--zonal flow system.Besides, we show that consideration of wave--flow resonance resolves the long-standing problem of how zonal flows (ZFs) saturate in the limit of weak or zero frictional drag and also determines the ZF scale directly from analysis. We show that resonant vorticity mixing, which conserves potential enstrophy, enables ZF saturation in the absence of drag, and so is effective at regulating the Dimits up-shift regime. Vorticity mixing is incorporated as a nonlinear, self-regulation effect in an extended 0D predator--prey model of drift--ZF turbulence.

This paper studies the effect of 3D magnetic perturbations (MPs) on fast-ion confinement, and its impact on the toroidal rotation velocity profile. Two low collisionality H-mode experiments carried out at the ASDEX Upgrade tokamak have been analysed. The two discharges feature different magnetic field helicity (q95), and differences in the velocity-space and level of fast-ion losses are observed. A new analysis technique has been developed that sheds light on the dependencies between fast-ion losses and toroidal rotation, providing for the first time correlation patterns resolved in radius and velocity space of the lost fast-ions. The correlation intensifies towards the plasma edge and is strongly dependent on the orbit topology of the lost fast-ions. The ASCOT orbit following code has been used to characterize the fast-ion resonant transport and beam driven torques, using the vacuum approach and including plasma response (PR). The change of the toroidal canonical momentum, which serves as figure of merit for resonant fast-ion transport, has been calculated with ASCOT. The beam geometry and q95 are found to have a strong impact on the fast-ion transport and losses. The fast-ion transport induced by the MPs affects the beam driven torques. The effect of the changes of the j × B and collisional torques on plasma rotation is analysed using the torques simulated by ASCOT and simple momentum balance calculations. For the low q95 = 3.8 discharge, which benefits from a resonant amplification, we find excellent agreement with the measured variation of the toroidal velocity. For the high q95 = 5.5 discharge, the inclusion of the PR improves the comparison with experimental data with respect to the vacuum estimation, but still some differences with experiments are observed. This suggests that other non-resonant effects could play a role for the determination of the toroidal rotation profile.

The prediction of plasma rotation is of high interest for fusion research due to the effects of the rotation upon magnetohydrodynamic (MHD) instabilities, impurities, and turbulent transport in general. In this work, an analysis method was studied and validated to reliably extract momentum transport coefficients from neutral beam injection (NBI) modulation experiments. To this end, a set of discharges was created with similar background profiles for the ion and electron temperatures, the heat fluxes, the electron density, and the plasma rotation that, therefore, should exhibit similar momentum transport coefficients. In these discharges, a range of temporal perturbations were imposed by modulating and varying the power deposition of the NBI, electron-cyclotron-resonance heating (ECRH), and ion-cyclotron-resonance heating (ICRH). The transport model including diffusion, convection, and residual stress was implemented within the ASTRA code. The Prandtl number P r = χ φ / χ i was assessed via the GKW code. A convective Coriolis pinch was fitted and the intrinsic torque from the residual stress was estimated. The obtained transport coefficients agree within error bars for sufficiently small imposed temperature perturbations, as would be expected, from the similar background profiles. This successful validation of the methodology opens the door to study the parametric dependence of the diffusive and convective momentum transport of the main ions of the plasma as well as the turbulent intrinsic torque in a future work.

The interaction of Alfvén eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25–250 kHz and toroidal mode numbers |n| < 20. During the 2019–2020 deuterium campaign, almost 7500 resonances and their frequencies f 0, net damping rates γ < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient’s effect on accessibility and ELM-related noise’s impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ion cyclotron resonance heating and neutral beam injection) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f 0, γ, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

Nonlinear 3D MHD asymmetric vertical displacement disruption simulations have been performed using JET equilibrium reconstruction initial data. Several experimentally measured quantities are compared with the simulation. These include vertical displacement, halo current, toroidal current asymmetry, and toroidal rotation. The experimental data and the simulations are in reasonable agreement. Also compared was the correlation of the toroidal current asymmetry and the vertical displacement asymmetry. The Noll relation between asymmetric wall force and vertical current moment is verified in the simulations. Also verified is the toroidal flux asymmetry. Although in many ways, JET is a good predictor of ITER disruption behavior, JET and ITER can be in different parameter regimes, and extrapolating from JET data can overestimate the ITER wall force.

The acceleration of the co-current toroidal rotations around resonant surfaces by resonant magnetic perturbations (RMPs) through turbulence is presented. These experiments were performed using a Langmuir probe array in the edge plasmas of the J-TEXT tokamak. This study aims at understanding the RMP effects on edge toroidal rotations and exploring its control method. With RMPs, the flat electron temperature T e profile, due to magnetic islands, appears around resonant surfaces (Zhao et al 2015 Nucl. Fusion 55 073022). When the resonant surface is closer to the last closed flux surface, the flat T e profile vanishes with RMPs. In both cases, the toroidal rotations significantly increase in the direction of the plasma current around the resonant surfaces with RMPs. The characteristics of turbulence are significantly affected by RMPs around the resonant surfaces. The turbulence intensity profile changes and the poloidal wave vector k θ increases with RMPs. The power fraction of the turbulence components in the ion diamagnetic drift direction increases with RMPs. The measurements of turbulent Reynolds stresses are consistent with the toroidal flows that can be driven by turbulence. The estimations of the energy transfer between the turbulence and toroidal flows suggest that turbulence energy transfers into toroidal flows. The result has the implication of the intrinsic rotation being driven by RMPs via turbulence.

The continuous spectrum of analytical toroidally rotating magnetically confined plasma equilibria is investigated analytically and numerically. In the presence of purely toroidal flow, the ideal magnetohydrodynamic equations leave the freedom to specify which thermodynamic quantity is constant on the magnetic surfaces. Introducing a general parametrization of this quantity, analytical equilibrium solutions are derived that still posses this freedom. These equilibria and their spectral properties are shown to be ideally suited for testing numerical equilibrium and stability codes including toroidal rotation. Analytical expressions are derived for the low-frequency continuous Alfvén spectrum. These expressions still allow one to choose which quantity is constant on the magnetic surfaces of the equilibrium, thereby generalizing previous results. The centrifugal convective effect is shown to modify the lowest Alfvén continuum branch to a buoyancy frequency, or Brunt–Väisälä frequency. A comparison with numerical results for the case that the specific entropy, the temperature, or the density is constant on the magnetic surfaces yields excellent agreement, showing the usefulness of the derived expressions for the validation of numerical codes.

We investigate by numerical simulation the rotational flows in a toroid confining a conducting magnetofluid in which a current is driven by the application of externally supported electric and magnetic fields. The computation involves no microscopic instabilities and is purely magnetohydrodynamic (MHD). We show how the properties and intensity of the rotations are regulated by dimensionless numbers (Lundquist and viscous Lundquist) that contain the resistivity and viscosity of the magnetofluid. At the magnetohydrodynamic level (uniform mass density and incompressible magnetofluids), rotational flows appear in toroidal, driven MHD. The evolution of these flows with the transport coefficients, geometry, and safety factor are described.

Strong suppression of impurity accumulation is observed in long pulse hydrogen discharges with high power NBI (neutral beam injection) heating (Pnbi > 10 MW) on the large helical device (LHD), even in the impurity accumulation window where the intrinsic impurities such as Fe and C are always accumulated into the plasma core. Density scan experiments in these discharges demonstrate to vanish the window and a new operational regime without impurity accumulation is found in steady state hydrogen discharges. Impurity pinch decreases with increasing ion temperature gradient and carbon Mach number. The peaking of the measured carbon profiles shows strong anti-correlations with the Mach number and its radial gradient. An external torque has a big impact on impurity transport and strong co-current rotation leads to an extremely hollow carbon profile, so-called 'impurity hole' observed in high ion temperature modes. Impurity pinch in the plasmas with net zero torque input (balanced NBI injection) is also strongly reduced by increasing ion temperature gradient, which can drive turbulent modes. The combination effect of turbulence and toroidal rotation plays an important role in the impurity transport.

Toroidal rotation can be produced by disruptions, as observed in several experiments. There is a concern that rotating asymmetric forces during an ITER disruption might resonate with the blanket and other structures surrounding the plasma. Here it is shown, both computationally using the M3D code, and analytically, that toroidal rotation is produced by magnetohydrodynamic turbulence. In particular, rotation is produced during an asymmetric vertical displacement event (AVDE) disruption. Toroidal and poloidal rotation are also produced during edge localized modes (ELMs), and may be consistent with a scaling law found for intrinsic toroidal rotation in H-mode tokamaks.

Fast ions are an inseparable part of fusion plasmas. They can be generated using electromagnetic waves or injected into plasmas as neutrals to heat the bulk plasma and to drive toroidal rotation and current. In future power plants fusion born fast ions deliver the main heating into the plasma. Understanding and controlling the fast ions is of crucial importance for the operation of a power plant. Furthermore, fast ions provide ways to probe the properties of the thermal plasma and get insight of its confinement properties. In this thesis, numerical code packages are used and developed to simulate JET experiments for a range of physics issues related to fast ions. Namely, the clamping fast ion distribution at high energies with RF heating, fast ion ripple torque generation and the toroidal momentum transport properties using NBI modulation technique are investi gated. Through a comparison of numerical simulations and the JET experimental data it is shown that the finite Larmor radius effects in ion cyclotron resonance heating are important and that they can prevent fast ion tail formation beyond certain energy. The identified mechanism could be used for tailoring the fast ion distribution in future experiments. Secondly, ASCOT simulations of NBI ions in a ripplefield showed that most of the reduction of the toroidal rotation that has been observed in the JET enhanced ripple experiments could be attributed to fast ion ripple torque. Finally, fast ion torque calculations together with momentum transport analysis have led to the conclusion that momentum transport in not purely diffusive but that a convective component, which increases monotonically in radius, exists in a wide range of JET plasmas. Using parameter scans, the convective transport has been shown to be insensitive to collisionality and q-profile but to increase strongly against density gradient.

The influence of toroidal plasma rotation on the existence of reversed shear Alfvén eigenmodes (RSAEs) near their minimum frequency is investigated analytically. An existence condition is derived showing that a radially decreasing kinetic energy density is unfavourable for the existence of RSAEs. The Coriolis effect is typically unfavourable for modes moving slower than the plasma or moving in the opposite direction. The generality of the analysis allows for the description of downwards sweeping RSAEs and also rotation-induced modes in regular shear plasmas.

Fast ions are an inseparable part of fusion plasmas. They can be generated using electromagnetic waves or injected into plasmas as neutrals to heat the bulk plasma and to drive toroidal rotation and current. In future power plants fusion born fast ions deliver the main heating into the plasma. Understanding and controlling the fast ions is of crucial importance for the operation of a power plant. Furthermore, fast ions provide ways to probe the properties of the thermal plasma and get insight of its confinement properties. In this thesis, numerical code packages are used and developed to simulate JET experiments for a range of physics issues related to fast ions. Namely, the clamping fast ion distribution at high energies with RF heating, fast ion ripple torque generation and the toroidal momentum transport properties using NBI modulation technique are investi gated. Through a comparison of numerical simulations and the JET experimental data it is shown that the finite Larmor radius effects in ion cyclotron resonance heating are important and that they can prevent fast ion tail formation beyond certain energy. The identified mechanism could be used for tailoring the fast ion distribution in future experiments. Secondly, ASCOT simulations of NBI ions in a ripplefield showed that most of the reduction of the toroidal rotation that has been observed in the JET enhanced ripple experiments could be attributed to fast ion ripple torque. Finally, fast ion torque calculations together with momentum transport analysis have led to the conclusion that momentum transport in not purely diffusive but that a convective component, which increases monotonically in radius, exists in a wide range of JET plasmas. Using parameter scans, the convective transport has been shown to be insensitive to collisionality and q-profile but to increase strongly against density gradient.

The continuous spectrum of analytical toroidally rotating magnetically confined plasma equilibria is investigated analytically and numerically. In the presence of purely toroidal flow, the ideal magnetohydrodynamic equations leave the freedom to specify which thermodynamic quantity is constant on the magnetic surfaces. Introducing a general parametrization of this quantity, analytical equilibrium solutions are derived that still posses this freedom. These equilibria and their spectral properties are shown to be ideally suited for testing numerical equilibrium and stability codes including toroidal rotation. Analytical expressions are derived for the low-frequency continuous Alfven spectrum. These expressions still allow one to choose which quantity is constant on the magnetic surfaces of the equilibrium, thereby generalizing previous results. The centrifugal convective effect is shown to modify the lowest Alfven continuum branch to a buoyancy frequency, or Brunt-Vaisala frequency. A comparison with numerical results for the case that the specific entropy, the temperature, or the density is constant on the magnetic surfaces yields excellent agreement, showing the usefulness of the derived expressions for the validation of numerical codes. (C) 2011 Elsevier Inc. All rights reserved.

The influence of toroidal plasma rotation on the existence of reversed shear Alfvén eigenmodes (RSAEs) near their minimum frequency is investigated analytically. An existence condition is derived showing that a radially decreasing kinetic energy density is unfavourable for the existence of RSAEs. The Coriolis effect is typically unfavourable for modes moving slower than the plasma or moving in the opposite direction. The generality of the analysis allows for the description of downwards sweeping RSAEs and also rotation-induced modes in regular shear plasmas.

Flux-driven global gyrokinetic codes are now mature enough to make predictions in terms of turbulence and transport in tokamak plasmas. Some of the recent breakthroughs of three such codes, namely GYSELA, ORB5 and XGC1, are reported and compared wherever appropriate. In all three codes, turbulent transport appears to be mediated by avalanche-like events, for a broad range of ρ* = ρi/a values, ratio of the gyro-radius over the minor radius. Still, the radial correlation length scales with ρi, leading to the gyro-Bohm scaling of the effective transport coefficient below ρ* ≈ 1/300. The possible explanation could be due to the fact that avalanches remain meso-scale due to the interaction with zonal flows, whose characteristic radial wavelength appears to be almost independent of the system size. As a result of the radial corrugation of the turbulence driven zonal and mean flows, the shear of the radial electric field can be significantly underestimated if poloidal rotation is assumed to be governed by the neoclassical theory, especially at low collisionality. Indeed, the turbulence contribution to the poloidal rotation increases when collisionality decreases. Finally, the numerical verification of toroidal momentum balance shows that both neoclassical and turbulent contributions to the Reynolds' stress tensor play the dominant role. The phase space analysis further reveals that barely passing supra-thermal particles mostly contribute to the toroidal flow generation, consistently with quasi-linear predictions.

As a result of toroidal rotation, sequences of new global modes are predicted to arise in magnetically confined plasmas. The frequencies of these Alfv\'en modes lie inside gaps of the continuous magnetohydrodynamic (MHD) spectrum that are created or enlarged by toroidal flow. The numerically obtained results are compared with an analytical investigation, yielding a useful criterion for mode existence. Because of their low frequencies, these modes may be easily destabilized by energetic particles. Because of their sensitivity to the Mach number however, these modes can provide a valuable extension to MHD spectroscopy by giving information on the rotational velocity.

As a result of toroidal rotation, sequences of new global modes are predicted to arise in magnetically confined plasmas. The frequencies of these Alfv\'en modes lie inside gaps of the continuous magnetohydrodynamic (MHD) spectrum that are created or enlarged by toroidal flow. The numerically obtained results are compared with an analytical investigation, yielding a useful criterion for mode existence. Because of their low frequencies, these modes may be easily destabilized by energetic particles. Because of their sensitivity to the Mach number however, these modes can provide a valuable extension to MHD spectroscopy by giving information on the rotational velocity.

The ideal magnetohydrodynamic stability is investigated of localized interchange modes in a large-aspect ratio tokamak plasma. The resulting stability criterion includes the effects of toroidal rotation and rotation shear and contains various well-known limiting cases. The analysis allows for a general adiabatic index, resulting in a stabilizing contribution from the convective effect. A further stabilizing effect from rotation exists when the angular frequency squared decreases radially more rapidly than the density. Flow shear, however, also decreases the stabilizing effect of magnetic shear through the Kelvin– Helmholtz mechanism. Numerical simulations reveal the merits and limitations of the performed local analysis.

The ideal magnetohydrodynamic stability is investigated of localized interchange modes in a large-aspect ratio tokamak plasma. The resulting stability criterion includes the effects of toroidal rotation and rotation shear and contains various well-known limiting cases. The analysis allows for a general adiabatic index, resulting in a stabilizing contribution from the convective effect. A further stabilizing effect from rotation exists when the angular frequency squared decreases radially more rapidly than the density. Flow shear, however, also decreases the stabilizing effect of magnetic shear through the Kelvin– Helmholtz mechanism. Numerical simulations reveal the merits and limitations of the performed local analysis.