Analytical Solution for the Flow of a Generalized Oldroyd-B Fluid in a Circular Cylinder

Author(s): Haitao Qi1, Nida Fatima2, Hassan Waqas3, Junaid Saeed3
1School of Mathematics and Statistics, Shandong University at Weihai, Weihai 264209, China.
2Division of Science and Technology, University of Education, Lahore-54590, Pakistan.
3Department of Mathematics, Government College University, Faisalabad, Pakistan.
Copyright © Haitao Qi, Nida Fatima, Hassan Waqas, Junaid Saeed. 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.

Abstract

The tangential stress and velocity field corresponding to the flow of a generalized Oldroyd-B fluid in an infinite circular cylinder will be determined by mean of Laplace and finite Hankel transform. The motion is produced by the cylinder, that after \(t=0^{+}\), begins to rotate about its axis, under the action of oscillating shear stress \(\Omega R \sin(\omega t)\) given on boundary. The solutions are based on an important remark regarding the governing equation for the non- trivial shear stress. The solutions that have been obtained satisfy all imposed initial and boundary conditions. The obtained solution will be presented under series form in term of generalized G-function. The similar solutions for the ordinary Oldroyd-B fluid, Maxwell, ordinary Maxwell and Newtonian fluids performing the same motion will be obtained as special cases of our general solutions.

Keywords: fractional oldroyd-B fluid, cylindrical Domain, unsteady rotating flow, Velocity field, shear stress; laplace and finite Hankel transforms.

1. Introduction

The Oldroyd-B fluid models is very important among the fluids of rate type due to its special behavior. Also, this model contains the Newtonian fluid model and Maxwell fluid model as special cases. The Oldroyd-B fluid model [1, 2] considered the memory effects and elastic effects exhibited by a large class of fluids such as the biological and polymeric liquids. The motion of a fluid in the neighborhood of a moving body is of great interest for industry. The flow between cylinders or through a rotating cylinder has applications in the food industry and being one of the most important and interesting problem of motion near rotating bodies. Exact solutions for some simple flows of Oldroyd-B fluids were presented by many authors, see for example, Rajagopal and Bhatnagar [3], Hayat et al. [4, 5]. The velocity distriution for different motions of Newtonian fluids through a circular cylinder is given in [6]. Wood [7] has considered the general case of helical flow of an Oldroyd-B fluid, due to the combined action of rotating cylinders(with constant angular velocities) and a constant axial pressure gradient. Accurate solutions regarding motions of Non-Newtonian fluids in cylindrical domains appear to be those of Ting [2], Srivastava [8] and Water and King [9] for second grade, Maxwell and Oldroyd-B fluids respectively. The most general solution corresponding to the helical flow of a second grade fluid seem to be those of Fetecau and Cornia Fetecau [10], in which the cylinder is rotating around its axis and sliding along the same axis with time-dependent velocities. There is a vast literature dealing with such fluids, but we shall recall here only a few of the most recent papers [12, 13, 14, 15, 16]. Most existing solutions in the literature correspond to problems with boundary conditions on the velocity. Though, all above mentioned papers incorporate motion problem in which velocity is given on the boundary. In [16], Renardy explained how well posed boundary value problems can be formulated using boundary conditions on stress. Water and King [17] were among the first specialists who used the shear stress on boundary to find exact solution for motions of rate type fluids. Our goal is to investigate analytical solution for the flow of a generalized Oldroyd-B fluid in a circular cylinder. We considered the boundary conditions on the shear stress. The flow of fluid is due to rotation of the cyliner around its axis, under the action of oscillating shear stress \(\Omega R \sin(\omega t)\) given on boundary. These solutions are obtained by mean of integral transforms. The obtained solution satisfy the all imposed initial and boundary conditions. Finally, solution of the ordinary Oldroyd-B fluid, Maxwell, ordinary Maxwell and Newtonian fluid flows are obtained as particular cases of our general results.

2. Mathematical formulation of the problem

