# Heat and Mass Transfer on Squeezing Unsteady MHD Nanofluid Flow between Parallel Plates with Slip Velocity Effect

Authors: Khilap Singh; Sawan K. Rawat; Manoj Kumar
Journal: Journal of Nanoscience
 (2016)

---

## Abstract

Heat and mass transfer behavior of unsteady flow of squeezing nanofluids between two parallel plates in the sight of uniform magnetic field with slip velocity effect is investigated. The governing equations representing fluid flow have been transformed into nonlinear ordinary differential equations using similarity transformation. The equations thus obtained have been solved numerically using Runge-Kutta-Fehlberg method with shooting technique. Effects on the behavior of velocity, temperature, and concentration for various values of relevant parameters are illustrated graphically. The skin-friction coefficient and heat and mass transfer rate are also tabulated for various governing parameters. The results indicate that, for nanofluid flow, the rates of heat and mass transfer are inversely proportional to nanoparticle volume fraction and magnetic parameter. The rate of mass transfer increases with increasing values of Schmidt number and squeeze number.


---

## Body

## 1. Introduction

The study of heat and mass transfer in unsteady squeezing viscous nanofluid flow between two parallel plates is a stimulating topic of exploration because of its industrial use and intense biological situations, some of which include processing of polymer, compression, power transmitting, lubricant system, transient loading of mechanical components and the squeezed films in power transmission, food processing, and cooling water, modeling of synthetics transportation inside living bodies, hydromechanical machinery, chemical processing equipment, and crop destruction due to freezing. The first work on the squeezing flow under lubrication approximation was studied by Stefan [1]. The flow analysis between two parallel plates of Cu-water squeezing nanofluid was investigated by Domairry and Hatami [2]. Pourmehran et al. [3] studied the unsteady flow of squeezing nanofluid between parallel plates. Gupta and Ray [4] proposed a problem of unsteady flow of a squeezing nanofluid between two parallel plates. The squeezing flow of Cu-water (or kerosene) nanofluid between two parallel plates under the effects of viscous dissipation and velocity slip was investigated by Khan et al. [5]. Dib et al. [6] obtained an approximate analytic solution of squeezing unsteady nanofluid flow.The word nanofluid represents the fluid in which particles of size with order of nanometer (diameter < 100 nm) are mixed in the base fluid. The nanoparticles used in nanofluids are generally made of metals (Al, Cu), oxides (Al2O3, CuO, TiO2, and SiO2), carbides (SiC), nitrides (AlN, SiN), and nonmetal (graphite, carbon nanotubes) and the base fluid is usually a conductive fluid, such as water or ethylene glycol. Other base fluids are toluene, oil, other lubricants, biofluids, and polymer solution. Nanoparticles are present up to 5% volume fraction in nanofluids. The conventional heat transfer fluids are poor conductors of heat. Nanofluids make an edge over them because they have high heat transfer capability. Since these heating/cooling fluids play a vital role in the development of energy efficient heat transfer equipment for energy supply, to raise the thermal conductivity of these fluids, nanosized conducting metal particles are added to them. Therefore, their proper understanding is a must to use them efficiently in modern industry. Applications of nanofluids include microelectronics, fuel cells, and pharmaceutical processes. Choi and Eastman [7] were the first to propose the term nanofluid that represents the fluid in which nanoscale particles are suspended in the base fluid with low thermal conductivity such as water, ethylene glycol, and oil. In recent years, many researchers have studied and reported nanofluid technology experimentally or numerically in the presence of heat transfer [8–23].Ibrahim and Shankar [24] studied MHD nanofluid flow and heat transfer over a stretching sheet in the presence of thermal radiation and slip conditions. Malvandi and Ganji [25] investigated effect of magnetic field on heat transfer of alumina/water nanofluid inside a circular microchannel. Ul Haq et al. [26] obtained the influence of thermal radiation and slip on MHD nanofluid flow passed over a stretching sheet. Govindaraju et al. [27] solved the problem of magneto-hydrodynamic nanofluid flow on entropy generation in a stretching sheet with slip velocity. The problem based on effects of Stefan blowing and the slips on bioconvection nanofluid flow over a horizontal plate in motion was numerically investigated by Uddin et al. [28]. Hsiao [29] studied the MHD mixed-convection stagnation flow of a nanofluid over stretching sheet in the presence of slip. Kameswaran et al. [30] investigated chemical reaction and viscous dissipation effects on nanofluid flow through stretching or shrinking sheet. Matin and Pop [31] studied heat and mass transfer flow of a nanofluid with chemical reaction in porous channel. Pal and Mandal [32] observed mixed-convection heat and mass transfer stagnation-point flow in nanofluids through stretching/shrinking sheet in a porous medium with thermal radiation. The nanofluid flow and heat transfer in porous medium in the presence of magnetic field and radiation were made by Zhang et al. [33]. Elshehabey and Ahmed [34] analyzed effect of mixed convection in nanofluid flow with sinusoidal distribution of temperature on the both vertical walls using Buongiorno’s nanofluid model.The novelty of the present study is to account for the slip velocity, magnetic field, and mass transfer on squeezing unsteady nanofluid flow and heat transfer between two parallel plates. In this study, authors have applied Runge-Kutta-Fehlberg fourth-fifth-order method with shooting technique to find the solution of nonlinear differential equations. The effects of governing parameters such as slip, magnetic, and squeeze number, Schmidt number, and nanoparticle volume fraction on velocity, temperature, and concentration as well as on skin-friction coefficient, Nusselt, and Sherwood number are investigated. To the best of our best knowledge such investigation is not studied in the scientific literature. Some analytical methods for squeezing unsteady nanofluid flow can be found in [4, 6].


