Turbulence Transition in Magnetically Confined Hydrogen and Deuterium Plasmas.

In this study, we discovered a turbulence transition in a large helical device. The turbulence level and turbulence-driven energy transport decrease to a specific transition density and increase above it. The ruling turbulences below and above the transition density were ion-temperature gradient (ITG) and resistive-interchange (RI) turbulences, consistent with the predictions of gyrokinetic theory and two-fluid MHD model, respectively. Isotope experiments on hydrogen (H) and deuterium (D) clarified the role of transitions. In the ITG regime, turbulence levels and energy transport were comparable in the H and D plasmas. In contrast, in the RI regime, they were clearly suppressed in the D plasma. The results provide crucial knowledge for understanding isotope effects and future optimization of stellarator and heliotron devices.

Impact of Magnetic Field Configuration on Heat Transport in Stellarators and Heliotrons.

We assess the magnetic field configuration in modern fusion devices by comparing experiments with the same heating power, between a stellarator and a heliotron. The key role of turbulence is evident in the optimized stellarator, while neoclassical processes largely determine the transport in the heliotron device. Gyrokinetic simulations elucidate the underlying mechanisms promoting stronger ion scale turbulence in the stellarator. Similar plasma performances in these experiments suggests that neoclassical and turbulent transport should both be optimized in next step reactor designs.

Transition between Isotope-Mixing and Nonmixing States in Hydrogen-Deuterium Mixture Plasmas.

The transition between isotope-mixing and nonmixing states in hydrogen-deuterium mixture plasmas is observed in the isotope (hydrogen and deuterium) mixture plasma in the Large Helical Device. In the nonmixing state, the isotope density ratio profile is nonuniform when the beam fueling isotope species differs from the recycling isotope species and the profile varies significantly depending on the ratio of the recycling isotope species, although the electron density profile shape is unchanged. The fast transition from nonmixing state to isotope-mixing state (nearly uniform profile of isotope ion density ratio) is observed associated with the change of electron density profile from peaked to hollow profile by the pellet injection near the plasma periphery. The transition from nonmixing to isotope-mixing state strongly correlates with the increase of turbulence measurements and the transition of turbulence state from TEM to ion temperature gradient is predicted by gyrokinetic simulation.

Isotope Effect on Energy Confinement Time and Thermal Transport in Neutral-Beam-Heated Stellarator-Heliotron Plasmas.

The isotope effect on energy confinement time and thermal transport has been investigated for plasmas confined by a stellarator-heliotron magnetic field. This is the first detailed assessment of an isotope effect in a stellarator heliotron. Hydrogen and deuterium plasmas heated by neutral beam injection on the Large Helical Device have exhibited no significant dependence on the isotope mass in thermal energy confinement time, which is not consistent with the simple gyro-Bohm model. A comparison of thermal diffusivity for dimensionally similar hydrogen and deuterium plasmas in terms of the gyroradius, collisionality, and thermal pressure has clearly shown robust confinement improvement in deuterium to compensate for the unfavorable mass dependence predicted by the gyro-Bohm model.

Suppression of Trapped Energetic Ions Driven Resistive Interchange Modes with Electron Cyclotron Heating in a Helical Plasma.

The resistive interchange mode destabilized by the resonant interaction with the trapped energetic ions is fully suppressed when the injected power of electron cyclotron heating exceeds a certain threshold. It is shown for the first time that the complete stabilization of the energetic-particle-driven mode without relaxing the energetic particle (EP) pressure gradient is possible by reducing the radial width of the eigenmodes δ_{w}, especially when δ_{w} narrows to a small enough value relative to the finite orbit width of EP.

Resistive interchange modes destabilized by helically trapped energetic ions in a helical plasma.

A new bursting m=1/n=1 instability (m,n: poloidal and toroidal mode numbers) with rapid frequency chirping down has been observed for the first time in a helical plasma with intense perpendicular neutral beam injection. This is destabilized in the plasma peripheral region by resonant interaction between helically trapped energetic ions and the resistive interchange mode. A large radial electric field is induced near the edge due to enhanced radial transport of the trapped energetic ions by the mode, and leads to clear change in toroidal plasma flow, suppression of microturbulence, and triggering an improvement of bulk plasma confinement.

