Next Article in Journal
Photoionization of Astrophysically Relevant Atomic Ions at PIPE
Next Article in Special Issue
Plasma Spectroscopy on an Aluminum-Pellet Ablation Cloud in an LHD Plasma with an Echelle Spectrometer
Previous Article in Journal
The 13th International Colloquium on Atomic Spectra and Oscillator Strengths for Astrophysical and Laboratory Plasmas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spectroscopic Measurement of Hydrogen Atom Density in a Plasma Produced with 28 GHz ECH in QUEST

by Satoru Mori 1,*, Taiichi Shikama 1,*, Kazuaki Hanada 2, Nao Yoneda 1, Arseniy Kuzmin 1, Masahiro Hasuo 1, Hiroshi Idei 2, Takumi Onchi 2, Akira Ejiri 3, Yuki Osawa 3, Yi Peng 3, Kyohei Matsuzaki 3, Shinichiro Kado 4, Keiji Sawada 5, Takeshi Ido 2, Kazuo Nakamura 2, Ryuya Ikezoe 2, Yoshihiko Nagashima 2, Makoto Hasegawa 2, Kengo Kuroda 2, Aki Higashijima 2, Takahiro Nagata 2 and Shun Shimabukuro 2add Show full author list remove Hide full author list
1
Department of Mechanical Engineering and Science, Graduate School of Engineering, Kyoto University, Kyoto 615-8540, Japan
2
Research Institute for Applied Mechanics, Kyushu University, Fukuoka 816-8580, Japan
3
Graduate School of Frontier Sciences, The University of Tokyo, Chiba 277-8561, Japan
4
Institute of Advanced Energy, Kyoto University, Kyoto 610-0011, Japan
5
Graduate School of Engineering, Shinshu University, Nagano 380-8553, Japan
*
Authors to whom correspondence should be addressed.
Submission received: 23 June 2020 / Revised: 29 July 2020 / Accepted: 31 July 2020 / Published: 18 August 2020
(This article belongs to the Special Issue Atomic and Molecular Spectra in Magnetically Confined Torus Plasmas)

Abstract

:
The spatial distribution of the hydrogen atom density was evaluated in a spherical tokamak (ST) plasma sustained only with 28 GHz electron cyclotron heating (ECH). The radially resolved Hδ emissivity was measured using multiple viewing chord spectroscopy and Abel inversion. A collisional-radiative (CR) model analysis of the emissivity resulted in a ground-state hydrogen atom density of 1015–1016 m−3 and an ionization degree of 1–0.85 in the plasma.

1. Introduction

Q-shu University Experiment with Steady-State Spherical Tokamak (QUEST) is a medium-sized spherical tokamak (ST) device operated for the understanding of plasma–wall interactions in steady-state discharges [1] and for the development of non-inductive ST start-up schemes [2,3]. This study contributes to the latter by investigating the effect of neutrals (atoms) on an initial plasma produced with electron cyclotron heating (ECH).
In the ST discharge produced with ECH, a cylindrical plasma in an open field equilibrium is initially produced around the resonance layer [2,4,5]. Then, by ramping up the plasma current, an ST plasma sustained only with ECH is formed [2,4,6]. Among the characteristics of the initial ST plasma is a large variation in ion toroidal velocity in the range of 0–20 km/s with the size and position of the last closed flux surface (LCFS) [7]. Since this characteristic could be used for the active control of the flow profile, we intended to further investigate the driving and damping mechanisms of the flow. In this study, we measured the radial distribution of the atom density and the degree of ionization (DOI) in order to evaluate the loss and transport of the ion momentum through the charge-exchange reaction (e.g., [8]), since the bulk electron temperature and density of the initial plasmas are typically less than 1 keV and 1 × 1019 m−3 [2,6,9], respectively, and these relatively small plasma parameters can result in a greater effect of atoms.

2. Experimental Setup

