
Citation: | Chenyao JIN, Chi-Shung YIP, Haigang SUN, Di JIANG, Wei ZHANG, Guosheng XU, Liang WANG. Presheath formation and area design limit satellite-based Langmuir probes[J]. Plasma Science and Technology, 2023, 25(12): 124001. DOI: 10.1088/2058-6272/ace676 |
In this article, the effect of the finite conductive surface area of a satellite on the use of satellite-based Langmuir probes is reviewed in light of the basic theory of asymmetric double Langmuir probes (ADLPs). Recent theoretical and experimental studies have discussed electron sheath/presheath formation and the electron Bohm criterion along with their implications for satellite-based Langmuir probes. The effects predicted by the latest theory of the electron Bohm criterion were not experimentally observed and the experimental results remain supportive of a critical area ratio (AL/AS)crit = (mi/(2.3me))1/2 between the probe area AS and the satellite area AL as conventionally believed. A satellite-based Langmuir probe must satisfy this criterion to physically act as a single Langmuir probe. However, experimental investigations also found that high-energy electrons adversely affect (AL/AS)crit and a Langmuir probe's signal quality by giving additional electron current to AL. Based on these results, a number of limitations of the maximum probe area are derived when designing satellite-based Langmuir probes, with consideration of both the aim of the satellite and the plasma where the satellite-based probe works. These proposed measures are expected to only partially alleviate the effect of the inadequate satellite surface area on the application of satellite-based Langmuir probes. Using a larger satellite to carry a Langmuir probe remains the most viable means to obtain precise space plasma parameters.
A Hall thruster is a type of power device used in the field of electric propulsion, and its main principle is to ionize and accelerate the propellant gas (Xenon in general) using electric and magnetic fields and to generate thrust [1, 2]. Hall thrusters have the advantage of being highly efficient, of having highly specific impulses and a simple structure. In addition, with the rapid development of the commercial aerospace industry, Hall thrusters will become more widely used in aerospace fields such as the propulsion of small or microsatellites [3, 4].
Hall thruster consists of a discharge channel, anode, magnetic circuit, and cathode. The magnetic field in the Hall thruster is formed by the magnetic circuit, mainly in the radial direction, and the electrons are magnetized, whereas the ions are not. Electrons are emitted from the cathode, they enter the channel, and are captured by the magnetic field. The neutral gas diffuses from the anode to the downstream of the channel. Ions are generated by the ionization collisions between atoms and electrons, and then accelerated by the self-consistent electric field to form thrust.
A multiple ionization effect exists in the working process of a Hall thruster, although the ionization of the propellant gas is dominated by a single ionization [5, 6]. Meanwhile, the multiple ionization process is closely related to the discharge condition of the thruster, such as discharge voltage, propellant flow rate, magnetic field configuration, and the amount of vacuum [6–10]. The geometry of the thruster also affects the characteristics of multiple ionization [11, 12]. The characteristics of multiple ionization affect erosion of the discharge chamber and efficiency, therefore it is worthy of further study [9, 11, 12].
The magnetically shielded Hall thruster has a significant advantage in extending the service life of a Hall thruster, therefore this technique is widely used [13–16]. Hofer et al used an E × B probe to measure the proportion of multiple charged ions on two 6 kW thrusters, one provided with magnetic shielding technology and the other without such technology, and found that the application of magnetic shielding significantly enhances the multiple ionization process [17]. The main differences between magnetically shielded thrusters and unshielded thrusters in design include, the position of the strong magnetic field region with a positive gradient, the curvature of the magnetic field line near the wall area, and the shape of the discharge channel [18, 19]. For unshielded thrusters, the magnetic field with a positive gradient is fully located inside the channel, which is a normal loading manner. For magnetically shielded thrusters, part of the magnetic field with a positive gradient is located outside the channel; specifically, the magnetic field peak is shifted from the exit plane to outside. Therefore, it can be called the aft-loading of the magnetic field [20]. The experiment confirmed that the ionization zone and acceleration zone move downstream in the magnetically shielded thrusters [18]. As a result, part or all of the acceleration zone is located outside the channel. The majority of the potential drop is located downstream from the exit of the channel, which in turn causes great changes in the electron energy conversion process [18]. Thus, it is guessed that the transit from a normally loaded magnetic field to an aft-loaded one [20] is the main reason for the enhancement of multiple ionizations in the magnetically shielded Hall thruster. The physical mechanism of it has yet to be clearly explained. Based on this, a numerical simulation is used to compare and analyze the mechanisms of the influence of the two different loading manners of magnetic field on the multiple ionization process, and thereby help to understand the characteristics of the multiple ionization process of a Hall thruster with aft-loaded magnetic field.
The remainder of this paper is arranged as follows. Section 2 introduces the numerical model and research methods applied. Section 3 describes and analyzes the simulation results. Section 4 provides some concluding remarks.
Particle-in-cell (PIC) simulation is a powerful method to study the physical processes of Hall thrusters [21–24]. Because a Hall thruster has an axisymmetric cylindrical structure, its discharge can be considered uniform along the azimuthal direction. Therefore, only the axial (z) and radial (r) directions are considered when constructing the calculation model. Our team has developed a mature PIC platform with a 2D3V framework, which has been used to solve many physical problems. A brief introduction of the PIC model is presented below. The details can be found in [25–30]. To investigate the topic of this study, the PIC platform is updated by considering the multiple ionization process.
In the previous model, only elastic collision, excitation collision and single ionization collisions between atoms and electrons were considered. In order to study the multiple ionization, a total of eight types of independent collisions between electrons and atoms or ions are considered, where the electron is the incident particle, including an elastic collision, excitation collision, and six types of ionization collisions, namely, Xe→Xe+, Xe→Xe2+, Xe→Xe3+, Xe+→Xe2+, Xe+→Xe3+, Xe2+→Xe3+. In addition, the movement of high-charged ions is also considered. The collision cross-section data of elastic collision, excitation collision and single ionization are the same as those used by Szabo [31] and shown in figure 1(a). The collision cross-section data of multiple ionization collisions are summarized in [32–38] and shown in figure 1(b). The calculation process of an electron collision is optimized using the Monte Carlo Collison (MCC) method and the NULL-collision method proposed by Vahedi [39]. After electrons collide, the collision form can be determined through a single estimate. In addition, a Bohm anomalous conduction is also considered in the model [40] and the Bohm collision frequency νB is used to equal this effect:
νB=CBeB/me | (1) |
where CB is the semi-empirical coefficient, e is the basic charge quantity, B is the strength of the magnetic field, and me is the electron mass. Coulomb and charge exchange collisions have yet to be considered.
The induced magnetic field is ignored in the model since its magnitude is much smaller than that of the applied one. Therefore, only the electrostatic model is considered, and the Poisson equation is adopted and numerically solved to obtain the electric field.
For particle dynamics, the motion of the atoms is processed according to the free molecular flow; the motion of electrons is controlled by an electric and magnetic field; the motion of ions is only influenced by the electric field.
Five typical boundary conditions are considered, as shown in figure 2. On the boundary of the anode, the neutral atoms are injected and the potential is set as the discharge voltage. Generally, the material of the discharge channel is BN ceramic. Therefore, the channel wall is regarded as a dielectric boundary. On the boundary of the dielectric wall, a normal electric field is determined by the charge density deposited on the wall. In addition, the secondary electron emission (SEE) effect of an electron-wall collision is considered in the model [41, 42]. The face of the pole boundary is considered a capacitance boundary, and its potential is determined by the amount of charge deposited on the surface. On the open boundary, electrical neutrality is maintained by injecting electrons, and the potential is set to zero. Moreover, on the symmetry boundary, the radial electric field is set to zero and the particles are specularly reflected.
In this study, a simulated thruster with a size similar to those of SPT-100 and PPS-1350 [43] with a power of 1.35 kW was selected, as shown in figure 2. The magnetic field used in the PIC simulation is calculated using the open-source software named Finite Element Method Magnetics (FEMM). According to the research purpose, two magnetic field configurations, aft-loaded magnetic field and normal loaded magnetic field were designed, as shown in figure 3. For the aft-loaded magnetic field, the position of the maximum magnetic field strength realized along the centerline is pushed from the channel exit to a site 4 mm outside the channel. The magnetic field strength at the channel exit on the channel centerline accounts for 86% of the maximum magnetic field, and this ratio in the 6 kW power magnetic-shielding Hall thruster developed by Jet Propulsion Laboratory (JPL) is 83% [13, 14, 19]. Thus, they are similar in magnetic field intensity distribution, which satisfies the characteristics of the aft-loaded magnetic field.
In the simulation, the discharge voltage is set as 300 V and the anode mass flow rate is set as 5.65 mg s-1. Regarding the Bohm coefficient CB in the semi-empirical formula, it was found that the CB value outside the channel is larger than that inside the channel [40, 44]. Considering that part of the positive-gradient region of the aft-loaded magnetic field is located outside the channel, different Bohm coefficients are used depending on the gradient of the magnetic field [45]. CB is selected as 1/310 for the positive-gradient region and 1/18 for the negative-gradient region. These Bohm coefficients are chosen to ensure that the simulated performance is closer to the experimental one [44, 46, 47].
In order to compare the proportions of differently charged ions in the discharge process, the ion current ratio and ion quantity ratio need to define respectively. It is shown in formulas (2)–(4), where Iia is ion current formed by Xea+, Iitotal is sum of Iia, Rcurrenta is ion current ratio for Xea+, and Rquantitya is ion quantity ratio for Xea+. a represents the charge state of the Xe ion.
Itotal i=3∑a=1Iai | (2) |
Racurrent =IaiItotal i×100%(a=1,2,3) | (3) |
Raquantity =Iai/a∑(Iai/a)×100%(a=1,2,3) | (4) |
Table 1 shows the different charged ion currents and their proportions of the total currents under different magnetic field loading manners. It can be seen that the change in loading manner of magnetic field from normal loading to aft-loading enhances the multiple ionization effect, where the proportion of Xe2+ and Xe3+ increases.
Magnetic field loading manner | Aft-loading | Normal loading | |
Xe+ | Ion current (A) | 3.01 | 3.31 |
Ion current ratio (%) | 72.91 | 80.10 | |
Ion quantity ratio (%) | 85.20 | 89.30 | |
Xe2+ | Ion current (A) | 0.9 | 0.74 |
Ion current ratio (%) | 21.83 | 17.80 | |
Ion quantity ratio (%) | 12.75 | 9.92 | |
Xe3+ | Ion current (A) | 0.22 | 0.09 |
Ion current ratio (%) | 5.26 | 2.10 | |
Ion quantity ratio (%) | 2.05 | 0.78 |
To further analyze the influence of the magnetic field loading manner on the ionization process, the contour of the ionization collision rate distribution of the particle collisions is given. It should be noted that, according to the introduction to the PIC numerical model in section 2.1, six types of ionization collision processes were considered in this study. Undoubtedly, the process of Xe→Xe+ plays a leading role, and the Xe→Xe+ ionization collision rate distributions are shown in figure 4. The proportion of Xe2+ in the product is second only to Xe+, therefore, the following content will mainly focus on the generation of Xe2+. The generation of Xe2+ comes from Xe+→Xe2+ and Xe→Xe2+, and their ionization collision rate distributions are shown in figures 5 and 6.
First, from the perspective of the ionization rate, the ionization rate of the Xe→Xe+ ionization process corresponding to the normal loaded magnetic field is higher, and the ionization rates of Xe→Xe2+ and Xe+→Xe2+ corresponding to the aft-loaded magnetic field are both higher. The ionization rate represents the probability of occurring of multiple ionization. This is consistent with the conclusion that the aft-loading of the magnetic field enhances the multiple ionization process. In addition, the results in figures 5 and 6 can illustrate that there are two ways to generate Xe2+ in both cases, namely Xe→Xe2+ and Xe+→Xe2+, and between the two, Xe+→Xe2+ is dominant. Similarly, by using the PIC for a calculation and comparison, it can be seen that Xe+→Xe3+ is dominant among the three ionization processes generated for Xe3+, namely, Xe→Xe3+, Xe+→Xe3+, and Xe2+→Xe3+. So it can be inferred multiply charged ions are more likely created from singly charged ions than directly from neutrals.
Second, two characteristics of the multiple ionization process can be concluded from the contour of the ionization rate distribution. On one hand, compared to the normal loading case, in the aft-loading case, multiple ionization collisions occur at the position closer to the channel downstream. On the other hand, for the generation of Xe2+, in both cases, compared to the ionization process of Xe→Xe2+, the ionization process of Xe+→Xe2+ occurs on closer to the channel downstream.
In a Hall thruster, the ionization rate Rion of the propellant is shown in formula (5), where ne is the electron density, na is the atom density, and βi is the ionization rate coefficient. Since ionization rate coefficient and electron temperature are directly related, the atom density and electron temperature are the main factors determining the magnitude of the ionization rate. The distributions of electric potential and electron temperature on the channel centerline are shown in figure 7.
Rion =nenaβi | (5) |
In addition, the magnetic field strength is one of the key factors that determine the spatial conductivity of electrons, so the magnetic field strength and its distribution are the key factors that determine spatial conductivity. Combining figures 3 and 7, compared to normal loading, aft-loading of magnetic field makes the positive-gradient magnetic field partly located outside the channel, causing significant changes in the conductivity inside and outside the channel. In the aft-loading, the conductivity inside the channel increases, and the potential drop in the channel decreases.
In the aft-loading case, the main potential drop is located outside the discharge channel, reaching 240 V on the channel centerline. In addition, the electrons emitted from the cathode rapidly reach a high electron temperature with a peak value of 54 eV under the strong heating of the electric field outside the channel. Such a high electron temperature enhances the ionization collision effect outside the channel, although compared to that inside the channel, the density of heavy particles outside the channel remains low, and thus the electrons are still at a high temperature of 32 eV when reaching the channel exit.
In the normal loading case, the main potential drop is located inside the channel, reaching 80 V on the channel centerline, and thus the heating effect of the electric field outside the channel on the electrons is poor, and the peak temperature of the electrons outside the channel is only 19 eV. As the electrons move toward the channel exit, the ionization energy loss is slightly higher than the heating effect of the electric field, and the temperature of the electrons at the channel exit is 18 eV.
Combining figures 4 and 5, when an electron emitted from the cathode moves toward the near-exit, it first encounters a large number of propellant atoms, and when propellant is freely diffused from the anode into the channel downstream, the gas becomes evenly dispersed within the entire channel between the inner and outer walls, and the gas density along the centerline experiences a significant drop from upstream to downstream. In aft-loading case, the temperature of the electron entering the channel is 32 eV. Such a high electron temperature causes an immediately strong ionization collision between the atom and electron. Considering figure 1, the probability of occurrence of a multiple ionization collision is high. In a normal loading case, the electrons entering the channel are at a low electron temperature of 18 eV. According to formula (4), such a propellant density and electron temperature environment fails to meet the condition of a large number of ionizations. Therefore, the electron continues to move upstream under the heating effect of the electric field within the channel, although because the wall surface energy loss of the electron in the channel is also large, the electron temperature is basically maintained at 20 eV, which is much lower than the electron temperature of 32 eV at the channel exit in aft-loading case. However, as the electron moves closer to the anode, the atom density of the propellant increases significantly. According to formula (4), there will be numerous ionization collisions between the electrons and atoms, causing an avalanche effect, and thus ionization is completed. Compared to the aft-loading case, the rate of double and even higher ionization collisions is low owing to the low electron temperature, and as the ionization area is closer to the anode, the propellant is distributed more intensively and the density becomes higher; thus, the rate of single ionization collisions increases.
Combining figures 3 and 6, the aft-loading of the magnetic field leads to the zone of ionization process of Xe→Xe+ to move downstream, and the generation of Xe2+ is mainly through process of Xe+→Xe2+. According to formula (4), the zone of ionization process of Xe+→Xe2+ must be downstream of the zone of ionization process of Xe→Xe+. Therefore, the zones of ionization process of Xe+→Xe2+ and Xe→Xe+ both move downstream. Compared to the normal loading case, in the aft-loading case, the main potential drop is outside the channel, where interaction between the channel and the electron is lost, the electron temperature is higher at the channel exit. The distance between the Xe→Xe+ ionization process zone and the exit of the channel is reduced, which leads to the weakening of the cooling effect of the channel and heavy particles on the electrons. Finally, electrons in the zone of ionization process of Xe+→Xe2+ have higher energy, therefore, the generation of Xe2+ is enhanced.
In section 3.1, the multiple ionization process was simulated under different magnetic field distributions but the same maximum magnetic field strength at the centerline of the channel. The change of the multiple ionization in different loading manners is due to the relative change of the spatial conductivity, and the conductivity is greatly affected by the strength of the magnetic field. Therefore, the influence of magnetic field strength on ionization process in aft-loading case needs to be considered. This part focuses on the influence of magnetic field strength on the ionization process in aft-loading case, because the strength and configuration of the magnetic field are both key factors for the ionization process. Based on previous experimental experience, the range of maximum magnetic field intensity of the channel centerline has been selected. The range of variation of the maximum magnetic field is 180–280 G under the conditions of the aft-loaded magnetic field. The simulation result is shown in figure 8. The proportion of Xe2+ and Xe3+ increases with the increase in magnetic field strength under the aft-loaded magnetic field.
This part still uses Xe2+ as the research object. It was pointed out in section 3.1 that the Xe+→Xe2+ ionization process is dominant in the ionization process with the generation of Xe2+, and thus the Xe+ density is a key factor determining the ionization rate of the Xe+→Xe2+ process.
In the positive gradient region of the magnetic field, the plasma instability is suppressed and the classical conduction is dominant; the electron conductivity is thus inversely proportional to the square of the magnetic field strength. The Bohm conduction is dominant in the negative gradient region of the magnetic field, where the electron conductivity is inversely proportional to the magnetic field strength. Therefore, in the aft-loading case, with the increase of the magnetic field strength, the ratio of the conductivity outside the channel to the inside of the channel increases, and the potential drop outside the channel decreases, as shown in figure 9(a). The value of electron temperature is the result of electric field heating and collision cooling. When the electron gains higher energy in the acceleration zone, the collision loss energy increases, and the electron temperature decreases in the ionization zone, as shown in figure 9(b). The production rate of Xe2+ is mainly determined by processes of Xe→Xe+ and Xe+→Xe2+, in which the electron temperature is the key factor. So higher electron temperature in the ionization zone is the main reason for the increase in production rate, as shown in figures 10(a) and (b). The results show that the proportion of Xe2+ increases as the magnetic field strength increases.
The variance of the proportion of Xe3+ with magnetic field strength has the same trend as that of Xe2+. It can be well-understood since both of them are generated mainly from the ionization of Xe+.
In this work, the mechanisms of the influence of two different loading manners of a magnetic field, normal loading and aft-loading, on the multiple ionization process of a Hall thruster were studied using a particle-in-cell simulation method. The simulation results proved the assumption that the shift of the strong magnetic field region with a positive gradient from the exit plane to outside indeed enhances the multiple ionization process. The increase of electron temperature in the ionization region is the main reason for the enhancement of multiple ionization. Compared with normal loading, in the aft-loading case, the position of the ionization zone moves downstream of the channel, resulting in the weakening of the cooling effect of the channel and the increase of the electron temperature. The increase of the magnetic field strength changes the relative conductivity of the ionization zone, reduces the collision loss in the acceleration zone, and leads to the increase of the electron temperature.
The authors gratefully acknowledge the late Dr Noah Hershkowitz, Irving Langmuir Professor Emeritus of the University of Wisconsin-Madison and our long-time collaborator, who has given invaluable advice and suggestions to the studies associated with sheath/presheath and diagnostics associated physics. This work is supported by National Natural Science Foundation of China (Nos. 12275305 and 12205334), the Chinese Academy of Science Hundred Youth Talent Program, China Postdoctoral Science Foundation (No. 2022M713188) and the Director's Fund of Hefei Institutes of Physical Science, Chinese Academy of Sciences (No. YZJJ2022QN19).
[1] |
Xu J C et al 2018 IEEE Trans. Plasma Sci.
46 1331 doi: 10.1109/TPS.2018.2794533
|
[2] |
Xu J C et al 2016 Rev. Sci. Instrum.
87 083504 doi: 10.1063/1.4960181
|
[3] |
Stangeby P C 1982 J. Phys. D: Appl. Phys.
15 1007 doi: 10.1088/0022-3727/15/6/011
|
[4] |
Godyak V A, Piejak R B and Alexandrovich B M 2002 Plasma Sources Sci. Technol.
11 525 doi: 10.1088/0963-0252/11/4/320
|
[5] |
Park J 2008 Adv. Space Res.
41 650 doi: 10.1016/j.asr.2007.04.089
|
[6] |
Lebreton J P et al 2006 Planet. Space Sci.
54 472 doi: 10.1016/j.pss.2005.10.017
|
[7] |
Ranvier S et al 2017 IEEE Trans. Plasma Sci.
45 2007 doi: 10.1109/TPS.2017.2700211
|
[8] |
Lee J C et al 2013 Curr. Appl Phys.
13 846 doi: 10.1016/j.cap.2012.12.011
|
[9] |
Lobbia R B and Beal B E 2017 J. Propul. Power
33 566 doi: 10.2514/1.B35531
|
[10] |
Liu C et al 2016 Acta Phys. Sin.
65 335(in Chinese)
|
[11] |
Liu C et al 2013 Chin. J. Space Sci.
33 331(in Chinese)
|
[12] |
Baalrud S D et al 2020 Plasma Sources Sci. Technol.
29 053001 doi: 10.1088/1361-6595/ab8177
|
[13] |
Baalrud S D, Hershkowitz N and Longmier B 2007 Phys. Plasmas
14 042109 doi: 10.1063/1.2722262
|
[14] |
Amemiya H and Fuchs G 1991 Jpn. J. Appl. Phys.
30 3531 doi: 10.1143/JJAP.30.3531
|
[15] |
Fish C et al 2012 Dice mission design, development, and implementation: success and challenges Proc. Small Satellite Conf. Logan 2012 https://digitalcommons.usu.edu/smallsat/2012/all2012/81/
|
[16] |
Hoang H et al 2019 Space Sci. Rev.
215 21 doi: 10.1007/s11214-019-0586-x
|
[17] |
Waets A, Cipriani F and Ranvier S 2019 IEEE Trans. Plasma Sci.
47 3689 doi: 10.1109/TPS.2019.2920136
|
[18] |
Langmuir I 1925 Phys. Rev.
26 585 doi: 10.1103/PhysRev.26.585
|
[19] |
Hershkowitz N 1989 How langmuir probes work ed O Auciello and D L Flamm Plasma Diagnostics(Boston: Academic) 113
|
[20] |
Sheridan T E 2000 Phys. Plasmas
7 3084 doi: 10.1063/1.874162
|
[21] |
Yee B T et al 2017 Plasma Sources Sci. Technol.
26 025009 doi: 10.1088/1361-6595/aa56d7
|
[22] |
Scheiner B 2016 et al Theory of sheaths near positively biased electrodes 2016 IEEE Int. Conf. on Plasma Science (ICOPS)(https://doi.org/10.1109/PLASMA.2016.7534341)
|
[23] |
Jin C Y et al 2022 Plasma Sources Sci. Technol.
31 115007 doi: 10.1088/1361-6595/ac9b8d
|
[24] |
Hershkowitz N 2005 Phys. Plasmas
12 055502 doi: 10.1063/1.1887189
|
[25] |
Scheiner B et al 2015 Phys. Plasmas
22 123520 doi: 10.1063/1.4939024
|
[26] |
Riemann K U 1991 J. Phys. D: Appl. Phys.
24 493 doi: 10.1088/0022-3727/24/4/001
|
[27] |
Chen F F 2006 Plasma Sources Sci. Technol.
15 773 doi: 10.1088/0963-0252/15/4/022
|
[28] |
Shimoyama M et al 2012 J. Phys. D: Appl. Phys.
45 075205 doi: 10.1088/0022-3727/45/7/075205
|
[29] |
Fang H K et al 2018 AIP Adv.
8 105220 doi: 10.1063/1.5052489
|
[30] |
Tejumola T R et al 2018 Acta Astronaut.
144 388 doi: 10.1016/j.actaastro.2018.01.016
|
[31] |
Suryana R et al 2021 J. Phys. Conf. Ser.
1876 012025 doi: 10.1088/1742-6596/1876/1/012025
|
[32] |
Bekkeng T A et al 2013 Ann. Geophys.
31 187 doi: 10.5194/angeo-31-187-2013
|
[33] |
Leon O et al 2022 IEEE Trans. Plasma Sci.
50 349 doi: 10.1109/TPS.2021.3137765
|
[34] |
Sun G Y et al 2022 Plasma Sci. Technol.
24 095401 doi: 10.1088/2058-6272/ac6aa7
|
[35] |
Yip C S and Hershkowitz N 2015 J. Phys. D: Appl. Phys.
48 395201 doi: 10.1088/0022-3727/48/39/395201
|
[36] |
Yip C S et al 2022 Plasma Sources Sci. Technol.
31 045002 doi: 10.1088/1361-6595/ac5c61
|
[37] |
Lim Y et al 2022 Plasma Sources Sci. Technol.
31 024001 doi: 10.1088/1361-6595/ac4d03
|
[38] |
Yip C S et al 2020 Plasma Sci. Technol.
22 085404 doi: 10.1088/2058-6272/ab7f3d
|
[39] |
Yip C S and Hershkowitz N 2015 Plasma Sources Sci. Technol.
24 034004 doi: 10.1088/0963-0252/24/3/034004
|
[40] |
Yip C S et al 2013 Plasma Sources Sci. Technol.
22 065002 doi: 10.1088/0963-0252/22/6/065002
|
[41] |
Godyak V A and Demidov V I 2011 J. Phys. D: Appl. Phys.
44 233001 doi: 10.1088/0022-3727/44/23/233001
|
[42] |
Liu Y Y et al 2019 J. Astronaut.
40 1115
|
[43] |
Chen F F 2001 Phys. Plasmas
8 3029 doi: 10.1063/1.1368874
|
[44] |
Diego P et al 2019 Ann. Geophys. Discuss. (https://doi.org/10.5194/angeo-2019-136)
|
[45] |
Golubkov G V et al 2019 Russ. J. Phys. Chem. B
13 874 doi: 10.1134/S1990793119050166
|
[46] |
Suvorova A V et al 2016 J. Geophys. Res. Space Phys.
121 5880 doi: 10.1002/2016JA022622
|
[47] |
Li P X et al 2020 Plasma Sources Sci. Technol.
29 025015 doi: 10.1088/1361-6595/ab69e5
|
[48] |
Oksuz L and Hershkowitz N 2004 Plasma Sources Sci. Technol.
13 263 doi: 10.1088/0963-0252/13/2/010
|
[49] |
Godyak V 2021 J. Appl. Phys.
129 041101 doi: 10.1063/5.0024258
|
[50] |
Sudit I D and Chen F F 1994 Plasma Sources Sci. Technol.
3 162 doi: 10.1088/0963-0252/3/2/006
|
[1] | Zhihui ZOU, Ping ZHU, Charlson C KIM, Wei DENG, Xianqu WANG, Yawei HOU. Frequency multiplication with toroidal mode number of kink/fishbone modes on a static HL-2A-like tokamak[J]. Plasma Science and Technology, 2022, 24(12): 124005. DOI: 10.1088/2058-6272/aca00a |
[2] | Yinxia GUAN (关银霞), Shiqiang WANG (王世强), Shiya TANG (唐诗雅), Hongxiang MU (牟洪祥), Quanzhen LIU (刘全桢). Effect of power parameters on micro-discharge induced by corona discharge[J]. Plasma Science and Technology, 2019, 21(5): 54007-054007. DOI: 10.1088/2058-6272/ab0336 |
[3] | Jiaqi DONG (董家齐). Kinetic micro-instabilities in the presence of impurities in toroidal magnetized plasmas[J]. Plasma Science and Technology, 2018, 20(9): 94005-094005. DOI: 10.1088/2058-6272/aad4f4 |
[4] | Shoujie HE (何寿杰), Peng WANG (王鹏), Jing HA (哈静), Baoming ZHANG (张宝铭), Zhao ZHANG (张钊), Qing LI (李庆). Effects of discharge parameters on the micro-hollow cathode sustained glow discharge[J]. Plasma Science and Technology, 2018, 20(5): 54006-054006. DOI: 10.1088/2058-6272/aab54b |
[5] | ZHENG Ting (郑婷), WU Bin (吴斌), XU Liqing (徐立清), HU Chundong (胡纯栋), ZANG Qing (臧庆), DING Siye (丁斯晔), LI Yingying (李颖颖), WU Xingquan (伍兴权), WANG Jinfang (王进芳), SHEN Biao (沈飙), ZHONG Guoqiang (钟国强), LI Hao (李昊), SHI Tonghui (石同辉), EAST Team. Fishbone Mode Excited by Deeply Trapped Energetic Beam Ions in EAST[J]. Plasma Science and Technology, 2016, 18(6): 595-600. DOI: 10.1088/1009-0630/18/6/03 |
[6] | LIU Chao(刘超), LIU Yue(刘悦), MA Zhaoshuai(马照帅). Effect of Equilibrium Current Profiles on External Kink Modes in Tokamaks[J]. Plasma Science and Technology, 2014, 16(8): 726-731. DOI: 10.1088/1009-0630/16/8/02 |
[7] | GAO Xiang (高翔), ZHANG Tao (张涛), HAN Xiang (韩翔), ZHANG Shoubiao (张寿彪), et al.. Observation of Pedestal Plasma Turbulence on EAST Tokamak[J]. Plasma Science and Technology, 2013, 15(8): 732-737. DOI: 10.1088/1009-0630/15/8/03 |
[8] | WANG Zhongtian (王中天), WANG Long (王龙), LONG Yongxing (龙永兴), DONG Jiaqi (董家齐), HE Zhixiong (何志雄), LIU Yu (刘宇), TANG Changjian (唐昌建). Shaping Effects of the E-Fishbone in Tokamaks[J]. Plasma Science and Technology, 2013, 15(1): 12-16. DOI: 10.1088/1009-0630/15/1/03 |
[9] | HE Feng (何锋), HE Shoujie (何寿杰), ZHAO Xiaofei (赵晓菲), GUO Bingang (郭滨刚), OUYANG Jiting (欧阳吉庭). Study of the Discharge Mode in Micro-Hollow Cathode[J]. Plasma Science and Technology, 2012, 14(12): 1079-1083. DOI: 10.1088/1009-0630/14/12/08 |
[10] | LI Jiquan, Y. KISHIMOTO. Wave-Number Spectral Characteristics of Drift Wave Micro-Turbulence with Large-Scale Structures[J]. Plasma Science and Technology, 2011, 13(3): 297-301. |
1. | Chen, Z., Huang, Z., Jiang, M. et al. J-TEXT achievements in turbulence and transport in support of future device/reactor. Plasma Science and Technology, 2024, 26(11): 114001. DOI:10.1088/2058-6272/ad663b |
2. | Tao, Y., Li, M.H., Qian, J.P. et al. Simulation of current density profile for EAST advanced scenario with METIS code. Fusion Engineering and Design, 2024. DOI:10.1016/j.fusengdes.2024.114177 |
3. | Liu, Z., Fu, G. A simple model for internal transport barrier induced by fishbone in tokamak plasmas. Journal of Plasma Physics, 2023, 89(6): 905890612. DOI:10.1017/S0022377823001344 |
4. | Zhang, W., Wu, T.Y., Li, Y.G. et al. Compact phase comparison system for the synthetic HCOOH laser diagnostic system in HL-2A. Journal of Instrumentation, 2022, 17(9): P09037. DOI:10.1088/1748-0221/17/09/P09037 |
5. | Chao, Y., Zhang, W., Hu, L. et al. Observations of mode frequency increase and the appearance of ITB during the m/n = 1/1 kink mode in EAST high electron temperature long pulse operation. Plasma Science and Technology, 2022, 25(2): 025107. DOI:10.1088/2058-6272/ac92d0 |
Magnetic field loading manner | Aft-loading | Normal loading | |
Xe+ | Ion current (A) | 3.01 | 3.31 |
Ion current ratio (%) | 72.91 | 80.10 | |
Ion quantity ratio (%) | 85.20 | 89.30 | |
Xe2+ | Ion current (A) | 0.9 | 0.74 |
Ion current ratio (%) | 21.83 | 17.80 | |
Ion quantity ratio (%) | 12.75 | 9.92 | |
Xe3+ | Ion current (A) | 0.22 | 0.09 |
Ion current ratio (%) | 5.26 | 2.10 | |
Ion quantity ratio (%) | 2.05 | 0.78 |