## 2. Mathematical Formulation

We consider an unsteady two-dimensional flow to observe heat and mass transfer of a squeezing nanofluid in the middle of two parallel plates extended infinitely and implanted in a system occupied with nanofluid (water as a base fluid) containing different types of nanoparticles, that is, copper (Cu), silver (Ag), alumina (Al2O3), and titanium oxide (TiO2) with slip velocity effect. The thermophysical properties of the nanofluids are given in Table 1. A transverse magnetic field of variable strength is imposed in direction perpendicular to both the plates. The distance between two plates is z=±l1-αt1/2=±ht, where l is the initial position (at time t=0). Flow is incompressible with no chemical reaction in system. Further, viscous dissipation effects are retained. The graphical model support to the present study has been given in Figure 1. The governing equations representing flow are as follows: (1)∂u∂x+∂v∂y=0,(2)ρnf∂u∂t+u∂u∂x+v∂u∂y=-∂p∂x+μnf∂2u∂x2+∂2u∂y2-σnfB02u,(3)ρnf∂v∂t+u∂v∂x+v∂v∂y=-∂p∂y+μnf∂2v∂x2+∂2v∂y2,(4)∂T∂t+u∂T∂x+v∂T∂y=knfρCpnf∂2T∂x2+∂2T∂y2+μnfρCpnf4∂u∂x2+∂u∂y+∂u∂x2,(5)∂C∂t+u∂C∂x+v∂C∂y=D∂2C∂x2+∂2C∂y2.The associated boundary conditions are given as (6)u=-L∂u∂y,v=vw=dhdt,T=Th,C=Chaty=ht,v=∂u∂y=∂T∂y=0,C=C0aty=0, where(7)ρnf=1-φρf+φρs,ρCpnf=1-φρCpf+φρCps,μnf=μf1-φ2.5Brinkman,knfkf=ks+2kf-2φkf-ksks+2kf+2φkf-ksMaxwell-Garnetts,σnf=1-φσf+φσs.Table 1
Thermophysical properties of pure water and nanoparticles.
ρ (Kg/m3) C p (J/kgK) k (W/mK) β × 1 0 5  ( K - 1 ) Pure water 997.1 4179 0.613 21 Copper (Cu) 8933 385 401 1.67 Silver (Ag) 10,500 235 429 1.89 Alumina (Al2O3) 3970 765 40 0.85 Titanium oxide (TiO2) 4250 686.2 8.9538 0.9Figure 1
Flow configuration and coordinate system.Equations (2)–(5) can be converted to a system of nonlinear ordinary differential equations via the following similarity variables:(8)η=yl1-αt1/2,u=αx21-αtf′η,v=-αl21-αt1/2fη,θ=T-T0Th-T0,ϕ=C-C0Ch-C0.The transformed equations are(9)fiv-SA11-φ2.5ηf′′′+3f′′+f′f′′-ff′′′-Mf′′=0θ′′+Pr⁡SA2A3θ′f-ηθ′+Pr⁡EcA31-φ2.5f′′2+4j2f′2=0,ϕ′′+ScSfϕ′-ηϕ′=0,where A1, A2, and A3 are dimensionless constants defined as follows:(10)A1=1-φ+φρsρf,A2=1-φ+φρCpsρCpf,A3=knfkf.The boundary conditions (6) in the terms of similarity variables (8) become(11)f0=0,f′′0=0,θ′0=0,ϕ0=0,f1=1,f′1=-δf′′1,θ1=1,ϕ1=1,where S=αl2/2νf is the squeeze number, Pr⁡=μfρCpf/ρfkf is the Prandtl number, Sc=υf/Dnf is Schmidt number, M=2σnfB02ht2/ρfμnf is the magnetic parameter, j=l/x is the reference length and δ=L/l1-αt1/2 is the velocity slip parameter, and Ec=ρf/ρCpfαx/21-αt2 is the Eckert number.The physical quantities of interest are the skin-friction coefficientCf, the Nusselt number Nux, and the Sherwood number Shx defined as(12)Cf=τwρnfvw2,Nux=lqwkfTh-T0,Shx=lmwDnfCh-C0,where(13)τw=μnf∂u∂yy=ht,qw=-knf∂T∂yy=ht,mw=-Dnf∂C∂yy=ht.Using (8) and (13) in (12), we get(14)Cf∗=x21-αtRexCfl2=f′′1A11-ϕ2.5,Nux∗=1-αtNux=-A3θ′1,Shx∗=1-αtShx=-ϕ′1,where Rex=αl5/2x31-αt1/2υf is the local Reynolds number.


