簡易檢索 / 詳目顯示
| 研究生： |
林庭億 Ting-Yi Lin |
|---|---|
| 論文名稱： |
一維羅倫茲電漿粒子模擬的動力學特性 Kinetic Properties of the Particle-in-Cell Simulation of one dimensional Lorentz Plasmas |
| 指導教授： |
陳仕宏
S. H. Chen |
| 口試委員: | |
| 學位類別： |
碩士 Master |
| 系所名稱： |
理學院 - 物理學系 Department of Physics |
| 畢業學年度： | 98 |
| 語文別： | 中文 |
| 論文頁數： | 85 |
| 中文關鍵詞： | 數值熱化 、羅倫茲電漿 、熱化 |
| 外文關鍵詞： | Lorentz plasma, numerical thermalization, Particle-in-Cell, PIC |
| 相關次數： | 點閱：9 下載：0 |
| 分享至: |
| 查詢本校圖書館目錄 查詢臺灣博碩士論文知識加值系統 勘誤回報 |
在粒子模擬法中數值熱化的現象(即速度分佈變成馬克士威爾)已經被廣泛的進行，在早期階段，其發展[1]已被認為很好地解釋。然而，最近在模擬弱游離氣體的模擬碰撞研究[2](這個碰撞模型描述電子與背景粒子(如原子和分子)的碰撞，當電子與中性背景氣體碰撞，電子將會被反彈。)發現當存在是一個額外的隨機力作用於粒子，對於之前對熱化時間的描述可能需要改變。在我們的工作中，我們感興趣的問題是，電子-離子在完全游離電漿中所發生的碰撞。我們研究在一維羅倫茲電漿(Lorentz plasma)[3]中粒子模擬法的熱化過程。所用到的小角度碰撞(small-angle collition)模型是使用蒙地卡羅演算法讓電子對固定的離子背景發生碰撞。我們的計算結果表明，熱化時間與 (每個德拜長度包含粒子的數目)成正比，而不是一般標準粒子模擬所認同的 。另一方面，我們的研究結果還指出，碰撞率對自我加熱機制(self-heating)的影響並不大，即使在發生強烈碰撞的環境也一樣。我們的研究結果可以補充那些先前被發現的研究[2]。
The phenomenon of numerical thermalization (i.e., relaxation of the velocity distribution toward a Maxwellian) in the standard particle-in-cell (PIC) simulation of Vlasov plasmas has been extensively studied at the early stage of its development[1] and was considered well understood. However, it was recently reported[2] that the well-established scaling law for the thermalization time could be compromised by the presence of an additional stochastic force acting on the particles, which has been used to simulate collisional processes in a weakly ionized gas. The collision model is described that the electrons collide with background particles, i.e., atoms and molecules. When electron collide with neutral background gas, the electron will be rebounded. In the our work, we are interested in the problem of electron-ion collisions in a fully ionized plasma. We examined the thermal relaxation processes in the PIC simulation of a Lorentz plasma in one dimension[3]. The small-angle collision of the electrons by the stationary ion background is modeled by a Monte-Carlo algorithm. Our numerical results show that the thermal relaxation time is proportional to ND (the number of particles per Debye length), but not ND2 as shown in the standard PIC simulations. On the other hand, our results also point out that the collision rate affects self-heating slightly even in a strong collision environment. Our results appear to complement those found by the previous study[2].
[1] C. K. Birdsall and A. B. Langdon, Plasma Physics via Computer Simulation (McGraw-Hill, New York, 1985).
[2] M. M. Turner, Phys. of Plasmas 13, 033506 (2006).
[3] R. Shanny, J. M. Dawson, and J. M. Greene, Phys. of Fluids 10, 1281 (1967).
[4] C. K. Birdsall and A. b. Langdon, Plasma Physics via Computer Simulation, IOP Publishing Ltd, London, 1991.
[5] R. W. Hockney and J. W. Eastwood, Computer Simulation Using Particles, IOP Publishing Ltd, Bristol, 1988.
[6] C. K. Birdsall, IEEE Trans. Plasma Science, vol. 19, p. 65, 1991.
[7] J. M. Dawson, Rev. Mod. Phys. Vol. 55, p. 403, 1983.
[8] C. K. Birdsall and N. Maron, J. Comp. Phys. Vol. 36, p. 1, 1980.
[9]J. M. Dawson, Phys. Fluids 7, 419 1964.
[10]J. Virtamo and H. Tuomisto, Phys. Fluids 22, 172 1979.
[12] F. F. Chen, Plasma Physics and Controlled Fusion (Plenum, New York, 1984).
[11]A. H. Boozer and G. Kuo-Petravic, Phys. Fluids 24, 851 (1981).
[13] Eldridge, O.C., Feix, M.R. (1963): Numerical Experiments with a Plasma Model. Phys.Fluids 6, 398-406.
[14] 等離子體粒子模擬－邵福球編著－北京：科學出版社，2002
