Search

A nonlinear simulation of the energetic particle driven instabilities in the Chinese First Quasi-axisymmetric Stellarator (CFQS) has been conducted for the first time. MEGA, a hybrid simulation code for energetic particles interacting with a magneto-hydrodynamic fluid, was used in the present work. Both the m/n = 3/1 energetic-particle-mode (EPM) like mode and the m/n = 5/2 toroidal Alfvén eigenmode (TAE) were found, where m is the poloidal mode number and n is the toroidal mode number. Four important results were obtained as follows. First, the instability in the CFQS in three-dimensional form was shown for the first time. Second, strong toroidal mode coupling was found for the spatial profiles, and it is consistent with the theoretical prediction. Third, the resonant condition caused by the absence of axial symmetry in CFQS was demonstrated for the first time. The general resonant condition is f mode = Nf ϕ − Lf θ , where f mode, f ϕ , and f θ are mode frequency, particle toroidal transit frequency, and particle poloidal transit frequency, respectively; N and L are arbitrary integers, represent toroidal and poloidal resonance numbers. For EPM-like mode, the dominant and subdominant resonant conditions are f mode = 3f ϕ − 7f θ and f mode = f ϕ − f θ , respectively. For TAE, the dominant and subdominant resonant conditions are f mode = 4f ϕ − 9f θ and f mode = 2f ϕ − 3f θ , respectively. On the one hand, the toroidal resonance numbers are different from the toroidal mode numbers by 2. This indicates that the two-fold rotational symmetry affects the resonance condition. On the other hand, the subdominant resonances satisfy N = n, which is expected for the axisymmetric plasmas and most of the toroidal plasmas including stellarators. Fourth, the nonlinear frequency chirpings in CFQS were demonstrated for the first time. Hole and clump structures were formed in the pitch angle and energy phase space, and the particles comprising the hole and clump were kept resonant with the modes during the mode frequencies chirping.

Beam modulation effects on Alfvén eigenmode stability have been investigated in a recent DIII-D experiment and show that variations in neutral beam modulation period can have an impact on the beam driven Alfvén eigenmode spectrum and resultant fast ion transport despite similar time-averaged input power. The experiment was carried out during the current ramp phase of L-mode discharges heated with sub-Alfvénic 50–80 kV deuterium neutral beams that drive a variety of Alfvén eigenmodes unstable. The modulation period of two interleaved beams with different tangency radii was varied from shot to shot in order to modify the relative time dependent mix of the beam pitch angle distribution as well as the persistence of a bump-on-tail feature near the injection energy (a feature confirmed by imaging neutral particle analyzer measurements). As the beam modulation period is varied from 7 ms to 30 ms on/off (typical full energy slowing down time of τ slow ≈ 50 ms at mid-radius), toroidicity-induced Alfvén eigenmodes (TAEs) located in the outer periphery of the plasma become intermittent and coincident with the more tangential beam. Core mode activity changes from reversed shear Alfvén eigenmodes (RSAEs) to a mix of RSAE and beta-induced Alfvén eigenmodes. Discharges with 30 ms on/off period do not have a persistent bump-on-tail feature, have the lowest average mode amplitude and least fast ion transport. Detailed analysis of an individual TAE using TRANSP kick modeling (Monte Carlo evolution of the distribution function with probabilistic ‘kicks’ by the AEs) and the resistive MHD code with kinetic fast ions, MEGA, find no strong role of energy gradient drive due to bump-on-tail features. Instead, the observed TAE modulation with interleaved beams is likely a pitch angle dependent result combined with slowing down of the tangential beam between pulses. For the conditions investigated, bump-on-tail contributions to TAE drive were found to be 5% or less of the total drive at any given time.

Spherical tokamak (ST) research in Japan has produced many innovative results: (i) plasma start-up to I p > 70 kA was achieved by electron cyclotron wave (ECW) with N ∥ = 0.75, while electron heating to T e > 500 eV was achieved with N ∥ = 0.26 on QUEST. (ii) The radiofrequency (RF)-induced transport model was combined with the x-ray emission model, and extended magnetohydrodynamics equilibrium with kinetic electrons was developed to interpret fast-electron-dominated lower hybrid wave sustained plasmas on TST-2. (iii) Density as high as 30 times the cutoff density was achieved by electron Berstein wave current drive combined with electron beam injection on LATE. (iv) Multiple plasmoids formed by tearing instability in the elongated current sheet were observed, and flux closure and ion heating by plasmoid-mediated fast magnetic reconnection were observed on HIST. (v) Optimization of ECW-assisted inductive start-up with a vertical field with positive decay index was performed on TST-2. (vi) Stabilization of the vertical displacement event by a set of upper and lower helical field coils was demonstrated on TOKASTAR-2. (vii) A 6 h discharge was achieved by cool-down of the center stack cover on QUEST, where the plasma duration limit was consistent with the wall saturation time estimated by modeling. (viii) Extension of ion heating by plasma merging was achieved on TS-3U, TS-4U, UTST, MAST, and ST40.

