Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 692728, 14 pages

http://dx.doi.org/10.1155/2014/692728

Research Article

Analytical Modelling of Three-Dimensional Squeezing Nanofluid Flow in a Rotating Channel on a Lower Stretching Porous Wall

Navid Freidoonimehr,¹ Behnam Rostami,¹ Mohammad Mehdi Rashidi,² and Ebrahim Momoniat³

¹Young Researchers & Elite Club, Islamic Azad University, Hamedan Branch, Hamedan 65181 15743, Iran

²Mechanical Engineering Department, Engineering Faculty of Bu-Ali Sina University, Hamedan 65178 38695, Iran

³Centre for Differential Equations, Continuum Mechanics and Applications, School of Computational and Applied Mathematics, University of the Witwatersrand, Private Bag 3, Johannesburg 2050, South Africa

Received 3 July 2014; Accepted 30 July 2014; Published 17 August 2014

Academic Editor: Sandile Motsa

Copyright © 2014 Navid Freidoonimehr et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A coupled system of nonlinear ordinary differential equations that models the three-dimensional flow of a nanofluid in a rotating channel on a lower permeable stretching porous wall is derived. The mathematical equations are derived from the Navier-Stokes equations where the governing equations are normalized by suitable similarity transformations. The fluid in the rotating channel is water that contains different nanoparticles: silver, copper, copper oxide, titanium oxide, and aluminum oxide. The differential transform method (DTM) is employed to solve the coupled system of nonlinear ordinary differential equations. The effects of the following physical parameters on the flow are investigated: characteristic parameter of the flow, rotation parameter, the magnetic parameter, nanoparticle volume fraction, the suction parameter, and different types of nanoparticles. Results are illustrated graphically and discussed in detail.

1. Introduction

Due to the vast number of applications in chemical as well as mechanical engineering processes such as the manufacture of thin plastic sheets, insulating materials, paper fabrication, and other various processes, this paper pays particular attention to the investigation of rotating flow over stretching surfaces [1].

Sakiadis [2, 3] initiated the study of boundary-layer flow over a continuous solid surface that moves with constant speed. Considering the effect of suction/injection, Erickson et al. [4] studied heat and mass transfer over a moving surface with constant surface velocity and temperature. Tsou et al. [5] studied heat transfer effects of a moving solid surface with constant velocity and temperature. Crane [6] investigated the two-dimensional flow of a viscous fluid over a stretching wall. Andersson [7] investigated MHD effects on boundary-layer flow of a viscoelastic fluid flow past a stretching sheet. Prasad et al. [8] employed a fourth order Runge-Kutta integration scheme to investigate the effect of variable fluid viscosity, magnetic parameter, Prandtl number, variable thermal conductivity, heat source/sink parameter, and thermal radiation parameter on MHD fluid flow over a stretching sheet. The generalized three-dimensional flow and heat transfer over a stretching sheet and in a channel bounded by the lower stretching plate and upper permeable wall were studied by Mehmood and Ali [9, 10]. Three-dimensional flow in a channel with a stretching wall was investigated by Borkakoti and Bharali [11]. Munawar et al. [12] analyzed the slip effect on the flow in a channel bounded by two stretching disks analytically.

Flow squeezed by two parallel plates has been investigated by various researchers. Chamkha et al. [13] determined an analytical solution for the problem of fully developed free convective flow of a micropolar fluid between two vertical parallel plates. Bhargava et al. [14] investigated the fully developed flow and heat transfer of an electrically conducting micropolar fluid (a strong cross magnetic field) between two parallel and porous plates in which the temperature has been considered to be dependent on the heat source, including the effect of frictional heating. The quasilinearization method is used to solve the governing system of ordinary differential equations. Ariel [15] presented solutions for two problems of laminar forced convection of a second-grade (viscoelastic) fluid through two parallel porous walls considering rectangular and cylindrical geometries. Applying similarity transformations to the governing nonlinear partial differential equations, Hayat and Abbas [16] derived nonlinear ordinary differential equations and studied the two-dimensional boundary-layer flow of an upper-convected Maxwell fluid in a channel with chemical reaction. Domairry and Aziz [17] obtained an analytic solution for unsteady MHD squeezing flow with suction and injection effects by the use of homotopy perturbation method.

