ELMy H-mode Plasmas Research Articles

Frequency linearization of voltage controlled oscillator (VCO) for 2 µs frequency modulated continuous wave (FMCW) reflectometer.

Frequency modulated continuous wave reflectometers have been widely used to measure plasma density profiles in many magnetic fusion devices. The frequency modulation (FM) time of the KSTAR reflectometer was 20 µs, that is, the FM rate was 50kHz. However, the edge density of the KSTAR tokamak fluctuates typically over the frequency range of 20-50kHz in the ELMy H-mode plasmas. Therefore, the density profile changes significantly during the FM time, causing significant distortion in the density profile measurements. The FM rate must be increased at least ten times faster than the density fluctuation frequency. A new voltage controlled oscillator (VCO) driver is designed based on the AD8067 operational amplifier to achieve a FM time of 2 µs. The VCO output frequency is linearly modulated by predistorting the tuning voltage of VCO using a 400 MSamples/s arbitrary waveform generator. The resulting output frequency of the VCO is measured by mixing the VCO output with a fixed frequency synthesizer signal and measuring the mixer intermediate frequency using a 2.5 GSamples/s digitizer. The tuning voltage is adjusted to minimize the amount by which the measured frequency deviates from the linearly modulated frequency. A simple linear adjustment proportional to the deviation does not effectively suppress the nonlinearity in FM. The response time of the VCO driver and the VCO input interface should be taken into account by solving the entire circuit equation. In this paper, the developed VCO driver circuit and the frequency linearization method are described.

Observation of internal transport barrier evolution in ELMy H-mode plasma in the EAST

Observation of internal transport barrier evolution in ELMy H-mode plasma in the EAST

Overview of recent experimental results on the EAST Tokamak

Since the last IAEA-FEC in 2021, significant progress on the development of long pulse steady state scenario and its related key physics and technologies have been achieved, including the reproducible 403 s long-pulse steady-state H-mode plasma with pure radio frequency (RF) power heating. A thousand-second time scale (∼1056 s) fully non-inductive plasma with high injected energy up to 1.73 GJ has also been achieved. The EAST operational regime of high β P has been significantly extended (H 98y2 > 1.3, β P ∼ 4.0, β N ∼ 2.4 and n e/n GW ∼ 1.0) using RF and neutral beam injection (NBI). The full edge localized mode suppression using the n = 4 resonant magnetic perturbations has been achieved in ITER-like standard type-I ELMy H-mode plasmas with q 95 ≈ 3.1 on EAST, extrapolating favorably to the ITER baseline scenario. The sustained large ELM control and stable partial detachment have been achieved with Ne seeding. The underlying physics of plasma-beta effect for error field penetration, where toroidal effect dominates, is disclosed by comparing the results in cylindrical theory and MARS-Q simulation in EAST. Breakdown and plasma initiation at low toroidal electric fields (<0.3 V m−1) with EC pre-ionization is developed. A beneficial role on the lower hybrid wave injection to control the tungsten concentration in the NBI discharge is observed for the first time in EAST suggesting a potential way toward steady-state H-mode NBI operation.

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy (E in) and angles (A in) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in, and W impurity transport.

Validated edge and core predictions of tungsten erosion and transport in JET ELMy H-mode plasmas

Predictive edge and core simulations of tungsten (W) erosion and transport in JET ITER-like wall plasmas are shown to be consistent with the experimentally inferred W density in the main plasma, within the uncertainty inherited from the measurements of the deuterium plasma conditions and from the W density measurements. The ERO2.0 code is applied to predicting the W erosion and edge transport, whereas JINTRAC predicts W transport from the pedestal top to the core plasma. The studied plasma scenarios range from L-mode to the highest-performance deuterium ELMy H-mode in JET.

Divertor power load investigations with deuterium and tritium in type-I ELMy H-mode plasmas in JET with the ITER-like wall

