Open access peer-reviewed chapter

Heat Transfer in a MHD Nanofluid Over a Stretching Sheet

Written By

Vikas Poply

Submitted: 26 July 2020 Reviewed: 06 December 2020 Published: 28 December 2020

DOI: 10.5772/intechopen.95387

From the Edited Volume

Heat Transfer

Edited by Miguel Araiz Vega

Chapter metrics overview

459 Chapter Downloads

View Full Metrics


The intention behind carrying out this research work is to analyze the heat transfer characteristics in a Magnetohydrodynamic (MHD) boundary layer nanofluid flow over a stretching sheet. Two phase representation of nanofluid studied the consequence of Brownian motion along with thermophoresis. The major purpose of study is to investigate the significant role of prominent fluid parameters thermophoresis, Brownian motion, Eckert number, Schmidt number and magnetic parameter on profile of velocity, temperature distribution and concentration. Runge–Kutta Fehlberg (RKF) method was adopted to numerically solve the non-linear governing equations and the linked boundary conditions by use of shooting technique. Over all the consequence of prominent fluid parameters are explained via graphs, whereas distinction of several valuable engineering quantities like skin friction coefficient, local Nusselt number and local Sherwood number are also tabulated. The finding of present study helps to control the rate of heat transportation as well as fluid velocity in any manufacturing processes and industrial applications to make desired quality of final product.


  • heat transfer
  • nanofluid
  • stretching sheet
  • MHD flow

1. Introduction

In various fields of science and technology rapid progress has urged the researchers to extend their study towards the regime of boundary layer flow over a stretching sheet. The boundary layer flow behavior towards a linearly or non-linearly stretching sheet plays a significant role for solving engineering problems and possess vast applications in manufacturing and production processes including metal spinning, rubber sheet manufacturing, production of glass fibers, wire drawing, extrusion of polymer sheets, petroleum industries, polymer processing etc. In these cases, the final product of desired characteristics depends on the rate of cooling in the process and the process of stretching. The dynamics of the boundary layer fluid flow over a stretching surface originated from the pioneering work of Crane [1] and he examined the incompressible steady boundary layer fluid flow caused by stretching sheet which moves in its own plane with linear velocity due to the uniform stress applications. This problem is particularly interesting as Crane [1] obtained the exact solution of 2D Navier–Stokes equations. Afterthat, Gupta and Gupta [2] extend Crane [1] work over different mathematical geometries. Yoon et al. [3] studied the theoretical and experimental results using Coulomb friction model by considering punch-sheet interface. Also, Sarma and Rao [4] examined the viscoelastic fluid flow by considering stretched sheet. In view of this, Vajravelu [5] studied flow and heat transfer in a viscous fluid over a nonlinear stretching sheet without using the impact of viscous dissipation. Cortell [6] examined heat and fluid flow transportation over a nonlinear stretching sheet for two different types of thermal boundary conditions, prescribed surface temperature (PST) and constant surface temperature (CST). The influence of heat transfer on the stagnation point flow of a third-order fluid over a shrinking surface has been studied by Nadeem et al. [7]. Recently, Prasad et al. [8] examined the mixed convection heat transfer aspects with variable fluid flow properties over a non-linear stretching surface.

Fluid heating and cooling are important in many industries such as power, manufacturing and transportation. Effective cooling techniques are greatly needed for cooling any sort of high energy device. Common heat transfer fluids such as water, ethylene glycol, and engine oil have limited heat transfer capabilities due to their low heat transfer properties. In contrast, metals have thermal conductivities up to three times higher than these fluids, so it is naturally desirable to combine the two substances to produce a heat transfer medium that behaves like a fluid, but has the thermal of a metal. Since last two decades, study of nanofluid has urged the researcher’s attention due to their heat transportation rate. Nanofluid comes in existence when we add a small quantity of nano-sized109107 particles to the base fluids. Low heat transportation fluids like fluorocarbons, glycol, deionized water, etc. have badly thermal conductivity and therefore deliberated necessary for heat transfer coefficient surrounded by heat transfer medium and surface. The nanoparticles are typically made up of metals (Al,Cu), nitrides (AlN,SiN), carbides (Sic), oxides (Al2O3), or nonmetals (carbon nanotubes, Graphite, etc.) and the base fluid (conductive fluid) is usually water or ethylene glycol. Also, it has been experimentally proved that rate of heat conduction of nanofluids is more than rate of heat conduction of the base fluids. The concept of nanofluid was initially proposed by Choi and Eastman [9] to indicate engineered colloids composed of nanoparticles dispersed in a base fluid. An MIT based comprehensive survey has been done by Buongiorno [10] for convective transportation in nanofluids by considering seven slip conditions that may produce a relative velocity within the base fluid and nanoparticles. Only two (Brownian motion and thermophoresis) out of these seven slip mechanisms were found to be important mechanisms. By adopting Buogiorno’s model, Kuznetsov and Nield [11] explored the nanofluid boundary layer uniform convecting fluid flow.

