Formation Of Transport Barrier Research Articles

Global flux-driven gyrofluid simulations of internal transport barrier formation by external torque in weak magnetic shear configuration

Using a global gyrofluid code, we conduct comprehensive flux-driven simulations incorporating external heat and momentum sources/sinks. They involve an external torque ramp while keeping a heating power constant. Simulations show the formation of internal transport barriers (ITBs) in ion heat and parallel momentum in weak magnetic shear configuration. The ITB formation is attributed to the reduction of turbulent transport accompanied by strong E×B flow shear generation, where the majority of E×B shear arises from plasma rotation shear. In contrast, only a minor confinement improvement is observed in strong magnetic shear configuration. There exists a substantial difference in the rotation shear between the weak and strong magnetic shear cases at the same amount of external torque. Transport analysis implies that the difference comes from intrinsic rotation generated by residual stress. Global linear simulations demonstrate that the effect of E×B shear on k∥ symmetry breaking is enhanced in weak magnetic shear configuration, explaining the origin of the difference in intrinsic rotation. The ITB is not achieved when the external torque is injected in counter-direction to intrinsic rotation. It is also shown that threshold torque for the ITB formation depends on power and the ratio of external torque to heating power should exceed a critical value to attain an ITB.

Statistical analysis of magnetic divertor configuration influence on H-mode transitions

DIII-D plasmas are compared for two upper divertor configurations: with the outer strike point on the small angle slot (SAS) divertor target and with the outer strike point on the horizontal divertor target (HT). Scanning the vertical distance between the magnetic null point and the divertor target over a range 0.10–0.16 m is shown to increase the threshold power, Pth , and edge plasma power, PLoss , for the low-to-high confinement (L–H) and H–L transitions respectively, by up to a factor of 1.4. The X-point height scans were performed at three L-mode core plasma line average electron densities, n¯e= 1.2, 2.2 and 3.6 ×1019m−3 , to investigate the density dependence of divertor magnetic configuration influence on Pth . The X-point height, Zx-pt , was further extended across the range 0.16–0.22 m with the more open HT divertor configuration, for which a clear decrease in Pth with increasing Zx-pt is observed. The dependence of Pth on divertor magnetic geometry is further investigated using a time-dependent probability density function (PDF) model and information geometry to elucidate the roles played by pedestal plasma turbulence and perpendicular velocity flows. The degree of stochasticity of the plasma turbulence is observed to be sensitive to the plasma heating rate. The calculated square of the information rate shows changes in the relative density fluctuations and perpendicular velocity PDFs begin 2–5 ms prior to the L–H transition for three plasmas; providing a crucial measurement of the dynamic timescale of external transport barrier formation. Additionally, both information length and rate provide potential predictors of the L–H transition for these plasmas.

Formation, propagation and conversion of transport barriers triggered by dynamical critical gradient in tokamak plasmas

In this work, we propose a reduced model with a dynamical critical gradient to study the formation, propagation, and conversion of the transport barriers. In contrast to the commonly adopted static critical gradient, an evolving critical gradient self-consistently softens the profile stiffness, so as to facilitate the generation of transport barriers. This is especially crucial to the internal transport barrier (ITB) formation. Numerically, we show that the inhomogeneity of turbulent and neoclassical diffusivities can induce the global wave front propagation of the transport barrier. When the heating power ramps quickly, the ITB propagates unidirectionally to the edge region and converts into an edge transport barrier. For slow power ramping, the propagating ITB will bifurcate into bidirectional wavefronts and finally convert into a steady double transport barrier state. Our model uncovers the vital role of a dynamical ‘profile-stiffness’ in depicting the global dynamics of the transport barrier.

Saturation of Fishbone Instability by Self-Generated Zonal Flows in Tokamak Plasmas.

Gyrokinetic simulations of the fishbone instability in DIII-D tokamak plasmas find that self-generated zonal flows can dominate the nonlinear saturation by preventing coherent structures from persisting or drifting in the energetic particle phase space when the mode frequency down-chirps. Results from the simulation with zonal flows agree quantitatively, for the first time, with experimental measurements of the fishbone saturation amplitude and energetic particle transport. Moreover, the fishbone-induced zonal flows are likely responsible for the formation of an internal transport barrier that was observed after fishbone bursts in this DIII-D experiment. Finally, gyrokinetic simulations of a related ITER baseline scenario show that the fishbone induces insignificant energetic particle redistribution and may enable high performance scenarios in ITER burning plasma experiments.