Kinetic-magnetohydrodynamic hybrid simulations were performed to investigate the linear growth and the nonlinear evolution of off-axis fishbone mode (OFM) destabilized by trapped energetic ions in tokamak plasmas. The spatial profile of OFM is mainly composed of m/n = 2/1 mode inside the q = 2 magnetic flux surface while the m/n = 3/1 mode is predominant outside the q = 2 surface, where m and n are the poloidal and toroidal mode numbers, respectively, and q is the safety factor. The spatial profile of the OFM is a strongly shearing shape on the poloidal plane, suggesting the nonperturbative effect of the interaction with energetic ions. The frequency of the OFM in the linear growth phase is in good agreement with the precession drift frequency of trapped energetic ions, and the frequency chirps down in the nonlinear phase. Two types of resonance conditions between trapped energetic ions and OFM are found. For the first type of resonance, the precession drift frequency matches the OFM frequency, while for the second type, the sum of the precession drift frequency and the bounce frequency matches the OFM frequency. The first type of resonance is the primary resonance for the destabilization of OFM. The resonance frequency which is defined based on precession drift frequency and bounce frequency of the nonlinear orbit for each resonant particle is analyzed to understand the frequency chirping. The resonance frequency of the particles that transfer energy to the OFM chirps down, which may result in the chirping down of the OFM frequency. A detailed analysis of the energetic ion distribution function in phase space shows that the gradient of the distribution function along the E′ = const. line drives or stabilizes the instability, where E′ is a combination of energy and toroidal canonical momentum and conserved during the wave–particle interaction. The distribution function is flattened along the E′ = const. line in the nonlinear phase leading to the saturation of the instability.

A transition from an interchange mode to a non-resonant mode is found in the nonlinear magnetohydrodynamic simulation for the partial collapse in a large helical device (LHD) plasma with a net toroidal current. This transition can occur when the magnetic shear is weak and the rotational transform is close to unity in the core region. In this transition, the mode number of the dominant Fourier component is reduced. As a result of the nonlinear evolution, the (m, n) = (1, 1) component can be dominant, where m and n are the poloidal and the toroidal mode numbers, respectively. This transition is considered to be a candidate to explain the observation in the LHD experiments with the net toroidal current that show partial collapses are caused by the (1, 1) mode.

Alfvén eigenmodes (AEs) driven by energetic electrons were investigated via hybrid simulations of an MHD fluid interacting with energetic electrons. The investigation focused on AEs with the toroidal number n = 4. Both energetic electrons with centrally peaked beta profile and off-axis peaked profile are considered. For the centrally peaked energetic electron beta profile case, a toroidal Alfvén eigenmode (TAE) propagating in the electron diamagnetic drift direction is found. The mode is mainly driven by deeply trapped energetic electrons. It is also found that a few passing energetic electrons spatially localized around rational surfaces can resonate with the mode. For the off-axis peaked energetic electron beta profile case, an AE propagating in the ion diamagnetic drift direction is found when a q − p r o f i l e with weak magnetic shear is adopted. The destabilized mode is an elliptical-Alfvén-eigenmode-type (EAE-type) mode which has a spatial profile peaking at the rational surface and a frequency close to the second Alfvén frequency gap. It is found that passing energetic electrons and barely trapped energetic electrons are responsible for this EAE-type mode destabilization. The saturation levels are compared for a TAE with the same linear growth rate among energetic electron driven mode and energetic ion driven mode with isotropic and anisotropic velocity space distributions. The saturation level of TAE driven by trapped energetic electrons is comparable to that driven by energetic electrons with isotropic velocity space distribution where the contribution of trapped particles is dominant. It is found that the trapped energetic ion driven TAE has a larger saturation level than the passing energetic ion driven TAE, which indicates the difference in particle trapping by the TAE between trapped and passing energetic ions.

Energetic-particle-driven geodesic acoustic modes (EGAMs) channeling in the Large Helical Device (LHD) plasmas are systematically investigated for the first time using MEGA code. MEGA is a hybrid simulation code for energetic particles interacting with a magnetohydrodynamic (MHD) fluid. In the present work, both the energetic particles and the bulk ions are described kinetically. The EGAM profiles in the three-dimensional form is illustrated. Then, EGAM channeling behaviors are analyzed under different conditions. During the EGAM activities without frequency chirping, EGAM channeling occurs in the linear growth stage but terminates in the decay stage after the saturation. During the EGAM activities with frequency chirping, EGAM channeling occurs continuously. Also, low-frequency EGAM makes the energy transfer efficiency (Eion/EEP) higher, and this is confirmed by changing the energetic particle pressure, energetic particle beam velocity, and energetic particle pitch angle. Moreover, higher bulk ion temperature makes the energy transfer efficiency higher. In addition, under a certain condition, the energy transfer efficiency in the deuterium plasma is lower than that in the hydrogen plasma.

