# Analytical Solutions for the Flow of a Fractional Second Grade Fluid due to a Rotational Constantly Accelerating Shear

Authors: M. Kamran; M. Imran; M. Athar
Journal: ISRN Mathematical Physics
 (2012)

---

## Abstract

Exact analytic solutions are obtained for the flow of a generalized second grade fluid in an annular region between two infinite coaxial cylinders. The fractional calculus approach in the governing equations of a second grade fluid is used. The exact analytic solutions are constructed by means of Laplace and finite Hankel transforms. The motion is produced by the inner cylinder which is rotating about its axis due to a constantly accelerating shear. The solutions that have been obtained satisfy both the governing equations and all imposed initial and boundary conditions. Moreover, they can be easily specialized to give similar solutions for second grade and Newtonian fluids. Finally, the influence of the pertinent parameters on the fluid motion, as well as a comparison between the three models, is underlined by graphical illustrations.


---

## Body

## 1. Introduction

The study of the non-Newtonian fluids has recently achieved much importance because of well-established applications in a number of processes that occur in industry. Such applications include the extrusion of polymer fluids, cooling of the metallic plate in a bath, animal bloods, foodstuffs, exotic lubricants and colloidal and suspension solutions. For these fluids, the classical Navier-Stokes theory is inadequate. Because of their complexity, there are several models of non-Newtonian fluids in the literature. One of the most popular models for non-Newtonian fluids is the model that is called second-grade fluid [1]. Although there are some criticisms regarding the applications of this model [2], it has been shown by Walters [3] that, for many types of problems in which the flow is slow enough in the viscoelastic sense, the results given using Oldroyd fluid will be substantially similar to those obtained for second grade fluid. Thus, if this is the manner of interpretation of the results, it is reasonable to use the second-grade fluid [4–6] to carry out the calculations. This is particularly so because of the fact that the calculations are generally simpler. This is true not only for exact analytic solutions but even for numerical solutions. The second-grade fluid is the simplest subclass of non-Newtonian fluids for which one can reasonably hope to obtain exact analytic solutions. Moreover, the exact analytic solutions are very important for several reasons. They provide a standard for checking the accuracies of many approximate solutions which can be numerical or empirical. These exact solutions can also be used as tests for verifying numerical schemes that are developed for studying more complex flow problems. Therefore, various researchers [7–9] are engaged in obtaining exact solutions.Recently, the fractional calculus has encountered much success in the description of complex dynamics. In particular, it has been proved to be a valuable tool for handling viscoelastic properties. The starting point of the fractional derivative model of non-Newtonian fluids is usually a classical differential equation which is modified by replacing the time derivative of an integer order by the so-called Riemann-Liouville fractional operator. This generalization allows one to define precisely noninteger order derivatives. The fractional calculus has been found to be quite flexible in describing viscoelastic behavior of fluids. In many different situations fractional calculus has been used to handle various rheological problems [10–21].The aim of this note is to provide exact solutions for the flow of a generalized second grade fluid in the annular region between two infinite coaxial circular cylinders. The motion is produced by the inner cylinder that applies a time-dependent couple to the fluid. More exactly, we would like to extend the results from [7, Section 5] to a larger class of fluids and to a time-dependent couple on the boundary. The general solutions, obtained by means of the integral transforms, will be easily specialized to give the similar solutions for Newtonian and ordinary second grade fluids. Finally, the influence of the pertinent parameters on the fluid motion, as well as a comparison between models, will be underlined by graphical illustrations.


## 2. Governing Equations