For an Oldroyd-B fluid constitutive equations is $$T=-p \textbf{I}+\textbf{S}\,\,\ ; \,\,\ S+\lambda\bigg({\frac{d \textbf{S}}{d t}}-\textbf{LS}-\textbf{S} \textbf{L}^{T}\bigg),$$
\begin{eqnarray} T=\mu{\textbf{A}+\mu\lambda_{r}\bigg(\frac{d\textbf{A}}{dt}-\textbf{LA}-\textbf{AL}^{T}\bigg)}, \end{eqnarray}
(1)
where T (Cauchy stress tensor), -pI (indeterminate spherical stress), S (stress tensor), L (velocity gradient), \(\mu\) (dynamic viscosity), \(\textbf{A}=\textbf{L}+\textbf{L}^{T}\) (first the Rivilin-Erickson tensor,) \(\lambda\) and \(\lambda_{r}\) \((0\leqslant\lambda_{r}< \lambda)\) are relaxation and retardation time. Assume an infinite circular cylinder rotate along z-axis with radius R. Cylinder is filled with an Oldroyd-B fluid which is at rest at time \(t=0\). After \(t=0^{+}\) the cylinder applies an oscillating rotational shear stress \(\Omega R \sin(\omega t)\) to the fluid, where \(\omega\) is the angular frequency. We assumed that velocity field and the extra shear stress are of the form
\begin{eqnarray} \textbf{V}=\textbf{V} (r,t)=W(r,t)e_{\theta} \,\,\,\,\,\,\,\,\ ; \,\,\,\,\,\,\,\ \textbf{S}=\textbf{S}(r,t), \end{eqnarray}
(2)
where \(e_{\theta}\) is unit vector in the \(\theta\)-direction of the cylindrical coordinate system.We assume that \(\textbf{S}\) and \(\textbf{V}\) is a function of time and radius only. At \(t=0\) there is no motion in fluid i.e; fluid is at rest then
\begin{eqnarray} \textbf{V}(r,0)=0 \,\,\,\,\,\,\,\,\,\,\ ; \,\,\,\,\,\,\,\,\,\ \textbf{S}(r,0)=0. \end{eqnarray}
(3)
Introducing Eqs.(2) in(1) and using (3) we get \(S_{rr}=S_{rz}=S_{z\theta}=S_{zz}=0\) and the meaning partial differential equation.
\begin{eqnarray} \bigg(1+\lambda\frac{\partial}{\partial t}\bigg) \tau(r,t) = \mu\bigg(1+\lambda_{r}\frac{\partial}{\partial t}\bigg) \bigg(\frac{\partial}{\partial r}-\frac{1}{r}\bigg) \omega(r,t), \end{eqnarray}
(4)
where \(\tau(r,t)=S_{r\theta}(r,t)\) is non-zero component of extra stress tensor. If we neglect the body force, then due to rotation symmetry the balance of linear momentum leads to the relevant equations.
\begin{eqnarray} \rho\frac{\partial }{\partial t}\omega(r,t) = \bigg(\frac{\partial}{\partial r}+\frac{2}{r}\bigg)\tau(r,t), \end{eqnarray}
(5)
where \(\rho\) is the constant density of the field. In order to obtained the governing equations for shear stress on boundary we eliminate \(\omega(r,t)\) between equation (4) and (5) after elimination \(\omega(r,t)\) we get,
\begin{eqnarray} \bigg(1+\lambda\frac{\partial}{\partial t}\bigg)\frac{\partial }{\partial t}\tau (r,t)=\nu\bigg(1+\lambda_{r}\frac{\partial}{\partial t}\bigg)\bigg(\frac{\partial^{2}}{\partial r^{2}}+\frac{1}{r}\frac{\partial}{\partial r}-\frac{4}{r^{2}}\bigg)\tau(r,t), \end{eqnarray}
(6)
where \(\nu=\frac{\mu}{\rho}\) is the Kinematic viscosity of the fluid. The governing model by using fractional derivative is shown as;
\begin{eqnarray} \bigg(1+ \lambda D_t^{\alpha}\bigg)\tau(r,t)=\mu\bigg(1+ \lambda_{r} D_t^{\beta}\bigg)\bigg(\frac{\partial}{\partial r}-\frac{1}{r}\bigg)\omega(r,t), \end{eqnarray}
(7)
\begin{eqnarray} \rho D_t^{\alpha} \omega(r,t) = \bigg(\frac{\partial}{\partial r}+\frac{2}{r}\bigg)\tau(r,t), \end{eqnarray}
(8)
\begin{eqnarray} \bigg(1+ \lambda D_t^{\alpha}\bigg)\frac{\partial}{\partial t}\tau (r,t)=\nu\bigg(1+ \lambda_{r} D_t^{\beta}\bigg)\bigg(\frac{\partial^{2}}{\partial r^{2}}+\frac{1}{r}\frac{\partial}{\partial r}-\frac{4}{r^{2}}\bigg)\tau(r,t). \end{eqnarray}
(9)
Where $$D_{t}^{\alpha}f(t)=\frac{1}{\Gamma(1-\alpha)}\frac{d}{dt}\int\limits_0^{t}\frac{f(\tau)}{(t-\tau)^{\alpha}}d\tau \,\,\,\ ; \,\,\,\ 0\leqslant\alpha< 1$$ is Caputo fractional derivative operator and $\Gamma(.)$ is the Euler integral of second kind or Gamma function.
At \(t=0\) fluid is at rest because there is no rotation in the cylinder, at \(t=0^{+}\) cylinder starts its rotation along its axis and boundary of cylinder applies shear stress on fluid and radius of cylinder is \(R\). Appropriate conditions i.e; initial and boundary conditions are,
\begin{eqnarray} \tau(r,t)|_{t=0}=\frac{\partial \tau(r,t)}{\partial t}|_{t=0} =0; \,\,\,\,\,\,\,\,\,\,\, \,\,\,\,\,\,\,\,\, r\epsilon[0,R] \end{eqnarray}
(10)
\begin{eqnarray} \tau(R,t)=\Omega R \sin(\omega t),: \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\, t\geqslant0, \end{eqnarray}
(11)
where \(\Omega\) is constant.

