Analysis of thin film flow of generalized Maxwell fluid confronting withdrawal and drainage on non-isothermal cylindrical surfaces

Abstract

This investigation is concerned with the study of thin film flow of a generalized Maxwell fluid along with slip conditions, confronting withdrawal and drainage on non-isothermal cylindrical surfaces. The governing equations have been formulated from the continuity equation, momentum equation, and energy equation. Analytical solutions for the velocity field, volume flow rate, average film velocity, tangential stress, and temperature are obtained in series form through the Binomial expansion technique in both withdrawal and drainage cases. The well-known solutions for a Newtonian fluid are regained as a particular case of our acquired general solutions in all flow cases. In addition, solutions for the power-law fluid model, executing alike motion, can be recovered as a limiting case of our acquired general solutions. The influence of different dimensionless parameters on all physical quantities (i.e. velocity, volume flow rate, average film velocity, tangential stress, and temperature profile) is examined and discussed graphically for both generalized Maxwell and Newtonian fluids.

Keywords

Thin film flow generalized Maxwell fluid binomial series method withdrawal and drainage non-isothermal cylindrical surfaces

Introduction

In recent time, the non-Newtonian fluids have gained astonishing interest by researchers and engineers in various branches of science and technology. Such interest is initiated by practical applications of these fluids in bioengineering, chemical industries, and material processing such as slurry fuels, colloidal and suspension solutions, exotic lubricants, extrusion of polymer fluids, and many more. Fluids belonging to this group including drilling mud, cement, shampoo, grease, ketchup, sludge, granular suspension, paints, aqueous foams, food products, plastics, paper pulp, and several others illustrate the characteristics of non-Newtonian fluids which are pretty different when equated with the linearly viscous fluids. These fluids are elucidated by a non-linear association between stress and the rate of deformation tensors, and hence, there is not a solo constitutive model which can forecast all the prominent features of non-Newtonian fluids because of their assorted physical structures. Generally, non-Newtonian fluid models are classified as (a) the rate type, (b) the differential type, and (c) the integral type. But the most prominent among them is the first one which has gained stunning importance in the field of research, especially in applied mathematics, engineering, and industry-related research problems. An imperative class of the rate type model of non-Newtonian fluids, that have the properties of both viscosity and elasticity, is viscoelastic fluids.^1,2 The study of viscoelastic fluids has many applications to industrial processes, including the cooling of metallic plates in a bath and food stuff.³ The first rate type model of viscoelastic fluids which is still extensively used to approximate the response of some dilute polymeric solutions is Maxwell fluids.⁴ This is a simplest subclass of rate type viscoelastic fluids which can designate the physical characteristics of relaxation time. Due to its important applications in industrial technology, it has gained considerable interest in the technical fields.^5,6 These applications comprise power generators, cooling of nuclear reactors, liquid metal flow control, micro pumps, high-temperature plasmas, biological transportation, drying processes, solidification of binary alloy, and so on.^3,7,8 Rajagopal and Srinivasa⁹ established a methodical thermodynamic structure for several viscoelastic fluids of rate type which is based on the influential work of Maxwell for different flow problems.

Thin liquid films are quite omnipresent in nature and technology, so an understanding of their mechanics, the amount of fluid withdrawn, and knowledge of their thickness is of great importance in numerous practical applications, such as coating of photographic films, enameling, hot tinning, lubrication of moving machine parts as in lubricating gears and bearings, chemical and nuclear reactor design, manufacturing of CDs and plastic sheets, painting, rinsing, dip coating, spreading of sauce on food, and pickling. A typical thin film flow comprises a span of fluid partly bounded by a solid surface, whereas the other surface is freely interrelated to another fluid, generally a gas and, most frequently, air in applications. For the formulation of these films, three basic conditions such as gravitational force, surface tension, and centrifugal force are required. A simple and clear example is the flow of a thin raindrop down a windowpane under the action of gravity. With the research advancement, thin film flow of non-Newtonian fluids has received significant approach by many intellectual workers and researchers^10–14 due to their fruitful applications in the vicinities of physical and biological sciences, particularly with the progress of polymer, petroleum, and other kinds of pulp industries.^15–19 Landau and Lifshitz²⁰ studied the drainage thin film flow phenomena of Newtonian fluids. Over the past few decades, significant attempts have been made by scientists^21–24 for the establishment of numerical solutions and analytical algorithms of thin film flows.^25–28 Hayat and colleagues^29–31 analyzed the thin film flow by considering differential type fluids in varieties of articles. They studied the influence of numerous physical parameters on flow fields. Siddiqui et al.^32–34 investigated the flow of thin films of third- and fourth-grade fluids down an inclined plane and vertical cylinder by means of the homotopy perturbation method. They also studied the flow behavior of non-Newtonian fluid films over a vertical moving belt. At present, the study of non-Newtonian fluid flow along with slip boundary conditions has become spirited due to the extensive usage of such fluids in polymer melt, power engineering, petroleum production, and food engineering.^35–37 Farooq et al.^38,39 explored the withdrawal and drainage on a vertical cylinder for generalized second-grade fluid with slip conditions. Gul and colleagues^40–42 probed the lifting and drainage for magnetohydrodynamic (MHD) thin film flow of third-grade fluid using constant and variable viscosities.

The purpose of this article is to venture further into the establishment of a generalized Maxwell fluid. More precisely, this article claims to examine the impact of slip condition on thin film flow of a generalized Maxwell fluid over a vertical upward moving cylinder and down a stationary vertical cylinder. Series solutions in both cases are acquired by using efficient methodical techniques, namely, binomial series method. Analytical solutions for velocity field, volume flow rate, average film velocity, tangential stress, and temperature are obtained in both withdrawal and drainage cases. The equivalent solutions for Newtonian fluid are also achieved as a particular case of our general solutions by taking the flow behavior index is equal to zero. Moreover, the influence of various physical parameters on velocity field, volume flow rate, average film velocity, tangential stress, and temperature is discussed and presented graphically in both cases under the state of generalized Maxwell and Newtonian fluids.

Governing equations

The main equations governing the motion of an incompressible generalized Maxwell fluid on non-isothermal cylindrical surfaces are^35,38

\nabla . v = 0

(1)

ρ D_{t} v + \nabla p = ρ b + \nabla . S

(2)

ρ ℘ D_{t} T = ε \nabla^{2} T + S . L

(3)

where $D_{t} = \partial / \partial t + v . \nabla$ is the material time derivative, $ρ$ is the fluid density, $S$ is the extra stress tensor, $b$ is the body force, $v$ is the velocity vector, $T$ is the temperature, $ε$ is the thermal conductivity, and $℘$ is the specific heat constant.^37,39

The Cauchy stress tensor $τ$ for a generalized Maxwell fluid is defined as³⁷

τ = - p I + S; S + λ [\overset{\cdot}{S} - LS - S L^{T}] = μ_{eff} A_{1}

(4)

where $- p I$ symbolizes the indeterminate spherical stress due to the constraint of incompressibility, $λ$ is the relaxation time, $A_{1} = L + L^{T}$ is the first Rivlin–Ericksen tensor,⁴³ $L = \nabla v$ is the velocity gradient, the superposed dot indicates material time derivative, the superscript $T$ designates the transpose operation, and $μ_{eff}$ is the effective viscosity which for a generalized Maxwell fluid is defined as^23,37

