International Scholarly Research Notices

International Scholarly Research Notices / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 513717 |

E. A. Ashmawy, "Unsteady Rotational Motion of a Slip Spherical Particle in a Viscous Fluid", International Scholarly Research Notices, vol. 2012, Article ID 513717, 8 pages, 2012.

Unsteady Rotational Motion of a Slip Spherical Particle in a Viscous Fluid

Academic Editor: A. Stefanov
Received21 Sep 2011
Accepted13 Oct 2011
Published24 Jan 2012


The unsteady rotational motion of a slip spherical particle with a nonuniform angular velocity in an incompressible viscous fluid flow is discussed. The technique of Laplace transform is used. The slip boundary condition is applied at the surface of the sphere. A general formula for the resultant torque acting on the surface of the sphere is deduced. Special fluid flows are considered and their results are represented graphically.

1. Introduction

The unsteady incompressible viscous fluid flows have been discussed extensively in the literature both analytically and numerically and have been reported in numerous studies. Basset [1] gave the force exerted on a sphere in an arbitrary axisymmetric time-dependent motion by integrating the result of Stokes [2]. Mazur and Bedeaux used the Faxen theorem to obtain the force on a moving sphere when moving with time-dependent velocity in [3].

Sano [4] utilized an asymptotic expansion to solve the problem of the low-Reynolds number unsteady fluid flow past a sphere when a constant rectilinear velocity is suddenly imparted to the sphere. Maxey and Riley [5] discussed the forces acting on a small rigid sphere in a non-uniform viscous fluid flow. In [6], Lawrence and Weinbaum used the linearized Navier-Stokes equations to obtain the force acting on an arbitrary axisymmetric body in oscillatory motion. Lawrence and Weinbaum [7] presented a more general analysis of the unsteady Stokes equations for the axisymmetric flow past a spheroidal body. Mei et al. examined in [8] the dependence of the unsteady drag on the frequency of the fluctuations at various Reynolds numbers. Lovalenti and Brady calculated the hydrodynamic force acting on a rigid spherical particle translating with arbitrary time-dependent motion in a time-dependent flowing fluid for small but finite values of the Reynolds number based on the particle’s slip velocity relative to the uniform flow in [9]. Feng and Joseph used in [10] a numerical simulation to examine the validity of the quasisteady method for solving unsteady viscous fluid flows. Ranger [11] studied time-dependent decay of the motion of a sphere translating and rotating in a viscous liquid. In [12], Coimbra and Rangel used a fractional-differential operator to obtain the general solution of the particle momentum equation for unsteady Stokes flows analytically. Asmolov [13] discussed the flow induced by a sphere undergoing unsteady motion in a Newtonian fluid at small Reynolds number. Venkatalaxmi et al. discussed a simple method of solution for the problem of an arbitrary unsteady Stokes flow in the presence of a sphere with no-slip condition in [14]. Avudainayagam and Ramakrisnna [15] used Saffman’s technique of integral identities to obtain the total couple acting on a sphere rotating in a viscous fluid with classical no-slip boundary condition.

The slip condition introduced by Navier [16] assumes that the surface friction force determined by the tangential velocity of fluid relative to the boundary at a point on its surface is proportional to the tangential stress acting at that point. The constant of proportionality between the fluid and the boundary is called the coefficient of sliding friction. It depends only on the nature of the fluid and the surface of the boundary [17]. While slip conditions are of actual interest for gases, recently it has been found that these conditions are also of interest for liquids as well [18–24]. Albano et al. [25] have obtained an expression to the drag force exerted on a sphere moving with a time-dependent velocity through an incompressible fluid in an unsteady inhomogeneous flow for the case of arbitrary slip on the surface of the sphere. Datta and Deo [26] investigated the forces experienced by randomly and homogeneously distributed parallel circular cylinders or spheres in uniform viscous flow with slip boundary condition under Stokes approximation using particle-in-cell model technique. Senchenko and Keh [27] have obtained analytical approximations for the resistance relations for a rigid, slightly deformed slip sphere in an unbounded Stokesian steady flow. The problem of slow steady rotation of two axisymmetric particles (separated by a certain distance) about their common axis of symmetry in an infinite viscous fluid with slip boundary conditions at their surfaces has been studied numerically by Tekasakul et al. [28]. Keh and Lee [29] discussed the quasisteady creeping flow caused by a spherical fluid or solid particle with a slip surface translating in a viscous fluid within a spherical cavity along the line connecting their centers in the limit of small Reynolds number.