β-induced Alfvén eigenmodes and fishbone driven by energetic electrons in EAST electron heating dominant plasmas

β-induced Alfvén eigenmodes (e-BAE) and fishbone (e-fishbone) driven by the energetic electrons have been identified on the experimental advanced superconducting tokamak in electron cyclotron resonance heating (ECRH), low hybrid wave (LHW) together with neutral beam injection plasma and in pure radio frequency heating plasma, respectively. The e-BAE is predominantly located near the pedestal region, and its frequency typically displays a chirping characteristic. The frequency is notably associated with the count of energetic electrons with energy between 30 keV and 80 keV driven by LHW and undergoes three distinct phases. A direct correlation is established between the frequency of e-BAE and the global parameters, including electron temperature, Alfvén velocity and βp . The e-fishbone emerges in the core and becomes noticeable shortly after ECRH initiation, in an electron heated dominant and core confinement enhanced discharge. Experimental analyses suggest that the e-fishbone ties closely with both the population and energy distribution of energetic electrons. These observations can provide valuable insights into the behavior of energetic electrons, which play a crucial role in the formation and sustainability of transport barriers within hybrid plasma, both in the internal and pedestal regions.

Role of the edge electric field in the resonant mode-particle interactions and the formation of transport barriers in toroidal plasmas

The impact of an edge radial electric field on the particle orbits and the orbital spectrum in an axisymmetric toroidal magnetic equilibrium is investigated using a guiding centre canonical formalism. Poloidal and bounce/transit-averaged toroidal precession frequencies are calculated, highlighting the role of the radial electric field. The radial electric field is shown to drastically modify the resonance conditions between particles with certain kinetic characteristics and specific perturbative non-axisymmetric modes, and to enable the formation of transport barriers. The locations of the resonances and the transport barriers that determine the particle, energy and momentum transport are shown to be accurately pinpointed in the phase space by employing the calculated orbital frequencies.

Detection of Internal Transport Barrier in the T-10 Tokamak Using Thomson Scattering Diagnostics

At the T-10 tokamak, the Thomson scattering diagnostics was upgraded in 2016. The new system is based on the Nd:YAG laser with the pulse repetition frequency of 100 Hz; the measurements can be performed every 10 ms during the entire plasma discharge at the spatial resolution of up to 5 mm. Using the upgraded diagnostics, the electron temperature measurements were performed during experimental campaigns of the T-10 tokamak in 2016−2018. In the regimes with the electron cyclotron resonance heating, it was demonstrated that the regions with increased temperature gradients form in the plasma, which was interpreted as the formation of internal transport barrier. The ASTRA-code-based simulations of the current profile time evolution made it possible to correlate the positions of the barrier and the rational surface q = 1.

On the transport of tracer particles in two-dimensional plasma edge turbulence

Shear flows in turbulent fluids have been known to act as transport barriers for some time. An example of a shear flow generating mechanism is the E×B shear in plasma, which has a substantial impact on the dynamics of magnetic confinement fusion devices. The influence of this may be seen in the scrape-off layer where blobs or filaments may be sheared and velocity impacted, and in the edge and core of the plasma, where the formation of transport barriers and suppression of turbulence is strongly associated with such shearing effects. A dynamical picture of transport through these effects has been elusive—the development of a reduced model would be beneficial. We consider the application of an “observational” random walk to such transport, in order to determine whether it is a suitable approach upon which to base the development of reduced models. The observational random walk is modification of the random walk approach, introducing an intrinsic time separating observations, which reproduces the basic results of previous random walk models given a Gaussian jump function, assuming spatially homogenous jump function. We demonstrate that the jump function can be inferred from the statistics of passive particles propagated by E×B drift on a synthetic turbulence field and that the transport equation found from the jump function matches the expected diffusive transport very well. We, then, consider passive particles on simulations of the classic and modified Hasagawa–Wakatani equations in a statistical steady state for a variety of adiabaticity values and find normal transport in the near-hydrodynamic limit. When zonal flows appear, we find jump functions with non-Gaussian features, which result in transport equations with fractional differential terms in addition to, or in place of, diffusion terms. We surmise that the non-local fractional terms are related to the zonal flows acting as transport barriers. Overall, we find that the approach developed is a suitable starting point for the development of reduced models.