Charge exchange spectroscopy using spatial heterodyne spectrometer in the large helical device.

In this study, the use of a spatial heterodyne spectrometer (SHS) to measure the toroidal flow velocity (Vf) and the ion temperature (TC6+ ) of the C6+ impurity ion by charge exchange spectroscopy was explored. The instrumental width (IW) of the SHS (aperture size = 16.77 mm2, etendue = 2.9867 mm2sr) was extrapolated to be 0.09 nm, which is half of the 0.17 nm IW extrapolated for a conventionally used dispersive spectrometer (DS) (aperture size = 2.6 mm2, etendue = 0.2605 mm2sr). The resulting Vf and TC6+ measurements were found to be in good agreement with those measured using the DS.

Receiver circuit improvement of dual frequency-comb ka-band Doppler backscattering system in the large helical device (LHD).

Doppler-backscattering (DBS) has been used in several fusion plasma devices because it can measure the perpendicular velocity of electron density perturbation v⊥, the radial electric field Er, and the perpendicular wavenumber spectrum S(k⊥) with high wavenumber and spatial resolution. In particular, recently constructed frequency comb DBS systems enable observation of turbulent phenomena at multiple observation points in the radial direction. A dual-comb microwave DBS system has been developed for the large helical device plasma measurement. Since it is desirable to control the gain of each frequency-comb separately, a frequency-comb DBS system was developed with a function to adjust the gain of the scattered signal intensity of each channel separately. A correction processing method was also developed to correct the amplitude ratio and the phase difference between the in-phase and quadrature-phase signals of the scattered signals. As a result, the error in Doppler-shift estimation required to observe vertical velocity and the radial electric field was reduced, which enables more precise measurements.

Preceding propagation of turbulence pulses at avalanche events in a magnetically confined plasma.

The preceding propagation of turbulence pulses has been observed for the first time in heat avalanche events during the collapse of the electron internal transport barrier (e-ITB) in the Large Helical Device. The turbulence and heat pulses are generated near the foot of the e-ITB and propagate to the peripheral region within a much shorter time than the diffusion timescale. The propagation speed of the turbulence pulse is approximately 10 km/s, which is faster than that of the heat pulse propagating at a speed of 1.5 km/s. The heat pulse propagates at approximately the same speed as that in the theoretical prediction, whereas the turbulence pulse propagates one order of magnitude faster than that in the prediction, thereby providing important insights into the physics of non-local transport.

Measurements of radial profile of isotope density ratio using bulk charge exchange spectroscopy.

A bulk charge exchange spectroscopy (BCXS) system using a grism (grating prism) spectrometer has been applied to measure the profile of the deuterium (D) fraction in deuterium and hydrogen (H) mixture plasma in the Large Helical Device. The observed spectrum can be fitted with four Gaussian functions successfully by reduction of free parameters for the least-squares fit. The plasma flow velocity and ion temperature profile measured by charge exchange spectroscopy using carbon impurity are used for estimation of the wavelength shift of hot components to reduce the free parameter. The ion temperature is used to estimate the apparent wavelength shift due to the energy dependent emission cross section only and is not used to set the Doppler width for H and D in the fitting. The sensitivity of the evaluated D fraction on the velocity is increased for a higher D fraction. The error of the D fraction is calculated from the error in the fitted parameter and sensitivity on the velocity of the hot component. The difference in the profile and time trace of the D fraction with D pellet and H pellet injection was observed clearly by BCXS using a grism spectrometer.

Analysis of the Motional Stark Effect (MSE) diagnostic to measure the rotational transform and current profile in the Large Helical Device.

