Ion Temperature Gradient Mode Research Articles

Core tungsten transport in WEST long pulse L-mode plasmas

Tungsten transport is investigated in WEST long pulse L-mode plasmas operated with the strike point on the actively cooled upper tungsten divertor. The pulses are mostly heated by lower hybrid waves. It is experimentally found that tungsten does not centrally accumulate throughout these ∼ 30 s reproducible discharges despite large normalised electron density gradients . To explain these observations, turbulent and neoclassical transport of electrons and tungsten ions are computed with GKW Peeters A.G. et al (2009 Computer Phys. Commnun. 180 2650) and NEO Belli E. and Candy J. (2008 Plasma Phys. Control. Fusion 50 095010), Belli E. and Candy J. (2012 Plasma Phys. Control. Fusion 54 015015) respectively. Additionally, interpretative integrated modelling simulations are also performed to keep data coherency despite the lack of measurements of some quantities such as the Ti profiles. Modelled are found consistent with interferometry inversions and the tungsten peaking factor remains comparable to due to dominant turbulent diffusivities inside . In the central region neoclassical W transport dominates but the convective velocities are several order of magnitudes lower compared to plasmas with toroidal rotation velocities induced by a neutral beam injection (NBI) torque. Finally, nitrogen is seeded in these pulses leading to an enhanced energy content which is consistent with stabilised ion temperature gradient modes from dilution.

On hybrid scenarios in KSTAR

We report the status of hybrid scenario experiments in Korea Superconducting Tokamak Advanced Research (KSTAR). The hybrid scenario is defined as stationary discharges with 2.4 and 2.0 at < 6.5 without or with very mild sawtooth activities in KSTAR. It is being developed towards reactor-relevant conditions. High performance of , and G-factor () 6 has been achieved and sustained for at ne/nGW ~ 0.7 with heating power of 5 MW. Some KSTAR hybrid discharges exhibit a unique feature of a slow transition from conventional H-mode to hybrid mode after the third neutral beam injection. The reason for the confinement enhancement is extensively studied in this transition period of a representative discharge exhibiting a common feature of KSTAR hybrid scenarios. 0D performance analysis with magnetohydrodynamic activities, 1D kinetic profile dynamics, power balance analysis, linear gyro-kinetic analysis and edge pedestal stability analysis were conducted. The enhancement is thought to be from both the core and the pedestal. The improvement in the core region of the ion energy channel is observed from the linear gyro-kinetic analysis considering the electromagnetic, the fast ion, the Shafranov shift, , and the magnetic shear effect. The electromagnetic finite stabilisation plays a role in the inner core region at together with the fast ion effect. The alpha stabilisation effect is also found at . , which could reduce the linear growth of the ion temperature gradient mode in the outer core region at with the highest contribution from the toroidal rotation. Regarding the improvement in the pedestal, Shafranov shift broadens the stability boundary of the pedestal in support of the diamagnetic effect. The pedestal height and width could be reproduced by the EPED model, while a realistic current profile is used to calculate the internal inductance for Shafranov shift. Based on these findings, a comprehensive confinement enhancement mechanism has been proposed by considering the core-edge interplay.

Effects of temperature gradient driven turbulence and core MHD instability on particle transport in HL-2A L-mode plasmas

Experiments in L-mode plasmas in the HL-2A tokamak show that the electron and trace Al impurity transport is related to the normalized electron temperature gradient with opposite trends. Increasing in the confinement zone (0.25 ⩽ ρ⩽ 0.5, ρ is the normalized minor radius) tends to pump out electrons but accumulate impurity, leading to slightly hollow electron density profiles meanwhile with an overall stronger inward convection of injected impurity. Moreover, after the injection of Al impurity by laser blow-off (LBO), the impurity density profiles remain deeply hollow due to the combined effect of the impurity expulsion by core magneto-hydrodynamic (MHD) activity in the plasma centre and the inward impurity convection driven by turbulence in the outer confinement region. Gyrokinetic simulation results reveal that the local change of affects the relative strength of the trapped electron mode (TEM) and the ion temperature gradient (ITG) mode, eventually determining the electron particle transport in a way consistent with the experimental observations. Comparison of simulations with and without collisionallity shows that the plasma collisionality can significantly reduce the growth rate of the TEM and the associated outward particle flux.

Finite β effects on short wavelength ion temperature gradient modes

The electromagnetic effect is studied on the short wavelength branch of the ion temperature gradient mode in the linear regime for the first time using a global gyrokinetic model. The short wavelength ion temperature gradient mode growth rate is found to be reduced in the presence of electromagnetic perturbations at finite plasma β. The effect on real frequency is found to be weak. The threshold value of ηi is found to increase for the mode as the magnitude of β is increased. The global mode structure of the short wavelength branch of the ion temperature gradient mode is compared with the conventional branch. The magnetic character of the mode, measured as the ratio of mode average square values of electromagnetic potential to electrostatic potential, is found to increase with increasing values of the plasma β. The mixing length estimate for flux shows that the maximum contribution still comes from the long wavelengths modes. The magnitude of the flux decreases with increasing β.

Ion Temperature Gradient Mode–Driven Solitary and Shock Waves in Electron-Positron-Ion Magnetized Plasma

The characteristics of linear and nonlinear waves of ion temperature gradient mode in nonuniform electron-positron-ion magnetized plasma are investigated. A modified linear dispersion relation and its reduction in various limits are discussed. In nonlinear regime, ion temperature gradient mode–driven solitary and shock wave structures are obtained and discussed. It is found that corresponding peak potential and width of these nonlinear structures strongly depend on the positron number density and temperature. Our study reveals that in electron-positron-ion plasma, the ion energy transportation in laboratory as well as in astrophysical plasmas can be affected because of positrons.

Effects of trapped electrons on the ion temperature gradient mode in tokamak plasmas with hollow density profiles

The ion temperature gradient (ITG) mode in the presence of impurity ions and trapped electrons (TEs) is numerically investigated in tokamak plasmas with hollow density profiles, using the gyrokinetic integral eigenmode equation. It is found that in the inverted density plasma, the increase of the ITG enhances the growth rate and frequency of the ITG, and the density gradient plays an important role in the ITG modes. For weak density gradient situations, the trapped electron effects increase the instability of the ITG, while the impurity has an obviously stabilizing effect. For the strong density gradient cases, both the impurities and trapped electrons enhance the ITG instabilities. In addition, it is shown that the growth rate of the ITG decreases with positive magnetic shear s while the real frequency increases with positive magnetic shear. The growth rate of the ITG increases with negative magnetic shear s while the real frequency decreases with negative magnetic shear. The length of the calculated mode structure in the positive and negative magnetic shear intervals is also presented.

A reduced MHD model for ITG-NTM interplay

A six-field reduced-MHD model is derived for plasma dynamics. The new model describes coherently both ion temperature gradient (ITG) mode and tearing mode and includes neoclassical effects. The model allows the construction of an energylike quantity with a linear pressure contribution that is conserved except for dissipative, finite Larmor radius, and neoclassical terms. This model may be used to study the nonlinear interaction between ITG microturbulence and neoclassical tearing mode, which is responsible for large-scale magnetic islands in tokamaks, and opens the way to a coherent description of turbulent impurity transport in magnetic islands.

Nonlinear gyrokinetic analysis of linear ohmic confinement to saturated ohmic confinement transition