In the present work, we discuss the unsteady slow motion of a rigid sphere with interfacial slip condition in a viscous fluid flow. A simple formula for the torque acting on the surface of the sphere is deduced. The case of steady flow and the case of classical no-slip condition are obtained as special cases of this work.

2. Formulation and Solution of the Problem

Assume that a rigid sphere of radius β€œπ‘Žβ€ rotates suddenly with a nonconstant angular velocity β€œΞ©(𝑑)” about an instantaneous axis of revolution in a viscous fluid which is otherwise at rest. For the creeping motion of a viscous fluid, the equation of motion takes the following form:πœŒπœ•πͺπœ•π‘‘=βˆ’βˆ‡π‘+πœ‡βˆ‡2πͺ,(2.1) where 𝜌 is the density of the fluid, πͺ represents the velocity vector of the fluid, 𝑝 is the pressure, and πœ‡ denotes the dynamical viscosity coefficient. Moreover, βˆ‡ is the usual Hamilton β€œnabla” operator.

The stress tensor can also be represented as follows:𝑑𝑖𝑗=βˆ’π‘π›Ώπ‘–π‘—+2πœ‡π‘’π‘–π‘—,with𝑒𝑖𝑗=12ξ€·π‘žπ‘–,𝑗+π‘žπ‘—,𝑖.(2.2) Working with the spherical polar coordinates (π‘Ÿ,πœƒ,πœ™),with origin at the centre of the sphere, the only non vanishing velocity component is the azimuthal πœ™-component denoted byπ‘žπœ™. Introducing the Laplace transform (denoted by an over bar) defined by ξ€œπ‘“(π‘Ÿ,πœƒ,𝑠)=∞0𝑓(π‘Ÿ,πœƒ,𝑑)π‘’βˆ’π‘ π‘‘π‘‘π‘‘,(2.3) the equation of motion then reduces to the following form:𝐸2βˆ’β„“2ξ€Έξ€½π‘Ÿsinπœƒπ‘žπœ™ξ€Ύ(π‘Ÿ,πœƒ,𝑠)=0,(2.4) where β„“2=𝑠/𝜈 and 𝜈represents the kinematical viscosity coefficient. Moreover the operator β€œπΈβ€ is defined by𝐸2=πœ•2πœ•π‘Ÿ2βˆ’cotπœƒπ‘Ÿ2πœ•+1πœ•πœƒπ‘Ÿ2πœ•2πœ•πœƒ2.(2.5)

The initial condition is taking the following form:π‘žπœ™ξ€·π‘Ÿ,πœƒ,0+ξ€Έ=ξ‚»Ξ©(0)π‘Žsinπœƒ,π‘Ÿ=π‘Ž,0,π‘Ÿ>π‘Ž.(2.6) Moreover, the following slip boundary condition is applied on the surface of the rigid sphere:π›½ξ€½π‘žπœ™ξ€Ύ(π‘Ž,πœƒ,𝑑)βˆ’Ξ©(𝑑)π‘Žsinπœƒ=π‘‘π‘Ÿπœ™||π‘Ÿ=π‘Ž,𝑑>0,(2.7) where 𝛽 is the slip parameter. It is a measure of the degree of tangential slip existing between the fluid and the boundary. This coefficient depends only on the nature of the fluid and solid surface. The classical case of no-slip can be recovered as a special case when π›½β†’βˆž.

