本发明公开了一种适用于粗颗粒流的气粒混合流动数值计算方法,针对单个粗颗粒占据多个网格单元的情况,通过建立适用的气体-颗粒耦合曳力模型,将计算流体力学与离散单元法相耦合;该方法以欧拉-拉格朗日法为框架,将空隙率及动量交换源项引入N-S方程,通过有限体积法及二阶迎风格式将求解方程离散,并采用传统的SIMPLE算法求解气相各个物理量;同时应用离散单元法求解颗粒的受力及运动;气体-颗粒间的相互作用通过动量交换源项及耦合曳力模型进行求解。该方法跟踪研究每个颗粒的受力及运动,能够充分给出粗颗粒流的微观动力学特性,同时克服了一般CFD-DEM耦合方法求解粗颗粒流精度低误差大的缺点,填补了CFD-DEM耦合模型在精确计算粗颗粒流上的空缺。