Divertor power load is a major challenge towards a burning plasma in a next-step tokamak. Here, the first results of a divertor power load characterisation in tritium plasmas in type-I ELMy H-mode, obtained in the JET deuterium-tritium campaign (DTE2) performed in 2021, are presented. It is demonstrated that both, transient loads due to type-I ELMs as well as the power fall-off length, do not exhibit an explicit ion mass dependence, with remarkably similar values in the tritium plasmas and in the deuterium references. This gives an improved credence to published scaling law predictions, solely based on deuterium plasma experiments. Moreover, the type-I ELM impact on the inner divertor target is studied in deuterium discharges. A slightly increased parallel energy fluence on the inner target with a factor of 1.08 compared to the outer target is observed. This is explained by the smaller major radius of the inner target.

Isotope physics of heat and particle transport with tritium in JET-ILW type-I ELMy H-mode plasmas

As part the DTE2 campaign in the JET tokamak, we conducted a parameter scan in T and D-T complementing existing pulses in H and D. For the different main ion masses, type-I ELMy H-modes at fixed plasma current and magnetic field can have the pedestal pressure varying by a factor of 4 and the total pressure changing from βN=1.0 to 3.0. We investigated the pedestal and core isotope mass dependencies using this extensive data set. The pedestal shows a strong mass dependence on the density, which influences the core due to the strong coupling between both plasma regions. To better understand the causes for the observed isotope mass dependence in the pedestal, we analysed the interplay between heat and particle transport and the edge localised mode (ELM) stability. For this purpose, we developed a dynamic ELM cycle model with basic transport assumptions and a realistic neutral penetration. The temporal evolution and resulting ELM frequency introduce an additional experimental constraint that conventional quasi-stationary transport analysis cannot provide. Our model shows that a mass dependence in the ELM stability or in the transport alone cannot explain the observations. One requires a mass dependence in the ELM stability as well as one in the particle sources. The core confinement time increases with pedestal pressure for all isotope masses due to profile stiffness and electromagnetic turbulence stabilisation. Interestingly, T and D-T plasmas show an improved core confinement time compared to H and D plasmas even for matched pedestal pressures. For T, this improvement is largely due to the unique pedestal composition of higher densities and lower temperatures than H and D. With a reduced gyroBohm factor at lower temperatures, more turbulent drive in the form of steeper gradients is required to transport the same amount of heat. This picture is supported by quasilinear flux-driven modelling using TGLF-SAT2 within Astra. With the experimental boundary condition TGLF-SAT2 predicts the core profiles well for gyroBohm heat fluxes >15 , however, overestimates the heat and particle transport closer to the turbulent threshold.

Parameter dependence of density pedestal width and correlation with pedestal turbulence in EAST type-I ELMy H-mode plasma

The density pedestal structures in a series of type-I ELMy H-mode discharges (0.6 MA/2.5 T) with upper single null (USN) divertor configuration have been studied in EAST. To obtain the pedestal width, density profiles are fitted via the modified hyperbolic tangent (mtanh) function and the two-line method. The pedestal widths from these two methods have a similar tendency. Significant impacts of pedestal poloidal beta and pedestal collisionality on the density pedestal width have been observed for the first time in EAST. The non-linear regressions suggest the density width scaling in type-I ELMy EAST plasma are wne∝c(νped⁎)0.27−0.33⋅(βp,ped)0.5−0.58. The density fluctuations in the bottom of the pedestal from the fluctuation reflectometry imply an enhanced density fluctuation (Pfluc) and an obvious spectral broadening in the high collisionality regime of νped⁎≥3. The experimental results suggest that the pedestal turbulence would play an important role in the evolution of the pedestal structure.

Comparison of MHD stability properties between QH-mode and ELMy H-mode plasmas by considering plasma rotation and ion diamagnetic drift effects

