Bead milling is a shear force-dominant method used for particle processing; and in designing and optimizing a bead mill process, it is desirable to apply numerical technologies. The discrete element method (DEM) coupled with computational fluid dynamics (CFD) is widely used to simulate bead motion in a bead mill. However, the existing Eulerian-Lagrangian method has a problem regarding the computation of bead motion in an arbitrarily shaped bead mill, while the Lagrangian-Lagrangian method involves high computational costs because of the calculation of fluid particles. In order to solve these problems, the advanced DEM-CFD method has been proposed. In this method, collision detection for the solid phase and boundary conditions for the fluid phase in a moving boundary system are modeled by signed distance functions and the immersed boundary method. In this study, the advanced DEM-CFD method was applied to bead milling. The simulation method was validated by comparing the simulated and experimental motion of beads. Since shear force is dominant in bead milling, the numerical study focused on friction parameters. The results suggest that the simulations are in good agreement with experiments, and the coefficients of friction are very important for bead mill simulation.