Plasma heating due to X-B mode conversion in a cylindrical ECR plasma system arxiv:physics/41112v2 [physics.plasm-ph] 6 Nov 24 Vipin K. Yadav and D. Bora Institute for Plasma Research, Bhat, Gandhinagar, Gujarat - 382 428, India Abstract Extra Ordinary (X) mode conversion to Bernstein wave near Upper Hybrid Resonance (UHR) layer plays an important role in plasma heating through cyclotron resonance. Wave generation at UHR and parametric decay at high power has been observed during Electron Cyclotron Resonance (ECR) heating experiments in toroidal magnetic fusion devices. A small linear system with ECR and UHR layer within the system has been used to conduct experiments on X-B conversion and parametric decay process as a function of system parameters. Direct probing in situ is conducted and plasma heating is evidenced by soft x-ray emission measurement. Experiments are performed with hydrogen plasma produced with 16-8 W microwave power at 2.45 GHz of operating frequency at 1 3 mbar pressure. The axial magnetic field required for ECR is such that the resonant surface (B = 875 G) is situated at the geometrical axis of the plasma system. Experimental results will be presented in the paper. I. INTRODUCTION Electron cyclotron resonance heating (ECRH) is a highly efficient, controllable and localized plasma heating [1] and current drive scheme in fusion machines such as tokamaks [2 5], 1
stellerators[6] and plasma systems such as tandem mirrors[7]. These devices, when opearted at high plasma densities (overdense plasma, ω pe ω ce ), direct ECR heating and current drive by the incident electromagnetic (EM) wave is not possible. Alternatively, linear conversion of the EM waves into the electrostatic (ES) EBW [8] at the UHR surface and subsequent damping at the cyclotron harmonics is being tried for overdense plasmas [9]. In ECRH, the launched EM wave in extraordinary (X) mode, ( k B and E 1 B, k is the wave vector, B is the magnetic field in the system and E 1 is the wave electric field), is absorbed at the UHR surface in plasma as shown in the X-mde dispersion relation c 2 k 2 = 1 ω2 pe ω 2 ω 2 ( ω 2 ωpe 2 ) ω 2 ω 2 uh (1) where, ω is the launched microwave frequency, ω pe is the electron plasma frequency and ω uh is the upper hybrid frequency given by ω 2 uh = ω 2 pe +ω 2 ce (2) here, ω ce is the electron cyclotron frequency. Near the UHR layer, the X-mode wavelength becomes very short and the amplitude of wave electric field becomes very large (high power-density). This large wave electric field gives rise to nonlinear phenomena of parametric excitation in which the incoming EM wave breaks up into EBW (X-B mode conversion) and lower hybrid (LH) wave [11]. The EBW propagates perpendicularly to the magnetic field towards the plasma centre and gets strongly damped at ECR layer as shown by the EBW dispersion relation given by 1+ ( kb v th ω pe ) 2 ( ) ( ) = e ( k2 B r2 L) I k 2 ω B rl 2 2 e ( k2 ωce 2 B r2 L) I q(k 2 ( B r2 L ) ) (3) q q2 ω2 ωce 2 here, v th is the electron thermal velocity, r L is the electron Larmour Radius and I q is the Bessel function of the first kind of imaginary argument, q = 1, 2... The absorption of EBW contribute in efficient plasma heating [12 14] The parametric decay spectrum consists of wave with frequency of the incident pump wave, ω andω±ω LH. Thecharacteristicsidebandsinthehigh-frequencyspectrumandtheassociated 2
low-frequency waves from the parametric excitation have been observed in tokamaks[15] and stellarators [16,17]. The EBW is observed experimentally in linear or toroidal plasma system [18,13]. The EBW generation via mode conversion is observed in a large diameter linear plasma system with filament aided discharge [12]. II. EXPERIMENTAL SET-UP The experimental setup consists of a cylindrical stainless steel vacuum chamber with internal radius 6.4 cm and axial length 1 cm. A base pressure of 6 1 6 mbar is achieved prior to the experiment using a diffusion pump. The magnetic field is produced by two identical circular coils rested on the axial ports at the two ends of the vacuum vessel. Each magnetic field coil has an electrical resistance of 3 Ω and require 5.3 A of dc current to produce magnetic field of 875 G at the axis of the vacuum vessel which is kept constant for more than 1 second. The magnetic field contours in the system are shown in figure 2 [19]. The microwave system consists of a 8 W, 2.45 ±.2 GHz magnetron, a three port circulator with water cooled dummy load to dissipate the reflected power, WR34 directional coupler, low barrier Schottky diode microwave detectors to measure forward and reflected power and a 1 mm thick, 11 mm 7 mm rectangular toughen glass sheet to serve as the waveguide window. In absence of plasma 8 % of the forward power gets reflected whereas the reflected power drops to 2 % of the forward power in presence of plasma. The plasma parameters are measured by a movable Langmuir probe [2] in the radial direction. The probe is made of tungsten with tip length and diameter of 5 mm and.5 mm respectively. A capacitive probe [21] as shown in figure 3, is used to measure the wave field. Probe signal is directly analysed with a spectrum analyser. Two capacitive probes [22] are used to measure the wave characteristics having tip length less than the expected wavelength to be detected, estimated theoretically as 7.5 mm. The signals from the capacitive probes are fed into two double balanced mixers (DBM). Each mixer has three ports : two input - 3
RF and LO (local oscillator) ports and one output - IF (intermediate frequency). III. EXPERIMENTAL RESULTS ECR plasma is produced in the experimental system with hydrogen as the filling gas at an operating pressure of 1 1 3 mbar and an input microwave power of 8 W. The maximum plasma density, n e in the system is 3.1 1 1 cm 3 at the geometrical axis of the system. Typical plasma temperature, n e in the system is 12 ev. Radial profiles of f ce, f pe and f uh are calculated and plotted for the measured density and magnetic field profiles as shown in figure 4. Plasma density during the discharge for microwave power of 8 W is such that the UHR layer during the experiment lies at R = 2 cm. The f uh variation with input power is plotted in figure 5 which shows that despite the change in input power, the UHR layer postion remains unchanged. TheparametricdecayspectrumatUHRlayerisshowninfigure7. TheresidualLHfrequency is observed separately as shown in figure 8 which is calculated as 23.4 MHz hence, it matches well with the measured peaks. The capacitive probe is moved in hydrogen plasma from R = cm towards the UHR layer at R = 2 cm. Figure 9 shows the capacitive probe output at the UHR layer as a function of microwave power that is increased gradually from 16 W to the maximum possible value of 8 W. As the power is increased, the spectrum is broadened on the lower frequency edge and at 3 W a peak starts appearing at 2.442 GHz which develops as the power is increased, Power threshold for excitation of the sideband is measured and the peak at 2.442 GHz appears only beyond 5 W of input power [21]. At the peak power, the signal is measured along the radius. As the capacitive probe is placed close to thr UHR layer, within R = 1 cm to R = 2 cm, the original spectrum splits to smaller peaks as shown in figure 1. At UHR, the peak at f = 2.442 GHz is distinct. The figure indicates the two ECR surfaces at R = and R = 6 where the frequency spectrum is clean without fluctuations. The spectrum at R =2 (UHR surface) shows the side-peaks with the launched frequency indicating the parametric decay there. The frequency spectrum at 4
the radial positions near UHR surface, R = 1 and R = 3 contains fluctuations suggesting the splitting of power at UHR and propagation in both forward and reverse directions [21]. The EBW is detected and characterized with two identical capacitive probes [22]. The axial capacitive probe is fixed at R = cm and the radial capacitive probe is moved upto UHR layer. This small gapof 2 cm between the UHR and ECR layers limits the spatial resolution. The EBW parameters are given in the table below. S.No. Parameter Value 1. Gas Hydrogen 2. f GHz 2.45 ±.2 3. f EBW GHz 2.438 4. λ B mm 6.56 5. v φ ms 1 1.6 1 7 6. v g ms 1 3.88 1 4 7. n 2 The mode converted EBW propagate towards the centre of the experimental system near which the first ECR layer is residing. The EBW is absorbed at the ECR surface and gives rise to the localized electron heating. This is evident from the time evolution of plasma temperature as shown in figure 11 where the temperature increases with time and then attains a saturated value 5 ev afterwards for 8 W of input power. This indicates the absorption of the launched EM wave in the system. The localized heating of plasma with mode converted EBW in this experimental plasma system is further verified by the observation of the emission of soft x-rays [24]. The EEDF obtained from the Langmuir probe characteristics at the centre of the plasma chamber is shown in figure 12. Here, (V p - V B ) is the accelerating voltage for the electrons in the plasma and di/dv signifies the number of such electrons. In the figure, the hump in the tail of the distribution indicates the presence of high energy electrons with energies more than 1 ev. These energetic electrons hit metal targets of tungsten, molybdenum, etc. in 5
the plasma and leads to emission of soft x-ray. X-ray radiation from plasma can be energy resolved by using filters in the form of thin Aluminium foils. In this experiment also the soft x-ray emission due to the generation of energetic electrons by the absorption of mode converted EBW is observed [24]. The soft x-ray emission signal with and without any target inside the plasma is shown in figure 13. This signal is very feeble compared to the one with Langmuir probe tip made of thoriated tungsten placed in the plasma. This indicates that the presence of a metal target increases the soft x-ray emission. The soft x-ray emission also depends on the nature of the metal target as shown in figure 14 where the signals for both the targets differ. IV. CONCLUSION Plasma heating due to X-B mode conversion is observed in a linear ECR plasma system. The X-B mode conversion is nonlinear in nature associated with parametric decay of the launched wave frequency. Parametric decay spectrum of the launched X-mode is observed at the UHR layer in system with capacitive probe. The decay spectrum contains ω ±ω LH about the launched frequency, ω (2.45 +.2 GHz) indicating the three-wave nonlinear interaction. The experimentally observed value of ω LH in the spectrum matches quite well with the estimated values. The parametric decay spectrum has an input power threshold ( 6 W) to occur. This is in accordance with the estimated input power values. The parametric decay peaks are observed only near the calculated UHR layer in the system. The shift in probe position from UHR vanishes the peaks. The EBW generated near the UHR surface as a result of the mode conversion is experimentally observed and characterized. The EBW then propagates towards the first harmonic cyclotron resonance layer residing near the centre of the geometrical axis of the system and is absorbed. This absorption leads to the localized electron heating at that position. Experimental results indicate the presence of high energy electrons. This is further verified by EEDF obtained from the I-V characteristic at the centre of the system. The soft x-ray emission is observed from the ECR plasma using 6
a VPD. The soft x-ray signal is feeble in the absence of a target and gets enhanced in the presence of a tungsten target. The signal is also observed to be a dependent function of the material of the target. 7
E-mail address: vipin@ipr.res.in REFERENCES E-mail address: dbora@ipr.res.in [1] Brian Lloyd, Plasma Phys. Cont. Fusion, 4, A119 (1998) [2] Wallace M. Manheimer, Infrared and Millimeter Waves Vol. 2 : Instrumentation, (Academic Press, New York, 1979), chapter 5, 299 [3] A.C. England, IEEE Trans. Plasma Sci., PS-12, 124 (1984) [4] O.C. Eldridge and A.C. England, Nucl. Fusion, 29, 1583 (1989) [5] V. Erckmann and U. Gasparino, Plasma Phys. Cont. Fusion, 36, 1869 (1994) [6] V.V. Abrakov, D.K. Akulina, Eh. D. Andryukhina, et. al., Nucl. Fusion, 37, 233 (1997) [7] M. Porkolab, L. Friedland and I.B. Bernstein, Nucl. Fusion, 21, 1463 (1981) [8] Ira B. Bernstein, Phys. Rev., 19, 1 (1958) [9] A.D. Piliya, A. Yu Popov and E.N. Tregubova, Plasma Phys. Control. Fusion, 45, 139 (23) [1] Thomas H. Stix, Phys. Rev. Lett., 15, 878 (1965) [11] Ya. N. Istomin and T.B. Leyser, Phys. Plasmas, 2, 284 (1995) [12] H. Sugai, Phys. Rev. Lett., 47, 1899 (1981) [13] S. Nakajima and H. Abe, Phys. Rev. A, 38, 4373 (1988) [14] A.K. Ram and S.D. Schultz, Phys. Plasmas, 7, 484 (2) [15] F.S. McDermott, G. Bekeffi, K.E. Hackett, J.S. Levine and M. Porkolab, Phys. Fluids, 25, 1488 (1992) [16] H.P. Laqua, H. Maassberg, N.B. Marushchenko, F. Volpe, A. Weller, W7-AS Team, W. 8
Kasparek and ECRH-Group, Phys. Rev. Lett., 9, 753-1 (23) [17] H.P. Laqua, V. Erckmann, H.J. Hartfuβ, H. Laqua, W7-AS Team and ECRH group, Phys. Rev. Lett., 78, 3467 (1997) [18] W. Gekelman and R.L. Stenzel, Rev. Sci. Instrum., 46, 1386 (1975) [19] Vipin K. Yadav and D. Bora, Pramana, 63, 563, (24) [2] I. Langmuir, Phys. Rev., 36, 954 (1929) [21] Vipin K. Yadav and D. Bora, Phys. Plasmas, 11, 349 (24) [22] Vipin K. Yadav and D. Bora, Phys. Plasmas, 11, 4582 (24) [23] C.V.S. Rao, Y. Shankara Joisa and C.J. Hansalia, A.K. Hui, Ratan Paul and Prabhat Ranjan, Rev. Sci. Instrum., 68, 1142 (1997) [24] Vipin K. Yadav and D. Bora, Plasma Sources Sci. Technol., 13, 231, (24) 9
FIGURES Circular Port Magnetic Coils Main Chamber Magnetic Coils Port 1 Rectangular Port Port 2 Magnetic Coils Magnetic Coils Circular Port FIG. 1. Main chamber. 6 438 4 6 75 2 1.6e+3 1e+3 1.6e+3 R (cm) 2 1.4e+3 1.2e+3 1e+3 875 875 75 1.2e+3 1.4e+3 4 6 6 438 4 6 8 1 Z (cm) FIG. 2. Magnetic field contours in the system. 8 mm SS tube Glass cap Graphite 1 Silver coating 1 1 Teflon ring 1111111111111111 11111111111111111 11111111111111111 11111111111111111 11111111111111111111111111111111 1 Probe tip 11 1111 11111111111111111 11111111111111111 11111111111111111 11111111111111111 1111111111111111 11 11 7 cm 3 mm Semi rigid co axial cable BALUN section 2.1 cm for teflon for 2.45 GHz FIG. 3. Capacitive probe. 1
3 f uhr 2.5 f RF f (GHz) 2 1.5 f ce 1 f pe.5 6 4 2 2 4 6 R (cm) FIG. 4. Radial variation of f ce, f ce and f uhr in hydrogen plasma. 2.8 2.6 2.4 f RF f uhr (GHz) 2.2 p = 1 3 mbar 2 P RF = 8 W : x 1.8 P RF = 6 W : + 1.6 P RF = 4 W : o 1.4 6 4 2 2 4 6 R (cm) FIG. 5. f uhr with input microwave power variation in hydrogen plasma. 2 1 < 2.47 1 2 3 4 FIG. 6. Power output from magnetron at forward port of directional coupler. 11
< 2.475 1 2 3 2.47 > 2.442 > < 2.49 < 2.58 4 2.3 2.35 2.4 2.45 2.5 2.55 2.6 FIG. 7. Parametric decay spectrum in hydrogen plasma. 28 Amplitude (A.u.) 3 32 34 < 37.9 36 2 3 4 5 6 7 Frequency (MHz) FIG. 8. The residual low frequency. 12
2 P = 8 W 2 P = 4 W 4 4 2 P = 7 W 2 P = 3 W 4 4 2 P = 6 W 2 P = 225 W 4 4 2 P = 5 W 2 P = 16 W 4 4 FIG. 9. Parametric decay spectrum with input microwave power variation. 13
2 r = 2 r = 4 4 4 2 r = 1 2 r = 5 4 4 2 r = 2 2 r = 6 4 4 2 4 r = 3 FIG. 1. Radial variation of parametric decay spectrum at 8 W input microwave power. 6 5 p = 1 3 mbar 4 P = 8 W T e (ev) 3 2 1 P = 4 W 5 1 15 2 25 Time (msec) FIG. 11. Time evolution of plasma temperature in hydrogen plasma. 14
.5.4 di/dv.3.2.1 5 1 15 [V p V B ], (V) FIG. 12. The electron energy distribution function (EEDF). Soft x ray signal (A.u.).4.3.2.1 with target without target Pressure = 5 X 1 5 mbar Power = 8 Watt.2.4.6.8 Time (sec) FIG. 13. Soft x-ray signal with and without metal target in plasma..3 (A.u.).2.1 Al wall.2.4.6.8 SXR signal.3.2.1 pressure = 5 X 1 5 mbar Power = 8 W SS wall.2.4.6.8 Time (sec) FIG. 14. Soft x-ray signal variation with Al and SS wall. 15