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 behavior (for systems with small wall separation H) to the three-dimensional behavior (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 buildup 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 behavior.

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.Our approach involves expanding the fluid velocity field into spherical and Cartesian fundamental sets of Stokes flows. The interaction of the fluid with the particles is described using the spherical basis fields; the flow scattered with the walls is expressed in terms of the Cartesian fundamental solutions. At the core of our method are transformation relations between the spherical and Cartesian basis sets. These transformations allow us to describe the flow field in a system that involves both the walls and particles.We used our accurate numerical results to test the single-wall superposition approximation for the hydrodynamic friction matrix. The approximation yields fair results for quantities dominated by single particle contributions, but it fails to describe collective phenomena, such as a large transverse resistance coefficient for linear arrays of spheres.

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 GreenÕs matrix that represents Stokes flow in the wall-bounded fluid in terms of a multipolar spherical basis. This GreenÕs 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.

Using a recently developed image representation for Stokes flow in a half-space bounded by a planar wall ͓Cichocki and Jones, Physica A 258, 273 ͑1998͔͒, the image system is constructed for the flow field produced by a force multipole in the space bounded by two parallel walls. The image singularities are expressed in terms of products of double-reflection matrices, and the expansion is simplified using symmetries of the double-reflection operation. Our analysis yields recurrence relations for the strengths of the image multipoles. The relations are solved explicitly, and a complete image system is obtained for an arbitrary source-force multipole. Applications of our image representation for evaluating the hydrodynamic friction and mobility matrices of particles interacting with two parallel planar walls are indicated.