The analysis method of the Motional Stark Effect (MSE) diagnostic to measure the rotational transform and current profiles in the Large Helical Device has been improved. This was done by using the Variational Moments Equilibrium Code to calculate an equilibrium database for various pressure profiles and current profiles. This method looks for the radial profile of the rotational transform in the equilibrium database that gives the best fit to the polarization angle profiles measured with the MSE diagnostic. This analysis improves the measurements of rotational transform, especially near the magnetic axis, where the sensitivity of the polarization angle measurements becomes low and the uncertainty due to error in the estimation of the Pfirsch-Schlüter current becomes large. The radial profiles of the rotational transform and current profiles for Electron Cyclotron Current Drive and Neutral Beam Current Drive are obtained in the new analysis method with a sufficiently high accuracy to discuss the discrepancy of the current density profiles between the measurements and the calculations.

A new multi-tracer pellet injection for a simultaneous study of low- and mid/high-Z impurities in high-temperature plasmas.

A new multi-tracer technique in the Tracer-Encapsulated Solid Pellet (TESPEL) method has been developed in order to acquire simultaneously the information about the behaviors of various impurities, i.e., to study concurrently the behaviors of low- and mid/high-Z impurities in magnetically confined high-temperature plasmas. In this new technique, an inorganic compound (for example, lithium titanate, Li2TiO3) is proposed to be used as a tracer embedded in the core of the TESPEL, instead of pure elements. The results of the proof-of-principle experiment clearly demonstrate the applicability of the new multi-tracer technique in the TESPEL method for the simultaneous study of behaviors of low- and mid/high-Z impurities in high-temperature plasmas.

W-band millimeter-wave back-scattering system for high wavenumber turbulence measurements in LHD.

A 90 GHz W-band millimeter-wave back-scattering system is designed and installed for measuring electron scale turbulence (k⊥ρs ∼ 40). A metal lens relay antenna is used for in-vessel beam focusing, and a beam diameter of less than 40 mm is achieved in the plasma core region. This antenna can be steered at an angle of 159° ± 6°, which almost covers the plasma radius. The estimated size of the scattering volume is ∼105 mm at the edge and 135 mm at the core, respectively. A 60 m corrugated waveguide is used to achieve a low transmission loss of ∼8 dB. A heterodyne detection system for millimeter-wave circuits with probing power modulation can distinguish the scattered signal from background noise.

Observation of reversed-shear Alfvén eigenmodes excited by energetic ions in a helical plasma.

Reversed-shear Alfvén eigenmodes were observed for the first time in a helical plasma having negative q0'' (the curvature of the safety factor q at the zero shear layer). The frequency is swept downward and upward sequentially via the time variation in the maximum of q. The eigenmodes calculated by ideal MHD theory are consistent with the experimental data. The frequency sweeping is mainly determined by the effects of energetic ions and the bulk pressure gradient. Coupling of reversed-shear Alfvén eigenmodes with energetic ion driven geodesic acoustic modes generates a multitude of frequency-sweeping modes.

Isotope effects in self-organization of internal transport barrier and concomitant edge confinement degradation in steady-state LHD plasmas.

The isotope effect, which has been a long-standing mystery in the turbulent magnetically confined plasmas, is the phenomena that the plasma generated with heavier hydrogen isotope show a mitigated transport. This is on the contrary to what is predicted with the simple scaling theory, in which the heavier ions easily diffuse because of its larger gyro-radius. Thanks to the newly developed analysis method and a comprehensive parameter scan experiment in the steady-state plasmas in the Large Helical Device (LHD), the isotope effect was clearly observed in the self-organized internal transport barrier (ITB) structure for the first time. Comparing the ITB intensity in deuterium (D) and hydrogen (H) plasmas, two distinct hydrogen isotope effects are found: stronger ITB is formed in D plasmas and a significant edge confinement degradation accompanied by the ITB formation emerges in H plasmas. This observation sheds light on a new aspect of the turbulent plasmas regarding how the basic properties of the fluid material affect the turbulent structure formation in the open-system.

Measurements of radial profile of hydrogen and deuterium density in isotope mixture plasmas using bulk charge exchange spectroscopy.

