《3维全电磁粒子模拟大规模并行程序NEPTUNE Massively parallel code named NEPTUNE for 3D fully electromagnetic and PIC simulations》.pdfVIP

第23巷第6期 强 激 光 与 粒 子 束 V01．23，No．6 2011年6月 HIGHP()WERI。ASERANDPARTlCI，EBEAMS Jun．，2011 文章编号： 100卜4322(2011)06—1607一09 董 烨， 陈 军， 杨温渊， 赵 强， 夏 芳， 肖 丽， 马 彦， 廖 丽， 孙会芳， 董志伟， 周海京， 陈 虹， 莫则尧 (北京应用物理与计算数学研究所。北京lO0088) 摘 要： 介绍了自主编制的3维全电磁粒子模拟大规模并行程序NEPTUNE的基本情况。该程序具备 对多种典型高功率微波源器件的3维模拟能力，可以在数百乃至上千个CPU上稳定运行。该程序使用时域有 粒子与电磁场的耦合关系．最后利用Boris方法求解泊松方程对电场散度进行修正。以确保计算精度。该程序 初步具备复杂几何结构建模能力。可以对典型高功率微波器件中常见的一些复杂结构。如任意边界形状的轴对 称几何体、正交投影面几何体．慢波结构、耦合孔洞、金属线和曲面薄膜等进行几何建模。该程序将理想导体边 界、外加波边界、粒子发射与吸收边界及完全匹配层边界等物理边界应用于几何边界上．实现了数值计算的封 闭求解。最后以算例的形式。介绍了使用NEPTUNE程序对磁绝缘线振荡器、相对论返波管、虚阴极振荡器及 相对论速调管等典型高功率微波源器件进行的模拟计算情况，验证了模拟计算结果的可靠性，同时给出了并行 效率的分布情况。 关键词： 高功率微波； 3维全电磁粒子模拟； 大规模并行程序l NEPTUNE程序 中图分类号： TN25 文献标志码： A doi：10．3788／HPLP1607 作用是非常强的非线性物理过程，无法运用纯理论方法精确解析求解；另一方面，HPM器件实验费用较为昂 贵且耗时耗力。基于以上原因及计算机技术的迅速发展和应用，利用计算机粒子模拟技术来研究和设计 HPM器件成为一个新的发展方向。随着HPM器件研究的深入及计算机并行技术的发展。人们有需求也有能 德国CST公司的MAFIA以及俄罗斯的KARAT等程序具备模拟HPM器件的能力“．5]。 为适应我国HPM源发展研究的需求，我们研制了3维全电磁粒子模拟大规模并行程序NEPTUNE[17‘181。 用Boris方法[1 2]求解泊松方程对电场散度进行修正，以确保计算精度)；程序模拟计算结果采用可视化工具 JavisL91进行显示和输出。 边界条件设计与实现请参考文献[15—16]，程序的整体并行架构设计请参照参考文献[17—18]。 l NEPTUNE程序核心数值算法 1．1 基本物理量的定义与说明 NEPTUNE程序采用笛卡尔坐标系，计算网格为标准的均匀长方体网格．数值计算中用到的六类基本物 理量，包括：电场矢量、磁场矢量、电荷密度、电流密度矢量、粒子位置及粒子动量矢量，这6类基本物理量基于 计算网格的定义方式如下所述。 ·收稿日期：2010—09一06·修订日期：2011一Ol-11 基金项目：国家高技术发展计划项Ij·中闷r程物理研究院科学技术发展基金项日(200980402046) 作者简介：董烨(1981一)．男．硕士．助理研究员．从事高功率微波技术方面的研究-dongye0682@sina．com． 万方数据 1608 强 激 光 与 粒 子 束 第23卷 如图1～2所示，电场矢量在笛卡尔坐标系中化为3个分量，E，E，和E。，分别定义在计算网格的棱心位 置，时刻定义在整数时刻点；磁场矢量化为3个分量，H，，H，和H：，分别定义在计算网格的面心位置，时刻定 义在半整数时刻点；电荷密度定义在计算网格的节点位置，时刻定义在整数时刻点；电流密度矢量化为3个分 量，．，，，J，与．，。，分别定义在计算网格的棱心位置，时刻定义在半整数时刻点；粒子位置坐标为(z，y，z)，可位 于计算网格的任意位置，时刻定义在整数时刻点；粒子动量矢量可以化为3个分量以，p，，户：，可以位于计算网 格的任意位置，时刻定义在半整数时刻点。