Next Article in Journal
Unique Fiber Morphologies from Emulsion Electrospinning—A Case Study of Poly(ε-caprolactone) and Its Applications
Previous Article in Journal
On Classification of Water-in-Oil and Oil-in-Water Droplet Generation Regimes in Flow-Focusing Microfluidic Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Slow Rotation of a Soft Colloidal Sphere Normal to Two Plane Walls

Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
*
Author to whom correspondence should be addressed.
Colloids Interfaces 2023, 7(1), 18; https://doi.org/10.3390/colloids7010018
Submission received: 6 February 2023 / Revised: 20 February 2023 / Accepted: 22 February 2023 / Published: 23 February 2023

Abstract

:
The creeping flow of a viscous fluid around a soft colloidal sphere rotating about a diameter normal to two planar walls at an arbitrary position between them is theoretically investigated in the steady limit of small Reynolds numbers. The fluid velocity outside the particle consists of the general solutions of the Stokes equation in circular cylindrical and spherical coordinates, while the fluid velocity inside the porous surface layer of the particle is expressed by the general solution of the Brinkman equation in spherical coordinates. The boundary conditions are implemented first on the planar walls by means of the Hankel transforms and then at the particle and hard-core surfaces by a collocation technique. The torque exerted on the particle by the fluid is calculated as a function of the ratio of the core-to-particle radii, ratio of the particle radius to the flow penetration length of the porous layer, and relative particle-to-wall spacings over the entire range. The wall effect on the rotating soft particle can be significant. The hydrodynamic torque exerted on the confined soft sphere increases as the relative particle-to-wall spacings decrease and stays finite even when the soft sphere contacts the plane walls. It is smaller than the torque on a hard sphere (or soft one with a reduced thickness or penetration length of the porous layer), holding the other parameters constant. For a given relative wall-to-wall spacing, this torque is minimal when the particle is situated midway between the walls and rises as it locates closer to either wall.

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 b has a hard core of radius a , covered by a permeable porous layer of thickness b a . 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 b (a hard core of radius a covered by a porous layer of thickness b a ) rotating with an angular velocity of Ω about a diameter in an unbounded fluid of viscosity η at low Reynolds numbers is [4]
T 0 = 8 π η λ 2 b Ω R λ a cosh ( λ b λ a ) + sinh ( λ b λ a ) ,
where,
R = ( λ 3 a b 2 3 λ b + 3 λ a ) cosh ( λ b λ a ) + ( λ 2 b 2 3 λ 2 a b + 3 ) sinh ( λ b λ a ) ,
and 1 / λ is the penetration length (square root of permeability) of fluid flow within the surface layer of the soft particle ( T 0 and Ω are in opposite directions). In the limiting case λ b , Equation (1) degenerates to the Stokes result for a hard sphere of radius b .
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 b rotating steadily with a constant angular velocity Ω about a diameter perpendicular to two large planar walls whose distances from the particle center are c and d , respectively ( c d is taken without loss of generality), and ( r , θ , φ ) and ( ρ , φ , z ) represent the spherical and cylindrical coordinate systems, respectively, originating from the particle center. The soft sphere comprises a permeable porous surface layer of thickness b a . Thus, the radius of its hard core is a . 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
[ 2 ρ 2 h ( r ) λ 2 ] v ϕ = 1 r 2 r ( r 2 v ϕ r ) + 1 r 2 θ [ 1 sin θ θ ( v ϕ sin θ ) ] h ( r ) λ 2 v ϕ = 0 ,
where v ϕ ( ρ , z ) in cylindrical coordinates or v ϕ ( r , θ ) in spherical coordinates is the ϕ (only nontrivial) component of the fluid velocity distribution, the continuity equation is satisfied, the dynamic pressure is constant everywhere, λ 1 is the penetration length (square root of permeability) of fluid flow within the surface layer, and h ( r ) equals unity as a < r < b 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,
r = a : v ϕ = 0 ,
r = b : v ϕ   and   τ r ϕ   are   continuous ,
z = c ,   d : v ϕ = Ω ρ ,
ρ   and c < z < d : v ϕ = Ω ρ ,
where τ r ϕ 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]
v ϕ = n = 1 ( λ r ) 1 / 2 [ C n I n + 1 / 2 ( λ r ) + D n K n + 1 / 2 ( λ r ) ] P n 1 ( cos θ )   if   a r b ,
v ϕ = Ω ρ + λ 2 0 [ X ( ω ) e ω z + Y ( ω ) e ω z ] ω J 1 ( ω ρ ) d ω + n = 1 A n ( λ r ) n 1 P n 1 ( cos θ )
if   r b   and c z d ,
where P n 1 is the associated Legendre function of the first kind of order n and degree 1, J n is the Bessel function of the first kind of order n , I n and K n are the modified Bessel functions of the first and second kinds, respectively, of order n , X ( ω ) , Y ( ω ) , A n , C n , and D n (all having the dimension of velocity) are the unknown functions and constants, respectively, to be determined. The parts of v ϕ involving P n 1 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 v ϕ involving J n 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
0 [ X ( ω ) e ω c + Y ( ω ) e ω c ] ω J 1 ( ω ρ ) d ω = λ 2 n = 1 A n α n ( ρ , c ) ,
0 [ X ( ω ) e ω d + Y ( ω ) e ω d ] ω J 1 ( ω ρ ) d ω = λ 2 n = 1 A n α n ( ρ , d ) ,
where
α n ( ρ , z ) = [ λ 2 ( ρ 2 + z 2 ) ] ( n + 1 ) / 2 P n 1 [ z ( ρ 2 + z 2 ) 1 / 2 ] .
The application of the Hankel transform on the variable ρ to Equations (10) and (11) yields
X ( ω ) e ω c + Y ( ω ) e ω c = λ 2 n = 1 A n 0 α n ( ρ , c ) ρ J 1 ( ω ρ ) d ρ ,
X ( ω ) e ω d + Y ( ω ) e ω d = λ 2 n = 1 A n 0 α n ( ρ , d ) ρ J 1 ( ω ρ ) d ρ .
The solution of Equations (13) and (14) leads to
X ( ω ) = n = 1 A n X n ( ω ) ,
Y ( ω ) = n = 1 A n Y n ( ω ) ,
where
X n ( ω ) = e ω c [ B n ( ω , c ) + e ω ( c + d ) B n ( ω , d ) ] 1 + e 2 ω ( c + d ) ,
Y n ( ω ) = e ω d [ e ω ( c + d ) B n ( ω , c ) B n ( ω , d ) ] 1 + e 2 ω ( c + d ) ,
and
B n ω , z = e ω z n 1 ! ( ω z λ z ) n 1 .
Substitution of Equations (15) and (16) back into Equation (9) results in
v ϕ = Ω ρ + n = 1 A n γ n ( r , θ ) if   r b   and c z d ,
where
γ n ( r , θ ) = λ 2 0 [ X n ( ω ) e ω r cos θ + Y n ( ω ) e ω r cos θ ] ω J 1 ( ω r sin θ ) d ω + ( λ r ) n 1 P n 1 ( cos θ ) ,
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
n = 1 [ C n I n + 1 / 2 ( λ a ) + D n K n + 1 / 2 ( λ a ) ]   ( λ a ) 1 / 2 P n 1 ( cos θ ) = 0 ,
n = 1 { [ C n I n + 1 / 2 ( λ b ) + D n K n + 1 / 2 ( λ b ) ] ( λ b ) 1 / 2 P n 1 ( cos θ ) A n γ n ( b , θ ) } = Ω b sin θ ,
n = 1 { [ C n { λ b I n 1 / 2 ( λ b ) + λ b I n + 3 / 2 ( λ b ) 3 I n + 1 / 2 ( λ b ) } D n { λ b K n 1 / 2 ( λ b )
+ λ b K n + 3 / 2 ( λ b ) + 3 K n + 1 / 2 ( λ b ) } ] ( λ b ) 1 / 2 P n 1 ( cos θ ) 2 A n γ n * ( b , θ ) } = 0 ,
where
γ n * ( r , θ ) = r 2 r [ γ n ( r , θ ) r ] = λ 2 r 0 { [ X n ( ω ) e ω r cos θ Y n ( ω ) e ω r cos θ ] J 1 ( ω r sin θ ) cos θ
[ X n ( ω ) e ω r cos θ + Y n ( ω ) e ω r cos θ ] J 2 ( ω r sin θ ) sin θ } ω 2 d ω ( n + 2 ) ( λ r ) n 1 P n 1 ( cos θ )
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 A n , C n , and D n . The collocation technique [29] permits these boundary conditions to be imposed at N points on the meridian semicircle of each surface and the infinite series in Equations (8) and (20) to be truncated after N terms, leading to 3 N simultaneous linear algebraic equations. These algebraic equations can be numerically solved for sufficiently large N to result in the 3 N constants A n , C n , and D n . 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]
T = 8 π η λ 2 A 1 ,
where η is the viscosity of the fluid and only the lowest-order constant A 1 makes a contribution.
When the surface layer of the soft sphere vanishes, it degenerates to a hard sphere of radius b = a , the constants C n = D n = 0 , Equations (5), (8), (22), (24), and (25) are trivial, and just Equation (23) is required to be solved for the N constants A n . When the hard core vanishes, the soft sphere reduces to an entirely porous sphere of radius b , the constants D n = 0 , Equations (4) and (22) become trivial, and only Equations (23) and (24) need to be solved for the 2 N constants A n and C n .

