等离子体粒子模拟.doc

等离子体粒子模拟 1 PIC/MCC模拟的理论基础 第一步，依据带电超粒子的速度与受力，由运动学方程计算下一时刻带电超粒子的位置与速度，进而计算带电超粒子的能量。 第二步，依据边界条件假设，处理运动到边界后的带电超粒子状态。 第三步，依据带电超粒子的能量，计算带电超粒子与中性粒子是否发生碰撞以及发生碰撞后的状态。 第四步，依据一定的权重化方法，将带电超粒子的电荷量累加到相应的网格点上，得到等离子体电荷密度分布。 第五步，依据等离子体密度分布和相应的边界条件，解泊松方程，得到电场分布。 第六步，带电超粒子的位置与电场分布，计算出带电超粒子的受力。至此完成一次循环，进入下一次循环。 2 MCC方法 MCC方法是产生一系列[0,1]之间均匀分布的随机数，通过比较随机数与由碰撞截面表示的碰撞几率，来得出粒子是否碰撞。如果碰撞，由碰撞几率决定发生何种碰撞，再结合能量和动量守恒等定律来确定碰撞后粒子的新状态。 MCC过程中，考虑了氩气放电中的主要碰撞过程，包括e-n:弹性、激发、电离碰撞；i-n:弹性、电荷交换碰撞。相应的碰撞截面计算公式列于表2-1中。 碰撞类型 碰撞反应 碰撞截面 电子弹性碰撞 e+Ar→e+Ar 电子激发碰撞 e+Ar→e+Ar* 电子电离碰撞 e+Ar→e++e 离子电荷交换碰撞 +Ar→Ar+ 离子弹性碰撞 +Ar→+Ar 注：E为；粒子能量，单位为eV 其中计算电子电离碰撞截面积的系数改成了 用MCC方法描述粒子间的碰撞，每个时间步长内带电粒子发生碰撞的概率为： 其中N为中性粒子的密度；和分别为带电粒子的能量、速率和总碰撞截面。 为电子与中性粒子的弹性碰撞截面； 为电子与中性粒子的激发碰撞截面； 为电子与中性粒子的电离碰撞截面； 为离子与中性粒子的弹性碰撞截面； 为离子与中性粒子的电荷交换碰撞截面。 所采用的各种碰撞截面均依赖于带电粒子的能量。 电子部分： 如果，则电子与中性粒子发生弹性碰撞； 如果，则电子与中性粒子发生激发碰撞； 如果，则电子与中性粒子发生电离碰撞。 离子部分： 如果，则电子与中性粒子发生弹性碰撞； 如果，则电子与中性粒子发生电荷交换碰撞。 发生碰撞后粒子的新状态依据具体碰撞类型决定。 电子与中性粒子的弹性碰撞 e+Ar→e+Ar 因为（M为中性粒子的质量），所以假设碰撞后电子的速率或能量保持不变，只是方向改变，电子的散射角可由下式求得： 式中为碰撞前电子的能量，R为[0,1]之间均匀分布的随机数。另外，电子散射方向与以入射方向为轴所成的极角（即散射角向角），通常看成均匀分布于[0,]范围内，可表示为： 其中，为[0,1]之间的均匀随机数。 散射后电子的运动法向与z轴的夹角 其中为碰撞前电子运动方向与z轴的夹角。 由相对论的动能关系：（没用这种方法，算出来的速度差三个数量级） 其中，e为电子电量，T为碰撞后粒子的动能（单位eV），为电子散射过程中的能量损耗，W为相对论下粒子的总能量，为粒子的静止质量，P为动量大小。 假设碰撞前电子各方向上的动量大小分别为（平行于图2-2中x方向），则有： 设分别为的归一化参量，则可导出散射后电子的三个方向归一化动量如下： 可得碰撞后电子动量的各分量分别为： 电子与中性粒子的激发碰撞 e+Ar→e+Ar* 碰撞后，电子失去该反应所需阈值能量，即电子的能量为 其中为电子激发阈值。分别计算和，再结合弹性碰撞处理方法即可求解电子碰撞激发后的速度。 碰撞电离过程 e+Ar→e++e 电离过程通常产生一组电子离子对，其能量守恒方程可写为： 其中，分别为散射、新产生以及入射电子能量，为新产生的离子能量，为中性粒子能量，为电离所需阈值能量。 由于电子与中性粒子质量比很小，假设入射电子的动量远小于中性粒子，即可看成：中性粒子脱落一个电子，然后自身转化为离子，碰撞后产生的离子能量等于碰撞前中性粒子的能量。因此： 散射电子与新电子能量将由入射电子损耗电离能后的能量随机分配或平均分配，分别计算和，进而得到各自的动量。 新电子的运动方向、方位角由随机数决定，运动方向与轴向夹角由随机数决定。 被散射电子的散射角和方位角由随机数决定。 新产生的离子通常给定一个相对较低的能量值，且其速度遵循Maxwellian分布。 离子与中性粒子碰撞 （1）粒子与中性粒子的弹性碰撞 +Ar→+Ar 散射后离子的能量为 其中，分别为散射和入射离子的能量。为能量损耗因子，其表达式为： 当和M大小相近时，有： 交换电荷碰撞过程 +Ar→Ar+ 碰撞后新生离子的能量由中性气体温度的Maxwellian分布取样。碰撞后新生离子的速度方向随机决定。 模拟初始化 2.1 一维模型 初始化包括空间位置初始化和速度初始化。假设模拟系统的长度为L，模拟的带电粒子数量为N，为第i个粒子的位置。 蒙特卡罗法，产生一个0到1之间的