Ion Internal Transport Barrier Research Articles

Identification of I-mode with ion ITB in NBI-heated plasmas on the HL-2A tokamak

Improved energy confinement mode (I-mode) has been observed on the HL-2A tokamak. The I-mode features an edge transport barrier (ETB) in electron temperature and a low confinement mode like edge density. A weakly coherent mode (WCM) is observed in the edge region with the frequency of kHz and the poloidal wavenumber of . The maximum of WCM amplitude is located near the top of electron temperature pedestal. A critical value of velocity shear for triggering the L–I transition has been found, and is much lower than that for triggering the L–H transition. Additionally, ion internal transport barrier (ITB) has been observed in the I-mode. The formation of ion ITB is due to the increase of velocity shear, leading to the suppression of turbulence. Transport analysis further confirms the existence of electron ETB and ion ITB. The coexistence of electron ETB and ion ITB leads to an improved plasma confinement, which is comparable to that in the high confinement mode, suggesting that I-mode with ITB regimes could be an interesting operation scenario for future fusion devices.

Explaining the lack of power degradation of energy confinement in wide pedestal quiescent H-modes via transport modeling

Wide pedestal quiescent H (WPQH)-mode is an attractive scenario for future burning plasmas as they operate without ELMs. WPQH is characterized by formation of a wider and higher pedestal (than quiescent H-mode), and broadband fluctuations in the pedestal. Unlike conventional H-modes, where the energy confinement time reduces with increasing heating power, the WPQH plasmas reported in this paper do not show power degradation of the energy confinement. As the injected neutral beam power was increased, reduced core (ρ ⩽ 0.45) transport calculated by transp, as well as increased core temperatures, pressure gradient and diamagnetic E × B shear rate were observed. The reduction in the heat transport and rapid decrease in the ion temperature gradient scale length suggest the formation of an ion internal transport barrier (ITB) that was accompanied by increased stored energy in the core. Quasilinear turbulent transport modeling using the trapped gyro Landau fluid (tglf) code was used to predict the ITB and its turbulence stability properties. By using profiles and equilibria produced by matching the transp transport fluxes with the tglf transport model within the tgyro transport solver, the energy confinement time captures the experimentally observed insensitivity to the increased P NBI. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by E × B shear, T i/T e ratio and Shafranov shift; the latter was found to have the strongest effect on the turbulence suppression at the highest heating level.

Predictive multi-channel integrated modeling of a reversed magnetic shear H-mode discharge with internal transport barrier in EAST