The flows to be here considered have the velocity field of the form [22, 23]
(2.1)v=v(r,t)=w(r,t)eθ,
where eθ is the unit vector along the θ-direction of the cylindrical coordinate system r, θ, and z. For such flows the constraint of incompressibility is automatically satisfied. The nontrivial shear stress τ(r,t)=Srθ(r,t) corresponding to such a motion of a second grade fluid is given by [24]
(2.2)τ(r,t)=(μ+α1∂∂t)[∂w(r,t)∂r-w(r,t)r],
where μ is the viscosity and α1 ia a material modulus. In the absence of a pressure gradient in the flow direction and neglecting the body forces, the balance of the linear momentum leads to the relevant equation [25, 26]
(2.3)ρ∂w(r,t)∂t=(∂∂r+2r)τ(r,t).Eliminatingτ(r,t) between (2.2) and (2.3), we get the governing equation
(2.4)∂w(r,t)∂t=(ν+α∂∂t)(∂2∂r2+1r∂∂r-1r2)w(r,t),
where ν=μ/ρ is the kinematic viscosity of the fluid, ρ is its constant density and α=α1/ρ.Generally, governing equations for generalized fluids with fractional derivatives are derived from those of the ordinary fluids by replacing the inner time derivatives of an integer order with the so-called Riemann-Liouville operator [11, 27]
(2.5)Dtβf(t)=1Γ(1-β)ddt∫0t‍f(τ)(t-τ)βdτ,0≤β<1,
where Γ(·) is the Gamma function.Consequently, the governing equations corresponding to the motion (2.1) of a generalized second grade fluid are (cf. [22, Equations (2) and (4)])
(2.6)∂w(r,t)∂t=(ν+αDtβ)(∂2∂r2+1r∂∂r-1r2)w(r,t);τ(r,t)=(μ+α1Dtβ)(∂∂r-1r)w(r,t),
where the new material constant α1 (for simplicity, we are keeping the same notation) goes to the initial α1 for β→1.In this paper, we are interested into the motion of a generalized second grade fluid whose governing equations are given by (2.6). The fractional partial differential equations (2.6), with adequate initial and boundary conditions, can be solved in principle by several methods, the integral transforms technique representing a systematic, efficient, and powerful tool. The Laplace transform will be used to eliminate the time variable and the finite Hankel transform to remove the spatial variable. However, in order to avoid the lengthy calculations of residues and contour integrals, the discrete inverse Laplace transform will be used.


## 3. Rotational Flow between Two Infinite Cylinders

Consider an incompressible generalized second grade fluid at rest in the annular region between two infinitely long coaxial cylinders. At timet=0+, let the inner cylinder of radius R1 be set in rotation about its axis by a time-dependent torque per unit length 2πR1ft and let the outer cylinder of radius R2 be held stationary. Owing to the shear, the fluid between cylinders is gradually moved, its velocity being of the form (2.1). The governing equations are given by (2.6) and the appropriate initial and boundary conditions are (see also [7, Equations (5.2) and (5.3)])
(3.1)w(r,0)=0;r∈[R1,R2],(3.2)τ(R1,t)=(μ+α1Dtβ)(∂w(r,t)∂r-w(r,t)r)|r=R1=ft,w(R2,t)=0;t≥0,
where f is a constant.


### 3.1. Calculation of the Velocity Field

