Document Type : Full Lenght Research Article
Authors
^{1} University of Gour Banga, Malda 732 103, WB, India
^{2} University of Gour Banga, Malda 732 103, India
^{3} Vidyasagar University, Midnapore 721 102, India
Abstract
Keywords
Main Subjects
History: Submitted 1 January 2016 Revised 3 April 2016 Accepted 22 May 2016
Keywords: slip flow couple stress fluid stretching sheet thermal radiation viscous dissipation 
Slip flow of an optically thin radiating nongray couple stress fluid past a stretching sheet
Sanatan Das*^{1}, Akram Ali^{1} and Rabindra Nath Jana^{2}
^{1}Department of Mathematics, University of Gour Banga, Malda 732 103, India ^{2}Department of Applied Mathematics, Vidyasagar University, Midnapore 721 102, India ^{ }

A B S T R A C T This study addresses the combined effects of couple stresses, thermal radiation, viscous dissipation and slip condition on a free convective flow of a viscous incompressible couple stress fluid induced by a vertical stretching sheet. Cogley VincentiGilles equilibrium model is employed to include the effects of thermal radiation in the study. The governing boundary layer equations are transformed into a system of nonlinear differential equations, and solved numerically using the RungeKutta fourthorder method with shooting technique. Numerical results are obtained for the fluid velocity, temperature as well as the shear stress and rate of heat transfer. Then, the effects of the pertinent parameters on these quantities are examined. It is found that both the fluid velocity and temperature reduce in the presence of thermal radiation. Increasing the values of the couple stress parameter thickens the momentum boundary layer. The slip parameter greatly influences the fluid flow and shear stress on the surface of the stretching sheet.
© 2016 Published by Semnan University Press. All rights reserved. 

Journal of Heat and Mass Transfer Research
Journal homepage: http://jhmtr.journals.semnan.ac.ir 

Address of correspondence author: Sanatan Das, University of Gour Banga, Malda 732 103, WB, India Email: tutusanasd@yahoo.co.in