This report summarizes the contributions presented at the 16th IAEA Technical Meeting on Energetic Particles in Magnetic Confinement Systems—Theory of Plasma Instabilities, held in Shizuoka, Japan, 3–6 September 2019. The meeting brought together about 100 experts from nuclear fusion research sites worldwide to discuss the physics of energetic particles and plasma instabilities, at the first joint meeting from the two scientific disciplines. The main topics of the meeting were: alpha particles physics; transport of energetic particles; effects of energetic particles in magnetic confinement fusion devices; collective phenomena such as Alfven eigenmodes, energetic particle modes and others; runaway electrons and disruptions; diagnostics for energetic particles; control of energetic particles confinement; multiscale physics and instabilities in burning plasmas.

Energetic electron effects on an energetic-ion driven toroidal Alfvén eigenmode (TAE) are investigated via hybrid simulations of an MHD fluid interacting with energetic particles. Both energetic electrons and energetic ions described by drift-kinetic equations are included in the present work. It is found that the TAE can be effectively stabilized by off-axis peaked energetic electrons which are located near the mode center, while the centrally peaked energetic electrons fail to stabilize the mode. It is confirmed that the spatially localized pressure profile of energetic electrons causes the stabilization of TAE. The stabilized TAE has a more localized mode structure accompanied by a significant reduction in the energetic ion driving rate. The small change of mode frequency and dissipation rate indicate the stabilization mechanism is different from the so-called pressure gradient stabilization that drives the TAE into continuum. The results suggest that the strong plasma non-uniformity induced by the energetic electron beta profile may be responsible for the change of mode structure. It is also found that this stabilizing effect is more effective for a high-n TAE. Moreover, it is numerically verified that the positive (negative) pressure gradient at the TAE center will increase (decrease) the mode frequency. The wave-particle interactions are also analysed for a case with energetic electrons peaked at the inner side of the TAE center. It is found that the power transfer to a resonant barely trapped energetic electron, which taps energy from the wave, can be comparable to the power transfer from a resonant energetic ion. This suggests that if a sufficient number of resonant barely trapped electrons are present, they might stabilize energetic-ion driven TAE through the wave-particle interaction.

Comprehensive computer simulations of the Alfven eigenmode burst, which is the synchronized sudden growth of multiple Alfven eigenmodes (AEs) interacting with energetic particles, were conducted with continuous neutral beam injection, collisions, and particle losses. It is found that the energetic-particle distribution in phase space reaches a 'critical distribution' with a stairway structure where a resonance overlap triggers the Alfven eigenmode burst. Before the burst, the gradual growth of the AEs associated with the beam injection broadens the resonant regions in phase space forming the distribution into a stairway shape. When the distribution reaches the 'critical distribution,' a resonance overlap triggers multiple resonance overlaps leading to the synchronized growth of AEs and the collapse of the distribution. For another run with the beam deposition power reduced to one-half, the energetic-particle distribution function just before the Alfven eigenmode burst is close to that for the original beam power. This result indicates that the critical distribution for the Alfven eigenmode burst is present.

Energetic particle driven geodesic acoustic modes (EGAMs) in the Large Helical Device (LHD) plasmas are investigated using MEGA code. MEGA is a hybrid simulation code for energetic particles interacting with a magnetohydrodynamic (MHD) fluid and in the present work, both the energetic particles and bulk ions are described by the kinetic equations. The low frequency EGAMs are reproduced. Also, the energy transfer is analyzed and the bulk ion heating during the EGAM activity is observed. The ions obtain energy when the energetic particles lose energy, and this indicates that an energy channel is established by EGAM. The EGAM channeling is reproduced by simulation with realistic parameters for the first time. The heating power to bulk ions is 3.4 kW/malt;supagt;3alt;/supagt;. It is found that the sideband resonance is dominant during the energy transfer from EGAM to the bulk ions, and the transit frequencies of resonant bulk ions are one-half of the EGAM frequency.

The linear stability analysis of resistive ballooning modes in Large Helical Device plasmas has been carried out by numerical simulation based on the kinetic magnetohydrodynamic (MHD) model where thermal ions are simulated with the drift-kinetic model. Due to the curvature drift and the gradient B drift, the orbit of trapped ions not only deviates from the magnetic surface but also has the precession motion in the poloidal and toroidal directions. Even when the trapped ions remain well inside the radial profile of the MHD mode, the response of the trapped ions is weakened when the precession drift frequency of the trapped ions with respect to the mode phase is close to or larger than the linear growth rate of the instability. This results in the suppression of the ion perpendicular pressure perturbation to the magnetic field leading to the reduction of the linear growth rate.