Magnetohydrodynamic (MHD) stability at tokamak edge pedestal in a quiescent H-mode (QH-mode) and type-I ELMy H-mode plasmas in DIII-D experiment was analyzed by considering plasma rotation and ion diamagnetic drift effects. QH-mode plasma is marginally stable to kink/peeling mode (K/PM), but ELMy H-mode one is almost unstable to peeling-ballooning mode (PBM). It was identified that there are three physics features responsible for the difference in the MHD stability properties between QH-mode plasma and ELMy H-mode one. These are the distance of pedestal foot from the last closed flux surface (LCFS), the amount of the ion diamagnetic drift frequency at pedestal, and impact of coupled rotation and ion diamagnetic drift effects. These features were confirmed through the numerical experiments that the stability properties of the QH-mode plasma can be changed to that of the ELMy H-mode one by shifting the plasma profiles inward in the radial direction and halving the ion diamagnetic drift frequency. The reasons of the change in the stability properties are thought as that K/PM is stabilized due to the inward shift of the bootstrap current profile, and PBM is destabilized due to the reduction of the coupled rotation and ion diamagnetic drift stabilizing effect. Importance of these features was validated through numerical experiments with experimental data of other QH-mode plasmas in DIII-D. All the results show that MHD stability properties of QH-mode plasma can be obtained in case that pedestal foot is close to LCFS, ion diamagnetic drift frequency is large due to high ion temperature, and strong rotation shear exists near pedestal.

Non-parametric inference of impurity transport coefficients in the ASDEX Upgrade tokamak

We present a non-parametric inference of impurity transport coefficients by using charge exchange recombination spectroscopy measurements of Ne X, Ne VIII, O VIII, and C VI lines. Due to their close atomic numbers, neon, oxygen and carbon impurity ions are assumed to have the same diffusion coefficient D and convection velocity v. Unlike conventional techniques that modulate or perturb the impurity contents, we employ a quasi-stationary plasma with static impurity profiles. Since the ratio of v to D only describes the equilibrated profile of the sum of all impurity charge states, steady-state measurements can still decouple D and v if different charge states are simultaneously observed. We have formulated a non-parametric analysis framework based on the Bayesian probability theory and conducted transport coefficient measurements for a Type III ELMy H-mode plasma at ASDEX Upgrade. The charge exchange reactions with the background neutrals, which are known to affect the impurity charge state balance, are taken into account by introducing additional free parameters. While D at the pedestal is close to the neoclassical level (1 m s−2), a large diffusion coefficient and a strong outward convection are inferred right inside the pedestal top.

Experimental investigation on divertor tungsten sputtering with neon seeding in ELMy H-mode plasma in EAST tokamak

Neon (Ne) seeding is used to cool the edge plasma by radiation to protect the divertor tungsten (W) target in the Experimental Advanced Superconducting Tokamak (EAST). The W sputtering in the outer divertor target with Ne seeding is assessed by the divertor visible spectroscopy system. It is observed that the W sputtering flux initially increases with Ne concentration in the divertor despite the decreasing plasma temperature. After reaching a maximum around 25 eV, the W sputtering rate starts to decrease, presenting a suppression effect. The effect on the divertor W sputtering flux and yield due to the competition between the increase of the Ne concentration and the decrease of the plasma temperature is discussed. The results show that enough Ne seeding is essential to effectively reduce the electron temperature and thus to suppress W sputtering. Moreover, ELM suppression is observed when Ne and W impurities enter the core plasma, which could be correlated to the enhanced turbulence transport in the pedestal.

Study of the mechanism of ITB formation and sustainment with optimized q profiles in ELMy H mode discharges on the EAST

Recently, improved confinement plasma with the internal transport barrier (ITB) on the Experimental Advanced Superconducting Tokamak has been achieved through extensive experiments in high β N scenario development. ELMy H-mode plasmas are heated by lower hybrid wave and neutral beam injection, with B T = 1.6 T, I p = 400/450 kA and q 95 in the range 3.7–4.6. ITB with optimized reversal q profiles and central flat q profiles are observed in these H-mode plasmas, which lead to better confinement. The electron temperature ITB, ion temperature ITB, and electron density ITB exhibit different characteristics. In this paper, the formation and sustainment of ITB on n e channel has great relationship with magnetohydrodynamics (MHD) behaviors. For T e channel, q profile has a dominating effect. Reversal magnetic shear and MHD behaviors are both effective in T i ITB formation. Suppression of turbulence in the ITB region is the common reason on all the channels. In addition, the power threshold of the ITB formation should be concentrated on electrons and ions separately. The results of this study can widen the understanding of ITB and be helpful in exploring better candidate operational scenarios for future fusion devices.

Dependence of upstream SOL density shoulder on divertor neutral pressure observed in L-mode and H-mode plasmas in the EAST superconducting tokamak

Upstream density profiles in the scrape-off layer (SOL) have been examined in low-confinement mode (L-mode) and high-confinement mode (H-mode) plasmas in the EAST superconducting tokamak. A weak density shoulder forms in the near SOL region in upper single-null configurations when the neutral pressure measured at the lower divertor exceeds a threshold value of 2 × 10−2 Pa in L-mode plasmas. When the neutral pressure is below this threshold, the weak density shoulder is absent and the sidebands of the lower hybrid waves associated with SOL parametric instabilities are reduced. Active detachment control with neon–deuterium seeding demonstrate that the weak density shoulder can form before the onset of the outer divertor detachment as long as the neutral pressure is above the threshold. Furthermore, no remarkable expansion of a shoulder is observed during divertor detachment, suggesting that divertor detachment is not a necessary condition for the formation or growth of a density shoulder. Through the increase in neutral pressure in the lower divertor by an order of magnitude, the weak shoulder was observed to expand into the far SOL and reach the leading edge of the limiter. The results in L-mode discharges identified the neutral pressure in the lower divertor as a primary factor for the formation of an SOL density shoulder in the upper single-null discharges. For the type-I ELMy H-mode plasmas, a similar density shoulder was detected during the inter-ELM phase when the neutral pressure in the lower divertor exceeded a threshold value of 4 × 10−2 Pa. On the other hand, the shoulder was absent when the divertor neutral pressure went below this threshold even though the plasma discharge was conducted with a higher core line-averaged density and divertor collisionality. This is consistent with the observations in L-mode plasmas. The neutral particle ionization of the working gas is thus believed to play a key role during the formation of the SOL density shoulder in the EAST tokamak.

Comparing spontaneous and pellet-triggered ELMs via non-linear extended MHD simulations

Injecting frozen deuterium pellets into an ELMy H-mode plasma is a well established scheme for triggering edge localized modes (ELMs) before they naturally occur. This paper presents non-linear simulations of spontaneous type-I ELMs and pellet-triggered ELMs in ASDEX Upgrade performed with the extended MHD code JOREK. A thorough comparison of the non-linear dynamics of these events is provided. In particular, pellet-triggered ELMs are simulated by injecting deuterium pellets into different time points during the pedestal build-up described in A Cathey et al (2020 Nuclear Fusion 60 124007). Realistic ExB and diamagnetic background plasma flows as well as the time dependent bootstrap current evolution are included during the build-up to accurately capture the balance between stabilising and destabilising terms for the edge instabilities. Dependencies on the pellet size and injection times are studied. The spatio-temporal structures of the modes and the resulting divertor heat fluxes are compared in detail between spontaneous and triggered ELMs. We observe that the premature excitation of ELMs by means of pellet injection is caused by a helical perturbation described by a toroidal mode number of n = 1. In accordance with experimental observations, the pellet-triggered ELMs show reduced thermal energy losses and a narrower divertor wetted area with respect to spontaneous ELMs. The peak divertor energy fluence is seen to decrease when ELMs are triggered by pellets injected earlier during the pedestal build-up.

Load-tolerant external matching unit based on a wideband hybrid coupler for the ICRH system of DTT

In the first phase of the Divertor Tokamak Test Facility (DTT), the goal of coupling a RF power of 3 MW to the plasma will be reached with a single ICRH module constituted by 4 generators, 2 or 4 External Matching Units (EMU) and 2 antennas. The operation with high levels of RF power into ELMy H-mode plasmas requires a robust and reliable Load-Tolerant EMU. Currently the most appropriate solution, ensuring a VSWR at the generators under the safety threshold of 1.5:1, will be based on a Wideband 2-section coaxial Hybrid Coupler operating in the frequency range of 60−90 MHz. HFSS and CST-MS simulation software was used for the optimization of the characteristic impedances and lengths of the transversal and longitudinal branches aimed to achieve excellent performance (amplitude and phase balance at the output ports, and low reflection at the input port) over the entire frequency range. The analysis has been carried out for Zo = 50 Ω and 30 Ω coaxial transmission lines in terms of S-parameters and electric field. Several good solutions with an amplitude imbalance about ±0.25 dB and a return loss lower than −20 dB in a large part of the frequency range will be presented. Then a simulation of the overall (end-to-end) EMU (Wideband Hybrid Coupler with shifter, stub and service stub in each of the 2 output branches) exploiting a suitable circuital solver of Ansys Electronics tool of HFSS, has been implemented. The resilience of the EMU has been tested determining the VSWR in input and the amplitude balance at the outputs for a typical resistive and reactive loading variation occurring during ELMs at f = 75 MHz and compared with the results at 65 MHz and 85 MHz.