## 3. Method of Solution

In this present paper, Runge-Kutta-Fehlberg fourth-fifth-order method (RKF45) has been employed to solve the system of nonlinear ordinary differential (9) with the boundary conditions given by (11) for different values of governing parameters. The RKF45 method has a procedure to determine if the suitable step size h is being used. The formula of fifth-order Runge-Kutta-Fehlberg method can be defined as follows:(15)zn+1=zn+16135k0+665612825k2+2856156430k3-950k4+255k5h,where the coefficients k0 to k5 are defined as follows:(16)k0=fxn,yn,k1=fxn+14h,yn+14hk0,k2=fxn+38h,yn+332k0+932k1h,k3=fxn+1213h,yn+19322197k0-72002197k1+72962197k2h,k4=fxn+h,yn+439216k0-8k1+3680513k2-8454104k3h,k5=fxn+12h,yn+-827k0+2k1-35442565k2+18594104k3-1140k4h.The computation of the error can be achieved by subtracting the fifth-order from the fourth-order method, (17)yn+1=yn+25216k0+14082565k2+21974104k3-15k4h.If the error goes beyond a specified antechamber, the results can be recalculated using a smaller step size. The approach to computing the new step size is shown as follows:(18)hnew=holdεhold2zn+1-yn+11/4.The variation of the dimensionless velocity, temperature, and concentration is ensured to be less than 10-6 between any two consecutive iterations for the convergence criterion.To solve the nonlinear differential equations (9) subject to the boundary conditions (11), first boundary conditions for η=1 are replaced by f(1)=1,f′1=-δf′′1,θ1=1, and ϕ1=1. We consider that f=f1,f′=f2,f′′=f3,f′′′=f4,θ=f5,θ′=f6,ϕ=f7, and ϕ′=f8. The nonlinear equations (9) are first converted into first-order ordinary linear differential equations as follows:(19)f′=f2,f2′=f3,f3′=f4,f4′=SA11-φ2.5ηf4+3f3+f2f3-f1f4+Mf3,f5′=f6,f6′=-Pr⁡SA2A3f6f1-ηf6-Pr⁡EcA31-φ2.5f32+4j2f22,f7′=f8,f8′=-ScSf1f8-ηf8,subject to the following initial conditions:(20)f10=0,f20=α1,f30=0,f40=α2,f50=α3,f60=0,f70=0,f80=α4.We ran the computer code written in the MATLAB for different values of step length ∇η. We found that there is no or only a negligible change in the physical quantities of interest like the skin-friction coefficient, the couple stress coefficient, Nusselt number, and Sherwood number for various values of ∇η>0.01. Therefore, in the present paper we have set step size ∇η=0.01. There are four initial conditions at η=0 and four conditions on boundary η=1. To get the solution of the problem, four more initial conditions at η=0, that is, values of α1, α2, α3, and α4, are required, which have been obtained by shooting technique. Finally the transformed initial value problem is solved by employing the Runge-Kutta-Fehlberg fourth-fifth-order method along with calculated boundary conditions.


