Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina

IMPACT FACTOR 2018: 1.005

CiteScore 2018: 1.01

SCImago Journal Rank (SJR) 2018: 0.237
Source Normalized Impact per Paper (SNIP) 2018: 0.541

ICV 2017: 162.45

Open Access
See all formats and pricing
More options …
Volume 16, Issue 1


Volume 13 (2015)

Unsteady flow of fractional Oldroyd-B fluids through rotating annulus

Madeeha Tahir / Muhammad Nawaz Naeem / Maria Javaid / Muhammad Younas / Muhammad Imran / Naeem Sadiq / Rabia Safdar
Published Online: 2018-04-18 | DOI: https://doi.org/10.1515/phys-2018-0028


In this paper exact solutions corresponding to the rotational flow of a fractional Oldroyd-B fluid, in an annulus, are determined by applying integral transforms. The fluid starts moving after t = 0+ when pipes start rotating about their axis. The final solutions are presented in the form of usual Bessel and hypergeometric functions, true for initial and boundary conditions. The limiting cases for the solutions for ordinary Oldroyd-B, fractional Maxwell and Maxwell and Newtonian fluids are obtained. Moreover, the solution is obtained for the fluid when one pipe is rotating and the other one is at rest. At the end of this paper some characteristics of fluid motion, the effect of the physical parameters on the flow and a correlation between different fluid models are discussed. Finally, graphical representations confirm the above affirmation.

Keywords: Fractional Oldroyd-B model; fractional differential operator; annulus; rotational flow; velocity field; shear stress; integral transforms

PACS: 47.10.A-; 47.27.nf

1 Introduction

Humans benefit from fluid flow but also need to understand the negative effects in nature, e.g. hurricanes and cyclones. The topographic conditions and the sufficient supply of rain flow is important in the cultivation of crops. Cultivation in deserts can also be achieved with artificial irrigation. Artificial breathing machines, artificial hearts and dialysis systems are designed using fluid dynamics. Hence fluids are vital.

The stress and rate of strain of a fluid is related: fluid is Newtonian if this relationship is linear, if not, the fluid is non-Newtonian. The non-Newtonian fluids, like polymers, plasma, paste and pulps are very important in industry. Therefore, the discussion of non-Newtonian fluids in fluid mechanics is important. All existing fluids cannot be described by the Navier-stocks equations. By categorizing different fluids, we have different mathematical models to represent fluids. Oldroyd-B fluid is one of these rate-type models.

The non-Newtonian fluids possessing shear dependent viscosity are known as viscoelastic fluids. The first rate type fluid model is a Maxwell model that is still in used. Fetecau et. al [1] derived the results of the velocity field and shear stress of an Oldroyd-B fluid on a plate, accelerating constantly. Hayat et. al [2] determined the solutions for MHD flow over an infinite oscillatory plate about an axis that is normal to the oscillatory plate. Some authors determined solutions for a viscoelastic fluid model by reproducing a kernel method (RKM) [3,4,5,6]. The unsteady incompressible Oldroyd-B fluid is analyzed by Jamil et. al [7] and special cases are also obtained for Maxwell, second grade and Newtonian fluids. Similarly, Burdujan [8] used this idea on Taylor-Couette flow and calculated the velocity field and shear stress of a rotating fluid between annuls of double cylinders. Numerical solutions of fractional fluids are also under attention [9, 10].

A lot of interest is given to a rotating circular domain corresponding to helical flow for second grade, Maxwell fluids and Oldroyd-B [11,12,13,14]. The flow of a polymer in a circular cylinder by considering pulsatile APG was analyzed by Barnes et. al [15, 16]. Davies et. al [17] and Phan-Thien [18] presented the same problem for a White-Metzner fluid. To solve the viscoelastic motion of different fluid models, one needs to learn about fractional calculus [19, 20]. Fractional calculus is very useful in fields of mathematics and physics [21,22,23]. Furthermore, the interested reader can consult the references [24,25,26,27,28,29,30].

In this proposal, we consider the fractional Oldroyd-B fluid model and study the flow due to the circulation of two pipes around its axial. Integral transforms are used to obtain the general solutions and the final result is in the form of hypergeometric functions. There are some limiting cases by applying the limits on physical parameters, i.e ς →1, λr →0 and λ → 0, solutions for the ordinary Oldroyd-B, fractional and ordinary Maxwell and Newtonian fluids are discussed respectively.

2 Development of the flow