In recent years, MHD fluid flow has gained researchers attention due to its controllable heat transfer rate. Magnetohydrodynamics (MHD) effect also play and influential role in controlling the rate of cooling as well as segregation of molten metal’s from various non-metallic impurities. Magnetohydrodynamic (MHD) fluid flow has enormous utilization in manufacturing processes, even in the industrial areas as well. The terminology “Magnetohydrodynamic” is combination of three elementary terms magneto that stands for magnetic field, hydro that stands for fluid/liquid and dynamics that stands for evolution of particles. The existence of external magnetic field gives rise to Lorentz drag force which acts on the fluid, so potentially altering the characteristics of fluid flow especially velocity, temperature and concentration. Grouping of electromagnetism Maxwell’s equation and fluid mechanics Navier’s stokes equations therefore provides Magnetohydrodynamic (MHD) relation [12, 13]. Hayat et al. [14] studied the MHD fluid flow transportation over stretching surfaces. Later, the influence of viscous and Ohmic dissipation (i.e. joule heating) in nanofluid has been presented by Hussain et al. [15]. Vajravelu and Canon [16] studied the flow behavior of fluid towards a non-linear stretching sheet. Further, Matin et al. [17] analyzed the entropy effect in MHD nanofluid flow over stretching surface. Shawky et al. [18] studied the Williamson nanofluid flow in porous medium and he acknowledged that enhancement in non-Newtonian parameter escalates skin friction coefficient along with the rate of heat transfer. Basir et al. [19] examined the consequences of Peclet and Schmidt number in existence of partial slip towards a stretching surface. After that, rate of heat transfer along with partial slip condition was generalized by Pandey and Kumar [20]. Recently, Vinita and Poply [21] discussed MHD slip fluid flow of nanofluid in the existence of free stream velocity or outer velocity towards a stretching surface. Vinita et al. [22] studied MHD fluid flow with variable slip conditions over non-linear stretching surface. Furthermore, non-linearity effect towards the stretching surface under different physical circumstances has been examined by researchers in [6, 23, 24, 25].

The outcomes of current study reveal that the results obtained is very significant in the formation of quality object in various manufacturing processes like polymer engineering, paper technology, wire and plastic industries. This study holds important industrial application, particularly in the field of extrusion where the fluid dispersed with particles is used to augment the strength and durability of the material.


2. Materials and methods

In present analysis, 2-D incompressible fluid flow in MHD nanofluid over linear stretching sheet has been considered. Linear behavior generates flow and sheet is stretched in both direction of x axis with stretching velocityuw=ax, wherea and x denotes a constant and stretching surface coordinate respectively. Tw=T+T0xm aty=0, whereT0 refers to the positive constant,T refers to the ambient temperature attained and m refers to the physical parameter known as surface temperature parameter. Also, by introducingm=0, we have a special case of constant surface temperature (CST). Figure 1 represents the physical model of the current study. The continuity, momentum, energy and concentration equations of the incompressible nanofluid boundary layer flow are as follows [10]

Figure 1.

Physical model and coordinate system.


Boundary conditions are given as:


Here horizontal and vertical velocities are represented byu andv, respectively. Alsoν denotes kinematic viscosity,ρ is the density of fluid, τ=ρcpρcf defines a proportion of heat capacities,DT reflects thermophoretic diffusion coefficient,B is the magnetic field intensity,DB denotes Brownian diffusion coefficient,σ represents electrical conductivity

The fundamental Eqs. (1)-(4) with boundary conditions (5) and (6) are transformed using similarity variables


Inserting Eq. (7) into Eqs. (2)(4), the governing Eqs. (1)(4) takes the form


The relevant boundary conditions are reduced to