To the best of authors’ knowledge, Choi and Eastman [18] were probably the first researchers who employed a mixture of nanoparticles and base fluid and called this mixture a “nanofluid.” A wide range of review papers have been published on nanofluids in recent years. Xuan and Li [19] considered the Reynolds number and volume fraction of nanoparticle influences in turbulent flows for nanofluids in tubes experimentally. Bachok et al. [20] employed a numerical, Keller-Box technique for steady nanofluid flow over a porous rotating disk. Khan and Pop [21] studied laminar flow for a nanofluid across a stretching flat surface using an implicit finite difference method. Abolbashari et al. [22] employed HAM to study the entropy analysis in an unsteady MHD nanofluid regime adjacent to an accelerating stretching permeable surface. Beg et al. [23] presented a comparative numerical solution for both single- and two-phase models for bionanofluid transport phenomena. Rashidi et al. [24] compared the two phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel numerically. Abu-Nada et al. [25] illustrated the impacts of variable properties in natural convection nanofluid flow. Rashidi et al. [26] showed how the second law of thermodynamics can be applied to MHD incompressible nanofluid flow over a porous rotating disk. The stagnation flow for a nanofluid over a stretching sheet was studied by Mustafa et al. [27] analytically. The interested reader is referred to the following papers for further reading on the application of nanoparticles in fluid flow [28–30].

In this paper we derive a coupled system of nonlinear ordinary differential equations to model the three-dimensional flow of a nanofluid in a rotating channel on a lower permeable stretching wall. The resulting system of equations is solved using the differential transform method (DTM) [31, 32]. The DTM method has been applied successfully to solve nonlinear differential equations without requiring linearization or discretization [33]. The present DTM code is benchmarked with numerical results based on a shooting technique and previously published results. The DTM method shows excellent correlation with results obtained using these other methods.

The paper is divided up as follows: in Section 2, the problem is formulated and a coupled system of nonlinear ordinary differential equations is derived. The DTM method is applied to solve the resulting system of nonlinear ordinary differential equations in Section 3. Results are discussed in Section 4. Concluding remarks are made in Section 5.

2. Mathematical Formulation of the Problem

Consider unsteady 3D rotating nanofluid flow of an incompressible electrically conducting viscous fluid between two infinite horizontal plane walls. The lower plane is placed at and is stretched with a time-dependent velocity in x-direction. The upper plane is also placed at a variable distance and the fluid is squeezed with a time-dependent velocity in the negative y-direction. The fluid and the channel are rotating around y-axis with an angular velocity and also the lower plate intakes the flow with the velocity . A magnetic field with density is applied along the y-axis about the system which is rotating. These velocities and magnetic fields are introduced to obtain similarity solutions by reducing governing equations into ordinary differential equations (ODEs). The physical model of the considered problem along with the coordinate system is illustrated in Figure 1. The governing equations of continuity and momentum of nanofluid flow in a rotating frame of reference are given as [34, 35]

where is the Cauchy stress tensor, the magnetic flux, and the current density. The above governing equations can be also described by the following set of Navier-Stokes equations [36, 37]:

where is the nanofluid density, is the nanofluid kinematic viscosity, where has been proposed by Brinkman [38], is the electrical conductivity, is the magnetic field, and is the characteristic parameter with dimension of (time)⁻¹ and . The above nanofluid constants are defined as follows:

where is the viscosity of the fluid fraction, is the nanoparticle volume fraction, and and are the densities of the fluid and of the solid fractions, respectively. The thermophysical properties of the base fluid (water) and different nanoparticles are given in Table 1 [39]. The appropriate boundary conditions are introduced as follows:

where is the stretching rate of the lower plate. The following appropriate similarity transformations are employed to convert the above governing equations (2)–(4) into a system of ordinary differential equations in terms of a stream function :

Substituting the similarity transformations (5) and (6) into (2)–(4), we obtain the following system of nonlinear ordinary differential equations:

where is the characteristic parameter of the flow, is the rotation parameter, is the magnetic parameter, and prime denotes differentiation with respect to . In order to squeeze the flow, we take , for which the upper plate moves downward with velocity . For the upper plate moves apart with respect to the plane and corresponds to the steady state case of the considered problem or stationary upper plate. In order to reduce the number of independent variables and to retain the similarity solution, (7) are simplified by cross differentiation and thus we obtain the following system of differential equations:

The transformed boundary conditions take the form