Applying the Laplace transform to (2.6)1 and (3.2), we get
(3.3)qw¯(r,q)=(ν+αqβ)(∂2∂r2+1r∂∂r-1r2)w¯(r,q),(3.4)τ¯(R1,q)=(μ+α1qβ)(∂∂r-1r)w¯(r,q)|r=R1=fq2;w¯(R2,q)=0,
where w¯(r,q) and τ¯(R1,q) are the Laplace transforms of the functions w(r,t) and τ(R1,t), respectively.We denote by [22, Equation (34)]
(3.5)w¯H(rn,q)=∫R1R2‍rw¯(r,q)B(rrn)dr,
the finite Hankel transform of the function w¯(r,q), where
(3.6)B(rrn)=J1(rrn)Y2(R1rn)-J2(R1rn)Y1(rrn),
where rn are the positive roots of the equation B(R2r)=0, and Jp(·), Yp(·) are the Bessel functions of the first and second kind of order p.The inverse Hankel transform ofw¯H(rn,q) is given by [22, Equation (35)]
(3.7)w¯(r,q)=π22∑n=1∞‍rn2J12(R2rn)B(rrn)J22(R1rn)-J12(R2rn)w¯H(rn,q).By means of (3.4)2 and of the identity
(3.8)J1(z)Y2(z)-J2(z)Y1(z)=-2πz,
we can easily prove that
(3.9)∫R1R2‍r(∂2∂r2+1r∂∂r-1r2)w¯(r,q)B(rrn)dr=-rn2w¯H(rn,q)+2πrn(∂∂r-1r)w¯(r,q)|r=R1.Combining (3.3), (3.4), and (3.9), we find that
(3.10)w¯H(rn,q)=2fπrn1q21ρq+α1qβrn2+μrn2.
Writing w¯H(rn,q) under the equivalent form
(3.11)w¯H(rn,q)=2fμπrn3[1q2-1+αrn2qβ-1q{q+(ν+αqβ)rn2}]=2fμπrn3[1q2-q-β-1+αrn2q-2(q1-β+αrn2)+νrn2q-β],
and applying the inverse Hankel transform and using the identities
(3.12)1(q1-β+αrn2)+νrn2q-β=∑k=0∞‍(-νrn2)kq-βk(q1-β+αrn2)k+1,(3.13)π∑n=1∞‍J12(R2rn)B(rrn)rn[J22(R1rn)-J12(R2rn)]=12(R1R2)2(r-R22r),
we find that
(3.14)w¯(r,q)=f2μ(R1R2)2(r-R22r)1q2-πfμ∑n=1∞‍J12(R2rn)B(rrn)rn[J22(R1rn)-J12(R2rn)]×∑k=0∞‍(-νrn2)k(q1-β+αrn2)k+1[q-βk-β-1+αrn2q-βk-2].Now applying the inverse Laplace transform to (3.14), we find for the velocity field the expression
(3.15)w(r,t)=f2μ(R1R2)2(r-R22r)t-πfμ∑n=1∞‍J12(R2rn)B(rrn)rn[J22(R1rn)-J12(R2rn)]×∑k=o∞‍(-νrn2)k[G1-β,-βk-β-1,k+1(-αrn2,t)+αrn2G1-β,-βk-2,k+1(-αrn2,t)],
where the generalized function Ga,b,c(d,t) is defined by [28, Equations (97) and (101)]
(3.16)Ga,b,c(d,t)=L-1{qb(qa-d)c}=∑j=0∞‍djΓ(c+j)Γ(c)Γ(j+1)t(c+j)a-b-1Γ[(c+j)a-b];Re(ac-b)>0,|dqa|<1.


### 3.2. Calculation of the Shear Stress

Applying the Laplace transform to (2.6)2, we find that
(3.17)τ¯(r,q)=(μ+α1qβ)(∂∂r-1r)w¯(r,q).
In order to get a suitable form for τ(r,t), we rewrite (3.10) under the equivalent form
(3.18)w¯H(rn,q)=2fπrn31q2(μ+α1qβ)-2fπrn31q(μ+α1qβ)(q+αqβrn2+νrn2).Applying the inverse Hankel transform to (3.18) and using (3.7) and the identity (3.13), we find that
(3.19)w¯(r,q)=f2(R1R2)2(r-R22r)1q2(μ+α1qβ)-πf∑n=1∞‍J12(R2rn)B(rrn)rn[J22(R1rn)-J12(R2rn)]1q(μ+α1qβ)(q+αqβrn2+νrn2).
Introducing (3.19) into (3.17), it results that
(3.20)τ¯(r,q)=(R1r)2f1q2+πf∑n=1∞‍J12(R2rn)B1(rrn)J22(R1rn)-J12(R2rn)1q(q+αqβrn2+νrn2),
or equivalently (see also (3.12))
(3.21)τ¯(r,q)=(R1r)2f1q2+πf∑n=1∞‍J12(R2rn)B1(rrn)J22(R1rn)-J12(R2rn)∑k=0∞‍(-νrn2)kq-βk-β-1(q1-β+αrn2)k+1,
where B1(rrn)=J2(rrn)Y2(R1rn)-J2(R1rn)Y2(rrn).Now taking the inverse Laplace transform of both sides of (3.21), we get
(3.22)τ(r,t)=(R1r)2ft+πf∑n=1∞‍J12(R2rn)B1(rrn)J22(R1rn)-J12(R2rn)∑k=0∞‍(-νrn2)kG1-β,-βk-β-1,k+1(-αrn2,t).