Consider a fractional Oldroyd-B fluid (OBFFD) between two coaxial pipes. The radii of the pipes are R1 and R2(>R1), after time t = 0+, both pipes and the fluid between them start rotating along their axis. For the considered problem, the z-axis is along the axis of the pipes, the r-axis is perpendicular to the axis of the pipe and θ is along the tangent to the boundary. Let vr, vθ and vz be the velocity components in the direction of the r, θ and z axes respectively.

The continuity and linear momentum equations for incompressible flow is



where V, t, ρ, b, d/dt is velocity, time, density, body force and the material time derivative respectively.

For an incompressible Oldroyd-B fluid [36]


where p is the hydrostatic pressure, I denotes the identity tensor, S is the extra-stress tensor, λ and λr are relaxation and retardation times, superposed dot indicates the material time derivative, L represents the velocity gradient, A = L + LT is the first Rivlin Ericsen tensor, μ is the dynamic viscosity of the fluid and the superscript T denotes the transpose operation. For the problem in consideration, the velocity field and the extra stress tensor has the following form [31]


For such flows, (1) is true.

The governing equations for a Oldroyd-B fluid, corresponding to the above defined motions are [31]



where ν is the kinematic viscosity and τ(r, t) is the non zero shear stress.

The governing equations for a fractional Oldroyd-B fluid is obtained by using a fractional differential operator Dtς, defined by [32, 33]


in governing equations (5) and (6) we have



3 Initial and boundary conditions

We need three conditions to solve a fractional partial differential equation (9). The relevant initial and boundary conditions are



where 1, 2 and κ are real numbers.

3.1 Analytic solutions for the velocity field

Using the Laplace transform Eqs. (9) and (11) implies that



where the Laplace transform of the function ω(r, t) is ω(r, q). The finite Hankel transform [34] of the function ω(r, q), can be expressed as




rn are the positive roots of the transcendental equation A(R1, r) = 0, and Jp(·) is Bessel functions of the first kind and Yp(·) is Bessel functions of second kind of order p.

Taking the Hankel transform of the differential equation (12) and keeping in mind the boundary condition (13) and following identities



we get


Writing the above equations in a suitable form


The inverse Hankel transform of ωH(rn, q) [34]


Using the inverse Hankel transform on Eq. (19), and the identity


we obtain


Using the identity


Eq. (21) can be further simplified to give


Now using inverse Laplace transform on Eq. (23) and taking into account (A1) [35] from the Appendix, we find that


3.2 Calculation of the shear stress

Using the Laplace transform, Eq. (5) implies that


Using Eqs. (21) in the above equation we get




Now, using the inverse Laplace transform on Eq. (26), finally the shear stress


4 Limiting cases

  1. Letting ς → 1 into Eqs. ((24)) and ((28)) the velocity field


    and the shear stress


    for the ordinary Oldroyd-B fluid.

  2. Making λr →0 in Eqs. ((24)) and ((28)), the velocity field


    and the shear stress


    corresponding to a fractional Maxwell fluid.

  3. Making ς → 1 in Eqs. ((31)) and ((32)), the velocity field is


    and the shear stress is


    for the ordinary Maxwell fluid.

    Using Appendix A2, the expressions (33) and (34) in the simplified form for κ = 1, in terms of expositional functions, are




  4. By now letting λ→ 0 in Eqs. (35) and (36) or ς → 1 and λ→ 0 in Eqs. (31) and (32) and using limλ01λkG1,b,k1λ,t=tb1Γ(b); b < 0, the velocity field is


    and the shear stress is


    corresponding to a Newtonian fluid.

5 Flow through a circular pipe

The finite Hankel transform of the function f(r, t) is


and the inverse Hankel transform of fH (rn) is


where J1(R2r) = 0 as rn are the positive roots of the equation.

By making J1(R2rn) = 0, and using the identities J0(z) + J2(z) = 2zJ1 (z) in Eq. (15) we get


Writing the Integral


Introducing Eq. (41) and taking R1 = 0, the above integral takes the form


as a result from (24) and (28) we recovered the expansions for the velocity field




through a cylinder already obtained by Kamran et al. [31] in Equation (21) and (27) for p = 1.

6 Conclusions and results

In this article, exact solutions for a fractional Oldroyd-B fluid between two rotating pipes are determined. Fluid motion is produced due to rotation of the pipes around their axis with time dependent angular velocities. The solutions determined by the use of integral transforms and presented in terms of Bessel functions and hypergeometric functions which are free of integrals. The final result satisfies the initial and boundary conditions. In the limiting cases, the corresponding results for Oldroyd-B fluid, Maxwell and Newtonian fluid, are obtained from general results. Moreover, the solution for the fluid through a circular pipe is obtained as a special case.