A simple model for internal transport barrier induced by fishbone in tokamak plasmas

Fishbone bursts have been observed to strongly correlate to internal transport barrier (ITB) formation in a number of tokamak devices. A simple model incorporating the fishbone dynamics and ion pressure gradient evolution is proposed in order to investigate the key physics parameters assisting the triggering of ITB. The time evolution of fishbone is described by the well-known predator–prey model. For each burst cycle, the energetic particles (EPs) resonantly interact with fishbone and are radially expelled from inner region leading to a radial current. A compensating bulk plasma return current and, hence, poloidal flow can be induced if the fishbone cycle frequency is greater than the poloidal flow damping rate. When the shear of the poloidal flow exceeds a critical value, the turbulent fluctuations are suppressed and the bulk ion pressure gradient transits to the high-confinement state. It is shown that this process is only sensitive to the deposition rate of the trapped EPs within the $q=1$ surface, but not sensitive to other parameters. A quantitative formula for the shearing rate of poloidal flow induced by fishbone bursts is derived and verified numerically.

Impact of avalanche type of transport on internal transport barrier formation in tokamak plasmas.

In magnetic fusion plasmas, a transport barrier is essential to improve the plasma confinement. The key physics behind the formation of a transport barrier is the suppression of the micro-scale turbulent transport. On the other hand, long-range transport events, such as avalanches, has been recognized to play significant roles for global profile formations. In this study, we observed the impact of the avalanche-type of transport on the formation of a transport barrier for the first time. The avalanches are found to inhibit the formation of the internal transport barrier (ITB) observed in JT-60U tokamak. We found that (1) ITBs do not form in the presence of avalanches but form under the disappearance of avalanches, (2) thesurface integral of avalanche-driven heat fluxe is comparable to the time rate change of stored energy retained at the ITB onset, (3) the mean E × B flow shear is accelerated via the ion temperature gradient that is not sustained under the existence of avalanches, and (4) after the ITB formation,avalanches are damped inside the ITB, while they remain outside the ITB.

A self-sustaining mechanism for internal transport barrier formation in HL-2A tokamak plasmas

The formation of Internal Transport Barrier (ITB) is studied in HL-2A plasmas by means of nonlinear gyrokinetic simulations. A new paradigm for the ITB formation is proposed in which different physics mechanisms play a different role depending on the ITB formation stage. In the early stage, fast ions, introduced by Neutral Beam Injection ion system, are found to stabilize the thermal-ion-driven instability by dilution, thus reducing the ion heat fluxes and finally triggering the ITB. Such dilution effects, however, play a minor role after the ITB is triggered as electromagnetic (EM) effects are dominant in the presence of established high pressure gradients. We define the concept of ITB self-sustainment, as the low turbulence levels found within the fully formed ITB are consequences of large scale zonal flows, which in turn are fed by a non-linear interplay with large scale high frequency EM perturbations destabilized by the ITB itself.

Turbulence stabilization in tokamak plasmas with high population of fast ions

This letter provides a new physical insight into the fast ion effects on turbulence in plasmas having a high fast ion fraction and peaked fast ion density profile. We elucidate turbulence stabilization mechanisms by fast ions that result in internal transport barrier formation in the fast ion regulated enhancement mode plasma. Both linear and nonlinear gyrokinetic simulations show that the dominant turbulence suppression mechanisms are the dilution effects. In particular, we find that turbulence can be sufficiently suppressed solely by an inverted main ion density gradient due to fast ions, for the first time. New physical findings reported here improve our understanding of fast ion effects on turbulence, essential for fusion energy production where . Moreover, they will open up a new methodology to control plasma turbulence applicable to a wide range of plasma confinement regimes.

Oscillation dynamics of m/n= 3/1 double tearing mode