μ_{eff} = β {| \sqrt{\frac{1}{2} tr A_{1}^{2}} |}^{q}

(5)

where $β$ is the flow consistency index, $q$ is the flow behavior index with $- 1 < q < 1$ . The model characterized by constitutive equation (4) contains the significant consequences such that for $q < 0$ , the fluid is shear thinning (pseudoplastic); for $q > 0$ , the fluid is shear thickening (dilatant); and for $q = 0$ , the classical Maxwell fluid is retrieved. Furthermore, for $λ \to 0$ , the generalized Maxwell fluid model will be reduced to the power-law fluid model.^17,22

Problem formulation for withdrawal case

Consider a container filled with an incompressible generalized Maxwell fluid. An infinite vertical cylinder of radius $R$ is lifted upward through this container with constant speed $u_{0}$ . As the cylinder moves, a thin layer of the fluid of constant thickness $X$ adheres to the cylindrical surface. This film has a tendency to drain down the cylinder under the influence of gravity. We use cylindrical coordinates system such that z-axis is alongside the axis of the cylinder in the upward direction and r-axis is normal to it as shown in Figure 1.^22,37,39 Suppose that the flow is laminar, uniform, and steady, the fluid wets the solid surface completely. The force of gravity is acting in vertically downward direction. As the flow is only in z-direction, $r$ and $θ$ components of velocity are zero. We find velocity and temperature profiles of the following form^23,35,39

v = [0, 0, u (r)] and T = T (r)

(6)

Figure 1.

Geometry of withdrawal case.

By making use of equation (6), equation (1) is satisfied identically and equation (2) is written as^35,37

r - component : \frac{2 β λ}{r} D_{r} [r {(D_{r} u)}^{2} {(D_{r} u)}^{q}] - \partial_{r} p = 0

(7)

θ - component : - \frac{1}{r} \partial_{θ} p = 0 \Rightarrow \partial_{θ} p = 0

(8)

z - component : \frac{β}{r} D_{r} [r (D_{r} u) (D_{r} u)^{q}] - ρ g - \partial_{z} p = 0

(9)

Suppose that the atmospheric pressure $p =$ constant $\Rightarrow \partial_{z} p = 0$ , then equation (9) becomes

\frac{β}{r} D_{r} [r (D_{r} u) {(D_{r} u)}^{q}] - ρ g = 0 \Rightarrow D_{r} [r {(D_{r} u)}^{q + 1}] = \frac{r ρ g}{β}

(10)

Using equations (4) and (6), equation (3) becomes^35,37

\begin{matrix} ε D_{r}^{2} T + \frac{ε}{r} D_{r} T + β [{(D_{r} u)}^{2} {(D_{r} u)}^{q}] = 0 \\ \Rightarrow ε [D_{r}^{2} T + \frac{1}{r} D_{r} T] = - β {(D_{r} u)}^{q + 2} \end{matrix}

(11)

Boundary conditions related with equations (10) and (11) are^35,39

\begin{matrix} Free space boundary conditions : \\ S_{rz} = 0, D_{r} T = 0 at r = D \end{matrix}

(12)

\begin{matrix} Slip boundary conditions : \\ u = u_{0} - α S_{rz}, T = T_{0} at r = R \end{matrix}

(13)

\begin{matrix} where α is the slip coefficient, \\ D = R + X and S_{rz} = β (D_{r} u)^{q + 1} \end{matrix}

(14)

By integration of equation (10) with respect to “ $r$ ” and consuming the boundary conditions given in equation (12), we have ^37,39

D_{r} u = (- 1)^{\frac{1}{q + 1}} {\frac{ρ g}{2 β}}^{\frac{1}{q + 1}} {[\frac{D^{2} - r^{2}}{r}]}^{\frac{1}{q + 1}}

(15)

In series form, using binomial expansion, equation (15) becomes

D_{r} u = (- 1)^{\frac{1}{q + 1}} {\frac{ρ g}{2 β}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) (- 1)^{k} r^{2 k - \frac{1}{q + 1}} D^{- 2 k + \frac{2}{q + 1}}

(16)

Equation (11) in view of equation (16) becomes

\begin{matrix} ε [D_{r}^{2} T + \frac{1}{r} D_{r} T] = - β (- 1)^{\frac{q + 2}{q + 1}} {\frac{ρ g}{2 β}}^{\frac{q + 2}{q + 1}} \\ \sum_{k = 0}^{\infty} (\begin{matrix} \frac{q + 2}{q + 1} \\ k \end{matrix}) (- 1)^{k} r^{2 k - \frac{q + 2}{q + 1}} D^{- 2 k + \frac{2 q + 4}{q + 1}} \end{matrix}

(17)

Solutions for withdrawal case

Solutions for the generalized Maxwell fluid and for Newtonian fluid are presented here in this section.

Solutions for generalized Maxwell fluid

For $q \neq 0$ , the following expressions are obtained for the generalized Maxwell fluid.

Velocity field

Integrating equation (16) with respect to “ $r$ ” and inserting the given boundary conditions, the expression for velocity field is obtained as

\begin{matrix} u = u_{0} + \frac{α ρ g}{2} [\frac{D^{2} - R^{2}}{R}] + (- 1)^{\frac{1}{q + 1}} {\frac{ρ g}{2 β}}^{\frac{1}{q + 1}} \\ \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \frac{{(- 1)}^{k} [r^{2 k + \frac{q}{q + 1}} - R^{2 k + \frac{q}{q + 1}}]}{(2 k + \frac{q}{q + 1}) D^{2 k - \frac{2}{q + 1}}} \end{matrix}

(18)

Volume flow rate

Volume flow rate $Q$ of the fluid passing through the surface of the cylinder is defined as^37,39

Q = \int \int ru (r) drd θ, R \leq r \leq D, 0 \leq θ \leq 2 π

(19)

Using equation (18), we acquire

\begin{array}{l} Q = π (D^{2} - R^{2}) [u_{0} + \frac{α ρ g}{2 R} (D^{2} - R^{2})] \\ + 2 π {(- 1)}^{\frac{1}{q + 1}} {\frac{ρ g}{2 β}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) {(- 1)}^{k} \\ \times [\frac{D^{2 k + \frac{3 q + 2}{q + 1}} - R^{2 k + \frac{3 q + 2}{q + 1}}}{(2 k + \frac{q}{q + 1}) D^{2 k - \frac{2}{q + 1}} (2 k + \frac{3 q + 2}{q + 1})} - \frac{R^{2 k + \frac{q}{q + 1}} (D^{2} - R^{2})}{2 (2 k + \frac{q}{q + 1}) D^{2 k - \frac{2}{q + 1}}}] \end{array}

(20)

Average film velocity

Average film velocity $U$ is defined as^37,39

U = \frac{Q}{π (D^{2} - R^{2})}

(21)

Inserting equation (20) into equation (21), the average film velocity is

