Open access peer-reviewed chapter

Thermal Conductivity in the Boundary Layer of Non-Newtonian Fluid with Particle Suspension

Written By

Rudraswamy N.G., Ganeshkumar K., Krishnamurthy M.R., Gireesha B.J. and Venkatesh P.

Submitted: 27 November 2017 Reviewed: 09 March 2018 Published: 05 September 2018

DOI: 10.5772/intechopen.76345

From the Edited Volume

Impact of Thermal Conductivity on Energy Technologies

Edited by Aamir Shahzad

Chapter metrics overview

1,354 Chapter Downloads

View Full Metrics


The present chapter is focused on studies concerned with three-dimensional flow and heat transfer analysis of Carreau fluid with nanoparticle suspension. The heat transfer analysis in the boundary was carried out with the fluid flow over a stretching surface under the influence of nonlinear thermal radiation, mixed convection and convective boundary condition. Suitable similarity transformations are employed to reduce the governing partial differential equations into coupled nonlinear ordinary differential equations. The equations in non-linear form are then solved numerically using Runge-Kutta-Fehlberg fourth fifth-order method with the help of symbolic algebraic software MAPLE. The results so extracted are well tabulated and adequate discussions on the parameters affecting flow and heat transfer analysis were carried out with the help of plotted graphs.


  • Carreau nano fluid
  • nonlinear thermal radiation
  • mixed convection
  • stretching sheet
  • convective boundary condition
  • numerical method

1. Introduction

Thermal radiation, the fundamental mechanism of heat transfer is an indispensable activity in rocket propulsion, plume dynamics, solar collector performance, materials processing, combustion systems, fire propagation and other industrial and technological processes at high temperatures. With the developments in computational dynamics, increasing attention has been diverted towards thermal convection flows with the significant radiative flux. Rayleigh initiated the theory of thermal convection, by deriving critical temperature gradient (Critical Rayleigh number). Importance of such radiations is intensified with absolute temperatures at higher level. Thus a substantial interest is driven towards thermal boundary layer flows with a strong radiation. Governing equation of radiative heat transfer with its integro-differential nature makes numerical solutions of coupled radiative-convective flows even more challenging. Multiple studies were conducted employing several models to investigate heat and mass transfer in boundary layer and fully-developed laminar convection flows. As a consequence several simultaneous multi-physical effects in addition to radiative heat transfer including gravity and pressure gradient effects [1], mhd flow of nanofluids [2], buoyancy effects [3, 4], ferrofluid dynamics [5], stretching surface flow [6, 7], time-dependent, wall injection and Soret/Dufour effects [8, 9, 10, 11].

These studies have however been confined to Newtonian flows. But industries related with fabrication of polymers and plastics at high temperatures show greater importance towards radiative flows of non-Newtonian fluids. The potential of non-Newtonian flows in ducts with radiative transfer were significantly developed after the studies on novel propellants for spacecraft [12]. The developments are extant and diversified the application of non-Newtonian fluid models. Most studies in this regard have employed the Rosseland model which is generally valid for optically-thick boundary layers. Recently, Kumar et al. [13] used such model to study melting heat transfer of hyperbolic tangent fluid over a stretching sheet with suspended dust particles. Cortell [14] and Batalle [15] have shown their earlier contribution towards radiative heat transfer of non-Newtonian fluids past stretching sheet under various circumstances. Relating to the studies Khan et al. [16] developed a numerical studies correlating MHD flow of Carreau fluid over a convectively heated surface with non-linear radiation. Appending to this studies Khan et al. [17] provided his results on hydromagnetic nonlinear thermally radiative nanoliquid flow with Newtonian heat along with mass conditions. Meanwhile, Rana and Bhargava [18] provided a numerical elucidation to study of heat transfer enhancement in mixed convection flow along a vertical plate with heat source/sink utilizing nanofluids. Hayat et al. [19] investigated the mixed convection stagnation-point flow of an incompressible non-Newtonian fluid over a stretching sheet under convective boundary conditions. Many diverse -physical simulations with and without convective and/or radiative heat transfer have been studied. Representative studies in this regard include [20, 21, 22, 23] with analogous to the property of radiation flow.

Endeavoring the complications in three dimensional flow analysis, Shehzad et al. [24] studied the effect of thermal radiation in Jeffrey nanofluid by considering the characteristics of thermophoresis and Brownian motion for a solar energy model. Hayat et al. [25] analyzed the effect non-linear thermal radiation over MHD three-dimensional flow of couple stress nanofluid in the presence of thermophoresis and Brownian motion. Rudraswamy et al. [26] observations on Soret and Dufour effects in three-dimensional flow of Jeffery nanofluid in the presence of nonlinear thermal radiation clearly showed that concentration and associated boundary layer thickness are enhanced by increasing Soret and Dufour numbers. Many such problems [27, 28, 29] were considered disclosing the feature of thermal radiation in three dimensional flow of non-Newtonian fluids.

