Application of the Multi-Particle Collision Dynamics Method to a Suspension of Magnetic Spherical Particles

Abstract

In order to apply the multi-particle collision dynamics (MPCD) method to a magnetic particle suspension, we have elucidated the dependence of the translational and rotational Brownian motion of magnetic particles on the MPCD parameters that characterize the MPCD simulation method. We here consider a two-dimensional system composed of magnetic spherical particles in thermodynamic equilibrium. The diffuse reflection model has been employed for treating the interactions between fluid and magnetic particles. In the diffuse reflection model, the interactions between fluid and magnetic particles are transferred into the translational motion more strongly than into the rotational motion of magnetic particles. The employment of relatively small simulation time steps gives rise to a satisfactory level of the translational Brownian motion. The activation level of the Brownian motion is almost independent of both the size of the unit collision cell and the number of fluid particles per cell. Larger values of the maximum rotation angle induce stronger translational and rotational Brownian motion, but in the present magnetic particle suspension the range between around π/4 and π/2 seems to be reasonable. We may conclude that the MPCD method with the simple diffuse reflection model is a feasible simulation technique as the first approximation for analyzing the behavior of magnetic particles in a suspension. If more accurate solutions regarding the aggregate structures of magnetic particles are required, the introduction of the scaling coefficient regarding the interactions between fluid and magnetic particles can yield more accurate and physically reasonable aggregate structures in both a qualitative and quantitative meanings.