3. Calculation Of Shear Stress

We shall use the Laplace transform and Finite Hankel transform to determine the exact analytical solution. Taking the Laplace transform of the equations (9) and (11) we have
\begin{eqnarray} q(1+ \lambda q^{\alpha})\overline{\tau}(r,q)=\nu(1+ \lambda_{r}q^{\beta})\bigg(\frac{\partial^{2}}{\partial r^{2}}+\frac{1}{r}\frac{\partial}{\partial r}-\frac{4}{r^{2}}\bigg)\overline{\tau}(r,q), \end{eqnarray}
(12)
\begin{eqnarray} \overline{\tau}(R,q)=\Omega R \frac{\omega}{q^{2}+\omega^{2}}, \end{eqnarray}
(13)
where \(\overline{\tau}(r,q)\) represent the Laplace transform of the function \(\tau(r,t)\). We can write equation (12) as,
\begin{eqnarray} \overline{\tau}(r,q)=\frac{\nu(1+\lambda_{r}q^{\beta})}{q+\lambda q^{\alpha +1}} \bigg(\frac{\partial^{2}}{\partial r^{2}}+\frac{1}{r}\frac{\partial}{\partial r}-\frac{4}{r^{2}}\bigg)\overline{\tau}(r,q). \end{eqnarray}
(14)
Finite Hankel transform of the function \(\overline{\tau}(r,q)\) defined as,
\begin{eqnarray} \overline{\tau}_{H}(r_{n},q)=\int\limits_0^{R} r J_{2}(rr_{n})\overline{\tau}(r,q)dr, \end{eqnarray}
(15)
identity which we used here is,
\begin{multline} \int\limits_0^{R} r J_{2}(rr_{n})\bigg(\frac{\partial^{2}}{\partial r^{2}}+\frac{1}{r}\frac{\partial}{\partial r}-\frac{4}{r^{2}} J_{2}(rr_{n})\bigg)\overline{\tau}(r,q)dr = -R\overline{\tau}(R,q)J’_{2}(Rr_{n})r_{n} -r^{2}_{n}\overline{\tau}_{H}(r_{n},q). \end{multline}
(16)
Multiplying equation (14) by \(rJ_{2}(rr_{n})\) then integrate from \(0\) to R with respect to r, we get $$\int\limits_0^{R} rJ_{2}(rr_{n})\overline{\tau}(r,q)dr=\frac{\nu (1+\lambda_{r}q^{\beta})}{q+\lambda q^{\alpha+1}}\int\limits_0^{R} r\bigg(\frac{\partial^{2}}{\partial r^{2}}+\frac{1}{r}\frac{\partial}{\partial r}-\frac{4}{r^{2}}\bigg)J_{2}(rr_{n})\overline{\tau}(r,q)dr.$$ using equations (13), (15) and (16) we get, $$\overline{\tau}_{H}(r_{n},q)=\frac{\nu (1+\lambda_{r} q^{\beta})}{q+\lambda q^{\alpha +1}}[-Rr_{n}J’_{2}(Rr_{n})\overline{\tau}(R,q)-r^{2}_{n}\overline{\tau}_{H}(r_{n},q)].$$ Now simplification for \(\overline{\tau}_{H}(r_{n},q)\) we get
\begin{eqnarray} \overline{\tau}_{H}(r_{n},q) = \frac{-R^{2}\Omega \omega J’_{2}(Rr_{n})}{(q^{2}+\omega^{2})}\frac{\nu r_{n}(1+\lambda_{r}q^{\beta})}{(q+\nu r^{2}_{n} +\lambda q^{\alpha+1}+\nu \lambda _{r} q^{\beta} r^{2}_{n})}. \end{eqnarray}
(17)
Separating the function in suitable form as, $$\overline{\tau}_{H}(r_{n},q) = \frac{-R^{2}\Omega \omega J’_{2}(Rr_{n})}{(q^{2}+\omega^{2}) r_{n}}+\frac{R^{2}\Omega \omega J’_{2}(Rr_{n})}{(q^{2}+\omega^{2})}\frac{q+\lambda q^{\alpha +1}}{r_{n}(q+\lambda q^{\alpha+1}+\nu r^{2}_{n}+\nu \lambda _{r} q^{\beta} r^{2}_{n})},$$
\begin{eqnarray} \overline{\tau}_{H}(r_{n},q)=\overline{\tau}_{1H}(r_{n},q)+\overline{\tau}_{2H}(r_{n},q), \end{eqnarray}
(18)
where,
\begin{eqnarray} \overline{\tau}_{1H}(r_{n},q)=\frac{-R^{2}\Omega \omega J’_{2}(Rr_{n})}{r_{n}(q^{2}+\omega^{2})}=-R^{2}\Omega J_{1}(Rr_{n})\frac{\omega}{q^{2}+\omega^{2}}, \end{eqnarray}
(19)
$$\overline{\tau}_{2H}(r_{n},q)=\frac{R^{2}\Omega \omega J’_{2}(Rr_{n})}{r_{n}(q^{2}+\omega^{2})}\frac{q+\lambda q^{\alpha +1}}{(q+\lambda q^{\alpha+1}+\nu r^{2}_{n}+\nu \lambda _{r} q^{\beta} r^{2}_{n})},$$
\begin{eqnarray} \overline{\tau}_{2H}(r_{n},q)=\frac{R^{2}\Omega \omega J_{1}(Rr_{n})}{q^{2}+\omega^{2}}\frac{q+\lambda q^{\alpha +1}}{(q+\lambda q^{\alpha+1}+\nu r^{2}_{n}+\nu \lambda_{r} q^{\beta} r^{2}_{n})}. \end{eqnarray}
(20)
Using the identity i.e; \begin{align*} &\frac{1}{q+\nu r^{2}_{n}+\lambda q^{\beta+1}+\nu r^{2}_{n}\lambda_{r}q^{\gamma}} \\&=\frac{1}{\lambda}\sum^{\infty}_{k=0} \sum^{k}_{m=0}\frac{k!}{(k-m)!m!}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \lambda_{r}^{m}\frac{q^{\gamma m-k-1}}{(q^{\beta}+\frac{1}{\lambda})^{k+1}}.\end{align*} Now equation (20) becomes,
\begin{multline} \overline{\tau}_{2H}(r_{n},q)=R^{2}\Omega J_{1}(Rr_{n})\frac{1}{\lambda}\sum^{\infty}_{k=0}\sum^{k}_{m=0}\frac{k!}{(k-m)!m!} \lambda_{r}^{m}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \frac{\omega}{q^{2}+\omega^{2}} \bigg[\frac{q^{\beta m-k}}{(q^{\alpha}+\frac{1}{\lambda})^{k+1}}+ \lambda \frac{q^{\alpha+\beta m-k}}{(q^{\alpha}+\frac{1}{\lambda})^{k+1}}\bigg]. \end{multline}
(21)
Taking the inverse Laplace transform of equation (19)
\begin{eqnarray} \tau_{1H}(r_{n},t)=-R^{2}\Omega J_{1}(Rr_{n})\sin(\omega t). \end{eqnarray}
(22)
Take the inverse Laplace transform and use convolution theorem of equation (21) we get
\begin{multline} \tau_{2H}(r_{n},t)=R^{2}\Omega J_{1}(Rr_{n})\frac{1}{\lambda}\sum^{\infty}_{k=0} \sum^{k}_{m=0}\frac{k!}{m!(k-m)!}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k}\lambda_{r}^{m}\bigg[\int\limits_0^{t}\sin\omega(t-s) G_{\alpha,\beta m-k,k+1}(-\lambda^{-1},s)ds+\lambda\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\alpha+\beta m-k,k+1}(-\lambda^{-1},s)ds\bigg]. \end{multline}
(23)
Where \(G_{a,b,c}(.,t)\) is the generalized G-function with \(£^{-1} \{\frac{q^{b}}{(q^{a}-d)^{c}}\} =G_{a,b,c}(d,t)\) , \(Re(ac-b)>0\) , \(Re(q)>0\) , \(\mid\frac{d}{q^{a}}\mid< 1\) and \(G_{a,b,c}(d,t)=\sum^{\infty}_{j=0}\frac{t^{(c+j)a-b-1}}{\Gamma ((c+j)a-b))}\frac{d^{j} \Gamma (c+j)}{\Gamma (c) \Gamma (j+1)}\) Taking inverse Laplace transform of equation (18) and using equation (22) and (23)
\begin{multline} \tau_{H}(r_{n},t)=-R^{2}\Omega J_{1}(Rr_{n})\sin(\omega t)+\frac{ R^{2}\Omega J_{1}(Rr_{n})}{\lambda}\sum^{\infty}_{k=0} \sum^{k}_{m=0}\frac{k!}{m!(k-m)!} \times \lambda_{r}^{m} \bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \bigg[\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\beta m-k,k+1}(-\lambda^{-1},s)ds +\lambda\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\alpha+\beta m-k,k+1}(-\lambda^{-1},s)ds\bigg]. \end{multline}
(24)
Apply the inverse Hankel transform to equation (24) and using the known formulae $$H(r^{2})=\int\limits_{0}^{R}r^{3}J_{2}(rr_{n})dr=\frac{-R^{3}J_1{Rr_{n}}}{r_{n}}$$ $$\tau(r,t)=2\sum^{\infty}_{n=1}\frac{J_{2}(rr_{n})}{[J’_{2}(r_{n})]^{2}}\tau_{H}(r_{n},t)$$
\begin{align} \tau(r,t)&=\frac{r^{2}r_{n}\Omega\sin(\omega t)}{R}+2R^{2}\Omega\sum^{\infty}_{n=1}\frac{J_{2}(rr_{n})J_{1}(Rr_{n})}{J_{1}^{2}(r_{n})}\frac{1}{\lambda}\nonumber \\ \times&\sum^{\infty}_{k=0} \sum^{k}_{m=0}\frac{k!}{(k-m)!m!}\lambda_{r}^{m}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k}\bigg[\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\beta m-k,k+1}(-\lambda^{-1},s)ds\nonumber \\&+\lambda\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\alpha+\beta m-k,k+1}(-\lambda^{-1},s)ds\bigg] \end{align}
(25)