Inspired by the above works, we put forth the studies on the effect of non-linear thermal radiation on three dimensional flow of Carreau fluid with suspended nanoparticles. Present studies even include the phenomenon of mixed convection and convective boundary conditions. A numerical approach is provided for the above flow problem by employing Runge-Kutta-fourth-fifth order method.


2. Mathematical formulation

A steady three-dimensional flow of an incompressible Carreau fluid with suspended nano particles induced by bidirectional stretching surface at z=0 has been considered. The sheet is aligned with the xy plane z=0 and the flow takes place in the domain z>0. Let u=uwx=ax and v=vwy=by be the velocities of the stretching sheet along x and y directions respectively. A constant magnetic field of strength B is applied in the z direction. Heat and mass transfer characteristics are taken in to account in the presence of Brownian motion and Thermophoresis effect. The thermo physical properties of fluid are taken to be constant.

Extra stress tensor for Carreau fluid is.


In which τ¯ij is the extra stress tensor, μ0 is the zero shear rate viscosity, Γ is the time constant, γ̇¯ is the power law index and is defined as.


Here Π is the second invariant strain tensor.

The governing boundary layer equations of momentum, energy and concentration for three-dimensional flow of Carreau nanofluid can be written as,


The boundary conditions for the present flow analysis are,


where ν is the kinematic viscosity of the fluid, μ is the coefficient of fluid viscosity, ρ is the fluid density, B is the magnetic field, σ is the electrical conductivity of the fluid, T is the fluid temperature, α is the thermal diffusivity of the fluid, k is the thermal conductivity. τ is the ratio of effective heat capacity of the nanoparticle material to heat capacity of the fluid, qr is the radiative heat flux, g is the gravitational acceleration, βT is thermal expansion coefficient of temperature, DB is the Brownian diffusion coefficient, hf is the heat transfer coefficient, DT is the thermophoretic diffusion coefficient, cp is the specific heat at constant pressure, Tf is the temperature at the wall, T is the temperatures far away from the surface. C is the concentration and C is the concentration far away from the surface. The subscript w denotes the wall condition.

Using the Rosseland approximation radiation heat flux qr is simplified as,


where σ and k are the Stefan-Boltzmann constant and the mean absorption coefficient respectively.

In view to Eq. (8), Eq. (4) reduces to.


The momentum, energy and concentration equations can be transformed into the corresponding ordinary differential equations by the following similarity variables,


where T=T1+θw1θη, θw=TfT, θw>1 being the temperature ratio parameter.

Then, we can see that Eq. (1) is automatically satisfied, and Eqs. (2)(7)are reduced to:


With the boundary conditions,


We=cUw2λ2ν is the Weissenberg number, M=σB2ρa is the magnetic parameter, c=ba is the ratio of stretching rates, Pr=να is Prandtl number, R=16σT33kk is the radiation parameter, Nb=τDBCwCν is the Brownian motion parameter, Nt=τDTTfTνT is the Thermophoresis parameter, λ=gβTTfTRex is the mixed convection parameter, Bi=hfkνa is the Biot number, Le=αDB is the Lewis number.

The local skin friction Cf, local Nusselt number Nux and local number Sherwood Shx are defined as,


The local skin friction, local Nusselt number and Sherwood number is given by,


where Rex=uwxν is the local Reynolds number based on the stretching velocity. uwx.


3. Numerical method

The non-linear ordinary differential Eqs. (11)(14) subjected to boundary conditions (15) has been solved using the Runge-Kutta-Fehlberg fourth-fifth order method with the help of symbolic algebraic software MAPLE. The boundary conditions for η= are replaced by fηmax=1,θηmax=0 and ϕηmax=0, where ηmax is a sufficiently large value of η at which the boundary conditions (15) are satisfied. Thus, the values of η=ηmax are taken to be 6. To validate the employed method, the authors have compared the results of f0 and g0 with the that of published works by Wang [27] and Hayat [30] for the different values stretching parameter. These comparisons are given in Table 1 and it shows that the results are in very good agreement.

cWang [27]Hayat et al. [30]Present result

Table 1.

Comparison of different values of c with Wang [27] and Hayat et al. [30].


4. Result and discussion

The purpose of this section is to analyze the effects of various physical parameters on the velocities, temperature and concentration fields. Therefore, for such objective, Figures 111 has been plotted. Observations over these data with plotted graphs are discussed below.

