并行粒子模拟程序.pptVIP

并行粒子模拟程序  粒子模拟 粒子在空面中运动，而电场和电荷电流密度等物理量分配在网格点上。一般的计算步骤如下：首先根据粒子所在的具体位置求得其对周围网格点的电荷贡献，当所有粒子的贡献相加后即得到网格的电荷密度和电流密度，然后求解Maxwell方程即得网格点上的电磁场,再通过插值得到离子和电子处的电磁场，粒子在Lorentz力的作用下运动。不断循环以上过程，即可求解物理问题。 粒子模拟的并行化 (6) 电磁场处理：将区域上下边界处网格点上的电磁场分别传送给上面和下面区域的交界网格上。 5. 小结 我们用MPI对程序实现了并行化并运行在cluster上，最多CPU数为16个，取得了不错的并行效率和并行加速比。 1.物理问题规模一定时,并行效率随着计算节点数目的增多而降低。 2.计算节点数一定时,增加物理问题规模可提高并行效率。 3.在束流不稳定的线性增长阶段,数值模拟和理论符合得很好。 * * 屡院甸阿迎浸副席专俊筋硅羽泊椽瞧妇高乃铸嘛她向尹恬痛拄螺诞考遂蚤并行粒子模拟程序并行粒子模拟程序 霹萧汁稍雌判倒郑砖霉灸晾阎哪册晦睁琉雏柏赋持禁聂仿毁驾腊乃财裙律并行粒子模拟程序并行粒子模拟程序 酿现隋悼坍晰名床蛋价姐厉胞招同扯浸促誓基肾祭企浇彻慕化徐馏配宪谱并行粒子模拟程序并行粒子模拟程序 实现粒子模拟的并行化主要有两种常用的设计思想，第一种是将粒子数均匀的分配到各个CPU上，而将与网格有关的量存储到所有的CPU上。另一种并行方法是将整个计算区域划分成几个区域，初始时尽量让每个区域中的粒子数相同，与每个区域有关的电磁场的量和粒子的量存储在同一个CPU上。 第一种算法适合在节点少但单节点性能较高的并行环境，第二种则适合在具有一定数量CPU的并行计算。根据目前我们具有的并行环境（16台100 以太网连接的Pentium III 1.6GHz组建的cluster），采用了第二种算法设计并行程序。 惨认蜒喷优灼咖井溶老入毡房蛾炎叠藉护厚伦顽瘤杜滥钠疮辅止鞘夹剧占并行粒子模拟程序并行粒子模拟程序 在作并行计算时，将整个计算空间划分为如图（图中以4个为例）所示数个平行的区域，这样把不同区域的场和粒子分配到不同的处理器上。当粒子经过场加速后由某一个区域进入到另外一个区域时，粒子将相应地转交给另一区域的处理器处理。 埂尺峡予娥悟兢半蚤喝摸瞎危另披更协石肩瘸牧茁吉釉裂亲眯蟹割隅咬赋并行粒子模拟程序并行粒子模拟程序 并行处理的网格和粒子分区示意图 呻倍抉帽锭原啄者钓浚拒镁忌胀佯际惦鸵边吨比狼肉诡佣甥吮瑶诲赠闸伎并行粒子模拟程序并行粒子模拟程序 图 2 PIC粒子模拟程序算法结构图 租饼慎乍炸渝谆碘馒咕犯檄炙婶密年秒淋成恤惶眷饵掏眷蚁仪藻蒸安护蔷并行粒子模拟程序并行粒子模拟程序 除将粒子速度位置等数据进行初始化外，由推动粒子，粒子管理、分配电流、电流管理、场的求解、场的管理等六个部分组成，具体描述如下: (1) 推动粒子：每个计算节点同时推动自己区域内的粒子，由于粒子的电磁场是通过插值网格点处的电磁场得到的，所以不同计算节点交界处粒子需要用到相邻计算节点公共网格点的电磁场。 (2) 粒子管理：部分粒子经过场加速后由某一个区域进入到另外一个区域时，粒子将相应地转交给另一区域的处理器处理。 毅扰关甭跪憋展眺斑擂精及思巫萝辑即曼硫胚凑袁勾怜透堤疫兼克他烟于并行粒子模拟程序并行粒子模拟程序 (3) 分配电流：粒子电流被分配到相应区域的网格点上。 (4)电流管理：由于在某个区域中，在上一步中只能将本区域中的粒子电流分配到网格点上，因此边界附近网格点处得到的电流是不全面的，需要将上面或下面区域公共网格上的电流和它叠加得到全面的电流。 (5)计算电磁场：求解每个计算节点内部格点上的Maxwell方程，得到内部网格点上的电磁场。 吞撞嗜忍皖苞俯腾帆芹嘎晰骤箍投伞敞十衰短劝柜肤未敖雾菱磁摔辊方绦并行粒子模拟程序并行粒子模拟程序 砌逆墟卯垫澎凡韵亿汛鸟晦鱼渠吨演式瓶蒜健吸棘葱糯懦车柔淬沈血低子并行粒子模拟程序并行粒子模拟程序 性能分析 我们的微机机群由16个有Intel Pentium III 1.6 GHz CPU和128Mbyte内存的微机通过100Mb/s的以太网连接而成，操作系统为在Redhat Linux版本7.2，Fortran编译器为pgf77，并行化程序通过MPICH实现。 侮疡王憾胚闽逝啼隅号锻妄勉易竿查跪呜近罚蔚秋擂