The steady two-dimensional MHD mixed convection boundary layer flow and heat transfer of a Jeffrey fluid over an exponentially stretched plate is investigated. The governing partial differential equations are first reduced to nonlinear ordinary differential equations, before being solved numerically using an implicit finite difference scheme. Local similarity solutions are obtained for some embedded parameters, such as Deborah number β, mixed convection parameter λ, Prandtl number Pr and Hartmann number H, are analyzed and discussed.
I. INTRODUCTION
In very-high-power output devices, the forced convection alone is not enough to dissipate all the heat. As such, combining natural convection with forced convection (mixed convection) will often give desired results.1 The phenomenon of mixed convection mainly occurs in many industrial and technical applications, such as the cooling of nuclear reactors during an emergency shutdown, a heat exchanger placed in a low-velocity environment, solar collectors, electronic devices cooled by fans and etc.2 Owing to its great practical applications, it is no wonder the subject has caught researchers’ attention over the last half century, leading to studies in the field involving diverse aspects and physical conditions. Kliegel3 obtained the theoretical solution of heat transfer by mixed convection about a vertical flat plate and Szewczyk4 analyzed the combination of forced and free convection laminar flow. Merkin5 considered the mixed convection problem, taking into account that heat was supplied to the fluid at the wall at a constant rate by using series expansion. Throughout the years, lots of studies have been conducted.6–12
The MHD boundary layer flow over stretching sheet is encountered numerously in many industrial and engineering processes. The quality of the final desired product very much depends on the heat transfer and stretching rate. For example, in metal extrusion processes, the metal sheets need to be “straightened” by stretching the sheets. Drawing the plates in an electrically-conducting fluid subjected to a transverse magnetic field is used to control the cooling rate, which improves the final output.13 The velocity of the stretching sheet may be linear and may not necessarily be linear.14 Nadeem et al.15 studied boundary layer flow of Jeffrey fluid over an exponentially stretching surface with thermal radiation effect, taking into account two cases of heat transfer analysis, i.e prescribed exponential order surface temperature (PEST) and prescribed exponential order heat flux (PEHF). Cortell16 performed numerical analysis of boundary layer flow induced by continuous stretching velocity uw(x) ∼ x1/3 with radiation effect. Mustafa et al.17 addressed a steady flow of Maxwell nanofluid induced by an exponentially stretching sheet subjected to convective heating, and Sharada and Shankar18 studied MHD mixed convection flow of a Casson fluid over an exponentially stretching sheet with the effects of Soret and Dufour, thermal radiation and chemical reaction. There are vast studies on the subject available in the literature.
Hence, motivated by the work of the above researchers, the present work attempts to study the MHD mixed convection Jeffrey fluid boundary layer flow and heat transfer over an exponentially stretching sheet. To the best knowledge of the authors, such a study has not been reported before.
II. ANALYSIS
Consider a steady two-dimensional incompressible Jeffrey fluid flow over an exponentially stretching sheet in the presence of a uniform magnetic field applied normally to the sheet. The x-axis is considered to be taken along the stretching sheet and the y-axis is perpendicular to it. The sheet is assumed to be stretched vertically with velocity and having surface temperature , where uo and To are constant, T∞ is the temperature far away from the stretching sheet and L is the reference length. It should be pointed out that Tw > T∞ and Tw < T∞ corresponds to the assisting and opposing flow, respectively. Under the Boussinesq and boundary layer approximations, the governing boundary layer equations and heat transfer for the problem are
where u and v are the velocity components in the x and y directions, respectively. λ1 is the ratio of the relaxation and retardation times, λ2 is the relaxation time and T is the fluid temperature. is the kinematic viscosity, where μ is the coefficient of fluid viscosity and ρ is the fluid density. g, βT, σ and α are gravitational acceleration, thermal expansion coefficient, electrical conductivity and thermal diffusivity, respectively. Here, B is the magnetic field which is assumed as where Bo is a constant magnetic field. Note that (2) reduced to viscous flow when λ1 = λ2 = 0. The corresponding boundary conditions for the flow problem are
The appropriate similarity variables are:18,19
where η is the similarity variable, f is the dimensionless stream function and θ is the dimensionless temperature. Substituting (5) into (1)–(3), Eq. (1) is automatically satisfied and Eqs. (2)-(3) are reduced to the following nonlinear ordinary differential equations
and the transformed boundary conditions can be written as
Here, prime denotes differentiation with respect to η and is the Deborah number. Since β depends on x, the resultant ordinary differential equation are not pure similarity equations for the whole flow domain but only represent local similarity for any given values of x.15 is the mixed convection parameter with and are the local Grashoff number and local Reynolds number, respectively. is the Hartmann number and is the Prandtl number. It should be pointed out that λ is a constant, with λ > 0 and λ < 0 represent assisting flow (heated plate) and opposing flow (cooled plate), respectively, while λ = 0 corresponds to forced convection and λ≠0 corresponds to mixed convection regime.
From an engineering point of view, there are two important physical quantities of interest; i.e. the skin friction coefficient Cf and the local Nusselt number Nux, which could be obtained once the values of and are found.
III. RESULTS AND DISCUSSION
Eqs. (6) and (7) subject to boundary condition (8) have been solved numerically using the method described in the book by Cebeci and Bradshaw,20 for some arbitrary values of Deborah number β, the mixed convection parameter λ, the Hartmann number H and the Prandtl number Pr. The ratio of the relaxation and retardation times λ1 is held fixed (=0) throughout the study. Here, we used a uniform grid of Δη = 0.01 and found it to be satisfactory for a convergence criterion of 10−5. The boundary layer thickness η chosen starts from 6 up to 15 to ensure the velocity and the temperature profiles approach zero, which validates the boundary condition (8).
The numerical results obtained for and when Pr = 7 for both assisting (λ > 0) and opposing flows (λ < 0) are depicted in figures 1 and 2, respectively. These figures suggest that an assisting buoyancy flow creates an increase in and , while an opposite phenomena is found to occur for the opposing buoyant flow. This is due to the fact that as λ ≫ 0, the flow is assisted by the buoyancy force to increase the velocity, which in turn increases the skin friction at the surface pictured by the term . Further, it should be pointed out that the solutions for both and are limited to certain values of λ when the flow is opposed by the buoyancy force, i.e when λ < 0.
The effect of the Hartmann number H on the and is also depicted on the same figures, i.e figures 1 and 2. It is seen that as H increases, both of and decrease. This is because of the Lorentz force, which acts against the flow and slows down the fluid velocity, decreasing friction at the surface. This is in line with the velocity profile plotted in figure 3. Even though the local heat transfer rate presented by decreases as H increases, it is found that the temperature increases as H increases, which results in a thicker thermal boundary layer, as displayed in figure 4.
Velocity profiles for different values of H and β when λ = 5 and Pr = 7.
Temperature profiles for different values of H and β when λ = 5 and Pr = 7.
Figures 1 and 2 also displayed the effect of the Deborah number β towards the flow and thermal fields. An interesting behaviour is found to happen when the flow is assisted by the buoyancy force, namely at λ = λc. At this particular point, the flow retains similar value of and respectively, regardless whether the flow is Newtonian or non-Newtonian. For λ < λc, the Deborah number is found to increase both and , and beyond this value (λ > λc), an opposite trend is observed. On top of that, the value of λc is found to arise at a higher value of the mixed convection parameter λ when H = 4 compared with H = 1.
The effect of Prandtl number when λ = 5 on the and are depicted in figures 5 and 6, respectively. These figures reveal that an increase in Pr is found to decrease both and . However, a significance drop is more pronounced for the as compared with the . The existence of Pr term in Eq. (5) is the main reason behind such a situation. Additionally, the effects of H and β are more noticeable for the flow field as compared with the thermal field. This is expected as H and β only take place in the flow field given by Eq. (6). As such, the results obtained for are almost similar, even with different values of H and β when Pr is fixed, as depicted in figure 6. Nonetheless, the close up numerical results obtained shows that the is slightly higher when H = 4 compared with H = 1.
IV. CONCLUSION
The present study considered the local similarity solution for MHD mixed convection Jeffrey fluid flow over an exponentially stretching sheet. The numerical results obtained for the and are found to exist for any values of assisting flow (λ > 0) and are limited for the opposing flow (λ < 0). As the Deborah number β increases, an interesting phenomenon occurred at a certain value of λc, mainly when the flow is being assisted.
ACKNOWLEDGMENTS
We would like to express our sincere gratitude to the Ministry of Higher Education, Malaysia for supported the study under the RAGS13-003-0066 grant.