where is the suction parameter.

Table 1: Thermophysical properties of the base fluid and different nanoparticles [39].

Physical properties
Fluid phase (water)
Cu
CuO
Ag
Al₂O₃
TiO₂

(kg/m³)
997.1
8933
6320
10500
3970
4250

Figure 1: Schematic diagram of the flow configuration and the coordinate system for the considered flow.

In this problem, the physical quantity of interest is the skin friction coefficients along the stretching wall at the lower and upper walls which are defined as [40]

Substituting (6) into (11) we obtain

where is the local Reynolds number.

In this section we have derived the coupled system of nonlinear ordinary differential equations (9) to model the three-dimensional flow of a nanofluid in a rotating channel on a lower permeable stretching wall. We next apply the DTM to solve the coupled system.

3. Analytical Approximation by Means of DTM

Taking the differential transform of (9), one can obtain (for more details, see [41, 42])

where and are differential transforms of and given by

where (16) is the transformed boundary conditions and , and are constants. By substituting (16) into (13), using recursion, and then substituting them in (14)-(15), we obtain the values of and given by

The constants , and can be determined by applying the remaining boundary conditions to (17). The number of required terms is determined by the convergence of the numerical values to one’s desired accuracy. We obtain the approximants using the computational software MATHEMATICA. The effect of different orders of the DTM solution on the convergence radius of the transverse velocity of is displayed in Figure 2.

Figure 2: The obtained results of transverse velocity of for different orders of DTM solution in comparison with the numerical solution when , , , , and .

In order to highlight the validity of the presented DTM solution, we compare some of our results with the obtained HAM results of [1]. A very good validation of the present analytical results has been achieved with the previously published study as shown in Table 2.

Table 2: Comparison results of shear stresses at lower and upper walls for different values of when , , and .

DTM Result
HAM result of [1]
DTM result
HAM result of [1]

0.0
+2.00373789
+2.0037379
−4.58270559
−4.5827056

0.3
+0.06127071
+0.0612707
−2.62278886
−2.6227888

0.6
−1.98233944
−1.9823394
−0.73129902
−0.7312990

0.9
−4.13494482
−4.1349448
+1.09241241
+1.0924124

1.2
−6.40476369
−6.4047636
+2.84916544
+2.8491655

4. Results and Discussion

In the previous section the nonlinear ordinary differential equations (9) are solved subject to the boundary conditions (10) via the DTM method for values of the five key parameters: characteristic parameter of the flow , rotation parameter , the magnetic parameter , nanoparticle volume fraction , the suction parameter , and different nanoparticles on the different velocity components. It should be stated that the copper nanoparticle is used in all of the cases in this section except for those figures which focus on the influences of the types of applied nanoparticles on the velocity component profiles. In addition, we assume that the values of the volume fraction parameter vary from 0 (regular Newtonian fluid) to 0.2. Representative values are employed to simulate the physically realistic flows. Tables 2 and 3 display the comparison between the DTM and numerical solution results based on a shooting technique for the shear stresses at lower and upper walls for different values of and .

Table 3: Comparison results of shear stresses at lower and upper walls or different values of and when , , and .

DTM result
Numerical result
DTM result
Numerical result

0.0
−3.0
−8.52098772
−8.52098769
+7.53714763
+7.53714759

−2.0
−6.76815369
−6.76815367
+5.55472891
+5.55472890

−1.0
−4.97050617
−4.97050616
+3.54801649
+3.54801648

0.0
−3.01634611
−3.01634612
+1.39835301
+1.39835302

+1.0
−0.89270759
−0.89270759
−0.90954349
−0.90954349

+2.0
+1.39749899
+1.39749899
−3.37418335
−3.37418335

+3.0
+3.84773079
+3.84773079
−5.99011133
−5.99011133

0.5
−3.0
−11.7761059
−11.7761055
+9.47810573
+9.47810571

−2.0
−9.51967384
−9.51967382
+7.18294804
+7.18294803

−1.0
−7.59674602
−7.59674601
+5.33619342
+5.33619342

0.0
−5.56301306
−5.56301307
+3.38177815
+3.38177815

+1.0
−3.37407009
−3.37407009
+1.26581626
+1.26581626

+2.0
−1.02490197
−1.02490197
−1.01698614
−1.01698614

+3.0
+1.48110397
+1.48110397
−3.46202219
−3.46202219