Figure 1 characterizes the influence of Weissenberg number We on velocity profiles for both x and y direction. It is found that increasing values of the Weissenberg number increases the momentum boundary layers in both directions. Physically, Weissenberg number is directly proportional to the time constant and reciprocally proportional to the body. The time constant to body magnitude relation is higher for larger values of Weissenberg number. Hence, higher Weissenberg number causes to enhance the momentum boundary layer thickness.

Figure 1.

Influence of We on velocity profiles of both fη and gη.

The developments of a magnetic field M on velocity profiles are circulated in Figure 2. We tend to discover depreciation within the velocity profile for ascent values of magnetic field parameter. Physically, the drag force will increase with a rise within the magnetic flux and as a result, depreciation happens within the velocity field.

Figure 2.

Influence of M on velocity profiles fη and gη.

Figure 3 designed the velocity profiles of f and g for various values of stretching parameter c. The velocity profiles and associated momentum boundary layer thickness decrease, once the stretching parameter will increase whereas velocity profile g’, exhibits the opposite behavior of f’. Figure 4 shows the velocity profiles for different values of mixed convection parameter λ. It depicts that the velocity field and momentum boundary layer thickness increases in both x and y direction by increasing mixed convection parameter.

Figure 3.

Influence of c on velocity profiles fη and gη.

Figure 4.

Influence of λ on velocity profiles fη and gη.

Figure 5 portraits the consequences of Brownian motion parameter on temperature and concentration profile. The Brownian motion parameter Nb will increase the random motion of the fluid particles and thermal boundary layer thickness conjointly will increase which ends up in an additional heat to provide. Therefore, temperature profile will increase however concentration profiles show opposite behavior.

Figure 5.

Influence of Nb on θη and ϕη profiles.

The development of the thermophoresis parameter Nt on temperature and concentration profiles is inspecting in Figure 6. Form this figure we observed that, the higher values of thermophoresis parameter is to increases both θη and ϕη profiles. Further, the thermal boundary layer thickness is higher for larger values of thermophoresis parameter. This is because, it’s a mechanism within which little particles area unit force off from the new surface to a chilly one. As a result, it maximizes the temperature and concentration of the fluid.

Figure 6.

Influence of Nt on θη and ϕη profiles.

Figure 7 describe the influences of Biot number Bi on temperature profile. One can observe form the figure, the larger values of Biot number cause an enhancing the temperature profile. This is because, the stronger convection leads to the maximum surface temperatures which appreciably enhance the thermal boundary layer thickness.

Figure 7.

Influence of Bi on temperature profile.

Figures 8 and 9 are sketched to analyze the effect of radiation parameter R and temperature ratio θw parameter on temperature profile. The above graphs elucidate that, the temperature profile and thermal boundary layer thickness area unit increased by ascent values of radiation parameter and temperature ratio. Larger values of thermal radiation parameter provide more heat to working fluid that shows an enhancement in the temperature and thermal boundary layer thickness.

Figure 8.

Influence of R on temperature profile.

Figure 9.

Influence of θw on temperature profile.

The effect of the Prandtl number Pr on θη is seen in Figure 10. Since Pr is that the magnitude relation of the viscous diffusion rate to the thermal diffusion rate, the upper worth of Prandtl number causes to scale back the thermal diffusivity. Consequently, for increasing values of Pr, the temperature profile gets decreases. The impact of Lewis number Le on nanoparticle concentration is plotted in Figure 11. It is evident that the larger values of Lewis number cause a reduction in nanoparticles concentration distribution. Lewis number depends on the Brownian diffusion coefficient. Higher Lewis number leads to the lower Brownian diffusion coefficient, which shows a weaker nanoparticle concentration.

Figure 10.

Influence of Pr on temperature profile.

Figure 11.

Influence of Le on concentration profile.

Table 2 presents the numerical values of skin friction for various physical values in the presence and absence We=n=0 of non-Newtonian fluid. It is observed that skin friction increase in both directions with increasing c for both presence and absence of non-Newtonian fluid. In the other hand, the skin friction coefficient decreases in both directions by increasing Bi. The skin friction is higher in the presence of non-Newtonian fluid than in the absence of non-Newtonian fluid.


Table 2.

Numerical result of skin friction coefficient for different physical parameter values for present and absence non Newtonian fluid.

Table 3 also elucidates that, the wall temperature for different physical parameter for linear as well as nonlinear radiation. It reveals that, the wall temperature increases for increasing values of Bi,R and c for both linear and nonlinear radiation but the wall temperature decreases by increasing Le,Nb,Nt and Pr. Further, it is noticed that the wall temperature is higher for nonlinear radiation than that linear radiation.


