# The Effect of Heat Transfer on MHD Marangoni Boundary Layer Flow Past a Flat Plate in Nanofluid

Authors: D. R. V. S. R. K. Sastry; A. S. N. Murti; T. Poorna Kantha
Journal: International Journal of Engineering Mathematics
 (2013)

---

## Abstract

The problem of heat transfer on the Marangoni convection boundary layer flow in an electrically conducting nanofluid is studied. Similarity transformations are used to transform the set of governing partial differential equations of the flow into a set of nonlinear ordinary differential equations. Numerical solutions of the similarity equations are then solved through the MATLAB “bvp4c” function. Different nanoparticles like Cu, Al2O3, and TiO2 are taken into consideration with water as base fluid. The velocity and temperature profiles are shown in graphs. Also the effects of the Prandtl number and solid volume fraction on heat transfer are discussed.


---

## Body

## 1. Introduction

The convection induced by the variations of the surface tension gradients is known as the Marangoni convection. This convection has received great consideration in view of its application in the fields of welding and crystal growth. Also this convection is necessary to stabilize the soap films and drying silicon wafers. During the study of the existence of the steady dissipative layers which occur along the liquid-liquid or liquid-gas interfaces, Napolitano [1] first called the boundary layer as the Marangoni boundary layer. Many researchers such as Okano et al. [2], Christopher and Wang [3], Pop et al. [4] and Magyari and Chamkha [5] have investigated the Marangoni convection in various geometries. Al-Mudhaf and Chamkha [6] obtained the similarity solution for the MHD thermosolutal Marangoni convection over a flat surface in the presence of heat generation or absorption with fluid suction and injection. Chen [7] investigated the flow and the heat transfer characteristics on the forced convection in a power law liquid film under an applied Marangoni convection over a stretching sheet. In recent years, the study on convective transport of nanofluids has become one of the popular topics of interest. Nanotechnology takes important part for the development of high performance, compact, and cost-effective liquid cooling systems. Moreover, nanofluids have effective applications in many industries such as electronics, transportation, biomedical, and many more. Nanotechnology has been an ongoing topic of discussion in public health as some of the researchers claimed that nanoparticles could present possible dangers in health and environment. Jang and Choi [8] have introduced nanosized particle in a base fluid, which is also termed nanofluid, for the first time. Arifin et al. [9] have examined the influence of nanoparticles on the Marangoni boundary layer flow using a model proposed by Tiwari and Das [10]. An extended work was done by Buongiorno [11], Daungthongsuk and Wongwises [12], Trisaksri and Wongwises [13], Wang and Mujumdar [14], and Kakaç and Pramuanjaroenkij [15]. Recently Hamid et al. [16] studied the radiation effects on the Marangoni boundary layer flow past a flat plate in nanofluid. In the present paper, we study a numerical solution of MHD heat transfer problem in nanofluid with nanoparticles Cu, Al2O3, and TiO2. We also observed the effects of the Prandtl number and solid volume fraction on the Nusselt number. The results are shown graphically.


## 2. Mathematical Formulation

