# Slow Rotation of Concentric Spheres with Source at Their Centre in a Viscous Fluid

Authors: Deepak Kumar Srivastava
Journal: Journal of Applied Mathematics
 (2009)

---

## Abstract

The problem of concentric pervious spheres carrying a fluid source at their centre and rotating slowly with different uniform angular velocitiesΩ1, Ω2 about a diameter has been studied. The analysis reveals that only azimuthal component of velocity exists, and the couple, rate of dissipated energy is found analytically in the present situation. The expression of couple on inner sphere rotating slowly with uniform angular velocity Ω1, while outer sphere also rotates slowly with uniform angular velocity Ω2, is evaluated. The special cases, like (i) inner sphere is fixed (i.e., Ω1=0), while outer sphere rotates with uniform angular velocity Ω2, (ii) outer sphere is fixed (i.e., Ω2=0), while inner sphere rotates with uniform angular velocity Ω1, and (iii) inner sphere rotates with uniform angular velocity Ω1, while outer sphere rotates at infinity with angular velocity Ω2, have been  deduced.


---

## Body

## 1. Introduction

Stokes flow is becoming increasingly important due to the miniaturization of fluid mechanical parts for example, in micromechanics as well as in nanomechanics. Slow rotation of spheroids (including the disc) in an infinite fluid was first solved by Jeffrey[1] using curvilinear coordinates. His approach was later extended to the spherical lens, torus, and other axisymmetric shapes. Proudman [2] and Stewartson [3] analyzed the dynamical properties of a fluid occupying the space between two concentric rotating spheres when the angular velocities of the spheres are slightly different, in other words, when the motion relative to a reference frame rotating with one of the spheres is due to an imposed azimuthal velocity which is symmetric about the equator. Rubinow and Keller [4] have considered the force on a spinning sphere which is moving through an incompressible viscous fluid by employing the method of matched asymptotic expansions to describe the asymmetric flow. Brenner [5] also obtained some general results for the drag and couple on an obstacle which is moving through the fluid. Childress [6] has investigated the motion of a sphere moving through a rotating fluid and calculated a correction to the drag coefficient. Wakiya[7] numerically evaluated the drag and angular velocity experienced by freely rotating spheres and compared with calculated from corresponding approximate formulae known before. Barrett [8] has tackled the problem of impulsively started sphere rotating with angular velocity Ω about a diameter. He modified the standard time-dependent boundary layer equation to give series solutions satisfying all the boundary conditions and gave solutions that are applicable at small times for non-zero Reynolds numbers. He found that the velocity components decay algebraically rather than exponentially at large distances. Pearson [9] has presented the numerical solution for the time-dependent viscous flow between two concentric rotating spheres. He governed the motion of a pair of coupled nonlinear partial differential equations in three independent variables, with singular end conditions. He also described the computational process for cases in which one (or both) of the spheres is given an impulsive change in angular velocity-starting from a state of either rest or uniform rotation. Majumdar [10] has solved, by using bispherical coordinates, the nonaxisymmetrical Stokes flow of an incompressible homogeneous viscous liquid in space between two eccentric spheres. It was proved that the resultant force acting upon the spheres is at right angles to the axis of rotation and the line of centres. The effect of the stationary sphere on the force and couple exerted by the liquid on the rotating sphere has been discussed, and the results are compared with those of the axisymmetrical case of Jeffrey [1]. Kanwal [11] has considered a disk performing simple harmonic rotary oscillations about its axis of symmetry in a nonconducting viscous fluid which is at rest at infinity. O’Neill and Majumdar [12, 13] have discussed the problem of asymmetrical slow viscous fluid motions caused by the translation or rotation of two spheres. The exact solutions for any values of the ratio of radii and separation parameters are found by them.Ranger [14] tackled the problem of axially symmetric flow past a rotating sphere due to a uniform stream of infinity. He has shown that leading terms for the flow consist of a linear superposition of a primary Stokes flow past a nonrotating sphere together with an anti symmetric secondary flow in the azimuthal plane induced by the spinning sphere. Philander [15] presented a note on the flow properties of a fluid between concentric spheres. This note concerns the flow properties of a spherical shell of fluid when motion is forced across the equator. The fluid under consideration is contained between two concentric spheres which rotate about a diameter with angular velocity Ω. The consequences of the forcing motion across the equator are explored in his work. Cooley [16] has investigated the problem of fluid motion generated by a sphere rotating close to a fixed sphere about a diameter perpendicular to the line of centres in the case when the motion is sufficiently slow to permit the linearization of the Navier-Stokes equations by neglecting the inertia terms. He used a method of matched asymptotic expansions to find asymptotic expressions for the forces and couples acting on the spheres as the minimum clearance between them tends to zero. In his paper, the forces and couples are shown to have the form a0 ln ε+a1+o(ε ln ε), where ε is the ratio of the minimum clearance between the spheres and the radius of the rotating sphere and where a0 and a1 are found explicitly. Munson and Joseph [17, 18] have obtained the high-order analytic perturbation solution for the viscous incompressible flow between concentric rotating spheres. In the second part of their analysis, they have applied the energy theory of hydrodynamic stability to the viscous incompressible flow of a fluid contained between two concentric spheres which rotate about a common axis with prescribed angular velocities. Riley and Mack [19] has discussed the thermal effects on slow viscous flow between rotating concentric spheres. Takagi [20] has considered the flow around a spinning sphere moving in a viscous fluid. He solved the Navier-Stokes equations, using the method of matched asymptotic expansions for small values of the Reynolds number. With the solution, the force and torque on the sphere are computed, and he found that the sphere experiences a force orthogonal to its direction of motion and that the drag is increased in proportion to the square of the spin velocity. Takagi [21] has studied the Stokes flow for the case in which two solid spheres in contact are steadily rotating with different angular velocities about their line of centres. For the case of two equal spheres, he found that, one of which is kept rotating with angular velocity ω while the other is left free, the latter will rotate with angular velocity ω/7. Munson and Menguturk [22] have studied the stability of flow of a viscous incompressible fluid between a stationary outer sphere and rotating inner sphere theoretically and experimentally. Wimmer [23] has provided some experimental results on  incompressible viscous fluid flow in the gap between two concentric rotating spheres. Takagi [24] further studied the problem of steady flow which is induced by the slow rotation of a solid sphere immersed in an infinite incompressible viscous fluid, on the basis of Navier-Stokes equations. He obtained the solution in the form of power series with respect to Reynolds number. Drew [25] has found the force on a small sphere translating relative to a slow viscous flow to order of the 1/2 power of Re for two different fluid flows far from the sphere, namely, pure rotation and pure shear. For pure rotation, the correction of this order to the Stokes drag consists of an increase in the drag. Kim [26] has calculated the torque and frictional force exerted by a viscous fluid on a sphere rotating on the axis of a circular cone of arbitrary vertex angle about an axis perpendicular to the cone axis in the Stokes approximation. Dennis et al. [28] have investigated the problem of viscous incompressible, rotationally symmetric flow due to the rotation of a sphere with a constant angular velocity about a diameter. The solutions of the finite-difference equations are presented for Reynolds number ranging from 1.0 to 5000. Davis and Brenner [29] have used the matched asymptotic expansion methods to solve the problem of steady rotation of a tethered sphere at small, nonzero Reynolds numbers. They obtained first-order Taylor number correction to both the Stokes-law drag and Kirchhoff’s law couple on the sphere for Rossby numbers of order unity. Gagliardi [30] has developed the boundary conditions for the equations of motion for a viscous incompressible fluid in a rotating spherical annulus. The solutions of the stream and circumferential functions were obtained in the form of a series of powers of the Reynolds number. Transient profiles were obtained for the dimensional torque, dimensionless angular velocity of the rotating sphere, and the dimensionless angular momentum of the fluid. Marcus and Tuckerman [31, 32] have computed numerically the steady and translation simulation of flow between concentric rotating spheres. O’Neill and Yano [33] derived the boundary condition at the surfactant and substrate fluids caused by the slow rotation of a solid sphere which is partially submerged in the substrate fluid. Yang et al. [34] have provided the numerical schemes for the problem of the axially symmetric  motion of an incompressible viscous fluid in an annulus between two concentric rotating spheres. Gagliardi et al. [35] reported the study of the steady state and transient motion of a system consisting of an incompressible, Newtonian fluid in an annulus between two concentric, rotating, rigid spheres. They solved the governing equations for the variable coefficients by separation of variables and Laplace Transform methods. They presented the results for the stream function, circumferential function, angular velocity of the spheres, and torque coefficient as a function of time for various values of the dimensionless system parameters. Ranger [36] has found an exact solution of the Navier-Stokes equations for the axi-symmetric motion (with swirl) representing exponentially time-dependent decay of a solid sphere translating and rotating in a viscous fluid relative to a uniform stream whose speed also decays exponentially with time. He also described a similar solution for the two-dimensional analogue where the sphere is replaced by a circular cylinder of infinite length. Tekasakul et al. [37] have studied the problem of the rotatory oscillation of an axi-symmetric body in an axi-symmetric viscous flow at low Reynolds numbers. They evaluated numerically the local stresses and torques on a selection of free, oscillating, axi-symmetric bodies in the continuum regime in an axi-symmetric viscous incompressible flow. Datta and Srivastava [38]  have tackled the problem of slow rotation of a sphere with fluid source at its centre in a viscous fluid. In their investigation, it was found that the effect of fluid source at the centre is to reduce the couple on slowly rotating sphere about its diameter. Kim and Choi [39]  conducted the numerical simulations for laminar flow past a sphere rotating in the streamwise direction, in order to investigate the effect of the rotation on the characteristics of flow over the sphere.Tekasakul and Loyalka [41] have investigated the rotary oscillations of several axi-symmetric bodies in axi-symmetric viscous flows with slip. A numerical method based on Green’s function technique is used, and analytic solutions for local stress and torque on spheres and spheroids as function of the frequency parameter and the slip coefficients are obtained. They have analysed that, in all cases, slip reduces stress and torque, and increasingly so with the increasing frequency parameter. Liu et al. [42] have developed a very efficient numerical method based on the finite difference technique for solving time-dependent nonlinear flow problems. They have applied this method to study the unsteady axisymmetric isotherm flow of an incompressible viscous fluid in a spherical shell with a stationary inner sphere and a rotating outer sphere. Ifidon [43] numerically investigated the problem of determining the induced steady axially symmetric motion of an incompressible viscous fluid confined between two concentric spheres, with the outer sphere rotating with constant angular velocity and the inner sphere fixed for large Reynolds numbers. Davis [44] obtained the expression for force and torque on a rotating sphere close to and within a fluid-filled rotating sphere. Romano [45] has introduced new exact analytic solutions for the rotational motion of a axially symmetric rigid body having two equal principal moments of inertia and subjected to an external torque which is constant in magnitude.In the present paper, the problem of slow rotation of concentric spheres, both assumed to be pervious, with a source at their centre has been tackled. If the strengthQ of the source was of the same order as the angular velocity Ω of  rotating spheres, the inertia terms could still be neglected and the total flow then consists of only the source solution superimposed on the Stokes solution. Therefore, in this case the Stokes drag and couple are not affected by the source. Also, if Q is large enough so that QΩ is not negligible, the inertia terms, being nonlinear, cannot be altogether omitted. The Navier-Stokes  equation can still be linearized by assuming that the velocity perturbation in the source flow on account of the Stokes flow is small so that the terms containing square of angular velocity (i.e., of order Ω2) can be neglected. This assumption is justifiable at least in the vicinity of the spheres where the Stokes approximation is valid. The present problem corresponds to the problem of Stokes flow past a sphere with source at its centre investigated by Datta [46] and slow rotation of sphere with source at its centre in a viscous fluid investigated by Datta and Srivastava [38], the results of which have found application in investigation of the diffusiophoresis target efficiency for an evaporating or condensing drop [47].


