Cite this Article
Mehmood Zaffar, Mehmood R., Iqbal Z.. Numerical Investigation of Micropolar Casson Fluid over a Stretching Sheet with Internal Heating. Communications in Theoretical Physics, 2017, 67(4): 443  
Permissions
Numerical Investigation of Micropolar Casson Fluid over a Stretching Sheet with Internal Heating
Mehmood Zaffar *, Mehmood R., Iqbal Z.
Department of Mathematics, HITEC University, Taxila, Pakistan

 

† Corresponding author. E-mail: zaffarilyas@gmail.com

Abstract
Abstract

This theoretical study investigates the microrotation effects on mixed convection flow induced by a stretching sheet. Casson fluid model along with microrotation is considered to model the governing flow problem. The system is assumed to undergo internal heating phenomenon. The governing physical problem is transformed into system of nonlinear ordinary differential equations using scaling group of transformations. These equations are solved numerically using Runge Kutta Fehlberg scheme coupled with shooting technique. Influence of sundry parameters for the case of strong and weak concentration of microelements on velocity, temperature, skin friction and local heat flux at the surface are computed and discussed. Lower skin friction and heat flux is observed for the case of weak concentration () compared to strong concentration of microelements () near the wall.

PACS: 44.25.+f;47.10.ad;47.50.-d
Keyword:micropolar;Casson fluid;mixed convection;internal heating;numerical solution
1 Introduction

Fluids with microstructures are termed as micropolar fluids. These fluids comprise of rigid, randomly oriented particles submerged in a glutinous medium. Micropolar fluids find tremendous applications in blood, foodstuffs, polymers, liquid metal and alloys, plasma and drilling of oil and gas wells etc. Such type of fluid model contains non-symmetric stress tensors. Eringen[1–2] introduced the theoretical explanations of micropolar fluids and discovered the effects of micro motion of fluid elements. He proposed a logical and significant overview of the classical Navier–Stokes model, covering, both in theory and applications, many more phenomena than the classical one. Moreover his generalization was well-designed and not too complex. Airman et al.[3,4] presented a detailed review on application of fluids experiencing micro rotation at particle level. Khonsari and Brewe[5] examined the effects of viscous dissipation on lubrication characteristics of micropolar fluids. They reported that existence of microstructure, according to the micropolar theory, tends to enhance the load-carrying capacity and friction coefficient. Rotating micropolar fluid between parallel plates with heat transfer under the influence of transverse magnetic field has been investigated by Rashid et al.[6] They concluded that the micro rotation is an increasing function of coupling parameter, magnetic field and Reynolds number for strong concentration and it is a decreasing function of viscosity parameter. Nazar et al.[7] analyzed the stagnation point flow of a micropolar fluid over a stretching sheet. They carried out numerical investigation by employing Keller Box method. Non-Newtonian fluids have been an intense topic of research for the past few decades. Much focus had been given to the modeling and analysis of non-Newtonian fluids with rheological characteristics because of usage of various non-Newtonian fluids such as lubricants in industry. The non-linearity can manifest itself in a variety of ways in many fields, such as in food processing, drilling operations and bio-engineering. In this regard Ellahi et al.[8] presented numerical analysis of MHD steady non-Newtonian flows in presence of heat transfer and nonlinear slip effects. Similarly, Makinde et al.[9] analyzed unsteady flow of a reactive variable viscosity non-Newtonian fluid through a porous saturated medium with asymmetric convective boundary conditions and observe that there is a transient increase in both fluid velocity and temperature with an increase in the reaction strength, viscous heating and fluid viscosity parameter. Among the class of several other non-Newtonian fluid models, Casson fluid is one such model with yield stress characteristics. Casson fluid falls in the category of dilatant fluids. It is assumed that Casson fluid has an infinite viscosity at zero shear rates. If the applied shear stress is less than the yield stress then fluid behave like a solid and when shear stress applied is greater than yield stress fluid starts to move. Casson fluid model best fit to rheological data for numerous materials such as jelly, sauce, honey, soup and concentrated fruit juices, etc.[10] The presence of protein, fibrinogen and globulin in aqueous base plasma, red blood cells makes human blood an ideal example of Casson fluid. Number of researchers used Casson fluid to mathematically model and examine the blood flow under a low shear rate in narrow arteries. Shahzad et al.[11] studied effects of mass transfer on generalized non-Newtonian fluid. They analyzed the influence of chemical reaction and suction on Casson fluid in presence of magnetic field. Nadeem et al.[12] observed analytically flow of Casson nanofluid. Mustafa et al.[13] discussed unsteady boundary layer flow of a Casson fluid due to an impulsively started moving flat plate. Similarly Bhattacharya et al.[14] presented analytic solution for magneto hydrodynamic boundary layer flow of Casson fluid over a stretching/shrinking sheet with wall mass transfer. Heat transfer of viscous non-Newtonian fluids past a stretching sheet is a considerable problem in fluid dynamics. The problems arise in the field of engineering and metallurgy depends on hydrodynamic flow and heat transfer rate. In polymer technology, stretching plastic sheets are used in manufacturing products. In electrically conducting fluids, strips are used to control the cooling process. Chen[15] worked out laminar mixed convection of stretching sheet adjacent to vertical wall. Ali et al.[16] examined laminar mixed convection boundary layers induced by a linearly stretching permeable surface. Ishak et al.[17] inspected mixed convection boundary layers in the stagnation-point flow toward a stretching vertical sheet. The thermal effects in fluid flows may cause heat transfer effects in manufacturing processes. In these processes, thermal buoyancy force arises due to heating of a surface that may be in rest or moving continuously under some circumstances. Shateyi et al.[18] analyzed the effects of thermal radiation, hall currents, soret and dufour on MHD flow and heat and mass transfer in a micro polar fluid by mixed convection over stretching surfaces in porous media. Vajravelu et al.[19] studied the heat transfer in a viscous fluid over a stretching sheet with viscous dissipation and internal heat generation. Some more studies related to the current topic can be found in Refs. [20–31].

