Processing math: 100%
Advanced Search+
Jing OU, Jiamin LONG. Estimates of required impurity fraction for EAST divertor detachment[J]. Plasma Science and Technology, 2025, 27(1): 015103. DOI: 10.1088/2058-6272/ad8ffb
Citation: Jing OU, Jiamin LONG. Estimates of required impurity fraction for EAST divertor detachment[J]. Plasma Science and Technology, 2025, 27(1): 015103. DOI: 10.1088/2058-6272/ad8ffb

Estimates of required impurity fraction for EAST divertor detachment

More Information
  • Author Bio:

    Jing OU: ouj@ipp.ac.cn

  • Corresponding author:

    Jing OU, ouj@ipp.ac.cn

  • Received Date: June 02, 2024
  • Revised Date: October 31, 2024
  • Accepted Date: November 06, 2024
  • Available Online: November 07, 2024
  • Published Date: December 30, 2024
  • During the EAST radiative divertor experiments, one of the key challenges was how to avoid the occurrence of disruptive events caused by excessive impurity seeding. To estimate the required impurity fraction for divertor detachment, we introduce a reduced edge plasma radiation model. In the model, based on the momentum conservation along the magnetic field line, the upstream pressure is determined by the plasma density and temperature at the divertor target, and then the impurity radiation loss is obtained by the balance of the heat and particle fluxes. It is found that the required impurity fraction shows a non-monotonic variation with divertor electron temperature (Td) when 0.1eV<Td<10eV. In the range of 0.1eV<Td<1eV, the position near the valley of required impurity fraction corresponds to strong plasma recombination. Due to the dependence of the volumetric momentum loss effect on the Td in the range of 1eV<Td<10eV, the required impurity fraction peaks and then decreases as Td is increased. Compared to neon, the usage of argon reduces the impurity fraction by about twice. In addition, for the various fitting parameters in the pressure–momentum loss model, it is shown that the tendency of required impurity fraction with Td always increases first and then decreases in the range of 1eV<Td<10eV, but the required impurity fraction decreases when the model that characterizes the strong loss in pressure momentum is used.

  • Neutral beam injection is a main auxiliary heating method with the clearest physical mechanism and the highest heating efficiency. It has been widely used in tokamak and other types of magnetic confinement fusion devices. The Experimental Advanced Superconducting Tokamak (EAST) is equipped with two beam lines, each equipped with two sets of 80 keV/50 A hot cathode bucket positive ion sources [13]. The designed injection power is 4 MW per beamline, and the total injection power can reach 8 MW [4, 5]. In order to support the physical research on the EAST, a new positive ion source with designed beam energy of 120 keV was proposed to be developed. Accelerator structure is one of the key components of the ion source. When the ion source is upgraded to 120 keV, the voltage applied between different grids is greatly increased. It is necessary to redesign a reasonable insulation support system to ensure the stability of the accelerator. The accelerator insulator must be made of materials with good sealing performance, high strength and good insulation performance. In 1981, Paterson et al designed the insulator structure of positive ion source accelerator of TFTR device, and gave the manufacturing process of Al2O3 ceramic titanium brazing rectangular insulator structure [6]. In 2008, Yamano et al analyzed and compared the vacuum discharge characteristics of fiber reinforced polymer (FRP) and Al2O3 ceramic insulator samples in JT-60U accelerator, and found that the flashover voltage of FRP samples was lower than that of Al2O3 ceramic samples [7]. In 2001, the Korean Atomic Energy Research Institute designed a round-hole type accelerator for the KSTAR neutral beam injection device at 120 keV, and determined the PEEK material insulator and its parameters [8]. Based on the above research, at present, in the engineering application of ion source, the relatively mature insulating materials that can meet the requirements of all aspects include FRP fiber-reinforced epoxy resin, PEEK and Al2O3 ceramics.

    In addition, scholars have carried out various studies on the insulation support system of ion sources in major fusion devices. In 2008, Ikeda et al measured the breakdown voltage of FRP insulator with and without grids in order to test the voltage holding capability of the 500 kV, 22 A on the three-grids accelerator in JT-60U [9]. In 2010, Tanaka et al tested the breakdown voltage at different gap lengths by changing the grid gap of the JT-60 negative ion source in order to determine the required grid gap length for stable voltage holding [10]. In 2014, Wei et al investigated the connection design of the insulator of the EAST ion source accelerator. Because the bonding technology has the advantages of simple operation and high strength, the flange and insulator were bonded by using epoxy adhesive, and then different samples were tested by tensile test, so as to determine the optimal bonding structure [11]. Based on the above research, scholars mainly focus on the electric field distribution around the grids and insulators of accelerator. In this paper, the design of the insulation support system of the 120 keV accelerator is introduced firstly, and then based on the structural model of the 120 keV positive ion source accelerator, the electrostatic analyses of the insulators and grids are carried out by using the method of finite element analysis. The voltage resistance characteristics of the insulators are analyzed, and the appropriate insulator material and corresponding design parameters are determined.

    The 80 keV hot cathode bucket positive ion source operated in EAST is composed of a bucket plasma generator with a multi-magnet pole field and a four-grids slit-type accelerator [4]. The accelerator is mainly composed of flanges, insulators, grid supports and four grids, including plasma grid (PG), gradient grid (GG), suppression grid (SG) and earth grid (EG). The insulation support system is made of three layers of FRP insulators and four layers of 304 stainless steel flanges bonded by epoxy adhesive, which serves to support the whole beam source system and voltage insulation.

    A key issue in the development of EAST ion source is the connection between the insulator and the support flange. The insulation support system needs to meet the requirements of high voltage holding, mechanical strength and vacuum sealing. Although the epoxy adhesive bonding technology used in the 80 keV positive ion source accelerator has the advantages of simple operation and low cost, it is very difficult to repair when the insulation support system leaks under long-term operation. At the same time, when the 80 keV positive ion source accelerator is optimized and upgraded, the disassembly device is also a big problem. Therefore, the insulation support system of 120 keV positive ion source accelerator was optimized and upgraded in these aspects. The 120 keV positive ion source uses a four-grids circular hole accelerator. The insulators and flanges of the 120 keV positive ion source accelerator are fixedly connected by M8 PEEK hex socket head cap screws, and sealed by sealing rings, which is more convenient for disassembly and maintenance. In addition, the 120 keV positive ion source accelerator uses PEEK material with better insulation performance than the FRP material used in the insulator of the 80 keV positive ion source accelerator. The 120 keV positive ion source accelerator has a length of 856 mm, a width of 470 mm and a height of 355 mm. The four layers of grids will be applied with voltages of 120 kV, 93.6 kV, −5 kV and 0 kV respectively. The grid plate of the 120 keV positive ion source accelerator is made of molybdenum material, which has higher hardness and melting point than the copper material used in the grid plate of the 80 keV positive ion source accelerator. The specific structural parameters are shown in the figures 1 and 2.

    Figure  1.  Schematic map of the 120 keV positive ion source.
    Figure  2.  The EG-SG insulator: (a) vertical view, (b) sectional view.

    The 80 keV positive ion source has breakdown and leakage issues during operation in the past few years, which may be due to the breakdown of vacuum gap between the accelerating grid and the grid support structure, as well as vacuum leakage near the insulators. In fact, obvious discharge marks were observed in these areas during the disassembly of the accelerator. Therefore, after completing the design phase, it is necessary to evaluate the insulation characteristics of the 120 keV accelerator and select appropriate insulation materials. There are currently three main types of insulation materials, namely FRP fiber-reinforced epoxy resin, PEEK, and Al2O3 ceramics.

    Based on the three-dimensional structural model of the 120 keV positive ion source accelerator, the study of the voltage holding performance of insulators with different heights, structures and materials was carried out by electrostatic simulations with COMSOL Multiphysics finite element software [12]. In order to ensure a certain safety margin, the voltage difference between grids was increased by 20% during simulations respectively, and the voltages of four layers of grids and flanges were set to be 144 kV (PG and flange), 112.32 kV (GG and flange), −6 kV (SG and flange), and 0 kV (EG and flange). The design heights of the three insulators were 55.00 mm (PG-GG insulator), 172.70 mm (GG-SG insulator) and 33.88 mm (SG-EG insulator). Under the condition of steady-state vacuum operation, different insulating materials were simulated by changing the relative dielectric constants of materials in the software, and then the changes of the distribution of electric fields around the insulators were observed by changing the height of the insulators.

    In order to save computing time, the accelerator was divided into three parts for electrostatic simulations. As shown in figure 3, the region near the insulators of each layer (blue-boxed region) is selected as the study object, in which the electric field distribution is most obviously affected by the two variables of insulator material and height, and then the maximum electric field (red-boxed region) is searched for in this region. As shown in figure 3(a), since the water-cooling tubes on GG flange are located between PG flange and GG flange, a layer of stainless steel cover plate is added on the surface of GG water-cooling tubes for protection. It can be clearly observed that near the PG-GG insulator, the electric field is the largest at the boundary of the GG flange water pipe protection cover plate, and when the height of the PG-GG insulator is 37 mm, which corresponds to a minimum gap of 9.15 mm in this region, the maximum electric field at this location is 7.4 kV/mm. Since the probability of breakdown increases significantly when the minimum gap between flanges is smaller than the grid gap, we started by increasing the height of the insulator corresponding to the minimum gap between flanges and the grid gap, and continued to increase until the design value is reached. We recorded the maximum electric field in the region at different heights. It can be seen from figure 4: (1) with the increasing of the insulator height, the maximum electric field on the surface decreases gradually, and the two are negatively correlated; (2) at the same insulator height, the PEEK material has the smallest maximum electric field and the best insulation performance among the three materials; (3) when the height of the insulator reaches the design value, the maximum electric field on the surface of PG-GG insulator and GG-SG insulator is around 4 kV/mm.

    Figure  3.  Distribution diagrams of electric fields around insulators: (a) PG-GG flanges and grids, (b) GG-SG flanges and grids, (c) EG-SG flanges and grids.
    Figure  4.  Plots of surface maximum electric fields as a function of the height of the insulator: (a) PG-GG insulator, (b) GG-SG insulator, (c) EG-SG insulator.

    In summary, PEEK material was ultimately chosen as the insulation material. Moreover, large ceramics are difficult to process, with poor dimensional accuracy and high weight. There are also some defects in the finished product such as cracks, porosity, brittle fracture, etc. PEEK material not only has better insulating properties, but also has better machinability and surface than FRP epoxy resin. The design height of the insulator can meet the withstand voltage requirement of 120 keV, at which time the maximum electric field in the insulator region between each flange is reduced to about 4 kV/mm or less than 4 kV/mm [13].

    During the four-year operation period from 1991 to 1994, 6% of the NBI system failure rate of the JT-60 device occurred on the ion source structure, which was mainly caused by the leakage of the cooling channels of the accelerator grids, the vacuum leakage of the insulators, and the voltage breakdown between the grids [14]. Therefore, the electric field concentration between grids is also an important factor affecting the insulation performance of accelerator. Through the beam optical analysis, the grid gap values were set to 3 mm (PG-GG), 10.5 mm (GG-SG), 3 mm (SG-EG), and the voltage values of four layers of grids and flanges were set to 120 kV (PG and flange), 93.6 kV (GG and flange), −5 kV (SG and flange), and 0 kV (EG and flange). Due to the small voltage difference between the SG and the EG, the maximum electric field was within the safe range. Therefore, the main simulations were conducted on the electric field distribution between the PG-GG and GG-SG with a large voltage difference under steady-state vacuum operation conditions.

    According to the Clump theory under non-uniform electric field, the local electric field E on the grid surface is an important factor affecting the breakdown voltage, as particles accelerate and bombard the opposite grid with energy QV proportional to the product of EV, where Q is the particle charge, E is the local electric field strength and V is the applied voltage [15]. Therefore, optimizing the geometric shape between grids to improve voltage retention capability is one of the solutions to alleviate strong field distribution. From figure 5, it can be seen that the electric field in the acceleration gap between PG and GG is not uniform, but concentrated in a certain position. The electric fields in positions 1 and 3 are too large, with a maximum value of 15.1 kV/mm, which may cause inter-electrode firing and breakdown in long-term operation. So the structural optimization is needed. However, in beam optical analysis, parameters such as grid gap and aperture have already been determined, so the strong field distribution can only be mitigated by changing the chamfer shape. Positions 1 and 3 between PG and GG are both controlled by the same chamfer. The chamfer value is continuously increased in step of 0.1 mm, and the maximum electric field simulation values of positions 1 and 3 are recorded respectively. As shown in table 1, as the chamfer values increase, the maximum electric fields of positions 1 and 3 decrease. However, when the chamfer value increases to 0.8 mm, the maximum electric field of position 1 increases instead, which is due to the superposition of the chamfers of the neighboring circular holes in the PG plate, resulting in a decrease of the PG-GG gap, so that the electric field begins to increase. It is finally determined that the chamfer values of the positions 1 and 3 between PG and GG are designed to be 0.7 mm.

    Figure  5.  Electric field distribution of PG-GG.
    Table  1.  The relationship between chamfer and electric field for positions 1 and 3.
    Chamfer (mm) Position 1 (kV/mm) Position 3 (kV/mm)
    0.415.114.2
    0.614.413.1
    0.714.212.7
    0.815.612.3
     | Show Table
    DownLoad: CSV

    As shown in figure 6, the voltage difference between the grids of GG-SG is the largest, with the electric fields at positions 1, 3 and 5 being too large, reaching a maximum value of 18.1 kV/mm. Structural optimization is needed for these three regions. Due to the limitation of the hole spacing of the GG plate, the chamfer value cannot be infinitely increased. Therefore, it is necessary to analyze these regions with larger electric fields within a reasonable range. The positions 1 and 3 between GG and SG are controlled by the same chamfer, and the chamfer value is continuously increased in 0.2 mm step to the maximum chamfer value of 1.6 mm with structural limitation, and the chamfer value of position 5 is continuously increased in 1.0 mm step. The maximum electric field simulation value of each point is shown in table 2. With the continuous increase of chamfer values, the maximum electric fields of positions 1, 3 and 5 are constantly decreasing. However, the chamfer value of position 5 cannot be too large, as it may cause the inner cooling water channel wall to be too thin, affecting the transportation of cooling water. The final determination is to design the chamfer values for positions 1 and 3 between GG and SG to be 1.6 mm, and for position 5 to be 5.0 mm.

    Figure  6.  Electric field distribution of GG-SG.
    Table  2.  The relationship between chamfer and electric field for positions 1, 3 and 5.
    Chamfer (mm) Position 1 (kV/mm) Position 3 (kV/mm) Chamfer (mm) Position 5 (kV/mm)
    1.0 14.4 13.5 2.0 18.1
    1.2 14.2 13.2 2.4 17.1
    1.4 14.1 13.0 3.0 16.0
    1.6 13.9 12.9 4.0 14.7
    - - - 5.0 13.7
    - - - 5.5 13.5
    - - - 6.0 12.9
     | Show Table
    DownLoad: CSV

    On the premise of ensuring the maximum beam density and the minimum beam divergence angle in the beam optical analysis [16], and according to the design reference of the 80 keV positive ion source accelerator, the maximum electric fields in the acceleration gaps between grids need to be less than 14 kV/mm. In summary, it is determined that the chamfer values for positions 1 and 3 between PG and GG are designed to be 0.7 mm, the chamfer values for positions 1 and 3 between GG and SG are designed to be 1.6 mm, and the chamfer value for position 5 is designed to be 5.0 mm. Therefore, the maximum electric fields between the four layers of grids are around 14 kV/mm, or less than 14 kV/mm, which can basically meet the 120 keV withstand voltage requirement. The finalized design values of all the chamfer values between the four layers of grids are shown in figure 7 and table 3. Table 3 includes the preliminary chamfer values determined by beam optical analysis and the final values determined by electrostatic simulation.

    Figure  7.  Schematic diagram of four grid holes.
    Table  3.  Chamfer data of four grid holes.
    Position Preliminary determined value (mm) Final value (mm) Position Preliminary determined
    value (mm)
    Final value (mm)
    R1 0.5 0.5 R7 1.0 1.6
    R2 0.5 0.5 R8 1.5 1.5
    R3 0.4 0.7 R9 0.5 0.5
    R4 1.0 1.0 R10 1.0 1.0
    R5 0.2 0.2 R11 1.0 1.0
    R6 1.4 1.4
     | Show Table
    DownLoad: CSV

    In order to support the physical research on the EAST, a new positive ion source with designed beam energy of 120 keV was proposed to be developed. A key issue in the development of the EAST ion source is the connection between insulator and support flange. The 120 keV positive ion source uses a four-grids circular hole accelerator, and PEEK insulators are fixed with 304 stainless steel flanges through bolts and sealed with sealing rings, which plays the role of supporting the entire beam source system and voltage insulation.

    The concentration of electric field is an important factor affecting the insulation performance of accelerator. In order to ensure the stability of the beam, a three-dimensional 120 keV positive ion source accelerator model was established, and the withstand voltage performance of insulations and grids was analyzed through electrostatic simulations. The electrostatic simulations of insulators show that: (1) the height of the insulator is negatively correlated with the maximum electric field of its surface in the region; (2) among the three materials, PEEK material has the best insulating property, so PEEK material is selected for insulators of the 120 keV accelerator; (3) when the heights of insulators are the design values, the maximum electric fields on the surface of the PG-GG insulator and GG-SG insulator are about 4 kV/mm, which is able to meet the pressure resistance requirement of 120 keV. In addition, the electric fields in the acceleration gaps between grids are not uniform, so grids need to be optimized in structure. Electrostatic simulations of grids show that: (1) the chamfer values between PG and GG in positions 1 and 3 are designed to be 0.7 mm; (2) the chamfer values for positions 1 and 3 between GG and SG are designed to be 1.6 mm, and the chamfer value for position 5 is designed to be 5.0 mm; (3) after optimizing the structures between grids, the maximum electric fields between the four layers of grids are around 14 kV/mm, or less than 14 kV/mm, which can meet the 120 keV withstand voltage requirement.

    This work was supported by National Natural Science Foundation of China (No. 12375227) and Innovation in Fusion Engineering Technology of Institute (No. E35QT1080C).

  • [1]
    Loarte A et al 2007 Nucl. Fusion 47 S203 doi: 10.1088/0029-5515/47/6/S04
    [2]
    Pitts R A et al 2019 Nucl. Mater. Energy 20 100696 doi: 10.1016/j.nme.2019.100696
    [3]
    Ye M Y et al 2019 Nucl. Fusion 59 096049 doi: 10.1088/1741-4326/ab2bd0
    [4]
    Zhang M Z et al 2024 Contrib. Plasma Phys. 64 e202300135 doi: 10.1002/ctpp.202300135
    [5]
    Kallenbach A et al 2010 Plasma Phys. Control. Fusion 52 055002 doi: 10.1088/0741-3335/52/5/055002
    [6]
    Maddison G P et al 2003 Nucl. Fusion 43 49 doi: 10.1088/0029-5515/43/1/306
    [7]
    Eldon D et al 2021 Nucl. Mater. Energy 27 100963. doi: 10.1016/j.nme.2021.100963
    [8]
    Ravensbergen T et al 2021 Nat. Commun. 12 1105 doi: 10.1038/s41467-021-21268-3
    [9]
    Koenders J T W et al 2023 Nucl. Fusion 63 026006 doi: 10.1088/1741-4326/aca620
    [10]
    Bernert M et al 2021 Nucl. Fusion 61 024001 doi: 10.1088/1741-4326/abc936
    [11]
    Wang L et al 2022 Nucl. Fusion 62 076002 doi: 10.1088/1741-4326/ac4774
    [12]
    Wu K et al 2023 Nucl. Mater. Energy 34 101398 doi: 10.1016/j.nme.2023.101398
    [13]
    Chen M W et al 2020 Nucl. Fusion 60 076009 doi: 10.1088/1741-4326/ab8c65
    [14]
    Wu K et al 2018 Nucl. Fusion 58 056019 doi: 10.1088/1741-4326/aab506
    [15]
    Tao Y Q et al 2023 Nucl. Fusion 63 076008 doi: 10.1088/1741-4326/acd014
    [16]
    Lengyel L L 1981 Analysis of radiating plasma boundary layers IPP Report 1/191 Garching: Max-Planck-Institut für Plasmaphysik
    [17]
    Kallenbach A et al 2016 Plasma Phys. Control. Fusion 58 045013 doi: 10.1088/0741-3335/58/4/045013
    [18]
    Goldston R J, Reinke M L and Schwartz J A 2017 Plasma Phys. Control. Fusion 59 055015 doi: 10.1088/1361-6587/aa5e6e
    [19]
    Reinke M L 2017 Nucl. Fusion 57 034004 doi: 10.1088/1741-4326/aa5145
    [20]
    Huber A and Chankin A V 2021 Nucl. Fusion 61 036049 doi: 10.1088/1741-4326/abdded
    [21]
    Stangeby P C 2018 Plasma Phys. Control. Fusion 60 044022 doi: 10.1088/1361-6587/aaacf6
    [22]
    Zhang X et al 2023 Nucl. Mater. Energy 34 101354 doi: 10.1016/j.nme.2022.101354
    [23]
    Krasheninnikov S I, Kukushkin A S and Pshenov A A 2022 Plasma Phys. Control. Fusion 64 125011 doi: 10.1088/1361-6587/ac9a6f
    [24]
    Krasheninnikov S I and Kukushkin A S 2017 J. Plasma Phys. 83 155830501 doi: 10.1017/S0022377817000654
    [25]
    Hiwatari R et al 2005 J. Nucl. Mater. 337–339 386 doi: 10.1016/j.jnucmat.2004.10.043
    [26]
    Ou J, Gan C Y and Ye L 2014 Plasma Sci. Technol. 16 907 doi: 10.1088/1009-0630/16/10/02
    [27]
    Chen J L, Jia G Z and Xiang N 2021 J. Fusion Energy 40 1 doi: 10.1007/s10894-021-00292-7
    [28]
    Wang M et al 2022 Plasma Sci. Technol. 24 015101 doi: 10.1088/2058-6272/ac320f
    [29]
    Wang F Q et al 2023 Plasma Sci. Technol. 25 115102 doi: 10.1088/2058-6272/ace026
    [30]
    Dragojlovic Z et al 2010 Fusion Eng. Des. 85 243 doi: 10.1016/j.fusengdes.2010.02.015
    [31]
    Eich T et al 2011 Phys. Rev. Lett. 107 215001 doi: 10.1103/PhysRevLett.107.215001
    [32]
    Stangeby P 2020 Plasma Phys. Control. Fusion 62 025013 doi: 10.1088/1361-6587/ab51d6
    [33]
    Goswami R et al 2021 Phys. Plasmas 8 857 doi: 10.1063/1.1342028
    [34]
    Mavrin A A 2017 J. Fusion Energy 36 161 doi: 10.1007/s10894-017-0136-z
    [35]
    Zhou Y L et al 2022 Plasma Phys. Control. Fusion 64 065006 doi: 10.1088/1361-6587/ac6827

Catalog

    Figures(9)

    Article views (21) PDF downloads (8) Cited by()

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return