Eligiusz Wajnryb

Friction and mobility of many spheres in Stokes flow

An efficient scheme is presented for the numerical calculation of hydrodynamic interactions of many spheres in Stokes flow. The spheres may have various sizes, and are freely moving or arranged in rigid arrays. Both the friction and mobility matrix are found from the solution of a set of coupled equations. The Stokesian dynamics of many spheres and the friction and mobility tensors of polymers and proteins may be calculated accurately at a modest expense of computer memory and time. The transport coefficients of suspensions can be evaluated by use of periodic boundary conditions.

Friction and mobility for colloidal spheres in Stokes flow near a boundary: The multipole method and applications

We obtain the many-body hydrodynamic friction and mobility matrices describing the motion in a fluid of N hard-spheres with stick boundary conditions in the presence of a planar hard wall or free surface using (1) a multipole expansion of the hydrodynamic force densities induced on the spheres and (2) an image representation to account for the fluid boundary. The coupled multipole equations may be truncated at any order to give positive definite approximations to the exact friction and mobility matrices. An extension of the Bossis–Brady lubrication correction to the friction matrix is also discussed and included. The resulting method for computing the mobility matrix may be used for the Stokesian or Brownian dynamics simulation of N spheres subject to interparticle and external forces and imposed shear flow. We illustrate the method by performing Stokesian dynamics simulation of particles near a hard wall. The simulations exhibit the rapid convergence of the multipole truncation scheme including lubricati...

Streaming potential studies of colloid, polyelectrolyte and protein deposition

Recent developments in the electrokinetic determination of particle, protein and polyelectrolyte monolayers at solid/electrolyte interfaces, are reviewed. Illustrative theoretical results characterizing particle transport to interfaces are presented, especially analytical formulae for the limiting flux under various deposition regimes and expressions for diffusion coefficients of various particle shapes. Then, blocking effects appearing for higher surface coverage of particles are characterized in terms of the random sequential adsorption model. These theoretical predictions are used for interpretation of experimental results obtained for colloid particles and proteins under convection and diffusion transport conditions. The kinetics of particle deposition and the structure of monolayers are analyzed quantitatively in terms of the generalized random sequential adsorption (RSA) model, considering the coupling of the bulk and surface transport steps. Experimental results are also discussed, showing the dependence of the jamming coverage of monolayers on the ionic strength of particle suspensions. In the next section, theoretical and experimental results pertaining to electrokinetics of particle covered surfaces are presented. Theoretical models are discussed, enabling a quantitative evaluation of the streaming current and the streaming potential as a function of particle coverage and their surface properties (zeta potential). Experimental data related to electrokinetic characteristics of particle monolayers, mostly streaming potential measurements, are presented and interpreted in terms of the above theoretical approaches. These results, obtained for model systems of monodisperse colloid particles are used as reference data for discussion of experiments performed for polyelectrolyte and protein covered surfaces. The utility of the electrokinetic measurements for a precise, in situ determination of particle and protein monolayers at various interfaces is pointed out.

Lubrication corrections for three-particle contribution to short-time self-diffusion coefficients in colloidal dispersions

It is shown that the standard treatment of lubrication effects in many-particle hydrodynamic interactions leads to divergent three-particle contributions to the short-time translational self-diffusion coefficient. To resolve the problem the improved method to account for lubrication is proposed. The translational and rotational self-diffusion coefficients of the Brownian semidilute suspension are then evaluated up to terms of the second order in volume fraction.

Hydrodynamic interactions of spherical particles in suspensions confined between two planar walls

Hydrodynamic interactions in a suspension of spherical particles confined between two parallel planar walls are studied under creeping-flow conditions. The many-particle friction matrix in this system is evaluated using our novel numerical algorithm based on transformations between Cartesian and spherical representations of Stokes flow. The cartesian representation is used to describe the interaction of the fluid with the walls and the spherical representation is used to describe the interaction with the particles. The transformations between these two representations are given in a closed form, which allows us to evaluate the coefficients in linear equations for the induced-force multipoles on particle surfaces. the friction matrix is obtained from these equations, supplemented with the superposition lubrication corrections. we have used our algorithm to evaluate the friction matrix for a single sphere, a pair of spheres, and for linear chains of spheres. The friction matrix exhibits a crossover from a quasi-two-dimensional behaviour (for systems with small wall separation H) to THE three-dimensional behaviour (when the distance H is much larger than the interparticle distance L). the crossover is especially pronounced for a long chain moving in the direction normal to its orientation and parallel to the walls. in this configuration, a large pressure build-up occurs in front of the chain for small values of the gapwidth H, which results in a large hydrodynamic friction force. a standard wall superposition approximation does not capture this behaviour.