Comprehensive magnetohydrodynamic (MHD) hybrid simulations with neutral beam injection and collisions were conducted to investigate the Alfven eigenmode (AE) bursts and the fast-ion losses in the large helical device (LHD) for the realistic conditions close to the experiments. It is found in the simulation of the slowing-down time scale that the AE bursts take place repetitively accompanied by fast-ion redistribution and losses leading to lower saturation levels of stored fast-ion energy than those in a classical calculation where the MHD perturbations are neglected. The fast-ion loss rate caused by the AE burst has the quadratic dependence on AE amplitude, which was observed in the LHD experiment. The majority of the lost fast ions are counter-passing particles whose velocity and pitch-angle are close to those of the beam injection. The second component of the lost fast ions is transit particles whose velocity is close to thermal velocity. The loss of the counter-passing particles occurs mainly during the AE bursts, while the transit particles are lost both during the AE bursts and the quiescent periods with larger loss rate than that in the classical calculation. The initial location of the lost counter-injected particles spreads from the plasma edge to the plasma center, while only the particles initially located in the peripheral region are lost for the co-injected beam.

This paper reports verification and validation of linear simulations of Alfvén eigenmodes in the current ramp phase of DIII-D L-mode discharge #159243 using gyrokinetic, gyrokinetic-MHD hybrid, and eigenvalue codes. Using a classical fast ion profile, all simulation codes find that reversed shear Alfvén eigenmodes (RSAE) are the dominant instability. The real frequencies from all codes have a coefficient of variation of less than 5% for the most unstable modes with toroidal mode number n = 4 and 5. The simulated RSAE frequencies agree with experimental measurements if the minimum safety factor is adjusted, within experimental errors. The simulated growth rates exhibit greater variation, and simulations find that pressure gradients of thermal plasmas make a significant contribution to the growth rates. Mode structures of the dominant modes agree well among all codes. Moreover, using a calculated fast ion profile that takes into account the diffusion by multiple unstable modes, a toroidal Alfvén eigenmode (TAE) with n = 6 is found to be unstable in the outer edge, consistent with the experimental observations. Variations of the real frequencies and growth rates of the TAE are slightly larger than those of the RSAE. Finally, electron temperature fluctuations and radial phase shifts from simulations show no significant differences with the experimental data for the strong n = 4 RSAE, but significant differences for the weak n = 6 TAE. The verification and validation for the linear Alfvén eigenmodes is the first step to develop an integrated simulation of energetic particles confinement in burning plasmas incorporating multiple physical processes.

Fast particles in fusion plasmas may drive Alfvén modes unstable leading to fluctuations of the internal electromagnetic fields and potential loss of particles. Such instabilities can have an impact on the performance and the wall-load of machines with burning plasmas such as ITER. A linear benchmark for a toroidal Alfvén eigenmode (TAE) is done with 11 participating codes with a broad variation in the physical as well as the numerical models. A reasonable agreement of around 20% has been found for the growth rates. Also, the agreement of the eigenfunctions and mode frequencies is satisfying. However, they are found to depend strongly on the complexity of the used model. © EURATOM 2018.

The confinement of energetic ions generated by neutral beam injection (NBI) and ion cyclotron resonance frequency heating is studied using GNET simulation code, in which the drift kinetic equation ...

Energetic ion-driven magnetohydrodynamic instabilities such as Alfven eigenmodes (AEs), energetic particle modes (EPMs), and their impacts on energetic ion confinement are being studied on the Larg...

This chapter reviews the progress accomplished since the redaction of the first ITER Physics Basis (1999 Nucl. Fusion 39 2137-664) in the field of energetic ion physics and its possible impact on burning plasma regimes. New schemes to create energetic ions simulating the fusion-produced alphas are introduced, accessing experimental conditions of direct relevance for burning plasmas, in terms of the Alfvénic Mach number and of the normalised pressure gradient of the energetic ions, though orbit characteristics and size cannot always match those of ITER. Based on the experimental and theoretical knowledge of the effects of the toroidal magnetic field ripple on direct fast ion losses, ferritic inserts in ITER are expected to provide a significant reduction of ripple alpha losses in reversed shear configurations. The nonlinear fast ion interaction with kink and tearing modes is qualitatively understood, but quantitative predictions are missing, particularly for the stabilisation of sawteeth by fast particles that can trigger neoclassical tearing modes. A large database on the linear stability properties of the modes interacting with energetic ions, such as the Alfvén eigenmode has been constructed. Comparisons between theoretical predictions and experimental measurements of mode structures and drive/damping rates approach a satisfactory degree of consistency, though systematic measurements and theory comparisons of damping and drive of intermediate and high mode numbers, the most relevant for ITER, still need to be performed. The nonlinear behaviour of Alfvén eigenmodes close to marginal stability is well characterized theoretically and experimentally, which gives the opportunity to extract some information on the particle phase space distribution from the measured instability spectral features. Much less data exists for strongly unstable scenarios, characterised by nonlinear dynamical processes leading to energetic ion redistribution and losses, and identified in nonlinear numerical simulations of Alfvén eigenmodes and energetic particle modes. Comparisons with theoretical and numerical analyses are needed to assess the potential implications of these regimes on burning plasma scenarios, including in the presence of a large number of modes simultaneously driven unstable by the fast ions. © 2007 IAEA, Vienna.