## 2. Formulation of the Problem

Let us consider two pervious spheres of radius “a” and “b” (where b>a) with source of strength “Q” at its centre generating radial flow around it in an infinite expense of incompressible fluid of density ρ and kinematic viscosity ν. The spheres are also made to rotate with small steady angular velocities Ω1 and Ω2 so that terms of an o(Ω2) may be neglected but terms of o(QΩ) retained. The motivation of this formulation has been taken from the author’s previous work [38] due to the fact that body geometry has not been changed, although the two concentric spheres are rotating slowly with different angular velocities instead of only one. The governed equations of motion will remain same and provide the new solutions under the defined boundary conditions.The motion is governed by Navier-Stokes equations(2.1)u·gradu=-(1ρ)gradp+ν∇2u
and continuity equation(2.2)divu=0,
together with no-slip boundary condition2.3au=aΩêx×êr,on the inner spherer=a,2.3bu=bΩêx×êr,on the outer spherer=b,
and the condition of vanishing of velocity at far off points(2.4)u=0asr→∞.
In the above equations symbols u,p,ρ,ν stand for velocity, pressure, density, and kinematic viscosity, and unit vectors  êx and êr are along x-axis and radial direction. It will be convenient to work in spherical polar coordinates (r,θ,φ) with x-axis as the polar axis. We nondimensionalize the space variables by a, velocity by a Ω, and pressure by  ρνΩ. Moreover, the symmetry of the problem and the boundary conditions ensure that velocity components are vr=vθ=0 and then we may express the nondimensional velocity vector u as(2.5)u=Qa2r2êr+vφ(r,θ)êφ
and pressure as(2.6)p=ρνΩ[p0(r)+p1(r,θ)].
By introducing the expressions (2.5) and (2.6) in (2.1), the azimuthal component vφ is seen to satisfy the equation(2.7)∇2vφ-vφr2sinθ=sr3∂∂r(rvφ),
where  s=Q/(νa) is the source parameter.The above equation is to be solved under the boundary conditions(2.8)vφ=sinθatr=1(nondimensional equation of spheres),vφ→0asr→∞.