\begin{matrix} U = u_{0} + \frac{α ρ g}{2 R} [D^{2} - R^{2}] + 2 {- \frac{ρ g}{2 β}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \\ \frac{{(- 1)}^{k} D^{- 2 k + \frac{2}{q + 1}}}{2 k + \frac{q}{q + 1}} [\frac{D^{2 k + \frac{3 q + 2}{q + 1}} - R^{2 k + \frac{3 q + 2}{q + 1}}}{(2 k + \frac{3 q + 2}{q + 1}) (D^{2} - R^{2})} - \frac{R^{2 k + \frac{q}{q + 1}}}{2}] \end{matrix}

(22)

Tangential stress

The tangential stress on cylindrical surface is given in equation (14). Using equation (4), we have

S_{rz} = - \frac{ρ g}{2} [\frac{D^{2}}{r} - r]

(23)

Consuming the slip boundary conditions given in equation (13), the tangential stress employed by the fluid on cylindrical surface in withdrawal case is

S_{rz} |_{r = R} = - \frac{ρ g}{2} [\frac{D^{2} - R^{2}}{R}]

(24)

Temperature profile

Solution of equation (17) for temperature profile by using boundary conditions defined in equations (12) and (13) is

\begin{matrix} T = T_{0} + \frac{β}{ε} (- 1)^{\frac{q + 2}{q + 1}} {\frac{ρ g}{2 β}}^{\frac{q + 2}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{q + 2}{q + 1} \\ k \end{matrix}) \\ \frac{{(- 1)}^{k} D^{- 2 k + \frac{2 q + 4}{q + 1}}}{2 k + \frac{q}{q + 1}} [\frac{R^{2 k + \frac{q}{q + 1}} - r^{2 k + \frac{q}{q + 1}}}{2 k + \frac{q}{q + 1}} + D^{2 k + \frac{q}{q + 1}} \ln {\frac{r}{R}}] \end{matrix}

(25)

Solutions for Newtonian fluid

When $q = λ = 0$ , the generalized Maxwell fluid diminishes to a Newtonian fluid which is the particular case of our obtained solutions for withdrawal case.

Velocity field

Putting $q = 0$ in equation (15) and consuming the given boundary condition, the expression for velocity field is

u = u_{0} - \frac{ρ g}{4 β} [2 D^{2} \ln {\frac{r}{R}} + (R^{2} - r^{2}) - 2 α β {\frac{D^{2} - R^{2}}{R}}]

(26)

Volume flow rate

Volume flow rate $Q$ for Newtonian fluid by using equation (26) in equation (19) is

\begin{matrix} Q = - \frac{π ρ g}{2 β} [D^{4} \ln {\frac{D}{R}}] + π (D^{2} - R^{2}) \\ [u_{0} + \frac{ρ g}{8 β} {2 D^{2} + [1 + \frac{4 α β}{R}] (D^{2} - R^{2})}] \end{matrix}

(27)

Average film velocity

By using equation (27) into equation (21), the average film velocity $U$ is

\begin{matrix} U = u_{0} - \frac{ρ g}{8 β} \\ {\frac{4 D^{4}}{D^{2} - R^{2}} \ln {\frac{D}{R}} - 2 D^{2} - [1 + \frac{4 α β}{R}] (D^{2} - R^{2})} \end{matrix}

(28)

Tangential stress

Tangential stress for Newtonian fluid is the same as it is for generalized Maxwell fluid, which is given by equation (24).

Temperature profile

Putting $q = 0$ and consuming boundary conditions given in equations (12) and (13), the solution of equation (11) for temperature profile is

T = T_{0} - \frac{β}{ε} {(\frac{ρ g}{8 β})}^{2} [4 D^{4} (4 \ln D - 3) \ln {\frac{R}{r}} + 8 D^{4} {{[\ln r]}^{2} - {[\ln R]}^{2}} + (r^{4} - R^{4}) + 8 D^{2} (R^{2} - r^{2})]

(29)

Solutions for withdrawal case in dimensionless form

Introducing the following dimensionless parameters^22,35,39

\begin{matrix} \hat{r} = \frac{r}{R}, \hat{u} = \frac{u}{u_{0}}, \hat{H} = \frac{D}{R}, \hat{T} = \frac{T - T_{0}}{T_{1} - T_{0}}, {\hat{μ}}_{eff} = β {(\frac{u_{0}}{R})}^{q}, \\ S_{t} = \frac{ρ g R^{2}}{u_{0} {\hat{μ}}_{eff}}, B_{r} = \frac{u_{0}^{2} {\hat{μ}}_{eff}}{ε (T_{1} - T_{0})}, {\hat{S}}_{rz} = \frac{1}{β} {[\frac{R}{u_{0}}]}^{q + 1} S_{rz} \end{matrix}

(30)

in differential equations (10) and (11) and in the boundary conditions given by equations (12) and (13). The new transformed differential equations along with boundary conditions in dimensionless form after dropping the “hats” for simplicity are³⁵

\frac{1}{r} D_{r} [r (D_{r} u) {(D_{r} u)}^{q}] - S_{t} = 0 \Rightarrow D_{r} [r {(D_{r} u)}^{q + 1}] = r S_{t}

(31)

\begin{matrix} D_{r}^{2} T + \frac{1}{r} D_{r} T + B_{r} [{(D_{r} u)}^{2} {(D_{r} u)}^{q}] = 0 \\ \Rightarrow D_{r}^{2} T + \frac{1}{r} D_{r} T = - B_{r} {(D_{r} u)}^{q + 2} \end{matrix}

(32)

\begin{matrix} Free space boundary condition : \\ D_{r} u = 0 and D_{r} T = 0 at r = H \end{matrix}

(33)

\begin{matrix} Slip boundary condition : \\ u = 1 - \frac{α}{u_{0}} A {(D_{r} u)}^{q + 1} and T = 0 at r = 1 \end{matrix}

(34)

where $S_{t}$ is the Stokes number, $B_{r}$ is the Brinkman number, $T_{1}$ is the reference temperature, $A = β (u_{0} / R)^{q + 1}$ , and $H$ is the fluid parameter.

By integration of equation (31) with respect to “r” and consuming boundary conditions defined in equation (33), we have

D_{r} u = (- 1)^{\frac{1}{q + 1}} {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} {[\frac{H^{2} - r^{2}}{r}]}^{\frac{1}{q + 1}}

(35)

In the above expression, signs of $D_{r} u$ and $r$ are always opposite to each other, since velocity of the fluid decreases as $r$ increases. In series form, using binomial expansion, equation (35) is written as³⁵

D_{r} u = (- 1)^{\frac{1}{q + 1}} {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) (- 1)^{k} r^{2 k - \frac{1}{q + 1}} H^{- 2 k + \frac{2}{q + 1}}

(36)

Equation (32) in view of equation (36) yields

\begin{matrix} D_{r}^{2} T + \frac{1}{r} D_{r} T = - B_{r} (- 1)^{\frac{q + 2}{q + 1}} {\frac{S_{t}}{2}}^{\frac{q + 2}{q + 1}} \\ \sum_{k = 0}^{\infty} (\begin{matrix} \frac{q + 2}{q + 1} \\ k \end{matrix}) (- 1)^{k} r^{2 k - \frac{q + 2}{q + 1}} H^{- 2 k + \frac{2 q + 4}{q + 1}} \end{matrix}

(37)

Inserting equation (35) into equation (34), we have

u = 1 + \frac{α S_{t}}{2 u_{0}} A [H^{2} - 1] and T = 0 at r = 1

(38)

Solutions for generalized Maxwell fluid

For $q \neq 0$ , the following expressions are obtained for generalized Maxwell fluid.

Velocity field

Integrating equation (36) with respect to “ $r$ ” and inserting boundary conditions given in equation (38), the expression for velocity field is acquired as

\begin{matrix} u = 1 + \frac{α S_{t}}{2 u_{0}} A [H^{2} - 1] + (- 1)^{\frac{1}{q + 1}} \\ {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \frac{{(- 1)}^{k} [r^{2 k + \frac{q}{q + 1}} - 1] H^{- 2 k + \frac{2}{q + 1}}}{(2 k + \frac{q}{q + 1})} \end{matrix}

(39)

Volume flow rate

Volume flow rate $Q$ of the fluid passing through cylindrical surface is defined as

Q = \int \int r u (r) drd θ, 1 \leq r \leq H, 0 \leq θ \leq 2 π

(40)

Inserting equation (39) into equation (40), the volume flow rate is obtained as³⁵

\begin{matrix} Q = π (H^{2} - 1) [1 + \frac{α S_{t}}{2 u_{0}} A (H^{2} - 1)] \\ + 2 π (- 1)^{\frac{1}{q + 1}} {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \frac{{(- 1)}^{k} H^{- 2 k + \frac{2}{q + 1}}}{(2 k + \frac{q}{q + 1})} \\ \times [\frac{H^{2 k + \frac{3 q + 2}{q + 1}} - 1}{2 k + \frac{3 q + 2}{q + 1}} - \frac{H^{2} - 1}{2}] \end{matrix}

(41)

Average film velocity

Average film velocity $U$ is given by Ullah³⁵

U = \frac{Q}{π (H^{2} - 1)}

(42)

Inserting equation (41) into equation (42), the average film velocity $U$ for withdrawal case is achieved as

\begin{matrix} U = 1 + \frac{α S_{t}}{2 u_{0}} A (H^{2} - 1) + 2 (- 1)^{\frac{1}{q + 1}} {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} \\ \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \frac{{(- 1)}^{k} H^{- 2 k + \frac{2}{q + 1}}}{(2 k + \frac{q}{q + 1}) [H^{2} - 1]} [\frac{H^{2 k + \frac{3 q + 2}{q + 1}} - 1}{2 k + \frac{3 q + 2}{q + 1}} - \frac{H^{2} - 1}{2}] \end{matrix}

(43)

Tangential stress

Tangential stress employed by the fluid over the cylindrical surface is

S_{rz} = {(D_{r} u)}^{q + 1}

(44)

Using equation (35) into equation (44), tangential stress is

S_{rz} |_{r = 1} = - \frac{S_{t} [H^{2} - 1]}{2}

(45)

Temperature profile

Solving equation (37) using boundary conditions given in equations (33) and (34), expression for temperature profile is

\begin{matrix} T = B_{r} (- 1)^{\frac{q + 2}{q + 1}} {\frac{S_{t}}{2}}^{\frac{q + 2}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{q + 2}{q + 1} \\ k \end{matrix}) \\ \frac{{(- 1)}^{k} H^{- 2 k + \frac{2 q + 4}{q + 1}}}{2 k + \frac{q}{q + 1}} [H^{2 k + \frac{q}{q + 1}} (\ln r) + \frac{1 - r^{2 k + \frac{q}{q + 1}}}{2 k + \frac{q}{q + 1}}] \end{matrix}

(46)

Solutions for Newtonian fluid

When $q = λ = 0$ , the generalized Maxwell fluid diminishes to Newtonian fluid which is the particular case of our general solutions. The exact solutions for Newtonian fluid in dimensionless form are presented below.

Velocity field

Putting $q = 0$ and consuming boundary conditions given in equation (38), solution of equation (35) for velocity field is

u = 1 + \frac{S_{t}}{2} [\frac{α β}{R} (H^{2} - 1) - H^{2} (\ln r) + \frac{r^{2} - 1}{2}]

(47)

Volume flow rate

Volume flow rate $Q$ is obtained from equation (40) by using equation (47), which is

\begin{matrix} Q = π (H^{2} - 1) [1 + \frac{S_{t}}{4} {\frac{2 α β}{R} (H^{2} - 1) - 1}] \\ - \frac{π S_{t}}{8} [4 H^{4} (\ln H) - (H^{2} - 1) (3 H^{2} + 1)] \end{matrix}

(48)

Average film velocity

By using equation (48) into equation (42), the average film velocity $U$ is

U = 1 + \frac{S_{t}}{8} [\frac{4 α β}{R} (H^{2} - 1) - \frac{4 H^{4} (\ln H)}{H^{2} - 1} + 3 H^{2} - 1]

(49)

Tangential stress

Using equation (35) into equation (44) and then solving for $q = 0$ , we obtain tangential stress for Newtonian fluid which is similar as it is calculated for generalized Maxwell fluid given by equation (45).

Temperature profile

Putting $q = 0$ and consuming boundary conditions given by equations (33) and (34), the solution of equation (32) for temperature profile is

\begin{matrix} T = \frac{- B_{r} S_{t}^{2}}{64} \\ [4 H^{4} {3 + 2 \ln r - 4 \ln H} \ln r + 8 H^{2} (1 - r^{2}) + r^{4} - 1] \end{matrix}

(50)

Problem formulation for drainage case

Consider a steady, uniform thin film of generalized Maxwell fluid draining down over a stationary vertical cylinder. In drainage case, z-axis is in the downward direction and the motion of fluid is due to the gravity as shown in Figure 2^21,22,37 The geometry and assumptions for drainage case are alike as in the preceding case, but the boundary conditions of velocity are changed. Therefore, equations (10) and (11) turn into the following form^23,35,39

\begin{matrix} \frac{β}{r} D_{r} [r (D_{r} u) {(D_{r} u)}^{q}] + ρ g = 0 \Rightarrow D_{r} [r {(D_{r} u)}^{q + 1}] \\ = - \frac{r ρ g}{β} \end{matrix}

(51)

\begin{matrix} ε D_{r}^{2} T + \frac{ε}{r} D_{r} T + β [{(D_{r} u)}^{2} {(D_{r} u)}^{q}] = 0 \\ \Rightarrow ε [D_{r}^{2} T + \frac{1}{r} D_{r} T] = - β {(D_{r} u)}^{q + 2} \end{matrix}

(52)

Figure 2.

Geometry of drainage case.

Boundary conditions connected with equations (51) and (52) are^35,37,39

\begin{matrix} Free space boundary condition : \\ S_{rz} = 0, D_{r} T = 0 at r = D \end{matrix}

(53)

\begin{matrix} Slip boundary condition : \\ u = - α S_{rz}, T = T_{0} at r = R \end{matrix}

(54)

where $D$ and $S_{rz}$ are given in equation (14). By integration of equation (51) with respect to “r” and then using equations (14) and (53), we have^37,39

{(D_{r} u)}^{q + 1} = \frac{ρ g}{2 β} [\frac{D^{2} - r^{2}}{r}] \Rightarrow D_{r} u = {\frac{ρ g}{2 β}}^{\frac{1}{q + 1}} {[\frac{D^{2} - r^{2}}{r}]}^{\frac{1}{q + 1}}

(55)

In series form, using binomial expansion, equation (55) becomes

D_{r} u = {\frac{ρ g}{2 β}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) (- 1)^{k} r^{2 k - \frac{1}{q + 1}} D^{- 2 k + \frac{2}{q + 1}}

(56)

Equation (52) in view of equation (56) yields

\begin{matrix} ε [D_{r}^{2} T + \frac{1}{r} D_{r} T] = - β {\frac{ρ g}{2 β}}^{\frac{q + 2}{q + 1}} \\ \sum_{k = 0}^{\infty} (\begin{matrix} \frac{q + 2}{q + 1} \\ k \end{matrix}) (- 1)^{k} r^{2 k - \frac{q + 2}{q + 1}} D^{- 2 k + \frac{2 q + 4}{q + 1}} \end{matrix}

(57)

Solutions for drainage case

Solutions for generalized Maxwell fluid and also for Newtonian fluid are presented here in this section.

Solutions for generalized Maxwell fluid

For $q \neq 0$ , the following expressions are obtained for generalized Maxwell fluid.

Velocity field

Solving equation (56) using the given boundary conditions, the expression for velocity field in drainage case is acquired as

\begin{matrix} u = - \frac{α ρ g}{2} [\frac{D^{2} - R^{2}}{R}] + {\frac{ρ g}{2 β}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \\ \frac{{(- 1)}^{k} [r^{2 k + \frac{q}{q + 1}} - R^{2 k + \frac{q}{q + 1}}]}{(2 k + \frac{q}{q + 1}) D^{2 k - \frac{2}{q + 1}}} \end{matrix}

(58)

Volume flow rate

Volume flow rate $Q$ in drainage case is obtained by using equation (58) into equation (19) as

\begin{matrix} Q = - \frac{π α ρ g}{2 R} {[D^{2} - R^{2}]}^{2} + 2 π {\frac{ρ g}{2 β}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \\ (- 1)^{k} [\frac{D^{2 k + \frac{3 q + 2}{q + 1}} - R^{2 k + \frac{3 q + 2}{q + 1}}}{(2 k + \frac{q}{q + 1}) D^{2 k - \frac{2}{q + 1}} (2 k + \frac{3 q + 2}{q + 1})} \\ - \frac{R^{2 k + \frac{q}{q + 1}} (D^{2} - R^{2})}{2 (2 k + \frac{q}{q + 1}) D^{2 k - \frac{2}{q + 1}}}] \end{matrix}

(59)

Average film velocity

Expression for average film velocity $U$ is attained by using equation (59) into equation (21) as under

\begin{matrix} U = - \frac{α ρ g}{2 R} [D^{2} - R^{2}] + 2 {\frac{ρ g}{2 β}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) (- 1)^{k} \\ [\frac{D^{2 k + \frac{3 q + 2}{q + 1}} - R^{2 k + \frac{3 q + 2}{q + 1}}}{(2 k + \frac{q}{q + 1}) D^{2 k - \frac{2}{q + 1}} (2 k + \frac{3 q + 2}{q + 1}) (D^{2} - R^{2})} \\ - \frac{R^{2 k + \frac{q}{q + 1}}}{2 (2 k + \frac{q}{q + 1}) D^{- 2 k + \frac{2}{q + 1}}}] \end{matrix}

(60)

Tangential stress

Tangential stress exerted on cylindrical surface in this case is

S_{rz} = \frac{ρ g}{2} [\frac{D^{2} - r^{2}}{r}] \Rightarrow S_{rz} |_{r = R} = \frac{ρ g}{2} [\frac{D^{2} - R^{2}}{R}]

(61)

Temperature profile

Solution of equation (57) for temperature profile, when the boundary conditions defined by equations (53) and (54) are used, is

\begin{matrix} T = T_{0} + \frac{β}{ε} {\frac{ρ g}{2 β}}^{\frac{q + 2}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{q + 2}{q + 1} \\ k \end{matrix}) \frac{{(- 1)}^{k} D^{- 2 k + \frac{2 q + 4}{q + 1}}}{2 k + \frac{q}{q + 1}} \\ [\frac{R^{2 k + \frac{q}{q + 1}} - r^{2 k + \frac{q}{q + 1}}}{2 k + \frac{q}{q + 1}} + D^{2 k + \frac{q}{q + 1}} \ln {\frac{r}{R}}] \end{matrix}

(62)

Solutions for Newtonian fluid

When $q = λ = 0$ , the generalized Maxwell fluid diminishes to a Newtonian fluid which is the particular case of our obtained solutions for drainage case.

Velocity field

Putting $q = 0$ and consuming boundary conditions defined in equation (53) and (54), solution of equation (51) for velocity field is

u = \frac{ρ g}{4 β} [2 D^{2} \ln {\frac{r}{R}} + (R^{2} - r^{2}) - 2 α β {\frac{D^{2} - R^{2}}{R}}]

(63)

Volume flow rate

The volume flow rate $Q$ for Newtonian fluid in drainage case is

\begin{matrix} Q = \frac{π ρ g}{2 β} [D^{4} \ln {\frac{D}{R}}] - \frac{π ρ g}{8 β} (D^{2} - R^{2}) \\ [2 D^{2} + {1 + \frac{4 α β}{R}} (D^{2} - R^{2})] \end{matrix}

(64)

Average film velocity

By utilizing equation (64) in equation (21), the average film velocity $U$ is

\begin{matrix} U = \frac{ρ g}{8 β} [\frac{4 D^{4}}{D^{2} - R^{2}} \ln {\frac{D}{R}} - 2 D^{2} - {1 + \frac{4 α β}{R}} (D^{2} - R^{2})] \end{matrix}

(65)

Tangential stress

Tangential stress for Newtonian fluid is the same as it is for generalized Maxwell fluid, which is given by equation (61).

Temperature profile

Putting $q = 0$ and consuming boundary conditions specified in equations (53) and (54), the solution of equation (52) for temperature profile is

\begin{matrix} T = T_{0} - \frac{β}{ε} {(\frac{ρ g}{8 β})}^{2} [4 D^{4} (4 \ln D - 3) \ln {\frac{R}{r}} \\ + 8 D^{4} {{(\ln r)}^{2} - {(\ln R)}^{2}} + (r^{4} - R^{4}) + 8 D^{2} (R^{2} - r^{2})] \end{matrix}

(66)

Solutions for drainage case in dimensionless form

Inserting the non-dimensional parameters defined in equation (30) into equations (51), (52), (53) and (54). After dropping “hats” for simplicity, the new converted dimensionless equations together with boundary conditions are³⁵

\frac{1}{r} D_{r} [r (D_{r} u) {(D_{r} u)}^{q}] + S_{t} = 0 \Rightarrow D_{r} [r {(D_{r} u)}^{q + 1}] = - r S_{t}

(67)

\begin{matrix} D_{r}^{2} T + \frac{1}{r} D_{r} T + B_{r} [{(D_{r} u)}^{2} {(D_{r} u)}^{q}] = 0 \\ \Rightarrow D_{r}^{2} T + \frac{1}{r} D_{r} T = - B_{r} {(D_{r} u)}^{q + 2} \end{matrix}

(68)

\begin{matrix} Free space boundary condition : \\ D_{r} u = 0 and D_{r} T = 0 at r = H \end{matrix}

(69)

\begin{matrix} Slip boundary condition : \\ u = - \frac{α}{u_{0}} A (D_{r} u)^{q + 1} and T = 0 at r = 1 \end{matrix}

(70)

By integration of equation (67) with respect to “r” and using the given boundary conditions, we acquire

D_{r} u = {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} {[\frac{H^{2} - r^{2}}{r}]}^{\frac{1}{q + 1}}

(71)

In series form, using binomial expansion, equation (71) is written as³⁵

D_{r} u = {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) (- 1)^{k} r^{2 k - \frac{1}{q + 1}} H^{- 2 k + \frac{2}{q + 1}}

(72)

Using equation (72), equation (68) yields

\begin{matrix} D_{r}^{2} T + \frac{1}{r} D_{r} T = - B_{r} {\frac{S_{t}}{2}}^{\frac{q + 2}{q + 1}} \\ \sum_{k = 0}^{\infty} (\begin{matrix} \frac{q + 2}{q + 1} \\ k \end{matrix}) (- 1)^{k} r^{2 k - \frac{q + 2}{q + 1}} H^{- 2 k + \frac{2 q + 4}{q + 1}} \end{matrix}

(73)

Equation (70) in view of equation (71) gives

u = - \frac{α S_{t}}{2 u_{0}} A [H^{2} - 1] and T = 0 at r = 1

(74)

Solutions for generalized Maxwell fluid

For $q \neq 0$ , the following expressions are obtained for generalized Maxwell fluid in drainage case.

Velocity field

Solving equation (72) using boundary conditions defined in equation (74), the expression for velocity field is

\begin{matrix} u = - \frac{α S_{t}}{2 u_{0}} A [H^{2} - 1] + {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \\ \frac{{(- 1)}^{k} [r^{2 k + \frac{q}{q + 1}} - 1] H^{- 2 k + \frac{2}{q + 1}}}{2 k + \frac{q}{q + 1}} \end{matrix}

(75)

Volume flow rate

Using equation (75) into equation (40), the volume flow rate $Q$ in drainage case is

\begin{matrix} Q = - \frac{π α S_{t}}{2 u_{0}} A {[H^{2} - 1]}^{2} + 2 π {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \\ \frac{{(- 1)}^{k} H^{- 2 k + \frac{2}{q + 1}}}{2 k + \frac{q}{q + 1}} [\frac{H^{2 k + \frac{3 q + 2}{q + 1}} - 1}{2 k + \frac{3 q + 2}{q + 1}} - \frac{H^{2} - 1}{2}] \end{matrix}

(76)

Average film velocity

Profile for average film velocity $U$ is attained by using equation (76) into equation (42) as under

\begin{matrix} U = - \frac{α S_{t}}{2 u_{0}} A [H^{2} - 1] + 2 {\frac{S_{t}}{2}}^{\frac{1}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{1}{q + 1} \\ k \end{matrix}) \\ \frac{{(- 1)}^{k} H^{- 2 k + \frac{2}{q + 1}}}{(2 k + \frac{q}{q + 1}) [H^{2} - 1]} [\frac{H^{2 k + \frac{3 q + 2}{q + 1}} - 1}{2 k + \frac{3 q + 2}{q + 1}} - \frac{H^{2} - 1}{2}] \end{matrix}

(77)

Tangential stress

Tangential stress exerted on cylindrical surface in the drainage case is

S_{rz} = \frac{S_{t}}{2} [\frac{H^{2} - r^{2}}{r}] \Rightarrow S_{rz} |_{r = 1} = \frac{S_{t}}{2} [H^{2} - 1]

(78)

Temperature profile

Solution of equation (73) for temperature profile, when the boundary conditions specified in equations (69) and (74) are used, is

\begin{matrix} T = B_{r} {\frac{S_{t}}{2}}^{\frac{q + 2}{q + 1}} \sum_{k = 0}^{\infty} (\begin{matrix} \frac{q + 2}{q + 1} \\ k \end{matrix}) \frac{{(- 1)}^{k} H^{- 2 k + \frac{2 q + 4}{q + 1}}}{2 k + \frac{q}{q + 1}} \\ [H^{2 k + \frac{q}{q + 1}} \ln (r) + \frac{1 - r^{2 k + \frac{q}{q + 1}}}{2 k + \frac{q}{q + 1}}] \end{matrix}

(79)

Solutions for Newtonian fluid

When $q = λ = 0$ , the generalized Maxwell fluid diminishes to a Newtonian fluid which is the particular case of our general solutions obtained for drainage case.

Velocity field

Putting $q = 0$ and consuming boundary conditions defined by equations (69) and (74), solution of equation (67) for velocity field is

u = \frac{S_{t}}{4} [2 H^{2} \ln (r) - r^{2} + 1 - \frac{2 α β}{R} (H^{2} - 1)]

(80)

Volume flow rate

Volume flow rate $Q$ in this case is obtained by using equation (80) in equation (40) as

Q = \frac{π S_{t}}{8} [4 H^{4} \ln H - (H^{2} - 1) {2 H^{2} + (1 + \frac{4 α β}{R}) [H^{2} - 1]}]

(81)

Average film velocity

By utilizing equation (81) into equation (42), the average film velocity $U$ is

U = \frac{S_{t}}{8} [\frac{4 H^{4} \ln H}{H^{2} - 1} - 2 H^{2} - {1 + \frac{4 α β}{R}} (H^{2} - 1)]

(82)

Tangential stress

Tangential stress for Newtonian fluid is same as it is for generalized Maxwell fluid given by equation (78).

Temperature profile

Putting $q = 0$ and consuming boundary conditions specified by equations (69) and (70), solution of the differential equation (68) for temperature profile is

\begin{matrix} T = \frac{- B_{r} S_{t}^{2}}{64} \\ [4 H^{4} {3 + 2 \ln r - 4 \ln H} \ln r + 8 H^{2} (1 - r^{2}) + r^{4} - 1] \end{matrix}

(83)

Results and discussion

In this article, we analyzed the thin film flow of generalized Maxwell fluid over a vertical cylinder under the impact of slip and free space boundary conditions. The problem has been modeled for withdrawal and drainage of velocity field and temperature profile. Furthermore, series solutions for volume flow rate, average film velocity, and tangential stress have been predicted by using the binomial series technique. The geometry of problem has been shown in Figures 1 and 2 for withdrawal and drainage, respectively. The influences of dimensionless model parameters, for instance, flow behavior index $q$ , Stocks number $S_{t}$ , Brinkman number $B_{r}$ , and fluid parameter $H$ , have been investigated on velocity, volume flow rate, average film velocity, tangential stress, and temperature for the proposed problem. The impacts of these pertinent parameters have been reflected in Figures 3 –24 which are plotted for real solutions.

Figure 3.

Impact of $S_{t}$ on velocity when $0 < q < 1$ : (a) withdrawal case and (b) drainage case.

Figure 4.

Effect of $S_{t}$ on velocity when $- 1 < q < 0$ : (a) withdrawal case and (b) drainage case.

Figure 5.

Impact of $S_{t}$ on velocity of Newtonian fluid: (a) withdrawal case and (b) drainage case.

Figure 6.

Impact of $H$ on velocity when $0 < q < 1$ : (a) withdrawal case and (b) drainage case.

Figure 7.

Effect of $H$ on velocity when $- 1 < q < 0$ : (a) withdrawal case and (b) drainage case.

Figure 8.

Impact of $H$ on velocity of Newtonian fluid: (a) withdrawal case and (b) drainage case.

Figure 9.

Impact of $q$ on velocity: (a) withdrawal case and (b) drainage case.

Figure 10.

Impact of $S_{t}$ on volume flow rate against $q$ : (a) withdrawal case and (b) drainage case.

Figure 11.

Impact of $S_{t}$ on volume flow rate against $H$ : (a) withdrawal case and (b) drainage case.

Figure 12.

Impact of $S_{t}$ on volume flow rate of Newtonian fluid: (a) withdrawal case and (b) drainage case.

Figure 13.

Impact of $q$ on volume flow rate against $H$ : (a) withdrawal case and (b) drainage case.

Figure 14.

Effect of $S_{t}$ on average film velocity against $q$ : (a) withdrawal case and (b) drainage case.

Figure 15.

Effect of $S_{t}$ on average film velocity against $H$ : (a) withdrawal case and (b) drainage case.

Figure 16.

Influence of $S_{t}$ on average film velocity of Newtonian fluid: (a) withdrawal case and (b) drainage case.

Figure 17.

Effect of $H$ on average film velocity against $q$ : (a) withdrawal case and (b) drainage case.

Figure 18.

Effect of $q$ on average film velocity against $H$ : (a) withdrawal case and (b) drainage case.

Figure 19.

Effect of $S_{t}$ on tangential stress: (a) withdrawal case and (b) drainage case.

Figure 20.

Effect of $H$ on tangential stress: (a) withdrawal case and (b) drainage case.

Figure 21.

Impacts of (a) S_t and (b) B_r, on temperature when 0 < q < 1.

Figure 22.

Impacts of (a) S_t and (b) B_r, on temperature when −1 < q < 0.

Figure 23.

Effect of $H$ on temperature: (a) $0 < q < 1$ and (b) $- 1 < q < 0$ .

Figure 24.

Impact of $q$ on temperature for withdrawal and drainage cases.

It has been investigated that, when $S_{t}$ increases, velocity decreases in withdrawal case, either the fluid is dilatant, pseudoplastic, or Newtonian, and this effect is shown graphically by Figures 3(a), 4(a), and 5(a), while Figures 3(b), 4(b), and 5(b) show that the magnitude of velocity increases in drainage case, respectively, for the above stated fluid types.

Figures 6 –8 show the effect of $H$ on velocity field for both cases. We diagnosed that velocity decreases when the value of $H$ increases in withdrawal case, either the fluid is shear thickening (dilatant), shear thinning (pseudoplastic), or Newtonian as this impact is given by Figures 6(a), 7(a), and 8(a), respectively, and in the drainage case, velocity rises with the rise in $H$ as shown by Figures 6(b), 7(b), and 8(b), respectively, for the above stated fluid types.

Figure 9(a) shows that velocity declines as the fluid is becoming thicker in withdrawal case, and it rises in the drainage case, as shown by Figure 9(b).

Figures 10 –12 display the influence of Stock’s number $S_{t}$ on the volume flow rate $Q$ for modeled problem. Figures 10(a), 11(a), and 12(a) indicate that with the increase in $S_{t}$ (Stock’s number), the volume flow rate decreases in withdrawal case and it increases in the drainage case as shown by Figures 10(b), 11(b), and 12(b).

Figure 13(a) shows that the volume flow rate rises as the fluid is becoming thicker in withdrawal case and it declines in the drainage case, as shown by Figure 13(b).

Figures 14 –16 display the influence of Stock’s number on average film velocity. Figures 14(a), 15(a), and 16(a) indicate that with the increase of Stock’s number, average film velocity decreases in withdrawal case and it increases in the drainage case as shown by Figures 14(b), 15(b), and 16(b).

Figures 17 and 18 present the effects of $H$ and $q$ on average film velocity. With the increase in $H$ , average film velocity decreases as shown by Figure 17(a) in withdrawal case, whereas in the drainage case, average film velocity increases and is presented by Figure 17(b). We examined from Figure 18 that average film velocity increases as the fluid is becoming shear thickening in withdrawal case and vice versa in the drainage case.

Figures 19 and 20 show the impacts of parameters $S_{t}$ and $H$ , respectively, on tangential stress. Figures 19(a) and 20(a) for withdrawal case indicate that tangential stress on the cylindrical surface decreases as $S_{t}$ and $H$ increased, and Figures 19(b) and 20(b) show that tangential stress rises with the rise in $S_{t}$ and $H$ for drainage case.

Figures 21 –23 indicate that temperature increases by rising the values of $S_{t}, B_{r}$ , and $H$ for both withdrawal and drainage cases, either the fluid is dilatant or pseudoplastic. Moreover, temperature increases as the fluid is becoming thicker in both withdrawal and drainage cases, and this effect is shown by Figure 22.

Conclusion

The important findings of this article are listed below.

The theoretical study of steady uniform thin film flow of incompressible generalized Maxwell fluid over vertical cylinder under the influence of non-isothermal effects for withdrawal and drainage cases is presented.

The problem is formulated in terms of differential equations along with appropriate initial and boundary conditions.

Exact methodical solutions are reported for velocity, volume flow rate, average film velocity, tangential stress, and temperature for both cases by applying binomial series technique which is the most efficient, effective, and reliable technique for the solutions of such problems.

Solutions for the Newtonian fluid are recovered as a particular case of our obtained solutions for generalized Maxwell fluid.

Tangential stress on cylindrical surface is similar for Maxwell and Newtonian fluids in withdrawal case. Moreover, the effect of tangential stress also remains invariable in drainage case for either type of the fluid.

Solutions for temperature profile in both cases remain the same for either type of the fluid.

The impacts of different dimensionless parameters on velocity field, volume flow rate, average film velocity, tangential stress, and temperature profile are examined and discussed graphically for both generalized Maxwell and Newtonian fluids.

The temperature of the fluid rises with a rise in the values of pertinent parameters used in the planed problem. This effect shows that there is a direct relation between these parameters and temperature profile. On the other hand, the magnitude of velocity decreases by rising the values of parameters in withdrawal case, while it increases in drainage case.

Footnotes

Acknowledgements

The authors express their heartfelt gratitude to the editor and reviewers for their precious comments, which have led to the improvement in clarity and quality of this paper.

Handling Editor: James Baldwin

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

ORCID iDs

Saif Ullah

Kareem Akhtar

References

Makris

Dargush

Constantinou

MC.

Dynamical analysis of generalized viscoelastic fluids. J Eng Mech 1993; 119: 1663–1679.

Liu

Grecov

Rheological and flow modelling of viscoelastic fluids between eccentric cylinders. Appl Math Model 2011; 35: 1603–1615.

Roylance

Engineering viscoelasticity. Cambridge, MA: Massachusetts Institute of Technology, 2001.

Maxwell

JC.

On the dynamical theory of gases. Phil Trans R Soc Lond A 1866; 157: 26–78.

Yang

Zhu

KQ.

Start-up flow of a viscoelastic fluid in a pipe with the fractional Maxwell model. Comput Math Appl 2010; 60: 2231–2238.

Shahid

Bhatti

Anwar Beg

, et al. Numerical study of radiative Maxwell viscoelastic magnetized flow from a stretching permeable sheet with the Cattaneo-Christov heat flux model. Neural Comput Appl 2018; 30: 3467–3478.

Rahbari

Abbasi

Rahimipetroudi

, et al. Heat transfer and MHD flow of non-Newtonian Maxwell fluid through a parallel plate channel: analytical and numerical solution. Mech Sci 2018; 9: 61–70.

Ersoy

HV.

Flow of a Maxwell fluid between two porous disks rotating about non-coincident axes. Adv Mech Eng 2014, https://journals-sagepub-com-s.web.bisu.edu.cn/doi/full/10.1155/2014/347196.

Rajagopal

Srinivasa

A thermodynamic framework for rate type fluid models. J Non-Newton Fluid Mech 2000; 88: 207–228.

10.

Sylvester

Tyler

Skelland

AHP

. Non-Newtonian film fluids: theory and experiment. Can J Chem Eng 1973; 51: 418–429.

11.

Hameed

Ellahi

Thin film flow of non-Newtonian MHD fluid on a vertically moving belt. Int J Numer Method Fluid 2011; 66: 1409–1419.

12.

Alam

Rahim

Islam

, et al. Adomian decomposition method for thin film flow of a non-Newtonian fluid on a vertical cylinder. J Frank Ins 2013; 350: 818–839.

13.

Anderson

Shang

DY.

An extended study of the hydrodynamics of gravity-driven film flow of Power-law fluids. Fluid Dyn Res 1998; 22: 345–357.

14.

Nadeem

Awais

Thin film flow of an unsteady shrinking sheet through porous medium with variable viscosity. Phys Lett A 2008; 372: 4965–4972.

15.

Asghar

Hayat

Kara

AH.

Exact solutions of thin film flows. Nonlin Dyn 2007; 50: 229–233.

16.

Khan

Islam

Gul

, et al. Thin film flow of a second grade fluid in a porous medium past a stretching sheet with heat transfer. Alex Eng J 2018; 57: 1019–1031.

17.

Wang

Pop

Analysis of the flow of a power-law fluid film on an unsteady stretching surface by means of homotopy analysis method. J Non-Newton Fluid Mech 2016; 138: 161–172.

18.

Faraz

Khan

Thin film flow of an unsteady Maxwell fluid over a shrinking/stretching sheet with variable fluid properties. Int J Numer Method Heat Fluid Flow 2018; 28: 1596–1612.

19.

Shah

Bonyah

Islam

, et al. Radiative MHD thin film flow of Williamson fluid over an unsteady permeable stretching sheet. Heliyon 2018; 4: 1–20.

20.

Landau

Lifshitz

EM.

Fluid mechanics. New York: Pergamon, 1987.

21.

Marinca

Herisanu

Nemes

Optimal homotopy asymptotic method with application to thin film flow. Cent Eur J Phys 2008; 6: 648–653.

22.

Memon

Islam

Siddiqui

, et al. Lift and drainage of electrically conducting power law fluid on a vertical cylinder. Appl Math Inf Sci 2014; 8: 45–55.

23.

Siddiqui

Akram

Memon

, et al. Withdrawal and drainage of thin film flow on a vertical cylinder. Sci Res Essay 2012; 41: 3554–3565.

24.

Gul

Islam

Shah

, et al. Analysis of thin film flow over a vertical oscillating belt with a second grade fluid. Eng Sci Tech Int J 2015; 18: 207–217.

25.

Alam

Islam

Thin film flow of magnetohydrodynamic (MHD) Johnson-Segalman fluid on vertical surfaces using the Adomian decomposition method. Appl Math Comput 2012; 219: 3956–3974.

26.

Ali

Islam

Gul

, et al. Magnetohydrodynamics thin film fluid flow under the effect of thermophoresis and variable fluid properties. Am Inst Chem Eng J 2017; 63: 5149–5158.

27.

Khan

Nie

Shah

Impact of thermal radiation and heat source/sink on MHD time-dependent thin-film flow of Oldroyed-B, Maxwell, and Jeffry fluids over a stretching surface. Processes 2019; 7: 1–18.

28.

Tahir

Gul

Islam

, et al. Flow of a nano-liquid film of Maxwell fluid with thermal radiation and magneto hydrodynamic properties on an unstable stretching sheet. Journal of Nanofluids 2017; 6: 1021–1030.

29.

Hayat

Sajid

On analytic solution for thin film flow of a fourth grade fluid down a vertical cylinder. Phys Lett A 2007; 361: 316–322.

30.

Hayat

Sajid

Thin film flow of an Oldroyd 8-constant fluid: an exact solution. Phys Lett A 2008; 372: 1827–1830.

31.

Sajid

Ali

Hayat

On exact solutions for thin film flows of a micropolar fluid. Commun Nonlin Sci Numer Simul 2009; 14: 451–461.

32.

Siddiqui

Mahmood

Ghori

QK.

Homotopy perturbation method for thin film flow of a third grade fluid down an inclined plane. Chaos Soliton Fract 2008; 35: 140–147.

33.

Siddiqui

Mahmood

Ghori

QK.

Homotopy perturbation method for thin film flow of a fourth grade fluid down a vertical cylinder. Phys Lett A 2006; 352: 404–410.

34.

Siddiqui

Ahmed

Ghori

QK.

Thin film flow of non-Newtonian fluids on a moving belt. Chaos Soliton Fract 2007; 33: 1006–1016.

35.

Ullah

. Withdrawal and drainage of generalized Oldroyd-B fluid on cylindrical surfaces. MPhil Thesis, 2015. Lahore, Pakistan: Department of Mathematics, Government College University.

36.

Khan

Mahmood

The influence of slip condition on the thin film flow of a third order fluid. Int J Nonlin Sci 2012; 13: 105–116.

37.

Azaz-ul-Haq

. Withdrawal and drainage of generalized Maxwell fluid on cylindrical surfaces. MPhil Thesis, 2015. Lahore, Pakistan: Department of Mathematics, Government College University.

38.

Farooq

Rahim

Islam

, et al. Series solutions of lifting and drainage problems of a nonisothermal modified second grade fluid using a vertical cylinder. J Appl Math 2014; 2014: 842650.

39.

Farooq

Rahim

Islam

, et al. Withdrawal and drainage of generalized second grade fluid on vertical cylinder with slip conditions. J Prime Res Math 2013; 9: 51–64.

40.

Gul

Mubashir Khan

, et al. Influence of slip condition on MHD thin film flow of a third grade fluid over a vertical belt with temperature dependent viscosity. J Appl Environ Biol Sci 2015; 5: 22–30.

41.

Gul

Shah

Islam

, et al. MHD thin film flow of a third grade fluid on a vertical belt with slip boundary conditions. J Appl Math 2013; 2013: 707286.

42.

Ghani

Gul

Islam

, et al. Unsteady magnetohydrodynamics thin film flow of a third grade fluid over an oscillating inclined belt embedded in a porous medium. Therm Sci 2017; 21: 875–887.

43.

Rivlin

Ericksen

JL.

Stress-deformation relations for isotropic materials. J Rational Mech Anal 1995; 4: 323–425.