The velocity profiles and shear stress are illustrated in Figures 27 for different values of parameters. Fig. 2 for the effect of t, Fig. 3 for the effect of ν, Fig. 4 for the effect of λ, Fig. 5 for the effect of λr and Fig. 6 shows the effect of the fractional parameter ς on the fluid motion. Clearly, the velocity as well as shear stress are increasing as a function of t and ν, The effect is qualitatively the same for fractional parameter ς and λ, more exactly the velocity ω(r, t) is a decreasing function with regards to ς.

Problem Geometry.
Figure 1

Problem Geometry.

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for R1 = 1, R2 = 7.4, ℧1 = 1, ℧2 = −1, λ = 4, λr = 1, ν = 0.002, ς = 0.1 and different values of time
Figure 2

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for R1 = 1, R2 = 7.4, ℧1 = 1, ℧2 = −1, λ = 4, λr = 1, ν = 0.002, ς = 0.1 and different values of time

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for t = 7, R1 = 1, R2 = 7.4, ℧1 = 1,℧2 = −1, λ = 20, λr = 5, ς = 0.9 and different values of ν
Figure 3

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for t = 7, R1 = 1, R2 = 7.4, ℧1 = 1,℧2 = −1, λ = 20, λr = 5, ς = 0.9 and different values of ν

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for t = 7, R1 = 1, R2 = 7.4,℧1 = 1, ℧2 = −1, ν = 0.2,λr = 1, ς = 0.5 and different values of λ
Figure 4

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for t = 7, R1 = 1, R2 = 7.4,℧1 = 1, ℧2 = −1, ν = 0.2,λr = 1, ς = 0.5 and different values of λ

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for t = 8, R1 = 1, R2 = 7.4, ℧1 = 1, ℧2 = −1,ν = 0.2,λ = 10,ς = 0.5 and different values of λr
Figure 5

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for t = 8, R1 = 1, R2 = 7.4, ℧1 = 1, ℧2 = −1,ν = 0.2,λ = 10,ς = 0.5 and different values of λr

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for t = 20, R1 = 1, R2 = 7.4,℧1 = 1, ℧2 = −1, ν = 0.001, λ = 4,λr = 1 and different values of ς
Figure 6

Profiles of the velocity field ω(r, t) given by Eqs. (24) and (28) for t = 20, R1 = 1, R2 = 7.4,℧1 = 1, ℧2 = −1, ν = 0.001, λ = 4,λr = 1 and different values of ς

Finally, Fig. 7 is the comparison of ω(r, t) corresponding to the motion of the Newtonian, fractional and ordinary Maxwell, and fractional and ordinary Oldroyd-B for the same values of time and for the same material and fractional parameters. From this comparison, the Oldroyd-B fluid is the slowest and the Newtonian fluid is the fastest on the whole flow domain. SI units are used for all figures, and the roots rn=nπ(R2R1)..

Profiles of the velocity field ω(r, t) corresponding to the Newtonian, Maxwell, Fractional Maxwell, Oldroyd-B and fractional Oldroyd-B fluids, for t = 8, R1 = 1, R2 = 7.4, ℧1 = 1, ℧2 = −1,ν = 0.001,λ = 3.6,λr = 3 and ς = 0.8
Figure 7