## 3. Solution

We take the trial solution as(3.1)vφ=rω(r)sinθ;
substituting this value of  vφ  in (2.7), we get, after some calculation  and adjustment,(3.2)ddr[r4dωdr-sr2ω]=0,
and the boundary conditions (2.8) in nondimensional form become(3.3)ω=1atr=1(i.e., on the surface),ω→0asr→∞.
The above boundary conditions may also be express in dimensional form as(3.4)ω=Ω1atr=a(i.e., on the inner sphere),ω→0asr→∞.
On integration of (3.2), we get the solution in nondimensional form as(3.5)ω(r)=-As3[s3r2-2sr+2]+Be-s/r
and in dimensional form as(3.6)ω(r)=-As3[s2a2r2-2sar+2]+Be-sa/r,
where A and B are constants of integration which can be obtained by applying boundary conditions (3.4) as(3.7)-As3=Ω1[s2-2s+2-2e-s]-1,B=-2Ω1[s2-2s+2-2e-s]-1.
Substituting the values of constants A and B in (3.6), we get the expression of angular velocity ω(r) in dimensional form as(3.8)ω(r)=-As3[s2a2r2-2sar+2]+Be-sa/r
or(3.9)ω(r)=Ω1[s2a2r2-2sar+2-e-sa/r][s2-2s+2-2e-s]-1
and consequently, with the help of (3.1), the expression for azimuthal component of velocity vφ comes out to be in dimensional form as(3.10)vφ=rω(r)sinθ=Ωrsinθ[s2a2r2-2sar+2-e-(sa/r)][s2-2s+2-2e-s].


