Application of the Brownian dynamics method to a rod-like hematite particle dispersion

Abstract

We have investigated various problems that arise in applying the Brownian dynamics method to a suspension composed of rod-like hematite particles, which have a magnetic moment normal to the particle axis direction. The accuracy and the deviation of simulation results from theoretical solutions have been discussed by comparison with the theoretical solutions that were obtained by solving the basic equations of the orientational distribution function. The characteristics of the negative viscosity are not observed to be dependent on a time interval unless a sufficiently short time interval is used. The present simulation results can satisfactorily reproduce the qualitative characteristics of the negative magneto-rheological effect that was predicted by the previous theoretical investigation. Good quantitative agreement is obtained in the range of small-applied magnetic fields, but agreement is not significantly good in the large magnetic field region. The deviation of the negative viscosity from the theoretical prediction cannot be improved by using a more accurate numerical algorithm such as moving from Euler to second-order or fourth-order Runge–Kutta. The results of the orientational distribution can well reproduce the characteristic features that the distribution has a gradual shape with low linear-like peak, which is in significant contrast to the sharp single-peak distribution of a ferromagnetic rod-like particle dispersion. The present orientational distributions are in significantly good agreement with those of the theoretical prediction in regard to the position and the height of a peak and the general shape of the overall profile. Good agreement of the present magnetisation curves with the theoretical prediction verifies that the spin rotational Brownian motion is activated at a physically reasonable level in the present simulations.