Most of researches have been carried out in micro polar fluids for characterizing the impact of magnetic field, heat transfer, mixed convection and viscous dissipation etc. To the best of our knowledge, micro polar fluids with rheological characteristics have not been discussed in past. Novelty of present study is to examine micro polar Casson fluid towards a stretching sheet influenced by internal heat generation. Recent research is a fresh contribution in this regard. Effects of sundry parameters on flow and heat transfer characteristics are examined and discussed in a physical manner.

2 Mathematical Formulation

Consider steady 2D flow of an incompressible micropolar Casson fluid towards a linear stretching convective sheet. Heat and mass transfer flow due to stretching of a heated or cooled surface of variable temperature and uniform ambient temperature is is considered. The governing equations of motion (i.e. the continuity, momentum, energy) in vector form for micropolar fluid with rheological characteristics are as follow:

Component form of Eqs. (1)–(4) under boundary layer approximations and considering buoyancy effects are given by:

with appropriate boundary conditions

where in above equations u and v are velocity components along coordinates axes, is velocity at wall, is stretching parameter, ρ is fluid density, ν is kinematic viscosity, k is vortex viscosity, β is Casson fluid parameter, N is micropolar rotation velocity, n is boundary concentration parameter of fluid, the case represents strong concentration, indicates weak concentration, Cp is the specific heat at constant pressure p, κ is thermal conductivity of the medium and T is fluid temperature. A stream of cold fluid at temperature is moving over sheet while the surface of sheet is heated from below by convection from hot fluid at temperature Tf which provides a heat transfer coefficient hf, g0 is acceleration due to gravity, m is coefficient of thermal expansion, Q0 is heat generation coefficient, γ is spin radiation viscosity defined as where is the micro-inertia density.

To convert above system of partial differential equations, we introduce following similarity transformations[32]

Equation (5) is automatically satisfied and Eqs. (6)–(8) become

and corresponding boundary conditions in Eqs. (9) and (10) take the form

where K is micropolar parameter, λ is thermal convective parameter, Pr is the Prandtl number, δ is heat generation parameter and Bi is Biot number which are defined as:

where Gr is Grashof number and is Reynold number are given by the relations