## 3.1. Calculation of the Velocity Field

Applying the Laplace transform to (2.6)1 and (3.2), we get
(3.3)qw¯(r,q)=(ν+αqβ)(∂2∂r2+1r∂∂r-1r2)w¯(r,q),(3.4)τ¯(R1,q)=(μ+α1qβ)(∂∂r-1r)w¯(r,q)|r=R1=fq2;w¯(R2,q)=0,
where w¯(r,q) and τ¯(R1,q) are the Laplace transforms of the functions w(r,t) and τ(R1,t), respectively.We denote by [22, Equation (34)]
(3.5)w¯H(rn,q)=∫R1R2‍rw¯(r,q)B(rrn)dr,
the finite Hankel transform of the function w¯(r,q), where
(3.6)B(rrn)=J1(rrn)Y2(R1rn)-J2(R1rn)Y1(rrn),
where rn are the positive roots of the equation B(R2r)=0, and Jp(·), Yp(·) are the Bessel functions of the first and second kind of order p.The inverse Hankel transform ofw¯H(rn,q) is given by [22, Equation (35)]
(3.7)w¯(r,q)=π22∑n=1∞‍rn2J12(R2rn)B(rrn)J22(R1rn)-J12(R2rn)w¯H(rn,q).By means of (3.4)2 and of the identity
(3.8)J1(z)Y2(z)-J2(z)Y1(z)=-2πz,
we can easily prove that
(3.9)∫R1R2‍r(∂2∂r2+1r∂∂r-1r2)w¯(r,q)B(rrn)dr=-rn2w¯H(rn,q)+2πrn(∂∂r-1r)w¯(r,q)|r=R1.Combining (3.3), (3.4), and (3.9), we find that
(3.10)w¯H(rn,q)=2fπrn1q21ρq+α1qβrn2+μrn2.
Writing w¯H(rn,q) under the equivalent form
(3.11)w¯H(rn,q)=2fμπrn3[1q2-1+αrn2qβ-1q{q+(ν+αqβ)rn2}]=2fμπrn3[1q2-q-β-1+αrn2q-2(q1-β+αrn2)+νrn2q-β],
and applying the inverse Hankel transform and using the identities
(3.12)1(q1-β+αrn2)+νrn2q-β=∑k=0∞‍(-νrn2)kq-βk(q1-β+αrn2)k+1,(3.13)π∑n=1∞‍J12(R2rn)B(rrn)rn[J22(R1rn)-J12(R2rn)]=12(R1R2)2(r-R22r),
we find that
(3.14)w¯(r,q)=f2μ(R1R2)2(r-R22r)1q2-πfμ∑n=1∞‍J12(R2rn)B(rrn)rn[J22(R1rn)-J12(R2rn)]×∑k=0∞‍(-νrn2)k(q1-β+αrn2)k+1[q-βk-β-1+αrn2q-βk-2].Now applying the inverse Laplace transform to (3.14), we find for the velocity field the expression
(3.15)w(r,t)=f2μ(R1R2)2(r-R22r)t-πfμ∑n=1∞‍J12(R2rn)B(rrn)rn[J22(R1rn)-J12(R2rn)]×∑k=o∞‍(-νrn2)k[G1-β,-βk-β-1,k+1(-αrn2,t)+αrn2G1-β,-βk-2,k+1(-αrn2,t)],
where the generalized function Ga,b,c(d,t) is defined by [28, Equations (97) and (101)]
(3.16)Ga,b,c(d,t)=L-1{qb(qa-d)c}=∑j=0∞‍djΓ(c+j)Γ(c)Γ(j+1)t(c+j)a-b-1Γ[(c+j)a-b];Re(ac-b)>0,|dqa|<1.


## 3.2. Calculation of the Shear Stress