1.0
−3.0
−17.7188652
−17.7188648
+14.8771382
+14.8771378

−2.0
−12.7797205
−12.7797203
+8.72969245
+8.72969242

−1.0
−10.6266507
−10.6266506
+6.92986184
+6.92986182

0.0
−8.46845956
−8.46845957
+5.16113154
+5.16113153

+1.0
−6.18186918
−6.18186918
+3.24194096
+3.24194096

+2.0
−3.74763053
−3.74763053
+1.14862333
+1.14862333

+3.0
−1.16400368
−1.16400368
−1.11763550
−1.11763550

The influence of nanoparticle volume fraction on the normal, axial, and transverse velocity components is shown in Figures 3, 4, and 5. The normal velocity component reduces for higher values of the nanoparticle volume fraction. The axial velocity profile decreases with the increase in nanoparticle volume fraction in the lower half channel while in the upper half of the channel, enhances with the increase in . In addition, transverse velocity profile increases near the lower surface for the large value of nanoparticle volume fraction, but the adverse trend happens on the other side (in the center and the upper half of the channel).

Figure 3: Effect of on the velocity component of when and .

Figure 4: Effect of on the velocity component of when and .

Figure 5: Effect of on the velocity component of when and .

The effects of the characteristic parameter of flow on the velocity component are plotted in Figure 6. As one can see, with the increase in , the velocity profile augments in the vicinity of the upper plane. It can be easily understood that it accumulates the squeezing effects on the flow. In Figure 7, the effect of characteristic parameter of the flow on the velocity component, parallel to x-axis, , is depicted and it is clearly obvious that the velocity increases with the increase in the value of . It should be noted that the higher values of amortize the reverse flow, but the negative values of support the reverse flow because of the squeezing and extricating effects of the upper wall. The same trend can be observed for transverse velocity component. Figure 8 shows that the transverse velocity component increases as increases.

Figure 6: Effect of on the velocity component of when , , and .

Figure 7: Effect of on the velocity component of when , , and .

Figure 8: Effect of on the velocity component of when , , and .

Figure 9 shows the effect of rotation parameter on the normal velocity component . The normal velocity component reduces for large values of the rotation parameter. The effect of the rotation parameter on the axial velocity component is presented in Figure 10. The results show that in the lower half channel, the velocity component parallel to x-direction decreases with the increase in rotation parameter, while in the upper half of the channel, increases with the increase in . This reverse trend occurs approximately in the center of the channel . The effect of rotation parameter on the transverse velocity is shown in Figure 11. It is clear that although the large rotation causes an increase in transverse velocity near lower surface, the adverse trend occurs on the other side. In the main nanoflow regime, two diverse trends can be noticed. In the lower half channel the transverse velocity increases, in the center and the upper half of the channel the transverse velocity decreases.

Figure 9: Effect of on the transverse velocity of when , , and .

Figure 10: Effect of on the transverse velocity of when , , and .

Figure 11: Effect of on the transverse velocity of when , , and .

The effect of magnetic parameter on the normal, axial, and transverse velocity components can be observed in Figures 12–14. From Figure 12, it is seen that the velocity decreases as magnetic parameter increases because Lorentz force tends to decrease the velocity profile . In Figure 13, two different trends can be seen. In the lower half of the channel, the velocity component parallel to the x-direction decreases with the increase in the magnetic parameter, while in the upper half of the channel, increases with the increase in . This reverse trend occurs approximately in the center of the channel. The transverse velocity shows increasing behavior with the increase in magnetic parameter in all region of the channel (Figure 14).

Figure 12: Effect of on the velocity component of when , , and .

Figure 13: Effect of on the velocity component of when , , and .

Figure 14: Effect of on the velocity component of when , , and .

The consequence of suction parameter on all the velocity profiles is illustrated in Figures 15, 16, and 17. Figure 15 depicts that as the suction parameter increases, the normal velocity component increases and the variation in the velocity profile confines in the vicinity of the upper plate. It is known that the large suction values cause to decrease, which results in the occurrence of reverse flow. The reverse flow is more prominent near the upper plate rather than the lower plate. The reverse flow near the lower plate is produced by the adverse pressure gradient because the large amount of fluid particles escapes from the lower wall. Figure 17 shows the same behavior for the transverse velocity component .