The physical quantities of interest are skin friction coefficient Cf and local Nusselt number Nux which are defined as

where the wall friction τw and heat transfer at wall qw, are expressed as

In view of Eq. (19), expressions described in (18) provide the skin friction and local Nusselt number as

in which is local Reynolds number.

3 The Numerical Solutions

Shooting method along with Runge Kutta fifth order technique was incorporated to tackle the system of nonlinear differential equations. Thus, solution of coupled nonlinear governing boundary layer Eqs. (12)–(14) together with boundary conditions in Eq. (15) are computed by means of shooting method along Runge Kutta fifth order technique. Initially higher order nonlinear differential equations (12)–(14) are converted into a system of first order differential equations and further transformed into initial value problem by labeling the variables as

Associated boundary conditions in Eq. (15) can be transformed as

Above nonlinear coupled ODEs along with initial conditions are solved using Runge Kutta method of order 5 integration techique. Appropriate values of unknown initial conditions S1, S2 and S3 are approximated through Newton’s method. Computations are carried out using mathematics software MATLAB. End of boundary layer region i.e., when to each group of parameters, is determined when the values of unknown boundary conditions at do not change to a successful loop with error less than 10−6 (see Refs. [33–34]).

4 Results and Discussion

This section is dedicated to examine the influence of sundry parameters on velocity , microrotation and temperature profile in the presence of strong and weak concentrations. Figure 1 is plotted to discover influence of micropolar parameter K on velocity profile for the case of weak concentration. It can be observed that velocity profile as well as corresponding momentum boundary layer thickness rises with the increasing behavior of micropolar parameter K for weak concentration. Figure 2 depicts that microrotation profile increases with micropolar parameter K near the wall for both cases of concentration but reverse behavior is observed away from surface. Figure 3 depicts the influence of mixed convection parameter λ on microrotation profile for strong as well as weak concentrations is positive. Influence of Casson fluid parameter β on microrotation profile is presented through Fig. 4. It is quite evident that microrotation profile is higher for the case of weak concentration as compared to strong concentration . The graph for various values of Biot number Bi for temperature profile is displayed in Fig. 5. It can be seen that temperature increases with attractive conduct of Biot number Bi. Figure 6 is plotted to examine the effect of δ on temperature profile for weak concentration . Here a grow in temperature and thermal boundary layer thickness is observed for mount in δ for weak concentration .

Fig. 1 Behavior of velocity profile against .
Fig. 2 Behavior of microrotation profile against .
Fig. 3 Behavior of microrotation profile against .
Fig. 4 Behavior of microrotaion profile aganist β for and .
Fig. 5 Behavior of temperature profile against Biot number Bi.
Fig. 6 Behavior of temperature profile against δ.

To investigate effects of parameters on skin friction coefficient and local Nusselt number we have demonstrated Figs. 7–10. From Fig. 7 it is illustrated that skin friction lessens by rising fluid parameter β for strong as well as weak concentrations. The behavior of micropolar parameter K on skin friction is seen in Fig. 8. It explains that skin friction is a decreasing function of β as micropolar parameter K increases. Moreover, Fig. 9 demonstrates the effect of micropolar parameter K on Nusselt number when plotted against fluid parameter β. Heat flux for weak and strong concentrations is displayed in Fig. 10. It is observed that heat flux as a function of micropolar parameter K falls with growing values of fluid parameter β. Effects for concerning parameters have similar behavior as in skin friction for both cases of concentration.

Fig. 7 Behavior of skin friction coefficient aganist K for different values of .
Fig. 8 Behavior of skin friction coefficient aganist β for different values of .
Fig. 9 Behavior of Nusselt number aganist K for different values of .
Fig. 10 Behavior of skin friction coefficient aganist β for different values of .
5 Concluding Remarks

The aim of this study is to investigate microrotation effects on mixed convective flow of a Casson fluid induced by a stretching sheet. The governing physical problem is tackled numerically using Runge Kutta Fehlberg scheme coupled with shooting. The core outcomes of this study are:

The velocity profile and microrotation profile depict opposite behavior against micropolar parameter K.