Profiles of the velocity field ω(r, t) corresponding to the Newtonian, Maxwell, Fractional Maxwell, Oldroyd-B and fractional Oldroyd-B fluids, for t = 8, R1 = 1, R2 = 7.4, ℧1 = 1, ℧2 = −1,ν = 0.001,λ = 3.6,λr = 3 and ς = 0.8







  • [1]

    Fetecau C., Prasad S.C., Rajagopal K.R., A note on the flow induced by a constantly accelerating plate in an Oldroyd-B fluid, Applied mathematical modelling, 2007, 31(4), 647-54. Web of ScienceCrossrefGoogle Scholar

  • [2]

    Hayat T., Hutter K., Asghar S., Siddiqui A.M., MHD flows of an Oldroyd-B fluid, Mathematical and Computer Modelling, 2002, 36(9-10), 987-95. CrossrefGoogle Scholar

  • [3]

    Akgul A., Khan Y., Akgul E.K., Baleanu D., Qurashi M.M., Solutions of nonlinear systems by reproducing kernel method, Journal of Nonlinear Sciences and Applications, 2017, 10(8), 4408-4417. CrossrefWeb of ScienceGoogle Scholar

  • [4]

    Akgül A., A new method for approximate solutions of fractional order boundary value problems. Neural, parallel & scientific computations, 2014, 22(1-2), 223-237. Google Scholar

  • [5]

    Akgül A., Baleanu D., Inc M., Tchier F., On the solutions of electrohydrodynamic flow with fractional differential equations by reproducing kernel method, Open Physics, 2016, 14(1), 685-689. Web of ScienceGoogle Scholar

  • [6]

    Akgül A., Inc M., Karatas E., Reproducing kernel functions for difference equations, Discrete & Continuous Dynamical Systems-S, 2015, 8, 6, 1055-1064. CrossrefWeb of ScienceGoogle Scholar

  • [7]

    Jamil M., Fetecau C., Rana M., Some exact solutions for Oldroyd-B fluid due to time dependent prescribed shear stress, Journal of Theoretical and Applied Mechanics, 2012, 50(2), 549-562. Google Scholar

  • [8]

    Burdujan L., The flow of a particular class of Oldroyd-B fluids, J. Series on Math and its Applications, 2011, 3, 23-45. Google Scholar

  • [9]

    Hashemi M.S., Inc M., Karatas E., Akgül A., A numerical investigation on burgers equation by mol-gps method, Journal of Advanced Physics, 2017, 6(3), 413-417. CrossrefWeb of ScienceGoogle Scholar

  • [10]

    Akgül A., Karatas E., Baleanu D., Numerical solutions of fractional differential equations of Lane-Emden type by an accurate technique, Advances in Difference Equations, 2015, 2015(1), 220. Web of ScienceCrossrefGoogle Scholar

  • [11]

    Fetecău C., Fetecau C., On the uniqueness of some helical flows of a second grade fluid, Acta mechanica, 1985, 57(3), 247-52. CrossrefGoogle Scholar

  • [12]

    Fetecau C., Fetecau C., Vieru D., On some helical flows of Oldroyd- B fluids, Acta Mechanica, 2007, 189(1-4), 53-63. Web of ScienceCrossrefGoogle Scholar

  • [13]

    Vieru D., Akhtar W., Fetecau C., Fetecau C., Starting solutions for the oscillating motion of a Maxwell fluid in cylindrical domains, Meccanica, 2007, 42(6), 573-83. Web of ScienceCrossrefGoogle Scholar

  • [14]

    Nadeem S., Asghar S., Hayat T., Hussain M., The Rayleigh Stokes problem for rectangular pipe in Maxwell and second grade fluid, Meccanica, 2008, 43(5), 495-504. CrossrefWeb of ScienceGoogle Scholar

  • [15]

    Barnes H.A, Townsend P., Walters K., On pulsatile flow of non-Newtonian liquids, Rheologica Acta, 1971, 10(4), 517-527. CrossrefGoogle Scholar

  • [16]

    Barnes H.A., Townsend P., Walters K., Flow of non-Newtonian liquids under a varying pressure gradient, Nature, 1969, 224(5219), 585-587. CrossrefGoogle Scholar

  • [17]

    DaVies J.M., Bhumiratana S., Bird R.B., Elastic and inertial effects in pulsatile flow of polymeric liquids in circular tubes, Journal of Non-Newtonian Fluid Mechanics, 1978, 3(3), 237-259. CrossrefGoogle Scholar

  • [18]

    Phan-Thien N., Dudek J., Pulsating flow revisited, Journal of Non-Newtonian Fluid Mechanics, 1982, 11(1-2), 147-161. Web of ScienceCrossrefGoogle Scholar

  • [19]

    Bagley R.L., Torvik P.J., On the fractional calculus model of viscoelastic behavior, Journal of Rheology, 1986, 30(1), 133-155. CrossrefGoogle Scholar

  • [20]

    Friedrich C.H., Relaxation and retardation functions of the Maxwell model with fractional derivatives, Rheologica Acta, 1991, 30(2), 151-158. CrossrefGoogle Scholar

  • [21]

    Inc M., Kiliç B., Karataş E., Akgăl A., Solitary Wave Solutions for the Sawada-Kotera Equation, Journal of Advanced Physics, 2017, 6(2), 288-293. CrossrefWeb of ScienceGoogle Scholar

  • [22]

    Guariglia E., Silvestrov S., Fractional-Wavelet Analysis of Positive definite Distributions and Wavelets on D0(C). In: Engineering Mathematics II, 2016, 337-353. Springer, Cham. Google Scholar

  • [23]

    Coussot C., Kalyanam S., Yapp R., Insana M.F., Fractional derivative models for ultrasonic characterization of polymer and breast tissue viscoelasticity, IEEE transactions on ultrasonics, ferroelectrics, and frequency control, 2009, 56(4). Web of ScienceGoogle Scholar

  • [24]

    Wenchang T., Wenxiao P., Mingyu X., A note on unsteady flows of a viscoelastic fluid with the fractional Maxwell model between two parallel plates, International Journal of Non-Linear Mechanics, 2003, 38(5), 645-650. CrossrefGoogle Scholar

  • [25]

    Nazar M., Fetecau C., Awan A.U., A note on the unsteady flow of a generalized second-grade fluid through a circular cylinder subject to a time dependent shear stress, Nonlinear Analysis: Real World Applications, 2010, 11(4), 2207-2214. CrossrefWeb of ScienceGoogle Scholar

  • [26]

    Athar M., Kamran M., Fetecau C., Taylor-Couette flow of a generalized second grade fluid due to a constant couple, Nonlinear Analysis: Modelling and Control, 2010, 15(1), 3-13. Google Scholar

  • [27]

    Mahmood A., Fetecau C., Khan N.A., Jamil M., Some exact solutions of the oscillatory motion of a generalized second grade fluid in an annular region of two cylinders, Acta Mechanica Sinica. 2010, 26(4), 541-550. CrossrefWeb of ScienceGoogle Scholar

  • [28]

    Kamran M., Imran M., Athar M., Exact solutions for the unsteady rotational flow of a generalized second grade fluid through a circular cylinder, Nonlinear Analysis: Modelling and Control, 2010, 15(4), 437-44. Google Scholar

  • [29]

    Fetecau C., Awan A.U., Athar M., A note on “Taylor-Couette flow of a generalized second grade fluid due to a constant couple”, Nonlinear Analysis: Modelling and Control, 2010, 15(2), 155-158. Google Scholar

  • [30]

    Athar M., Kamran M., Imran M., On the unsteady rotational flow of a fractional second grade fluid through a circular cylinder, Meccanica, 2012, 47(3), 603-611. Web of ScienceCrossrefGoogle Scholar

  • [31]

    Kamran M., Imran M., Athar M., Exact solutions for the unsteady rotational flow of an Oldroyd-B fluid with fractional derivatives induced by a circular cylinder, Meccanica, 2013, 48(5), 1215-1226. CrossrefWeb of ScienceGoogle Scholar

  • [32]

    Podlubny I., Fractional Differential Equations Academic Press, San Diego, 1999. Google Scholar

  • [33]

    Guariglia E., Fractional Derivative of the Riemann Zeta Function, Fractional Dynamics, Chapter: 21, De Gruyter, Cattani C., Srivastava H., Yang X.J., (2015), 357-368. Google Scholar

  • [34]

    Sneddon I.N., Functional Analysis in: Encyclopedia of Physics, Vol. II, Springer Berlin, 1955 Google Scholar

  • [35]

    Lorenzo C.F., Hartley T.T., Generalized functions for the fractional calculus, National Aeronautics and Space Administration, Glenn Research Center, 1999. Google Scholar

  • [36]

    Hayat T., Siddiqui A.M., Asghar S., Some simple flows of an Oldroyd-B fluid, International Journal of Engineering Science, 2001, 39(2), 135-147. CrossrefGoogle Scholar

About the article

Received: 2017-09-18

Accepted: 2018-01-15

Published Online: 2018-04-18

Funding: This research is supported by the Government College University, Faisalabad, Pakistan and the Higher Education Commission Pakistan.

Competing interestCompeting interests: The authors declare that they have no competing interests.

Citation Information: Open Physics, Volume 16, Issue 1, Pages 193–200, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0028.

Export Citation

© 2018 Madeeha Tahir et al., published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Citing Articles

Here you can find all Crossref-listed publications in which this article is cited. If you would like to receive automatic email messages as soon as this article is cited in other publications, simply activate the “Citation Alert” on the top of this page.

Bo Wang, Madeeha Tahir, Muhammad Imran, Maria Javaid, and Chahn Yong Jung
IEEE Access, 2019, Volume 7, Page 72482

Comments (0)

Please log in or register to comment.
Log in