where prime denotes derivative with respect toξ and the key crucial parameters are defined by:


HereM is the magnetic parameter,Nt is the thermophoresis parameter,Sc is the Schmidt number,Pr is Prandtl number,Ec is the Eckert number andNb is the Brownian motion parameter. Also, the physical quantities of interest skin friction coefficient, local Nusselt number and local Sherwood number are respectively defined as:


whereτw, qw and qm are wall shear stress, local heat flux and local mass flux at the stretching surface serially given as:


3. Results and discussion

Present study finds numerical solution of differential Eqs. (8)(10) subjected to the boundary conditions (11) and (12) that are computed using RKF method by applying shooting technique. The main reason behind to solve the present problem are to determine the impact of prominent fluid parameters namely Eckert numberEc, thermophoresisNt, Brownian motion parameterNb, Schmidt numberSc and magnetic parameterM onfδ0,θδ0 andϕδ0. Table 1 demonstrate the impact of fluid parametersNb,Nt,Ec,Sc and M on skin friction coefficientfδ0, local Nusselt numberθδ0 and local Sherwood numberϕδ0 by taking fixed entries of fluid parameters Prandtl number, Pr as 5.0 and surface temperature parameter m as 1.0.


Table 1.

Values of skin friction coefficient fδ0, local Nusselt number θδ0 and local Sherwood number ϕδ0 for crucial fluid parameters Nb,Nt,Ec,Sc and M with fixed entries of Pr=5.0 and m=1.0.

Figure 2 manifests variation in fluid velocity against magnetic parameter M0.81.01.2. This figure shows that existence of magnetic parameter M resists the fluid particle to move freely and main reason behind the resistance is that magnetic parameter M produces Lorentz force and this magnetism behavior can be adopted for controlling the fluid movement. Thus, enhancement in the value of magnetic parameter M causes the declination of velocity distribution

Figure 2.

Impact of magnetic parameter M on velocity profile fδξ.

Figure 3 examines temperature distribution variation against the fluid parameter Brownian motion parameter Nb0.050.100.15. The striking of atoms or molecules of the fluid particles with each other will create an arbitrary motion called Brownian motion of suspended (pendulous) particles and that will enhances width of boundary layer. Hence, fluid temperature increases for higher Brownian motion parameter Nb and in consequence local Nusselt number decreases.

Figure 3.

Impact of Brownian motion parameter Nb on temperature profile θδξ.

Figure 4 deliberates the impact of fluid temperature under the consequence of thermophoresis parameter Nt0.20.30.4. Temperature gradient falls down for higher values of thermophoresis parameter Nt that result in reduction of conduction of nanoparticles. Thus, width of boundary layer enhances due to reallocation of ultrafine particles from hotter to colder part and hence, temperature enhances for higher thermophoresis parameter Nt that can be seen in Figure 4.

Figure 4.

Impact of thermophoresis parameter Nt on temperature profile θδξ.

Figure 5 demonstrate fluid temperature variation against Eckert number Ec0.00.10.2. A dimensionless quantity Ec is the fraction of advective transportation and heat dissipation potential. As Eckert number Ec enhances, thermal buoyancy effect raises that results in increasing temperature and that is the main reason behind the conversion of kinetic energy into thermal energy. Hence, fluid temperature enhances because of this conversion effect. Consequently, declination in Nusselt number Nux is noticed that can be seen via Table 1.

Figure 5.

Impact of Eckert number Ec on temperature profile θδξ.

Figure 6 reflects variation of temperature distribution against magnetic parameter M. With an increases in magnetic parameter M, velocity profile decreases because of generation of Lorentz force that consequently intensify the boundary thickness and rate of heat transportation and hence fluid temperature enhances as shown via Figure 6.

Figure 6.

Impact of magnetic parameter M on temperature profile θδξ.

Figure 7 manifests the impact of Brownian motion parameter Nb0.050.100.15 on nanoparticle concentration ϕδξ. With an increase in the value of Brownian motion parameter Nb, fluid particles collides with each other with higher speed which results in increase in the nanoparticle concentration and consequently, local Sherwood number reduces as depicted in the Table 1.

Figure 7.

Impact of Brownian motion parameter Nb on Concentration profile ϕδξ.

