1. Introduction
The low-Reynolds-number translational and rotational motions of colloidal particles in incompressible Newtonian fluids have attracted wide attention from researchers in the fields of chemical, biomedical, mechanical, civil, and environmental engineering. These motions are practical and fundamental in numerous processes such as agglomeration, sedimentation, centrifugation, microfluidics, aerosol technology, and rheology of suspensions. The theoretical investigation of this topic began with Stokes’ studies [
1,
2] on the creeping motions of hard spherical particles in unbounded viscous fluids. Masliyah et al. [
3] and Keh and Chou [
4] extended this analysis to the translation and rotation, respectively, of a soft sphere.
A soft particle of radius
has a hard core of radius
, covered by a permeable porous layer of thickness
. Polystyrene latices with surface layers [
5] and biological cells with surface attachments [
6] are examples of soft particles. To sterically stabilize colloidal dispersions, polymers are deliberately adsorbed by particles to form permeable layers [
7]. When the porous layers of soft spheres disappear, the particles revert to hard spheres. When the hard cores of soft spheres vanish, they become fully porous spheres (like permeable colloid flocs and polymer coils) [
8].
The hydrodynamic torque on a soft sphere of radius
(a hard core of radius
covered by a porous layer of thickness
) rotating with an angular velocity of
about a diameter in an unbounded fluid of viscosity
at low Reynolds numbers is [
4]
where,
and
is the penetration length (square root of permeability) of fluid flow within the surface layer of the soft particle (
and
are in opposite directions). In the limiting case
, Equation (1) degenerates to the Stokes result for a hard sphere of radius
.
In real situations of the rotation of particles, the surrounding fluid is bounded by solid walls [
9,
10,
11,
12]. Thus, it is necessary to know whether the proximity of boundary walls meaningfully affects particle rotation. The slow rotations of a hard sphere confined by adjacent boundaries, such as in a spherical cavity [
13,
14,
15,
16,
17], in a circular cylinder [
18,
19,
20], and near one or two planar walls [
13,
21,
22,
23], were analyzed. Alternatively, the low-Reynolds-number rotations of a soft or porous spherical particle in a spherical cavity [
4,
24,
25,
26,
27] and in a cylinder [
28] were also theoretically investigated. These studies show that the effect of boundaries on the rotation of particles can be very substantial and interesting.
In the general theories of stirred vessels and rotational viscometers for highly viscous liquids, it is important to understand the variation of torque as the confinement boundary approaches. The objective of this paper is to analyze the rotation of a soft colloidal sphere (having a porous layer of arbitrary thickness and permeability) about its diameter normal to one or two plane walls at an arbitrary position between them at a low Reynolds number. The fluid velocity was found by solving the Stokes and Brinkman equations using the boundary collocation method, and semianalytical results were obtained for the hydrodynamic torque acting on the particle for various values of the relevant parameters (the core-to-particle radius ratio, shielding parameter of the porous surface layer, and relative separation distances from the walls), with excellent convergence over the entire range.
2. Analysis
As illustrated in
Figure 1, we studied the creeping flow of a constant-property fluid around a soft spherical particle of radius
rotating steadily with a constant angular velocity
about a diameter perpendicular to two large planar walls whose distances from the particle center are
and
, respectively (
is taken without loss of generality), and
and
represent the spherical and cylindrical coordinate systems, respectively, originating from the particle center. The soft sphere comprises a permeable porous surface layer of thickness
. Thus, the radius of its hard core is
. The fluid velocity inside the porous layer is finite, while the external fluid far from the particle is at rest. The objective is to find the correction to Equation (1) for the particle rotation caused by the confining plane walls.
The creeping flow is governed by the Stokes and Brinkman equations, yielding
where
in cylindrical coordinates or
in spherical coordinates is the
(only nontrivial) component of the fluid velocity distribution, the continuity equation is satisfied, the dynamic pressure is constant everywhere,
is the penetration length (square root of permeability) of fluid flow within the surface layer, and
equals unity as
and zero otherwise. The boundary conditions require that the fluid is no slip at the hard-core surface and plane walls, and that both velocity and stress are continuous at the particle surface. Thus,
where
is the nontrivial shear stress at the particle surface. Equations (3)–(7) take the reference frame rotating with the particle.
The fluid velocity can be expressed in the form [
23,
27]
where
is the associated Legendre function of the first kind of order
and degree 1,
is the Bessel function of the first kind of order
,
and
are the modified Bessel functions of the first and second kinds, respectively, of order
,
,
,
,
, and
(all having the dimension of velocity) are the unknown functions and constants, respectively, to be determined. The parts of
involving
in the previous equations are separable solutions to Equation (3) in spherical coordinates that represent the disturbance generated by the particle and the part of
involving
in Equation (9) is a Fourier-Bessel integral solution to Equation (3) in cylindrical coordinates representing the disturbance produced by the planar walls. Note that Equation (9), which is a superposition of the general solutions in cylindrical and spherical coordinates due to the linearity of Equation (3), satisfies Equation (7) immediately.
Substitution of boundary condition (6) into Equation (9) leads to
where
The application of the Hankel transform on the variable
to Equations (10) and (11) yields
The solution of Equations (13) and (14) leads to
where
and
Substitution of Equations (15) and (16) back into Equation (9) results in
where
in which the integral can be calculated numerically.
The remaining boundary conditions to be fulfilled are those at the particle and hard core surfaces. Substituting Equations (8) and (20) into Equations (4) and (5) yields
where
The satisfaction of boundary conditions (22)–(24) at the inner and outer surfaces of the porous layer of the soft sphere requires solutions of the constants
,
, and
. The collocation technique [
29] permits these boundary conditions to be imposed at
points on the meridian semicircle of each surface and the infinite series in Equations (8) and (20) to be truncated after
terms, leading to
simultaneous linear algebraic equations. These algebraic equations can be numerically solved for sufficiently large
to result in the
constants
,
, and
. The details of the boundary collocation scheme were given in an early paper [
30] for a hard sphere translating perpendicular to two parallel plane walls.
The hydrodynamic torque acting on the soft sphere is [
27]
where
is the viscosity of the fluid and only the lowest-order constant
makes a contribution.
When the surface layer of the soft sphere vanishes, it degenerates to a hard sphere of radius , the constants , Equations (5), (8), (22), (24), and (25) are trivial, and just Equation (23) is required to be solved for the constants . When the hard core vanishes, the soft sphere reduces to an entirely porous sphere of radius , the constants , Equations (4) and (22) become trivial, and only Equations (23) and (24) need to be solved for the constants and .
3. Results and Discussion
The numerical solutions for the hydrodynamic torque
acting on a soft spherical particle rotating about its diameter perpendicular to two plane walls as a function of the ratio of the particle radius to the porous layer penetration length
, ratio of the core-to-particle radii
, particle-wall spacing parameter
, and relative position parameter
obtained from the boundary collocation method, are provided in
Table A1 and
Table A2 in
Appendix A for the distinct case of
(a fully porous sphere) and the general case, respectively. The torque
, given by Equation (1) for the soft sphere in the unbounded fluid is used to normalize
. The accuracy and convergence behavior of the collocation technique depends upon the relevant parameters. All the results obtained converge to at least six significant figures. For the most difficult case, the number of collocation points,
, is sufficiently large to achieve this convergence. These results are the same as those obtained for a hard sphere [
23] in the limiting case of
or
. Obviously,
is the limit
, regardless of other parameters. The wall effects on the rotational motion of the soft sphere can be significant.
The normalized torque
for a fully porous sphere rotating about its diameter perpendicular to two planar walls is plotted against the parameters
,
,
in
Figure 2,
Figure 3 and
Figure 4, respectively, over the entire range. For fixed values of
and
, as expected and shown in
Table A1 and
Figure 2a,b,
Figure 3a and
Figure 4b,
is a monotonically increasing function of the shielding parameter
(decreasing function of the permeability) for the fluid in the porous particle from unity (with
) at
to a larger finite value as
. When
is smaller than unity, the variation of
with
and
is weak (<1.4%).
of a porous sphere with low permeability (say,
) in general is close to that of a hard one (with
), though their difference can be noticeable when the particle is very close to a wall (
).
For the given values of
and
, as indicated in
Table A1 and
Figure 2b,
Figure 3a,b and
Figure 4a, the normalized torque
, acting on the confined porous sphere, is an increasing function of the particle-to-wall spacing parameter
from unity at
to a greater finite value at
(note that
is still finite even for the limit that the particle touches the plane walls), since the hydrodynamic hindrance caused by the plane walls is stronger when they locate closer to the particle. The dependence of
on
is robust when
is large but vanishes in the limit
. The supposition that the two-wall effect on the rotation of a particle can be viewed as a sum of single-wall effects will overestimate the hydrodynamic torque exerted on the particle. That is, the increase in
from unity for the two-equidistant-wall case
is less than twice that for the corresponding single-wall case
, which can be seen in
Table A1 and
Figure 2a,
Figure 3b and
Figure 4a,b.
For specified values of
and
, the normalized torque
of the porous sphere increases with an increase in the parameter
(denoting the relative position of the porous sphere between the walls) from a finite value at
(the case of a single wall) to a greater one at
(the case of two equally distant walls). Namely, the nearness of a second wall will enhance the torque acting on the particle close to the first wall. The variation of
with
can be significant when the value of
is large, though it disappears in the limits
and
. For a given value of
(the ratio of the particle diameter to the wall-to-wall distance), as revealed by the dashed curves in
Figure 4a, the torque is minimum when the particle locates in the middle between the two walls [
] and increases monotonically as the particle approaches either wall.
After understanding the hydrodynamic effect of two parallel plane walls on the axially symmetric rotation of a porous particle, we study the general case of that on a rotating soft particle. The results of the normalized torque
on a soft sphere rotating about its diameter perpendicular to two planar walls for different values of the core-to-particle radius ratio
, shielding parameter in the porous layer
, dimensionless spacing parameter
, and relative position parameter
are presented in
Figure 5,
Figure 6,
Figure 7 and
Figure 8 (together with
Table A2), respectively, over the entire ranges. Again,
increases as
increases from unity at
to a finite value at
and increases as
increases from a finite value at
to another at
, keeping the other parameters unchanged. Also,
is a monotonic increasing function of
from a constant (equal to zero for the entirely porous limit
) at
(the porous surface layer is completely permeable) to a great one as
(the surface layer is impermeable).
For fixed values of
,
, and
,
Table A2 and
Figure 5,
Figure 6,
Figure 7 and
Figure 8 show that the normalized torque
for the confined soft spherical particle undergoing rotation increases monotonically with an increase in the ratio of the core-to-particle radii
, where the limits
and
denote a porous sphere and an impermeable sphere, respectively. That is, if the porous layer is thicker for specified permeability, particle size, and separation from walls, the torque exerted on the particle will be less. All results for the soft spherical particle fall between the upper and lower bounds of
and
, respectively. For the circumstance where the surface layer has low to mediate permeability (e.g.,
), as shown in
Figure 5a and
Figure 8,
on the particle with
smaller than about 0.8 can be well approximated by that for a fully porous particle of the same size, permeability, and distances from walls. In this case, the relative motion of the fluid is barely felt by the hard core of the soft sphere, and its hindrance to the flow is negligible. However, this approximation is not valid for the porous layer with high permeability.
Recently, collocation results were obtained for the normalized hydrodynamic torque
of a soft sphere of radius
rotating about a diameter on the axis of a circular cylinder of radius
[
28]. Similar to the currently considered case of axisymmetric rotation of the particle perpendicular to two equidistant plane walls (i.e., at the center of a slit),
is a monotonically growing function of the shielding parameter
(from a value at
to a higher one at
), particle-wall spacing parameter
(from unity at
to a greater constant at
), and core-to-particle radius ratio
, holding other parameters constant. The particle in the circular cylinder bears much more torque than the particle in the slit does. This result manifests that the retardation to the particle rotation caused by the confinement walls is freed in both principal lateral directions of the slit, though only in an axial direction of the cylinder.
4. Concluding Remarks
The slow rotational motion of a soft spherical particle in a viscous fluid about its diameter perpendicular to one or two planar walls is semianalytically studied using the method of boundary collocation. Convergent numerical results for the torque exerted on the particle by the fluid were obtained for various values of the ratio of the particle radius to the flow penetration length of the porous layer , the ratio of the core-to-particle radii , particle-wall spacing parameter , and relative position parameter . The wall effect on the rotating soft particle can be significant. The normalized torque, , acting on the confined particle increases with an increase in from unity as (the particle is far from the walls) and remains finite even at the contact limit . This torque is smaller than that on a hard sphere (or soft one with larger or ), keeping the other parameters’ constant. For a given ratio of the particle diameter to the wall-to-wall distance , is minimal when the particle is midway between the two walls [] and increases as it locates closer to either wall [the value of decreases]. Experimental data of the normalized torque for the slow rotation of a soft particle near one or two plane walls would be needed to confirm the validity of our semianalytical collocation results at various ranges of , , , and .
Section three provides results for a resistance problem, where the hydrodynamic torque
on a particle rotating normal to two planar walls is considered for a given angular velocity
[equal to
according to Equation (1)]. In a mobility problem, the torque,
(equal to
), imposed to the particle is assumed and the boundary-corrected angular velocity,
, is considered. For a soft sphere rotating normal to two plane walls dealt with here, the normalized angular velocity
for the mobility problem is equal to
, as given in
Table A1 and
Table A2 and
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6,
Figure 7 and
Figure 8 for the resistance problem.