Applying the Laplace transform to (2.6)2, we find that
(3.17)τ¯(r,q)=(μ+α1qβ)(∂∂r-1r)w¯(r,q).
In order to get a suitable form for τ(r,t), we rewrite (3.10) under the equivalent form
(3.18)w¯H(rn,q)=2fπrn31q2(μ+α1qβ)-2fπrn31q(μ+α1qβ)(q+αqβrn2+νrn2).Applying the inverse Hankel transform to (3.18) and using (3.7) and the identity (3.13), we find that
(3.19)w¯(r,q)=f2(R1R2)2(r-R22r)1q2(μ+α1qβ)-πf∑n=1∞‍J12(R2rn)B(rrn)rn[J22(R1rn)-J12(R2rn)]1q(μ+α1qβ)(q+αqβrn2+νrn2).
Introducing (3.19) into (3.17), it results that
(3.20)τ¯(r,q)=(R1r)2f1q2+πf∑n=1∞‍J12(R2rn)B1(rrn)J22(R1rn)-J12(R2rn)1q(q+αqβrn2+νrn2),
or equivalently (see also (3.12))
(3.21)τ¯(r,q)=(R1r)2f1q2+πf∑n=1∞‍J12(R2rn)B1(rrn)J22(R1rn)-J12(R2rn)∑k=0∞‍(-νrn2)kq-βk-β-1(q1-β+αrn2)k+1,
where B1(rrn)=J2(rrn)Y2(R1rn)-J2(R1rn)Y2(rrn).Now taking the inverse Laplace transform of both sides of (3.21), we get
(3.22)τ(r,t)=(R1r)2ft+πf∑n=1∞‍J12(R2rn)B1(rrn)J22(R1rn)-J12(R2rn)∑k=0∞‍(-νrn2)kG1-β,-βk-β-1,k+1(-αrn2,t).


## 4. The Special Caseβ→1

Makingβ→1 into (3.15) and (3.22), we obtain the similar solutions
(4.1)wSG(r,t)=f2μ(R1R2)2(r-R22r)t-πfμ∑n=1∞‍J12(R2rn)B(rrn)rn[J22(R1rn)-J12(R2rn)]×(1+αrn2)∑k=0∞‍(-νrn2)kG0,-(k+2),k+1(-αrn2,t),τSG(r,t)=(R1r)2ft+πf∑n=1∞‍J12(R2rn)B1(rrn)J22(R1rn)-J12(R2rn)∑k=0∞‍(-νrn2)kG0,-(k+2),k+1(-αrn2,t),
corresponding to a second grade fluid performing the same motion.Now, in view of the identity(4.2)∑k=0∞‍(-νrn2)kG0,-(k+2),k+1(-αrn2,t)=1νrn2[1-exp⁡(-νrn2t1+αrn2)],
equation (4.1) can be written under the simplified forms
(4.3)wSG(r,t)=f2μ(R1R2)2(r-R22r)t-πfμν∑n=1∞‍J12(R2rn)B(rrn)rn3[J22(R1rn)-J12(R2rn)]×(1+αrn2){1-exp⁡(-νrn2t1+αrn2)},(4.4)τSG(r,t)=(R1r)2ft+πfν∑n=1∞‍J12(R2rn)B1(rrn)rn2[J22(R1rn)-J12(R2rn)][1-exp⁡(-νrn2t1+αrn2)].
The velocity field can be also processed to give the equivalent form
(4.5)wSG(r,t)=f2μ(R1R2)2(r-R22r)(t-α1μ)-πfμν∑n=1∞‍J12(R2rn)B(rrn)rn3[J22(R1rn)-J12(R2rn)]×{1-(1+αrn2)exp⁡(-νrn2t1+αrn2)}.Makingα1 and then α→0 into (4.3) and (4.4), the velocity field
(4.6)wN(r,t)=f2μ(R1R2)2(r-R22r)t-πfμν∑n=1∞‍J12(R2rn)B(rrn)rn3[J22(R1rn)-J12(R2rn)][1-exp⁡(-νrn2t)]
and the associated shear stress
(4.7)τN(r,t)=(R1r)2ft+πfν∑n=1∞‍J12(R2rn)B1(rrn)rn2[J22(R1rn)-J12(R2rn)][1-exp⁡(-νrn2t)],
corresponding to a Newtonian fluid, are obtained.


