拟颗粒模拟自1996年提出以来在颗粒流体系统的离散模拟上得到了很多重要的应用，但串行计算限制了模拟的规模，为了将拟颗粒模型应用于更大规模体系的模拟，并行化成为必然选择。　　 本文在一维区域分割气固拟颗粒模拟的并行程序基础上实现了二维区域分割的并行算法。在本算法中，粒子在每一个时间步内独立运动，运动后移出进程区域边界的粒子经Shift通信模式传递到相应的目标进程。对于粒子碰撞的计算，首先进行进程区域内部粒子碰撞的计算；然后经过边界信息交换，相邻进程共享边界粒子信息后，再进行边界部分粒子碰撞的计算。针对基于硬球作用模型的拟颗粒模拟，本算法保证了在时驱算法下边界粒子在相邻的各进程内碰撞顺序的一致性。算法的实现使并行模拟的拓扑分区不再受流场区域分割数目的限制，适应于更大舰模体系的模拟。　　 同时，本文针对上述并行算法进行了效率分析和优化，通过对程序中主要函数的耗时分析，发现进程分配的粒子数不均匀和粒子链表结构复杂是影响程序运行效率的重要原因。在程序的优化上，选择IntelR icc编译器且参数选择“-fast”时，可以使所测算例的程序耗时减少约20％。同时，通过在程序中使用动态负载均衡和数据链表重排，可以使程序运行时间在相同计算规模上减少了30％-50％。最后提出了进一步优化的方向。　　 作为应用，本文分别模拟了含11550个固体粒子、2.43×106个气体拟颗粒和含102400个固体粒子、1.8×107个气体拟颗粒的流态化体系，得到了与文献相一致的结果。最后利用所得模拟结果对能量最小多尺度模型(EMMS)稳定性条件成立的空间尺度进行了初步探索。