Applying the boundary conditions (2.6)-(2.7), in the Laplace transform domain, to the solution of the differential equation (2.4), we obtain the velocity π‘žπœ™ in the following form:π‘žπœ™=Ω(𝑠)π›Όβ„“π‘Ž3β„“2π‘Ž2+ξƒͺξ‚€1(β„“π‘Ž+1)(𝛼+3)π‘Ÿ+1β„“π‘Ÿ2ξ‚π‘’βˆ’β„“(π‘Ÿβˆ’π‘Ž)sinπœƒ,(2.8) where𝛼=π‘Žπ›½πœ‡.(2.9) This later dimensionless parameter 𝛼 is known as Navier number. Moreover, the ratio of πœ‡/𝛽 is called the slip length.

Also, the shear stress has the following form:π‘‘π‘Ÿπœ™ξƒ©=βˆ’Ξ©(𝑠)π›Όπ‘Ž3πœ‡β„“2π‘Ž2+ξƒͺξ‚΅β„“(β„“π‘Ž+1)(𝛼+3)2π‘Ÿ+3β„“π‘Ÿ2+3π‘Ÿ3ξ‚Άπ‘’βˆ’β„“(π‘Ÿβˆ’π‘Ž)sinπœƒ.(2.10) The torque acting on the surface of the rigid sphere can be obtained, in the Laplace transform domain, by 𝑇𝑧=2πœ‹π‘Ž3ξ€œπœ‹0π‘‘π‘Ÿπœ™||π‘Ÿ=π‘Žsin2πœƒπ‘‘πœƒ.(2.11) This can be simplified to the following form:𝑇𝑧=βˆ’8πœ‹π‘Ž3πœ‡π›ΌβŽ›βŽœβŽœβŽ(𝛼+3)Ξ©(𝑠)+𝛼𝑠Ω(𝑠)3ξ‚€βˆšπ‘ +((𝛼+3)/π‘Ž)π‘ πœˆ+(𝛼+3)𝜈/π‘Ž2ξ‚βŽžβŽŸβŽŸβŽ .(2.12) Whenπ›Όβ†’βˆž, we recover the results obtained for the case of no-slip given by [15]𝑇𝑧=βˆ’8πœ‹π‘Ž3πœ‡βŽ›βŽœβŽœβŽπ‘ŽΞ©(𝑠)+2𝑠Ω(𝑠)3ξ‚ƒπ‘Žβˆšξ‚„βŽžβŽŸβŽŸβŽ π‘ πœˆ+𝜈.(2.13) Taking the inverse Laplace transform of (2.12), we arrive at𝑇𝑧(𝑑)=βˆ’8πœ‹π‘Ž3πœ‡π›Όξƒ©π›Ό(𝛼+3)Ξ©(𝑑)+3πΏβˆ’1𝑠Ω(𝑠)𝑠+(𝛼+3)π‘Žβˆšπ‘ πœˆ+(𝛼+3)πœˆπ‘Ž2ξ‚Άβˆ’1ξƒ°ξƒͺ,(2.14) where πΏβˆ’1represents the inverse Laplace operator.

By using the convolution theorem together with the complex Laplace inversion formula with the aid of contour integration, we can obtain the inverse Laplace transform of the formula (2.14) in the following form:𝑇𝑧(𝑑)=βˆ’8πœ‹π‘Ž3πœ‡π›Όξ‚΅ξ€œ(𝛼+3)Ξ©(𝑑)+𝑑0ξ‚€π‘‘Ξ©ξ‚ξ‚Άπ‘‘πœ+Ξ©(0)𝐹(π‘‘βˆ’πœ)π‘‘πœ,(2.15) whereπ‘ŽπΉ(𝑑)=3√3πœ‹πœˆπœˆξ€œβˆž0√π‘₯π‘’βˆ’π‘₯𝑑1βˆ’π‘₯π‘Ž2ξ€Έ/𝜈(𝛼+3)2+π‘₯π‘Ž2/πœˆπ‘‘π‘₯.(2.16) The general formula (2.15) can be employed to obtain the total torque acting on a sphere rotating in a Stokes viscous fluid flow with a nonuniform angular velocity.