### 3.1. Couple on Inner Sphere Rotating with Uniform Angular Velocity Ω1 (When Outer Sphere Is Also Rotating with Different Uniform Angular Velocity Ω2)

If  there exists an external concentric pervious sphere of radiusb(b>a), rotating with small angular velocity Ω2, that is, the boundary conditions for this situation will be(3.11)ω=Ω2atr=b(at outer surface),ω=Ω1atr=a(at inner surface),
by using the above boundary conditions, in (3.6), the constant of integration A and B comes out to be(3.12)As3=Ω1-Ω2e-s(1-a/b)[(s2a2/b2-2(sa/b)+2)e-s(1-a/b)+(-s2+2s-2)],B=esa/b[Ω1{s2a2/b2-2(sa/b)+2}+Ω2{-s2+2s-2}e-s(1-a/b){s2a2/b2-2(sa/b)+2}+{-s2+2s-2}].
The expression of angular velocity ω(r) can be written with the help of (3.6)(3.13)ω(r)=-As3[s2a2r2-2sar+2]+Be-sa/r,
where A and B are given in (3.12). On differentiating the function ω(r), we have(3.14)dωdr=-As3[-2s2a2r3+2sar2]+Be-sa/r(sar2),
and the value of  dω/dr  at  r=a  can be written as(3.15)(dωdr)r=a=1a[2As3(s2-s)+Be-s].
The moment of force pφ is pφ·rsinθ, where  pφ=μ·rsinθ·(dω/dr) is the only nonvanishing component of force p. If N is the couple on the sphere of radius a, then by using (3.15), we have(3.16)N=∫0π(pφ·rsinθ)r=adS=∫0π(μrsinθdωdrrsinθ)r=a·(2πasinθ·adθ)=83πa3μ[2As3(s2-s)+Be-s]=83πa3μ[s2a2b22(s2-s){Ω1-Ω2e-s(1-a/b)}+e-s(1-a/b)×{Ω1(s2a2b2-2sab+2)+Ω2(-s2+2s-2)}]×[{e-s(1-a/b)(s2a2b2-2sab+2)+(-s2+2s-2)}]-1.
The rate of dissipated energy is given by NΩ1, where the value of N is given in (3.16).