Consider a steady two-dimensional Marangoni boundary layer flow past a permeable flat plate in a water-based nanofluid containing different types of nano particles like Cu (Copper), Al2O3 (Aluminium Oxide), and  TiO2 (Titanium dioxide). Assume that the fluid is incompressible and the flow is laminar. Also it is assumed that the base fluid and the particles are in thermal equilibrium and no slip occurs between them. The thermophysical properties of nanoparticles are given in the Table 1. Further, we consider a Cartesian coordinate system (x,y), where x and y are the coordinates measured along the plate and normal to it, respectively, and the flow takes place at y≥0. Assume that the temperature of the plate is Tw(x) and that of the ambient fluid is T∞.Table 1
Thermophysical properties of pure water and nanoparticles (Oztop and Abu-Nada [17]).
Physical property
Pure water
Cu
Al2O3
TiO2
ρ (kg/m3)
997.1
8933
3970
4250
C
p (J/kg K)
4179
385
765
686.2
k (W/m K)
0.613
401
40
8.9538We further assume that the surface tensionσ is to vary linearly with temperature as
(1)σ=σ0[1-γ(T-T∞)],
where  σ0 is the surface tension at the interface and γ is the rate of change of surface tension with temperature (a positive fluid property). It is also assumed that a uniform magnetic field, H0 is imposed in the direction normal to the surface (Figure 1). Then, the steady state boundary layer equations for a nanofluid in the Cartesian coordinates are given by
(2)∂u∂x+∂v∂y=0,(3)u∂u∂x+v∂u∂y=μnfρnf∂2u∂y2-σρnfH02u,(4)u∂T∂x+v∂T∂y=αnf∂2T∂y2,
together with the boundary conditions
(5)v=0,T=T∞+ax2,μnf∂u∂y=∂σ∂T∂T∂xaty=0,u=0,T=T∞asy⟶∞.Figure 1
Schematic diagram of the problem.Hereu and v are the components of velocity along the x- and y-axes, respectively. T is the temperature, αnf is the thermal diffusivity, ρnf is the effective density, knf is the effective thermal conductivity, and μnf is the effective viscosity of the nanofluid. Moreover, a is the coefficient of temperature gradient. Consider the following:
(6)αnf=knf(ρCp)nf,ρnf=(1-ϕ)ρf+ϕρs,μnf=μf(1-ϕ)2.5,(ρCp)nf=(1-ϕ)(ρCp)f+ϕ(ρCp)s,knfkf=(ks+2kf)-2ϕ(kf-ks)(ks+2kf)+ϕ(kf-ks),a=ΔTL2,
where ϕ is the solid volume fraction of the nanofluid, ρf is the reference density of the fluid fraction, ρs is the reference density of the solid fraction, μf is the viscosity of the fluid fraction, kf is the thermal conductivity of the fluid, ks is the thermal conductivity of the solid, and (ρCp)nf is the heat capacity of the nanofluid. L is the length of the surface, and ΔT is the constant characteristic temperature.A similarity solution of (2)–(5) is obtained by introducing an independent variable η and a dependent variable f in terms of the stream function ψ as
(7)ψ=C1xf(η),η=C2y,
with  u=∂ψ/∂yand v=-∂ψ/∂x.The constantsC1 and C2 are given by
(8)C1=(σ0γaμfρf2)1/3,C2=(σ0γaρfμf2)1/3.Further, the dimensionless temperatureθ  is given by
(9)θ(η)=T-T∞ax2.Substituting (6), (7), (8), and (9) into (3) and (4), we obtain a set of nonlinear ordinary differential equations:
(10)f′′′=(1-ϕ)2.5[(1-ϕ)+ϕρsρf](f′2-ff′′)+M(Pr)1/3(kf)1/3(Cpf)-1/3(ρf)-2/3f′,θ′′=[(1-ϕ)+ϕ(ρsCps/ρfCpf)]knf/kf(2f′θ-fθ′)Pr,
and the boundary conditions become
(11)f(0)=0,θ(0)=1,1(1-ϕ)2.5f′′(0)=-2,f′(∞)=0,θ(∞)=0,
where the magnetic field parameter M=σ1/3H02/(γa)2/3.Also one can define the surface velocity and the local Nusselt number, respectively, as(12)uw(x)=(σ0γa)2ρfμf3xf′(0),(13)Nux=xqw(x)kf[T(x,0)-T(x,∞)],
where qw(x) is the heat flux from the surface of the plate and is given by
(14)qw(x)=-knf(∂T∂y)y=0.Using the above nondimension quantities, one can obtain the local Nusselt number as(15)Nux=-knfkfC2θ′(0).Based on the average temperature difference between the temperature of the surface and the ambient fluid temperature we define(16)NuL=-knfkf(MaLPr)1/3θ′(0),
where MaL is the Marangoni based on L and is defined as
(17)MaL=(∂σ/∂T)(ΔT)Lαfμf.


## 3. Results and Discussion