The boundary layer flow over a stretching sheet in a uniform stream of fluid has been studied extensively in fluid mechanics. A large number of studies has been done in the area of boundary layer flow over a continuous stretching surface with regard to its numerous industrial and engineering applications. It occurs frequently in manufacturing involving hot metal rolling, wire drawing, producing glassfiber, producing paper, drawing of plastic films, manufacturing foods, crystal growing, liquid films in condensation process and metal spinning, as well as metal and polymer extrusion processes (Fisher [3]). The flow and heat transfer over a stretching sheet have tremendous industrial applications, for instance, in the extrusion of a polymer sheet from a die. In the process of manufacturing such sheets, the melt issues from a slit are subsequently stretched. The rates of stretching and cooling have a significant influence on the quality of the final product with desired characteristics. The aforementioned processes involve cooling of a molten liquid by drawing it into a cooling system. The desired properties for the product of such process mainly depend on two characteristics:the first is the cooling liquid used and the other is the rate of stretching. The liquids of nonNewtonian characteristics can be chosen for as a cooling liquid as their flow and hence the heat transfer rate can be regulated through some external means. The optimal rate of stretching is important, as rapid stretching results in sudden solidification, thereby destroying the properties expected from the product.
After the pioneering study of Sakiadis [4] and Crane [5], a large number of studies on a stretching sheet have been published by considering various governing parameters such as suction/injection, porosity, magnetic field parameter, and radiation with different types of fluids such as Newtonian, nonNewtonian, polar, and couple stress fluids. However, the abundant literature on the boundary layer flow over a stretching sheet is limited to Newtonian and some nonNewtonian fluids flow with traditional noslip flow boundary condition over various stretching geometries such as linear and nonlinear stretching sheet and a little attention is given to stretching sheet with slip boundary condition. However, the fluids with microscale or nanoscale dimensions have flow behavior that greatly differs from the traditional fluid flow characteristics and belongs to the slip flow regime. For the flow in the slip regime, the fluid motion still obeys the NavierStoke's equations, but with slip velocity, temperature and concentration boundary conditions. For instance, the flow in many applications of micro/nano systems such as hard disk drive, micropump, microvalve and micronozzles is in slip transition regime, which is characterized by slip boundary at the wall.
Several researchers have considered various aspects of momentum and heat transfer characteristics in boundary layer flow over a stretching boundary (Gupta and Gupta [6], Rajagopal et al. [7], Siddappa and Abel [8], Andersson [9], Kumaran and Ramanaiah [10], Wang [11] and Cortell [12]). Flow of elasticoviscous fluid induced by a stretching sheet with slip has been presented by Ariel et al. [13]. Akyildiz et al. [14] have discussed the diffusion of chemically reactive species in a porous medium over a stretching sheet. Wang [15] have made an analysis of viscous flow due to a stretching sheet with surface slip and suction. Fang et al. [1617] have obtained the exact solution of an MHD slip flow over a stretching sheet. The thermal boundary layers over a shrinking sheet have been analytically studied by Fang and Zhang [18]. The heat transfer characteristics in a viscoelastic boundary layer flow over a stretching sheet have been analyzed by Arnold el al.[19]. Shantha et al. [20] have presented the free convection flow of a conducting couple stress fluid in a porous medium. Srinivasacharya and Kaladhar [21] have conducted a study of a mixed convection flow of couple stress fluid in a nondarcy porous medium with Soret and Dufour effects. The secondorder slip flow and heat transfer over a stretching sheet with nonlinear Navier boundary condition have been numerically investigated by Nandeppanavar et al.[22]. Singh and Makinde [23] have studied an MHD slip flow of viscous fluid over an isothermal reactive stretching sheet. Hayat et al. [24] have presented the stagnationpoint flow of couple stress fluid with melting heat transfer. Turkyilmazoglu [25] has obtained an exact solution for twodimensional laminar flow over a continuously stretching or shrinking sheet in an electrically conducting quiescent couple stress fluid. The MHD flow and heat transfer of an exponential stretching sheet in a BoussinesqStokes suspension have been examined by Siddheshwar et al. [26]. Salema et al. [27] have studied the hydromagnetic flow of Cu water nanofluid past a moving wedge with viscous dissipation. Zhu et al. [28] have presented the secondorder slip MHD flow and heat transfer of nanofluids with thermal radiation and chemical reaction. Sheikholeslami et al. [29] have investigated the magnetic nanofluid forced convective heat transfer with variable magnetic field using twophase model. The effect of nonuniform magnetic field on forced convection heat transfer of Fe O  water nanofluid has been examined by Sheikholeslami et al. [30]. Kandelousi [31] has studied the effect of spatially variable magnetic field on ferrofluid flow and heat transfer considering constant heat flux boundary condition. Sheikholeslami and Ganji [32] have presented the nanofluid flow and heat transfer between parallel plates considering Brownian motion using DTM. The effect of thermal radiation on a magnetohydrodynamic nanofluid flow and heat transfer has been examined by Sheikholeslamia et al. [33]. Sheikholeslami [34] has presented the forced convective heat transfer in a semi annulus under the influence of a variable magnetic field.
In the present study, the free convective slip flow induced by a vertical linearly stretching sheet in couple stress fluid is presented in the presence of thermal radiation by adopting the Cogley VincentiGilles equilibrium model introduced by Cogley et al. [35]. The sheet surface is isothermal. The partial differential equations governing the flow are transformed into ordinary differential equations using the similarity transformation and are numerically solved. The interplay between pertinent parameters is explained using graphs.
1. Mathematical Formulation
A steady slip flow of a viscous incompressible couple stress fluid over an isothermal vertical linearly stretching sheet is considered. A Cartesian coordinate system with xaxis along the sheet is chosen and the yaxis is normal to it (see Fig.1). The sheet stretches in its own plane with the velocity, being a constant known as stretching rate. The sheet surface is heated by the temperature T_{w} while the temperature of the ambient cold fluid isT_{∞}. All the fluid properties are assumed to be constant except the density in the buoyancy term of the momentum equation.
The Stokes constitutive model for the couple stresses is used. Under the usual Boussinesq's approximation, the free convection flow of a radiating couple stress fluid is governed by the following system of equations:
(1) 