3. Results and Discussion

The numerical solutions for the hydrodynamic torque T 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 λ b , ratio of the core-to-particle radii a / b , particle-wall spacing parameter b / c , and relative position parameter c / ( c + d ) obtained from the boundary collocation method, are provided in Table A1 and Table A2 in Appendix A for the distinct case of a = 0 (a fully porous sphere) and the general case, respectively. The torque T 0 , given by Equation (1) for the soft sphere in the unbounded fluid is used to normalize T . 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, N = 46 , 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 λ b or a = b . Obviously, T / T 0 = 1 is the limit b / c = 0 , regardless of other parameters. The wall effects on the rotational motion of the soft sphere can be significant.
The normalized torque T / T 0 for a fully porous sphere rotating about its diameter perpendicular to two planar walls is plotted against the parameters λ b , b / c , c / ( c + d ) in Figure 2, Figure 3 and Figure 4, respectively, over the entire range. For fixed values of b / c and c / ( c + d ) , as expected and shown in Table A1 and Figure 2a,b, Figure 3a and Figure 4b, T / T 0 is a monotonically increasing function of the shielding parameter λ b (decreasing function of the permeability) for the fluid in the porous particle from unity (with T = T 0 = 0 ) at λ b = 0 to a larger finite value as λ b . When λ b is smaller than unity, the variation of T / T 0 with b / c and c / ( c + d ) is weak (<1.4%). T / T 0 of a porous sphere with low permeability (say, λ b > 100 ) in general is close to that of a hard one (with λ b ), though their difference can be noticeable when the particle is very close to a wall ( b / c 1 ).
For the given values of λ b and c / ( c + d ) , as indicated in Table A1 and Figure 2b, Figure 3a,b and Figure 4a, the normalized torque T / T 0 , acting on the confined porous sphere, is an increasing function of the particle-to-wall spacing parameter b / c from unity at b / c = 0 to a greater finite value at b / c = 1 (note that T / T 0 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 T / T 0 on b / c is robust when λ b is large but vanishes in the limit λ b = 0 . 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 T / T 0 from unity for the two-equidistant-wall case c / ( c + d ) = 1 / 2 is less than twice that for the corresponding single-wall case c / ( c + d ) = 0 , which can be seen in Table A1 and Figure 2a, Figure 3b and Figure 4a,b.
For specified values of λ b and b / c , the normalized torque T / T 0 of the porous sphere increases with an increase in the parameter c / ( c + d ) (denoting the relative position of the porous sphere between the walls) from a finite value at c / ( c + d ) = 0 (the case of a single wall) to a greater one at c / ( c + d ) = 1 / 2 (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 T / T 0 with c / ( c + d ) can be significant when the value of λ b is large, though it disappears in the limits λ b = 0 and b / c = 0 . For a given value of 2 b / ( c + d ) (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 [ c / ( c + d ) = 1 / 2 ] 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 T / T 0 on a soft sphere rotating about its diameter perpendicular to two planar walls for different values of the core-to-particle radius ratio a / b , shielding parameter in the porous layer λ b , dimensionless spacing parameter b / c , and relative position parameter c / ( c + d ) are presented in Figure 5, Figure 6, Figure 7 and Figure 8 (together with Table A2), respectively, over the entire ranges. Again, T / T 0 increases as b / c increases from unity at b / c = 0 to a finite value at b / c = 1 and increases as c / ( c + d ) increases from a finite value at c / ( c + d ) = 0 to another at c / ( c + d ) = 1 / 2 , keeping the other parameters unchanged. Also, T / T 0 is a monotonic increasing function of λ b from a constant (equal to zero for the entirely porous limit a / b = 0 ) at λ b = 0 (the porous surface layer is completely permeable) to a great one as λ b (the surface layer is impermeable).
For fixed values of λ b , b / c , and c / ( c + d ) , Table A2 and Figure 5, Figure 6, Figure 7 and Figure 8 show that the normalized torque T / T 0 for the confined soft spherical particle undergoing rotation increases monotonically with an increase in the ratio of the core-to-particle radii a / b , where the limits a / b = 0 and a / b = 1 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 a / b = 1 and a / b = 0 , respectively. For the circumstance where the surface layer has low to mediate permeability (e.g., λ b 10 ), as shown in Figure 5a and Figure 8, T / T 0 on the particle with a / b 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 T / T 0 of a soft sphere of radius b rotating about a diameter on the axis of a circular cylinder of radius c [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), T / T 0 is a monotonically growing function of the shielding parameter λ b (from a value at λ b = 0 to a higher one at λ b ), particle-wall spacing parameter b / c (from unity at b / c = 0 to a greater constant at b / c = 1 ), and core-to-particle radius ratio a / b , 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 λ b , the ratio of the core-to-particle radii a / b , particle-wall spacing parameter b / c , and relative position parameter c / ( c + d ) . The wall effect on the rotating soft particle can be significant. The normalized torque, T / T 0 , acting on the confined particle increases with an increase in b / c from unity as b / c = 0 (the particle is far from the walls) and remains finite even at the contact limit b / c = 1 . This torque is smaller than that on a hard sphere (or soft one with larger a / b or λ b ), keeping the other parameters’ constant. For a given ratio of the particle diameter to the wall-to-wall distance 2 b / ( c + d ) , T / T 0 is minimal when the particle is midway between the two walls [ c / ( c + d ) = 1 / 2 ] and increases as it locates closer to either wall [the value of c / ( c + d ) 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 λ b , a / b , b / c , and c / ( c + d ) .
Section three provides results for a resistance problem, where the hydrodynamic torque T on a particle rotating normal to two planar walls is considered for a given angular velocity Ω [equal to ( T 0 / 8 π η λ 2 b R ) [ λ a cosh ( λ b λ a ) + sinh ( λ b λ a ) ] according to Equation (1)]. In a mobility problem, the torque, T (equal to 8 π η λ 2 b Ω 0 R / [ λ a cosh ( λ b λ a ) + sinh ( λ b λ a ) ] ), 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 Ω / Ω 0 for the mobility problem is equal to ( T / T 0 ) 1 , 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.

Author Contributions

Conceptualization, H.J.K.; methodology, H.J.K. and C.L.C.; investigation, H.J.K. and C.L.C.; writing—original draft preparation, H.J.K. and C.L.C.; writing—review and editing, H.J.K.; supervision, H.J.K.; funding acquisition, H.J.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology, Taiwan (Republic of China) grant number MOST 110-2221-E-002-017-MY3.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

The collocation solutions for the normalized torque T / T 0 acting on a soft sphere 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 λ b , ratio of the core-to-particle radii a / b , particle-wall spacing parameter b / c , and relative position parameter c / ( c + d ) are presented in Table A1 and Table A2 for the limiting case of a = 0 (a fully porous sphere) and general case, respectively.
Table A1. Normalized torque T / T 0 for a porous sphere ( a = 0 ) rotating about a diameter perpendicular to two parallel planar walls.
Table A1. Normalized torque T / T 0 for a porous sphere ( a = 0 ) rotating about a diameter perpendicular to two parallel planar walls.
c / ( c + d ) b / c T / T 0
λ b = 1 λ b = 10 λ b = 100 λ b = 600
00.11.000011.000091.000121.00012
0.31.000211.002471.003291.00337
0.51.000951.011561.015441.01585
0.61.001651.020211.027141.02786
0.71.002621.032701.044321.04554
0.81.003921.050311.069371.07144
0.91.005591.075261.107661.11143
0.951.006591.091831.136411.14210
0.991.007461.108041.170651.18075
0.9951.007581.110311.176471.18807
0.9991.007671.112181.181651.19516
1/40.11.000011.000091.000121.00013
0.31.000211.002511.003331.00342
0.51.000971.011731.015671.01607
0.61.001671.020501.027531.02826
0.71.002661.033171.044951.04619
0.81.003971.051031.070341.07243
0.91.005671.076311.109071.11289
0.951.006681.093071.138111.14386
0.991.007571.109461.172611.18280
0.9951.007691.111751.178491.19017
0.9991.007781.113631.183711.19734
1/20.11.000011.000161.000221.00022
0.31.000371.004461.005941.00609
0.51.001721.021011.028131.02886
0.61.002971.036901.049721.05106
0.71.004731.060071.081801.08410
0.81.007081.093081.129181.13312
0.91.010111.140401.202631.20992
0.951.011921.172121.258471.26959
0.991.013521.203371.325581.34554
0.9951.013731.207751.337101.36004
0.9991.013891.211351.347331.37420
Table A2. Normalized torque T / T 0 for a soft sphere with λ b = 1 rotating about a diameter perpendicular to two parallel planar walls.
Table A2. Normalized torque T / T 0 for a soft sphere with λ b = 1 rotating about a diameter perpendicular to two parallel planar walls.
c / ( c + d ) b / c T / T 0
a / b = 0.5 a / b = 0.8 a / b = 0.95
00.11.000021.000061.00011
0.31.000551.001751.00290
0.51.002531.008161.01362
0.61.004381.014201.02387
0.71.006981.022821.03880
0.81.010461.034691.06020
0.91.014981.050791.09159
0.951.017691.060881.11354
0.991.020091.070191.13645
0.9951.020411.071451.13982
0.9991.020661.072461.14263
1/40.11.000021.000071.00011
0.31.000551.001771.00294
0.51.002561.008271.01381
0.61.004441.014401.02421
0.71.007081.023151.03935
0.81.010611.035191.06104
0.91.015201.051511.09283
0.951.017941.061731.11501
0.991.020381.071171.13814
0.9951.020701.072441.14154
0.9991.020961.073471.14437
1/20.11.000041.000121.00019
0.31.000981.003161.00525
0.51.004571.014791.02477
0.61.007921.025841.04365
0.71.012641.041721.07145
0.81.018991.063781.11175
0.91.027271.093971.17167
0.951.032251.113041.21400
0.991.036671.130741.25853
0.9951.037251.133131.26511
0.9991.037721.135071.27059

References

  1. Stokes, G.G. On the theories of the internal friction of fluids in motion and of the equilibrium and motion of elastic solids. Trans. Camb. Phil. Soc. 1845, 8, 287–319. [Google Scholar]
  2. Stokes, G.G. On the effect of the internal friction of fluids on the motion of pendulums. Trans. Camb. Phil. Soc. 1851, 9, 8–106. [Google Scholar]
  3. Masliyah, J.H.; Neale, G.; Malysa, K.; van de Ven, T.G.M. Creeping flow over a composite sphere: Solid core with porous shell. Chem. Eng. Sci. 1987, 42, 245–253. [Google Scholar] [CrossRef]
  4. Keh, H.J.; Chou, J. Creeping motion of a composite sphere in a concentric spherical cavity. Chem. Eng. Sci. 2004, 59, 407–415. [Google Scholar] [CrossRef]
  5. Anderson, J.L.; Solomentsev, Y. Hydrodynamic effects of surface layer on colloidal particles. Chem. Eng. Commun. 1996, 148–150, 291–314. [Google Scholar] [CrossRef]
  6. Wunderlich, R.W. The effects of surface structure on the electrophoretic mobilities of large particles. J. Colloid Interface Sci. 1982, 88, 385–397. [Google Scholar] [CrossRef]
  7. Napper, D.H. Polymeric Stabilization of Colloidal Dispersions; Academic Press: London, UK, 1983. [Google Scholar]
  8. Neale, G.; Epstein, N.; Nader, W. Creeping flow relative to permeable spheres. Chem. Eng. Sci. 1973, 28, 1865–1874. [Google Scholar] [CrossRef]
  9. Malysa, K.; van de Ven, T.G.M. Rotational and translational motion of a sphere parallel to a wall. Int. J. Multiph. Flow 1986, 12, 459–468. [Google Scholar] [CrossRef]
  10. Liu, Q.; Prosperetti, A. Wall effects on a rotating sphere. J. Fluid Mech. 2010, 657, 1–21. [Google Scholar] [CrossRef]
  11. Daddi-Moussa-Ider, A.; Lisicki, M.; Gekle, S. Slow rotation of a spherical particle inside an elastic tube. Acta Mech. 2018, 229, 149–171. [Google Scholar] [CrossRef] [Green Version]
  12. Romanò, F.; des Boscs, P.-E.; Kuhlmann, H.C. Forces and torques on a sphere moving near a dihedral corner in creeping flow. Eur. J. Mech. B Fluids 2020, 84, 110–121. [Google Scholar] [CrossRef]
  13. Jeffery, G.B. On the steady rotation of a solid of revolution in a viscous fluid. Proc. Lond. Math. Soc. 1915, 14, 327–338. [Google Scholar] [CrossRef]
  14. Keh, H.J.; Chang, J.H. Boundary effects on the creeping-flow and thermophoretic motions of an aerosol particle in a spherical cavity. Chem. Eng. Sci. 1998, 53, 2365–2377. [Google Scholar] [CrossRef]
  15. Lee, T.C.; Keh, H.J. Slow motion of a spherical particle in a spherical cavity with slip surfaces. Int. J. Eng. Sci. 2013, 69, 1–15. [Google Scholar] [CrossRef]
  16. Papavassiliou, D.; Alexander, G.P. Exact solutions for hydrodynamic interactions of two squirming spheres. J. Fluid Mech. 2017, 813, 618–646. [Google Scholar] [CrossRef] [Green Version]
  17. Chou, C.Y.; Keh, H.J. Slow rotation of a spherical particle in an eccentric spherical cavity with slip surfaces. Eur. J. Mech. B Fluids 2021, 86, 150–156. [Google Scholar] [CrossRef]
  18. Brenner, H.; Sonshine, R.M. Slow viscous rotation of a sphere in a circular cylinder. Quart. J. Mech. Appl. Math. 1964, 17, 55–63. [Google Scholar] [CrossRef]
  19. Greenstein, T.; Schiavina, G.L. Torque exerted on a slowly rotating eccentrically positioned sphere within an infinitely long circular cylinder. Int. J. Multiph. Flow 1975, 2, 353–355. [Google Scholar] [CrossRef]
  20. Lee, M.C.; Keh, H.J. Slow axisymmetric rotation of a sphere in a circular tube with slip surfaces. Fluid Dyn. Res. 2021, 53, 065502. [Google Scholar] [CrossRef]
  21. Dean, W.R.; O’Neill, M.E. A slow motion of viscous liquid caused by the rotation of a solid sphere. Mathematika 1963, 10, 13–24. [Google Scholar] [CrossRef]
  22. Chen, P.Y.; Keh, H.J. Slow motion of a slip spherical particle parallel to one or two plane walls. J. Chin. Inst. Chem. Eng. 2003, 34, 123–133. [Google Scholar]
  23. Liao, J.C.; Keh, H.J. Slow rotation of a sphere about its diameter normal to two planes with slip surfaces. Fluid Dyn. Res. 2022, 54, 035502. [Google Scholar] [CrossRef]
  24. Srinivasacharya, D.; Krishna Prasad, M. Steady rotation of a composite sphere in a concentric spherical cavity. Acta Mech. Sin. 2012, 28, 653–658. [Google Scholar] [CrossRef]
  25. Prakash, J.; Raja Sekhar, G.P. Slow motion of a porous spherical particle with a rigid core in a spherical fluid cavity. Meccanica 2017, 52, 91–105. [Google Scholar] [CrossRef]
  26. Sherief, H.H.; Faltas, M.S.; Saad, E.I. Stokes resistance of a porous spherical particle in a spherical cavity. Acta Mech. 2016, 227, 1075–1093. [Google Scholar] [CrossRef]
  27. Chou, C.Y.; Keh, H.J. Low-Reynolds-number rotation of a soft particle inside an eccentric cavity. Eur. J. Mech. B Fluids 2022, 91, 194–201. [Google Scholar] [CrossRef]
  28. Jhuang, L.J.; Keh, H.J. Slow axisymmetric rotation of a soft sphere in a circular cylinder. Eur. J. Mech. B Fluids 2022, 95, 205–211. [Google Scholar] [CrossRef]
  29. Ganatos, P.; Weinbaum, S.; Pfeffer, R. A strong interaction theory for the creeping motion of a sphere between plane parallel boundaries Part 2 Parallel motion. J. Fluid Mech. 1980, 99, 755–783. [Google Scholar] [CrossRef]
  30. Chang, Y.C.; Keh, H.J. Slow motion of a slip spherical particle perpendicular to two plane walls. J. Fluids Struct. 2006, 22, 647–661. [Google Scholar] [CrossRef]
Figure 1. Geometrical sketch of a soft spherical particle rotating about a diameter normal to two planar walls.
Figure 1. Geometrical sketch of a soft spherical particle rotating about a diameter normal to two planar walls.
Colloids 07 00018 g001
Figure 2. Normalized torque T / T 0 for a porous sphere ( a = 0 ) rotating about a diameter perpendicular to two planar walls vs. the shielding parameter λ b : (a) b / c = 9 / 10 ; (b) c / ( c + d ) = 1 / 2 .
Figure 2. Normalized torque T / T 0 for a porous sphere ( a = 0 ) rotating about a diameter perpendicular to two planar walls vs. the shielding parameter λ b : (a) b / c = 9 / 10 ; (b) c / ( c + d ) = 1 / 2 .
Colloids 07 00018 g002
Figure 3. Normalized torque T / T 0 for a porous sphere ( a = 0 ) rotating about a diameter perpendicular to two planar walls vs. the spacing parameter b / c : (a) c / ( c + d ) = 1 / 2 ; (b) λ b = 10 .
Figure 3. Normalized torque T / T 0 for a porous sphere ( a = 0 ) rotating about a diameter perpendicular to two planar walls vs. the spacing parameter b / c : (a) c / ( c + d ) = 1 / 2 ; (b) λ b = 10 .
Colloids 07 00018 g003
Figure 4. Normalized torque T / T 0 for a porous sphere ( a = 0 ) rotating about a diameter perpendicular to two planar walls vs. the relative position parameter c / ( c + d ) : (a) λ b = 10 , where the dashed curves show given values of 2 b / ( c + d ) (the ratio of the particle diameter to the wall-to-wall distance); (b) b / c = 9 / 10 .
Figure 4. Normalized torque T / T 0 for a porous sphere ( a = 0 ) rotating about a diameter perpendicular to two planar walls vs. the relative position parameter c / ( c + d ) : (a) λ b = 10 , where the dashed curves show given values of 2 b / ( c + d ) (the ratio of the particle diameter to the wall-to-wall distance); (b) b / c = 9 / 10 .
Colloids 07 00018 g004
Figure 5. Normalized torque T / T 0 for a soft sphere rotating about a diameter perpendicular to two planar walls vs. the ratio of the core-to-particle radii a / b : (a) b / c = 9 / 10 ; (b) λ b = 1 . The solid and dashed curves denote cases c / ( c + d ) = 1 / 2 and c / ( c + d ) = 0 , respectively.
Figure 5. Normalized torque T / T 0 for a soft sphere rotating about a diameter perpendicular to two planar walls vs. the ratio of the core-to-particle radii a / b : (a) b / c = 9 / 10 ; (b) λ b = 1 . The solid and dashed curves denote cases c / ( c + d ) = 1 / 2 and c / ( c + d ) = 0 , respectively.
Colloids 07 00018 g005
Figure 6. Normalized torque T / T 0 for a soft sphere rotating about a diameter perpendicular to two planar walls vs. the shielding parameter λ b with b / c = 9 / 10 . The solid and dashed curves denote cases of c / ( c + d ) = 1 / 2 and c / ( c + d ) = 0 , respectively.
Figure 6. Normalized torque T / T 0 for a soft sphere rotating about a diameter perpendicular to two planar walls vs. the shielding parameter λ b with b / c = 9 / 10 . The solid and dashed curves denote cases of c / ( c + d ) = 1 / 2 and c / ( c + d ) = 0 , respectively.
Colloids 07 00018 g006
Figure 7. Normalized torque T / T 0 for a soft sphere rotating about a diameter perpendicular to two planar walls vs. the spacing parameter b / c with c / ( c + d ) = 1 / 2 . The solid and dashed curves denote cases of λ b = 1 and λ b = 3 , respectively.
Figure 7. Normalized torque T / T 0 for a soft sphere rotating about a diameter perpendicular to two planar walls vs. the spacing parameter b / c with c / ( c + d ) = 1 / 2 . The solid and dashed curves denote cases of λ b = 1 and λ b = 3 , respectively.
Colloids 07 00018 g007
Figure 8. Normalized torque T / T 0 for a soft sphere rotating about a diameter perpendicular to two planar walls vs. the relative position parameter c / ( c + d ) with λ b = 1 . The solid and dashed curves denote cases of b / c = 7 / 10 and b / c = 9 / 10 , respectively.
Figure 8. Normalized torque T / T 0 for a soft sphere rotating about a diameter perpendicular to two planar walls vs. the relative position parameter c / ( c + d ) with λ b = 1 . The solid and dashed curves denote cases of b / c = 7 / 10 and b / c = 9 / 10 , respectively.
Colloids 07 00018 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chang, C.L.; Keh, H.J. Slow Rotation of a Soft Colloidal Sphere Normal to Two Plane Walls. Colloids Interfaces 2023, 7, 18. https://doi.org/10.3390/colloids7010018

AMA Style

Chang CL, Keh HJ. Slow Rotation of a Soft Colloidal Sphere Normal to Two Plane Walls. Colloids and Interfaces. 2023; 7(1):18. https://doi.org/10.3390/colloids7010018

Chicago/Turabian Style

Chang, Chia L., and Huan J. Keh. 2023. "Slow Rotation of a Soft Colloidal Sphere Normal to Two Plane Walls" Colloids and Interfaces 7, no. 1: 18. https://doi.org/10.3390/colloids7010018

APA Style

Chang, C. L., & Keh, H. J. (2023). Slow Rotation of a Soft Colloidal Sphere Normal to Two Plane Walls. Colloids and Interfaces, 7(1), 18. https://doi.org/10.3390/colloids7010018

Article Metrics

Back to TopTop