Table 3.

Numerical result of Nusselt number for different physical parameter values for linear and non nonlinear radiation.

Table 4 clearly shows the numerical values of skin friction, Nusselt number and Sherwood number for various physical parameters values. It reveals that, numerical values of wall temperature θ0 increase by increasing Bi,θw,R and c. In the other hand Nusselt number decreases by increasing. Le,M,Nb,Nt and Pr. From this table, the skin friction coefficient increases by increasing Bi and m. Further, the Sherwood number increases by increasing Bi,θw,R,Pr and We.


Table 4.

Numerical result of local skin friction coefficient, Sherwood number and Nusselt number for different physical parameter.


5. Conclusions

In the present study, influence of nonlinear radiation on three dimensional flow of an incompressible non-Newtonian Carreau nanofluid has been obtained. The obtained results are presented in tabulated and graphical form with relevant discussion and the Major findings from this study are:

The velocity profiles increase in x directions and decrease in the y direction by increasing the stretching parameter.

Concentration profile increase by increasing the values Nb but in case of Nt concentration profile decreases.

Nb and Nt parameter shows the increasing behavior for temperature profile.

Effects of Le nanoparticle fraction ϕη show the decreasing behavior.

Magnetic parameter reduces the velocity profiles in both x and y directions.

Temperature and thermal boundary layer thickness are decreased when the Pr and tl number increases.

Nonlinear thermal radiation should be kept low to use it as a coolant factor.

The rate of heat transfer increases with the increases in parameters Rd and θw.