Numerical solutions were obtained for the effect of the Prandtl number and solid volume fraction on the Marangoni heat transfer in a nanofluid. In this paper, we considered three different nanoparticles whose thermophysical properties were given in Table1. The nonlinear ordinary differential equations (10) subject to the boundary conditions (11) were solved numerically using the MATLAB “bvp4c” routine. We considered the range of nanoparticles volume fraction ϕ as 0≤∅≤0.3 and the Prandtl number Pr as 2≤Pr≤8 (for the base fluid (water) Pr=6.2). The influences of the magnetic field parameter (M), the nanoparticles volume fraction (ϕ) on velocity and, temperature and also the influence of the Prandtl number (Pr) and solid volume fraction (ϕ) on the Nusselt number are presented in graphs.Figures2, 3, and 4 display the velocity profiles, and Figures 5, 6, and 7 display the temperature profiles of Cu-water, Al2O3-water, and TiO2-water, respectively, for different values of magnetic field parameter M. It is observed from the figures that the velocity in the boundary layer decreases and temperature increases as the Magnetic field parameter increases; this is due to the resistive force, called the Lorentz force, which is produced by the induced magnetic field within the boundary layer.Figure 2
Velocity profile for Cu nanoparticles for variousM.Figure 3
Velocity profile for Al2O3 nanoparticles for various M.Figure 4
Velocity profile for TiO2 nanoparticles for various M.Figure 5
Temperature profile for Cu nanoparticles for variousM.Figure 6
Temperature profile for Al2O3 nanoparticles for various M.Figure 7
Temperature profile for TiO2 nanoparticles for various M.Figure8 depicts the influence of volume fraction on the velocity profile of the nanofluid particles. It is observed near the wall that velocity decreases with an increase in the volume fraction ϕ.  Also it is observed that the velocity of TiO2 nanoparticles is higher than that of Cu nanoparticles. From Figure 9, it is clear that an increase in the value of volume fraction enhances the temperature profile, and Cu nanoparticles exhibit more temperature than that of the other nanoparticles. It is also known from Figure 10 that temperature decreases with an increase in the Prandtl number. This is because of a decrease in thermal diffusivity with an increase in the Prandtl number (Pr).Figure 8
Velocity profile for differentϕ.Figure 9
Temperature profile for differentϕ.Figure 10
Temperature profile for different Pr.Figures11 and 12 depict the influence of the Prandtl number and volume fraction on heat transfer, respectively. It is observed that the Nusselt number increases with an increase in Prandtl number and decreases with an increase in the volume fraction.Figure 11
Heat transfer effect against the Prandtl number.Figure 12
Heat transfer effect against the volume fraction.From Table2, it is observed that skin friction decreases with the increase in the volume fraction.Table 2
The numeric values of skin friction and the Nusselt number for various values ofϕ.
ϕ
f
′
′
(
0
)
-
θ
′
(
0
)
0
2
4.68429586
0.1
1.53686694
3.40361315
0.2
1.1448668
2.58644151


## 4. Conclusion

In the present paper, we studied the effect of heat transfer on the Marangoni boundary layer flow past a flat plate in nanofluids in presence of transverse magnetic field. With the similarity transformation, the governing equations of motion together with boundary conditions were transformed to a set of nonlinear ordinary differential equations. The numerical solutions are then obtained for these equations by the help of MATLAB “bvp4c” programming tool. Different types of nanoparticles like Cu, Al2O3, and TiO2 were taken into consideration with H2O as base fluid. The effects of magnetic field parameter M, solid volume fraction of the nanofluid ϕ on the velocity and temperature fields for different nano particles, and the Prandtl number Pr on temperature field were plotted and analyzed. Also the effects of the Prandtl number Pr and solid volume fraction ϕ on local the Nusselt number for the different nanoparticles were discussed for a fixed value of magnetic field parameter M.  It is found that the inclusion of the magnetic field parameter on the flow increased the temperature and decreased the velocity fields in all types of nanofluids. A similar profile was observed on the inclusion of solid volume fraction of the nanoparticles. It was noted that presence of the Prandtl number reduced the temperature field. Also it was observed that for a fixed Prandtl number and other parameters, the rate of heat transfer is more in TiO2-H2O.



---

*Source: 581507-2013-05-25.xml*