等离子体模拟改进电磁粒子网格算法

摘要

为了消除传统电磁粒子网格算法在交错推进电场和磁场、粒子动量和位置时存在的半时间步长差,推导出一种电磁场和粒子非交错推进的改进电磁粒子网格算法.该算法通过引入一个中间过渡时刻,对麦克斯韦方程组和牛顿-洛仑兹方程的差分格式进行了变换,使算法实现了同时将所有物理量从n时刻推进到n+1时刻.该算法能得到时间上同步的参数,有利于计算结果的数据处理以及在模拟过程中对电磁场和粒子进行实时诊断.应用改进后的电磁粒子网格算法对无箔二极管实例进行了模拟,得到了与实验基本一致的结果.