The case of no-slip can be recovered in the limiting case of π›Όβ†’βˆž to take the following form:𝑇𝑧(𝑑)=βˆ’8πœ‹π‘Ž3πœ‡ξƒ©π‘ŽΞ©(𝑑)+3√3πœ‹πœˆπœˆξ€œπ‘‘0𝑑Ω(𝜏)ξ‚Άξ€œπ‘‘πœ+Ξ©(0)∞0√π‘₯π‘’βˆ’π‘₯(π‘‘βˆ’πœ)1+π‘₯π‘Ž2ξƒͺ/πœˆπ‘‘π‘₯π‘‘πœ.(2.17) To recover the result of Basset [1] for the torque acting on a rigid sphere, of radius π‘Ž, rotating steadily in a viscous fluid with a slip-flow boundary condition, we put Ξ©(𝑑)=Ξ©π‘œ(constant)andπ‘‘β†’βˆž into (2.15) to get𝑇𝑧(𝑑)=βˆ’8πœ‹π‘Ž3πœ‡π›ΌΞ©π‘œ(𝛼+3).(2.18)

2.1. Some Special Flows

In this section we will employ the obtained general formula (2.15) to some non-uniform viscous fluid flows.

Case 1 (Damping oscillation). Here we assume that the rigid sphere starts to move with the angular velocity: Ξ©(𝑑)=Ξ©π‘œπ‘’βˆ’πœ”π‘‘sin(πœ”π‘‘).(2.19) In this case, the torque acting on the sphere is given by 𝑇𝑧(𝑑)=βˆ’8πœ‹π‘Ž3πœ‡π›ΌΞ©π‘œξ‚»π‘’(𝛼+3)βˆ’π‘€π‘‘ξ€œsin(πœ”π‘‘)+πœ”π‘‘0π‘’βˆ’πœ”πœξ‚Ό(cos(πœ”πœ)βˆ’sin(πœ”πœ))𝐹(π‘‘βˆ’πœ)π‘‘πœ,(2.20) where Ξ©π‘œ is the characteristic angular velocity and πœ” is the oscillation angular velocity.

Case 2 (Accelerating velocity). In the second special case we assume that the sphere starts to rotate with the angular velocity of the following form: Ξ©(𝑑)=Ξ©π‘œ{πœ”π‘‘+𝑐},(2.21) where 𝑐 is an arbitrary constant.
The torque in this case takes the following form: 𝑇𝑧(𝑑)=βˆ’8πœ‹π‘Ž3πœ‡π›ΌΞ©π‘œξƒ―(𝛼+3)(πœ”π‘‘+𝑐)+πœ”π‘Ž3√3πœ‹πœˆπœˆξ€œπ‘‘0ξ€œβˆž0√π‘₯π‘’βˆ’π‘₯(π‘‘βˆ’πœ)ξ€·1βˆ’π‘₯π‘Ž2ξ€Έ/𝜈(𝛼+3)2+π‘₯π‘Ž2ξƒ°/πœˆπ‘‘π‘₯π‘‘πœ.(2.22)

Case 3 (Impulsive rotation). Assume that the sphere rotates impulsively by the angular velocity: Ξ©(𝑑)=Ξ©o𝐻(𝑑),(2.23) where 𝐻(𝑑) is the Heaviside step function.
The total Torque is found to be𝑇𝑧(𝑑)=βˆ’8πœ‹π‘Ž3πœ‡π›ΌΞ©π‘œξƒ―π‘Ž(𝛼+3)𝐻(𝑑)+3√3πœ‹πœˆπœˆξ€œβˆž0√π‘₯π‘’βˆ’π‘₯𝑑1βˆ’π‘₯π‘Ž2ξ€Έ/𝜈(𝛼+3)2+π‘₯π‘Ž2ξƒ°/πœˆπ‘‘π‘₯.(2.24)

3. Conclusion

