In the safety analysis of sodium-cooled fast reactors, numerical simulations of various thermal-hydraulic phenomena with multicomponent and multiphase flows in core disruptive accidents (CDAs) are regarded as particularly difficult. In the material relocation phase of CDAs, core debris settle down on a core support structure and/or an in-vessel retention device and form a debris bed. The bed’s shape is crucial for the subsequent relocation of the molten core and heat removal capability as well as re-criticality. In this study, a hybrid numerical simulation method, coupling the multi-fluid model of the three-dimensional fast reactor safety analysis code SIMMER-IV with the discrete element method (DEM), was applied to analyze the sedimentation and bed formation behaviors of core debris. Three-dimensional simulations were performed and compared with results obtained in a series of particle sedimentation experiments. The present simulation predicts the sedimentation behavior of mixed particles with different properties as well as homogeneous particles. The simulation results on bed shapes and particle distribution in the bed agree well with experimental measurements. They demonstrate the practicality of the present hybrid method to solid particle sedimentation and bed formation behaviors of mixed as well as homogeneous particles.