4. Calculation for Velocity Field

Rewrite equation (8) and used equation (25) we get the non-integer order differential equation for velocity. $$\rho D_t^{\alpha} \omega(r,t) = \bigg(\frac{\partial}{\partial r}+\frac{2}{r}\bigg)\tau(r,t)$$ where $$\tau(r,t)=\tau_{1}(r,t)+\tau_{2}(r,t)$$ So above equation becomes;
\begin{eqnarray} \rho D_t^{\alpha} \omega(r,t) = \bigg(\frac{\partial}{\partial r}+\frac{2}{r}\bigg)\tau_{1}(r,t)+\bigg(\frac{\partial}{\partial r}+\frac{2}{r}\bigg)\tau_{2}(r,t) \end{eqnarray}
(26)
where $$\tau_{1}(r,t)=\frac{r^{2}r_{n}\Omega\sin(\omega t)}{R}$$ $$\tau_{2}(r,t)=2R^{2}\Omega\sum^{\infty}_{n=1}\frac{J_{2}(rr_{n})J_{1}(Rr_{n})}{J_{1}^{2}(r_{n})}\frac{1}{\lambda}\sum^{\infty}_{k=0}\sum^{k}_{m=0}\frac{k!}{(k-m)!m!}\lambda_{r}^{m}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k}$$ \begin{eqnarray*}\times\bigg[\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\beta m-k,k+1}(-\lambda^{-1},s)ds\\+\lambda\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\alpha+\beta m-k,k+1}(-\lambda^{-1},s)ds\bigg]. \end{eqnarray*} So,
\begin{eqnarray} \bigg(\frac{\partial}{\partial r}+\frac{2}{r}\bigg)\tau_{1}(r,t)=\frac{4rr_{n}\Omega\sin(\omega t)}{R}, \end{eqnarray}
(27)
and
\begin{multline} \bigg(\frac{\partial}{\partial r}+\frac{2}{r}\bigg)\tau_{2}(r,t)=\frac{2 R^{2}\Omega}{\lambda}\sum^{\infty}_{n=1}\frac{J_{1}(rr_{n})J_{1}(Rr_{n})r_{n}}{J_{1}^{2}(r_{n})}\sum^{\infty}_{k=0}\sum^{k}_{m=0} \frac{k!}{(k-m)!m!}\\\times\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k}\lambda_{r}^{m} \bigg[\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\beta m-k,k+1}(-\lambda^{-1},s)ds+\lambda\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\alpha+\beta m-k,k+1}(-\lambda^{-1},s)ds\bigg]. \end{multline}
(28)
Use equation (27) and (28) in equation (26) we get,
\begin{multline} \rho D_t^{\alpha} \omega(r,t) =\frac{4rr_{n}\Omega\sin(\omega t)}{R}+\frac{2 R^{2}\Omega}{\lambda}\sum^{\infty}_{n=1}\frac{J_{1}(rr_{n})J_{1}(Rr_{n})r_{n}}{J_{1}^{2}(r_{n})} \\ \times\sum^{\infty}_{k=0} \sum^{k}_{m=0}\frac{k!}{(k-m)!m!}\lambda_{r}^{m}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \bigg[\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\beta m-k,k+1}(-\lambda^{-1},s)ds+\lambda\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\alpha+\beta m-k,k+1}(-\lambda^{-1},s)ds\bigg]. \end{multline}
(29)
The Laplace transform of equation (29) is
\begin{multline} \varpi(r,q)=\frac{4rr_{n}\Omega}{\rho R}\frac{1}{q^{\alpha}}\frac{\omega}{q^{2}+\omega^{2}}+\frac{2 R^{2}\Omega}{\rho\lambda}\sum^{\infty}_{n=1}\frac{J_{1}(rr_{n})J_{1}(Rr_{n})r_{n}}{J_{1}^{2}(r_{n})} \sum^{\infty}_{k=0}\sum^{k}_{m=0}\frac{k!}{(k-m)!m!}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \lambda_{r}^{m}\bigg[\frac{\omega}{q^{2}+\omega^{2}}\frac{q^{\beta m-\alpha-k}}{(q^{\alpha}+\frac{1}{\lambda})^{k+1}} +\lambda\frac{\omega}{q^{2}+\omega^{2}}\frac{q^{\beta m-k}}{(q^{\alpha}+\frac{1}{\lambda})^{k+1}}\bigg]. \end{multline}
(30)
Now Apply inverse Laplace transform to equation (30) and using the Convolution theorem
\begin{multline} \omega(r,t)=\frac{4rr_{n}\Omega}{\rho R}\int\limits_{0}^{t}\sin\omega(t-\tau)\frac{t^{\alpha-1}}{\Gamma(\alpha)}d\tau+\frac{2R^{2}\Omega}{\rho \lambda}\sum^{\infty}_{n=1}\frac{J_{1}(rr_{n})J_{1}(Rr_{n})r_{n}}{J_{1}^{2}(r_{n})} \\ \times\sum^{\infty}_{k=0}\sum^{k}_{m=0}\frac{k!}{(k-m)!m!}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \lambda_{r}^{m}\bigg[\int\limits_{0}^{t}\sin\omega(t-s)G_{\alpha,\beta m-\alpha-k,k+1}(-\lambda^{-1},s)ds + \lambda\int\limits_{0}^{t}\sin\omega(t-s)G_{\alpha,\beta m-k,k+1}(-\lambda^{-1},s)ds \bigg]. \end{multline}
(31)