(2) 

(3) 
Where and are the velocity components along the x and ydirections, respectively, T is the temperature of the fluid, gis the acceleration due to gravity, μ: the dynamic viscosity of the fluid, : the kinematic viscosity, : the couple stress viscosity, : the fluid density, k: the thermal conductivity, : the thermal expansion coefficient,c_{P}: the specific heat at constant pressure andq_{r}: the radiative heat flux. The asymptotical case as which corresponds to the purely Newtonian case. The boundary conditions arethe following:
Where L denotes the slip coefficient and Tw is the temperature of the sheet. WhenL = 0, the noslip condition is recovered.
It has been shown by Cogley et al. that in the optically thin limit for a nongray gas near equilibrium, the following relation holds:
(4) 
(5) 
Where is the absorption coefficient, is the wave length, is the Planck's function and the subscript indicates that all quantities have been evaluated at the temperature which is the temperature of the stretching sheet. Thus, our study is limited to small difference of surface temperature to the fluid temperature. Greif et al. [36]showthat for an optically thin limit, the fluid does not absorb its own emitted radiation, it means that there is no selfabsorption, but the fluid does absorb radiation emitted by the boundaries. Thetreatments to the radiative heating are either in the limit where photon mean free paths are very small, called optically thick or very long, called optically thin. At high temperature the presence of thermal radiation alters the distribution of temperature in the boundary layer, which in turn affects the heat transfer at the channel walls.
On the use of equation (5), the equation (3) becomes
(6)

Where
(7) 
The continuity equation (1) is satisfied by introducing a stream function such as
(8) 
The following similarity variables are also introduced
(9) 
Where is the independent similarity variable, is the nondimensional stream function and is the nondimensional fluid temperature. Substituting the equation (9) into equations (2) and (6), we obtain the following ordinary differential equations:
(10) 

