电浆模拟.DOCVIP

電漿模擬 文/藍永強、邱行偉 電漿的動態行為非常複雜，主要透過自洽與外加的電磁場來交互作用，常涉及多重時間與空間耦合的非線性過程，因之，以電腦模擬來研究電漿的複雜現象與作用機制是必要也實際有效的方法。本文將介紹電漿研究中經常使用的主要模擬方法，包括粒子模擬以及電漿磁流體模擬。 (phase space)分佈函數隨時間的變化需滿足波茲曼傳輸方程式(Boltzmann transport equation)或是伏拉索夫方程式(Vlasov equation，在無碰撞電漿系統中)的描述[1]，在假設電漿滿足特定的分佈函數下，對波茲曼傳輸方程式作moment展開可以推導出流體力學方程式，電漿模擬方法中的流體模型(fluid model)即是根據流體力學方程式進行電漿的模擬[2-4]，然而還有一種更為基本的模擬方法，也就是根據Maxwell電磁方程式直接模擬電漿分佈函數隨著時間的變化，電漿模擬方法中的粒子模型即是以模擬電漿粒子的運動來得到電漿的分佈函數。 電漿的粒子式模擬方法最早是由Dawson、Buneman、HockneyBirdsall等學者於1950代中期起逐漸發展而來[5-8]，這種方法是直接模擬空間中荷電粒子的運動以研究電漿的特性以及電漿與電磁波的交互作用，由於電漿系統中的荷電粒子非常多(一個理想的電漿系統其空間尺度必須遠大於電漿的Debye長度，即，而在一個Debye球內電漿粒子數要足夠多以產生遮蔽效應，即)，如果要模擬每一個荷電粒子的運動，即使採用現今的大型電腦主機與平行計算也無法在合理的時間內獲得足夠的電漿分佈訊息，取而代之的是以模擬巨粒子(macro particles)的運動來獲得電漿的統計資料，每一顆巨粒子都代表模擬空間中一群真實的荷電粒子，而且巨粒子為一具有固定大小的粒子雲。能夠使用巨粒子來取代真實電漿的原因在於，在電漿系統中我們有興趣的物理機制通常是電漿的集合運動行為(collective behavior)，而非個別電漿粒子的運動行為，當電漿的集合運動行為其特徵尺度大於電漿的Debye長度時，個別電漿粒子的運動行為將變得不重要，因此可以使用巨粒子來模擬電漿粒子的集合運動，也正因為電漿粒子的這項特性，我們才得以用有限的電腦資源與時間模擬電漿粒子的動態行為；至於採用固定大小粒子雲的原因則在於，如果以點狀粒子模擬荷電粒子的運動，當兩個粒子距離非常小時，粒子間的庫倫作用力會非常大()，另一方面，在電漿系統中由於Debye遮蔽效應，電漿間發生小角度碰撞的機率遠大於大角度的庫倫散射，採用固定大小的電子雲可以消除電漿在近距離發生庫倫作用力發散的狀況並避免大角度的庫倫散射。 粒子式模型還必須能自洽地(self-consistently)模擬電漿與電磁場的交互作用。在電磁場的計算部份，粒子式模型通常採用有限差分時域(finite-difference time-domain)模型與Yee格網(mesh)以計算空間中電磁場隨時間的變化，而在Maxwell方程式的電流密度項中，除了外加電流密度外，也要加入因電漿粒子的運動所產生的電流密度，近年來也有學者採用有限元素(finite-element)模型來計算電磁場[9]。在巨粒子的模擬部份，通常選取粒子雲的大小與有限差分模型格網的大小相同，以粒子雲的中心位置來代表粒子的位置，當粒子移動時整個粒子雲一起跟著移動(如圖一所示)，並採用一階(線性)或高階加權(weighting)法將電荷加權至鄰近格網的格點上，因此這種方法被成為particle-in-cell(PIC)模型，其中採用一階加權的粒子式模型又稱為cloud-in-cell(CIC)模型[5]，採用高階加權可以降低電漿密度與電磁場計算的雜訊，但會大幅增加計算量，因此粒子模型通常採用一階加權；巨粒子的移動會產生電流，通過格網與格網介面之電流密度的計算亦是採用加權法算出每一個時間步驟(time step)通過格網介面的電荷，將此電荷除以通過格網介面的截面積與時間，即可得到格網介面上的電流密度，這樣的作法可以確保粒子運動仍滿足粒子連續方程式(continuity equation)的要求，最後此電流密度再加入Maxwell方程式，並根據電磁場的有限差分迭代以得到新的電場與磁場；在粒子所受電磁力以及推動粒子移動的計算方面，由於有限差分模型所得的電場與磁場都在格網的格線與介面上，因此在計算粒子所受到的電磁力時也必須採用加權法算出在粒子所在位置上的電場與磁場，同時描述荷電粒子運動的Newton-Lorentz方程式： 也採用有限差分時域模型，並根據粒子原來的位置、速度與所受的電磁力迭代計算出粒子的新位置與新速度。 圖一： PIC模型粒子移動示意圖 (time-domain)的計算中，粒子式模