Figure 15: Effect of on the velocity component of when , , and .

Figure 16: Effect of on the velocity component of when , , and .

Figure 17: Effect of on the velocity component of when , , and .

Figures 18, 19, and 20 demonstrate the effect of different nanoparticle materials on the normal velocity component , axial velocity component , and also the transverse velocity profile . As the results show, the maximum amount of normal velocity component belongs to nanoparticle. Moreover, the effects of nanoparticle volume fraction and the types of nanoparticles on the shear stresses at lower and upper walls are depicted in Table 4.

Table 4: Effect of nanoparticle volume fraction and nanoparticles material on the shear stresses at lower and upper walls when and .

Ag
Cu
CuO
TiO₂
Al₂O₃

0.00
−9.66151
+5.987272
−9.66151
+5.987272
−9.66151
+5.987272
−9.66151
+5.987272
−9.66151
+5.987272

0.05
−7.72463
+4.511947
−8.08637
+4.787926
−8.79156
+5.326106
−9.46511
+5.840246
−9.5657
+5.917038

0.10
−6.85916
+3.85395
−7.34313
+4.222859
−8.37683
+5.011777
−9.49215
+5.863511
−9.67086
+5.999997

0.15
−6.45632
+3.548727
−7.00493
+3.966591
−8.24921
+4.916332
−9.71233
+6.034163
−9.95938
+6.222928

0.20
−6.31569
+3.44302
−6.91686
+3.900868
−8.3391
+4.986844
−10.1216
+6.349316
−10.435
+6.588886

Figure 18: Effect of nanoparticle types on the velocity component of when , , and .

Figure 19: Effect of nanoparticle types on the velocity component of when , , and .

Figure 20: Effect of nanoparticle types on the velocity component of when , , and .

5. Conclusions

In this paper, we have used DTM to solve a coupled system of nonlinear ordinary differential equations for three-dimensional flow of a nanofluid in a rotating channel on a lower permeable stretching wall. We have considered water the base fluid and four different types of nanoparticles, copper, copper oxide, aluminum oxide, and titanium dioxide have been examined in this simulation. The upper wall is moving along the direction normal to the surface with time dependent velocity. The transformed dimensionless equations have been formulated and solved with robust boundary conditions. Important physical parameters have been investigated graphically. These parameters include the flow parameter, the rotation parameter, the magnetic parameter, the nanoparticle volume fraction, the suction parameter, and also different types of nanoparticles. Results show that the vertical motion of the upper plate interrupts the velocity in the channel remarkably. The downward motion of the upper plate augments the forward flow, whereas the upward motion reverses the flow.

Nomenclature

| Stretching rate of the lower plate
: | External uniform magnetic field
: | Constant magnetic flux density
: | Skin friction coefficient
: | Self-similar velocities
: | Upper plane distance
: | Magnetic flux
: | Pressure
: | Time
: | Cauchy stress tensor
: | Velocity component in the x direction
: | Velocity component in the y direction
: | Velocity component in the z direction.

Dimensionless Parameters

: | Characteristic parameter of the flow
: | Rotation parameter
: | Magnetic parameter
: | Suction parameter
: | Local Reynolds number .

Greek Symbols

: | DTM constant parameters
: | Electrical conductivity
: | Density
: | Dynamic viscosity
: | A scaled boundary-layer coordinate
: | Kinematic viscosity
: | Characteristic constant parameter
: | Stream function
: | Constant angular velocity
: | Nanoparticle volume fraction.

Subscripts

: | Fluid phase
: | Nanofluid
: | Solid phase.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

Ebrahim Momoniat acknowledges support from the National Research foundation of South Africa under Grant number 76868.

References

S. Munawar, A. Mehmood, and A. Ali, “Three-dimensional squeezing flow in a rotating channel of lower stretching porous wall,” Computers and Mathematics with Applications, vol. 64, no. 6, pp. 1575–1586, 2012.

B. C. Sakiadis, “Boundary layer behaviour on continuous solid surfaces: I Boundary layer equations for two dimensional and axisymmetric flow,” AIChE Journal, vol. 7, pp. 26–28, 1961.

B. C. Sakiadis, “Boundary layer behaviour on continuous solid surfaces: II boundary layer on a continuous flat surface,” AIChE Journal, vol. 7, pp. 221–225, 1961.