The interaction between pressure driven modes and magnetic islands generated by a resonant magnetic perturbation (RMP) in the large helical device (LHD) is numerically analyzed. In this analysis, three-dimensional treatment is essential in the equilibrium and dynamics calculations, because the equilibrium pressure profile is deformed by the RMP. The deformation changes the linear mode structure from the interchange type to the ballooning-like type that is localized around the X-point of the island in the equilibrium magnetic field including the RMP. This mode causes a pressure collapse in the nonlinear evolution, which spreads from the X-point to the core. Therefore, the spatial phase of the collapse is fixed to the island geometry. The fixed phase agrees with the LHD experimental results with a natural error field.

Characteristics of the MHD instabilities for high beta LHD (large helical device) plasmas in the inward shifted configurations have been investigated by numerical simulations. On the condition that the magnetic Reynolds number is lower than the experimental value, the time evolution of the MHD instabilities is as follows: (1) resistive ballooning modes, the toroidal mode number of which is higher than the LHD's toroidal pitch (), are destabilized in the plasma peripheral region; (2) low toroidal modes typified by are destabilized by the nonlinear mode coupling of the resistive ballooning modes the toroidal mode number of which is adjacent to each other where n is the toroidal mode number. In particular, the velocity of the mode is finite at the magnetic axis so that there is a large velocity directed to the magnetic axis; and (3) when the high n modes are saturated, the destabilized low n modes, which have a global mode structure and shift the magnetic axis, also begin to be saturated. The self mode coupling of the low n modes induces the core crush and transports the plasma from the core region to the peripheral region. Since the nonlinear MHD phenomena are dominated by the resistive modes, the MHD phenomena for the experimental high magnetic Reynolds number are considered to be milder than our numerical results. Thus, when the experimental high magnetic Reynolds number is taken in the simulation, the core crush may be suppressed so that the numerical results are expected to become close to the experimental results where the stable high beta plasmas are obtained.

As the finalization of a hydrogen experiment towards the deuterium phase, the exploration of the best performance of hydrogen plasma was intensively performed in the large helical device. High ion and electron temperatures, Ti and Te, of more than 6 keV were simultaneously achieved by superimposing high-power electron cyclotron resonance heating onneutral beam injection (NBI) heated plasma. Although flattening of the ion temperature profile in the core region was observed during the discharges, one could avoid degradation by increasing the electron density. Another key parameter to present plasma performance is an averaged beta value $\left\langle \beta \right\rangle $ . The high $\left\langle \beta \right\rangle $ regime around 4% was extended to an order of magnitude lower than the earlier collisional regime. Impurity behaviour in hydrogen discharges with NBI heating was also classified with a wide range of edge plasma parameters. The existence of a no impurity accumulation regime, where the high performance plasma is maintained with high power heating >10 MW, was identified. Wide parameter scan experiments suggest that the toroidal rotation and the turbulence are the candidates for expelling impurities from the core region.

© 2014 IAEA. A multi-phase simulation that is a combination of classical simulation and hybrid simulation for energetic particles interacting with a magnetohydrodynamic (MHD) fluid is developed to simulate the nonlinear dynamics on the slowing down time scale of the energetic particles. The hybrid simulation code is extended with realistic beam deposition profile, collisions and losses, and is used for both the classical and hybrid phases. The code is run without MHD perturbations in the classical phase, while the interaction between the energetic particles and the MHD fluid is simulated in the hybrid phase. In a multi-phase simulation of DIII-D discharge #142111, the stored beam ion energy is saturated due to Alfvén eigenmodes (AE modes) at a level lower than in the classical simulation. After the stored fast ion energy is saturated, the hybrid simulation is run continuously. It is demonstrated that the fast ion spatial profile is significantly flattened due to the interaction with the multiple AE modes with amplitude v/vA∼ δB/B ∼ O(10-4). The dominant AE modes are toroidal Alfvén eigenmodes (TAE modes), which is consistent with the experimental observation at the simulated moment. The amplitude of the temperature fluctuations brought about by the TAE modes is of the order of 1% of the equilibrium temperature. This is also comparable with electron cyclotron emission measurements in the experiment.

