-
With the detailed consideration of electrochemical reactions and collision relations, a direct numerical simulation model of helicon plasma discharge with three-dimensional fluid-dynamic equations is proposed in the present work. It can improve the precision of results and widen the model applicability by discarding the small perturbation theory in previous helicon models which are partially analytical in essence. Under the assumption of weak ionization, the Maxwell equations coupled with the plasma parameters are directly solved in the whole computational domain. Thus the energy deposited from electromagnetic wave to plasma can be then easily calculated. The values of plasma parameters which include electron density, mean electron energy and heavy species density are obtained by solving a set of drift-diffusion equations. Meanwhile, seven kinds of chemical reactions in the plasma and two kinds of surface reactions on the wall are taken into account. All of the partial differential equations are solved by the finite element solver of COMSOL MultiphysicsTM with the full coupled method.#br#The results of numerical cases employing argon as the working medium show that there exists a sharp density jump from a low to high value as the radiofrequency power is raised. The density jump phenomenon is in accordance with the experimental results of Toki (Toki K, Shinohara S, Tanikawa T, Shamrai K P 2006 Thin Solid Films 506-507 597) and Chen (Chen F F 2007 Plasma Sources Sci. Technol. 16 593). The electron temperature decreases with an increase of the gas pressure, which is similar to Toki's (Toki K, Shinohara S, Tanikawa T, Shamrai K P 2006 Thin Solid Films 506-507 597) measurement by a RF compensation probe. In comparison with the classical sheath theory, the simulation result demonstrates that the distribution of parameters such as particle number density, the Deby length, electric potential and electron temperature can be solved exactly. In addition, the phenomenon of low-field density peak in helicon discharge was studied in the work. Previous research by Chen (Chen F F 2003 Phys. Plasmas 10 2586) suggests that this peak is caused by constructive interference from the reflected wave. The effect of length of the discharge chamber on the relation of electron density and background magnetic field is investigated numerically. The results validate the mechanism of wave interference reflected by endplates of the discharge chamber. Furthermore, the time-averaged magnetic energy density has more than one peak on the axial direction. Comparing the distribution of the magnetic energy density to that of the dimensionless amplitude of the helicon wave and the TG wave in the one-dimensional undamped condition, it found that the length of peak to peak of the helicon wave is just as twice as that of the magnetic energy density, which indicates that the substance of wave interference is involved in the standing wave generated by the helicon wave and its reflected wave from endplates.
-
Keywords:
- plasma /
- helicon wave /
- direct numerical simulation /
- wave interference
[1] Chen F F 1991 Plasma Phys. Controlled Fusion 33 339
[2] Suwon C 1996 Phys. Plasmas 3 4268
[3] Suwon C 1997 Phys. Plasmas 4 4167
[4] Arnush D, Chen F F 1998 Phys. Plasmas 5 1239
[5] Arnush D 2000 Phys. Plasmas 7 3042
[6] Chen F F 2008 IEEE Trans. Plasma Sci. 36 2095
[7] Curreli D, Chen F F 2011 Phys. Plasmas 18 113501
[8] Chen F F, Curreli D 2013 Phys. Plasmas 20 057102
[9] Ahedo E 2009 Phys. Plasmas 16 113503
[10] Ahedo E, Jaume N C 2013 Phys. Plasmas 20 043512
[11] Cheng Y G, Cheng M S, Wang M G, Li X K 2014 Acta Phys. Sin. 63 035203 (in Chinese)[成玉国, 程谋森, 王墨戈, 李小康2014 63 035203]
[12] Cheng Y G, Cheng M S, Wang M G, Li X K, Yang X 2014 Plasma Sources Sci. Technol. 16 1111
[13] Boswell R W 1984 Plasma Phys. Controlled Fusion 26 1147
[14] Toki K, Shinohara S, Tanikawa T, Shamrai K P 2006 Thin Solid Films 506-507 597
[15] Chen F F 2007 Plasma Sources Sci. Technol. 16 593
[16] Cheng Y G 2015 Ph. D. Dissertation (Changsha:National University of Defense Technology) (in Chinese)[成玉国2015博士学位论文(长沙:国防科技大学)]
[17] Hagelaar G J M, Pitchford L C 2005 Plasma Sources Sci. Technol. 14 722
[18] Dimitris P L, Demetre J E 1995 J. Res. Nat. Inst. Stand. Technol. 100 473
[19] Chen F F, Torreblanca H 2007 Plasma Phys. Controlled Fusion 49 81
[20] Chen F F 1992 J. Vac. Sci. Technol. A 10 1389
[21] Chen F F 2003 Phys. Plasmas 10 2586
[22] Chen F F, Arnush D 1997 Phys. Plasmas 4 3411
-
[1] Chen F F 1991 Plasma Phys. Controlled Fusion 33 339
[2] Suwon C 1996 Phys. Plasmas 3 4268
[3] Suwon C 1997 Phys. Plasmas 4 4167
[4] Arnush D, Chen F F 1998 Phys. Plasmas 5 1239
[5] Arnush D 2000 Phys. Plasmas 7 3042
[6] Chen F F 2008 IEEE Trans. Plasma Sci. 36 2095
[7] Curreli D, Chen F F 2011 Phys. Plasmas 18 113501
[8] Chen F F, Curreli D 2013 Phys. Plasmas 20 057102
[9] Ahedo E 2009 Phys. Plasmas 16 113503
[10] Ahedo E, Jaume N C 2013 Phys. Plasmas 20 043512
[11] Cheng Y G, Cheng M S, Wang M G, Li X K 2014 Acta Phys. Sin. 63 035203 (in Chinese)[成玉国, 程谋森, 王墨戈, 李小康2014 63 035203]
[12] Cheng Y G, Cheng M S, Wang M G, Li X K, Yang X 2014 Plasma Sources Sci. Technol. 16 1111
[13] Boswell R W 1984 Plasma Phys. Controlled Fusion 26 1147
[14] Toki K, Shinohara S, Tanikawa T, Shamrai K P 2006 Thin Solid Films 506-507 597
[15] Chen F F 2007 Plasma Sources Sci. Technol. 16 593
[16] Cheng Y G 2015 Ph. D. Dissertation (Changsha:National University of Defense Technology) (in Chinese)[成玉国2015博士学位论文(长沙:国防科技大学)]
[17] Hagelaar G J M, Pitchford L C 2005 Plasma Sources Sci. Technol. 14 722
[18] Dimitris P L, Demetre J E 1995 J. Res. Nat. Inst. Stand. Technol. 100 473
[19] Chen F F, Torreblanca H 2007 Plasma Phys. Controlled Fusion 49 81
[20] Chen F F 1992 J. Vac. Sci. Technol. A 10 1389
[21] Chen F F 2003 Phys. Plasmas 10 2586
[22] Chen F F, Arnush D 1997 Phys. Plasmas 4 3411
Catalog
Metrics
- Abstract views: 7877
- PDF Downloads: 391
- Cited By: 0