This work presents a nonlinear gyrokinetic analysis, which addresses one of the perennial conundrums in ohmically heated fusion plasmas. The widely observed linear ohmic confinement (LOC) to saturated ohmic confinement (SOC) transition of the energy confinement time with increasing density is successfully reproduced from nonlinear gyrokinetic simulations for the first time. Suppression of trapped electron turbulence by collisional detrapping due to increasing density is found to be responsible for the transition. The microturbulence transition from trapped electron modes to ion temperature gradient modes is found to coincide, but is not a necessary condition for the LOC–SOC transition. Damping of nonlinearly generated zonal flow by increasing collisionality with density can be responsible for the energy confinement degradation in the SOC regime.

Extended investigations of isotope effects on ECRH plasma in LHD

Isotope effects of ECRH plasma in LHD were investigated in detail. A clear difference of transport and turbulence characteristics in H and D plasmas was found in the core region, with normalized radius ρ < 0.8 in high collisionality regime. On the other hand, differences of transport and turbulence were relatively small in low collisionality regime. Power balance analysis and neoclassical calculation showed a reduction of the anomalous contribution to electron and ion transport in D plasma compared with H plasma in the high collisionality regime. In core region, density modulation experiments also showed more reduced particle diffusion in D plasma than in H plasma, in the high collisionality regime. Ion scale turbulence was clearly reduced at ρ < 0.8 in high collisionality regime in D plasma compared with H plasma. The gyrokinetic linear analyses showed that the dominant instability ρ = 0.5 and 0.8 were ion temperature gradient mode (ITG). The linear growth rate of ITG was reduced in D plasma than in H plasma in high collisionality regime. This is due to the lower normalized ITG and density gradient. More hollowed density profile in D plasma is likely to be the key control parameter. Present analyses suggest that anomalous process play a role to make hollower density profiles in D plasma rather than neoclassical process. Electron scale turbulence were also investigated from the measurements and linear gyrokinetic simulations.

Development and testing of an unstructured mesh method for whole plasma gyrokinetic simulations in realistic tokamak geometry

In this work, we have formulated and implemented a mixed unstructured mesh-based finite element–Fourier decomposition scheme for gyrokinetic simulations in realistic tokamak geometry. An efficient particle positioning (particle-triangle mapping) scheme for the charge deposition and field scattering using an intermediate grid as the search index for triangles has been implemented, and a significant speed-up by a factor of ∼30 is observed as compared with the brute force scheme for a medium-size simulation. The TRIMEG (TRIangular MEsh based Gyrokinetic) code has been developed. As an application, the ion temperature gradient (ITG) mode is simulated using the simplified gyrokinetic Vlasov-Poisson model. Our simulation and that using the ORB5 code for the DIII-D Cyclone case show reasonable agreement. As an additional application, ITG simulations using an ASDEX Upgrade equilibrium have been performed with density and temperature gradient profiles similar to the Cyclone case. Capabilities of the TRIMEG code for simulations with realistic experimental equilibria in the plasma core and in the whole plasma volume with open field lines are demonstrated.

Re-splitting [formula omitted] method for electro-magnetic gyrokinetic particle-in-cell (PIC) simulation of tokamak plasmas

A new method based on a re-splitting δf technique is proposed to mitigate the cancellation problem in electro-magnetic gyrokinetic PIC simulation. This new scheme, together with the pull-back mitigation (PBM) scheme (Mishchenko et al., 2014), have recently been implemented in the GEM code (Chen et al., 2003). These two algorithms are systematically compared with the original GEM algorithm for the simulation of the ion temperature gradient (ITG) mode, the kinetic ballooning mode (KBM), the toroidal Alfvén eigenmode (TAE) and the energetic particle mode (EPM). The three algorithms agree well for all of these modes. Moreover, both the re-splitting method and PBM can mitigate the cancellation problem and improve the computational efficiency.

The dominant micro-turbulence instabilities in the lower q95 high βp plasmas on DIII-D and predict-first extrapolation