## 4. Results and Discussion

In order to validate the numerical results obtained, we compare our results with those reported by Gupta and Ray [4] as shown in Table 2, and they are found to be in a good agreement. The effects of the volume fraction of solid nanoparticles, magnetic parameter, velocity slip parameter, squeeze number, and Schmidt number are inspected for different kinds of nanoparticles when the base fluid is water, Ec=0.01,Pr⁡=6.2, and j=0.01.Table 2
Comparison offη and θ(η) for Cu-water nanofluid when S=1,Pr⁡=6.2,Ec=j=0.01,φ=0.02, and M=Sc=δ=0.
η Gupta and Ray [4] Present result f ( η ) θ ( η ) f ( η ) θ ( η ) 0 0 1.03206637 0 1.02996637 0.1 0.14135866 1.03206407 0.14135879 1.02996428 0.2 0.28066605 1.03203202 0.28066639 1.02993531 0.3 0.41578075 1.03189263 0.41578137 1.02980932 0.4 0.54437882 1.03150811 0.54437979 1.02946175 0.5 0.66385692 1.03066423 0.66385837 1.02869896 0.6 0.77122923 1.02903983 0.77123132 1.02723066 0.7 0.86301562 1.02615205 0.86301853 1.02462037 0.8 0.93511971 1.02125916 0.93512364 1.02019765 0.9 0.98269524 1.01318608 0.98270044 1.01290031 1.0 1.00000000 1.00000000 1.00000000 1.00000000Figures2 and 3 illustrate the behavior of the velocity f′(η) and temperature θ(η) of the Cu-water nanofluid for different values of magnetic parameter M. Figure 2 shows that, with an increase in the values of M, the velocity decreases near the lower plate surface, but after a certain distance it increases. It is noticed from Figure 3 that the temperature decreases monotonically with increasing values of M.Figure 2
Effects ofM on velocity f′(η).Figure 3
Effects ofM on temperature θ(η).The variations of the velocity, temperature, and concentration for different values of slip parameterδ are shown in Figures 4–6. It is noted from Figure 4 that near the wall the value of velocity f′(η) increases with rising values of δ when 0≤η≤0.6 and for η≥0.6 the velocity decreases as δ increases. Figure 5 depicts that a rise in the values of slip parameter δ increases temperature monotonically. It is evident from Figure 6 that the change in the concentration ϕ(η) with rising values of slip parameter δ is negligible.Figure 4
Effects ofδ on velocity f′(η).Figure 5
Effects ofδ on temperature θ(η).Figure 6
Effects ofδ on concentration ϕ(η).The effects of the squeeze numberS on velocity, temperature, and concentration profiles are depicted in Figures 7–9. Physically the squeeze number (S) describes the movement of the plates (S>0 corresponds to the plates moving apart, while S<0 corresponds to the plates moving together). It can be easily seen from Figure 7 that the value of velocity f′(η) near the lower plate surface decreases regularly with the increase in the value of S, and as we move away from lower plate surface this value increases. Figure 8 shows that increasing value of squeeze number S decreases the temperature θ(η). Figure 9 depicts that the value of ϕ(η) goes up with increasing η and high squeeze number implies slight drop in concentration ϕ(η).Figure 7
Effects ofS on velocity f′(η).Figure 8
Effects ofS on temperature θ(η).Figure 9
Effects ofS on concentration ϕ(η).Figures10 and 11 show the effects of the volume fraction φ on the velocity and temperature profiles. Figure 10 exhibits initially an increase in the values of nanoparticles volume fraction φ, the velocity f′(η) decreases, and after a fixed distance from lower plate surface it increases slightly. From Figure 11, it is observed that the temperature decreases with increasing values of φ.Figure 10
Effects ofφ on velocity f′(η).Figure 11
Effects ofφ on temperature θ(η).Figures12 and 13 illiterate the effect of the different nanoparticles on velocity and temperature profile when the base fluid is water. It is observed from Figure 12 that the different nanofluids have different velocities and also noted that titanium oxide (TiO2) has higher velocity when compared to other nanoparticles such as Al2O3, Cu, and Ag for 0≤η≤0.5, while for η>0.5 the results get reversed. It is noted from Figure 13 that the order of nanofluid for decreasing value of θ(η) when η change from 0 to 1 is TiO2-water, Al2O3-water, Cu-water, and Ag-water nanofluid. Figure 14 depicts the effect of the Schmidt number Sc on concentration profile. The values of the temperature increase with the increase in the value of Sc.Figure 12
Velocity of different nanofluids.Figure 13
Temperature of different nanofluids.Figure 14
Effects of Sc on concentrationϕ(η).The effects of the squeeze numberS on the skin-friction coefficient Cf∗, the Nusselt number Nux∗, and the Sherwood number Shx∗ are given in Table 3. From Table 3, it is obvious that the skin-friction coefficient and the Nusselt number are inversely proportional to S, whereas the Sherwood number is directly proportional to S.Table 3
Variation ofCf∗,Nux∗, and Shx∗ with different values of Squeeze number S for Cu-water nanofluid when M=1,Sc=1,δ=0.1, and φ=0.02.
S C f ∗ Nu x ∗ Sh x ∗ - 1.0 - 1.34295 0.187977 - 1.0694254 - 0.5 - 1.77787 0.127657 - 1.0293233 0.0 - 2.09419 0.104693 - 1.0000000 0.5 - 2.34148 0.094096 - 0.9770841 1.0 - 2.54426 0.088577 - 0.9583107Table4 displays the effects of the skin-friction coefficient, the Nusselt number, and Sherwood number for different values of the magnetic parameter M and slip parameter δ. It is noticed from the table that the effect of increasing values of M is to increase the skin-friction coefficient Cf∗, the heat transfer rate Nux∗, and mass transfer rate Shx∗. Further, from Table 4 it is concluded that the increasing value of δ decreases the skin-friction coefficient and increases the heat and mass transfer rate.Table 4
Variation ofCf∗,Nux∗, and Shx∗ with different values of magnetic parameter M and slip parameter δ for Cu-water nanofluid when S=1,Sc=1, and φ=0.02.
M δ C f ∗ Nu x ∗ Sh x ∗ 0 0.1 - 2.46084 0.089198 - 0.95676 1 0.1 - 2.54426 0.088577 - 0.95831 2 0.1 - 2.6223 0.088075 - 0.95974 1 0.2 - 2.95257 0.118598 - 0.95158 1 0.5 - 3.26534 0.144432 - 0.94633The effects of the nanoparticle volume fractionφ and Schmidt number Sc on the skin-friction coefficient Cf∗, Nusselt number (the heat transfer rate) Nux∗, and Sherwood number (the mass transfer rate) Shx∗ are given in Table 5. From this table it is concluded that the increasing value of φ increases the skin-friction coefficient and decreases the heat and mass transfer rate. In addition, the rate of mass transfer increases with the increase in Sc, but there is no effect of Sc on the skin-friction coefficient and the heat transfer rate.Table 5
Variation ofCf∗,Nux∗, and Shx∗ with different values of volume fraction φ and Schmidt number Sc for Cu-water nanofluid when M=1,S=1, and δ=0.1.
φ Sc C f ∗ Nu x ∗ Sh x ∗ 0.00 1 - 2.76550 0.092885 - 0.95773 0.02 1 - 2.54426 0.088577 - 0.95831 0.05 1 - 2.31844 0.082295 - 0.95911 0.02 0 - 2.54426 0.088577 - 1.00000 0.02 10 - 2.54426 0.088577 - 0.63318 0.02 50 - 2.54426 0.088577 - 0.05952 0.02 100 - 2.54426 0.088577 - 0.00182The comparison for metallic nanoparticles (Ag, Cu) and nonmetallic nanoparticles (Al2O3, TiO2) is done in Table 6. It is observed form Table 6 that value of the skin-friction coefficient is more for the metallic nanoparticles than nonmetallic nanoparticles, whereas the nonmetallic nanoparticles have higher heat and mass transfer rate in comparison with metallic nanoparticles.Table 6
Variation ofCf∗,Nux∗, and Shx∗ for different nanoparticles when M=1,S=1,Sc=1,δ=0.1, and φ=0.02.
Nanoparticle C f ∗ Nu x ∗ Sh x ∗ Ag - 2.48678 0.088535 - 0.95851 Cu - 2.54426 0.088577 - 0.95831 TiO2 - 2.73607 0.088742 - 0.95777 Al2O3 - 2.74807 0.088794 - 0.95772