L. E. Erickson, L. T. Fan, and V. G. Fox, “Heat and mass transfer on a moving continuous flat plate with suction or injection,” Industrial and Engineering Chemistry Fundamentals, vol. 5, no. 1, pp. 19–25, 1966.

F. K. Tsou, E. M. Sparrow, and R. J. Goldstein, “Flow and heat transfer in the boundary layer on a continuous moving surface,” International Journal of Heat and Mass Transfer, vol. 10, no. 2, pp. 219–235, 1967.

L. J. Crane, “Flow past a stretching plate,” Zeitschrift für angewandte Mathematik und Physik ZAMP, vol. 21, no. 4, pp. 645–647, 1970.

H. I. Andersson, “MHD flow of a viscoelastic fluid past a stretching surface,” Acta Mechanica, vol. 95, no. 1–4, pp. 227–230, 1992.

K. V. Prasad, D. Pal, V. Umesh, and N. S. P. Rao, “The effect of variable viscosity on MHD viscoelastic fluid flow and heat transfer over a stretching sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 2, pp. 331–344, 2010.

A. Mehmood and A. Ali, “Analytic solution of generalized three-dimensional flow and heat transfer over a stretching plane wall,” International Communications in Heat and Mass Transfer, vol. 33, no. 10, pp. 1243–1252, 2006.

A. Mehmood and A. Ali, “Analytic homotopy solution of generalized three-dimensional channel flow due to uniform stretching of the plate,” Acta Mechanica Sinica, vol. 23, no. 5, pp. 503–510, 2007.

A. K. Borkakoti and A. Bharali, “Hydromagnetic flow and heat transfer between two horizontal plates, the lower plate being a stretching sheet,” Quarterly of Applied Mathematics, vol. 40, no. 4, pp. 461–467, 1983.

S. Munawar, A. Mehmood, and A. Ali, “Effects of slip on flow between two stretchable disks using optimal homotopy analysis method,” Canadian Journal of Applied Science, vol. 1, pp. 50–68, 2011.

A. J. Chamkha, T. Groşan, and I. Pop, “Fully developed free convection of a micropolar fluid in a vertical channel,” International Communications in Heat and Mass Transfer, vol. 29, no. 8, pp. 1119–1127, 2002.

R. Bhargava, L. Kumar, and H. S. Takhar, “Numerical solution of free convection MHD micropolar fluid flow between two parallel porous vertical plates,” International Journal of Engineering Science, vol. 41, no. 2, pp. 123–136, 2003.

P. D. Ariel, “On exact solutions of flow problems of a second grade fluid through two parallel porous walls,” International Journal of Engineering Science, vol. 40, no. 8, pp. 913–941, 2002.

T. Hayat and Z. Abbas, “Channel flow of a Maxwell fluid with chemical reaction,” Zeitschrift für angewandte Mathematik und Physik, vol. 59, no. 1, pp. 124–144, 2008.

G. Domairry and A. Aziz, “Approximate analysis of MHD dqueeze flow between two parallel disks with suction or injection by homotopy perturbation method,” Mathematical Problems in Engineering, vol. 2009, Article ID 603916, 19 pages, 2009.

S. U. S. Choi and J. A. Eastman, “Enhancing thermal conductivity of fluids with nanoparticles,” Materials Science, vol. 231, pp. 99–105, 1995.

Y. Xuan and Q. Li, “Investigation on convective heat transfer and flow features of nanofluids,” Journal of Heat Transfer, vol. 125, no. 1, pp. 151–155, 2003.

N. Bachok, A. Ishak, and I. Pop, “Flow and heat transfer over a rotating porous disk in a nanofluid,” Physica B: Condensed Matter, vol. 406, no. 9, pp. 1767–1772, 2011.

W. A. Khan and I. Pop, “Boundary-layer flow of a nanofluid past a stretching sheet,” International Journal of Heat and Mass Transfer, vol. 53, no. 11-12, pp. 2477–2483, 2010.

M. H. Abolbashari, N. Freidoonimehr, F. Nazari, and M. M. Rashidi, “Entropy analysis for an unsteady MHD flow past a stretching permeable surface in nano-fluid,” Powder Technology, vol. 267, pp. 256–267, 2014.

O. A. Beg, M. M. Rashidi, M. Akbari, and A. Hosseini, “Comparative numerical study of single-phase and two-phase models for bio-nanofluid transport phenomena,” Journal of Mechanics in Medicine and Biology, vol. 14, Article ID 1450011, 31 pages, 2014.

