89b Multiphase Modeling of Electro- and Magnetorheological Suspensions

David Kittipoomwong1, Yevgeny Yurkovetsky2, and Daniel J. Klingenberg1. (1) University of Wisconsin, 1415 Engineering Drive, Madison, WI 53706, (2) Levich Institute, City College of CUNY, 140th St. & Convent Ave., New York, NY 10031

Application of an electric or magnetic field to suspensions of polarizable or magnetizable particles, respectively, produces an increase of the apparent viscosity of the suspensions in shear flow. The viscosity increase is rapid, occurring over 1-100 milliseconds. For large applied field strengths, this initial, rapid transient response is often followed by a much slower transient increase in the apparent suspension viscosity. The slow transient response has been associated with the field-induced formation of lamellar structures.

A two-fluid continuum model has been employed to investigate the formation of the lamellar structures in these materials. In this model, particle migration is related to the particle contribution to the normal stresses induced by the applied field and by the shear flow. Using naive approximations for these stresses, the model successfully predicts the formation of lamellar structures, and presumably the onset of a slow transient rheological response, for field strengths above a "critical" value. However, the naive approximations result in several other predictions that do not agree with experimental observations. Most notably, the predicted shear rheology is predicted to be insensitive to the applied field strength and the volume fraction dependence of the critical field strength does not agree with that observed experimentally.

We employ particle-level simulations to determine more accurate representations of the particle contributions to the stress, as functions of the applied field strength, shear rate, particle concentration, and shear history. Results for the normal stress components provide improved information to describe both the onset of lamella formation, as well as the long-time structure. Results for the deviatoric stress components provide information about the conventional field-induced rheology as well as its transient behavior.