5. Limiting case

5.1. Ordinary Oldroyd-B fluid

Letting \(\alpha\rightarrow1\) , \(\beta\rightarrow1\) into equations (25) and (31)we get the result of shear stress and velocity field respectively for ordinary Oldroyd-B fluid.
\begin{multline} \tau_{OB}(r,t)=\frac{r^{2}r_{n}\Omega\sin(\omega t)}{R}+2R^{2}\Omega\sum^{\infty}_{n=1}\frac{J_{2}(rr_{n})J_{1}(Rr_{n})}{J_{1}^{2}(r_{n})} \frac{1}{\lambda}\sum^{\infty}_{k=0}\sum^{k}_{m=0}\frac{k!}{(k-m)!m!} \times \bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k}\lambda_{r}^{m} \bigg[\int\limits_0^{t}\sin\omega(t-s)G_{1,m-k,k+1}(-\lambda^{-1},s)ds+\lambda\int\limits_0^{t}\sin\omega(t-s)G_{1,1+ m-k,k+1}(-\lambda^{-1},s)ds\bigg], \end{multline}
(32)
\begin{multline} \omega_{OB}(r,t)=\frac{4rr_{n}\Omega}{\rho R}\int\limits_{0}^{t}\sin\omega(t-\tau)d\tau+\frac{2R^{2}\Omega}{\rho \lambda}\sum^{\infty}_{n=1}\frac{J_{1}(rr_{n})J_{1}(Rr_{n})r_{n}}{J_{1}^{2}(r_{n})} \times \sum^{\infty}_{k=0}\sum^{k}_{m=0}\frac{k!}{(k-m)!m!}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \lambda_{r}^{m}\bigg[\int\limits_{0}^{t}\sin\omega(t-s)G_{1, m-k-1,k+1}(-\lambda^{-1},s)ds\\+\lambda\int\limits_{0}^{t}\sin\omega(t-s)G_{1,m-k,k+1 }(-\lambda^{-1},s)ds \bigg]. \end{multline}
(33)

