Abstract The motion behavior of nonmetallic inclusions at the interface of molten steel and slag fundamentally affects the removal of inclusions. Therefore, from an analysis of forces, this study constructed a mathematical model of inclusion movement. Compared with other models that only consider the forces acting on nonmetallic inclusions at the interface, the proposed model considers not only cases in the inclusions which enter the slag interior and rebound into the molten steel, but also the effect of fluid flow containing the inclusions with different Re numbers on the drag force. The application of this established model has not taken Reynolds number of fluid flow into consideration. The model can predict the motion of inclusions at the interface and in nearby areas and provide a curve of inclusion displacement vs time. The mathematical model was verified with a physical model, with the curve of displacement vs time obtained from physical experiment being consistent with the calculated curve. The preliminary calculation results show that inclusions having liquid film at their surfaces are rebounded into the steel when they have size within a certain range but enter the slag phase directly when they are beyond that size range.