This study systematically investigates the roles of the plasma viscosity and resistivity in the oscillation dynamics during the decay phase of the m/n = 3/1 double tearing mode using the Ci-Liu-Ti (CLT) code. The primary objectives of this research are to examine the driving and suppressing mechanisms of the oscillation. The oscillation and steady-state are the result of the competition between the external injection and the reconnection annihilation of magnetic flux during the decay phase. In a regime with a higher viscosity (or a lower resistivity), the steady-state arises from the significant damping (weak generation) of plasma flows, resulting in the formation of saturated islands. In a regime with a lower viscosity (or a higher resistivity), the suppression of the oscillation amplitude can be attributed to a strong residual flow that quickly takes the injected magnetic flux away toward to the reconnection region, which caused no enough accumulated magnetic flux to drive oscillations and the system evolves toward a steady-state configuration. The steady-state condition results in the generation of a narrow radial vortex region which promotes formation of internal transport barriers. The upper threshold of the resistivity within the low-resistivity regime to achieve a steady-state decreases as the viscosity increases.

Self-generated vortex flows in a tokamak magnetic island with a background flow

We present a gyrokinetic theory of self-generated E × B vortex flows in a magnetic island in a collisionless tokamak plasma with a background vortex flow. We find that the long-term evolution of the self-generated vortex flows can be classified into two regimes by the background vortex flow potential Φ, with an asymptotic criterion given by , where T is temperature, ε is the inverse aspect ratio and r is the radial coordinate. We find that the background vortex flow above the criterion significantly weakens the toroidal precession-induced long-term damping and structure change of the self-generated vortex flows. That is, the finite background vortex flow is beneficial to maintain the self-generated vortex flows, favorable to an internal transport barrier formation. Our result indicates that the island boundary region is a prominent location for triggering the transition to an enhanced confinement state of the magnetic island.

Study of turbulent transport in magnetized plasmas with flow using symplectic maps

Turbulent transport in a magnetized plasma of the kind found in tokamaks is modeled by a 2D wave spectrum that allows reduction to a symplectic map. The properties of particle transport when chaos sets in are analyzed in various circumstances including finite Larmor radius (FLR) effects and a background plasma flow. For large wave amplitudes, regular particle orbits become chaotic, which represents a type of Lagrangian turbulence. When chaos becomes global, it leads to the loss of particle confinement. Poloidal flows tend to decrease the chaos in some regions, and they can give rise to the formation of transport barriers. FLR effects not only reduce chaos but also give rise to non-local behavior. Thus, when the particles have a thermal distribution of Larmor radii, a non-Gaussian particle distribution function in space is obtained. However, the transport preserves its diffusive scaling when there is no flow. Previous results about the dependence of the diffusion coefficient with amplitude are re-derived analytically and numerically taking into account FLR effects. In the presence of general poloidal flows, the transport has to be described by a two-step map. They modify the nature of transport in the direction of the flow from diffusive to ballistic to super-ballistic, depending on the type of flow. The transverse transport, in turn, shows suppression of the oscillations with wave amplitude that are present in the absence of flow. When the plasma flow varies linearly with radius, the transport can be studied with a similar single-step map, and the transverse diffusion coefficient is reduced while parallel transport can become super-ballistic. For non-monotonic flows, there are accelerating modes that can produce ballistic-like particles while the bulk of the particles behaves diffusively.

Global gyrofluid simulations of turbulence in tokamak plasmas

We develop a tokamak plasma turbulence simulation code based on a fluid model using the BOUT++ framework. The model corresponds to a global extension of the local gyrofluid model (Beer and Hammett, 1996 [11]), where external sources and neoclassical equilibrium flows are included but some higher-order kinetic effect terms are ignored for the global nonlinear simulation. Hyper-viscous damping terms are introduced to compensate for the stabilizing effect of the missing kinetic terms. With an appropriate hyper-viscosity, we obtain reasonable linear and nonlinear results on ion temperature gradient modes compared to global gyrokinetic simulations. We perform global flux-driven simulations with self-consistently evolving profiles with external heat sources/sinks. A spectral analysis of heat flux in flux-driven turbulence exhibits a broad 1/f (f is the frequency) spectrum observed in gyrokinetic simulations, indicating the prevalence of self-organized criticality-like heat avalanches. We also investigate the effect of an externally imposed vorticity source on turbulence suppression and transport barrier formation. It is shown that internal transport barrier can form even for a monotonic safety factor profile when a sufficiently strong vorticity source is applied.

