We study an ad hoc extension of the Lattice-Boltzmann method that allows the simulation of non-Newtonian fluids described by generalized Newtonian models. We extensively test the accuracy of the method for the case of shear-thinning and shear-thickening truncated power-law fluids in the parallel plate geometry, and show that the relative error compared to analytical solutions decays approximately linear with the lattice resolution. Finally, we also tested the method in the reentrant-flow geometry, in which the shear-rate is no-longer a scalar and the presence of two singular points requires high accuracy in order to obtain satisfactory resolution in the local stress near these points. In this geometry, we also found excellent agreement with the solutions obtained by standard finite-element methods, and the agreement improves with higher lattice resolution.

The dynamics of macroscopically homogeneous sheared suspensions of neutrally buoyant, non-Brownian spheres is investigated in the limit of vanishingly small Reynolds numbers using Stokesian dynamics. We show that the complex dynamics of sheared suspensions can be characterized as a chaotic motion in phase space and determine the dependence of the largest Lyapunov exponent on the volume fraction ϕ. We also offer evidence that the chaotic motion is responsible for the loss of memory in the evolution of the system and demonstrate this loss of correlation in phase space. The loss of memory at the microscopic level of individual particles is also shown in terms of the autocorrelation functions for the two transverse velocity components. Moreover, a negative correlation in the transverse particle velocities is seen to exist at the lower concentrations, an effect which we explain on the basis of the dynamics of two isolated spheres undergoing simple shear. In addition, we calculate the probability distribution function of the transverse velocity fluctuations and observe, with increasing ϕ, a transition from exponential to Gaussian distributions.The simulations include a non-hydrodynamic repulsive interaction between the spheres which qualitatively models the effects of surface roughness and other irreversible effects, such as residual Brownian displacements, that become particularly important whenever pairs of spheres are nearly touching. We investigate, for very dilute suspensions, the effects of such a non-hydrodynamic interparticle force on the scaling of the particle tracer diffusion coefficients Dy and Dz, respectively, along and normal to the plane of shear, and show that, when this force is very short-ranged, both are proportional to ϕ2 as ϕ → 0. In contrast, when the range of the non-hydrodynamic interaction is increased, we observe a crossover in the dependence of Dy on ϕ, from ϕ2 to ϕ as ϕ → 0. We also estimate that a similar crossover exists for Dz but at a value of ϕ one order of magnitude lower than that which we were able to reach in our simulations.

We performed macroscopic experiments on the motion of a sphere through an array of obstacles that highlight the deterministic nature of the lateral displacements that lead to particle separation in microfluidic systems. The motion of the spheres is irreversible and displays directional locking. The locking directions can be predicted with a single parameter that distinguishes between reversible and irreversible particle-obstacle collisions. These results stress the need to incorporate irreversible interactions to predict the movement of a non-Brownian sphere passing through a periodic array.