A bulk charge exchange spectroscopy system has been applied to measure the radial profiles of the hydrogen (H) and deuterium (D) density ratio in the isotope mixture plasma in a large helical device. Charge exchange lines of Hα and Dα are fitted by 4 Gaussian of H and D cold components and H and D hot components with 5 parameters by combining the measurement of plasma toroidal rotation velocity with carbon charge exchange spectroscopy. The radial profiles of the relative density of hydrogen and deuterium ions are derived from H and D hot components measured and the beam density calculated from beam attenuation calculation. A proof-of-principle experiment is performed by the H pellet and the D pellet injections into the H-D mixture plasma.

Measurement of derivative of ion temperature using high spatial resolution charge exchange spectroscopy with space modulation optics.

A new technique to measure the first and second derivatives of the ion temperature profile has been developed by using a charge exchange spectroscopy system with space modulation optics. The space observed is scanned up to +/-3 cm with a cosine wave modulation frequency up to 30 Hz by shifting the object lens in front of the optical fiber bundle by 0.5 mm with a piezoelement. The first and second derivatives of ion temperature are derived from the modulation component of the ion temperature measured by using Fourier series expansion.

Trigger mechanism for the abrupt loss of energetic ions in magnetically confined plasmas.

Interaction between a quasi-stable stationary MHD mode and a tongue-shaped deformation is observed in the toroidal plasma with energetic particle driven MHD bursts. The quasi-stable stationary 1/1 MHD mode with interchange parity appears near the resonant rational surface of q = 1 between MHD bursts. The tongue-shaped deformation rapidly appears at the non-resonant non-rational surface as a localized large plasma displacement and then collapses (tongue event). It curbs the stationary 1/1 MHD mode and then triggers the collapse of energetic particle and magnetic field reconnection. The rotating 1/1 MHD mode with tearing parity at the q = 1 resonant surface, namely, the MHD burst, is excited after the tongue event.

A motional Stark effect diagnostic analysis routine for improved resolution of iota in the core of the large helical device.

A new Motional Stark Effect (MSE) analysis routine has been developed for improved spatial resolution in the core of the Large Helical Device (LHD). The routine was developed to reduce the dependency of the analysis on the Pfirsch-Schlüter (PS) current in the core. The technique used the change in the polarization angle as a function of flux in order to find the value of diota/dflux at each measurement location. By integrating inwards from the edge, the iota profile can be recovered from this method. This reduces the results' dependency on the PS current because the effect of the PS current on the MSE measurement is almost constant as a function of flux in the core; therefore, the uncertainty in the PS current has a minimal effect on the calculation of the iota profile. In addition, the VMEC database was remapped from flux into r/a space by interpolating in mode space in order to improve the database core resolution. These changes resulted in a much smoother iota profile, conforming more to the physics expectations of standard discharge scenarios in the core of the LHD.

Novel analysis technique for measuring edge density fluctuation profiles with reflectometry in the Large Helical Device.

A new method for measuring density fluctuation profiles near the edge of plasmas in the Large Helical Device (LHD) has been developed utilizing reflectometry combined with pellet-induced fast density scans. Reflectometer cutoff location was calculated by proportionally scaling the cutoff location calculated with fast far infrared laser interferometer (FIR) density profiles to match the slower time resolution results of the ray-tracing code LHD-GAUSS. Plasma velocity profile peaks generated with this reflectometer mapping were checked against velocity measurements made with charge exchange spectroscopy (CXS) and were found to agree within experimental uncertainty once diagnostic differences were accounted for. Measured density fluctuation profiles were found to peak strongly near the edge of the plasma, as is the case in most tokamaks. These measurements can be used in the future to inform inversion methods of phase contrast imaging (PCI) measurements. This result was confirmed with both a fixed frequency reflectometer and calibrated data from a multi-frequency comb reflectometer, and this method was applied successfully to a series of discharges. The full width at half maximum of the turbulence layer near the edge of the plasma was found to be only 1.5-3 cm on a series of LHD discharges, less than 5% of the normalized minor radius.