Magnetic separation is an important method in the processing process, and its essence is the targeted dispersion of the mineral processing slurry pulp in the magnetic field space. The slurry is a complex multiphase fluid system with continuous phase carrying a large number of discrete phase particles, in which the magnetic particles agglomerate, migrate, and disperse under the dominance of magnetic force. In this process, there is nonlinear and unstable dynamic coupling between the continuous phase (liquid) and the discrete phase (solid particles) and between the discrete phases. In this paper, a dynamic cyclic multi-dipole magnetic moment algorithm with a higher calculation accuracy is innovatively proposed to calculate the magnetic interaction force between particles. Moreover, the P-E magnetization model suitable for a two-dimensional uniform magnetic field is further improved and optimized to make it applicable to a three-dimensional gradient magnetic field. Finally, based on the coupling of the Finite Element Method (FEM), Computational Fluid Dynamics (CFD), and Discrete Element Method (DEM), a dynamic coupling model capable of accurately simulating the magnetic separation process is developed. This model can be used to study the separation behavior of particles under a multiphase flow and multi-force field and to explore the motion behavior of magnetic particles.

High-gradient magnetic separation (HGMS) is a conventional and effective method for processing weak magnetic materials. A multi-field dynamic coupling simulation method integrating the Finite Element Method (FEM), Computational Fluid Dynamics (CFD), and the Discrete Element Method (DEM) was employed to investigate the separation behavior in HGMS. The dynamic deposition process of magnetic particles under the interactions of magnetic fields, fluid flow fields, and particle–particle forces was simulated using a two-way fluid–solid coupling algorithm based on the FEM-CFD-DEM coupling approach. Experimental results demonstrated that the particle deposition profiles predicted by the double-wire medium model were in good agreement with the measured data. The research findings indicated that the separation process could be divided into three distinct stages—the adsorption stage, the closure stage, and the clogging stage—each characterized by unique dynamic behaviors and pressure-drop evolution patterns. Additionally, the effects of key parameters such as the feeding velocity and medium filling ratio on the separation process were analyzed, providing theoretical foundations and technical support for the optimization of HGMS processes and the enhancement of separation efficiency.