In this paper, we have applied the slip condition and obtained a general formula that can be used to evaluate the torque on a sphere rotating with a non-uniform angular velocity in a viscous fluid flow. The classical no-slip case can be recovered as a special case of our work in the limiting case ofπ›Όβ†’βˆž. Moreover, the couple acting on a sphere moving steadily in a viscous fluid flow can be obtained from Case 3 of this paper whenπ‘‘β†’βˆž. From Figure 1, we conclude that the torque in the case of damping oscillation goes to zero after short time. Also, the torque in the case of accelerating velocity represented in Figure 2 is increasing linearly with the time. Figure 3 shows that after very short time the torque tends to the steady state as expected. From Figures 1, 2, and 3 we note that the torque is identically zero in the case of perfect slip, that is, when 𝛼=0. Also, we observe that the values of the total couple increase with the increase of the slip parameter.


  1. A. B. Basset, A Treatise on Hydrodynamics, Dover Publications, New York, NY, USA, 1961.
  2. G. G. Stokes, β€œOn the theories of internal friction of fluids in motion,” Transactions of the Cambridge Philosophical Society, vol. 8, pp. 287–305, 1845. View at: Google Scholar
  3. P. Mazur and D. Bedeaux, β€œA generalization of Faxén's theorem to nonsteady motion of a sphere through an incompressible fluid in arbitrary flow,” Physica, vol. 76, no. 2, pp. 235–246, 1974. View at: Google Scholar
  4. T. Sano, β€œUnsteady flow past a sphere at low Reynolds number,” Journal of Fluid Mechanics, vol. 112, pp. 433–441, 1981. View at: Google Scholar
  5. M. R. Maxey and J. J. Riley, β€œEquation of motion for a small rigid sphere in a nonuniform flow,” Physics of Fluids, vol. 26, no. 4, pp. 883–889, 1983. View at: Google Scholar
  6. C. J. Lawrence and S. Weinbaum, β€œThe force on an axisymmetric body in linearized time dependent motion,” Journal of Fluid Mechanics, vol. 171, pp. 209–218, 1986. View at: Google Scholar
  7. C. J. Lawrence and S. Weinbaum, β€œThe unsteady force on a body at low Reynolds number; the axisymmetric motion of a spheroid,” Journal of Fluid Mechanics, vol. 189, pp. 463–489, 1988. View at: Google Scholar
  8. R. Mei, C. J. Lawrence, and R. J. Adrian, β€œUnsteady drag on a sphere at finite Reynolds number with small fluctuations in the free-stream velocity,” Journal of Fluid Mechanics, vol. 233, pp. 613–631, 1991. View at: Google Scholar
  9. P. M. Lovalenti and J. F. Brady, β€œThe hydrodynamic force on a rigid particle undergoing arbitrary time- dependent motion at small Reynolds number,” Journal of Fluid Mechanics, vol. 256, pp. 561–605, 1993. View at: Google Scholar
  10. J. Feng and D. D. Joseph, β€œThe unsteady motion of solid bodies in creeping flows,” Journal of Fluid Mechanics, vol. 303, pp. 83–102, 1995. View at: Google Scholar
  11. K. B. Ranger, β€œTime-dependent decay of the motion of a sphere translating and rotating in a viscous liquid,” Quarterly Journal of Mechanics and Applied Mathematics, vol. 49, no. 4, pp. 621–633, 1996. View at: Google Scholar
  12. C. F. M. Coimbra and R. H. Rangel, β€œGeneral solution of the particle momentum equation in unsteady Stokes flows,” Journal of Fluid Mechanics, vol. 370, pp. 53–72, 1998. View at: Google Scholar
  13. E. S. Asmolov, β€œFlow past a sphere undergoing unsteady rectilinear motion and unsteady drag at small Reynolds number,” Journal of Fluid Mechanics, vol. 446, pp. 95–119, 2001. View at: Google Scholar
  14. A. Venkatalaxmi, B. S. Padmavathi, and T. Amaranath, β€œA general solution of unsteady Stokes equations,” Fluid Dynamics Research, vol. 35, no. 3, pp. 229–236, 2004. View at: Publisher Site | Google Scholar
  15. A. Avudainayagam and J. Ramakrisnna, β€œSlow unsteady motion of a sphere in a viscous incompressible fluid,” Journal of Mathematical and Physical Sciences, vol. 16, pp. 359–365, 1982. View at: Google Scholar
  16. C. L . M. H. Navier, β€œMemoirs de l'Academie,” Royale des Sciences de l'Institut de France, vol. 1, pp. 414–416, 1823. View at: Google Scholar
  17. J. Happel and H. Bernner, Low Reynolds Number Hydrodynamics, Noordhoff International Publishing, The Netherlands, 1973.
  18. J. L. Barrat and L. Bocquet, β€œLarge slip effect at a nonwetting fluid-solid interface,” Physical Review Letters, vol. 82, no. 23, pp. 4671–4674, 1999. View at: Google Scholar
  19. R. Pit, H. Hervet, and L. Léger, β€œDirect experimental evidence of slip in hexadecane: solid interfaces,” Physical Review Letters, vol. 85, no. 5, pp. 980–983, 2000. View at: Publisher Site | Google Scholar
  20. D. C. Tretheway, L. Zhu, L. Petzold, and C. D. Meinhart, β€œExamination of the slip boundary condition by μ-PIV and Lattice Boltzmann simulations,” in Proceedings of the ASME International Mechanical Engineering Congress and Exposition, pp. 551–556, New Orleans, La, USA, November 2002. View at: Google Scholar
  21. M. E. O'Neill, K. B. Ranger, and H. Brenne, β€œSlip at the surface of a translating–rotating sphere bisected by a free surface bounding a semi-infinite viscous fluid, removal of the contact-line singularity,” Physics of Fluids, vol. 29, pp. 913–924, 1986. View at: Google Scholar
  22. E. E. Michaelides and Z. G. Feng, β€œThe equation of motion of a small viscous sphere in an unsteady flow with interface slip,” International Journal of Multiphase Flow, vol. 21, no. 2, pp. 315–321, 1995. View at: Google Scholar
  23. B. S. Padmavathi, T. Amaranath, and S. D. Nigam, β€œStokes flow past a sphere with mixed slip-stick boundary conditions,” Fluid Dynamics Research, vol. 11, no. 5, pp. 229–234, 1993. View at: Google Scholar
  24. D. Palaniappan and P. Daripa, β€œExterior Stokes flows with stick-slip boundary conditions,” Zeitschrift für Angewandte Mathematik und Physik, vol. 53, no. 2, pp. 281–307, 2002. View at: Publisher Site | Google Scholar | MathSciNet
  25. A. M. Albano, D. Bedeaux, and P. Mazur, β€œOn the motion of a sphere with arbitrary slip in a viscous incompressible fluid,” Physica A, vol. 80, no. 1, pp. 89–97, 1975. View at: Google Scholar
  26. S. Datta and S. Deo, β€œStokes flow with slip and Kuwabara boundary conditions,” Proceedings of the Indian Academy of Sciences, vol. 112, no. 3, pp. 463–475, 2002. View at: Google Scholar
  27. S. Senchenko and H. J. Keh, β€œSlipping Stokes flow around a slightly deformed sphere,” Physics of Fluids, vol. 18, no. 8, Article ID 088104, pp. 1–4, 2006. View at: Publisher Site | Google Scholar
  28. P. Tekasakul, R. V. Tompson, and S. K. Loyalka, β€œA numerical study of two coaxial axi-symmetric particles undergoing steady equal rotation in the slip regime,” Zeitschrift für Angewandte Mathematik und Physik, vol. 50, no. 3, pp. 387–403, 1999. View at: Google Scholar
  29. H. J. Keh and T. C. Lee, β€œAxisymmetric creeping motion of a slip spherical particle in a nonconcentric spherical cavity,” Theoretical and Computational Fluid Dynamics, vol. 24, no. 5, pp. 497–510, 2010. View at: Publisher Site | Google Scholar

Copyright © 2012 E. A. Ashmawy. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles