Modeling a high power fusion plasma reactor-type ion source: Applicability of particle methods

被引:29
作者
Fubiani, G. [1 ,2 ]
Hagelaar, G. J. M. [1 ,2 ]
Boeuf, J. P. [1 ,2 ]
Kolev, S. [3 ]
机构
[1] Univ Toulouse, UPS, INPT, LAPLACE Lab Plasma & Convers Energie, F-31062 Toulouse 9, France
[2] CNRS, LAPLACE, F-31062 Toulouse, France
[3] Univ Sofia, Fac Phys, BG-1164 Sofia, Bulgaria
关键词
ASYMPTOTIC PRESERVING SCHEME; TIME INTEGRATION; SIMULATION;
D O I
10.1063/1.3696037
中图分类号
O35 [流体力学]; O53 [等离子体物理学];
学科分类号
070204 ; 080103 ; 080704 ;
摘要
A simple two-dimensional (2D) Cartesian, three-velocities (3 V), particle-in-cell Monte Carlo collisions (PIC MCC) model of a tandem type inductively coupled plasma (ICP) discharge is presented. The conditions are similar to those of negative ion sources for fusion applications, i.e., a high absorbed power (on the order of 100 kW) and a high density plasma (typically 10(18) m(-3)) in a large volume vessel with a magnetic field barrier. We show that the plasma transport properties may be calculated with sufficient accuracy by performing simulations at low plasma density (low absorbed power) and rescaling the results to large plasma density by assuming linear variations of plasma density with power (in case of non-linear plasma chemistry, appropriate scale of the cross-sections or reaction rates is also necessary). This is equivalent to performing the simulations with the real absorbed power but using a larger value of the vacuum permittivity in Poisson's equation, to force quasineutrality. Rescaled results from explicit PIC MCC simulations compare well with those from direct implicit PIC MCC simulations (performed with real power and plasma density). Implicit simulations have less stringent constraints on time step and grid spacing but numerical heating may be more difficult to control and verification of energy conservation during the calculation is necessary. Rescaled explicit simulations may be employed for non-turbulent plasmas provided that the sheath length is small with respect to the plasma dimensions. The simulation results are used to discuss the physics of plasma transport across the magnetic filter and the role of the drift currents (E x B drift and diamagnetic drift). (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3696037]
引用
收藏
页数:13
相关论文
共 46 条
[1]   A Langmuir probe system for high power RF-driven negative ion sources on high potential [J].
不详 .
PLASMA SOURCES SCIENCE & TECHNOLOGY, 2009, 18 (01)
[2]  
[Anonymous], 2005, PRINCIPLES PLASMA DI, DOI [10.1002/0471724254, DOI 10.1002/0471724254]
[3]  
[Anonymous], THESIS U TOULOUSE FR
[4]  
[Anonymous], 1984, CHEN
[5]  
[Anonymous], 1983, Introduction to Plasma Theory
[6]  
[Anonymous], THESIS MIT BOSTON
[7]  
Birdsall C.K., 1998, Plasma Physics via Computer Simulation
[8]   Model of an inductively coupled negative ion source: II. Application to an ITER type source [J].
Boeuf, J. P. ;
Hagelaar, G. J. M. ;
Sarrailh, P. ;
Fubiani, G. ;
Kohen, N. .
PLASMA SOURCES SCIENCE & TECHNOLOGY, 2011, 20 (01)
[9]   AN IMPLICIT METHOD FOR ELECTROMAGNETIC PLASMA SIMULATION IN 2 DIMENSIONS [J].
BRACKBILL, JU ;
FORSLUND, DW .
JOURNAL OF COMPUTATIONAL PHYSICS, 1982, 46 (02) :271-308
[10]  
Briggs W.L., 2000, A Multigrid Tutorial, V2nd