Large-radius internal transport barriers (ITB) are the signature of high scenarios on DIII-D. Previous studies show that a large Shafranov shift, rather the shear, suppresses the turbulence and helps in the formation of the large-radius ITB. New gyrokinetic simulations suggest that the remaining micro-instabilities in lower q95 (<7.0), high ITB plasmas are drift wave instabilities, including the collisionless trapped electron mode in the core and ITB peak gradient region, the electron temperature gradient mode in the ITB peak gradient region and at the ITB foot and the ion temperature gradient mode at the ITB foot. Gyrokinetic simulation results qualitatively agree with the density fluctuation analysis from beam emission spectroscopy, which suggests the existence of a low-k ion mode at the ITB foot. The gyrokinetic simulations also predict that a larger Shafranov shift can overwhelm the driving sources for turbulence from the profile gradients at higher , leading to stronger turbulence suppression and stronger ITBs in lower q95, high plasmas.

Scaling study for positive magnetic shear ELMy H-mode plasmas in JT-60U

Scaling of the density peaking for positive magnetic shear ELMy H-mode plasmas in JT-60U has been developed. Although the density peaking generally depends on collisionality, a variation of the density peaking factor for the same collisionality exists, and the variation is different between plasmas with co-directed and ctr-directed toroidal rotation(co-VT plasmas and ctr-VT plasmas). The variation of the co-VT plasmas can be explained by particle source profile. As a result of scaling, the peaking factor of the co-VT plasmas depends on the particle source rate profile from neutral beams (NBs) as much as collisionality. The dataset of the ctr-VT plasmas has a larger variation of the peaking factor compared to that of the co-VT plasmas. The larger variation stems from that the peaking factor of the ctr-VT plasmas also depends on the normalized ion temperature gradient at the edge region. As a result of scaling, the normalized ion temperature gradient influences the peaking factor as much as or a little larger than collisionality. The parameter regime of the ctr-VT plasmas ranges from the trapped electron mode (TEM) to the ion temperature gradient mode (ITG mode). The plasmas with the positive correlation between the normalized density gradient and the normalized ion temperature gradient are dominated by the TEM. On the other hand, the plasmas with the negative correlation are dominated by the ITG mode.

Characterisation of highly radiating neon seeded plasmas in JET-ILW

To study the impact of strong impurity radiation on the energy confinement and the discharge stability at JET-ILW, dedicated high-density, highly heated experiments with neon seeding have been performed. In these experiments an increase of radiation in the core and especially the pedestal region plus a characteristic strong X-point radiator inside the confined region have been observed. The increased radiation inside the separatrix had no impact on the energy confinement. Only at the highest neon puff rates and highest heating powers () a weak H-mode without back-transitions to L-mode (M-mode) could be achieved, while at lower heating powers or seeding levels L-modes and oscillations between H-mode and L-mode could be observed. In H-mode, at highest neon seeding levels the targets are in complete detachment. A degradation of the pedestal density profile was observed when neon was seeded. A correlation between this density pedestal degradation, the appearance of an X-point radiator and the detachment of the targets was found in neon seeded H-mode pulses.In ASTRA-TGLF core transport simulations it could be shown that the increase of the electron temperature in the core with neon seeding is attributed to a partial stabilisation of ion temperature gradient modes.

Gyrokinetic simulations of nonlinear interactions between magnetic islands and microturbulence

A conservative scheme of kinetic electrons for gyrokinetic simulations in the presence of magnetic islands has been implemented and verified in the gyrokinetic toroidal code, where zonal and non-zonal components of all perturbed quantities are solved together. Using this new conservative scheme, linear simulation of kinetic ballooning mode has been successfully benchmarked with the electromagnetic hybrid model. Simulations of nonlinear interactions between magnetic islands and the ion temperature gradient (ITG) mode in a tokamak show that the islands rotate at the electron diamagnetic drift velocity. The linear ITG structure shifts from the island O-point toward the X-point due to the pressure flattening effect inside the islands, and the nonlinear ITG structure peaks along the magnetic island separatrix because of the increased pressure gradient there.