The internal transport barrier formation on EAST tokamak during the fishbone instability

The internal transport barrier (ITB) which is related to the fishbone instability has been observed on the Experimental Advanced Superconducting Tokamak (EAST) in ELMy H-mode discharges. An interpretation of the formation of the ITB on EAST tokamak is provided, based on both analytical and numerical calculations. The fishbone instability induces the redistribution of fast ions and leads to the accumulation of fast ions in a local region where the ITB is going to appear. Correspondingly, the gradients of fast ions are enhanced, where the ion temperature gradient (ITG) mode exists. Fast ions can interact with the ITG mode through the dilution, Shafranov shift and wave-particle resonance mechanisms. It is found that the ITG mode is stabilized by fast ions and the stabilizing effects are determined mainly by the density, temperature and their gradients of fast ions. The enhanced density and temperature gradients of fast ions lead to a stronger stabilizing effect on ITG mode. Compared with the stabilizing effect before the appearance of fishbone instability, the stabilization on ITG mode is enhanced after the fishbone instability, which is beneficial to the formation of the ITB.

Experimental study of the core instability before and after internal transport barrier formation in EAST

In a reversed shear discharge in the Experimental Advanced Superconducting Tokamak (EAST), an internal transport barrier (ITB) was formed. For the first time, the turbulence in the plasma core before and after the formation of the ITB in EAST was investigated by microwave reflectometry. It was found that during the formation of the ITB, the turbulence inside the barrier was not completely suppressed. The density fluctuation level decreased in the early phase of ITB and significantly increased later. It was found that the density fluctuation level increased with the density gradient after the appearance of reversed shear Alfvén eigenmodes (RSAEs). The change in turbulence, measurable by the reflectometer, did not affect the rate of increase in core density. Through the measurement of RSAEs, it was found that the formation of this ITB may be related to the minimum q (q min). q min was close to 2.

Characteristics of electron temperature profile stiffness in electron-heated plasmas on EAST

A very high core electron temperature (T e0 ∼ 10 keV) plasma has been established and stably sustained by applying both lower hybrid wave (LHW) and on-axis electron cyclotron resonance heating (ECRH) in the Experimental Advanced Superconducting Tokamak (EAST). In this work, power balance analysis shows that the increase of ECRH power can increase the normalized T e gradient significantly at the plasma core region (ρ < 0.6), but does not change the T e profile stiffness in the low-density L-mode plasmas. This has been considered to be due to a strong synergistic effect between ECRH and LHW. Furthermore, three distinguishable stages characterized by different T e profile stiffnesses can be identified from the density ramp-up in the electron-heated plasma on EAST. A stronger T e profile stiffness at ρ = 0.3 has been observed in the Stage-II, where the LHW power deposition gradually moves away from the plasma core region, following the electron density increases. Furthermore, the formation of an internal plasma density transport barrier inside ρ ∼ 0.6, accompanied by a sudden drop in core T e and a rise in both core plasma density and ion temperature, has been observed for the first time during the transition from the Stage-II to the Stage-III when the central line-averaged plasma density reaches a threshold of 2.2 × 1019 m−3. This finding strongly affects further development of high-performance gas-fueled electron-heated plasma scenarios in EAST and suggests an advanced operational regime with a wide internal plasma density transport barrier.

Electron thermal internal transport barriers triggered by the effect of ion shielding

The formation of electron thermal internal transport barriers (eITBs) is investigated by critical gradient threshold analysis for recent EAST experiments with dominant RF heating and low torque injection. The ratio of electron temperature to ion temperature, τ = T e/T i, is identified to be the key parameter that triggers eITB in the hot electron mode through electron turbulence suppression by the effect of ion shielding. The critical gradient of electron-temperature-gradient turbulence can be greatly increased with τ due to strong electron heating and weak electron–ion energy coupling. The formation and evolution of eITB is found to be determined by the critical temperature gradient of electron turbulence.