5.2. Generalized Maxwell Fluid

By placing \(\lambda_{r}\rightarrow0\) , \(\beta\rightarrow0\) into equations (25) and (31)we get the results of shear stress and velocity field for generalized Maxwell fluid respectively.
\begin{multline} \tau_{GM}(r,t)=\frac{r^{2}r_{n}\Omega\sin(\omega t)}{R}+2R^{2}\Omega\sum^{\infty}_{n=1}\frac{J_{2}(rr_{n})J_{1}(Rr_{n})}{J_{1}^{2}(r_{n})}\frac{1}{\lambda}\sum^{\infty}_{k=0}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \bigg[\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,-k,k+1}(-\lambda^{-1},s)ds +\lambda\int\limits_0^{t}\sin\omega(t-s)G_{\alpha,\alpha -k,k+1}(-\lambda^{-1},s)ds\bigg], \end{multline}
(34)
\begin{multline} \omega_{GM}(r,t)=\frac{4rr_{n}\Omega}{\rho R}\int\limits_{0}^{t}\sin\omega(t-\tau)\frac{t^{\alpha-1}}{\Gamma(\alpha)}d\tau+\frac{2R^{2}\Omega}{\rho \lambda}\sum^{\infty}_{n=1}\frac{J_{1}(rr_{n})J_{1}(Rr_{n})r_{n}}{J_{1}^{2}(r_{n})} \times \sum^{\infty}_{k=0}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \bigg[\int\limits_{0}^{t}\sin\omega(t-s)G_{\alpha,-\alpha-k,k+1}(-\lambda^{-1},s)ds\\+\lambda\int\limits_{0}^{t}\sin\omega(t-s)G_{\alpha, -k,k+1}(-\lambda^{-1},s)ds \bigg]. \end{multline}
(35)

