We report an extension of the smoothed profile method (SPM)[Y. Nakayama, K. Kim, and R. Yamamoto, Eur. Phys. J. E {bf 26}, 361(2008)], a direct numerical simulation method for calculating the complex modulus of the dispersion of particles, in which w
e introduce a temporally oscillatory external force into the system. The validity of the method was examined by evaluating the storage $G(omega)$ and loss $G(omega)$ moduli of a system composed of identical spherical particles dispersed in an incompressible Newtonian host fluid at volume fractions of $Phi=0$, 0.41, and 0.51. The moduli were evaluated at several frequencies of shear flow; the shear flow used here has a zigzag profile, as is consistent with the usual periodic boundary conditions.