Figure 8 portraits variation for nanoparticle volume fraction ϕδξ against thermophoresis parameter Nt0.20.30.4. This graph shows that with an increase in thermophoresis parameter, nanoparticle concentration increases. Basically, in case of thermophoresis force applied by a particle on the other particle will generates the movement of particles from hotter to colder part and hence fluid moves from hotter to colder region and hence intensification in the nanoparticle volume fraction is observed via Figure 8.

Figure 8.

Impact of thermophoresis parameter Nt on Concentration profile ϕδξ.

Figure 9 portraits the impact of Schmidt number Sc1.11.41.7 on profile of nanoparticle concentration. Intensification in the value of physical parameter Sc, declination in mass diffusivity is observed. Due to this effect nanoparticle concentration decreases.

Figure 9.

Impact of Schmidt number Sc on Concentration profile ϕδξ.

Figure 10 reflects the variation for nanoparticle concentration ϕδξ against the magnetic parameter M0.81.01.2. With increase in magnetic parameterM, rate of mass transportation decreases that consequently increase nanoparticle concentration and hence reduction in the value of local Sherwood number is notice as seen in Table 1.

Figure 10.

Impact of magnetic parameter M on Concentration profile ϕδξ.


4. Conclusions

Present study reflects the heat, mass and flow transportation of Magnetohydrodynamic (MHD) nanofluid towards a sheet which is stretched linearly. Key findings of current analysis are summarized as:

  1. Skin friction coefficient elevates with increment in magnetic parameter M due to produced Lorentz force that ultimately improves local Sherwood number along with Nusselt number for higher magnetic parameter M.

  2. Fluid temperature enhances for greater values of physical parameters Eckert number Ec, Brownian motion parameter Nb and thermophoresis parameter Nt.

  3. An enhancement in the profile of nanoparticle concentration is noticed for greater values of thermophoresis parameter Nt,. Whereas, it declines for Brownian motion parameter Nb and Schmidt number Sc.




Cartesian coordinates


Magnetic field intensity


Positive constant


Surface temperature parameter


Schmidt number




Nanoparticle volume fraction


Mass flux


Prandtl number


Nusselt number


Thermophoresis parameter


Ambient temperature attained


Sheerwood number


Vertical velocity


Brownian diffusion coefficient


Brownian motion parameter




Stretching velocity


Thermophoresis diffusion coefficient


Temperature at the sheet


Hass flux


Ambient nanoparticle volume fraction


Horizontal velocity


Magnetic parameter

Greek symbolsν

Kinematic viscosity


Casson fluid parameter


Electrical conductivity


Similarity variable


Thermal diffusivity


Ratio of heat capacities


Non-dimensional temperature