5.3. Ordinary Maxwell Fluid

By putting \(\alpha\rightarrow1\) in equation (34) and (35) we get expressing for the shear stress and velocity field for ordinary Maxwell fluid respectively.
\begin{multline} \tau_{OM}(r,t)=\frac{r^{2}r_{n}\Omega\sin(\omega t)}{R}+2R^{2}\Omega\sum^{\infty}_{n=1}\frac{J_{2}(rr_{n})J_{1}(Rr_{n})}{J_{1}^{2}(r_{n})}\frac{1}{\lambda}\sum^{\infty}_{k=0}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \bigg[\int\limits_0^{t}\sin\omega(t-s)G_{1,-k,k+1}(-\lambda^{-1},s)ds+\lambda\int\limits_0^{t}\sin\omega(t-s)G_{1,1 -k,k+1}(-\lambda^{-1},s)ds\bigg], \end{multline}
(36)
\begin{multline} \omega_{OM}(r,t)=\frac{4rr_{n}\Omega}{\rho R}\int\limits_{0}^{t}\sin\omega(t-\tau)d\tau+\frac{2R^{2}\Omega}{\rho \lambda}\sum^{\infty}_{n=1}\frac{J_{1}(rr_{n})J_{1}(Rr_{n})r_{n}}{J_{1}^{2}(r_{n})}\sum^{\infty}_{k=0}\bigg(\frac{-\nu r^{2}_{n}}{\lambda}\bigg)^{k} \bigg[\int\limits_{0}^{t}\sin\omega(t-s)G_{1,-1-k,k+1}(-\lambda^{-1},s)ds+\lambda\int\limits_{0}^{t}\sin\omega(t-s)G_{1, -k,k+1}(-\lambda^{-1},s)ds\bigg]. \end{multline}
(37)

5.4. Newtonian Fluid

Let \(\lambda_{r}\rightarrow0\) and using the limit, i.e; $${\lim_{\lambda\rightarrow 0}}\frac{1}{\lambda^{k}}G_{1,b,k}\bigg(\frac{-1}{\lambda},t\bigg)=\frac{t^{-b-1}}{\Gamma(-b)},$$ \(b< 0\) in (25) and (31) we can get results for Newtonian fluid as,
\begin{multline} \tau_{NF}(r,t)=\frac{r^{2}r_{n}\Omega\sin(\omega t)}{R}+\frac{2R^{2}\Omega}{\omega^{2}+\nu^{2}r^{4}_{n}} \sum^{\infty}_{n=1}\frac{J_{2}(rr_{n})J_{1}(Rr_{n})}{J_{1}^{2}(r_{n})}\{\omega^{2}\sin(\omega t)+ \nu r^{2}_{n}\omega\cos(\omega t)-\omega \nu r^{2}_{n} e^{-\nu r^{2}_{n} t}\}, \end{multline}
(38)
\begin{multline} \omega_{NF}(r,t)=\frac{4rr_{n}\Omega}{\rho R}\int\limits_{0}^{t}\sin\omega(t-\tau)\frac{t^{\alpha-1}}{\Gamma(\alpha)}d\tau+\frac{2R^{2}\Omega}{\rho (\omega^{2}+\nu^{2}r^{4}_{n})}\sum^{\infty}_{n=1}\frac{J_{1}(rr_{n})J_{1}(Rr_{n})r_{n}}{J_{1}^{2}(r_{n})} \bigg[-\omega \nu r^{2}_{n}\int\limits_{0}^{t} e ^{-\nu r^{2}_{n}(t-s)}\frac{\tau^{\alpha-1}}{\Gamma(\alpha)}ds +\omega^{2}\int\limits_{0}^{t}\sin\omega(t-s)\frac{\tau^{\alpha-1}}{\Gamma(\alpha)}ds +\nu r^{2}_{n} \omega \int\limits_{0}^{t}\cos\omega(t-s)\frac{\tau^{\alpha-1}}{\Gamma(\alpha)}ds \bigg]. \end{multline}
(39)