(11) 
Fig. 1 The geometry of the problem 
Where is the radiation parameter, is the couple stress parameter, is the Eckert number, is the Grashof number that approximates the ratio of the buoyancy force to the viscous force acting on a fluid and is the Prandtl number which measures the ratio of momentum diffusivity to the thermal diffusivity. The corresponding boundary conditions are
(12) 
Where is the velocity slip parameter. The noslip case is recovered for . In theabsence of thermal radiation ( ) and neglecting viscous dissipation ( ) and slip condition ( ), the present problem reduces to the problem considered by Siddheshwar et al. [26]. It is observed that the trends of the velocity profiles are similar in the absence of Chandrasekhar's number ( ) as well as heat source/sink ( ) in the case of the prescribed surface temperature.
2. Numerical Solution
It is wellknown that most boundary layer transmission problems are described by a set of nonlinear partial differential equations. However, due to the strongly nonlinear and unconventional nature of these problems, the solving processes are extraordinarily complex, and the effective solutions are hardly obtained. Recently, various methods have been tried to solve these problems. The nondimensional ordinary differential equations (10) and (11) with boundary conditions (12) are solved numerically. The highly nonlinear momentum boundary layer equation and the thermal boundary layer equation are converted into similarity equations and then solved numerically using the modified fourthorder RungKutta method with shooting technique [37]. RungKutta method is the fourthorder method meaning that the local truncation error is of order whereas the total accumulated error is of order . This is a usual case that, the truncation error also includes discretization error, which is the error that arises from taking the finite number of steps in the computation to approximate an infinite process. The resulting higher order ordinary differential equations are reduced to firstorder differential equations by letting
(13) 
Thus, the corresponding coupled higher order nonlinear differential equations become:
(14) 
With the boundary conditions:
(15) 
Whereb, cand d are unknown which are to be determined as a part of the numerical solution. In the shooting method, the unspecified initial conditions b, cand d are assumed and equation (14) is integrated numerically as an initially valued problem to a given terminal point. The accuracy of the assumed missing initial condition is checked by comparing the calculated value of the dependent variable at the terminal point with its given value there. If a difference exists, the improved values of the missing initial conditions must be obtained and the process is repeated. The numerical computations are done by MATLAB bvp4c routine. The stepsize is taken as . The process is repeated until we get the results correct up to the desired accuracy of level.
3. Results and Discussion
In order to get a physical insight of the problem, we illustrate the effects of the pertinent parameters, namely couple stress parameter , Grashof number , radiation parameterR, Prandtl number , Eckert number and slip parameter on the fluid velocity, temperature, shear stress and rate of heat transfer. Prandtl number ( ) is chosen as the values ranging from which are the most encountered fluids used in industries. presents no Joule and viscous heating. For a large value of slip parameter which corresponds to a very small at the leading edge, the boundary layer assumption is not appropriate, i.e. the boundary layer equations become inaccurate. Also, the Knudsen number is greater than 0.1 for a large and hence the NavierStokes equation fails to model the flow regime. We therefore limit the discussion in this study to a relatively small range of from 0 to 2.5 as this range exhibits the slip flow region. The default values of the other parameters are mentioned in the description of the respective figures.
4.1. Effects of parameters on stream function
Figs.23 present the stream function for several values of Grashof number and Eckert number . It is seen from these figures that the stream function increases with an increase in either or . Figs.46 reveals that the stream function decreases with an increase in either couple stress parameter or radiation parameter or Prandtl number .
4.2. Effects of parameters on velocity profiles
Figs.78 show that the fluid velocity increases with an increase in either Grashof number or Eckert number . Physically speaking, the Grashof number signifies the relative importance of buoyancy force to the viscous hydrodynamic force.
Fig. 2 The nondimensional stream function for different when , , , , and 
Fig. 3 The nondimensional stream function for different when , , , and 
Fig. 4 The nondimensional stream function for different when , , , , and 
Fig.5 The nondimensional stream function for different when , , , and 
Larger Grashof number indicates lower viscous effects in momentum equation. For higher values of , the fluid velocity overshoots the sheet velocity in the regions close to the surface of the sheet. Physically, means heating of the fluid by cooling the boundary surface. It is also revealed that the momentum boundary layer thickness reduces as the Grashof number increases and it increases astheEckert number increases. The Eckert number expresses the relationship between the kinetic energy in the flow and the enthalpy. It embodies the conversion of kinetic energy into internal energy by work done against the viscous fluid stresses. Greater viscous dissipative heat causes a rise in the fluid temperature and hence the fluid velocity enhances. Fig.9 reveals that the fluid velocity decreases near the sheet surface and increases away from the surface of the sheet with an increase in the couple stress parameter . This means that the increasing values of result in thickening of the momentum boundary layer. Expectedly, the couple stresses enhance the fluid's intermolecular cohesion. This increases the fluid's resistance to shear stress, resulting in adecrease in velocity. Fig.10 depicts the effects of the radiation parameter on the velocity profile. The fluid velocity decreases for increasing the values of . The thermal radiation has increased the momentum boundary layer thickness. Fig.11 presents the fluid velocity profiles for several values of Prandtl number . The thickness of the momentum boundary layer is also greater. The reason for such a behavior is that an increase in the Prandtl number is due to an increase in the fluid viscosity, which makes the fluid thick, thus decreasing its velocity. It is seen from Fig.12 that the fluid velocity decreases with an increase in the slip parameter . Physically, when slip occurs, the slipping fluid shows a decrease in the surface skinfriction between the fluid and the stretching sheet because not all the pulling force of the stretching sheet can be transmitted to the fluid. So, increasing the value of decreases the fluid velocity in the region of the boundary layer.
4.3. Effects of parameters on temperature profiles
Figs. 1315 illustrate the temperature field against thecoordinate . The fluid temperature is the highest at the sheet and decreases asymptotically to zero far away from the sheet satisfying the boundary conditions. Figs. 13 and 14 present the profiles of the fluid temperature for several values of radiation parameter and Prandtl number .
Fig.6 The nondimensional stream function for different when , , , and 


Fig.7 The velocity profile for different when , , , and 