Microrotation profile rises with mixed convective parameter λ while it decreases with Casson fluid parameter β for strong as well as weak concentration.

Skin friction coefficient and local Nusselt number rise with micropolar parameter K while decrease with Casson fluid parameter β.

Higher skin friction coefficient and local Nusselt number are observed for the case of strong concentration compared to weak concentration .

Reference
[1] Eringen A. C. Int. J. Engg. Sci. 2 1964 205
[2] Eringen A. C. J. Math. Mech. 16 1966 1
[3] Ariman T. Turk M. A. Sylvester N. D. Int. J. Engg. Sci. 11 1973 905
[4] Ariman T. Turk M. A. Sylvester N. D. Int. J. Engg. Sci. 12 1974 273
[5] Khonsari M. M. Brewe D. E. Acta Mech. 105 1994 57
[6] Mehmood R. Nadeem S. Masood S. J. Magn. Magn. Mat. 401 2016 1006
[7] Nazar R. Amin N. Filip D. Pop I. Int. J. Nonlinear Mech. 39 2004 1227
[8] Ellahi R. Hameed M. Int. J. Num. Meth. Heat Fluid Flow 22 2012 24
[9] Makinde O. D. Chinyoka T. Rundora L. Comp. Math. Appl. 62 2011 3343
[10] Casson N. Rheology of Disperse Systems in Flow Equation for Pigment Oil Suspensions of the Printing Ink Type, Rheology of disperse systems Mill C. C. Pergamon Press London, U.K 1959 84 102
[11] Shahzad S. A. Hayat T. Qasim M. Asghar S. J. Chem. Engg. 30 2013 187
[12] Nadeem S. Mehmood R. Akbar N. S. Int. J. Thermal Sci. 78 2014 90
[13] Mustafa M. Hayat T. Pop I. Aziz A. Heat Transf. –Asian Res. 40 2011 563
[14] Bhattacharyya K. Hayat T. Alsaedi A. Chin. Phys. B 22 2013 024702
[15] Chen C. H. J. Heat Mass Transf. 33 1998 471
[16] Ali M. E. Al-Yousef F. Int. J. Heat Mass Transf. 45 2002 4241
[17] Ishak A. Nazar R. Pop I. Meccanica 41 2006 509
[18] Shateyi S. Motsa S. S. Sibanda P. Math. Prob. Eng. 2010 2010 1
[19] Vajravelu K. Hadjinicolaou A. Int. Commun. Heat Mass Transfer 20 1993 417
[20] Chen C. H. Int. J. Non-Linear Mech. 44 2009 596
[21] Rashidi M. M. Abbasbandy S. Commun. Nonlin. Sci. Num. Simul. 16 2011 1874
[22] Eldabe N. T. M. Salwa M. G. E. J. Phys. Soc. Jpn. 64 1995 41
[23] Nadeem S. Mehmood R. Akbar N. S. J. Comput. Theor. Nano Sci. 11 2014 1422
[24] Pop I. Ingham D. B. Convective heat transfer Mathematical and computational modeling of viscous fluids and porous media pergamon Oxford 2001
[25] Makinde O. D. Int. J. Phys. Sci. 5 2010 700
[26] Ishak A. Appl. Math. Comput. 217 2010 837
[27] Makinde O. D. Olanrewaju P. O. Trans. ASME J. Fluid Eng. 132 2010 044502
[28] Nandeppanavar M. M. Vajravelu K. Subhas Abel M. Siddalingappa M. N. Meccanica 48 2013 1451
[29] Lin T. F. Chang C. J. Yan W. M. J. Heat Transf. 110 1988 337
[30] Karwe M. V. Jaluria Y. Int. J. Heat Mass Transfer 35 1992 493
[31] Yacob N. A. Ishak A. Pop I. Vajravelu K. Nanoscale Research Lett. 6 2011 1
[32] Qasim M. Khan I. Shafie S. Plos ONE 8 2013 49045
[33] Mehmood Z. Iqbal Z. J. Mol. Liq. 224 2016 1083
[34] Azhar E. Iqbal Z. Maraj E. N. Zeitschrift fr Naturforschung A 71 2016 837848