Swapping trajectories : a new wall-induced cross-streamline particle migration mechanism in a dilute suspension of spheres

Binary encounters between spherical particles in shear flow are studied for a system bounded by a single planar wall or two parallel planar walls under creeping flow conditions. We show that wall proximity gives rise to a new class of binary trajectories resulting in cross-streamline migration of the particles. The spheres on these new trajectories do not pass each other (as they would in free space) but instead they swap their cross-streamline positions. To determine the significance of the wall-induced particle migration, we have evaluated the hydrodynamic self-diffusion coefficient associated with a sequence of uncorrelated particle displacements due to binary particle encounters. The results of our calculations quantitatively agree with the experimental value obtained by Zarraga & Leighton ( Phys. Fluids , vol. 14, 2002, p. 2194) for the self-diffusivity in a dilute suspension of spheres undergoing shear flow in a Couette device. We thus show that the wall-induced cross-streamline particle migration is the source of the anomalously large self-diffusivity revealed by their experiments.

Many-particle hydrodynamic interactions in parallel-wall geometry: Cartesian-representation method

This paper describes the results of our theoretical and numerical studies of hydrodynamic interactions in a suspension of spherical particles confined between two parallel planar walls, under creeping-flow conditions. We propose a novel algorithm for accurate evaluation of the many-particle friction matrix in this system—no such algorithm has been available so far.

Hydrodynamic interactions and collision efficiencies of spherical drops covered with an incompressible surfactant film

A theory is developed for the hydrodynamic interactions of surfactant-covered spherical drops in creeping flows. The surfactant is insoluble, and flow-induced changes of surfactant concentration are small, i.e. the film of adsorbed surfactant is incompressible. For a single surfactant-covered drop in an arbitrary incident flow, the Stokes equations are solved using a decomposition of the flow into surface-solenoidal and surface-irrotational components on concentric spherical surfaces. The surface-solenoidal component is unaffected by surfactant; the surface-irrotational component satisfies a slip-stick boundary condition with slip proportional to the surfactant diffusivity. Pair hydrodynamic interactions of surfactant-covered bubbles are computed from the one-particle solution using a multiple-scattering expansion. Two terms in a lubrication expansion are derived for axisymmetric near-contact motion. The pair mobility functions are used to compute collision efficiencies for equal-size surfactant-covered bubbles in linear flows and in Brownian motion. An asymptotic analysis is presented for weak surfactant diffusion and weak van der Waals attraction

Far-field approximation for hydrodynamic interactions in parallel-wall geometry

A complete analysis is presented for the far-field creeping flow produced by a multipolar force distribution in a fluid confined between two parallel planar walls. We show that at distances larger than several wall separations the flow field assumes the Hele-Shaw form, i.e. it is parallel to the walls and varies quadratically in the transverse direction. The associated pressure field is a two-dimensional harmonic function that is characterized by the same multipolar number m as the original force multipole. Using these results we derive asymptotic expressions for the Greens matrix that represents Stokes flow in the wall-bounded fluid in terms of a multipolar spherical basis. This Greens matrix plays a central role in our recently proposed algorithm [Physica A 356 (2005) 294] for evaluating many-body hydrodynamic interactions in a suspension of spherical particles in the parallel-wall geometry. Implementation of our asymptotic expressions in this algorithm increases its efficiency substantially because the numerically expensive evaluation of the exact matrix elements is needed only for the neighboring particles. Our asymptotic analysis will also be useful in developing hydrodynamic algorithms for wall-bounded periodic systems and implementing acceleration methods by using corresponding results for the two-dimensional scalar potential.

Generalization of the Rotne-Prager-Yamakawa mobility and shear disturbance tensors

The Rotne–Prager–Yamakawa approximation is one of the most commonly used methods of including hydrodynamic interactions in modelling of colloidal suspensions and polymer solutions. The two main merits of this approximation are that it includes all long-range terms (i.e. decaying as