The present paper deals with study of free convection in two-dimensional magnetohydrodynamic (MHD) boundary layer flow of an incompressible, viscous, electrically conducting, and steady nanofluid. The governing equations representing fluid flow are transformed into a set of simultaneous ordinary differential equations by using appropriate similarity transformation. The equations thus obtained have been solved numerically using adaptive Runge-Kutta method with shooting technique. The effects of physical parameters like magnetic parameter, temperature buoyancy parameter on relative velocity and temperature distribution profile, shear stress profile, and temperature gradient profile were depicted graphically and analyzed. Significant changes were observed due to these parameters in velocity and temperature profiles.
1. Introduction
Applications of heat transfer through moving material in a moving fluid medium have wide range of applications in real world. Various researchers studied the concepts of moving sheet in a flowing fluid medium. Numerical and analytical solution for momentum and energy in laminar boundary layer flow over continuously moving sheet was discussed by Zheng and Zhang [1]. The authors found that the effect of velocity ratio parameter and other parameters on heat transfer were significant. Al-Sanea [2] studied the flow and thermal characteristic of a moving vertical sheet of extruded material close to and far downstream from the extrusion slot. Regimes of forced, mixed, and natural convection have been delineated, in buoyancy flow, as a function of Reynolds and thermal Grashof numbers for various values of Prandtl number and buoyancy parameter. Seddeek [3] analyzed the effect of magnetic field on the flow of micropolar fluid past a continuously moving plate. Patel et al. [4] suggested a new model for thermal conductivity of nanofluids which is found to agree excellently with a wide range of experimental published data. The momentum and heat transfer flow of an incompressible laminar fluid past a moving sheet based on composite reference velocity was studied by Cortell [5]. The author found that the direction of the wall shear stress changes in such an interval and an increase of the relative velocity parameter yields an increase in temperature. Ishak et al. [6] studied the magnetohydrodynamic boundary layer flow due to a moving extensible surface. It was found that the dual solutions exist for the flow near y-axis, where the velocity profiles show a reversed flow. Bachok et al. [7] investigated steady, laminar boundary layer flow of a nanofluid past a moving semi infinite flat plate. The results indicate that dual solutions exist when the plate and the free stream move in the opposite directions. Habib and El-Zahar [8] proposed a model for determining the heat transfer between a moving sheet and flowing fluid. The authors found that the heat transfer depends on the relative velocity between the moving fluid and the moving sheet to certain values. Energy conversion for conjugate convection, conduction, and radiation analysis have been performed for free convective heat transfer with radiation by Hsiao [9]. The results show that the free convection effect will produce a larger heat transfer effect better than the forced convection. Raguraman et al. [10] investigated the effects of some important design parameters for coal-water slurry in agitated vessel used in coal gasification such as stirrer speed, location of stirrer, D/d ratio, and coal-water ratio. Confirmation test verified that the Taguchi method achieved optimization of heat transfer coefficient in agitated vessel. The objective of the present study is to analyze the effect of magnetic field and temperature buoyancy effect on the heat transfer flow between a moving sheet and moving alumina water nanofluid environment. Hsiao [11, 12] studied two-dimensional conjugate heat transfer with Ohmic dissipation mixed convection of an incompressible Maxwell fluid on a stagnation point and also investigated energy conversion problems of conjugate conduction, convection, and radiation heat and mass transfer with viscous dissipation and magnetic effects. The results have shown that it should produce greater heat transfer effect with larger values of Prandtl number, free convection parameter, and conduction-convection number. On the other hand, the Eckert number and magnetic parameter will reduce the heat transfer effect. For heat convection energy conversion aspect, some importance parameters applied to the system, such as buoyancy parameters, radiative energy parameter, boundary proportional parameter, and Prandtl number which can produce positive effects for larger values of those parameters. Idusuyi et al. [13] formulated a computational model for the heat generation and dissipation in a disk brake during braking and the following release period. The results show similarity in thermal behaviour at the contact surface for the asbestos and aramid brake pad materials with a temperature difference of 1.8 K after 10 seconds.
The objective of the present study is to analyze the effect of magnetic field and temperature buoyancy effect on the heat transfer flow between a moving sheet and moving alumina water nanofluid environment.
2. Mathematical Description
The graphical model of the problem considered has been given to depict the flow configuration and coordinate system. The system deals with two-dimensional MHD boundary layer, incompressible, viscous, electrically conducting, and steady fluid flow (Figure 1). The governing equations representing flow are the following:
(1)∂u∂x+∂v∂y=0(2)u∂u∂x+v∂u∂y=-1ρnf∂p∂x+υnf∂2u∂y2+gβTnfT-T∞-σnfB2(x)ρnfu(3)u∂T∂x+v∂T∂y=αnf∂2T∂y2,
where u and v are the velocity components in x and y direction, respectively, υnf, βTnf, ρnf,σnf, and αnf are the kinematic viscosity, thermal expansion coefficient, density, electrical conductivity, and thermal diffusivity of the nanofluid. g is the gravitational acceleration and T, Tw, and T∞ are the temperature of the nanofluid inside the thermal boundary layer, the plate temperature, and the nanofluid temperature at free stream, respectively. Bx=1/xB0 denotes the magnetic field imposed in the transverse direction to the fluid.
Flow configuration and coordinate system.
The appropriate boundary conditions for the abovementioned problem are as follows (Table 1):(4)u=Us,v=0,T=Ts,aty=0,u⟶U∞,T⟶T∞,asy⟶∞.
To solve the governing equations (1), (2), and (3) with the boundary conditions (4), the following similarity transformation has been introduced:
(5)ψ=υnfU∞-Usx1/2fη,θ=T-T∞Ts-T∞,η=U∞-Usυnfx1/2y.
Here U∞ and Us are free stream velocity and moving sheet velocity, respectively. ψ is the stream function which satisfies (1) with u=∂ψ/∂y and v=-∂ψ/∂x. Thus, (5) gives
(6)u=U∞-Usf′η,v=-12υnfU∞-Usx1/2fη-ηf′η.
Now using (5) and (6) in (2) and (3), we get the following nonlinear differential equations:
(7)f′′′η+12fηf′′η-Mf′η+gtθη=0θ′′+12Prfηθ′η=0;
the boundary conditions (4) will become
(8)f0=0,f′0=11-α,f′′0=unknown,θ0=1,θ′0=unknown,f′η⟶α1-α,θη⟶0,asη⟶∞,
where α=0 represent the state of moving surface in stationary fluid (U∞=0) and 0<α<∞ is the condition of moving surface in moving fluid with U∞>Us. The dimensionless parameters introduced in (7) are defined as follows.
M=σnfB02/ρnf(U∞-Us) is the magnetic parameter, Gr=gβTnfTs-T∞x3/υnf2 is the local temperature Grashof number, Re=U∞x/υnf is the local Reynolds number, gt=Gr/Re2 is the temperature buoyancy parameter, Pr=υnf/αnf is the Prandtl number, and α=U∞/Us is the relative velocity parameter.
To solve the set of nonlinear differential equations (7) with subject to the boundary conditions (8), adaptive Runge-Kutta method with shooting technique has been applied. This method is based on the discretization of the problem domain and the calculation of unknown boundary conditions. The domain of the problem is discretized and the boundary conditions for η=∞ are replaced by f′ηmax=α/abs(1-α), and θηmax=0 where ηmax is sufficiently large value of ηcomparing to step size at which the boundary conditions (8) for f′(η) and θ(η) are satisfied. On account of the consistency and to fulfill stability criteria ηmax=10 and step size Δη=0.01 have been taken. To solve the problem the nonlinear equations (7) are first converted into first order ordinary linear differential equations as follows:
(9)f′η=y2;f′′η=y3;f′′′η=-12y1y3+My2-gty4;θ′=y5;θ′′η=-12Pry1y5.
There are three conditions on the boundary η=0 and two conditions at η=∞ as given in (8). Shooting technique has been used to find required missing initial conditions. The value of unknowns f′′(0) and θ′0 has been tabulated in Table 2. Adaptive Runge-Kutta method estimates the truncation error at each integration step and automatically adjusts the step size to keep the error within prescribed limits. The formula used to measure the actual conservative error is defined by eh=1/5∑i=15Ei2(h)1/2; here Ei(h) is the measure of error in the dependent variable yi. Per-step error control is achieved by adjusting the increment h by hi+1=0.9hi(ε/e(hi))1/5, i=1,2,3,4,5, so that it is approximately equal to the prescribed tolerance ε=10-6. The assumed initial step size is h=0.01.
Values of f′′(0) and θ′(0) for Pr=6.6556406 and ∅=0.05 obtained by shooting method.
Magnetic parameter (M)
Temperature buoyancy parameter (gt)
Wall shear stresses f′′(0)
Skin friction coefficient θ′(0)
0.0
0.2
-0.664300
-1.97611
0.2
0.0
-1.297800
-1.90991
0.2
0.5
-1.157890
-1.91991
0.2
1.0
-1.016890
-1.92921
0.2
3.0
-0.479190
-1.96421
0.5
0.2
-1.768590
-1.85411
1.0
0.2
-2.269990
-1.79411
3.0
0.2
-3.605449
-1.65431
4. Result and Discussion
The computations have been made for velocity, shear stress, concentration, and concentration gradient profiles corresponding to fixed values of magnetic parameter (M), α, and temperature buoyancy parameter (gt). The value of Pr=6.6556406 is taken corresponding to 0.05 volume fraction of the alumina water nanofluid at 20°C. The parameter gt has been introduced to represent the temperature buoyancy effect on flow field. Figures 2–6 exhibit the velocity, shear stress, concentration, and concentration gradient profiles, respectively, for various values of magnetic parameter (i.e., M=0, 0.2, 0.5, 1, 3). It is noticed that in Figures 2 and 3 that the velocity decreases and the shear stress decreases up to (2.93290777, −0.0311476673) and then increases with an increase in the magnetic parameter. The temperature increases with an increase in the magnetic parameter; temperature gradient increases up to (0.524497748, −0.835081441) and then decreases with an increment in the values of magnetic parameter as shown in Figures 4 and 5. Figures 6 and 7 show the effect of convection parameter on velocity, shear stress, concentration, and concentration gradient profiles, respectively, for gt=0,0.2,0.5,1. It has been observed that the velocity increases with an increase in the convection parameter as shown in Figure 6, and shear stress increases up to (0.433856945, −1.03343467) with an increase in the convection parameter and then decreases as shown in Figure 7. The effect of convection parameter on temperature gradient and temperature also has been studied through Figures 8 and 9. It is noticed that the temperature and temperature gradient show significant changes with an increase in the temperature buoyancy parameter.
Effect of magnetic parameter on velocity.
Effect of magnetic parameter on shear stress.
Effect of magnetic parameter on temperature.
Effect of magnetic parameter on temperature gradient.
Effect of temperature buoyancy parameter on velocity.
Effect of temperature buoyancy parameter on shear stress.
Effect of temperature buoyancy parameter on temperature gradient.
Effect of temperature buoyancy parameter on temperature.
5. Conclusion
Free convective heat transfer in an electronically conducting nanofluid over a moving sheet with magnetic field is considered. The nonlinear partial differential equations are transformed into a system of nonlinear ordinary differential equations by using similarity transformation and then solved numerically using the Runge-Kutta fifth order method along with the shooting technique. From the numerical calculations of the skin friction coefficient and wall shear stresses the following is concluded.
An increase in magnetic parameter decreases the velocity and wall shear stress but enhances the temperature, skin friction coefficient. The shear stress decrease up to a fixed value of magnetic parameter and after that goes to increases. For the magnetic parameter, temperature gradient shows the opposite behaviour to the shear stress.
As the temperature buoyancy parameter increases, it increases the velocity, wall shear stress and decreases skin friction coefficient and temperature gradient.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
ZhengL. C.ZhangX. X.Skin friction and heat transfer in power-law fluid laminar boundary layer along a moving surface200245132667267210.1016/S0017-9310(01)00363-5ZBL1101.763232-s2.0-0037030215Al-SaneaS. A.Convection regimes and heat transfer characteristics along a continuously moving heated vertical plate200324688890110.1016/j.ijheatfluidflow.2003.10.0012-s2.0-0345330172SeddeekM. A.Flow of a magneto-micropolar fluid past a continuously moving plate2003306425525710.1016/S0375-9601(02)01513-X2-s2.0-17144434395PatelH. E.SundararajanT.PradeepT.DasguptaA.DasguptaN.DasS. K.A micro-convection model for thermal conductivity of nanofluids20056558638692-s2.0-3084445294110.1007/BF02704086CortellR.Flow and heat transfer in a moving fluid over a moving flat surface200721643544610.1007/s00162-007-0056-zZBL1161.764592-s2.0-35148866524IshakA.NazarR.PopI.MHD boundary-layer flow due to a moving extensible surface2008621233310.1007/s10665-007-9169-zMR24299262-s2.0-48249136943BachokN.IshakA.PopI.Boundary-layer flow of nanofluids over a moving surface in a flowing fluid20104991663166810.1016/j.ijthermalsci.2010.01.0262-s2.0-77955432252HabibH. M.El-ZaharE. R.Mathematical modeling of heat-transfer for a moving sheet in a moving fluid2013633693732-s2.0-84881527002HsiaoK.-L.Energy conversion conjugate conduction-convection and radiation over non-linearly extrusion stretching sheet with physical multimedia effects20135949450210.1016/j.energy.2013.06.0412-s2.0-84883309334RaguramanC. M.RagupathyA.SivakumarL.Experimental determination of heat transfer coefficient in stirred vessel for coal-water slurry based on the taguchi method20132013871929610.1155/2013/719296HsiaoK.-L.Conjugate heat transfer for mixed convection and Maxwell fluid on a stagnation point20143964325433210.1007/s13369-014-1065-zMR3219590HsiaoK.-L.Nanofluid flow with multimedia physical features for conjugate mixed convection and radiation20141041810.1016/j.compfluid.2014.08.001MR3260148IdusuyiN.BabajideI.AjayiO. K.OlugasaT. T.A computational study on the use of an aluminium metal matrix composite and aramid as alternative brake disc and brake pad material20142014649469710.1155/2014/494697EinsteinA.1956New York, NY, USADoverMR0077443DrewD. A.PassmanS. L.1999New York, NY, USASpringerXuanY.RoetzelW.Conceptions for heat transfer correlation of nanofluids200043193701370710.1016/S0017-9310(99)00369-5ZBL0963.760922-s2.0-0343192359MamutE.Characterization of heat and mass transfer properties of nanofluids200651512