Research article
Open Access
Numerical Method for the Analysis of Thermal Radiation
on Heat Transfer in Nanofluid
J.I. Oahimire^{1*}, F.E Bazuaye^{2}, Taylor S. Harry^{2}
^{1}Department of Mathematics, Michael Okpara University of Agriculture, Umudike, Nigeria
^{2}Department of Mathematics and Statistics, University of Port Harcourt, Port Harcourt, Nigeria
^{2}Department of Mathematics and Statistics, University of Port Harcourt, Port Harcourt, Nigeria
*Corresponding author: Sergei Department of Mathematics, Michael Okpara University of Agriculture, Umudike, Nigeria, Email:
@
Received: July 02, 2016; Accepted: July 23, 2016; Published: July 30, 2016
Citation: Oahimire JI, Bazuaye FE, Harry TS (2016) Numerical Method for the Analysis of Thermal Radiation on Heat Transfer in
Nanofluid. Nanosci Technol 3(1): 14. DOI: http://dx.doi.org/10.15226/23748141/3/1/00135
AbstractTop
This present study investigates the effects of thermal radiation on
heat transfer in nanofluid. The governing equations are formulated
base on already existing model and are transformed to ordinary
differential equations using stream function and similarity variables.
The resulting dimension less equations are then solved numerically
by RungeKuta Fehlberg method with shooting techniques using
Maple software. With the help of graphs and tables, influences of
the governing parameters are discussed. The result shows that the
radiation has significance influences on heat transfer.
Keywords: Nanofluid; Runge KutaFehlberg method; Shooting technique; Heat transfer; Radiation
Keywords: Nanofluid; Runge KutaFehlberg method; Shooting technique; Heat transfer; Radiation
Nomenclature
a = Constant; g = Acceleration due to gravity; k = Thermal
Conductivity; P_{r} = Prandtl Number; T = Fluid Temperature;
T_{w} = Surface Temperature; T_{∞} = Free Stream Temperature;
u,v = Velocity Components; x, y = Cartesian Coordinates; f(x)
= Dimensionless Stream Function; G_{r} = Grashof Number; q_{r} =
Heat Flux Radiation; B_{o} = Magnetic Field of Constant Strength;
R = Radiation Parameter; K_{s} = Rosseland Mean Absorption
Coefficient; K = Thermal Conductivity Coefficient
Greek Symbols
β = Thermal Expansion Coefficient; μ = Dynamic Coefficient
of Viscosity; θ (η) = Dimensionless Temperature; η = Similarity
Variable; ρ = Fluid Density; ψ = Stream Function; σ' = Stefan
Bottzman Constant
Introduction
Nanofluid is formed when Nanoparticles such as Aluminium
Oxide (Al2O3), Copper (Cu), Copper Oxide (CuO), Gold (Au), Silver
(Ag), Silica particles, e.t.c are mixed with base fluids such as
water, oil, acetone, ethylene, etc. The discovery of nanofluids in
enhancing heat transfer in industrial processes has drawn the
attention of both scholars and industrialists to make researches
into this relatively new area.
Wang and Choi, [1] studied the thermal condutivity of nanoparticle fluid mixture containing Al_{2}O_{3} and CuO nanoparticles and showed that the thermal conductivity of nanofluids increased with increasing volume fraction of the nanoparticles. Steve (2006) found out that a small size radiator filled with well prepared nanofluid containing monosized nanoparticles with 2nm nominal diameter will be okay for heavy vehicle instead of a very big radiator and this will save cost and space. And he also discovered that fuel efficiency will be reduced by 5% in these Heavy Trucks when nanofluids are introduced into their cooling system. John and Baldev, [2] carried out experiment and they concluded that by controlling the linear aggregation length from nano to macron scales, the thermal conductivity of the nanofluid can be tuned from a low to very high value and that under repeated magnetic cycling, the thermal conductivity is reversible.
Khan and Ipop, [3] studied the boundarylayer flow of a nanofluid past a stretching sheet and their result shows that the heat transfer at fixed values of the Lewis number (Le), the Brownian motion parameter (Nb) and the thermophoresis parameter (Nt) increases with the Prandtl number (Pr). This is because a higher Pr fluid has relatively lower thermal conductivity which reduces conduction and thereby increases the heat transfer rate at the surface of the sheet.
Godson, et al. [4] in separate papers agreed that Nanofluids are the most option to enhance heat transfer and showed that a very small amount of nanoparticles suspended stably in base fluids can provide impressive improvement in the thermal property of such fluids.
Zeinali and Salehi, [5] studied experimentally, the effect of magnetic field of various strengths on the thermal performance of Silver/water nanofluid and their results showed that the thermal efficiency in the presence of magnetic field increased.
Yohannes and Shankar, [6] examine the boundary layer flow of heat and mass transfer in MHD flow of nanofluids through porous media with viscous dissipation and chemical reaction. The governing boundary layer equations were formulated and reduced to a set of ordinary differential equations using similarity transformations and then solved numerically by an explicit finite differential scheme known as the Keller box method. WanMohd, [7] applied Rungekutta Fehlberg method and shooting technique for solving classical Blasius equation [7]. The wellknown Blasius equation is governed by the third order nonlinear ordinary differential equation and then solved numerically using rungekutta Fehlberg method with shooting techniques. Hamad, [8] found analytical solution to the differential equations that model heat transfer in nanofliud. We extended the work of Hamad, [8] by incorporating radiation term to have differential equations that model the effect of thermal radiation on heat transfer of a nanofliud.
Wang and Choi, [1] studied the thermal condutivity of nanoparticle fluid mixture containing Al_{2}O_{3} and CuO nanoparticles and showed that the thermal conductivity of nanofluids increased with increasing volume fraction of the nanoparticles. Steve (2006) found out that a small size radiator filled with well prepared nanofluid containing monosized nanoparticles with 2nm nominal diameter will be okay for heavy vehicle instead of a very big radiator and this will save cost and space. And he also discovered that fuel efficiency will be reduced by 5% in these Heavy Trucks when nanofluids are introduced into their cooling system. John and Baldev, [2] carried out experiment and they concluded that by controlling the linear aggregation length from nano to macron scales, the thermal conductivity of the nanofluid can be tuned from a low to very high value and that under repeated magnetic cycling, the thermal conductivity is reversible.
Khan and Ipop, [3] studied the boundarylayer flow of a nanofluid past a stretching sheet and their result shows that the heat transfer at fixed values of the Lewis number (Le), the Brownian motion parameter (Nb) and the thermophoresis parameter (Nt) increases with the Prandtl number (Pr). This is because a higher Pr fluid has relatively lower thermal conductivity which reduces conduction and thereby increases the heat transfer rate at the surface of the sheet.
Godson, et al. [4] in separate papers agreed that Nanofluids are the most option to enhance heat transfer and showed that a very small amount of nanoparticles suspended stably in base fluids can provide impressive improvement in the thermal property of such fluids.
Zeinali and Salehi, [5] studied experimentally, the effect of magnetic field of various strengths on the thermal performance of Silver/water nanofluid and their results showed that the thermal efficiency in the presence of magnetic field increased.
Yohannes and Shankar, [6] examine the boundary layer flow of heat and mass transfer in MHD flow of nanofluids through porous media with viscous dissipation and chemical reaction. The governing boundary layer equations were formulated and reduced to a set of ordinary differential equations using similarity transformations and then solved numerically by an explicit finite differential scheme known as the Keller box method. WanMohd, [7] applied Rungekutta Fehlberg method and shooting technique for solving classical Blasius equation [7]. The wellknown Blasius equation is governed by the third order nonlinear ordinary differential equation and then solved numerically using rungekutta Fehlberg method with shooting techniques. Hamad, [8] found analytical solution to the differential equations that model heat transfer in nanofliud. We extended the work of Hamad, [8] by incorporating radiation term to have differential equations that model the effect of thermal radiation on heat transfer of a nanofliud.
Mathematical formulation
Consider the steady laminar twodimensional flow of an
incompressible viscous nanofluid past a linearly semiinfinite
stretching sheet under the influence of a constant magnetic field
of strength B0 which is applied normally to the sheet. x'and y' are
the coordinates along and perpendicular to the sheet while u' and
v' are the velocity components in the directions respectively as
shown in the figure below:
Figure 1:
A water based nanofluid containing different types of
nanoparticles: Al_{2}
O_{3}, Cu, TiO_{2} and Ag is used with the assumption
that both the fluid and the nanoparticles are in thermal
equilibrium. Based on the already existing model of Hamad,
[8], the governing boundary layer equations of continuity,
momentum and heat equations formulated are:
$$\frac{\partial u\text{'}}{\partial x\text{'}}+\frac{\partial v\text{'}}{\partial y\text{'}}=0\left(1\right)$$
$${\rho}_{nf}\left[u\text{'}\frac{\partial {u}^{\prime}}{\partial {x}^{\prime}}+v\text{'}\frac{\partial {v}^{\prime}}{\partial {y}^{\prime}}\right]={\mu}_{nf}\frac{{\partial}^{2}u\text{'}}{\partial y{\text{'}}^{2}}\sigma {B}_{0}{u}^{\prime}+g{\beta}_{t}\left({T}^{\prime}{{T}^{\prime}}_{\infty}\right)\left(2\right)$$
$${\left(\rho {c}_{p}\right)}_{nf}\left[u\text{'}\frac{\partial T\text{'}}{\partial {x}^{\prime}}+v\text{'}\frac{\partial T\text{'}}{\partial {y}^{\prime}}\right]={K}_{nf}\frac{{\partial}^{2}T\text{'}}{\partial y{\text{'}}^{2}}\frac{\partial {q}_{r}}{\partial y\text{'}}\left(3\right)$$
The boundary conditions of the equations are
$${u}^{\prime}={u}_{w}^{\text{'}}\left({x}^{\prime}\right)=a{x}^{\prime},{v}^{\prime}=0,{T}^{\text{'}}={T}_{w}^{\text{'}}at{y}^{\prime}=0$$
$${u}^{\prime}\to 0,{T}^{\prime}\to {T}_{\infty}^{\text{'}},y\text{'}\to \infty \left(4\right)$$
Where qr is the radiative heat flux, T′ is the temperature of the
fluid, x' and y' are the coordinates along and perpendicular to the
sheet while u' and v' are the velocity components in the x' and y'
directions respectively and a is a constant.
The effective density (ρnf), effective dynamic viscosity (μnf),
heat capacitance (ρCp)nf and the effective thermal conductivity
(knf) of the nanofluid, in that order, are given as
$$\begin{array}{l}{\rho}_{nf}=(1A){\rho}_{f}+A{\rho}_{s}\\ {\mu}_{nf}=\frac{{\mu}_{f}}{{(1A)}^{2.5}}\\ {(}_{\rho}=(1A){(}_{\rho}+A{(}_{\rho}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(\text{}5\right)\\ k{}_{nf}={k}_{f}\left\{\frac{{k}_{s}+2{k}_{f}2A({k}_{f}{k}_{s})}{{k}_{s}+2{k}_{f}+2A({k}_{f}{k}_{s})}\right\}\end{array}$$
Where A is the solid volume fraction (A ≠ 1), μf is the dynamic
viscosity of the base fluid, while ρ_{f} and ρ_{s} are the densities of
the pure fluid and the nanoparticle respectively. The constants
k_{f} and k_{s} are the thermal conductivities of the base fluid and the
nanoparticle respectively.
Following Rosseland approximation Brewstar [9] the radiative heat flux q_{r} is modeled as
$${q}_{r}=\frac{4{\sigma}^{\prime}}{3{k}^{\prime}}\frac{\partial {T}^{\text{'}4}}{\partial {y}^{\prime}}\left(6\right)$$
Following Rosseland approximation Brewstar [9] the radiative heat flux q_{r} is modeled as
Where σ^{'} is the StefanBoltzman constant and k^{'} is the
mean absorption coefficient. Assuming that the difference in
temperature within the flow is such that T^{ '4} can be expressed
as a linear combination of the temperature, we expand T ^{'4 }in
Taylor's series about T_{∞} as follows:
$${T}^{\text{'}4}={T}_{\infty}^{\text{'}4}+4{T}_{\infty}^{\text{'}3}\left({T}^{\prime}{T}_{\infty}^{\text{'}}\right)+6{T}_{\infty}^{\text{'}2}{\left({T}^{\prime}{T}_{\infty}^{\text{'}}\right)}^{2}+\dots \mathrm{..}\left(7\right)$$
and neglecting higher order terms beyond the first degree in
(
T'  T'_{∞} ), we have
$${T}^{\text{'}4}\approx 3{T}_{\infty}^{;4}+4{T}_{\infty}^{\text{'}3}{T}^{\prime}\left(8\right)$$
Differentiating equation (6) with respect to y' and using
equation (8) to obtain
$$\frac{\partial {q}_{r}}{\partial {y}^{\prime}}=\frac{16{T}_{\infty}^{3}{\sigma}^{\prime}}{3k\text{'}}\frac{{\partial}^{2}{T}^{\prime}}{\partial {y}^{\text{'}2}}\left(9\right)$$
Then equation (3) becomes
$${\left(\rho {c}_{p}\right)}_{nf}\left[u\text{'}\frac{\partial T\text{'}}{\partial {x}^{\prime}}+v\text{'}\frac{\partial T\text{'}}{\partial {y}^{\prime}}\right]={K}_{nf}\frac{{\partial}^{2}T\text{'}}{\partial y{\text{'}}^{2}}+\frac{16{T}_{\infty}^{3}{\sigma}^{\prime}}{3k\text{'}}\frac{{\partial}^{2}{T}^{\prime}}{\partial {y}^{\text{'}2}}\left(10\right)$$
By the introduction of the following variables
$$u=\frac{u\text{'}}{\sqrt{a{v}_{f}}}v=\frac{v\text{'}}{\sqrt{a{v}_{f}}},\theta =\frac{{T}^{\prime}{T}_{\infty}^{\text{'}}}{{T}_{w}^{\text{'}}{T}_{\infty}^{\text{'}}},x=\frac{x\text{'}}{\sqrt{\frac{{v}_{f}}{a}}}.y=\frac{y\text{'}}{\sqrt{\frac{{v}_{f}}{a}}}\left(11\right)$$
Equation (11) transform equation (1), (2) and (10) to the
followings
$$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0\left(12\right)$$
$$u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\frac{1}{(1A){\rho}_{f}+{\rho}_{s}}(\frac{1}{{(1A)}^{2.5}}\frac{{\partial}^{2}u}{\partial {y}^{2}}Mu+{G}_{r}\theta )\left(13\right)$$
$$u\frac{\partial \theta}{\partial x}+v\frac{\partial \theta}{\partial y}=\frac{1}{{p}_{r}}\frac{1}{(1A){(\rho {C}_{p})}_{f}+A{(\rho {C}_{p})}_{s}}\left(\frac{{k}_{nf}}{{k}_{f}}+{R}_{d}\right)\frac{{\partial}^{2}\theta}{\partial {y}^{2}}\left(14\right)$$
Where $M=\frac{\sigma {\beta}_{0}}{a}$
is the magnetic field parameter,${p}_{r}=\frac{{v}_{f}}{{K}_{f}}$
is
the Prandtl number,${R}_{d}=\frac{16\sigma \text{'}{T}_{\infty}^{\text{'}3}}{3{K}^{*}{K}_{f}}$
is the radiation parameter,${G}_{r}=\frac{g{\beta}_{t}\left({T}_{w}^{\text{'}}{T}_{\infty}^{;}\right)x}{a}$
Is the Grashof number
And the boundary conditions becomes
$$\begin{array}{l}u=x,v=0,\text{\hspace{0.17em}}\theta =1\text{\hspace{0.17em}}\text{\hspace{0.17em}}at\text{\hspace{0.17em}}y=0\\ u\to 0,\text{\hspace{0.17em}}\theta \to 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}as\text{\hspace{0.17em}}\text{\hspace{0.17em}}y\to \infty \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{(15)}\end{array}$$
By introducing the stream function ψ, which can be defined as
$$\text{u=}\frac{\partial \text{\psi}}{\partial \text{y}}\text{,v=}\frac{\text{}\partial \text{\psi}}{\partial \text{\hspace{0.17em}}\text{x}}.$$
and using the similarity transforma t i o n s
$$\eta =y,\psi =xf(\eta ),\theta =\theta (\eta )\left(16\right)$$
We have
$${f}^{\u2034}+{(1A)}^{2.5}\left\{f{f}^{\u2033}{(f\text{'})}^{2}\right\}\left[(1A){\rho}_{f}+{\rho}_{s})(M{f}^{\prime}+{G}_{r}\theta )\right]=0\left(17\right)$$
$$\frac{1}{{P}_{r}}\frac{1}{(1A){(\rho {C}_{p})}_{f}+A{(\rho {C}_{p})}_{s}}\left[\frac{{k}_{nf}}{{k}_{f}}+{{\displaystyle R}}_{d}\right]{\theta}^{\u2033}(\eta )+f(\eta ){\theta}^{\prime}(\eta )=0\left(18\right)$$
And the boundary conditions now becomes
$$\begin{array}{l}f(0)=0,\text{}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{f}^{\prime}(0)=1\text{\hspace{0.17em}}\text{\hspace{0.17em}}at\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\eta =0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\\ f\to 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}at\text{\hspace{0.17em}}\text{\hspace{0.17em}}\eta \to \infty .\\ \theta (0)=1\text{\hspace{0.17em}}\text{\hspace{0.17em}}and\text{\hspace{0.17em}}\text{\hspace{0.17em}}\theta (\infty )=0,\end{array}\}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}(19)$$
Numerical Solution
Table 1: Thermo physical properties of water and nanoparticles Hamad
[8].
Compound 
ρ (kg/ m^{3}) 
C_{p}(J/ kgK) 
k (W/ mK 
Pure water 
997.1 
4179 
0.613 
Copper (Cu) 
8933 
385 
401 
Alumina (Al_{2}O_{3 }) 
3970 
765 
40 
Silver (Ag) 
10500 
235 
429 
Titanium Oxide (TiO_{2}) 
4250 
686.2 
8.9538 
An algorithm was written and run with Maple codes to plot graphs and
generate tables for analysis.
The thermo physical properties of pure water and those of
the nanoparticles as given in table 1 by Hamad [8] will now be
substituted in the transformed coupled ordinary differential
equations. These equations will now be solved using Runge
KuttaFehlberg (RKF) method with shooting technique. The
algorithm of RungeKuttaFehlberg (RKF) method is order 5 and
it gives a better approximation of the solution than the general
Runge – Kutta method of order 4.
An algorithm was written and run with Maple codes to plot graphs and generate tables for analysis.
An algorithm was written and run with Maple codes to plot graphs and generate tables for analysis.
Results
Results are shown in figures 2 3 and table 2
From derivation, the the Nuselt number can be written as follows;
$$Nu={\theta}^{\prime}(0)\text{\hspace{0.17em}}\text{\hspace{0.17em}}$$
From derivation, the the Nuselt number can be written as follows;
Discussion of Results
Figure 2:Effect of M on velocity distribution profiles.
Figure 3:Effect of Rd on Temperature distribution profiles.
Figure 4:
Table 2: Values of −θ′(0) for various Rd, A When Gr = 0.1 , M = 0.5
and Pr = 6.2 .
Rd 
A 
θ'(O) 

Cu 
Ag 
Al_{2}CO_{3}_{ } 
TiO_{2} 

0.5 
0.30 
1.31430 
1.23109 
1.28489 
1.28961 

0.35 
1.33022 
1.23834 
1.29732 
1.30086 
0.40 
1.34938 
1.24650 
1.31206 
1.31409 

0.45 
1.37892 
1.25349 
1.33263 
1.33158 

2 
0.30 
0.85353 
0.79483 
0.83119 
0.82915 

0.35 
0.86695 
0.80360 
0.84275 
0.84011 
0.40 
0.87137 
0.80646 
0.84654 
0.84371 

0.45 
0.87576 
0.80929 
0.85030 
0.84727 

4 
0.30 
0.60425 
0.55941 
0.58668 
0.58351 
We have formulated the differential equations modeling the
effect of radiation on heat transfer in nanofluid based on the
already existing model. Numerical evaluation of the numerical
solutions reported in the previous sections was performed and
the results are presented in graphical and tabular forms. This was
done to illustrate the influence of some parameter involved.
The effect of magnetic field parameter on velocity distribution profiles across the boundary layer is presented in Figure 1. It is obvious that the effect of increasing values of the magnetic field parameter M results in a decreasing velocity distribution. This is due to the fact that the introduction of transverse magnetic field normal to the flow direction has a tendency to create a drag force due to Lorentz force and hence results in retarding the velocity profiles Figure 2 shows the influence of the radiation parameter on the temperature distribution profiles. The temperature increases as the radiation increases. Figure 3 presents the temperature distribution profiles for different values of the Prandtl number (Pr). The results show that the effect of increasing values of the Prandtl number results in a decrease in the temperature.
In table 2, as the radiation parameter increases, the nuselt number of the nanoparticles decreases, that is to say the effect of increasing the radiation parameters is to decrease the rate of heat transfer of the nanoparticles.Also, as the fraction volume parameter (A) increases, nuselt number increases.
The effect of magnetic field parameter on velocity distribution profiles across the boundary layer is presented in Figure 1. It is obvious that the effect of increasing values of the magnetic field parameter M results in a decreasing velocity distribution. This is due to the fact that the introduction of transverse magnetic field normal to the flow direction has a tendency to create a drag force due to Lorentz force and hence results in retarding the velocity profiles Figure 2 shows the influence of the radiation parameter on the temperature distribution profiles. The temperature increases as the radiation increases. Figure 3 presents the temperature distribution profiles for different values of the Prandtl number (Pr). The results show that the effect of increasing values of the Prandtl number results in a decrease in the temperature.
In table 2, as the radiation parameter increases, the nuselt number of the nanoparticles decreases, that is to say the effect of increasing the radiation parameters is to decrease the rate of heat transfer of the nanoparticles.Also, as the fraction volume parameter (A) increases, nuselt number increases.
Conclusion
Numerical solution was applied to the differential equation
that model the effect of thermal radiation on heat transfer on
nanofluid using rungekutta Fehlberg method. It was found that
the effect of high radiation is to decrease the rate of heat transfer.
ReferencesTop
 Wang X, Xu X, Choi SUS. Thermal conductivity of nanoparticlefluid mixture. J. Therm. Physics Heat Transfer. 1999;13(4):47480.
 John P, Shima PD and Baldev R. Nanofluid with tunable thermal properties. Appl. Phys. Lett. 2008;92: 043108.
 Khan WA and Pop I. Boundary Layer flow of a nanofluid past a streching sheet. International Journal of Heat and Mass Transfer. 2010;53(1112):24772483.
 Godson L, Raja B, Mohan LD and Wongwises S. Enhancement of Heat Transfer Using nanofluids. Renewable and Sustainable energy Reviews. 2010; 14(2): 629 – 641.
 Zeinali SH, Salehi H, Noie SH. The effect of Magnetic field and nanofluid on thermal performance of twophase closed Thermosyphon (TPCT). Int. Journal of Physical Sciences. 2012;7(4):534 – 543. DOI: 10.5897/ IJPS11.1019.
 Yohannes KY and Shankar B. Heat and mass transfer in MHD flow of nanofliuds through a porous media due to a stretching sheet with viscious disspation and chemical reaction effects. Caribbean journal of science and technology. 2013;1:001017.