The postulated core disruptive accidents (CDAs) are regarded as particular difficulties in the safety analysis of liquid-metal fast reactors (LMFRs). in CDAs, the motions and interactions of solid particles, such as refrozen fuels, disrupted pellets, etc., not only dominate fundamental behaviors of multiphase flows, but also drastically influence the process of CDAs. The fast reactor safety analysis code, SiMMER-iii, which is a 2D, multi-velocity-field, multiphase, multicomponent, Eulerian, fluid dynamics code coupled with a fuel-pin model and a space- and energy-dependent neutron kinetics model, was successfully applied to a series of CDA assessments. However, strong interactions among solid particles as well as particle characteristics in multiphase flows with rich solid particles were not taken into consideration for fluid-dynamics models of SiMMER-iii. in this article, a hybrid method for multiphase flows analysis is developed by coupling the discrete element method (DEM) with the multi-fluid model of SiMMER-iii. in the coupling algorithm, motions of liquid and gas phases are solved by a time-factorization (time-splitting) method. For the solid phases, contacts among particles and interactions with fluid phases are considered through DEM. Numerical simulations of dam-break behavior with rich solid particles show reasonable agreements with corresponding experimental results. it is expected that SiMMER-iii coupled with DEM could provide a promising and useful computational tool for complicated multiphase-flow phenomena with high concentration of solid particles.