We performed experiments using a discharge sustained only with 28 GHz ECH. This discharge is described in detail elsewhere [2]. Briefly, the plasma was mainly heated by second harmonic resonance ECH at R = 0.32 m with a power of 120 kW, where R is the major radius. We used two discharges performed with an identical setup, i.e., one for Thomson scattering (#41176) and the other for spectroscopy (#41178) measurements. The reproducibility of plasma parameters for these discharges is within their uncertainties. Figure 1a shows the temporal evolutions of the plasma current, electron temperature T e , and density n e at the core [9], and the chord-integrated Hα intensity measured on a radial viewing chord on the midplane. We conducted all the measurements in the time period t = 2.7–3.0 s denoted by the shaded area in Figure 1a. The period was in the flattop phase, where the magnetic field coil currents were kept constant. The reconstructed equilibrium magnetic flux surfaces are shown in Figure 1b; these surfaces were obtained from the measured magnetic flux and plasma current using a modification of the method described in [10]. The plasma minor radius was reduced to approximately 0.1 m to increase the volume power density of the heating. The magnetic axis was located at R = 0.32 m and z = −0.1 m, and the last closed flux surface (LCFS) was located at R = 0.25 and 0.4 m on the midplane.
Figure 2 is a schematic illustration of the visible spectroscopic system. We measured the Hδ emission line spectrum (p = 2–6), where p is the principal quantum number, using 24 fan-shaped viewing chords on the midplane. The distance between the viewing chord and the torus axis, which we define as η , ranged between 0.24 and 1.01 m. The diameters of the viewing spots were approximately 30 mm at the positions closest to the torus axis.
The emission was collected using two sets of a collimator (Edmund Optics Japan, Japan, TS achromatic lens; 15 mm focal length, ϕ 14 mm) and an optical fiber bundle (Mitsubishi Cable Industries, Japan, STU230D; 230 and 250 μm core and cladding diameters, respectively, 0.2 NA). The collected emission was transferred to a Czerny–Turner spectrometer (Acton Research, USA, AM-510; 1 m focal length, F/8.7, 1800 grooves/mm grating) equipped with astigmatism compensation optics. The spectra were recorded with a CCD (Andor Technology, UK, DU440-BU2; 2048 × 512 pixels, 13.5 μm pixel size, 16 bit). We set the entrance slit width of the spectrometer at 50 μm and the central output wavelength at 412 nm. The instrumental width at the full width at half maximum and the reciprocal linear dispersion were 40 pm and 6.5 pm/pixel, respectively. The sensitivity of the system was absolutely calibrated using a standard lamp while separately considering the window transmittance.

3. Results

Figure 3a shows a Hδ spectrum measured on the viewing chord at η = 0.51 m, which is indicated by the bold line in Figure 2; this position corresponds to the outer scrape-off-layer (SOL). We evaluated the chord-integrated intensity I A from the area of the spectrum. The evaluated I A is marked in Figure 3b. The error bars consist of uncertainties in the sensitivity calibration coefficient and the zero level of the spectrum. Note that the pedestal of the OII emission line at 410.32 nm is superposed on that of Hδ, but its effect on I A is small and we neglected it. Additionally, we omitted I A obtained on the chord at η = 0.84 m owing to the large noise produced by the incidence of X-rays to the CCD.
From the evaluated I A , we calculated the radially resolved emissivity ε ( R ) using Abel inversion assuming the toroidal symmetry of the emission. In a tokamak discharge with similar plasma parameters and wall conditions, the variation in Hα (p = 2–3) intensity in the toroidal direction was reported to be small except for an approximately 10 times local increase at the poloidal limiter [11]. Since the present discharge is limited by the center stack, the validity of the assumption of the toroidal symmetry could degrade near the center stack.
For the numerical calculation of the Abel inversion, the discrete data of I A was interpolated with a smooth spline function in the range η = 0.24–1.01 m as shown by the red curve in Figure 3b. ε ( R ) was then calculated as follows:
ε ( R )   =   1 π R R end 1 η 2 R 2 d I A ( η ) d η d η ,
where R end is the outermost radius of Hδ emission, which we determined to be 1.01 m because I A obtained on viewing chords at η   1.01 m was less than the noise level. The calculated ε ( R ) is shown in Figure 3c. We evaluated the error bar by the procedure described in [12] using 500 pairs of data. A small peak was found in ε ( R ) around R = 0.7 m owing to an increase in I A at η = 0.71 m. We confirmed the reproducibility of this tendency in identical discharges, but we have not identified its origin.
The hydrogen atom density was calculated from ε ( R ) by collisional-radiative (CR) model analysis [13]. The atom density in an excited state p , n H , p ( R ) , is expressed from the equation of continuity and the quasi steady-state approximation as follows:
n H , p ( R ) = r 0 p ( R ; T e , n e ) n H + ( R ) n e ( R ) + r 1 p ( R ; T e , n e ) n H ( R ) n e ( R )   + r 2 p ( R ; T e , n e ) n H 2 ( R ) n e ( R ) ,
where n H + , n H , and n H 2 are the ground-state hydrogen ion, atom, and molecule densities, respectively. r 0 p , r 1 p , and r 2 p are called population coefficients whose values are calculated with a CR model code [13] as a function of T e and n e .
Figure 4 shows the calculated population coefficients for p = 6, where the green, red, and blue lines represent r 06 , r 16 , and r 26 , and the solid, dashed, and dash-dotted lines are those at n e = 1016, 1017, and 1018 m−3, respectively. T e and n e of the present plasma are in the ranges of 10 T e 400 eV and 1016 n e 1018 m−3, respectively (see Figure 5a), and the former is denoted by the shaded area in Figure 4. Considering the relative magnitudes of the population coefficients and the approximated charge neutrality condition n H + n e , we can neglect the first term of Equation (2). We can also neglect the third term by assuming that n H 2 < n H , which is reasonable except for near the wall. Equation (2) is then approximated as follows:
n H , 6 ( R )     r 16 ( R ; T e , n e ) n H ( R ) n e ( R ) .
The excited hydrogen atom density n H , 6 ( R ) is obtained from ε ( R ) using the following:
ε ( R )   =   h ν 4 π n H , 6 ( R ) A ,
where h is the Planck constant and ν and A are the frequency and spontaneous emission coefficient of the Hδ transition, respectively. Then, n H ( R ) is obtained from Equations (3) and (4) as follows:
n H ( R )   =   4 π h ν A r 16 ( R ; T e , n e ) n e ( R ) ε ( R ) .
We interpolated T e and n e measured by Thomson scattering [9] in the same way as I A as shown in Figure 5a. In addition, the data was linearly extrapolated on the inboard side (R = 0.28−0.34 m) to evaluate n H around the resonance layer. The extrapolated curves are shown by dashed lines.
n H ( R ) evaluated with the above procedures is shown in Figure 5b. The density is in the range of 1015–1016 m−3 and monotonically decreases toward the ECH resonance layer. The increase around R = 0.65 m is physically implausible and may have been caused by an error in ε ( R ) . The evaluated density is less than the hydrogen molecule density of approximately 1017 m−3 at the vacuum chamber wall (R = 1.3 m) obtained from the pressure measured by a fast ionization gauge assuming a temperature of 300 K.
The accuracy of the CR model analysis can be affected by the opacity of Lyman series emission lines (p = 1–q), where q is the upper-state principal quantum number [14]. When the opacity is significant, the reabsorption of Lyman series emission lines increases the excited hydrogen atom density, and ignoring the effect results in the overestimation of n H ( R ) . Since the present CR model analysis is mainly affected by the Lε line (p = 1–6) and the secondary effect, namely, the increase in n H , 6 ( R ) owing to the absorption of the other Lyman series emission lines, is small [14,15], we evaluated the effect of the opacity by calculating the optical depth of the Lε line at its line center [16] under the assumption of a uniform plasma as follows:
τ = e 2 4 m e ϵ 0 m H 2 π k B T H f 16 ν 16 n H L ,
where e is the elementary charge, ε0 is the vacuum permittivity, kB is the Boltzmann constant, me is the electron mass, mH is the hydrogen atom mass, f16 (=7.8035 × 10−3) and ν16(=3.1967 × 1015 Hz) are the absorption oscillator strength and transition frequency, respectively, nH is the density of ground-state atoms, and L is the absorption length. We assumed an identical temperature TH for the excited and ground-state atoms. Assuming TH = 300 K, nH = 1 × 1016 m−3, and L = 2 m gives the upper limit of the optical depth τ as 1.5 × 10−2. The reabsorption with an optical depth of this magnitude is negligible in the CR model analysis.
Finally, as a measure of the effect of atoms on the plasma, we evaluated the degree of ionization (DOI). The DOI is defined as n H + / ( n H + + n H + n H 2 ) and we approximated it to n e / ( n e + n H ) on the basis of the above-mentioned assumptions, where we used the interpolated curve of n e in Figure 5a. The evaluated DOI shown in Figure 5b reached nearly unity over the entire core region, suggesting that the effect of atoms on the core plasma is small. The DOI, however, decreased in the SOL region and this could induce the momentum loss of the ions trapped in orbits that intersect the LCFS.

4. Conclusions

In this study, in order to evaluate the ion momentum loss and transport through the charge-exchange reactions in a QUEST plasma sustained only with 28 GHz ECH, we measured the radial distribution of the hydrogen atom density. The obtained density was of the order of 1015–1016 m−3, and considering the fact that it is less than the hydrogen molecule density evaluated at the vacuum chamber wall and the observed tendency of a monotonic decrease toward the ECH resonance layer, the result is plausible. However, the effect of the violation of the assumed toroidally symmetric emission should be addressed to further improve the reliability of the result. The DOI evaluated using the atom density and n e measured by Thomson scattering suggested that the effect of atoms is small in the core region.

Author Contributions

Conceptualization, T.S.; methodology, S.M., T.S., S.K., M.H. (Masahiro Hasuo), and K.S.; investigation, S.M., T.S., K.H., N.Y., A.K., and M.H. (Masahiro Hasuo); QUEST operation, K.H., H.I., T.O., T.I., K.N., R.I., Y.N., M.H. (Makoto Hasegawa), K.K., A.H., T.N., and S.S.; Thomson scattering measurements, A.E., Y.O., Y.P., and K.M., writing—original draft preparation, S.M.; writing—review and editing, T.S.; funding acquisition, T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by JSPS KAKENHI (18K03576), bilateral (NIFS16KUTR114, NIFS19KUTR140), and LHD project research (NIFS16KOAP031) collaboration programs of the National Institute for Fusion Science, and the collaborative research program of the Research Institute for Applied Mechanics, Kyushu University (30FP-38).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hanada, K.; Yoshida, N.; Hasegawa, M.; Hatayama, A.; Okamoto, K.; Takagi, I.; Hirata, T.; Oya, Y.; Miyamoto, M.; Oya, M.; et al. Particle balance investigation with the combination of the hydrogen barrier model and rate equations of hydrogen state in long duration discharges on an all-metal plasma facing wall in QUEST. Nucl. Fusion 2019, 59, 076007. [Google Scholar] [CrossRef]
  2. Idei, H.; Kariya, T.; Imai, T.; Mishra, K.; Onchi, T.; Watanabe, O.; Zushi, H.; Hanada, K.; Qian, J.; Ejiri, A.; et al. Fully non-inductive second harmonic electron cyclotron plasma ramp-up in the QUEST spherical tokamak. Nucl. Fusion 2017, 57, 126045. [Google Scholar] [CrossRef]
  3. Kuroda, K.; Raman, R.; Hanada, K.; Hasegawa, M.; Onchi, T.; Ono, M.; Nelson, B.A.; Jarboe, T.R.; Nagata, M.; Mitarai, O.; et al. Initial results from solenoid-free plasma start-up using transient CHI on QUEST. Plasma Phys. Control. Fusion 2018, 60, 115001. [Google Scholar] [CrossRef] [Green Version]
  4. Maekawa, T.; Yoshinaga, T.; Uchida, M.; Watanabe, F.; Tanaka, H. Open field equilibrium current and cross-field passing electrons as an initiator of a closed flux surface in EC-heated toroidal plasmas. Nucl. Fusion 2012, 52, 083008. [Google Scholar] [CrossRef] [Green Version]
  5. Yoneda, N.; Shikama, T.; Zushi, H.; Hanada, K.; Fujikawa, A.; Onchi, T.; Kuroda, K.; Nii, K.; Hasuo, M.; Hasegawa, M.; et al. Spectroscopic measurements of impurity ion toroidal and poloidal flow velocities and their dependence on vertical magnetic field in QUEST toroidal ECR plasmas. Plasma Fusion Res. 2018, 13, 3402087. [Google Scholar] [CrossRef]
  6. Yoshinaga, T.; Uchida, M.; Tanaka, H.; Maekawa, T. Spontaneous formation of closed-field torus equilibrium via current jump observed in an electron-cyclotron-heated plasma. Phys. Rev. Lett. 2006, 96, 125005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Yoneda, N.; Shikama, T.; Hanada, K.; Mori, S.; Onchi, T.; Kuroda, K.; Hasuo, M.; Idei, H.; Nakamura, K.; Nagashima, Y.; et al. Toroidal flow measurements of impurity ions in QUEST ECH plasmas using multiple viewing chords emission spectroscopy. Nucl. Mater. Energy 2020, submitted. [Google Scholar]
  8. Versloot, T.W.; de Vries, P.C.; Giroud, C.; Brix, M.; von Hellermann, M.G.; Lomas, P.J.; Moulton, D.; O’Mullane, M.; Nunes, I.M.; Salmi, A.; et al. Momentum losses by charge exchange with neutral particles in H-mode discharges at JET. Plasma Phys. Control. Fusion 2011, 53, 065017. [Google Scholar] [CrossRef] [Green Version]
  9. Yamaguchi, T.; Ejiri, A.; Hiratsuka, J.; Hasegawa, M.; Nagashima, Y.; Narihara, K.; Takase, Y.; Zushi, H.; The QUEST Group. Electron Temperature Measurement on QUEST Spherical Tokamak by Thomson Scattering System. Plasma Fusion Res. 2015, 10, 1202082. [Google Scholar] [CrossRef] [Green Version]
  10. Yoshinaga, T.; Uchida, M.; Tanaka, H.; Maekawa, T. A current profile model for magnetic analysis of the start-up phase of toroidal plasmas driven by electron cyclotron heating and current drive. Nucl. Fusion 2007, 47, 210. [Google Scholar] [CrossRef]
  11. Zushi, H.; Nakamura, K.; Hanada, K.; Sato, K.N.; Sakamoto, M.; Idei, H.; Hasegawa, M.; Iyomasa, A.; Kawasaki, S.; Nakashima, H.; et al. Steady-state tokamak operation, ITB transition and sustainment and ECCD experiments in TRIAM-1M. Nucl. Fusion 2005, 45, S142. [Google Scholar] [CrossRef]
  12. Ueda, A.; Shikama, T.; Teramoto, T.; Higashi, T.; Iida, Y.; Hasuo, M. Spectroscopic measurement of the degree of ionization in a helium electron cyclotron discharge in a simple cusp field. Appl. Phys. Lett. 2017, 111, 074101. [Google Scholar] [CrossRef]
  13. Sawada, K.; Fujimoto, T. Effective ionization and dissociation rate coefficients of molecular hydrogen in plasma. J. Appl. Phys. 1995, 78, 2913. [Google Scholar] [CrossRef] [Green Version]
  14. Sawada, K.; Goto, M.; Ezumi, N. Spectroscopic Diagnostic of Helium-Hydrogen RF Plasma under the Influence of Radiation Trapping. Plasma Fusion Res. 2011, 6, 1401010. [Google Scholar] [CrossRef] [Green Version]
  15. Ueda, A.; Shikama, T.; Teramoto, T.; Higashi, T.; Iida, Y.; Hasuo, M. Helium atom line-intensity ratios as an integrated diagnostic tool for low-pressure and low-density plasmas. Phys. Plasmas 2018, 25, 054508. [Google Scholar] [CrossRef] [Green Version]
  16. Shikama, T.; Ogane, S.; Iida, Y.; Hasuo, M. Measurement of the helium 23S metastable atom density by observation of the change in the 23S–23P emission line shape due to radiation reabsorption. J. Phys. D Appl. Phys. 2016, 49, 025206. [Google Scholar] [CrossRef]
Figure 1. (a) Temporal evolutions of the plasma current, electron temperature, and density at the core (R = 0.34 m, z = 0 m) (#41176), and chord-integrated Hα intensity (#41178). (b) Equilibrium magnetic flux surfaces in the time period t = 2.7–3.0 s denoted by the shaded area in (a). The last closed flux surface (LCFS) is shown by the red line.
Figure 1. (a) Temporal evolutions of the plasma current, electron temperature, and density at the core (R = 0.34 m, z = 0 m) (#41176), and chord-integrated Hα intensity (#41178). (b) Equilibrium magnetic flux surfaces in the time period t = 2.7–3.0 s denoted by the shaded area in (a). The last closed flux surface (LCFS) is shown by the red line.
Atoms 08 00044 g001
Figure 2. Schematic illustrations of the cross-section of the Q-shu University Experiment with Steady-State Spherical Tokamak (QUEST) midplane and the visible spectroscopic system.
Figure 2. Schematic illustrations of the cross-section of the Q-shu University Experiment with Steady-State Spherical Tokamak (QUEST) midplane and the visible spectroscopic system.
Atoms 08 00044 g002
Figure 3. (a) Hδ spectrum measured on the viewing chord at η = 0.51 m (#41178). (b) Chord-integrated Hδ intensities and the interpolated curve. (c) Radial distribution of Hδ emissivity obtained using Abel inversion.
Figure 3. (a) Hδ spectrum measured on the viewing chord at η = 0.51 m (#41178). (b) Chord-integrated Hδ intensities and the interpolated curve. (c) Radial distribution of Hδ emissivity obtained using Abel inversion.
Atoms 08 00044 g003
Figure 4. Population coefficients for p = 6. The green, red, and blue lines represent r 06 , r 16 , and r 26 , respectively. The solid, dashed, and dash-dotted lines represent those at n e = 1016, 1017, and 1018 m−3, respectively.
Figure 4. Population coefficients for p = 6. The green, red, and blue lines represent r 06 , r 16 , and r 26 , respectively. The solid, dashed, and dash-dotted lines represent those at n e = 1016, 1017, and 1018 m−3, respectively.
Atoms 08 00044 g004
Figure 5. (a) T e and n e   measured by Thomson scattering and their interpolated and extrapolated curves. (b) Evaluated n H and the degree of ionization (DOI).
Figure 5. (a) T e and n e   measured by Thomson scattering and their interpolated and extrapolated curves. (b) Evaluated n H and the degree of ionization (DOI).
Atoms 08 00044 g005

Share and Cite

MDPI and ACS Style

Mori, S.; Shikama, T.; Hanada, K.; Yoneda, N.; Kuzmin, A.; Hasuo, M.; Idei, H.; Onchi, T.; Ejiri, A.; Osawa, Y.; et al. Spectroscopic Measurement of Hydrogen Atom Density in a Plasma Produced with 28 GHz ECH in QUEST. Atoms 2020, 8, 44. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms8030044

AMA Style

Mori S, Shikama T, Hanada K, Yoneda N, Kuzmin A, Hasuo M, Idei H, Onchi T, Ejiri A, Osawa Y, et al. Spectroscopic Measurement of Hydrogen Atom Density in a Plasma Produced with 28 GHz ECH in QUEST. Atoms. 2020; 8(3):44. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms8030044

Chicago/Turabian Style

Mori, Satoru, Taiichi Shikama, Kazuaki Hanada, Nao Yoneda, Arseniy Kuzmin, Masahiro Hasuo, Hiroshi Idei, Takumi Onchi, Akira Ejiri, Yuki Osawa, and et al. 2020. "Spectroscopic Measurement of Hydrogen Atom Density in a Plasma Produced with 28 GHz ECH in QUEST" Atoms 8, no. 3: 44. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms8030044

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop