Processing math: 100%
Advanced Search+
Qingrui ZHOU, Yanjie ZHANG, Chaofeng SANG, Jiaxian LI, Guoyao ZHENG, Yilin WANG, Yihan WU, Dezhen WANG. Simulation of tungsten impurity transport by DIVIMP under different divertor magnetic configurations on HL-3[J]. Plasma Science and Technology, 2024, 26(10): 104003. DOI: 10.1088/2058-6272/ad6817
Citation: Qingrui ZHOU, Yanjie ZHANG, Chaofeng SANG, Jiaxian LI, Guoyao ZHENG, Yilin WANG, Yihan WU, Dezhen WANG. Simulation of tungsten impurity transport by DIVIMP under different divertor magnetic configurations on HL-3[J]. Plasma Science and Technology, 2024, 26(10): 104003. DOI: 10.1088/2058-6272/ad6817

Simulation of tungsten impurity transport by DIVIMP under different divertor magnetic configurations on HL-3

More Information
  • Author Bio:

    Chaofeng SANG: sang@dlut.edu.cn

  • Corresponding author:

    Chaofeng SANG, sang@dlut.edu.cn

  • Received Date: February 21, 2024
  • Revised Date: July 24, 2024
  • Accepted Date: July 25, 2024
  • Available Online: July 26, 2024
  • Published Date: September 08, 2024
  • Tungsten (W) accumulation in the core, depending on W generation and transport in the edge region, is a severe issue in fusion reactors. Compared to standard divertors (SDs), snowflake divertors (SFDs) can effectively suppress the heat flux, while the impact of magnetic configurations on W core accumulation remains unclear. In this study, the kinetic code DIVIMP combined with the SOLPS-ITER code is applied to investigate the effects of divertor magnetic configurations (SD versus SFD) on W accumulation during neon injection in HL-3. It is found that the W concentration in the core of the SFD is significantly higher than that of the SD with similar total W erosion flux. The reasons for this are: (1) W impurities in the core of the SFD mainly originate from the inner divertor, which has a short leg, and the source is close to the divertor entrance and upstream separatrix. Furthermore, the W ionization source (SW0) is much stronger, especially near the divertor entrance. (2) The region overlap of SW0 and FW,TOT pointing upstream promote W accumulation in the core. Moreover, the influence of W source locations at the inner target on W transport in the SFD is investigated. Tungsten impurity in the core is mainly contributed by target erosion in the common flux region (CFR) away from the strike point. This is attributed to the fact that the W source at this location enhances the ionization source above the W ion stagnation point, which sequentially increases W penetration. Therefore, the suppression of far SOL inner target erosion can effectively prevent W impurities from accumulating in the core.

  • Power heat load on the divertor target is one of the most urgent issues for future fusion devices. The detachment scenario is believed to be a way to achieve stable H-mode operation [1, 2]. The physics of detachment is very complex. Plasma confinement varies with different gas species, materials of the first wall and divertor configurations in detached H-mode plasma [38]. Experimental results show that after detachment, plasma performance with relatively closed divertor configurations is better than that with relatively open divertor configurations in DIII-D [8, 9]. This means that exploring closed divertor physics is of great significance. The HL-2A tokamak has a very closed divertor configuration [10], which has advantages for the study of detachment physics.

    Recently, it has been reported that L-mode plasma with high power shows good confinement properties and no ELM activities via Ar seeding in ASDEX-Upgrade experiments [11], indicating the importance of L-mode detachment. It was proposed that the study of edge turbulence becomes important due to its effects of confinement degradation. In addition, the plasma profiles (such as radial profiles of electron density and temperature) redistribute after detachment [38], which has strong effects on the evolution of plasma confinement and edge turbulence. In turn, edge turbulence affected by detachment plays a role in the redistribution of plasma profiles. Hence, the study of edge turbulent transport evolution in the process of detachment is crucial. The radial profiles of relative density fluctuation at the plasma boundary are measured using reflectometry in different states of L-mode detachment in ASDEX-Upgrade [12], which shows an increasing tendency to change from the attached to the pronounced detached state, but turbulent transport is not reported. More studies of the evolution of edge turbulence and its induced transport under detachment are required, especially in plasma with a closed divertor.

    In this work, edge turbulent transport in L-mode plasma during the process of detachment is investigated in detail in the HL-2A tokamak. Based on the characteristics of detachment and edge turbulence evolution, the process of attachment to detachment is divided into five stages: the attached state, the pre-partially detached state, the partially detached state, the pre-pronounced detached state and the pronounced detached state. Plasma confinement barely varies in this process. The rest of this paper is as follows. In section 2, the experimental setup is described in detail. The experimental investigation of compatibility of confinement with detachment in L-mode plasma is presented in section 3, which includes observation of L-mode detachment, variation of plasma confinement and evolution of edge turbulent transport. The last two sections provide discussion and conclusions.

    Experiments were conducted in L-mode deuterium plasmas on the HL-2A tokamak, which has a major radius of R = 1.65 m, minor radius of α ≈ 0.4 m and toroidal magnetic field of Bt ≈ 1.3 - 2.7 T. Typical discharge parameters of the HL-2A tokamak can be found in [13]. Various diagnostics are employed in the experiments to investigate the process from attachment to detachment, part of which is shown in figure 1(a). Plasma radiation is estimated by two 16-channel arrays of absolute extreme ultraviolet (AXUV) detectors in the main chamber and the divertor regions with a sampling rate of 20 kHz. Passive visible spectroscopy can measure the line emission from impurity, covering almost the whole cross-section with a spatial resolution of 2.7 cm and a temporal resolution of 10 ms [14]. Meanwhile, high-speed vacuum ultraviolet (VUV) spectroscopy can measure the impurity line emission in the mid-plane with a time resolution of about 17 μs [15]. The Dα detector in the divertor has a time resolution of 1 μs. Two infrared (IR) cameras are installed at different toroidal positions to measure the heat flux onto the outer target. The target probe arrays consisting of 274 tips in total, are installed on both the inner and outer target plates [16], which have a spatial resolution of 1.2 cm and a temporal resolution of 1 μs.

    Figure  1.  (a) Skeleton of parts of the diagnostics employed in the main plasma chamber and divertor region and (b) a schematic diagram of a probe array with θ denoting the poloidal direction, φ denoting the toroidal direction and r denoting the radial direction.

    A fast reciprocating Langmuir probe array [17] is developed as the main tool to characterize the edge turbulence, which will be carefully discussed in this paper. It is located at the outer mid-plane at the low-field side, as shown in figure 1(a) and zoomed-in figure 1(b). It can move from the far SOL region to the edge region with a speed of about 1 m s−1, covering a radial distance of 8 cm. Radial profiles of plasma parameters can be obtained during discharge with a sampling rate of 1 MHz. The probe arrays consist of three probe heads, marked as Probes A, B and C in figure 1(b). Probe B is located at the mid-plane and the other two are separated poloidally with a central distance of 5 cm between the nearby probe heads. These three probe heads are also radially separated with a 0.5 cm distance between Probes B and C, and a 1 cm distance between B and A. Seven probe tips can be installed on each probe head, which are listed on the three radial steps with a step distance of 3 mm. Two of the seven tips on Probe C are left open in the experiments presented in this paper due to the limitation of the 19-tip vacuum socket. The height and diameter of each tip are 2.5 and 2 mm respectively. The poloidal separation of the adjacent tips is 5 mm. The two tips highlighted in green are to measure V- and V+, and the others are all designed to diagnose floating potential Vf The fluctuation of poloidal electric field ˜Eθ is derived from the poloidal floating potential difference neglecting electron temperature fluctuation, where symbol ~ denotes the fluctuation of parameters. The fluctuated radial velocity, turbulent particle flux and radial electric field are calculated as ˜Vr=˜Eθ/Bt, Γr=˜ne˜Vr and Er=rVp, respectively, where denotes a time ensemble average. In addition, the diagnostic system of Doppler backward scattering (DBS) is utilized to measure the relative density fluctuation as well as the poloidal velocity with a sampling rate of 1.25 MHz [18]. The 14-channel DBS measures plasma parameters in the plane that locates 13 cm lower than the mid-plane. The measured poloidal wave number kθ is in a range of 4–22 cm-1 which differs from the reciprocating probe arrays (kθ < 3 cm-1).

    In addition, gas is injected into the divertor chamber and the gas puffing signal denotes the voltage applied to solenoid valves on the pipes. The gas injection can be either continuous or in pulses, determined by the voltage applied to the valve. In the case of pulse injection, the duty cycle and the amplitude of the voltage applied to the valve control the amount of injected gas. In the case of continuous injection (which lasts for about 100–400 ms), the voltage is usually small to avoid too much gas puffing. The gas puffing in the inner and outer target can be set differently. It usually takes about 50–150 ms for the injected gas to move through the pipes into the divertor and the main chambers.

    In HL-2A L-mode plasma, detachment has been achieved in various ways, such as different impurity gases, toroidal injection locations and duration of pulses. In this paper, shot #38008 is selected as a typical shot of L-mode detachment. The main parameters of this shot are plasma current Ip165kA, toroidal magnetic field Bt1.3 T, NBI power PNBI300 kW, Ohmic heating voltage VL1.8 V, line-averaged density ¯ne(1.52.3)×1019 m3 with its ratio to Greenwald density limit ˉnenG0.370.57, edge safety factor qa3.1, and the stored energy We1218 kW for the duration of NBI heating.

    In this shot, mixed gas of N2(60% in volumn) +D2(60%) is injected into the divertor chamber in pulses. The pulse sequence starts at 1010 ms with different time intervals, and the injection location varies in the inner and outer divertor target, as shown in figure 2(b). The method of injection is designed to achieve different states of detachment and avoid excessive gas injection, which can lead to disruption. In the outer divertor, the duration of the four pulses is 20, 20, 15 and 15 ms, respectively, with an interval of 50 ms. In the inner divertor, the duration of the four pulses is 7, 7, 5 and 5 ms, respectively, with an interval of 50 ms. In figure 2, two detaching processes due to the first two gas pulses, highlighted in magenta, correspond to the time period of 1080–1135 and 1150–1270 ms. The two processes present similar performances and the second one is described in detail. With the increase in line-average density and radiation measured by VUV and AXUV in the main chamber, the radiation measured by AXUV and Dα in the divertor chamber decreases. The ion saturation current density measured by target probe arrays and the heat flux measured by the IR camera in the outer target show a dramatic reduction. In addition, the CCD camera pictures confirm stronger radiation close to the X-point after detachment in figure 3. In addition to the X-point, figure 3 shows that the visible radiation at the plasma boundary becomes stronger in the detached state than that in the attached state. These are all typical characteristics of detachment [2]. In fact, the time interval of detachment can be longer than the energy confinement time, which is about tens of milliseconds, but not ten times longer than energy confinement time due to many reasons, such as limited duration of NBI heating and discharge control in HL-2A in the last campaign. It is difficult to have slow evolution of the detachment states, which may result in transient rather than stationary regimes.

    Figure  2.  Time evolution of relevant plasma parameters in the main and divertor chambers, among which (a) is the NBI power PNBI and plasma current Ip, (b) the gas puffing in the inner and outer divertor, (c) the plasma radiation in the main chamber and in the divertor measured by AXUV, (d) the line emission of nitrogen measured by VUV, (e) the line-average density ¯ne and stored energy We (f) the Dα radiation in the divertor, (g) the heat flux on the outer target measured by the IR camera and (h) the ion saturation current density Jsat in the outer target measured by the target probe arrays. Two magenta boxes are the two processes of attachment to detachment.
    Figure  3.  CCD camera shows the distribution of plasma radiation in the main chamber (a) in the attached state (after puffing) and (b) in the pronounced detached state.

    Figure 4(a) illustrates the ion saturation current density Jsat measured by the target probe arrays on the inner divertor target, which shows an apparent decrease from 5 to about 2 A cm-2 when reaching detachment. The maximum electron temperature in the poloidal direction (Z in figure 4(b)) shows a sharp reduction from 20 to 3 eV. The heat flux measured by the IR camera decreases from 0.5 to about 0.2 MW m-2 which decreases by 60% in the outer target (figure 4(c)). These results from the divertor target demonstrate further that detachment has been realized in the HL-2A tokamak. In addition, it is found that the inner target detaches earlier than the outer target. The inner target starts detaching at about 1150 ms, while the outer target starts the process of detachment at about 1158 ms.

    Figure  4.  Spatial-temporal of (a) ion saturation current density by target probe arrays at the inner divertor target, (b) electron temperature from target probe arrays and (c) heat flux measured by the IR camera at the outer divertor target. Two magenta boxes are the two processes of attachment to detachment.

    Figure 5 presents the spatiotemporal information of NII line emission in the main chamber visible spectrometry. There is a pronounced peak around the time 1100 ms. Another NII emission peak shapes during the time of 1150–1200 ms and then decays gradually without keeping injecting N2 The interval between peaks is attributed to the pulse of injection. The NII emission concentrates in the lower half plane (-35 cm < Z < -15 cm), which is close to the X-point. In addition, this indicates a higher concentration level of nitrogen at 15 cm < Z < 35 cm in the upper half plane but a lower concentration in the mid-plane. The results are consistent with CCD camera measurements that the radiation near the X-point is much stronger in detachment.

    Figure  5.  Spatial-temporal profiles of N Ⅱ emission in log-scale in the main chamber shown by visible spectroscopy. Two magenta boxes are the two processes of attachment to detachment.

    As shown in figure 2(e), core confinement is not degraded after detachment on HL-2A L-mode plasma. In order to investigate how divertor detachment affects confinement, 32 shots with different gas species, gas injection methods and locations as well as various heating schemes are analyzed and listed in tables 1 and 2. In table 1, injected gas is either mixed gas (N2 + D2) or N2 only from the inner and/or outer target. Injection of other gases is not discussed here since the plasma behavior may be different. 'NW' means gas is injected only from the northwest sector, while 'NW + SO' means gas is injected from both the northwest and southeast sectors. 'Continuous' means gas is injected continuously for 100–400 ms. In table 2 the Ohmic heating power is not included since it is almost the same in different shots, which is about 300 kW (the plasma current is about 160 kA times, the loop voltage about 1.8 V). In order to obtain reliable results, the strategy is to make an average of stored energy of 50 ms (comparable to the energy confinement time) before detachment as the value of stored energy before in the attached state, and an average of stored energy within tens of milliseconds as the value of stored energy in the detachment state during each shot. Here, the detached state can be either the partially detached state or the pronounced detached state.

    Table  1.  Information of gas injection.
    Number of shots Gas species Injection ways Injection locations
    5 N2 + D2 pulse NW
    4 N2 continuous NW
    23 N2 continuous NW + SO
     | Show Table
    DownLoad: CSV
    Table  2.  Information of heating scheme.
    Number of shots NBI (kW) ECRH (kW) LHCD (kW) In total (kW)
    11 300–600 0 0 300–600
    3 440–550 670–980 0 1200–430
    9 220–620 0 370–900 590–1520
    8 200–600 370–1000 450–710 1400–220
    1 0 0 540 540
     | Show Table
    DownLoad: CSV

    Table 3 shows the increment of stored energy after detachment. In 75% of shots, the stored energy presents a low amplitude of increment (less than 10%); and in the other shots, the stored energy increases by 10%–20%. The average increment of stored energy after detachment is about 6% in these 32 shots. It can be concluded that the stored energy is not degraded but shows a slight increment after detachment in HL-2A L-mode plasma.

    Table  3.  Increment of stored energy after detachment in 32 shots.
    Increment Number of shots Ratio
    < 0 ~0 ~0%
    0–0.05 14 ~44%
    0.05–0.10 10 ~31%
    0.10–0.15 5 ~16%
    0.15–0.20 2 ~6%
    >0.2 0 ~0
    Average 6%
     | Show Table
    DownLoad: CSV

    As illustrated above, plasma radiation in the main chamber increases to a great extent (figures 2(c) and 4), while the total energy remains the same or increases a little (figure 2(e)). There must be a change of some other physical quantities to conserve the energy. Since the plasma profiles, especially edge profiles, change significantly during detachment, edge turbulence and transport, which are critical to plasma confinement, are investigated in the following sections. Since turbulence in plasma affects confinement, it is necessary to investigate the evolution of edge turbulent transport during detachment.

    Turbulence and induced transport near the last closed flux surface (LCFS) before and after detachment in the outer mid-plane are carefully analyzed. The Langmuir probe stays at a fixed position, which is very close to the LCFS during a time interval from 1140 to 1200 ms. Plasma performance is similar in process from attachment to detachment, so we take the second process in typical shot #38008 as an example. As shown in figure 6, the process of attachment to detachment is divided into five elaborate states according to plasma parameters on the divertor target and edge turbulent transport: the attached state (AB zone), the pre-partially detached state (BC), the partially detached state (CD), the pre-pronounced detached state (DE) and the pronounced detached state (EF). In the two pre-state states, plasma parameters evolve rapidly in the outer mid-plane. The classification of each state in this paper is identified according to plasma parameter evolution in the divertor as well as in the mid-plane, which is different to the classification in AUG [19, 20]. The phenomenology reported here does not correspond to what is typically performed in detachment experiments.

    Figure  6.  Time evolutions of (a) Jsat in the outer target, (b) ion saturation current Is, (c) electron temperature Te (d) electron density ne, (e) electron pressure pe measured near the LCFS by Langmuir probes in the outer mid-plane.

    In the attached state, Dα radiation in the divertor region is strong and Jsat is about 4 A cm−2 in the outer divertor target. The plasma parameters, such as the electron density ne and electron temperature Te and their radial gradients, are stable. In the pre-partially detached state, the inner and outer targets are still in the attached state. In addition to ne, which increases slightly, Te and electron pressure pe and their radial gradients increase rapidly. In the partially detached state, Jsat in the outer divertor target usually decreases and then remains at a low amplitude, while edge plasma parameters in the outer mid-plane are maintained at high levels. The Dα signal shows the same tendency as Jsat. It is noted that it is the sawtooth (t ~ 1148 ms) that causes a peak in the signals. Meanwhile, the radiation close to the X-point becomes stronger. The pre-pronounced detachment states are characterized by a dramatic reduction in edge plasma parameters in the outer mid-plane in a timescale of about 15 ms, particularly Te, which decreases from about 46 to 12 eV, by about 80%. The ne in the edge shows a smaller reduction since the line-average density increases rapidly, as shown in figure 2(e). An obvious and interesting fact is that the gradients of these parameters become smaller and smaller, and at the end of this pre-pronounced detached state they all approach 0, as shown in figures 7(a)(c). In this state, Dα in the divertor chamber and Jsat in the outer target decrease to very low amplitudes. In the pronounced detached state, Dα in the divertor chamber and Jsat in the outer target are maintained at very low amplitudes, while visible radiation is very strong near the X-point and it becomes stronger at the plasma boundary, as shown in figure 4. The strong radiation cooled the plasma boundary, resulting in the main plasma parameters, especially Te and pe, being at very low amplitudes. Their gradients are all near 0, implying a very flat edge plasma profile.

    Figure  7.  Time evolutions of the radial gradient of (a) electron density ∇rne, (b) electron temperature ∇rTe, (c) electron pressure ∇rpe, (d) poloidal velocity by time delay estimate (TDE) Vθ, TDE, and (e) the radial electric field Er measured near the LCFS by Langmuir probes in the outer mid-plane.

    The physics of the plasma parameter evolutions in the outer mid-plane before pronounced detachment can be understood by means of examining edge radial electric field Er and poloidal velocity Vθ, which are almost synchronously evolving. Er remains at −2000 V m−1 in the attached state, and then varies rapidly to about 0 V m−1 in the pre-partially detached state and keeps increasing to as large as 4000 V m−1 in the partially detached state, which results in a strong poloidal velocity Vθ shear layer around the LCFS. Er decreases from 3500 V m−1 sharply to a few hundred V m−1 in the pre-pronounced state and remains at a very low amplitude in the pronounced detached state. The evolutions of Er and Vθ imply the dramatic change of edge turbulence and flow shear. More information about edge turbulence and transport is required, which necessitates further study.

    Edge turbulence and transport play an important role in plasma confinement. Understanding how edge turbulence transport evolves and its relation to detachment will be helpful for detachment operation for tokamaks. Figure 8 shows the auto-spectra of ˜Is in the attached, partially detached and pronounced detached states around the LCFS measured by Langmuir probe arrays in the outer mid-plane. The auto-spectra of ˜Is in the pre-partial and pre-pronounced detached states are not discussed here because plasma parameters evolve rapidly during both states. According to the auto-spectra, the characterized bands of turbulence are in the frequency range of about 20–100 kHz. It is clear that turbulence intensity increases in the partially detached state. When the plasma reaches the pronounced detached state, the turbulence intensity decreases dramatically, but is much smaller than that of the attached state.

    Figure  8.  Auto-spectra of ˜Is in the attached, partially detached and pronounced detached states around the LCFS measured by Langmuir probe arrays in the outer mid-plane.

    The normalized probability distribution functions (PDFs) of fluctuated radial velocity fluctuations and electron density fluctuations are shown in figure 9. One can see that the PDFs of radial velocity fluctuations are very similar in the attached and partially detached states, lying in a range from −3000 to 3000 m s−1, while they are concentrated in a range from −1500 to 1500 m s−1 in the pronounced detached state. It is clear that the width of half height of the PDF in the pronounced state is much smaller than those in the attached and partially detached states. The kurtosis values are 3.6, 3.9 and 4.5 in the attached, partially detached and pronounced detached states, respectively. The larger kurtosis value agrees with the smaller width of half height of the PDF in the pronounced detached state, indicating the strong reduction in radial velocity fluctuations in this state. The values of skewness of radial velocity fluctuations in the three states are about −0.4, 0 and −0.2, respectively. Negative skewness value means that the PDF is negatively deviated. The PDF of fluctuated electron density varies slightly. Its skewness values are 0.3, 0.5 and 0.6 in the three states, respectively, meaning that distribution of density fluctuation is positively deviated and increases slightly.

    Figure  9.  Normalized PDFs of the (a) radial velocity fluctuation ˜Vr and (b) electron density fluctuation ˜ne.

    Turbulent particle flux is studied in detail in figure 10. ˜nrmse denotes the root mean square (rms) of density fluctuation, and ˜Vrmsr represents the rms of fluctuated radial velocity. The correlation coefficient between the density and radial velocity fluctuation is about 0.75–0.78, which varies slightly from the attached to the pronounced detached state. The cross-phase of fluctuated density and radial velocity is shown in figure 10(c). Density fluctuation varies moderately in the whole process. The radial velocity fluctuation level is about 1200 m s−1 in the attached state, and dramatically increases to 1800 m s−1 in the pre-partially detached state, and then remains at about 1500 m s−1 in the partially detached state, followed by a sharp decrease to about 800 m s−1 in the pre-pronounced detached state, and finally remains at a magnitude of a few hundred m s−1 in the pronounced detached state. The cross-phase is about negative 10° in the attached and pre-partially detached states, and then rises to about positive 20° and maintains this in the partially detached state. It decreases to about 0 in the pre-pronounced detached state and remains the same in the pronounced state.

    Figure  10.  Time evolution of the rms of (a) electron density fluctuation and (b) radial velocity fluctuation, (c) cross-phase between fluctuated density and radial velocity and (4) turbulent particle flux in the process from attached to pronounced detached state.

    In the attached state, the time-averaged turbulent particle flux is very small, due to the negative cross-phase between density and radial velocity fluctuation, which may be convective cells locating around the LCFS. Turbulent particle flux increases dramatically in the pre-partially detached state, mainly due to the sharp increase in turbulent radial velocity. However, the radial electric field is very small in this state. Turbulent particle flux decreases apparently and then remains in the partially detached state, which is consistent with the strong E × B shear in this state. In the pre-pronounced state, turbulent particle flux increases at first, corresponding to both the increase in density fluctuation and reduced Er; then turbulent particle flux decreases mainly due to the decreased radial fluctuated velocity as well as contributions from the cross-phase between electron density fluctuation and radial velocity fluctuation. The evolution of turbulent particle transport is consistent with the variation in Er. In the pronounced detached state, due to the redistribution of plasma profiles in the edge, radial gradients of electron temperature and electron pressure are approaching 0. This means that the free energy source for turbulence is very small, which is a major reason for the very low intensity of turbulence and transport in the pronounced detached state.

    Figure 11 shows the time-frequency spectra of the fluctuated density, radial velocity and turbulent particle flux in log-scale. It can be seen that the main frequency band is in the range of 20–100 kHz, and decays quickly in the above frequency band before the pronounced detached state. However, turbulence above 30 kHz is strongly suppressed in the pre-pronounced detached state. The major frequency band and intensity of turbulence of radial velocity fluctuation vary obviously in different states. The frequency band becomes broader in the pre-partially detached state and then becomes a little smaller in the partially detached state, and then decreases sharply in the pre-detached state. It decreases further in the pronounced state, whose main turbulence frequency band can hardly be seen. The evolution tendency of turbulent particle flux is similar to that of the fluctuated radial velocity. Here, the spectral function of the turbulent particle flux is the real part of the cross-spectral density (CSD) of density fluctuation and radial velocity fluctuation, as expressed in [21]. Hence, it is concluded that with the higher frequency band (f > 30 kHz) of fluctuated density, the whole frequency bands of fluctuated radial velocity and turbulent particle flux are strongly suppressed in the pronounced detached state. The results shown in figures 811 are consistent in terms of edge turbulence and transport.

    Figure  11.  (a) Time-frequency PSD of the density fluctuation, (b) time-frequency PSD of radial velocity fluctuation and (c) spectral function of turbulent particle flux.

    The relative density fluctuation and poloidal velocity are also measured by the DBS near the LCFS, as shown in figure 12. In figure 12(a), the intensity of relative density fluctuation increases in the partially detached state and decreases in the pronounced state, which agrees with the results in figure 8 measured by Langmuir probe arrays. The poloidal velocity increases in the partially detached state and approaches 0 in the pronounced detached state, which is also consistent with the results of the poloidal velocity measured by Langmuir probe arrays at the outer mid-plane.

    Figure  12.  (a) Relative density fluctuations and (b) spectrogram of the Doppler backscattering signal measured at the edge plasma together with the derived poloidal velocity. Poloidal velocity is obtained from the estimated Doppler shift (black line) times a constant 0.014.

    The reason for non-degradation of plasma confinement under pronounced detachment may be interpreted as follows: although the radiation in the main plasma and especially in the X-point are strong, the intensity of edge turbulent transport is very weak due to flat plasma profiles at the edge; electron density in the core is increasing; electron temperature does not change much in the core plasma. As a consequence, the confinement does not degrade in HL-2A L-mode detached plasma. The effects of the closed divertor configuration on the detachment and edge turbulence transport require further study.

    The change in the radial electric field Er in the pronounced detached state is discussed here. As stated in [22], the general tendency is that the better confinement is coupled with more negative Er and lower edge turbulent transport in H-mode plasma. However, in this work, Er is approaching zero in the pronounced state, while the plasma confinement is not degraded simultaneously mainly due to lower edge turbulent transport in detached L-mode plasma in HL-2A. The mechanisms behind the reduced edge turbulent transport in general H-mode plasma (due to stronger shear flow) and the detached L-mode plasma in HL-2A (due to edge cooling via radiation) are different. In addition to turbulent transport, one more significant channel of power exhaust radiation should be considered for confinement in the detached plasma. The fact that Er is approaching zero in the detached L-mode plasma measured by Langmuir probe can be validated by the DBS data in figure 12, where a black line denotes the poloidal velocity approaching zero. This result is different to the general tendency in [22].

    It is not very clear why the electron temperature and pressure increase rapidly in the pre-partially detached state. However, the soft X-ray signals are analyzed, and it is found that signals in the center channels decrease while signals in the outside channels increase. This may imply that the plasma energy may propagate from the core to the plasma edge in this case. There are similar observations from other shots, including the shots that have continuous injection up to 300 ms. This may contribute to the increment in electron temperature in the pre-partially detached state, but this requires further study in the near future.

    The paper presents the characteristics of L-mode detachment, together with its induced evolutions of edge turbulent transport and plasma confinement in NBI-heated plasma in the HL-2A tokamak. Partially detached and pronounced detached states have been achieved in L-mode plasma. The evolution of plasma parameters and edge turbulence in the outer mid-plane during detachment are carefully discussed to understand the detachment process. The stored energy is maintained in the process from attachment to pronounced detachment. In the partially detached state, edge turbulence and its transport have increased obviously compared to the attached state. In the pronounced detachment state, the redistribution of the plasma profiles due to detachment leads to low amplitude of electron temperature and pressure, as well as very weak edge turbulence and transport. Despite the strong plasma radiation in the pronounced detached state, the reduced edge turbulent transport contributes to the retention of stored energy in detached L-mode plasma in HL-2A. Detachment plays a significant role in the redistribution of plasma profiles, which should be studied further in each detached state.

    This work was supported by National Natural Science Foundation of China (Nos. 12235002 and 12122503), National Key R&D Program of China (No. 2018YFE0301101), Dalian Science & Technology Talents Program (No. 2022RJ11), and Xingliao Talent Project (No. XLYC2203182).

  • [1]
    Pitts R A et al 2019 Nucl. Mater. Energy 20 100696 doi: 10.1016/j.nme.2019.100696
    [2]
    Ryutov D D and Soukhanovskii V A 2015 Phys. Plasmas 22 110901 doi: 10.1063/1.4935115
    [3]
    Zhang Y J et al 2020 Nucl. Mater. Energy 25 100803 doi: 10.1016/j.nme.2020.100803
    [4]
    Zheng G Y et al 2016 Nucl. Fusion 56 126013 doi: 10.1088/0029-5515/56/12/126013
    [5]
    Stangeby P C and Elder J D 1992 J. Nucl. Mater. 196–198 258 doi: 10.1016/S0022-3115(06)80042-5
    [6]
    Wang J et al 2023 Phys. Plasmas 30 043905 doi: 10.1063/5.0137801
    [7]
    Nichols J H et al 2021 Nucl. Fusion 61 096018 doi: 10.1088/1741-4326/ac14e6
    [8]
    Guo J et al 2021 Plasma Phys. Control. Fusion 63 125006 doi: 10.1088/1361-6587/ac2afa
    [9]
    Xu G L et al 2022 Nucl. Mater. Energy 33 101259 doi: 10.1016/j.nme.2022.101259
    [10]
    Kumpulainen H A et al 2020 Nucl. Mater. Energy 25 100784 doi: 10.1016/j.nme.2020.100784
    [11]
    Wang H et al 2022 Nucl. Fusion 62 126018 doi: 10.1088/1741-4326/ac8fa4
    [12]
    Ma X et al 2020 Phys. Scr. T171 014072 doi: 10.1088/1402-4896/ab4a39
    [13]
    Sinclair G et al 2022 Nucl. Fusion 62 106024 doi: 10.1088/1741-4326/ac8b95
    [14]
    Abrams T et al 2021 Phys. Scr. 96 124073 doi: 10.1088/1402-4896/ac3c5f
    [15]
    Maeker J B et al 2022 Nucl. Mater. Energy 33 101309 doi: 10.1016/j.nme.2022.101309
    [16]
    Zamperini S A et al 2022 Nucl. Fusion 62 026037 doi: 10.1088/1741-4326/ac3fe7
    [17]
    Zamperini S et al 2023 Fusion Sci. Technol. 79 36 doi: 10.1080/15361055.2022.2082791
    [18]
    Gao S L et al 2021 AIP Adv. 11 025233 doi: 10.1063/5.0037381
    [19]
    Gao S L et al 2022 Plasma Sci. Technol. 24 075104 doi: 10.1088/2058-6272/ac608f
    [20]
    Wang Y L et al 2021 Plasma Phys. Control. Fusion 63 085002 doi: 10.1088/1361-6587/ac0351
    [21]
    Wang Y L et al 2023 Nucl. Fusion 63 096024 doi: 10.1088/1741-4326/aceb09
    [22]
    Gallo A et al 2020 Phys. Scr. T171 014013 doi: 10.1088/1402-4896/ab4308
    [23]
    Ciraolo G et al 2021 Nucl. Fusion 61 126015 doi: 10.1088/1741-4326/ac2439
    [24]
    Di Genova S et al 2021 Nucl. Fusion 61 106019 doi: 10.1088/1741-4326/ac2026
    [25]
    Schneider R et al 2006 Contrib. Plasma Phys. 46 3 doi: 10.1002/ctpp.200610001
    [26]
    Makkonen T et al 2011 J. Nucl. Mater. 415 S479 doi: 10.1016/j.jnucmat.2010.08.023
    [27]
    Elder J D et al 2005 J. Nucl. Mater. 337–339 79 doi: 10.1016/j.jnucmat.2004.10.138
    [28]
    Xu Y C et al 2021 Plasma Phys. Control. Fusion 63 095003 doi: 10.1088/1361-6587/ac11b3
    [29]
    Zhao X L et al 2022 Nucl. Fusion 62 126071 doi: 10.1088/1741-4326/ac9b77
    [30]
    Zhang Y J et al 2022 Nucl. Fusion 62 106006 doi: 10.1088/1741-4326/ac8564
    [31]
    Petrie T W et al 2008 Nucl. Fusion 48 045010 doi: 10.1088/0029-5515/48/4/045010
    [32]
    Warrier M, Schneider R and Bonnin X 2004 Comput. Phys. Commun. 160 46 doi: 10.1016/j.cpc.2004.02.011
    [33]
    Cupak C et al 2021 Appl. Surf. Sci. 570 151204 doi: 10.1016/j.apsusc.2021.151204
    [34]
    Ding R et al 2016 Nucl. Fusion 56 016021 doi: 10.1088/0029-5515/56/1/016021
    [35]
    Krieger K et al 1999 J. Nucl. Mater. 266–269 207 doi: 10.1016/S0022-3115(98)00890-3
    [36]
    Schmid K et al 2010 Nucl. Fusion 50 105004 doi: 10.1088/0029-5515/50/10/105004
    [37]
    Zhou Q R et al 2020 Nucl. Mater. Energy 25 100849 doi: 10.1016/j.nme.2020.100849
    [38]
    Wang L et al 2022 Nucl. Fusion 62 076002 doi: 10.1088/1741-4326/ac4774
    [39]
    Hitzler F et al 2020 Plasma Phys. Control. Fusion 62 085013 doi: 10.1088/1361-6587/ab9b00
    [40]
    Zhao X L et al 2020 Plasma Phys. Control. Fusion 62 055015 doi: 10.1088/1361-6587/ab831b
    [41]
    Guo J et al 2023 Nucl. Fusion 63 126033 doi: 10.1088/1741-4326/ad00cc
    [42]
    Stangeby P C 2000 The Plasma Boundary of Magnetic Fusion Devices (Boca Raton: CRC Press
    [43]
    Wang H et al 2024 Nucl. Fusion 64 046009 doi: 10.1088/1741-4326/ad29bc
    [44]
    Dux R et al 2011 Nucl. Fusion 51 119501 doi: 10.1088/0029-5515/51/11/119501
  • Related Articles

    [1]X T DING (丁玄同), W CHEN (陈伟). Review of the experiments for energetic particle physics on HL-2A[J]. Plasma Science and Technology, 2018, 20(9): 94008-094008. DOI: 10.1088/2058-6272/aad27a
    [2]Zhongbing SHI (石中兵), Wulyu ZHONG (钟武律), Min JIANG (蒋敏). Progress of microwave diagnostics development on the HL-2A tokamak[J]. Plasma Science and Technology, 2018, 20(9): 94007-094007. DOI: 10.1088/2058-6272/aad27b
    [3]Wulyu ZHONG (钟武律), Xiaolan ZOU (邹晓岚), Zhongbing SHI (石中兵), Xuru DUAN (段旭如), Min XU (许敏), Zengchen YANG (杨曾辰), Peiwan SHI (施培万), Min JIANG (蒋敏), Guoliang XIAO (肖国梁), Xianming SONG (宋显明), Jiaqi DONG (董家齐), Xuantong DING (丁玄同), Yong LIU (刘永), HL-A team (HL-A团队). Dynamics of oscillatory plasma flows prior to the H-mode in the HL-2A tokamak[J]. Plasma Science and Technology, 2017, 19(7): 70501-070501. DOI: 10.1088/2058-6272/aa6538
    [4]JIANG Chunyu (蒋春雨), CAO Jing (曹靖), JIANG Xiaofei (蒋小菲), ZHAO Yanfeng (赵艳凤), SONG Xianying (宋先瑛), YIN Zejie (阴泽杰). Real-Time Bonner Sphere Spectrometry on the HL-2A Tokamak[J]. Plasma Science and Technology, 2016, 18(6): 699-702. DOI: 10.1088/1009-0630/18/6/19
    [5]XIA Donghui(夏冬辉), ZHOU Jun(周俊), RAO Jun(饶军), HUANG Mei(黄梅), LU Zhihong(陆志鸿), WANG He(王贺), CHEN Gangyu(陈罡宇), WANG Chao(王超), LU Bo(卢波), ZHUANG Ge(庄革). Design of the Transmission Lines for 140 GHz ECRH System on HL-2A[J]. Plasma Science and Technology, 2014, 16(3): 267-272. DOI: 10.1088/1009-0630/16/3/17
    [6]WANG Chao (王超), ZHOU Jun (周俊), HUANG Mei (黄梅), WANG He (王贺), CHEN Gangyu (陈罡宇), RAO Jun (饶军). ECRH Launcher for Four-Beam Injection on HL-2A Tokamak[J]. Plasma Science and Technology, 2013, 15(5): 476-479. DOI: 10.1088/1009-0630/15/5/16
    [7]HUANG Xianli (黄贤礼), SHI Zhongbing (石中兵), CUI Zhengying (崔正英), ZHONG Wulv (钟武律), DONG Yunbo (董云波), CHEN Chengyuan (陈程远), FENG Beibin (冯北滨), YAO Lianghua (姚良骅), LIU Zetian (刘泽田), DING Xuantong (丁玄同), et al. Heat Transport During H-Mode in the HL-2A Tokamak[J]. Plasma Science and Technology, 2013, 15(3): 221-224. DOI: 10.1088/1009-0630/15/3/06
    [8]CHEN Wenguang (陈文光), RAO Jun (饶军), LI Bo (李波), LEI Guangjiu (雷光玖), CAO Jianyong (曹建勇), WANG Mingwei (王明伟), KANG Zihua (康自华), FENG Kun (冯鲲), HL-A NBI Group. Technical Design of Arc-Discharge and Deceleration Power Supply for MW Level NBI System on HL-2A Tokamak[J]. Plasma Science and Technology, 2012, 14(10): 936-940. DOI: 10.1088/1009-0630/14/10/15
    [9]HE Zhixiong, DONG Jiaqi, HE Hongda, JIANG Haibin, GAO Zhe, ZHANG Jinhua. MHD Equilibrium Configuration Reconstructions for HL-2A Tokamak[J]. Plasma Science and Technology, 2011, 13(4): 424-430.
    [10]YAO Lianghua (姚良骅), ZHAO Dawei (赵大为), FENG Beibin (冯北滨), CHEN.Chengyuan (陈程远), ZHOU Yan(周艳), HAN Xiaoyu (韩晓玉), LI Yonggao (李永高), Jerome BUCALOSSI, Duan Xuru (段旭如). Comparison of Supersonic Molecular Beam Injection from both low field side and high field side of HL-2A[J]. Plasma Science and Technology, 2010, 12(5): 529-534.
  • Cited by

    Periodical cited type(6)

    1. Wu, T., Xu, M., Wang, Z. et al. Compatibility of pronounced detachment with improved confinement on HL-2A tokamak. Nuclear Fusion, 2025, 65(2): 026022. DOI:10.1088/1741-4326/ad9e04
    2. Hao, G., Xu, J., Sun, Y. et al. Summary of the 11th Conference on Magnetic Confined Fusion Theory and Simulation. Plasma Science and Technology, 2024, 26(10): 101001. DOI:10.1088/2058-6272/ad5d8a
    3. Wu, N., Cheng, J., Yi, K.Y. et al. Facilitated core-edge integration through divertor nitrogen seeding in the HL-2A tokamak. Nuclear Fusion, 2024, 64(9): 096007. DOI:10.1088/1741-4326/ad5e96
    4. Zhou, Y., Dudson, B., Wu, T. et al. One-dimensional simulation and validation of divertor detachment induced through nitrogen seeding on HL-2A. Plasma Physics and Controlled Fusion, 2024, 66(5): 055005. DOI:10.1088/1361-6587/ad30fc
    5. Long, T., Ke, R., Wu, T. et al. Studies of edge poloidal rotation and turbulence momentum transport during divertor detachment on HL-2A tokamak | [HL-2A 托卡马克偏滤器脱靶时边缘极向旋转和湍流动量输运]. Wuli Xuebao/Acta Physica Sinica, 2024, 73(8): 088901. DOI:10.7498/aps.73.20231749
    6. Wu, T., Diamond, P.H., Nie, L. et al. How turbulent transport broadens the heat flux width: local SOL production or edge turbulence spreading?. Nuclear Fusion, 2023, 63(12): 126001. DOI:10.1088/1741-4326/acf5d9

    Other cited types(0)

Catalog

    Figures(12)  /  Tables(2)

    Article views (24) PDF downloads (3) Cited by(6)

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return