Open access peer-reviewed chapter - ONLINE FIRST

Heat Transfer in a MHD Nanofluid Over a Stretching Sheet

By Vikas Poply

Submitted: July 26th 2020Reviewed: December 6th 2020Published: December 28th 2020

DOI: 10.5772/intechopen.95387

Downloaded: 43


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-sized109107particles 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, whereaand xdenotes a constant and stretching surface coordinate respectively. Tw=T+T0xmaty=0, whereT0refers to the positive constant,Trefers 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 byuandv, respectively. Alsoνdenotes kinematic viscosity,ρis the density of fluid, τ=ρcpρcfdefines a proportion of heat capacities,DTreflects thermophoretic diffusion coefficient,Bis the magnetic field intensity,DBdenotes 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:


HereMis the magnetic parameter,Ntis the thermophoresis parameter,Scis the Schmidt number,Pris Prandtl number,Ecis the Eckert number andNbis 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, qwand qmare 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 numberScand magnetic parameterMonfδ0,θδ0andϕδ0. Table 1 demonstrate the impact of fluid parametersNb,Nt,Ec,Scand Mon skin friction coefficientfδ0, local Nusselt numberθδ0and local Sherwood numberϕδ0by taking fixed entries of fluid parameters Prandtl number, Pras 5.0 and surface temperature parameter mas 1.0.


Table 1.

Values of skin friction coefficient fδ0, local Nusselt number θδ0and local Sherwood number ϕδ0for crucial fluid parameters Nb,Nt,Ec,Scand Mwith fixed entries of Pr=5.0and 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 Mresists the fluid particle to move freely and main reason behind the resistance is that magnetic parameter Mproduces Lorentz force and this magnetism behavior can be adopted for controlling the fluid movement. Thus, enhancement in the value of magnetic parameter Mcauses 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 Nband 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 Ntthat 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 Ntthat 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 Ecis the fraction of advective transportation and heat dissipation potential. As Eckert number Ecenhances, 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 Nuxis 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.15on 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.7on 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 Mdue 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 Nband 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 Nband 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

Download for free

chapter PDF

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Vikas Poply (December 28th 2020). Heat Transfer in a MHD Nanofluid Over a Stretching Sheet [Online First], IntechOpen, DOI: 10.5772/intechopen.95387. Available from:

chapter statistics

43total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us