Fig.8 The velocity profile for different when , , , and 
Fig. 9 The velocity profile for different when , , , and 
Fig.10: The velocity profile for different when , , , and 
Fig.11 The velocity profile for different when , , , and 
Fig.12 The velocity profile for different when , , , and 
An increase in the thermal radiation causes a decrease in the fluid temperature within the boundary layer and consequently the thermal boundary layer thickness decreases as shown in Fig. 13. Thetemperature decreases with anincrease in the radiation parameter as therelease of heat energy from the flow region is also increased, resulting in adecrease in the fluid temperature. It is seen from Fig.14 that the fluid temperature decreases with an increase in the Prandtl number . From this figure it is noticed that an increase in the Prandtl number results in a decrease in temperature. The reason is that thesmaller values of the Prandtl number are equivalent to increasing thermal conductivity and therefore heat is capable of diffusing away from the heated surface more rapidly than at thehigher values of . Thus, the temperature falls more rapidly for water than for air and electrolyte solution. An increase in the causes a decrease in the thermal boundary layer thickness. Fig. 15 illustrates the influence of Eckert number on temperature in the boundary layer. The fluid temperature increases with an increase in theEckert number . TheEckert number is the ratio of the kinetic energy of the flow to the boundary layer enthalpy differences. It represents the conversion of the kinetic energy into the internal energy by work done against the viscous fluid stresses. The positive Eckert number means cooling of the surface of the stretching sheet, i.e loss of heat from the surface of the stretching sheet to the fluid. Hence, greater viscous dissipative heat causes a rise in the fluid temperature. Furthermore, it can be observed that the thermal boundary layer thickness becomes thicker for the increased Eckert number. The maximum temperature occurs in the vicinity of the surface of the stretching sheet and then the temperature asymptotically approaches to zero in the freestream region.
4.4. Effects of parameters on shear stress and rate of heat transfer
For engineering purposes, one is usually interested in the values of the shear stress (or skin friction) and the heat transfer rate. The shear stress is an important parameter in the heat transfer studies, since it is directly related to the heat transfer coefficient. The increased shear stress is generally viewed as a disadvantage in the technical applications, while the increased heat transfer can be exploited in some applications. Numerical values of the rate of heat transfer and the shear stress at the surface of the sheet are entered in the Table 1 for several values of , , , , and . The shear stress reduces with an increase in , and while it enhances with an increase in either or or .
This may be attributed to a rise in the velocity gradient at the surface of the sheet due to the buoyancy force. An increase intheshear stress is observed with an increase in the thermal radiation. As expected, a fall in the fluid temperature due to a rise in thethermal radiation enhances the fluid viscosity, leading to an increase in the shear stress. A decrease in shear stress is noticed as Eckert number increases. This may be explained by stating that,anincrease in thefluid temperature due to the viscous heating may decrease thefluid viscosity, since viscosity is temperature dependent. It is also concluded from Table 1 that both friction factor and surface heat flux reduce with an increase in the value of . This is expected, since an increase in the value of decreases thethermal diffusivity of the fluid and it also reduces the magnitude of thefrictional force between the viscous layers. Consequently,the shear stress reduces. Physically, the negative value of means the stretching sheet exerts a drag force on the fluid and the positive value means the opposite.
Fig.13 The temperature profile for different when and 
Fig.14 Thetemperature profile for different when and 
Fig.15 The temperature profile for different when and 
It is seen from the Table 1 that the rate of heat transfer increases for increasing values of or or while it decreases for increasing values of or or . Due to the influence of thermal radiation, the temperature gradient increases which ultimately enhances therate of heat transfer. Therefore, thermal radiation plays an important role in augmenting heat transfer rate or in other words it acts as a heat source. Thus, increasing the couple stresses and buoyancy force implies more intermolecular cohesion which generates more heat leading to an increase in the temperature which in turn causes an increase in the rate of heat transfer. On the other hand, means the heat transfer takes places from the surface of the sheet to fluid.
Table 1. The rate of heat transfer and shear stress at the sheet . 



0.5 
2 
5 
0.72 
0.5 
0.1 
0.72691 
1.55029 
1.0 
2 
5 
0.72 
0.5 
0.1 
0.85015 
1.53812 
1.5 
2 
5 
0.72 
0.5 
0.1 
0.88136 
1.53456 



0.5 
1 




0.70732 
1.27084 
0.5 
1.5 




0.71752 
1.41428 
0.5 
2 




