Research Article  Open Access
MHD Thin Film Flows of a Third Grade Fluid on a Vertical Belt with Slip Boundary Conditions
Abstract
The problem of heat transfer analysis is considered in electrically conducting thin film flows with slip boundary conditions. The flow is assumed to be obeying the nonlinear rheological constitutive equation of a third grade fluid. We have solved the governing nonlinear equations of present problems using the traditional Adomian decomposition method (ADM). Particular attention is given to the combined effect of heat and MHD on the velocity field. The results include the profile of velocity, volume flux, skin friction, average velocity, and the temperature distribution across the film. The effects of model parameters on velocity, skin friction and temperature variation have been studied. Optimal homotopy asymptotic method (OHAM) is also used for comparison. The numerical results and absolute errors are derived in tables.
1. Introduction
The flow and heat transfer inside thin films are ubiquitous in civil, environmental sciences, mechanical engineering, biological sciences, geophysics, and elsewhere. This is due to their several applications at large scale such as wire and fiber coating, reactor fluidization, paper production, different food stuffs like ketchup, sauce, and honey, transpiration cooling, gaseous diffusion, drilling mud, oil wills, heat pipes, and fluid cells. The problem of chambers for chemical and biological detection systems like fluid of many chemicals was considered by Lavrik et al. [1]. Heat transfer inside thin film flow with variable pressure was discussed by Khaled and Vafai [2]. Thin film unsteady flow with variable viscosity was investigated by Nadeem and Awais [3]. Munson and Young [4] discussed the thin film flow of Newtonian fluids. Alam et al. investigated the thin film flow of JohnsonSegalman fluids for lifting and drainage problems [5]. The MHD thin film flow was discussed by Hameed and Ellahi for a nonNewtonian fluid on a vertical moving belt [6]. Due to complexity of nonNewtonian fluids, it becomes difficult to suggest a single model which exhibits all properties of nonNewtonian fluids; therefore various empirical and semiempirical models have been imposed. The use of heat transfer together with the MHD and nonNewtonian fluids under the influence of slip boundary conditions is of particular interest in chemical processing. Relevant and interesting work about present work may be found in [7, 8]. One of the established models amongst nonNewtonian fluids is class of third grade fluids which have their constitutive equations based on strong theoretical foundations, where the relation between the stress and strain is not linear. To solve real world problems, different approximate techniques have been used in mathematics, fluid mechanics, and engineering sciences [9]. Some of the common methods are VIM [10–14], DTM [15, 16], HPM [17–19], HAM, and OHAM [20–24]. These methods deal with the nonlinear problems effectively. The work under various configurations on the thin film flows was discussed by Siddiqui et al. [25]. In [26] they examined the thin film flows of Sisko and Oldroyd6 constant fluids on a moving belt. The heat transfer analysis of thin film is also discussed by Chakraborty and Som [27]. The main aim of the present work is to study heat transfer into a thin film of a third grade fluid on a vertical belt under the influence of transverse magnetic field with slip boundary conditions using ADM. In 1992 Adomian [28, 29] introduced the ADM for the approximate solutions for linear and nonlinear problems. Wazwaz [30, 31] used ADM for the reliable treatment of Bratutype and EmdenFowler equations.
2. Basic Equations
The governing equations of an incompressible isothermal and electrically conducting third grade fluid are where is the constant density, is body force per unit mass, is velocity vector of the fluid, is temperature, is thermal conductivity, is specific heat, , denote material time derivative, is Cauchy stress tensor, and is the shear stress. A uniform magnetic field is imposed transversely on the belt. The Lorentz force per unit volume is given by
Shear stress tensor is given by where denotes spherical stress and is defined as
Here and are the material constants, and , , and are the kinematical tensors given by
3. Formulation of the Lift Problem
Consider a wide flat belt moving vertically upward at constant speed , through a large bath of third grade liquid as shown in Figure 1. The belt carries with it a layer of liquid of constant thickness, . For analysis coordinate system is chosen, in which the xaxis is taken parallel to the surface of the belt and yaxis is perpendicular to the belt. Uniform magnetic field is applied transversely to the belt. Assume the flow is steady and laminar after a small distance above the liquid surface layer. The external pressure is atmospheric everywhere.
Velocity and temperature fields are
Boundary conditions are
Inserting the velocity field given in (9) in (1) and (5)–(7), the continuity equation (1) satisfies identically and (5) gives the following components of stress tensor
Using (11), the momentum and energy equations reduce to
Introducing the following nondimensional variables: where is the gravitational parameter, is magnetic parameter, is nonNewtonian effect, is slip parameter, is heat dimensionless number, is the local Reynolds number, and is the nondimensional variable using the above dimensionless variables in (10) and in (12) and dropping bars we obtain
4. Adomian Decomposition Method (ADM)
A general description of the method is as follows. Begin with an equation , where represents a general nonlinear ordinary differential operator involving both linear and nonlinear terms and is a source term. The linear term is decomposed into , where is easily invertible and is the remainder of the linear operator. For convenience, may be taken as the highest order derivative which avoids difficult integrations which result when complicated Green’s functions are involved. Thus the equation may be written as follows: Solving for , where is nonlinear operator and is a source term, since is easily invertible.
Equation (20) can be written as follows: We use depending on the order of the differential equation. If differential equation is an initialvalue problem, if it is desired for boundary value problem as well, integrations are used. The constants of integration are evaluated from the given initial and boundary conditions. can also be treated as definite integral from ( to ). Solving (21) for we obtained
The nonlinear term , is defined as follows: Here are special polynomials called Adomian polynomials and will be equated to . The initial velocity is identified as follows: Comparison for different components of velocity profile is as follows: The Adomian polynomials depend on the velocity components , which play a flourishing role in the rapid convergence of the series. In the above series depends only on , depends on and , depends on , , and , and so forth. Relevant discussion about ADM can be seen in [16, 17].
4.1. The ADM Solution of Lifting Problem
Using the inverse operator of the Adomian decomposition method on (14), we obtained For series solution we use summation on (26) Adomian polynomials are defined from (27) as follows:
From (28) when , the Adomian polynomials in components form are The series solution becomes
The velocity components are obtained by comparing both sides of (33).
4.1.1. Zero Component Problem
Consider From (16) the boundary condition, for , is By making use of (32) in (31), after simplification we obtain
4.1.2. First Component Problem
Consider From (17) the boundary condition, for , is
The solution is
4.1.3. Second Component Problem
Consider From (17) the boundary condition, for , is
The solution is
The series solution up to the second component is Using (33), (36), and (39) in (40), we have The constants () and are given in the appendix.
The dimensionless shear stress is The shear rate becomes
The coefficient of skin friction is defined as follows:
By use of (43), we have
5. Volume Flow Rate and Average Velocity
Volume flow rate in nondimensional form is as follows:
Using velocity field from (41), volume flow rate is obtained as follows.
Average velocity in dimensionless form is given by
6. Temperature Distribution in Case of Lift Problem
On substituting the series solution for velocity field given in (41) in (15) and solving corresponding to the boundary condition given in (18) for fixed values of ; ; ; ; ; , we obtained
7. Formulation of Drainage Problem
The geometry and assumptions of the problem are the same as those in the previous problem. Consider a film of nonNewtonian liquid draining at volume flow rate down the vertical belt, as shown in Figure 2. The belt is stationary and the fluid drain down the belt due to gravity. The coordinate system is selected in the same way as that in the previous case. Assuming the flow is steady, laminar and external pressure is neglected. Consider that fluid shear forces balance gravity and the film thickness remain constant.
Boundary condition for electrically conducting drainage problem is as follows: Using nondimensional variables the slip boundary conditions for drainage problem become
7.1. Solution of the Drainage Problem by ADM
Using ADM on (14), the Adomian polynomials in (29) for both problems are the same. The different velocity components are obtained as follows.
7.1.1. Zero Component Problem
Consider The solution is
7.1.2. First Component Problem
Consider For different velocity components, the Adomian polynomials mentioned in (29) are used:
7.1.3. Second Component Problem
Consider The solution is The series solution of the velocity field is Substituting (53), (55), and (57) in (58), The constants () and are defined and listed in the appendix.
Substituting (59) the shear rate is obtained as
8. Volume Flow Rate and Average Velocity of Thin Film Flow
Using (59) in (47) we obtain the following.
The average velocity is defined as
9. Heat Distribution for Drainage Problem
Solving (15) with boundary conditions given in (18), after making use of (59) for fixed values of , , , , and , we obtain
10. Results and Discussion
The effects of magnetic parameter , nonNewtonian parameter , gravitational parameter , slip parameter , and the dimensionless number , for both lift and drainage velocity profiles, are discussed in Figures 5–14. Figures 1 and 2 show the geometry of lift and drainage velocity profiles. Figures 3 and 4 show the comparison of OHAM and ADM for lift and drainage velocity profiles. Numerical results and absolute error for both problems are shown in Tables 1 and 2, respectively. Figures 5 and 10 show that the rise in the nonNewtonian parameter increases the speed of flow. For small values of , the velocity profile differs little from the Newtonian one; however when is increased, these profiles become more flattened showing the shearthinning effect. Behavior of the velocity field , for different values of , by fixing other physical parameters, is shown in Figures 6 and 11 for lifting and drainage of fluid, respectively. Here, it can be seen that the boundary layer thickness is reciprocal to the transverse magnetic field and the velocity decreases as one progresses towards the surface of the fluid. We note that the velocity of fluid is maximum at the surface of the belt and minimum at the surface. Moreover, it is to be noted that for large values of , the velocity increases rapidly as compared to small values. When the gravitational parameter increases, the velocity decreases in lifting flow and increases in case of drainage. This can be seen from Figures 7 and 12, respectively. Due to friction force the gravitational effect seems to be smaller near the belt. It can be seen that there is a point in the domain, where the velocity of the fluid becomes approximately the same for different values of gravitational parameter. The reason is that the friction of the belt becomes negligible at this point. On increasing after this point in lifting flow, the velocity decreases due to negligible friction but in drainage flow after this point, the velocity of the fluid increases. The effect of slip parameter can be observed from Figures 8 and 13. It is noticed that the speed of the fluid near the belt is greater than the speed at the surface. When we increase the slip parameter, the velocity of the fluid increases and comparatively this increase can be seen more clearly between and because the friction goes on decreasing. Figures 9 and 14 indicate the dimensionless temperature distribution for different values of . It can be seen that temperature distribution increases as the dimensionless parameter increases and becomes more flattened for large values of . Figure 15 shows the effect of local Reynolds number versus skin friction. This figure shows that the Reynolds number decreases the skin friction. For large values of Reynolds number, the skin friction vanishes. Figure 16 shows the effect of nonNewtonian parameter versus skin friction.


11. Conclusion
The constitutive equation governing the flow of a third grade fluid for lifting and drainage of fluid with slip boundary conditions is solved analytically by using Adomian decomposition method. Expression for velocity field, volume flow rate, skin friction, and temperature distribution is derived and sketched. It is concluded that velocity increases as the gravitational parameter decreases in lifting case while velocity increases as this parameter increases in drainage. For small values of , the velocity profile differs little from the Newtonian one; however when is increased, these profiles become more flattened showing the shearthinning effect. It can be seen that the boundary layer thickness is reciprocal to the transverse magnetic effect and the velocity decreases as it progresses towards the surface of the fluid. On increasing the slip parameter, the velocity of the fluid increases and comparatively this increase can be seen more clearly between and because the friction goes on decreasing. According to the best of our knowledge there is no previous literature about discussed problem; this is our first attempt to handle this problem with slip boundary condition. Also this problem is more general when compared to linear viscous model and second grade fluid model.