Isotope dependence of the type I ELMy H-mode pedestal in JET-ILW hydrogen and deuterium plasmas

The pedestal structure, edge transport and linear MHD stability have been analyzed in a series of JET with the ITER-like wall hydrogen (H) and deuterium (D) type I ELMy H-mode plasmas. The pedestal pressure is typically higher in D than in H at the same input power and gas rate, with the difference mainly due to lower density in H than in D (Maggi et al (JET Contributors) 2018 Plasma Phys. Control. Fusion 60 014045). A power balance analysis of the pedestal has shown that higher inter-ELM separatrix loss power is required in H than in D to maintain a similar pedestal top pressure. This is qualitatively consistent with a set of interpretative EDGE2D-EIRENE simulations for H and D plasmas, showing that higher edge particle and heat transport coefficients are needed in H than in D to match the experimental profiles. It has also been concluded that the difference in neutral penetration between H and D leads only to minor changes in the upstream density profiles and with trends opposite to experimental observations. This implies that neutral penetration has a minor role in setting the difference between H and D pedestals, but higher ELM and/or inter-ELM transport are likely to be the main players. The interpretative EDGE2D-EIRENE simulations, with simultaneous upstream and outer divertor target profile constraints, have indicated higher separatrix electron temperature in H than in D for a pair of discharges at low fueling gas rate and similar stored energy (which required higher input power in H than in D at the same gas rate). The isotope dependence of linear MHD pedestal stability has been found to be small, but if a higher separatrix temperature is considered in H than in D, this could lead to destabilization of peeling-ballooning modes and shrinking of the stability boundary, qualitatively consistent with the reduced pedestal confinement in H.

Pedestal magnetic turbulence measurements in ELMy H-mode DIII-D plasmas by Faraday-effect polarimetry

Internal magnetic fluctuation measurements are utilized to identify turbulence associated with micro-tearing modes (MTM) in the DIII-D Edge-Localized-Mode (ELM)-y H-mode pedestal. Using a Faraday-effect polarimeter, magnetic turbulence (150–500 kHz) is directly observed with a typical line-averaged fluctuation amplitude of ∼0.8 G at peak frequency (250 kHz) and ∼15 G integrated over the spectrum from 150 to 500 kHz. Frequency, poloidal wavenumber, and propagation direction of the magnetic turbulence all serve to identify as MTM. Magnetic turbulence amplitude non-monotonically correlates with collision frequency, peaks off mid-plane, and correlates with electron temperature gradient evolution between ELMs, consistent with MTM features identified from theory and gyro-kinetic simulation. The magnetic turbulence growth correlates with confinement degradation in ELMy H-mode plasmas during a slow density ramp. These internal measurements provide unique constraints toward developing physics understanding and validating models of the H-mode pedestal for future devices.

Validation of EDGE2D-EIRENE and DIVIMP for W SOL transport in JET

Abstract Tungsten sputtering rates and density profiles predicted using the edge plasma codes EDGE2D-EIRENE and DIVIMP are found to agree within a factor of 4 with measurements of neutral and singly-ionized W spectral line emission in the JET low-field side (LFS) divertor, and within a factor of 2 with SXR, VUV, and bolometric calculations of the W density in the main plasma. The edge plasma W predictions are extended to the core plasma using JINTRAC integrated core-edge modelling. Prompt redeposition of W is identified as the primary reason for the discrepancy between predicted and measured W emission in the divertor. The studied plasmas include attached divertor conditions in L-mode and type-I ELMy H-mode plasmas typical for JET. To more accurately reproduce the spectroscopically inferred W sputtering rates in EDGE2D-EIRENE, imposing the experimentally observed Be concentration of order 0.5% in the divertor is necessary. However, the W density in the main plasma is predicted to be insensitive to whether or not W is sputtered by Be at the divertor targets. Instead, the majority of the predicted core W originated in L-mode from sputtering due to fast D charge-exchange atoms at the W-coated tiles above the LFS divertor, and in H-mode due to D and W ions at the targets during ELMs.