0.72691 
1.55029 





3 



0.82251 
1.54067 


4 



0.77469 
1.54565 


5 



0.72691 
1.55029 






0.72 


0.72691 
1.55029 



2 


0.72016 
1.47688 



3 


0.71505 
1.42013 







0.0 

0.74004 
1.66408 




0.5 

0.72691 
1.55029 




1.0 

0.71367 
1.43867 







0.05 
0.76583 
1.54306 





0.5 
0.51824 
1.58209 





1.5 
0.30328 
1.60229 
4. Conclusion
The boundary layer slip flow of couple stress fluid past a vertical stretching sheet in the presence of thermal radiation has been numerically investigated. Using similarity transformations the governing equations of the problems are transformed into nonlinear ordinary differential equations and solved for local similar solutions using shooting iteration technique together with fourthorder RungeKutta integration. The main results of the present study are summarized as follows:
Acknowledgments
We thank the respected anonymous reviewer/s who contributedtowards improvement of this paper.
References
[1]. Stokes, V. K. (1966). Couple stresses in fluid. The Physics of Fluids 9(9), 17091715 .
[2]. Stokes, V. K. (1984). Theories of Fluids with Microstructure: An Introduction, Springer Verlag, New York.
[3]. Fischer, E.G. (1976). Extrusion of plastics. Wiley, New York, 1976.
[4]. Sakiadis, B.C.(1961). Boundary layer behaviour on continuous solid surfaces: I boundary layer on a continuous flat surface. AICHE J. 7, 221–5.
[5]. Crane, L.(1970). Flow past a stretching plate. Z. Angew. Math. Phys. 21, 6457.
[6]. Gupta, P.S., Gupta, A.S. (1977). Heat and mass transfer on a stretching sheet with suction or blowing. Canad. J. Chem. Engg. 55, 744746.
[7]. Rajagopal, K.R., Na, T.Y., Gupta, A.S.(1984). Flow of viscoelastic fluid due to a stretching sheet. Rheol. Acta. 23, 213215.
[8]. Siddappa, B., Abel, M. S.(1985). Non Newtonian flow past a stretching surface, Z. Angew. Math. Phys. 36, 890892.
[9]. Andersson, H.I.(1992). MHD flow of a viscoelastic fluid past a stretching surface. Acta Mech. 95, 227230.
[10]. Kumaran, V., Ramanaiah, G.(1996) A note on the flow over a stretching sheet. Acta Mech. 116, 229 233.
[11]. Wang, C.Y. (2002). Flow due to a stretching boundary with partial slip. An exact solution of the Navier–Stokes equations. Chem. Eng. Sci. 57, 3745–7.
[12]. Cortell, R. (2007). Viscous flow and heat transfer over a nonlinearly stretching sheet. Appl. Math. Comput. 184 (2), 864873.
[13]. Ariel, P.D., Hayat, T., Asghar, S. (2006). The flow of an elasticoviscous fluid past a stretching sheet with slip. Acta Mech. 187, 29–36.
[14]. Akyildiz, F.T., Bellout, H., Vajravelu, K.(2006). Diffusion of chemically reactive species in a porous medium over a stretching sheet. J. Math. Anal, Appl. 320, 322–39.
[15]. Wang, C.Y.(2009). Analysis of viscous flow due to a stretching sheet with surface slip and suction. Nonlinear Anal Real World Appl. 10, 375–80.
[16]. Fang, T., Zhang, J., Yao, S. (2009). Slip MHD viscous flow over a stretching sheet An exact solution. Comm. Nonlinear Num. Simu. 14, 37313737.
[17]. Fang, T., Zhang, J., Yao, S. (2010). Slip MHD viscous flow over a permeable shrinking sheet. Chinese Phys. Lett. 27, 124702.
[18]. Fang, T.G., Zhang, J.(2010). Thermal boundary layers over a shrinking sheet: an analytical solution. Acta Mech. 209, 325–43.
[19]. Arnold, J.C., Asir, A.A.,Somasundaram, S., Christopher, T. (2010). Heat transfer in a viscoelastic boundary layer flow over a stretching sheet. Int. J. Heat Mass Transfer 53, 1112–8.
[20]. Shantha, G., Shanker, B. (2010). Free convection flow of a conducting couple stress fluid in a porous medium. Int. J. Numer. Methods Heat Fluid Flow 20, 250264.
[21]. Srinivasacharya, D., Kaladhar, K. (2012). Mixed convection flow of couple stress fluid in a nondarcy porous medium with Soret and Dufour effects. J. Appl. Sci. Eng 15, 415 422.
[22]. Nandeppanavar, M.M., Vajravelu, K., Abel, M.S., Siddalingappa, M. N.(2012). Second order slip flow and heat transfer over a stretching sheet with nonlinear Navier boundary condition. Int. J. Therm. Sci. 58, 143–50.
[23]. Singh, G., Makinde, O.D. (2013). MHD slip flow of viscous fluid over an isothermal reactive stretching sheet. ANNALS of Faculty Engineering Hunedoara – Int. J. Engn. Tome XI, 4146.
[24]. Hayat, T., Mustafa, M., Iqbal, Z., Alsaedi, A. (2013). Stagnationpoint flow of couple stress fluid with melting heat transfer. Appl. Math. Mech.Eng. Ed 34, 167176.
[25]. Turkyilmazoglu, M. (2014). Exact solutions for twodimensional laminar flow over a continuously stretching or shrinking sheet in an electrically conducting quiescent couple stress fluid. Int.J. Heat Mass Transfer 72, 18.
[26]. Siddheshwar, P.G., Sekhar, G.N., A. S. Chethan, A.S. (2014). MHD Flow and heat transfer of an exponential stretching sheet in a BoussinesqStokes suspension. J. Appl. Fluid Mech. 7(1), 169176.
[27]. Salem, A. M., Ismail, G. Fathy, R. (2014). Hydromagnetic flow of Cu water nanofluid past a moving wedge with viscous dissipation. Chin. Phys. B 23(4), 044402.
[28]. Zhu, J., Zheng, L., Zheng, L., Zhang, X. (2015). Secondorder slip MHD flow and heat transfer of nanofluids with thermal radiation and chemical reaction. Appl. Math. Mech.  Engl. Ed., 36(9), 11311146.
[29]. Sheikholeslami, M., Rashidi, M.M. , Ganji, D.D. (2015). Numerical investigation of magnetic nanofluid forced convective heat transfer in existence of variable magnetic field using two phase model. J. Molecular Liquids 212, 117126.
[30]. Sheikholeslami, M., Rashidi, M.M., Ganji, D.D. (2015). Effect of nonuniform magnetic field on forced convection heat transfer of Fe3O4 water nanofluid. Comput. Methods Appl. Mech. Engrg. 294, 299312
[31]. Kandelousi, M. S. (2014). Effect of spatially variable magnetic field on ferrofluid flow and heat transfer considering constant heat flux boundary condition. The European Phys. J. Plus, 129 248.
[32]. Sheikholeslami, M., Ganji, D.D.(2015). Nanofluid flow and heat transfer between parallel plates considering Brownian motion using DTM. Comput. Methods Appl. Mech. Eng. 283, 651663.
[33]. Sheikholeslamia, M., Ganji, D. D., Javed, M. Y., Ellahi, R. (2015). Effect of thermal radiation on magnetohydrodynamics nanofluid flow and heat transfer by means of two phase model. J. Magnetism and Magnetic Materials 374, 3643.
[34]. Sheikholeslami, M., Vajravelu, K., Rashidi, M. M. (2016). Forced convection heat transfer in a semi annulus under the influence of a variable magnetic field. Int. J. Heat and Mass Transfer 92, 339348
[35]. A.C. Cogley, W.C. Vincentine, S.E. Gilles, A differential approximation for radiative transfer in a nongray gas near equilibrium, AIAA Journal 6 (1968) 551555.
[36]. R. Grief, I. S. Habib, J. C. Lin, Laminar convection of a radiating gas in a vertical channel, J. Fluid Mech. 46 (1970) 513520.
[37]. T. Y. Na, Computational Method in Engineering Boundary Value Problems. Academic Press, New York, 197