## 5. Conclusions

The present paper deals with the numerical solution of combined heat and mass transfer effects of an unsteady MHD laminar two-dimensional flow of incompressible viscous nanofluids in middle of two parallel plates extended infinitely with slip velocity effect. The relevant nonlinear partial differential equations were transformed to a set of ordinary differential equations and then are solved numerically using the Runge-Kutta-Fehlberg fourth-fifth-order method along with shooting technique. From the above discussion the eye-catching results are as follows:(i)
Temperature drops in Cu-water nanofluid with increase in the magnetic field strength.(ii)
The velocity of the nanofluid with nonmetallic nanoparticles Al2O3 and TiO2 is greater than metallic nanoparticles Cu and Ag initially when water is the base fluid, but nature gets reversed after η=0.5. The temperature of the metallic nanoparticles is lower than the nonmetallic nanoparticles.(iii)
When plates start moving apart, then temperature starts decreasing.(iv)
Initially, with the increasing values of magnetic field strength and squeeze number, the velocity of Cu-water slows down but the nature in both the case gets reversed afterη=0.5.(v)
The Nusselt number and skin-friction coefficient decrease as plates move apart for Cu-water nanofluid.(vi)
Nusselt number and Sherwood number for nonmetallic nanoparticles are higher than the metallic nanoparticles.(vii)
Varying value of slip parameter has negligible effect on concentration and as values of squeeze number increase, the value of concentration profile also gets slightly higher.(viii)
Considerable amount of enhancement in concentration can be seen as Schmidt number rises.



---

*Source: 9708562-2016-11-30.xml*