## 5. Conclusions

The purpose of this note is to provide exact analytic solutions for the velocity fieldw(r,t) and the shear stress τ(r,t) corresponding to the unsteady rotational flow of a generalized second grade fluid between two infinite coaxial cylinders, the inner cylinder being set in rotation about its axis by a constantly accelerating shear. The solutions that have been obtained, presented under series form in terms of usual Bessel (J1(·), J2(·), Y1(·), and Y2(·)) and generalized Ga,b,c(·,t) functions, satisfy all imposed initial and boundary conditions. They can be easily specialized to give the similar solutions for ordinary second grade and Newtonian fluids. Furthermore, in view of some recent results [29, Equation (3.15)], our velocity field (4.3) is in accordance with that obtained in [7, Equation (5.17)] by a different technique.Now, in order to reveal some relevant physical aspects of the obtained results, the diagrams of the velocityw(r,t) and of the shear stress τ(r,t) are depicted against r for different values of the time t and of the pertinent parameters. From Figures 1(a) and 1(b), containing the diagrams of the velocity and the shear stress at several times, it clearly results in the influence of the rigid boundary on the fluid motion. The velocity is an increasing function of t. For the same values of the parameters, the shear stress, in absolute value, is also an increasing function of t. The influence of the kinematic viscosity ν on the fluid motion is shown in Figures 2(a) and 2(b). The velocity is a decreasing function of ν. The shear stress, in absolute value, on the first part of the flow domain, near the moving cylinder, is a decreasing function of ν. It is an increasing function of ν in the neighborhood of the stationary cylinder. Figure 3 shows the influence of the parameter α on the flow motion. Both the velocity and the shear stress, in absolute value, on the first part of the flow domain, near the moving cylinder, are increasing functions of α. They are decreasing functions of α in the neighborhood of the stationary cylinder. The influence of the fractional parameter β on the fluid motion is shown in Figure 4. Its effect on the fluid motion is qualitatively opposite to that of parameter α.Profiles of the velocityω(r,t) and the shear stress τ(r,t) given by (3.15) and (3.22), for f=-1,  R1=0.3, R2=0.5, ν=0.000259, μ=0.246, α=0.00153, β=0.3, and different values of t.
(a)
(b)Profiles of the velocityω(r,t) and the shear stress τ(r,t) given by (3.15) and (3.22), for f=-1, R1=0.3, R2=0.5, α=0.00153, β=0.3  t=10 s and different values of ν.
(a)
(b)Profiles of the velocityω(r,t) and the shear stress τ(r,t) given by (3.15) and (3.22), for f=-1, R1=0.3, R2=0.5, ν=0.000259, μ=0.246, β=0.3 and different values of α.
(a)
(b)Profiles of the velocityω(r,t) and the shear stress τ(r,t) given by (3.15) and (3.22), for f=-1, R1=0.3, R2=0.5, ν=0.000259, μ=0.246, α=0.00153 and different values of β.
(a)
(b)Finally, for comparison, the profiles ofw(r,t) and τ(r,t) corresponding to the motion of the three models (Newtonian, ordinary second grade, and generalized second grade) are together depicted in Figure 5, for the same values of t and of the common material parameters. In all the cases the velocity of the fluid is a decreasing function with respect to r, and the Newtonian fluid is the swiftest, while the generalized second grade fluid is the slowest in the region near the moving cylinder. The units of the material constants are SI units within all figures, and the roots rn have been approximated by 2(n-1)π/[2(R2-R1)].Profiles of the velocityω(r,t) and the shear stress τ(r,t) corresponding to Newtonian, second grade and generalized second grade fluids, for f=-1, R1=0.3, R2=0.5, ν=0.000259, μ=0.246, α=0.00153  β=0.3 and t=10 s.
(a)
(b)



---

*Source: 374670-2012-08-09.xml*