Recurrent bursts of toroidicity induced Alfven eigenmodes (TAEs) are reproduced with a Fokker-Planck-MHD simulation where a fast ion source and slowing down are incorporated self-consistently. The bursts take place at regular time intervals and the behaviours of all the TAEs are synchronized. The fast ion transport due to TAE activity spatially broadens the classical fast ion distribution and significantly reduces its peak value. Only a small change of the distribution takes place with each burst, leading to loss of a small fraction of the fast ions. The system stays close to the marginal stability state established through the interplay of the fast ion source, slowing down and TAE activity.

Remarkable progress has been made in diagnosing energetic particle instabilities on present-day machines and in establishing a theoretical framework for describing them. This overview describes the much improved diagnostics of Alfven instabilities and modelling tools developed world-wide, and discusses progress in interpreting the observed phenomena. A multi-machine comparison is presented giving information on the performance of both diagnostics and modelling tools for different plasma conditions outlining expectations for ITER based on our present knowledge., 37 pages, 20 figures

Recurring bursts of chirping Alfvén modes that were observed in JT-60U tokamak plasmas driven by negative-ion-based neutral beams (N-NB) are reproduced in first-principle simulations performed with an extended version of the hybrid code MEGA. This code simulates the interactions between gyrokinetic fast ions and magnetohydrodynamic (MHD) modes in the presence of a realistic fast ion source and collisions, so that it self-consistently captures dynamics across a wide range of time scales (0.01–100 ms). The simulation confirms that the experimentally observed phenomena known as 'fast frequency sweeping (fast FS) modes' are caused by bursts of energetic particle modes (EPM) with dominant toroidal mode number n = 1. On the long time scale (1–10 ms), the simulation reproduces the chirping range (40–60 kHz), the burst duration (few ms) and intervals (5–10 ms). On the short time scale (0.01–0.1 ms), it reproduces pulsations and phase jumps, which we interpret as the result of beating between multiple resonant wave packets. Having reproduced at multiple levels of detail the dynamics of low-amplitude long-wavelength Alfvén modes driven by N-NB ions, the next goal is to reproduce and explain abrupt large-amplitude events (ALE) that were seen in the same experiments at longer time intervals (10–100 ms).

The sensitivity of the growth and nonlinear evolution of fast-ion-driven modes is examined with respect to the choice of particle boundary conditions, diffusion coefficients, fast ion gyroradii and bulk compressibility. The primary purpose of this work is to justify the choice of parameters to be used in the self-consistent long-time simulations of fast ion dynamics using global MHD-kinetic hybrid codes that include fast ion sources and collisions. The present study is conducted for a scenario based on the N-NB-driven JT-60U shot E039672, which is subject to abrupt large events (ALE). We use realistic geometry, a realistic fast ion distribution, and focus on experimentally observed harmonics with low toroidal mode numbers n = 1, 2, 3. The use of realistic boundary conditions and finite Larmor radii for the fast ions is shown to be essential. The usual values μ 0 η = ν = χ ∼ 10 − 6 v A 0 R 0 used for resistivity, viscosity and thermal diffusivity, and Γ = 5 / 3 used for the specific heat ratio (controlling the effect of compressibility) are shown to be reasonable choices. Our method for performing the parameter scans around the threshold for the onset of convective amplification is proposed as a strategy for nonlinear benchmark studies.

Fast ion pressure profiles flattened by multiple Alfvén eigenmodes (AEs) are investigated for various neutral beam deposition powers in a multi-phase simulation, which is a combination of classical simulation and hybrid simulation for energetic particles interacting with a magnetohydrodynamic fluid. Monotonic degradation of fast ion confinement and fast ion profile stiffness is found with increasing beam deposition power. The confinement degradation and profile stiffness are caused by a sudden increase in fast ion transport flux brought about by AEs for fast ion pressure gradients above a critical value. The critical pressure gradient and the corresponding beam deposition power depend on the radial location. The fast ion pressure gradient stays moderately above the critical value, and the profiles of the fast ion pressure and fast ion transport flux spread radially outward from the inner region, where the beam is injected. It is found that the square root of the MHD fluctuation energy is proportional to the beam deposition power. Analysis of the time evolutions of the fast ion energy flux profiles reveals that intermittent avalanches take place with contributions from the multiple eigenmodes. Surface of section plots demonstrate that the resonance overlap of multiple eigenmodes accounts for the sudden increase in fast ion transport with increasing beam power. The critical gradient and critical beam power for the profile stiffness are substantially higher than the marginal stability threshold.