We also noticed that the velocity profile and its associated boundary layer thickness are increases by increasing the values of We.


  1. 1. Anwar Bég O, Ghosh SK, Narahari M, Bég TA. Mathematical modelling of thermal radiation effects on transient gravity-driven optically-thick gray convection flow along an inclined plate with pressure gradient. Chemical Engineering Communications. 2011;198:1-15
  2. 2. Rahman MM, Eltayeb IA. Radiative heat transfer in a hydromagnetic nanofluid past a non-linear stretching surface with convective boundary condition. Meccanica. 2013;48:601-615
  3. 3. Rashidi MM, Ganesh NV, Abdul Hakeem AK, Ganga B. Buoyancy effect on MHD flow of nanofluid over a stretching sheet in the presence of thermal radiation. Journal of Molecular Liquids. October 2014;198:234-238
  4. 4. Srinivas S, Muthuraj R. Effects of thermal radiation and space porosity on MHD mixed convection flow in a vertical channel using homotopy analysis method. Communications in Nonlinear Science and Numerical Simulation. 2010;15:2098-2108
  5. 5. Zeeshan A, Majeed A, Ellahi R. Effect of magnetic dipole on viscous ferro-fluid past a stretching surface with thermal radiation. Journal of Molecular Liquids. 2016;215:549-554
  6. 6. Cortell R. Fluid flow and radiative nonlinear heat transfer over stretching sheet. Journal of King Saud University Science. 2013;26:161-167
  7. 7. Ferdows M, Khan MS, Anwar Bég O, Azad MAK, Alam MM. Numerical study of transient magnetohydrodynamic radiative free convection nanofluid flow from a stretching permeable surface. Proceedings of Institution of Mechanical E-Part E: Journal of Process Mechanical Engineering. 2014;228(3):181-196
  8. 8. Anwar Bég O, Ferdows M, Bég A, Tasveer T, Ahmed MW, Alam MM. Radiative optically-dense magnetized transient reactive transport phenomena with cross diffusion and dissipation effects. Numer. Simul. Journal of Taiwan Institute of Chemical Engineers. 2016:15
  9. 9. Pal D, Mondal H. Effects of Soret Dufour, chemical reaction and thermal radiation on MHD non-Darcy unsteady mixed convective heat and mass transfer over a stretching sheet. Communications in Nonlinear Science and Numerical Simulation. 2011;16(4):1942-1958
  10. 10. Murthy PVSN, Mukherjee S, Srinivasacharya D. Combined radiation and mixed convection from a vertical wall with suction/injection in a non–Darcy porous medium. Acta Mechanica. 2004;168(3–4):145
  11. 11. Sheikholeslami M, Ganji DD, Javed MY, Ellahi R. Effect of thermal radiation on magnetohydrodynamics nanofluid flow and heat transfer by means of two phase model. Journal of Magnetism and Magnetic Materials. 2015;374:36-43
  12. 12. S. Snyder, N. Arockiam, P.E. Sojka, Secondary Atomization of Elastic Non-Newtonian Liquid Drops, In: AIAA Joint Propulsion Conference, Nashville, Tennessee, USA, 2010
  13. 13. Kumar KG, Gireesha BJ, Rudraswamy NG, Gorla RSR. Melting heat transfer of hyperbolic tangent fluid over a stretching sheet with fluid particle suspension and thermal radiation. Communications in Numerical Analysis. 2017;2017(2):125-140
  14. 14. Cortell R. Effects of viscous dissipation and radiation on the thermal boundary layer over a nonlinearly stretching sheet. Physical Letters A. 2008;372(5):631-636
  15. 15. Batalle RC. Viscoelastic fluid flow and heat transfer over a stretching sheet under the effects of a non-uniform heat source, viscous dissipation and thermal radiation. International Journal of Heat and Mass Transfer. 2007;50(15–16):3152-3162
  16. 16. Khan M, Hussain M, Azam M. Magnetohydrodynamic flow of Carreau fluid over a convectively heated surface in the presence of non-linear radiation. JMMM. 2016;412:63-68
  17. 17. Khan MI, Alsaedi A, Shehzad SA, Hayat T. Hydromagnetic nonlinear thermally radiative nanoliquid flow with Newtonian heat and mass conditions. Results Physics. 2017;7:2255-2260
  18. 18. Rana P, Bhargava R. Numerical study of heat transfer enhancement in mixed convection flow along a vertical plate with heat source/sink utilizing nanofluids. Communication in Nonlinear Science Numerical Simulation. 2011;16(11):4318-4334
  19. 19. Hayat T, Shehzad SA, Alsaedi A, Alhothuali MS. Mixed convection stagnation point flow of Casson fluid with convective boundary conditions. Chinese Physics Letters. 2012;29(11):114704.
  20. 20. Laxmi T, Shankar B. Effect of nonlinear thermal radiation on boundary layer flow of viscous fluid over nonlinear stretching sheet with injection/suction. Journal of Applied Mathematical Physics. 2016;4:307-319
  21. 21. Kumar KG, Rudraswamy NG, Gireesha BJ, Manjunatha S. Non-linear thermal radiation effect on Williamson fluid with particle-liquid suspension past a stretching surface. Results in Physics. 2017;7:3196-3202
  22. 22. Kumar KG, Gireesha BJ, Manjunatha S, Rudraswamy NG. Effect of nonlinear thermal radiation on double-diffusive mixed convection boundary layer flow of viscoelastic nanofluid over a stretching sheet. International Journal of Mechanical Material Engineering. 2017;12(1):18
  23. 23. Mabood F, Imtiaz M, Alsaedi A, Hayat T. Unsteady convective boundary layer flow of Maxwell fluid with nonlinear thermal radiation: A numerical study. International Journal of Nonlinear Science Numerical Simulation. 2016;17(5):221-229
  24. 24. Shehzad SA, Hayat T, Alsaedi A, Obid MA. Nonlinear thermal radiation in three-dimensional flow of Jeffrey nanofluid: A model for solar energy. Applied Mathematics and Computation. 2016;248:273-286
  25. 25. Hayat T, Muhammad T, Alsaedi A, Alhuthali MS. Magnetohydrodynamic three-dimensional flow of viscoelastic nanofluid in the presence of nonlinear thermal radiation. JMMM. 2015;385:222-229
  26. 26. Rudraswamy NG, Kumar KG, Gireesha BJ, Gorla RSR. Soret and Dufour effects in three-dimensional flow of Jeffery nanofluid in the presence of nonlinear thermal radiation. Journal of Nanoengineering Nanomanufacturing. 2016;6(4):278-287
  27. 27. Wang CY. The three-dimensional flow due to a stretching sheet. Physics of Fluids. 1984;27:1915-1917
  28. 28. Hayat T, Muhammad T, Shehzad SA, Alsaedi A. Similarity solution to three dimensional boundary layer flow of second grade nanofluid past a stretching surface with thermal radiation and heat source/sink. AIP Advances. 2015;5:017107.
  29. 29. Shehzad SA, Hayat T, Alsaedi A. MHD three dimensional flow of viscoelastic fluid with thermal radiation and variable thermal conductivity. A Journal of Central South University. 2014;21(10):3911-3917
  30. 30. Hayat T, Shehzad SA, Alsaedi A. Three-dimensional stretched flow of Jeffrey fluid with variable thermal conductivity and thermal radiation. Applied Mathematics and Mechanics. 2013;34:823-832

Written By

Rudraswamy N.G., Ganeshkumar K., Krishnamurthy M.R., Gireesha B.J. and Venkatesh P.

Submitted: 27 November 2017 Reviewed: 09 March 2018 Published: 05 September 2018