M. M. Rashidi, A. Hosseini, I. Pop, S. Kumar, and N. Freidoonimehr, “Comparative numerical study of single and two-phase models of nanofluid heat transfer in wavy channel,” Applied Mathematics and Mechanics, vol. 35, pp. 831–848, 2014.

E. Abu-Nada, Z. Masoud, H. F. Oztop, and A. Campo, “Effect of nanofluid variable properties on natural convection in enclosures,” International Journal of Thermal Sciences, vol. 49, no. 3, pp. 479–491, 2010.

M. M. Rashidi, S. Abelman, and N. Freidoonimehr, “Entropy generation in steady MHD flow due to a rotating porous disk in a nanofluid,” International Journal of Heat and Mass Transfer, vol. 62, no. 1, pp. 515–525, 2013.

M. Mustafa, T. Hayat, I. Pop, S. Asghar, and S. Obaidat, “Stagnation-point flow of a nanofluid towards a stretching sheet,” International Journal of Heat and Mass Transfer, vol. 54, no. 25-26, pp. 5588–5594, 2011.

M. M. Rashidi, N. Freidoonimehr, A. Hosseini, O. A. Bég, and T. K. Hung, “Homotopy simulation of nanofluid dynamics from a non-linearly stretching isothermal permeable sheet with transpiration,” Meccanica, vol. 49, pp. 469–482, 2014.

M. Sheikholeslami and D. D. Ganji, “Three dimensional heat and mass transfer in a rotating system using nanofluid,” Powder Technology, vol. 253, pp. 789–796, 2014.

M. Sheikholeslami, M. Gorji-Bandpy, and D. D. Ganji, “Lattice Boltzmann method for MHD natural convection heat transfer using nanofluid,” Powder Technology, vol. 254, pp. 82–93, 2014.

F. Ayaz, “Applications of differential transform method to differential-algebraic equations,” Applied Mathematics and Computation, vol. 152, no. 3, pp. 649–657, 2004.

F. Ayaz, “Solutions of the system of differential equations by differential transform method,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 547–567, 2004.

M. M. Rashidi and E. Erfani, “A new analytical study of MHD stagnation-point flow in porous media with heat transfer,” Computers and Fluids, vol. 40, no. 1, pp. 172–178, 2011.

M. Kurosaka, “The oscillatory boundary layer growth over the top and bottom plates of a rota ting channel,” Journal of Fluids Engineering, Transactions of the ASME, vol. 95, no. 1, pp. 68–74, 1973.

H. S. Takhar, A. J. Chamkha, and G. Nath, “MHD flow over a moving plate in a rotating fluid with magnetic field, hall currents and free stream velocity,” International Journal of Engineering Science, vol. 40, no. 13, pp. 1511–1527, 2002.

K. Vajravelu and B. V. R. Kumar, “Analytical and numerical solutions of a coupled non-linear system arising in a three-dimensional rotating flow,” International Journal of Non-Linear Mechanics, vol. 39, no. 1, pp. 13–24, 2004.

G. Domairry and A. Aziz, “Approximate analysis of MHD squeeze flow between two parallel disks with suction or injection by homotopy perturbation method,” Mathematical Problems in Engineering, vol. 2009, Article ID 603916, 19 pages, 2009.

H. C. Brinkman, “The viscosity of concentrated suspensions and solutions,” The Journal of Chemical Physics, vol. 20, no. 4, p. 571, 1952.

H. F. Oztop and E. Abu-Nada, “Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids,” International Journal of Heat and Fluid Flow, vol. 29, no. 5, pp. 1326–1336, 2008.

M. Sheikholeslami and D. D. Ganji, “Heat transfer of Cu-water nanofluid flow between parallel plates,” Powder Technology, vol. 235, pp. 873–879, 2013.

M. M. Rashidi, “The modified differential transform method for solving MHD boundary-layer equations,” Computer Physics Communications, vol. 180, no. 11, pp. 2210–2217, 2009.

M. M. Rashidi and N. Freidoonimehr, “Series solutions for the flow in the vicinity of the equator of an MHD boundary-layer over a porous rotating sphere with heat transfer,” Thermal Science, vol. 18, pp. S527–S537, 2012.