### 3.2. Couple on Outer Sphere Rotating with Uniform Angular Velocity Ω2 (When Inner Sphere is Fixed, that is, 1=0Ω)

The expression for angular velocityω(r) is given by (3.6):(3.17)ω(r)=-As3[s2a2r2-2sar+2]+Be-sa/r.
Now we use the following boundary conditions:(3.18)ω(r)=Ω2on surfacer=b,ω(r)→0asr→∞.
Under these boundary conditions, the values of constant A and B can be obtained as follows:(3.19)As3=Ω2[s2a2b2+2sab+2e-sa/b-2]-1,B=2Ω2[s2a2b2+2sab+2e-sa/b-2]-1.
Now, the expression for derivative of angular velocity at r=b comes out to be(3.20)[ddrω(r)]r=b=[-As3(-2s2a2r3+2sar2)+B(sar2)e-sa/r]r=b=1b[2As3(s2a2b2-sab)+Bsabe-sa/b],
which reduces in final form by (3.19):(3.21)=2Ω2b[s2a2b2-sab+asbe-sa/b][s2a2b2+2sab+2e-sa/b-2]-1.
Hence, the couple N on the outer sphere in the presence of inner sphere is(3.22)N=∫0π(μ·rsinθ·dωdr·rsinθ)r=b(2πbsinθ·bdθ)=83πb3μ(ddrω(r))r=b,
and by using (3.21), it reduces to(3.23)=163πb3μΩ2[s2a2b2-sab+sabe-sa/b][s2a2b2+2sab+2e-sa/b-2]-1.
The expression for rate of dissipated energy will be NΩ2, where N is given in (3.23).


## 3.1. Couple on Inner Sphere Rotating with Uniform Angular Velocity Ω1 (When Outer Sphere Is Also Rotating with Different Uniform Angular Velocity Ω2)

If  there exists an external concentric pervious sphere of radiusb(b>a), rotating with small angular velocity Ω2, that is, the boundary conditions for this situation will be(3.11)ω=Ω2atr=b(at outer surface),ω=Ω1atr=a(at inner surface),
by using the above boundary conditions, in (3.6), the constant of integration A and B comes out to be(3.12)As3=Ω1-Ω2e-s(1-a/b)[(s2a2/b2-2(sa/b)+2)e-s(1-a/b)+(-s2+2s-2)],B=esa/b[Ω1{s2a2/b2-2(sa/b)+2}+Ω2{-s2+2s-2}e-s(1-a/b){s2a2/b2-2(sa/b)+2}+{-s2+2s-2}].
The expression of angular velocity ω(r) can be written with the help of (3.6)(3.13)ω(r)=-As3[s2a2r2-2sar+2]+Be-sa/r,
where A and B are given in (3.12). On differentiating the function ω(r), we have(3.14)dωdr=-As3[-2s2a2r3+2sar2]+Be-sa/r(sar2),
and the value of  dω/dr  at  r=a  can be written as(3.15)(dωdr)r=a=1a[2As3(s2-s)+Be-s].
The moment of force pφ is pφ·rsinθ, where  pφ=μ·rsinθ·(dω/dr) is the only nonvanishing component of force p. If N is the couple on the sphere of radius a, then by using (3.15), we have(3.16)N=∫0π(pφ·rsinθ)r=adS=∫0π(μrsinθdωdrrsinθ)r=a·(2πasinθ·adθ)=83πa3μ[2As3(s2-s)+Be-s]=83πa3μ[s2a2b22(s2-s){Ω1-Ω2e-s(1-a/b)}+e-s(1-a/b)×{Ω1(s2a2b2-2sab+2)+Ω2(-s2+2s-2)}]×[{e-s(1-a/b)(s2a2b2-2sab+2)+(-s2+2s-2)}]-1.
The rate of dissipated energy is given by NΩ1, where the value of N is given in (3.16).


## 3.2. Couple on Outer Sphere Rotating with Uniform Angular Velocity Ω2 (When Inner Sphere is Fixed, that is, 1=0Ω)