Non-dimensional nanoparticle concentration


  1. 1. L. J. Crane, “Flow past a stretching plate,” Zeitschrift fr angewandte Mathematik und Physik ZAMP, vol. 21, pp. 645–647, July 1970.
  2. 2. P. S. Gupta and A. S. Gupta, “Heat and mass transfer on a stretching sheet with suction or blowing,” The Canadian Journal of Chemical Engineering, vol. 55, pp. 744–746, Dec. 1977.
  3. 3. B. B. Yoon, R. S. Rao, and N. Kikuchi, “Sheet stretching: A theoretical-experimental comparison,” International Journal of Mechanical Sciences, vol. 31, pp. 579–590, Jan. 1989.
  4. 4. M. S. Sarma and B. N. Rao, “Heat transfer in a viscoelastic fluid over a stretching sheet,” Journal of Mathematical Analysis and Applications, vol. 222, pp. 268–275, 1998.
  5. 5. K. Vajravelu, “Viscous flow over a nonlinearly stretching sheet,” Appl. Math. Comput., vol. 124, pp. 281–288, Dec. 2001.
  6. 6. R. Cortell, “Viscous flow and heat transfer over a nonlinearly stretching sheet,” Applied Mathematics and Computation, vol. 184, pp. 864–873, Jan. 2007.
  7. 7. S. Nadeem, A. Hussain, and K. Vajravelu, “Effects of heat transfer on the stagnation flow of a third order fluid over a shrinking sheet,” Zeitschrift fr Naturforschung A, vol. 65, pp. 969–994, Nov. 2010.
  8. 8. K. Prasad, K. Vajravelu, and P. Datti, “Mixed convection heat transfer over a non-linear stretching surface with variable fluid properties,” International Journal of Non-Linear Mechanics, vol. 45, pp. 320– 330, Apr. 2010.
  9. 9. S. U. S. Choi and J. A. Eastman, “Enhancing thermal conductivity of fluids with nanoparticles,” ASME International Mechanical Engineering Congress and Exposition, pp. 1-9, Nov. 1995.
  10. 10. J. Buongiorno, “Convective transport in nanofluids,” Journal of Heat Transfer, vol. 128, pp. 240–250, March 2006.
  11. 11. A. Kuznetsov and D. Nield, “Natural convective boundary-layer flow of a nanofluid past a vertical plate,” International Journal of Thermal Sciences, vol. 49, pp. 243–247, Feb. 2010.
  12. 12. C. Mamaloukas, S. Spartalis, and H. P. Mazumdar, “MHD flow of a Newtonian fluid over a stretching sheet: an approximate solution,” International journal of Computational and Numerical Analysis and Applications, vol. 1, pp. 299–310, 2002.
  13. 13. R. Tamizharasi and V. Kumaran, “Pressure in MHD/Brinkman flow past a stretching sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, pp. 4671–4681, Dec. 2011.
  14. 14. T. Hayat, M. Qasim, and S. Mesloub, “MHD flow and heat transfer over permeable stretching sheet with slip conditions,” International Journal for Numerical Methods in Fluids, vol. 66, pp. 963–975, July 2011.
  15. 15. A. Hussain, M. Malik, T. Salahuddin, S. Bilal, and M. Awais, “Combined effects of viscous dissipation and joule heating on MHD sisko nanofluid over a stretching cylinder,” Journal of Molecular Liquids, vol. 231, pp. 341–352, Apr. 2017.
  16. 16. K. Vajravelu and J. Cannon, “Fluid flow over a nonlinearly stretching sheet,” Applied Mathematics and Computation, vol. 181, pp. 609–618, Oct. 2006.
  17. 17. M. H. Matin, M. R. H. Nobari, and P. Jahangiri, “Entropy Analysis in Mixed Convection MHD flow of Nanofluid over a Non-linear Stretching Sheet,” Journal of Thermal Science and Technology, vol. 7, no. 1, pp. 104–119, 2012.
  18. 18. H. M. Shawky, N. T. M. Eldabe, K. A. Kamel, and E. A. Abd-Aziz, “MHD flow with heat and mass transfer of Williamson nanofluid over stretching sheet through porous medium,” Microsystem Technologies, vol. 25, pp. 1155–1169, Apr. 2019.
  19. 19. M. F. Md Basir, M. J. Uddin, A. I. Md. Ismail, and O. A. Bg, “Nanofluid slip flow over a stretching cylinder with schmidt and pclet number effects,” AIP Advances, vol. 6, pp. 1-15, May 2016.
  20. 20. A. K. Pandey and M. Kumar, “Boundary layer flow and heat transfer analysis on cu-water nanofluid flow over a stretching cylinder with slip,” vol. 56, no. 4, pp. 671–677, Dec. 2017.
  21. 21. V. Vinita and V. Poply, “Impact of outer velocity MHD slip flow and heat transfer of nanofluid past a stretching cylinder,” Materials Today: Proceedings, Vol. 26, pp. 3429-3435, 2020.
  22. 22. Vinita, Poply V, Goyal R, Sharma N., “Analysis of the velocity, thermal, and concentration MHD slip flow over a nonlinear stretching cylinder in the presence of outer velocity,” Heat Transfer, 2020; 1–27.
  23. 23. Goyal R, Vinita, Sharma N, Bhargava R., “GFEM analysis of MHD nanofluid flow toward a power-law stretching sheet in the presence of thermodiffusive effect along with regression investigation,” Heat Transfer, 2020;1–23.
  24. 24. Poply V, Singh P, Choudhary KK, “Analysis of laminar boundary layer flow along a stretching cylinder in the presence of thermal radiation,” WSEAS Trans Fluid Mech., vol. 8, pp. 159-164, 2013.
  25. 25. P. Rana and R. Bhargava, “Flow and heat transfer of a nanofluid over a nonlinearly stretching sheet: A numerical study,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, pp. 212– 226, Jan. 2012.

Written By

Vikas Poply

Submitted: 26 July 2020 Reviewed: 06 December 2020 Published: 28 December 2020