Isotope dependence of energy, momentum and particle confinement in tokamaks

The isotope dependence of plasma transport will have a significant impact on the performance of future D-T experiments in JET and ITER and eventually on the fusion gain and economics of future reactors. In preparation for future D-T operation on JET, dedicated experiments and comprehensive transport analyses were performed in H, D and H-D mixed plasmas. The analysis of the data has demonstrated an unexpectedly strong and favourable dependence of the global confinement of energy, momentum and particles in ELMy H-mode plasmas on the atomic mass of the main ion species, the energy confinement time scaling as ${\tau _E}\sim {A^{0.5}}$ (Maggi et al., Plasma Phys. Control. Fusion, vol. 60, 2018, 014045; JET Team, Nucl. Fusion, vol. 39, 1999, pp. 1227–1244), i.e. opposite to the expectations based only on local gyro-Bohm (GB) scaling, ${\tau _E}\sim {A^{ - 0.5}}$ , and stronger than in the commonly used H-mode scaling for the energy confinement (Saibene et al., Nucl. Fusion, vol. 39, 1999, 1133; ITER Physics Basis, Nucl. Fusion, vol. 39, 1999, 2175). The scaling of momentum transport and particle confinement with isotope mass is very similar to that of energy transport. Nonlinear local GENE gyrokinetic analysis shows that the observed anti-GB heat flux is accounted for if collisions, E × B shear and plasma dilution with low-Z impurities (9Be) are included in the analysis (E and B are, respectively the electric and magnetic fields). For L-mode plasmas a weaker positive isotope scaling ${\tau _E}\sim {A^{0.14}}$ has been found in JET (Maggi et al., Plasma Phys. Control. Fusion, vol. 60, 2018, 014045), similar to ITER97-L scaling (Kaye et al., Nucl. Fusion, vol. 37, 1997, 1303). Flux-driven quasi-linear gyrofluid calculations using JETTO-TGLF in L-mode show that local GB scaling is not followed when stiff transport (as is generally the case for ion temperature gradient modes) is combined with an imposed boundary condition taken from the experiment, in this case predicting no isotope dependence. A dimensionless identity plasma pair in hydrogen and deuterium L-mode plasmas has demonstrated scale invariance, confirming that core transport physics is governed, as expected, by the 4 dimensionless parameters ρ*, ν*, β, q (normalised ion Larmor radius, collisionality, plasma pressure and safety factor) consistently with global quasi-linear gyrokinetic TGLF calculations (Maggi et al., Nucl. Fusion, vol. 59, 2019, 076028). We compare findings in JET with those in different devices and discuss the possible reasons for the different isotope scalings reported from different devices. The diversity of observations suggests that the differences may result not only from differences affecting the core, e.g. heating schemes, but are to a large part due to differences in device-specific edge and wall conditions, pointing to the importance of better understanding and controlling pedestal and edge processes.

Investigation of the dithering-like cycles induced by deuterium pellet injection during ELMy H-mode plasmas on EAST

By injecting cryogenic deuterium pellets into H-mode plasmas with edge localized modes (ELMs) from the low-field side (LFS) on EAST tokamak, the impact of these pellets on the edge plasma can be investigated. In some cases of pellet injection (PI), the intermediate phase (I-phase), i.e, the dithering-like cycles, have been observed. Experiments show that during one dithering-like cycle, three phases can be distinguished, which are closely related to the variations in the pedestal density profiles and density gradient after PI. Phase 1 is characterized by rapid particle loss and a density pedestal crash, phase 2 has a relatively low density gradient and enhanced edge density fluctuation and phase 3 corresponds to the relatively fast recovery of the density pedestal and reduced edge density fluctuation. Through an experimental analysis, it is shown that the plasma will enter the I-phase when the edge line-averaged density increases beyond a certain threshold value after PI. H-I-H transitions are more likely to be observed when the edge line-averaged density is greater than the core line-averaged density.