A multi-phase simulation, which is a combination of classical simulation and hybrid simulation for energetic particles interacting with a magnetohydrodynamic (MHD) fluid including neutral beam injection, slowing-down, and pitch angle scattering, is applied to DIII-D discharge #142111 where the fast ion spatial profile is significantly flattened due to multiple Alfven eigenmodes (AEs). The large fast ion pressure profile flattening observed experimentally is successfully reproduced by these first of a kind comprehensive simulations. Temperature fluctuations due to three of the dominant toroidal Alfven eigenmodes in the simulation results are compared in detail with electron cyclotron emission measurements in the experiment. It is demonstrated that the temperature fluctuation profile and the phase profile are in very good agreement with the measurement, and the amplitude is also in agreement within a factor of two. This level of agreement validates the multi-phase hybrid simulation for the prediction of AE activity and alpha particle transport in burning plasmas.

This conference report summarizes the contributions to, and discussions at, the 4th Asia-Pacific Transport Working Group Meeting held at Kyushu University, Japan, during 10–13 June 2014. The topics of the meeting were organized under five main headings: turbulence suppression and transport barrier formation, effect of magnetic topology on MHD activity and transport, non-diffusive contribution of momentum and particle transport, non-local transport and turbulence spreading and coupling, energetic particles and instability. The Young Researchers' Forum which was held in this meeting is also described in this report.

Dynamics of fast ions and shear Alfven waves are simulated using MEGA, a global nonlinear hybrid code. The scenario considered is based on JT-60U shot E039672, driven by strong negative-ion-based neutral beams (N-NB), just before the onset of a so-called abrupt large event (ALE). It is found that modes with toroidal mode numbers n = 2, 3, 4 can be destabilized, besides the n = 1 mode studied previously. The properties of the modes with n > 1 are sensitive to the value of the plasma beta and the form of the fast ion distribution, so simulation conditions are set up as realistically as presently possible. When the fast ion drive exceeds a certain threshold, the n = 3 mode is enhanced through a convective amplification process while following fast ions were displaced either by the n = 3 mode itself or by the n = 1 mode. The fast ion transport in several cases, simulated with single- or multiple-n modes, is analysed and implications of the results for the explanation of ALEs are discussed.

The progress of physical understanding as well as parameter improvement of net-current-free helical plasma is reported for the Large Helical Device since the last Fusion Energy Conference in Daejeon in 2010. The second low-energy neutral beam line was installed, and the central ion temperature has exceeded 7 keV, which was obtained by carbon pellet injection. Transport analysis of the high-T-i plasmas shows that the ion-thermal conductivity and viscosity decreased after the pellet injection although the improvement does not last long. The effort has been focused on the optimization of plasma edge conditions to extend the operation regime towards higher ion temperature and more stable high density and high beta. For this purpose a portion of the open helical divertors are being modified to the baffle-structured closed ones aimed at active control of the edge plasma. It is compared with the open case that the neutral pressure in the closed helical divertor increased by ten times as predicted by modelling. Studies of physics in a three-dimensional geometry are highlighted in the topics related to the response to a resonant magnetic perturbation at the plasma periphery such as edge-localized-mode mitigation and divertor detachment. Novel approaches of non-local and non-diffusive transport have also been advanced.

The hybrid code MEGA is used to simulate an abrupt large-amplitude event (ALE) in JT-60U shot E039672 that was strongly driven by a pair of negative-ion-based neutral beams (N-NBs). This configuration is known to be unstable to energetic particle modes (EPMs) with toroidal mode number n = 1. The purpose of this study is to look for a threshold with respect to the on-axis fast ion beta value, βh0, beyond which the EPM undergoes convective amplification (CA) and causes significant fast ion transport. This is motivated by the hypothesis that such a threshold may work as a trigger mechanism for relaxation events, such as ALE. In order to facilitate quantitative comparisons with the experiment, a realistic geometry and bulk pressure is used. Furthermore, MEGA is initialized with a fast ion distribution computed for JT-60U by an orbit-following Monte-Carlo code, and a passive vacuum region allows particles to travel without encountering artificial loss boundaries. Consistently with the experiment, the simulation predicts a burst time of about 0.3 ms and peak magnetic fluctuation levels around δBϑ/B

