介绍了NEPTUNE软件采用的一些并行计算方法：采用“块-网格片”二层并行区域分解方法，使计算规模能够扩展到上千个处理器核。基于复杂几何特征采用自适应技术并行生成结构网格，在原有规则区域的基础上剔除无效网格，大幅降低了存储量和并行执行时间。在经典的Boris和SOR迭代方法基础上，采用红黑排序和几何约束，提出了非规则区域上的Poisson方程并行求解方法。采用这些方法后，当使用NEPTUNE软件模拟MILO器件时，可在1 024个处理器核上获得51.8%的并行效率。

介绍了自主编制的3维全电磁粒子模拟大规模并行程序NEPTUNE的基本情况。该程序具备对多种典型高功率微波源器件的3维模拟能力，可以在数百乃至上千个CPU上稳定运行。该程序使用时域有限差分(FDTD)方法更新计算电磁场，采用Buneman-Boris算法更新粒子运动状态，运用质点网格法(PIC)处理粒子与电磁场的耦合关系，最后利用Boris方法求解泊松方程对电场散度进行修正，以确保计算精度。该程序初步具备复杂几何结构建模能力，可以对典型高功率微波器件中常见的一些复杂结构，如任意边界形状的轴对称几何体、正交投影面几何体，慢波结构、耦合孔洞、金属线和曲面薄膜等进行几何建模。该程序将理想导体边界、外加波边界、粒子发射与吸收边界及完全匹配层边界等物理边界应用于几何边界上，实现了数值计算的封闭求解。最后以算例的形式，介绍了使用NEPTUNE程序对磁绝缘线振荡器、相对论返波管、虚阴极振荡器及相对论速调管等典型高功率微波源器件进行的模拟计算情况，验证了模拟计算结果的可靠性，同时给出了并行效率的分布情况。