6. Conclusion

The idea presented in this paper is to find a formula which is useful to find out exact solutions for shear stress and velocity field of any Oldroyd-B fluid which is present in rotationally oscillating cylinders. We used two transformation i.e; Hankel transform and Laplace transform. At time \(t=0^{+}\), cylinder starts its rotation about its axis. To obtain the solutions we used the finite Hankel and Laplace transforms. We express our results in the form of generalized G-function, which satisfy the governing equations and fulfilled all imposed initial and boundary conditions. Furthermore, the results for Ordinary Oldroyd-B fluid, fractional Maxwell fluid, ordinary Maxwell fluid and classical Newtonian fluid are obtained as limiting cases.

Competing Interests

The author(s) do not have any competing interests in the manuscript.

References:

  1. Batchelor, C. K., & Batchelor, G. K. (1967). An introduction to fluid dynamics. Cambridge university press. [Google Scholor]
  2. Ting, T. W. (1963). Certain non-steady flows of second-order fluids. Archive for Rational Mechanics and Analysis, 14(1), 1-26. [Google Scholor]
  3. Rajagopal, K. R., & Bhatnagar, R. K. (1995). Exact solutions for some simple flows of an Oldroyd-B fluid. Acta Mechanica, 113(1-4), 233-239. [Google Scholor]
  4. Hayat, T., Siddiqui, A. M., & Asghar, S. (2001). Some simple flows of an Oldroyd-B fluid. International Journal of Engineering Science, 39(2), 135-147. [Google Scholor]
  5. Hayat, T., Khan, M., & Ayub, M. (2004). Exact solutions of flow problems of an Oldroyd-B fluid. Applied Mathematics and Computation, 151(1), 105-119.[Google Scholor]
  6. Jamil, M., Fetecau, C., & Imran, M. (2011). Unsteady helical flows of Oldroyd-B fluids. Communications in Nonlinear Science and Numerical Simulation, 16(3), 1378-1386. [Google Scholor]
  7. Wood, W. P. (2001). Transient viscoelastic helical flows in pipes of circular and annular cross-section. Journal of non-newtonian fluid mechanics, 100(1-3), 115-126. [Google Scholor]
  8. Srivastava, P. N. (1966). Non-steady helical flow of a visco-elastic liquid(Nonsteady helical flow of viscoelastic liquid contained in circular cylinder, noting occurrence of oscillations in fluid decaying exponentially with time). Archiwum Mechaniki Stosowanej, 18(2), 145-150. [Google Scholor]
  9. Waters, N. D., & King, M. J. (1971). The unsteady flow of an elastico-viscous liquid in a straight pipe of circular cross section. Journal of Physics D: Applied Physics, 4(2), 204-211. [Google Scholor]
  10. Fetecau, C., & Fetecau, C. (1985). On the uniqueness of some helical flows of a second grade fluid. Acta mechanica, 57(3-4), 247-252. [Google Scholor]
  11. Hayat, T., Nadeem, S., & Asghar, S. (2004). Periodic unidirectional flows of a viscoelastic fluid with the fractional Maxwell model. Applied Mathematics and Computation, 151(1), 153-161.[Google Scholor]
  12. Khan, M., Nadeem, S., Hayat, T., & Siddiqui, A. M. (2005). Unsteady motions of a generalized second-grade fluid. Mathematical and Computer Modelling, 41(6-7), 629-637. [Google Scholor]
  13. Shaowei, W., & Mingyu, X. (2006). Exact solution on unsteady Couette flow of generalized Maxwell fluid with fractional derivative. Acta mechanica, 187(1-4), 103-112. [Google Scholor]
  14. Mahmood, A., Parveen, S., Ara, A., & Khan, N. A. (2009). Exact analytic solutions for the unsteady flow of a non-Newtonian fluid between two cylinders with fractional derivative model. Communications in Nonlinear Science and Numerical Simulation, 14(8), 3309-3319. [Google Scholor]
  15. Vieru, D., Akhtar, W., Fetecau, C., & Fetecau, C. (2007). Starting solutions for the oscillating motion of a Maxwell fluid in cylindrical domains. Meccanica, 42(6), 573-583. [Google Scholor]
  16. Renardy, M. (1990). An alternative approach to inflow boundary conditions for Maxwell fluids in three space dimensions. Journal of Non-Newtonian Fluid Mechanics, 36, 419-425. [Google Scholor]
  17. Waters, N. D., & King, M. J. (1970). Unsteady flow of an elastico-viscous liquid. Rheologica Acta, 9(3), 345-355. [Google Scholor]