A reversed magnetic shear high plasma performance H-mode discharge with internal transport barrier (ITB) has been realized and modeled on experimental advanced superconducting tokamak (EAST) in deuterium (D) plasmas with an ITER-like tungsten (W) upper divertor. In this shot (#71326), the dimensionless parameter G () can reach to 0.25 (H 89 ⩽ 2, β N ⩽ 2, q 95 ∼ 4, I p ⩽ 450 kA, B T ⩽ 1.6 T). The ITB has been observed in the channels of particle (electron density n e), electron and ion temperature (T e and T i) and toroidal rotation velocity (V t). However, in the formation process of V t and n e ITBs, the heavy impurities accumulate in the core plasmas, which may be one of the most important reasons for the limitation of the improvement of plasma performance. A time slice (t = 4.65 s) of the high plasma performance phase with n e, T i and V t ITBs has been analyzed. It is found that the position of minimum safety factor (q min) is at about normalized radius ρ = 0.4 where may be the location of ion ITB foot, during equilibrium reconstruction from EFIT using external magnetic and internal polarimeter-INTferometer (POINT) measurement constraints. Based on this equilibrium, the source is self-consistently calculated by ONETWO and NUBEAM. Then the n e, T e, T i and V t have been modeled predictively and simultaneously by TGYRO utilizing the reconstructed equilibrium, the self-consistent source and the experimental radiated power for the first time in EAST. This modeling well reproduces the experimental n e and T i profiles in the reversed magnetic shear plasmas. In addition, it is shown that reversed magnetic shear rather than shear of radial electric field (E r) mainly produced by V t plays a key role to the formation of n e and T i ITBs. The reasons for central accumulation of W impurity are also clarified by computing both neoclassical and turbulent transport components of W. It turns out that the neoclassical transport dominates over turbulent transport for W and neoclassical pinch of W causes its central accumulation.

Characterization of isotope effect on ion internal transport barrier and its parameter dependence in the Large Helical Device

In this paper, the background physics of the isotope effects in the ion internal transport barrier (ITB) are discussed in detail. An heuristic criterion for the ITB strength is defined based on the nonlinear dependence of the ion thermal diffusivity on the local ion temperature in the L-mode phase. Comparing deuterium plasmas and hydrogen plasmas, two isotope effects on the ion ITB are clarified: stronger ITBs formed in the deuterium plasmas and an ITB concomitant edge confinement degradation in the hydrogen plasmas. Principal component analysis reveals that the ion ITB becomes strong when a high input power normalized by the line averaged electron density is applied and electron density profile is peaked. A gyrokinetic simulation suggests that the ITB profile is determined by the ion temperature gradient driven turbulence, while the way the profile saturates in L-mode plasmas is unknown. In the electron density turbulence behavior, a branch transition is observed, where the increasing trend in turbulence amplitude against the ITB strength is flipped to a decreasing trend across the ITB formation. The radial electric field structure is measured by the charge exchange recombination spectroscopy system. It is found that the radial electric field shear plays a minor role in determining the ITB strength.

Isotope effects on transport in LHD

Isotope effects are one of the most important issues for predicting future reactor operations. Large helical device (LHD) is the presently working largest stellarator/helical device using super conducting helical coils. In LHD, deuterium experiments started in 2017. Extensive studies regarding isotope effects on transport have been carried out. In this paper, the results of isotope effect studies in LHD are reported. The systematic studies were performed adjusting operational parameters and nondimensional parameters. In L mode like normal confinement plasma, where internal and edge transport barriers are not formed, the scaling of global energy confinement time (τ E) with operational parameters shows positive mass dependence (M 0.27; where M is effective ion mass) in electron cyclotron heating plasma and no mass dependence (M 0.0) in neutral beam injection heating plasma. The non-negative ion mass dependence is anti-gyro-Bohm scaling. The role of the turbulence in isotope effects was also found by turbulence measurements and gyrokinetic simulation. Better accessibility to electron and ion internal transport barrier (ITB) plasma is found in deuterium (D) plasma than in hydrogen (H). Gyro kinetic non-linear simulation shows reduced ion heat flux due to the larger generation of zonal flow in deuterium plasma. Peaked carbon density profile plays a prominent role in reducing ion energy transport in ITB plasma. This is evident only in plasma with deuterium ions. New findings on the mixing and non-mixing states of D and H particle transports are reported. In the mixing state, ion particle diffusivities are higher than electron particle diffusivities and D and H ion density profiles are almost identical. In the non-mixing state, ion particle diffusivity is much lower than electron diffusivity. Deuterium and hydrogen ion profiles are clearly different. Different turbulence structures were found in the mixing and non-mixing states suggesting different turbulence modes play a role.

Progress of Experimental Studies in the HL-2A Tokamak

During the last several years, the HL-2A experiment has made significant progress in the following areas: (1) lower-hybrid wave (LHW) heating and current drive, (2) plasma confinement and turbulent transport, (3) magnetohydrodynamic (MHD) instabilities and energetic particle physics and (4) H-mode and edge localized mode (ELM) control. The results show that the LHW system working in the co-current mode can reach higher driving efficiency and full non-inductive lower-hybrid current drive (LHCD) has been achieved. The intrinsic poloidal torque characterized by the divergence of the residual stress is deduced from synthesis for the first time. The dynamics of spectral symmetry breaking in drift wave turbulence is in good agreement with the development of the poloidal torque to drive the edge poloidal flow. The influence of the cross-phase dynamics on turbulent stress was also investigated. The ion internal transport barrier has been observed in the NBI-heated plasma, and inside the barrier the ion thermal transport is reduced to the neoclassical level. Besides, micro-turbulence is modulated by the rotation frequency of the magnetic island, and this modulation effect is related to a critical island width. Strong E × B shear is found at the island boundary. Three kinds of axisymmetric modes, beta-induced Alfven eigenmode (BAE), toroidal Alfven eigenmode (TAE) and the ellipticity-induced Alfven eigenmode (EAE), are found to be driven unstable by nonlinear mode coupling between Alfven eigenmodes and tearing mode which is well explained by the nonlinear gyrokinetic theory. The fishbone and tearing modes were actively controlled by the electron cyclotron resonance heating (ECRH). The dynamics of the edge plasma flows and turbulence during the L–I–H transition have been dedicatedly investigated. The geodesic acoustic mode (GAM) and limit cycle oscillation (LCO) coexist for a short time and disappear in the H-mode plasma with the increasing of E × B shear flow before the I–H transition, which plays an important role in the turbulence suppression. Different techniques, such as LHW, ECRH, resonant magnetic perturbation (RMP), and impurity seeding by the laser blow-off (LBO) and supersonic molecular beam injection (SMBI), have been successfully applied to control the large ELMs. It has been found that pedestal turbulence enhancement might be responsible for the observed mitigation effect.

Transport characteristics of deuterium and hydrogen plasmas with ion internal transport barrier in the Large Helical Device

A remarkable extension of the high-ion-temperature (high-Ti) regime was obtained in deuterium plasma experiments in the Large Helical Device. In order to clarify transport characteristics in the ion internal transport barrier (ITB) formation with an isotope effect, a dataset of pure deuterium (nD/ne > 0.8) and pure hydrogen (nH/ne > 0.8) plasmas in the high-Ti regime were analyzed, and two mechanisms of transport improvement were characterized. A significant reduction of ion heat transport in the core of both deuterium and hydrogen plasmas was observed, indicating ion ITB formation. The dependence of the ion heat diffusivity on temperature ratio (Te/Ti) and normalized Ti-gradient (R/LTi = −(R/Ti)(dTi/dr)) was investigated in the core region, in which gyrokinetic simulations with the GKV code predict the destabilization of ion temperature gradient (ITG) modes. The Te/Ti dependence shows ITG-like property, while a significant deviation from the ITG-like property is found in the R/ dependence, indicating suppression of the ITG mode in the large R/ regime and resultant ion ITB formation. In the comparison between deuterium plasma and hydrogen plasma, the lower transport in the deuterium plasma is observed in both ion and electron heat diffusivities, indicating a significant isotope effect. It was found with the nonlinear turbulent transport simulation with GKV that the zonal flow enhancement contributes to the ITG suppression in the deuterium plasma.

The isotope effect on impurities and bulk ion particle transport in the Large Helical Device

The isotope effect on impurities and bulk ion particle transport is investigated by using the deuterium, hydrogen, and isotope mixture plasma in the Large Helical Device (LHD). A clear isotope effect is observed in the impurity transport but not the bulk ion transport. The isotope effects on impurity transport and ion heat transport are observed as a primary and a secondary effect, respectively, in the plasma with an internal transport barrier (ITB). In the LHD, an ion ITB is always transient because the impurity hole triggered by the increase of ion temperature gradient causes the enhancement of ion heat transport and gradually terminates the ion ITB. The formation of an impurity hole becomes slower in the deuterium (D) plasma than the hydrogen (H) plasma. This primary isotope effect on impurity transport contributes the longer sustainment of the ion ITB state because the low ion thermal diffusivity can be sustained as long as the normalized carbon impurity gradient R/Ln,c, where , is above the critical value (∼−5). Therefore, the longer sustainment of the ITB state in the deuterium plasma is considered to be a secondary isotope effect due to the mitigation of the impurity hole. The radial profile of H and D ion density is measured using bulk charge exchange spectroscopy inside the isotope mixture plasma. The decay time of H ion density after the H-pellet injection and the decay time of D ion density after D-pellet injection are almost identical, which demonstrates that there is no significant isotope effect on ion particle transport.

Realization of high Ti plasmas and confinement characteristics of ITB plasmas in the LHD deuterium experiments

The deuterium (D) operation was initiated in the LHD in 2017. In the first campaign of the D experiments, we successfully extended the high temperature regime in the LHD. The new record of the ion temperature (Ti) of 10 keV associated with the ion internal transport barrier (ITB) was achieved due to several operational optimization. The thermal confinement characteristics of ITB plasmas were compared between hydrogen and D discharges. The effective ion thermal diffusivity of the ion-ITB plasmas was found to be smaller in the D discharges than that in the H discharges. The profiles of the Ti, the electron density, and the impurity of the high Ti plasmas strongly depended on the magnetic configuration and these profiles tended to peaked in the inward-shifted configuration. It was also found that the electron thermal confinement of the electron-ITB plasmas was clearly improved in the deuterium case. The GKV simulation showed the linear growth rate of TEM/ITG reduced in the plasmas with D both for the ion ITB and the electron ITB plasmas and qualitatively agreed with the tendency of the change in the thermal diffusivity obtained from the power balance analysis.

Carbon impurities behavior and its impact on ion thermal confinement in high-ion-temperature deuterium discharges on the Large Helical Device

The behavior of carbon impurities in deuterium plasmas and its impact on thermal confinement were investigated in comparison with hydrogen plasmas in the Large Helical Device (LHD). Deuterium plasma experiments have been started in the LHD and high-ion-temperature plasmas with central ion temperature (T i) of 10 keV were successfully obtained. The thermal confinement improvement could be sustained for a longer time compared with hydrogen plasmas. An isotope effect was observed in the time evolution of the carbon density profiles. A transiently peaked profile was observed in the deuterium plasmas due to the smaller carbon convection velocity and diffusivity in the deuterium plasmas compared with the hydrogen plasmas. The peaked carbon density profile was strongly correlated to the ion thermal confinement improvement. The peaking of the carbon density profile will be one of the clues to clarify the unexplained mechanisms for the formations of ion internal transport barrier and impurity hole on LHD. These results could also lead to a better understanding of the isotope effect in the thermal confinement in torus plasma.

Collisionality dependence and ion species effects on heat transport in He and H plasma, and the role of ion scale turbulence in LHD

Surveys of the ion and electron heat transports of neutral beam (NB) heating plasma were carried out by power balance analysis in He and H rich plasma at LHD. Collisionality was scanned by changing density and heating power. The characteristics of the transport vary depending on collisionality. In low collisionality, with low density and high heating power, an ion internal transport barrier (ITB) was formed. The ion heat conductivity (χi) is lower than electron heat conductivity (χe) in the core region at ρ < 0.7. On the other hand, in high collisionality, with high density and low heating power, χi is higher than χe across the entire range of plasma. These different confinement regimes are associated with different fluctuation characteristics. In ion ITB, fluctuation has a peak at ρ = 0.7, and in normal confinement, fluctuation has a peak at ρ = 1.0. The two confinement modes change gradually depending on the collisionality. Scans of concentration ratio between He and H were also performed. The ion confinement improvements were investigated using gyro-Bohm normalization, taking account of the effective mass and charge. The concentration ratio affected the normalized χi only in the edge region (ρ ~ 1.0). This indicates ion species effects vary depending on collisionality. Turbulence was modulated by the fast ion loss instability. The modulation of turbulence is higher in H rich than in He rich plasma.

Extension of operational regime in high-temperature plasmas and effect of ECRH on ion thermal transport in the LHD

A simultaneous high ion temperature (Ti) and high electron temperature (Te) regime was successfully extended due to an optimized heating scenario in the LHD. Such high-temperature plasmas were realized by the simultaneous formation of an electron internal transport barrier (ITB) and an ion ITB by the combination of high power NBI and ECRH. Although the ion thermal confinement was degraded in the plasma core with an increase of Te/Ti by the on-axis ECRH, it was found that the ion thermal confinement was improved at the plasma edge. The normalized ion thermal diffusivity at the plasma edge was reduced by 70%. The improvement of the ion thermal confinement at the edge led to an increase in Ti in the entire plasma region, even though the core transport was degraded.

Overview of recent HL-2A experiments

Since the last Fusion Energy Conference, significant progress has been made in the following areas. The first high coupling efficiency low-hybrid current drive (LHCD) with a passive–active multi-junction (PAM) antenna was successfully demonstrated in the H-mode on the HL-2A tokamak. Double critical impurity gradients of electromagnetic turbulence were observed in H-mode plasmas. Various ELM mitigation techniques have been investigated, including supersonic molecular beam injection (SMBI), impurity seeding, resonant magnetic perturbation (RMP) and low-hybrid wave (LHW). The ion internal transport barrier was observed in neutral beam injection (NBI) heated plasmas. Neoclassical tearing modes (NTMs) driven by the transient perturbation of local electron temperature during non-local thermal transport events have been observed, and a new type of non-local transport triggered by the ion fishbone was found. A long-lasting runaway electron plateau was achieved after argon injection and the runaway current was successfully suppressed by SMBI. It was found that low-n Alfvénic ion temperature gradient (AITG) modes can be destabilized in ohmic plasmas, even with weak magnetic shear and low-pressure gradients. For the first time, the synchronization of geodesic acoustic mode (GAM) and magnetic fluctuations was observed in edge plasmas, revealing frequency entrainment and phase lock. The spatiotemporal features of zonal flows were also studied using multi-channel correlation Doppler reflectometers.

Ion internal transport barrier in neutral beam heated plasmas on HL-2A

Ion internal transport barriers (iITBs) are first observed in neutral beam injection (NBI) heated plasmas at the HL-2A tokamak. The position of the barrier foot, in the stationary state, coincides with the q = 1 surface within its uncertainty of measurement. iITBs can develop more easily at the beginning of NBI heating. Also, iITBs are unstable for the sawtooth plasma. Simulations reveal that the thermal diffusivity of ions (χi) inside the barrier can be as low as the neoclassical level. It is observed that the flow shear in the stationary iITB state reaches the level required for suppressing the ion temperature gradient mode instability, which indicates the important role of flow shear in sustaining the iITB.

Comparison of Ion Internal Transport Barrier Formation between Hydrogen and Helium Dominated Plasmas

Comparison of Ion Internal Transport Barrier Formation between Hydrogen and Helium Dominated Plasmas

Integrated discharge scenario for high-temperature helical plasma in LHD

The discharge scenario of high temperature plasma with a helical configuration has significantly progressed. The increase of central ion temperature due to the reduction of wall recycling was clearly observed. The peaking of the ion heating profile and the reduction of charge exchange loss of energetic ions play an important role for further improvement of ion heat transport in the ion internal transport barrier (ITB) core. The ion ITB and electron ITB have been successfully integrated due to the superposition of centrally focused electron cyclotron heating to the ion ITB plasma, and the high temperature regime of the ion temperature comparable to the electron temperature (Ti ~ Te) has been significantly extended. The width of the ion ITB formed with electron ITB is wider than the width of electron ITB. The positive radial electric field was observed in the integrated ITB plasma by a heavy ion beam probe, while the negative radial electric field was observed in ion ITB plasmas. The ion temperature gradient decreases with the increase of the temperature ratio (Te /Ti).

Overview of transport and MHD stability study: focusing on the impact of magnetic field topology in the Large Helical Device

The progress in the understanding of the physics and the concurrent parameter extension in the large helical device since the last IAEA-FEC, in 2012 (Kaneko O et al 2013 Nucl. Fusion 53 095024), is reviewed. Plasma with high ion and electron temperatures (Ti(0) ∼ Te(0) ∼ 6 keV) with simultaneous ion and electron internal transport barriers is obtained by controlling recycling and heating deposition. A sign flip of the nondiffusive term of impurity/momentum transport (residual stress and convection flow) is observed, which is associated with the formation of a transport barrier. The impact of the topology of three-dimensional magnetic fields (stochastic magnetic fields and magnetic islands) on heat momentum, particle/impurity transport and magnetohydrodynamic stability is also discussed. In the steady state operation, a 48 min discharge with a line-averaged electron density of 1 × 1019 m−3 and with high electron and ion temperatures (Ti(0) ∼ Te(0) ∼ 2 keV), resulting in 3.36 GJ of input energy, is achieved.

Quasilinear carbon transport in an impurity hole plasma in LHD

Comprehensive electrostatic gyrokinetic linear stability calculations for ion-scale microinstabilities in an LHD plasma with an ion-internal transport barrier (ITB) and carbon “impurity hole” are used to make quasilinear estimates of particle flux to explore whether microturbulence can explain the observed outward carbon fluxes that flow “up” the impurity density gradient. The ion temperature is not stationary in the ion-ITB phase of the simulated discharge, during which the core carbon density decreases continuously. To fully sample these varying conditions, the calculations are carried out at three radial locations and four times. The plasma parameter inputs are based on experimentally measured profiles of electron and ion temperature, as well as electron and carbon density. The spectroscopic line-average ratio of hydrogen and helium densities is used to set the density of these species. Three ion species (H,He,C) and the electrons are treated kinetically, including collisions. Electron instability drive does enhance the growth rate significantly, but the most unstable modes have characteristics of ion temperature gradient modes in all cases. As the carbon density gradient is scanned between the measured value and zero, the quasilinear carbon flux is invariably inward when the carbon density profile is hollow, so turbulent transport due to the instabilities considered here does not explain the observed outward flux of impurities in impurity hole plasmas. The stiffness of the quasilinear ion heat flux is found to be 1.7–2.3, which is lower than several estimates in tokamaks.

Dynamic transport study of heat and momentum transport in a plasma with improved ion confinement in the Large Helical Device

Dynamic transport study taking account of the slowing-down effect on the neutral beam injection heating is applied to a high ion temperature plasma with an ion internal transport barrier (ITB) obtained by carbon pellet injection, which records the highest ion temperature of around ∼7 keV in the Large Helical Device. The transient increase in ion heating is clarified during the density decay phase just after the carbon pellet injection by considering the slowing-down effect. The dynamic transport analysis also includes the change in heat flux due to the change in kinetic energy inside the plasma with the time scale of the global energy confinement time, which is required to investigate the heat and momentum transport during the transient phase more exactly. The characteristics of the ion heat and momentum transport improvement during the ion ITB formation phase are clarified by the dynamic transport study.

Role of external torque in the formation of ion thermal internal transport barriers

We present an analytic study of the impact of external torque on the formation of ion internal transport barriers (ITBs). A simple analytic relation representing the effect of low external torque on transport bifurcations is derived based on a two field transport model of pressure and toroidal momentum density. It is found that the application of an external torque can either facilitate or hamper bifurcation in heat flux driven plasmas depending on its sign relative to the direction of intrinsic torque. The ratio between radially integrated momentum (i.e., external torque) density to power input is shown to be a key macroscopic control parameter governing the characteristics of bifurcation.