Dynamics of electrorheological suspensions subjected to spatially nonuniform electric fields

A numerical method based on the distributed Lag range multiplier method (DLM) is developed for the direct simulation of electrorheological (ER) liquids subjected to spatially nonuniform electric field. The flow inside particle boundaries is constrained to be rigid body motion by the distributed Lag range multiplier method and the electrostatic forces acting on the particles are obtained using the point-dipole approximation. The numerical scheme is verified by performing a convergence study which shows that the results are independent of mesh and time step sizes. The dynamical behavior of ER suspensions subjected to nonuniform electric field depends on the solids fraction, the ratio of the domain size and particle radius, and four additional dimensionless parameters which respectively determine the importance of inertia, viscous, electrostatic particle-particle interaction and dielectrophoretic forces. For inertia less flows a parameter defined by the ratio of the dielectrophoretic and viscous forces, determines the time duration in which the particles collect near either the local maximums or local minimums of the electric field magnitude, depending on the sign of the real part of the Clausius-Mossotti factor. In a channel subjected to a given nonuniform electric field, when the applied pressure gradient is smaller than a critical value, the flow assists in the collection of particles at the electrodes, but when the pressure gradient is above this critical value the particles are swept away by the flow.