The expression for angular velocityω(r) is given by (3.6):(3.17)ω(r)=-As3[s2a2r2-2sar+2]+Be-sa/r.
Now we use the following boundary conditions:(3.18)ω(r)=Ω2on surfacer=b,ω(r)→0asr→∞.
Under these boundary conditions, the values of constant A and B can be obtained as follows:(3.19)As3=Ω2[s2a2b2+2sab+2e-sa/b-2]-1,B=2Ω2[s2a2b2+2sab+2e-sa/b-2]-1.
Now, the expression for derivative of angular velocity at r=b comes out to be(3.20)[ddrω(r)]r=b=[-As3(-2s2a2r3+2sar2)+B(sar2)e-sa/r]r=b=1b[2As3(s2a2b2-sab)+Bsabe-sa/b],
which reduces in final form by (3.19):(3.21)=2Ω2b[s2a2b2-sab+asbe-sa/b][s2a2b2+2sab+2e-sa/b-2]-1.
Hence, the couple N on the outer sphere in the presence of inner sphere is(3.22)N=∫0π(μ·rsinθ·dωdr·rsinθ)r=b(2πbsinθ·bdθ)=83πb3μ(ddrω(r))r=b,
and by using (3.21), it reduces to(3.23)=163πb3μΩ2[s2a2b2-sab+sabe-sa/b][s2a2b2+2sab+2e-sa/b-2]-1.
The expression for rate of dissipated energy will be NΩ2, where N is given in (3.23).


## 4. Particular Cases

Case 1.
We consider the outer spherical surface to be fixed, that is,Ω2=0; then in this case, by (3.16), we have couple on the inner sphere rotating with angular velocity Ω1 as
(4.1)N=83πa3μΩ1[2(s2-s)+(s2a2b2-2sab+2)e-s(1-a/b)]×[(-s2+2s-2)+(s2a2b2-2sab+2)e-s(1-a/b)]-1.
Now, on shifting the solid outer spherical body having radius b(b>a) to infinity that is, b→∞, then e-s(1-a/b)→e-s, and by (4.1), we can have the expression for couple on slowly rotating sphere of radius “a” alone and given by (4.1) as
(4.2)N=163πμa3Ω1[s2-s+se-s][-s2+2s-2+2e-s],
which matches with the expression of couple obtained by Datta and Srivastava [38] for slowly rotating pervious sphere of radius “a” rotating with slow uniform angular velocity Ω1 and further reduces to classical one  8πμa3Ω1 for  s=0 ( i.e., in the absence of source at the center).Case 2.
We consider the inner spherical surface to be fixed, that is,Ω1=0; then in this case, by (3.16), we have the couple on the outer sphere rotating with angular velocity Ω2 as
(4.3)N=83πμa3Ω2[-3s2+4s-2][(s2a2/b2-2(sa/b)+2)+(-s2+2s-2)es(1-a/b)].Case 3.
If the inner sphere rotates with uniform angular velocityΩ1, while outer rotates with uniform angular velocity Ω2 at infinity, that is, b→∞, then by expression (3.16), we have the couple on inner sphere as
(4.4)N=83πμa3[2Ω1(s2-s+e-s)+Ω2e-s(-3s2+4s-2)][-s2+2s-2+2e-s].
For large value of source parameter “s”, (4.4) reduces in to the form
(4.5)N≈163πμa3Ω1[1+(1s-6s3)],
which gives2/3M0, where M0=8πμa3Ω1, in the limit as  s→∝, again in good agreement with the result of couple on rotating sphere with uniform angular velocity for large source parameter obtained by Datta and Srivastava [38].Case 4.
If we consider the limiting situation asb→a  and  Ω2→Ω1, then we have the expression for couple on slowly rotating sphere having radius  “a”  by (3.23)
(4.6)N=163πa3μΩ1[s2-s+se-s][-s2+2s+2e-s-2]-1,
which matches with the result existed in the paper of Datta and Srivastava [38] and further reduces to the classical one 8πμa3Ω1  for  s=0 ( i.e., in the absence of source at the center).
The expressions for couple in general case ((3.16), (3.23)) and in cases ((4.1) to (4.6)) are expected to be new and never seen in literature. It was concluded there that the effect of fluid source at the center of sphere is to reduce the couple.



---

*Source: 740172-2009-11-16.xml*