
Citation: | Shanwei HOU, Tingfeng MING, Qiqi SHI, Wenmin ZHANG, Feifei LONG, Bang LI, Rongjing DENG, Qing ZHUANG, Xiaoju LIU, Ling ZHANG, Ang TI, Long ZENG, Guoqiang LI, Naoko ASHIKAWA, Satoshi OHDACHI, Xiang GAO. Experimental observations of naturally occurring dust using a high-speed vacuum ultraviolet imaging system in EAST[J]. Plasma Science and Technology, 2023, 25(5): 055101. DOI: 10.1088/2058-6272/acaaef |
In the ELMy H-mode experiment, naturally occurring dust originating at the high-field side is clearly observed using the high-speed vacuum ultraviolet imaging system developed on the Experimental Advanced Superconducting Tokamak (EAST). The main ablation cloud shape is similar to the classical shape observed in pellet fueling experiments. However, during the dust penetration, an erupted secondary cloudlet with a bent 'cigar' shape is observed and moves upwards along the direction perpendicular to the magnetic field line, which is different to the obviation in the pellet fueling experiments. This may be due to the ion diamagnetic drift effect. The velocities of the secondary cloudlet are estimated to be 50‒80 m s−1. In addition, a significant degradation of the plasma confinement is observed during the dust penetration.
The initiation of tokamak plasmas under different heating scenarios has been explored through radio-frequency heating [1, 2] and neutral beam injection [3, 4]. As a fundamental technique, Ohmic heating (OH) has emerged as an effective method to start up the plasma currents in tokamaks [5, 6]. Since the inception of tokamak research, central solenoid (CS) coils featuring an iron core have been employed to apply high toroidal-loop voltages. To save on the central space and improve the cost-effectiveness, modern tokamak designs tend to incorporate air-core coils as the CS. If the CS coils are manufactured using superconductive materials, the application of low loop voltages to the toroidal plasma is rendered feasible.
In particular, the CS space required for a spherical tokamak is less than that for conventional tokamaks; hence, air-core coils with reduced volt-seconds are installed as the CS. To overcome this limitation, which focuses on spherical tokamak, current start-up methods to reduce the dependence on OH have been established. From this perspective, co-axial helicity injection (CHI) [7–9], local helicity injection (LHI) [10], and plasma merging [11] have accomplished effective results in achieving plasma current start-up. As these methods consume high levels of power transiently, they result in the generation of high-density plasmas. In particular, radio-frequency heating is advantageous for maintaining the tokamak equilibrium across a longer time scale than the transient methods. Experiments with lower hybrid resonance heating and electron cyclotronheating (ECH) have verified that energetic electrons play an essential role in initiating non-inductive plasma [12].
A combination of OH and other heating methods would be the fundamental scenario to generate plasmas stably in advanced spherical tokamaks [13, 14]. ECH is established as the main heating source in the QUEST spherical tokamak [15–17]. In a previous experiment using a single flux swing, the plasma current reached 100 kA [12]. With a higher plasma current, a higher electron density could be achieved efficiently. Under the condition that the electron density exceeds the cut-off density for the EC wave, the electron Bernstein wave (EBW) transformed from the EC waves can be excited, thereby heating of the over-dense plasmas [18]. Notably, a steady-state heating and current-drive experiment is planned for the QUEST project using an 8.56 GHz EC-heating system with 250 kW RF power [19]. Additionally, this heating system is used for EBW experiments. Although energetic electrons are generated through RF start-up, the EBW heats bulk electrons. In this regard, a high-density discharge through CHI or LHI may help prevent the initiation of an excessive number of energetic electrons.
A target plasma is required to conduct an experiment on the heating and current-drive via EBW in QUEST. If the electron density exceeds
The remainder of this paper is organized as follows. The design of the new OH system is introduced in section 2. This section presents a conceptual circuit, a current start-up scenario, and a control system of the OH system. An interlock system based on a field programmable gate array (FPGA) is designed to protect the upgraded OH system. Section 3 illustrates the experimental results corresponding to the non-inductive and inductive current start-ups. Finally, the conclusions are drawn in section 4.
The power supply system is established to flexibly energize the CS coil. The CS coil is located in the center of the QUEST consisting of three coils, named PF4-1, PF4-2, and PF4-3, in series. PF4-1 and PF4-3 each has 72 turns and PF4-2 has 144 turns. In OH experiment, the total number of turns is 216 with PF4-1, PF4-3, and half of PF4-2 connected in series PF4-1 and PF4-3, suited on the top and bottom. Each has a radius of 163 mm and length of 590 mm, while PF4-2 in the middle, has a radius of 172 m and length of 792 mm. The CS coil can be considered as a combination of an inductance, LCS=1.88 mH, and a resistance,
The two DC power supplies, the positive supply (PS) and negative supply (NS), are powered by three-phase ACs through the control of six thyristor stacks. Using both the PS and NS power supplies, bipolar voltages can be applied to the CS coil, provided as they are seamlessly controlled. Additionally, the two IGBT stacks installed between the power supplies and the CS coil function as switches.
The bipolar power supply is controlled by the QUEST Central Control System (QCCS), while the IGBT stacks are controlled by the IGBT control system, which mainly comprises the FPGA, isolation amplifiers, transmitters, and receivers. The IGBT gate signals, transmitted by the FPGA, are used to switch on/off the IGBTs. These signals are transferred by transmitters and receivers through light cables to the IGBT drivers. Accordingly, an equivalent circuit is established, as illustrated in figure 1.
Lp=μ0R(ln8Ra+li2-2), | (1) |
Rp=(1ene·veie)2πRπa2=η2πRπa2, | (2) |
where the major radius is denoted by R, the minor radius is symbolized by a, and internal inductance values
η=(meene·veie)=1.65×10-9×Zeff×lnΛTe32, | (3) |
where Coulomb logarithm
As shown by the arrows depicted in figure 1, the operation controlling PS and NS in OH is conducted in two periods. In the first period, the PS and IGBT stack 1 are turned on, and the current in the CS loop flows clockwise, as indicated by the red arrow. Conversely, in the second period, the PS and IGBT stack 1 are turned off, and the NS and IGBT stack 2 are subsequently turned on. In this case, the CS current flows through IGBT 2 and NS counterclockwise, as shown by the blue arrows.
The equivalent circuit is simulated using a power electronics simulator (PSIM; Powersim Inc.), as indicated in figure 2. The plasma resistance decreases with increasing temperature. Considering both the thermal and energetic electrons,
The PS and NS supply negative and positive voltages in this circuit setup. On applying -375 V to the CS, the current can increase to 8 kA. A loop voltage through the first swing is supplied as the current rapidly decreases to 0 A. Furthermore, the plasma current starts increasing, following the working principle of the transformer; it increases up to 60 kA at 2.52 s. The NS provides a voltage of up to 375 V to the CS coil, and the current increases to 8 kA from 2.52 to 2.62 s. Finally, the plasma current attains a maximum value of approximately 120 kA through a double flux swing.
As the maximum CS current reaches 8 kA during the OH operation, this working condition entails a risk of damage to the power-supply if the IGBT stacks are not appropriately controlled. For example, the current via the PS may flow through IGBT stack 2 to the closed circuit with NS if both IGBT stacks 1 and 2 are turned on simultaneously. Therefore, to protect the circuit components and power supplies, an interlock system was designed using an FPGA configuration; this system detects the IGBT signals and yields enabling signals for the power supplies and IGBT stacks.
This FPGA interlock system comprises two parts:
(1) Timer
The logical control of the PS timer, utilized to control and drive other FPGA interlock signals, is presented as a block diagram in figure 3. A timer signal is produced at the trigger timing by the QCCS.
When a trigger is received from the QCCS, a tick count of the timer records the present time (
(2) Generator of IGBT gate.
The generator of the IGBT gate serves to receive a signal from the CS circuit. This gate signal controls the on-and-off operation of the IGBT stack by analyzing these analog signals. The detailed logic of the IGBT gate is presented in figure 4.
Figure 4 illustrates a switch connected with an AND gate and an OR gate controlled through the IGBT gate. In the case of PS, if both the IGBTs are in the 'Off' state, the output signal of the switch is produced through Port 1. If the answerbacks from PS IGBT and NS IGBT represent the 'On' and 'Off' states, respectively, then the switch output signal is produced through Port 2.
The AND gate of Port 1 protects the circuit from the arbitrary switching-on of the IGBTs. The FPGA receives the signal value of the CS coil current (
Figure 5 plots the conceptual logic waveforms of the IGBT interlock system. The timer signal PS ends at
The interlock system is tested, and the corresponding waveforms are plotted in figure 6. The QCCS generates signals of a sequence when the PS and NS are in the working state. At t=2.36 s, the sequence signal of PS ends but
In this study, an experiment was conducted to verify the effect of the upgraded OH and its control system. The results from shot numbers (Nos.) 47762 and 48052 were selected as the representative discharge of a non-inductive current start-up and an OH start-up with double flux swing, respectively. The waveforms of the plasma current and CS current are presented in figure 7.
In this experiment, 28 GHz- ECH was used for ionization and current start-up in both discharges. In shot No. 47762, the plasma current is driven mainly via EC heating. Although the CS current is low, it is necessary for maintaining equilibrium. When
As the power supplied through the heating is insufficient to maintain such a high current, the plasma current decreases after attaining the peak value. A simulation result of the double flux swing, which has the comparable value of CS current with the experiment, is plotted with dotted green line. The maximum plasma current is about 70 kA that is lower compared with experimental result which is heated by ECH. The simulation result of the non-inductive discharge is plotted by dotted blue line, and it is lower than experimental result because the simulation is only concluded the OH part. In the experiment, non-inductive plasma start-up is mainly realized by ECH. At 2.8 s, EC heating is halted. Thereafter, the plasma current gradually decreases to 0 A because the energetic electrons are confined optimally.
The waveform of the loop voltage,
(4) |
where N (=1) represents the turn of a flux loop. In shot No. 47762, a non-inductive current start-up, the flux consumption only reaches -0.01 Wb at t=2.8 s. Specifically, owing to the double flux swing, the flux approaches -0.045 Wb through the PS and exceeds -0.1 Wb through the NS in shot No. 48052.
The relationship between the flux and the maximum plasma current, corresponding to three The relationship between the flux and maximum plasma current types of start-up methods, is plotted in figure 9. For the non-inductive start-up, the flux swing can merely attain a value of 0.03 Wb, and the maximum plasma current is approximately 60 kA. The single flux swing yields a more significant flux (0.07 Wb) and higher plasma current (100 kA) than those via the non-inductive approach. Notably, the double flux swing yields the maximal flux and plasma current in these three types. Ideally, neglecting the plasma resistance, the as-recorded loop voltage is solely affected by the plasma inductance,
(5) |
is plotted by the black dashed-line in figure 9. Due to the 28 GHz-ECH, which generates the initial plasma current in advance, the plasma current is greater than the ideal value.
The bulk electron temperature and density measured via Thomson scattering are displayed in figure 10. The horizontal axis represents the tangent radial position. At t=2.6 s, the electron density exceeds the cut-off density of
OH, through a double flux swing, was applied to the QUEST spherical tokamak. Furthermore, IGBT switching was used to combine the two high-power current sources, and a logical control system with an FPGA was installed in this bipolar power supply. In the experiments, plasma currents were produced through OH.
The authors gratefully acknowledge the EAST groups. This work is supported in part by the National Key R & D Program of China (Nos. 2019YFE03080200, 2022YFE03030001 and 2022YFE03050003), National Natural Science Foundation of China (Nos. 11975271, 12075284, 12075283 and 12175277), the Chinese Academy of President's International Fellowship Initiative (No. 2021VMA0022) and the Post-CUP program, JSPS-CAS Bilateral Joint Research Projects, 'Control of wall recycling on metallic plasma facing materials in fusion reactor', 2019–2022, (No. GJHZ201984).
[1] |
Krasheninnikov S I et al 2008 Plasma Phys. Control. Fusion 50 124054 doi: 10.1088/0741-3335/50/12/124054
|
[2] |
Smirnov R D et al 2015 Phys. Plasmas 22 012506 doi: 10.1063/1.4905704
|
[3] |
Shoji M et al 2015 Nucl. Fusion 55 053014 doi: 10.1088/0029-5515/55/5/053014
|
[4] |
Smirnov R D et al 2009 Plasma Phys. Control. Fusion 51 055017 doi: 10.1088/0741-3335/51/5/055017
|
[5] |
Brochard F et al 2017 Nucl. Fusion 57 036002 doi: 10.1088/1741-4326/aa4e56
|
[6] |
Flanagan J C et al 2014 Plasma Phys. Control. Fusion 57 014037 doi: 10.1088/0741-3335/57/1/014037
|
[7] |
Autricque A et al 2017 Nucl. Mater. Energy 12 599 doi: 10.1016/j.nme.2016.11.012
|
[8] |
Tang M et al 2011 J. Nucl. Mater. 415 S1094 doi: 10.1016/j.jnucmat.2011.01.001
|
[9] |
Yu J H et al 2009 J. Nucl. Mater. 390 216 doi: 10.1016/j.jnucmat.2009.01.053
|
[10] |
Saito K et al 2007 J. Nucl. Mater. 363–365 1323 doi: 10.1016/j.jnucmat.2007.01.269
|
[11] |
Huang Z H et al 2015 J. Nucl. Mater. 463 856 doi: 10.1016/j.jnucmat.2014.11.021
|
[12] |
Ming T F et al 2012 Rev. Sci. Instrum. 83 10E513 doi: 10.1063/1.4734494
|
[13] |
Zhou F et al 2017 Rev. Sci. Instrum. 88 073505 doi: 10.1063/1.4991856
|
[14] |
Gao X et al 2017 Nucl. Fusion 57 056021 doi: 10.1088/1741-4326/aa626c
|
[15] |
Li B et al 2021 Plasma Sci. Technol. 23 035104 doi: 10.1088/2058-6272/abda9d
|
[16] |
Craciunescu T et al 2010 J. Nucl. Mater. 400 205 doi: 10.1016/j.jnucmat.2010.03.005
|
[17] |
Wang Z H et al 2020 Phys. Plasmas 27 033703 doi: 10.1063/1.5134787
|
[18] |
Terry J L et al 1992 Rev. Sci. Instrum. 63 5191 doi: 10.1063/1.1143849
|
[19] |
Cseh G et al 2017 Nucl. Fusion 57 016022 doi: 10.1088/0029-5515/57/1/016022
|
[20] |
Chang C T 1991 Phys. Rep. 206 143 doi: 10.1016/0370-1573(91)90053-O
|
[21] |
Kuteev B V 1995 Nucl. Fusion 35 431 doi: 10.1088/0029-5515/35/4/I05
|
[22] |
Neuhauser J and Wunderlich R 1989 A Model for Periodic Pellet Ablation in Toroidal Plasma Confinement Experiments (Garching: Max-Planck-Institut für Plasmaphysik)
|
[23] |
Zhu X et al 2017 J. Instrum. 12 C12060 doi: 10.1088/1748-0221/12/12/C12060
|
[24] |
Xu Z et al 2016 Rev. Sci. Instrum. 87 11D429 doi: 10.1063/1.4961294
|
[25] |
Yang Y et al 2017 Plasma Phys. Control. Fusion 59 085003 doi: 10.1088/1361-6587/aa6f2a
|
[26] |
Duan Y M et al 2013 J. Nucl. Mater. 438 S338 doi: 10.1016/j.jnucmat.2013.01.063
|
[27] |
Zhang L et al 2015 Rev. Sci. Instrum. 86 123509 doi: 10.1063/1.4937723
|
[28] |
Pütterich T et al 2008 Plasma Phys. Control. Fusion 50 085016 doi: 10.1088/0741-3335/50/8/085016
|
[1] | Yuan YAO, Yao YANG, Ang TI, Yang SONG, Jiamin ZHANG, Yan WANG, Yao ZHANG, Haiqing LIU, Yinxian JIE. A fringe jump counting method for the phase measurement in the HCN laser interferometer on EAST and its FPGA-based implementation[J]. Plasma Science and Technology, 2024, 26(4): 045601. DOI: 10.1088/2058-6272/ad0dec |
[2] | Haojie MA, Huasheng XIE, Bo LI. Simulation of ion cyclotron wave heating in the EXL-50U spherical tokamak based on dispersion relations[J]. Plasma Science and Technology, 2024, 26(2): 025105. DOI: 10.1088/2058-6272/ad0d53 |
[3] | Tongyu WU (吴彤宇), Wei ZHANG (张伟), Haoxi WANG (王浩西), Yan ZHOU (周艳), Zejie YIN (阴泽杰). Research on the phase adjustment method for dispersion interferometer on HL-2A tokamak[J]. Plasma Science and Technology, 2018, 20(6): 65601-065601. DOI: 10.1088/2058-6272/aaaa19 |
[4] | Jun WU (吴军), Jian WU (吴健), Haisheng ZHAO (赵海生), Zhengwen XU (许正文). Analysis of incoherent scatter during ionospheric heating near the fifth electron gyrofrequency[J]. Plasma Science and Technology, 2017, 19(4): 45301-045301. DOI: 10.1088/2058-6272/aa58db |
[5] | CHANG Pengxiang (常鹏翔), WU Bin (吴斌), WANG Jinfang (王进芳), LI Yingying (李颖颖), WANG Xiaoguang (王小光), XU Handong (徐旵东), WANG Xiaojie (王晓洁), LIU Yong (刘永), ZHAO Hailin (赵海林), HAO Baolong (郝宝龙), YANG Zhen (杨振), ZHENG Ting (郑婷), HU Chundong (胡纯栋). The Influence of Neutral Beam Injection on the Heating and Current Drive with Electron Cyclotron Wave on EAST[J]. Plasma Science and Technology, 2016, 18(11): 1064-1068. DOI: 10.1088/1009-0630/18/11/02 |
[6] | GAO Min (高敏), CHEN Shaoyong (陈少永), TANG Changjian (唐昌建). Electron Cyclotron Harmonic Wave Heating in Tokamak Plasmas with Different Polarization Modes[J]. Plasma Science and Technology, 2013, 15(4): 313-317. DOI: 10.1088/1009-0630/15/4/02 |
[7] | PENG Jianfei (彭建飞), XUAN Weimin (宣伟民), WANG Haibing (王海兵), LI Huajun (李华俊), WANG Yingqiao (王英翘), WANG Shujin (王树锦). Study on Matching a 300 MVA Motor Generator with an Ohmic Heating Power Supply in HL-2M[J]. Plasma Science and Technology, 2013, 15(3): 300-302. DOI: 10.1088/1009-0630/15/3/22 |
[8] | DENG Yongfeng(邓永锋), TAN Chang(谭畅), HAN Xianwei(韩先伟), TAN Yonghua(谭永华). Numerical Simulation of the Self-Heating Effect Induced by Electron Beam Plasma in Atmosphere[J]. Plasma Science and Technology, 2012, 14(2): 89-93. DOI: 10.1088/1009-0630/14/2/01 |
[9] | YANG Yao, GAO Xiang, the EAST team. Energy Confinement of both Ohmic and LHW Plasma on EAST[J]. Plasma Science and Technology, 2011, 13(3): 312-315. |
[10] | ZHU Xueguang(朱学光). Influence of the Phase of the Antenna Current Standing Wave on the Power Flux in Ion Cyclotron Heating[J]. Plasma Science and Technology, 2010, 12(5): 543-546. |
1. | Zhang, Y., Onchi, T., Nakamura, K. et al. A versatile power supply system for the central solenoid of the QUEST spherical tokamak. Fusion Engineering and Design, 2023. DOI:10.1016/j.fusengdes.2023.113648 |
1. | Zhang, Y., Onchi, T., Nakamura, K. et al. A versatile power supply system for the central solenoid of the QUEST spherical tokamak. Fusion Engineering and Design, 2023. DOI:10.1016/j.fusengdes.2023.113648 |