A boundary element model for particle dynamics simulation of magnetic fluids in magnetic fields is employed to calculate the field distributions and interactions in a particulate system.Compared with the finite element method,the boundary element method calculates sharper and more variable boundary area of particles.In the area other than boundary areas,the two methods show good agreement.With increasing node numbers,the magnetic force obtained by the boundary element method closes to the virtual force obtained by the finite element method.