Gyrokinetic GENE simulations of DIII-D near-edge L-mode plasmas

We present gyrokinetic simulations with the GENE code addressing the near-edge region of an L-mode plasma in the DIII-D tokamak. At radial position ρ = 0.80, simulations with the ion temperature gradient (ITG) increased by 40% above the nominal value give electron and ion heat fluxes that are in simultaneous agreement with the experiment. This gradient increase is consistent with the combined statistical and systematic uncertainty σ of the charge exchange recombination spectroscopy measurements at the 1.6σ level. Multiscale simulations are carried out with a realistic mass ratio and geometry for the first time in the near-edge. These multiscale simulations suggest that the highly unstable ion temperature gradient (ITG) modes of the flux-matched ion-scale simulations suppress electron-scale transport, such that ion-scale simulations are sufficient at this location. At radial position ρ = 0.90, nonlinear simulations show a hybrid state of ITG and trapped electron modes, which was not expected from linear simulations. The nonlinear simulations reproduce the total experimental heat flux with the inclusion of E × B shear effects and an increase in the electron temperature gradient by ∼23%. This gradient increase is compatible with the combined statistical and systematic uncertainty of the Thomson scattering data at the 1.3σ level. These results are consistent with previous findings that gyrokinetic simulations are able to reproduce the experimental heat fluxes by varying input parameters close to their experimental uncertainties, pushing the validation frontier closer to the edge region.

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.

Effects of trapped electrons and impurity ions on ITG modes in reversed-field pinch plasmas

The ion temperature gradient (ITG) driven electrostatic modes in reversed-field pinch (RFP) plasma are studied with inclusion of trapped electrons and impurity ions in the gyrokinetic integral equation. The effects of magnetic shear, safety factor, electron density gradient, electron to ion temperature ratio, temperature gradient of trapped electrons are studied. As is well known, the ITG modes in RFPs are more stable than that in tokamaks, and are hard to be excited. However, it is demonstrated that the impurity ions with density gradients opposite to those of primary ions enhance the ITG instabilities. The trapped electrons also have a destabilization effect on ITG modes. In addition, trapped electron modes are found unstable in certain parameter regimes when the ion temperature gradient does not exist.

Global gyrokinetic simulation of turbulence driven by kinetic ballooning mode

Nonlinear evolution of the kinetic ballooning mode (KBM) is investigated by extending the global toroidal gyrokinetic simulation code (GKNET) to an electromagnetic regime. It is found that the saturation process of KBM, which is unstable at high normalized pressure β, is significantly different from the ion temperature gradient (ITG) mode, which is unstable at low β. The KBMs get saturated by producing zonal flows and zonal magnetic fields. The production of zonal flow is weak in the initial saturation phase of KBM, which is in contrast to the ITG mode which produces strong zonal flows in the initial saturation phase. However, strong zonal flows are produced in the subsequent evolution of KBM, and a quasisteady state of KBM turbulence is established. In addition to the zonal flows, some low toroidal number modes, which are linearly stable against the KBM, dominate the KBM turbulence. The strong zonal magnetic field is also produced by the KBM. These zonal modes regulate the KBM turbulence.

Persistence of Ion Temperature Gradient Turbulent Transport at Finite Normalized Pressure.

Plasma β dependence of electromagnetic turbulent transport is investigated by means of gyrokinetic simulations with self-consistent change of the equilibrium magnetic field. It is found that energy transport due to ion-temperature-gradient (ITG) driven turbulence does not decrease with increasing β; that is, the ion energy diffusivity does not decrease, and the electron energy diffusivity increases with β. This is because magnetic fluctuations are significantly influenced by the background magnetic field structure change with β by the Pfirsch-Schluter current. The magnetic field change weakens the suppression effect of magnetic perturbations on the growth of the ITG mode, and it also suppresses nonlinear zonal flow production. The influence of the magnetic field change is significant as the global magnetic shear increases.