本发明属于高压开关设备绝缘设计与数值仿真领域，具体涉及一种GIS设备内部多金属颗粒流动特性的三维多场耦合仿真的方法。通过编写CFDEM的开源代码完成了同轴圆柱结构电场下颗粒流动的欧拉‑拉格朗日模型和几何建模，同时实现了FVM和DEM算法的数值耦合，计算在综合考虑了颗粒与颗粒、颗粒与流体、颗粒与电场、颗粒与电极和外壳的碰撞这一系列相互作用、以及颗粒与电极或者颗粒间发生碰撞时的荷电和电荷转移的物理过程的基础上完成，得到了颗粒流动特性的特征参数，如带电量、速度、位移和受到的作用力，再通过数形结合的后处理，最终观察到了多颗粒在直流和交流两种电压下的动态分布以及线形颗粒的翻滚状态。