基于Matlab的二维PIC/MCC模型的实现

Realization of Two-Dimensional PIC/MCC Model based on Matlab

为了实现限定条件下气体放电的二维模拟,基于Matlab通过编程手段开发模拟程序,并用此模拟程序对直流低压条件下氩气的放电过程进行验证。对电子速度进行麦克斯韦初始化,并进行结果验证。提取50个时间步长内的电子数变化趋势,验证了氩气放电过程中电子崩的发生;并且通过100个时间步长结束时的电子位置分布,可以观察到阳极鞘层的产生;根据解泊松方程至200步和300步结果的对比研究,验证在本文的模拟条件下200步解泊松方程是足够的;对电子速度均方根处理,得到温度分布示意图,验证了直流氩气放电时阳极板附近温度高的实验现象。至此,基于Matlab开发的二维粒子/蒙特卡洛碰撞的耦合(PIC/MCC)程序,能够进行气体放电模拟。

In order to achieve the two-dimensional simulation of gas discharge under limited conditions, a simulation program based on Matlab has been developed in the paper, and utilized to verify the discharge process of argon under direct current and low pressure. Electronic speed was initialized through Maxwell equations and the results were verified. The occurrence of electron avalanche in the process of argon discharge was verified by extracting the change trend of the number of electrons in the 50 time steps, and the production of anode sheath was observed through the distribution of electron position at the end of the 100 time step. According to the solution of Poisson＇s equation to the 200 step and 300 step, the solution of Poisson＇s equation in the 200 step was sufficient under the simulated conditions in this paper. The diagram of temperature distribution given by the root mean square treatments of the electron velocity successfully verified the experimental phenomena that the temperature near the anode plate was higher in the DC argon discharge. Thus, the two-dimensional PIC/MCC program based on Matlab can satisfactorily simulate the gas discharge.