蒙特卡罗碰撞模型在尤普软件中的应用研究.pdfVIP

蒙特卡罗碰撞模型在尤普软件中的应用 王洪广，刘纯亮，李永东 西安交通大学电子物理与器件研究所，710049 摘要：本文研究了蒙特卡罗碰撞模型在等离子体模拟程序——尤普(UNPIC)软件 中的应用，将模型扩展到高能量粒子碰撞范围，通过对充中性气体的相对论返波管进 行了数值模拟，对模型进行了验证。 关键词：蒙特卡罗碰撞等离子体粒子模拟尤普软件相对论返波管 1 已I吉 ．L· ’，I口 随着气体放电等离子体的广泛应用，迫切需要对气体放电现象作更深入的理解。 CarloCollision) 程序——尤普软件(uNIPIC)中加入蒙特卡罗碰撞(MCC，Monte 模型，来模拟研究气体放电现象。 一般的蒙特卡罗碰撞是采用碰撞时间随机的方法，而PIC方法中，粒子推进和场 推进的时间步长是固定的。因此在PIC方法中，采用了在一个时间步长内随机决定粒 子间是否发生碰撞的方法来实现蒙特卡罗碰撞【l】。加入蒙特卡罗碰撞的等离子体粒子 模拟(PIC—MCC)的流程图如图1所示。本文主要介绍尤普软件中的MCC模型，最 后通过充中性气体相对论返波管进行了模拟验证。 粒子推进H粒孟鎏界 计算粒子 蒙特卡罗 受力 碰撞 场推进H燃 圆1 pTr、．h斤rr谛辉圆 蒙特卡罗碰撞模型 2．1粒子碰撞过程的处理 碰撞模型中的粒子分为两类：耙粒子和源粒子。假设第i个源粒子和耙粒子之间 可发生N种碰撞，设第J种(．『∈(1，N))碰撞的碰撞截面为仃，(￡；)。那么，总碰撞截 面为所有类型的碰撞截面之和： 仃r(￡f)=仃I(￡f)+仃2(￡f)+…+仃Ⅳ(岛) (1) 则第i个源粒子与耙粒子发生碰撞(包括所有N种碰撞)的几率——总碰撞几率 只=1-exp(q(Ei)，z，(xi)viAt) (2) 其中：x，、v，为第i个源粒子的位置和速度，n,(xi)为第i个源粒子位置处的耙 粒子密度，出是PIC推进时间步长。 为了避免在每次推进中对每个源粒子进行大量的计算，向PIC—MCC方法中引入 了空碰撞【2】的概念，即什么都不发生的碰撞。引入空碰撞之后，总碰撞几率可以被设 定为一个常数P咖ll’每一次推进时只需随机抽取部分源粒子进行碰撞处理。随机抽取 480 源粒子发生碰撞。适当选取缸，使Pnullo．0l，从而使计算量大大减小。然后再由产 生的随机数RE(R2∈[0，1])确定碰撞类型。 各种碰撞类型碰撞后各粒子的位置、速度和方向处理： (1)电子与中性原子的弹性碰撞【3] 碰撞后电子的运动方向是随机的，散射角Z为和子午面角≯为： 驴：2万足，cosz：竺坐燮 (3) 。 一 E v·：v／1一T2memi(1-cosz) (4) ‘V (m。+mi)2 式中：弛、mi分别为电子、离子质量，v为入射电子速度。 (2)电子与中性原子的激发碰撞 按照先激发再弹性碰撞的方法来处理，设中性原子的激发能为＆，那么经过激 发碰撞后，入射电子的动能E’=E—ee，速度大小改变为√2El／m，即为弹性碰撞的新 入射速度，然后再用弹性碰撞来求取碰撞后的速度。 (3)电子与中性原子的电离碰撞 按照先电离再弹性碰撞的方法来处理对于新产生的电子，其初始位置可近似为 入射电子的位置，初始动能为