This paper examines nonlinear magneto-hydrodynamic effects on the energetic particle driven toroidal Alfvén eigenmode (TAE) for lower dissipation coefficients and with higher numerical resolution than in the previous simulations (Todo et al 2010 Nucl. Fusion 50 084016). The investigation is focused on a TAE mode with toroidal mode number n = 4. It is demonstrated that the mechanism of mode saturation involves generation of zonal (n = 0) and higher-n (n ⩾ 8) sidebands, and that the sidebands effectively increase the mode damping rate via continuum damping. The n = 0 sideband includes the zonal flow peaks at the TAE gap locations. It is also found that the n = 0 poloidal flow represents a balance between the nonlinear driving force from the n = 4 components and the equilibrium plasma response to the n = 0 fluctuations. The spatial profile of the n = 8 sideband peaks at the n = 8 Alfvén continuum, indicating enhanced dissipation due to continuum damping.

A hybrid simulation code for nonlinear magnetohydrodynamics (MHD) and energetic-particle dynamics has been extended to simulate recurrent bursts of Alfvén eigenmodes by implementing the energetic-particle source, collisions and losses. The Alfvén eigenmode bursts with synchronization of multiple modes and beam ion losses at each burst are successfully simulated with nonlinear MHD effects for the physics condition similar to a reduced simulation for a TFTR experiment (Wong et al 1991 Phys. Rev. Lett. 66 1874, Todo et al 2003 Phys. Plasmas 10 2888). It is demonstrated with a comparison between nonlinear MHD and linear MHD simulation results that the nonlinear MHD effects significantly reduce both the saturation amplitude of the Alfvén eigenmodes and the beam ion losses. Two types of time evolution are found depending on the MHD dissipation coefficients, namely viscosity, resistivity and diffusivity. The Alfvén eigenmode bursts take place for higher dissipation coefficients with roughly 10% drop in stored beam energy and the maximum amplitude of the dominant magnetic fluctuation harmonic δB m/n /B ∼ 5 × 10−3 at the mode peak location inside the plasma. Quadratic dependence of beam ion loss rate on magnetic fluctuation amplitude is found for the bursting evolution in the nonlinear MHD simulation. For lower dissipation coefficients, the amplitude of the Alfvén eigenmodes is at steady levels δB m/n /B ∼ 2 × 10−3 and the beam ion losses take place continuously. The beam ion pressure profiles are similar among the different dissipation coefficients, and the stored beam energy is higher for higher dissipation coefficients.

Nonlinear magnetohydrodynamic (MHD) effects on Alfvén eigenmode evolution were investigated via hybrid simulations of an MHD fluid interacting with energetic particles. The investigation focused on the evolution of an n = 4 toroidal Alfvén eigenmode (TAE) which is destabilized by energetic particles in a tokamak. In addition to fully nonlinear code, a linear-MHD code was used for comparison. The only nonlinearity in that linear code is from the energetic-particle dynamics. No significant difference was found in the results of the two codes for low saturation levels, δB/B ∼ 10−3. In contrast, when the TAE saturation level predicted by the linear code is δB/B ∼ 10−2, the saturation amplitude in the fully nonlinear simulation was reduced by a factor of 2 due to the generation of zonal (n = 0) and higher-n (n ⩾ 8) modes. This reduction is attributed to the increased dissipation arising from the nonlinearly generated modes. The fully nonlinear simulations also show that geodesic acoustic mode is excited by the MHD nonlinearity after the TAE mode saturation.

Remarkable progress in the physical parameters of net-current free plasmas has been made in the Large Helical Device (LHD) since the last Fusion Energy Conference in Chengdu, 2006 (Motojima et al 2007 Nucl. Fusion 47 S668). The beta value reached 5% and a high-beta state beyond 4.5% from the diamagnetic measurement has been maintained for longer than 100 times the energy confinement time. The density and temperature regimes have also been extended. The central density has exceeded 1 × 1021 m−3 due to the formation of an internal diffusion barrier. The ion temperature has reached 5.2 keV at the density of 1.6 × 1019 m−3, which is associated with the suppression of ion heat conduction loss. Although these parameters have been obtained in separated discharges, each fusion-reactor relevant parameter has elucidated the potential of net-current free heliotron plasmas. Diversified studies in recent LHD experiments are reviewed in this paper.

Bursting energetic particle mode (EPM) and fast ion loss induced by the EPM have been observed and experimentally investigated in the Compact Helical System (CHS). A directional probe method was newly applied to neutral beam heated plasmas for fast ion measurement, and it successfully revealed fast ion behaviour inside the last closed flux surface (LCFS) during bursting EPMs. A convective oscillation was observed only inside the LCFS, and drives little loss of fast ions. An anomalous loss of fast ions produced by the phase difference between the fast ion oscillation and the EPM with m = 3/n = 2 was observed in the outward shifted configurations. Another anomalous loss of fast ions due to EPMs with m = 2/n = 1 was also observed in the inward shifted configuration. The characteristics of the fast-ion-loss flux are different from those observed in the outward shifted configurations, and the loss mechanism seems to be related to a mode mixture between the EPM and a low frequency mode.