Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 861065, 12 pages
http://dx.doi.org/10.1155/2015/861065

Research Article

Heat and Mass Transfer for MHD Viscoelastic Fluid Flow over a Vertical Stretching Sheet with Considering Soret and Dufour Effects

Mohammad Mehdi Rashidi,1,2,3 Mohamed Ali,⁴ Behnam Rostami,³ Peyman Rostami,⁵ and Gong-Nan Xie⁶

¹Shanghai Automotive Wind Tunnel Center, Tongji University, 4800 Caoan Road, Jiading, Shanghai 201804, China
²ENN-Tongji Clean Energy Institute of Advanced Studies, Shanghai, China
³Mechanical Engineering Department, Engineering Faculty, Bu-Ali Sina University, Hamedan, Iran
⁴Mechanical Engineering Department, College of Engineering, King Saud University, P.O. Box 800, Riyadh 11421, Saudi Arabia
⁵Department of Mechanical Engineering, Isfahan University of Technology, Isfahan 84156 83111, Iran
⁶School of Mechanical Engineering, Northwestern Polytechnical University, Xi’an, Shaanxi, China

Received 22 March 2014; Revised 17 June 2014; Accepted 6 September 2014

Academic Editor: Haochun Zhang

Copyright © 2015 Mohammad Mehdi Rashidi 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.

The homotopy analysis method (HAM) with two auxiliary parameters is employed to examine heat and mass transfer in a steady two-dimensional magneto hydrodynamic viscoelastic fluid flow over a stretching vertical surface by considering Soret and Dufour effects. The two-dimensional boundary-layer governing partial differential equations are derived by considering the Boussinesq approximation. The highly nonlinear ordinary differential forms of momentum, energy, and concentration equations are obtained by similarity transformation. These equations are solved analytically in the presence of buoyancy force. The effects of different involved parameters such as magnetic field parameter, Prandtl number, buoyancy parameter, Soret number, Dufour number, and Lewis number on velocity, temperature, and concentration profiles are plotted and discussed. The effect of the second auxiliary parameter is also illustrated. Results show that the effect of increasing Soret number or decreasing Dufour number tends to decrease the velocity and temperature profiles (increase in Sr cools the fluid and reduces the temperature) while enhancing the concentration distribution.

1. Introduction

The analysis of the flow field in a boundary-layer near a stretching sheet is an important part in fluid dynamics and heat transfer. This type of flow occurs in a number of engineering processes such as extrusion of plastic sheets, polymer processing, and metallurgy [1, 2].

Some researchers neglect the Dufour and Soret effects on heat and mass transfer according to Fourier’s and Fick’s laws [3]; however, when density differences exist in the flow regime, these effects are important and cannot be neglected [4]. Afify [5] has shown that when heat and mass transfer occurred in a moving fluid, the energy flux can be generated by a composition gradient, namely, the Dufour or diffusion-thermo effect, and the mass fluxes developed by the temperature gradients are called the Soret or thermal-diffusion effect. In their numerical study they have used the Soret and Dufour effects of a steady flow due to a rotating disk in the presence of viscous dissipation and ohmic heating. Heat and mass transfer with hydrodynamic slip over a moving plate in porous media was investigated by Hamad et al. [6] via Runge-Kutta-Fehlberg fourth-fifth order method. The heat transfer of mixed convection of vertically moving surface in an ambient stagnant fluid was reported by Ali and Al-Yousef [7, 8] and the effect of variable viscosity of mixed convection was studied by Ali [9].

Das et al. [10] considered the effect of heat and mass transfer on a free convective flow of an incompressible electrically conducting fluid past a vertical porous plate. Chen [11] employed finite difference method in order to study the heat and mass transfer in MHD free convective flow with ohmic heating and viscous dissipation. Noor et al. [12] examined the MHD flow over an inclined surface with heat source/sink effects by shooting method. Abreu et al. [13] solved the boundary-layer flow with Dufour and Soret effects in both forced and natural convection. The effects of thermal radiation and first order chemical reaction on unsteady MHD convective flow past a semiinfinite vertical plate under oscillatory suction and heat source in slip-flow regime were taken into account by Pal and Talukdar [14]. Gbadeyan et al. [15] studied heat and mass transfer of a mixed convection boundary-layer flow considering porous medium over a stretching vertical surface. A vertical plate in a non-Darcy porous medium was selected to investigate the thermodiffusion and diffusion-thermo effects numerically using the Keller-box method by Prasad et al. [16]. Pal and Mondal [17–19] analyzed the effects of thermal diffusion and diffusion thermo on steady and unsteady MHD non-Darcy flow over a stretching sheet in a porous medium considering thermal radiation, nonuniform heat source/sink, variable viscosity, viscous dissipation, and first order chemical reaction using Runge-Kutta-Fehlberg integration method. Mansour et al. [20] analyzed the effects of chemical reaction and thermal stratification over a vertical stretching surface in a porous medium by Runge-Kutta scheme with considering Soret and Dufour numbers. Bég et al. [21] used Keller-box implicit method to analyze the heat and mass transfer micropolar fluid flow from an isothermal sphere with Soret and Dufour effects. Furthermore, Alam et al. [22], Tai and Char [23], Mahdy [24, 25], Pal and Sewli [26], and Tsai and Huang [27] have studied the effect of Soret and Dufour effects in their analyses for different aspects of heat and mass transfer flows.

One of the most effective and reliable methods in order to solve the high nonlinear problems is the homotopy analysis method. Homotopy analysis method (HAM) was firstly employed by Liao to offer a general analytic method for nonlinear problems [28, 29]. Rashidi et al. [30] analyzed the effect of partial slip, diffusion thermo, and thermal diffusion on MHD fluid flow in a rotating disk via HAM and discussed the effect of various slip parameters, magnetic field parameter, Schmidt number, and other important variables. Mustafa et al. [31] considered the effects of Brownian motion and thermophoresis in stagnation point flow of a nanofluid towards a stretching sheet. Rashidi and Pour [32] employed HAM for unsteady boundary-layer flow and heat transfer on a stretching sheet. Abbas et al. [33] studied the mixed convection of an incompressible Maxwell fluid flow over a vertical stretching surface by HAM. Dinarvand et al. [34] employed HAM to investigate unsteady laminar MHD flow near forward stagnation point of a rotating and translating sphere. Hayat et al. [35] illustrated the thermal-diffusion and diffusion-thermo effects on two-dimensional MHD axisymmetric flow of a second grade fluid in the presence of Joule heating and first order chemical reaction. The non-linear Brinkman equation for the stagnation-point flow was studied via HAM by Ziabakhsh et al. [36]. Analytical and numerical solutions of a radial stagnation flow over a stretching cylinder have been recently reported by Weidman and Ali [37] where aligned and nonaligned flow were studied. Rashidi et al. [38, 39] employed HAM to obtain the analytical solutions over stretching and shrinking sheets in the presence of buoyancy parameter.

The objective of this paper is to study the steady two-dimensional MHD viscoelastic fluid flow over a vertical stretching surface in the presence of the Soret and Dufour effects analytically via HAM. The effects of different involved parameters such as magnetic field parameter, Prandtl number, buoyancy parameter, Soret number, Dufour number, and Lewis number on the fluid velocity, temperature, and concentration distributions are plotted and discussed.

2. Flow Analysis

Consider a steady two-dimensional heat and mass transfer flow of an incompressible electrically conducting viscoelastic fluid over a stretching vertical surface with a variable magnetic field normally applied to the surface. Keeping the origin fixed, two equal and opposite forces are applied along the -axis. It is assumed that the stretching velocity is in the form of , where and are constants. The induced magnetic field is neglected in comparison to the applied magnetic field and the viscous dissipation is small. The governing equations subject to Boussinesq approximation, the boundary-layer assumptions, and the above assumptions can be written as (for more details see [41])

where and are velocity components in the directions of and along and normal to the surface, respectively (as shown in Figure 1). is the kinematic viscosity, is the viscoelasticity parameter, is the electrical conductivity, is the fluid density, is the acceleration due to gravity, is the coefficient of thermal expansion, is the coefficient of thermal expansion with concentration, is the thermal diffusivity, is the thermal diffusion ratio, is the concentration susceptibility, is the specific heat at constant pressure, is the coefficient of mass diffusivity, is the fluid temperature, is the fluid concentration, and is the mean fluid temperature.

Figure 1: The schematic diagram of the stretching sheet problem.

The corresponding boundary conditions are as follows:

We assume that and , where and are constants. Introducing stream function , and similarity variable , [42] the continuity equation is satisfied and the momentum, energy, and concentration equations are transformed into ordinary differential equations as follows:

where superscript denotes the derivative with respect to , is the viscoelasticity parameter (when the viscoelastic parameter takes the form of similar to the viscoelastic parameter obtained by Hayat et al. [41]), is the magnetic field parameter, is the buoyancy parameter, where is the Grashof number, is the Reynolds number, is the constant dimensionless concentration buoyancy parameter, is the Prandtl number, is the Lewis number, is the Soret number, and is the Dufour number. The corresponding boundary conditions are as follows:

where is the suction/injection parameter ( for suction and for injection). In this paper the suction parameter has been considered because the primary assumption in boundary-layer definition says that the boundary-layer thickness is supposed to be very thin and we are not allowed to increase it, so we do not present the injection parameters that may lead to enlarging the boundary-layer thickness and contravening the boundary-layer assumption presented by Prandtl in 1904.

3. HAM Solution

We choose the initial approximations to satisfy the boundary conditions. The appropriate initial approximations are as follows:

where is the second auxiliary parameter. The linear operators , and are

with the following properties:

where are arbitrary constants and the nonlinear operators are

The auxiliary functions are introduced as

The th order deformation equations (see (10)) can be solved by the symbolic software MATHEMATICA

where is the auxiliary nonzero parameter

For more information about the HAM solution, see [28, 29].

In Figure 2,  -curve is figured, obtained via 20th order of HAM solution. The averaged residual errors are defined as (12) to acquire optimal values of auxiliary parameters:

Figure 2: The -curves of , , and obtained by the 20th order approximation of the HAM solution when , , , , , , , , , and .

To check the accuracy of the method, the residual errors of (12) are illustrated in Figures 3 and 4. The residual errors are reduced when we use the second auxiliary parameter and this justifies why we use the second auxiliary parameter. In Figure 3, the effect of considering is to decrease the order of residual errors than at (without the second auxiliary parameter) in Figure 4 which improves the accuracy of the HAM method. The velocity profiles presented in Figure 5 show an excellent agreement between our results and [40].

Figure 3: The residual errors when , , , , , , , , , and .

Figure 4: The residual errors when , , , , , , , , , and .

Figure 5: Verification of obtained by the 20th order of HAM solution with pervious published paper [40] when , , , , and .

4. Results and Discussion

In this paper the MHD two-dimensional steady heat and mass transfer flow of an incompressible viscoelastic fluid over a stretching vertical surface with considering the effects of Soret and Dufour numbers is investigated. Applying numerical values to the problem parameters, we can discuss their effects on the velocity , temperature , and concentration distributions. Graphical illustration of the results is very useful and practical to discuss the effect of different parameters. In this analysis, it has been considered that [41]. Negative (thermal and concentration buoyancy forces oppose each other) induces a slight increase in both fluid temperature and concentration [43]. In this paper the value of is considered to be 0.5. The effect of magnetic parameter on the velocity is plotted in Figure 6. Transverse magnetic field parameter creates a drag force, namely, Lorentz force that resists the flow and slows down the flow and causes to decrease the velocity. In Figure 7 the effect of magnetic field parameter on temperature profiles is illustrated. Magnetic field parameter causes skin-frictional heating and so the wall temperature increases and the thickness of thermal boundary-layer increases. The effect of is to increase the concentration profile (Figure 8). The governing equations are coupled together only with the buoyancy parameters. When increases, the Grashof number accelerates the fluid so the velocity and the boundary-layer thickness increases, as shown in Figure 9. The effect of on temperature and concentration profiles is shown in Figures 10 and 11. Both the thermal and concentration boundary-layer thicknesses decrease with the increase in the value of buoyancy parameter. The effects of Prandtl number on velocity, temperature, and concentration distributions are illustrated in Figures 12–14, respectively. Increase in leads to increase in kinematic viscosity and velocity decreases. It is clearly shown that with the increase in the velocity profiles descends (Figure 12). With the increase in Prandtl number the thermal diffusion decreases, so the thermal boundary-layer becomes thinner and temperature decreases. A fluid with larger and higher heat capacity increases the heat transfer [40] (Figure 13). The Pr reduces the concentration distribution just the same as its effect on temperature profile (Figure 14). The Soret effect is a mass flux due to a temperature gradient and the Dufour effect is enthalpy flux due to a concentration gradient and appears in the energy equation. The effects of Soret and Dufour numbers on velocity, temperature, and concentration profiles are plotted in Figures 15, 16, and 17, respectively. We considered the effects of and so that their product remains constant at 0.05. As one can see the increase in the value of or decrease in descends the velocity and temperature profiles and enhances the concentration distribution. Increase in Soret number cools the fluid and reduces the temperature [43]. Lewis number is the ratio of thermal diffusivity to mass diffusivity. The Lewis number can also be expressed as the ratio of the Schmidt number to the Prandtl number (), where is the Schmidt number. Figure 18 displays the effect of Lewis number on the velocity profile. The effect of increasing the value of on the velocity is as the same as the effect of decreasing the value of and it can be easily understood that with the enhancement of the velocity distribution increases. The effect of on temperature profile is presented in Figure 19. The temperature decreases with the increase in Lewis number similar to the results presented by Hayat et al. [41]. With the increase in the mass diffusivity decreases and the concentration descends (Figure 20). It should be noticed that permits complete similarity solutions of the equations, where and are constants and not . However in this problem must be constant and is selected equal to 0.5 in order to reach the local similarity solution.

Figure 6: The effect of on velocity profile when , , , , , , , and .

Figure 7: The effect of on temperature profile when , , , , , , , and .

Figure 8: The effect of on concentration profile when , , , , , , , and .

Figure 9: The effect of on velocity profile when , , , , , , , and .

Figure 10: The effect of on temperature profile when , , , , , , , and .

Figure 11: The effect of on concentration profile when , , , , , , , and .

Figure 12: The effect of on velocity profile when , , , , , , , and .

Figure 13: The effect of on temperature profile when , , , , , , , and .

Figure 14: The effect of on concentration profile when , , , , , , , and .

Figure 15: The effect of , on velocity profile when , , , , , , and .

Figure 16: The effect of , on temperature profile when , , , , , , and .

Figure 17: The effect of , on concentration profile when , , , , , , and .

Figure 18: The effect of on velocity profile when , , , , , , , and .

Figure 19: The effect of on temperature profile when , , , , , , , and .

Figure 20: The effect of on concentration profile when , , , , , , , and .

5. Conclusion

In the present investigation, an analysis is carried out in order to study the steady magneto hydrodynamic incompressible viscoelastic fluid flow over a stretching surface in the presence of the Soret and Dufour effects analytically via HAM with two auxiliary parameters. Analytical solutions are obtained using the homotopy analysis method and its residual was reduced by using the second auxiliary parameter. These analytical solutions show excellent agreement with the data available in the literature (Figures 3–5). The effect of Mn is to decrease the velocity while increasing the thermal boundary-layer. The effect of increasing the buoyancy parameter is to reduce both the thermal and concentration boundary-layer thicknesses. The effect of increasing Sr or decreasing Du tends to decreases the velocity and temperature profiles while enhancing the concentration distribution. The temperature profiles are not sensitive to increasing Le however the concentration profiles are very sensitive.

Nomenclature

: | Constant values [–]
: | Magnetic field [kg s⁻² A⁻¹]
: | Arbitrary constant [–]
: | Concentration [kg m⁻²]
: | Specific heat at constant pressure [J kg⁻¹ K⁻¹]
: | Coefficient of mass diffusivity [m² s⁻¹]
: | Dufour number ( [–]
: | Auxiliary nonzero parameter
: | Auxiliary function
: | Auxiliary linear operator
: | Lewis number [–]
: | Magnetic field parameter [–]
: | Nonlinear operator
: | Constant dimensionless concentration buoyancy parameter
: | Prandtl number [–]
: | Reynolds number [–]
: | Soret number () [–]
: | Mean fluid temperature [K].

Greek Letters

: | Thermal diffusivity [m² s⁻¹]
: | Coefficient of thermal expansion [K⁻¹]
: | Coefficient of thermal expansion with concentration [kg⁻¹ m³]
: | Dimensionless fluid concentration ( [–]
: | The second auxiliary parameter
: | Similarity variable ( [–]
: | Dimensionless fluid temperature ( [–]
: | Density [kg m⁻³]
: | Fluid electrical conductivity [S m⁻¹]
: | Buoyancy parameter () [–]
: | Fluid kinematic viscosity [m² s⁻¹]
: | Stream function.

Subscripts

: | Wall condition
: | Infinity condition.

Superscript

′:
Differentiation with respect to .

Conflict of Interests

All the authors have no conflict of interests to report.

Acknowledgments

The authors express their gratitude to the anonymous referees for their constructive reviews of the paper and for helpful comments. The authors extend their appreciation to the Deanship of Scientific Research at King Saud University for funding this work through the research group Project no. RGP-VPP-080.

References

S. Abel, K. V. Prasad, and A. Mahaboob, “Buoyancy force and thermal radiation effects in MHD boundary layer visco-elastic fluid flow over continuously moving stretching surface,” International Journal of Thermal Sciences, vol. 44, no. 5, pp. 465–476, 2005.

R. Tamizharasi and V. Kumaran, “Pressure in MHD/Brinkman flow past a stretching sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 12, pp. 4671–4681, 2011.

T. Hayat and F. A. Hendi, “Thermal-diffusion and diffusion-thermo effects on MHD three-dimensional axisymmetric flow with Hall andion-slip currents,” Journal of American Science, vol. 8, pp. 284–294, 2012.

S. P. A. Devi and R. U. Devi, “Soret and Dufour effects on MHD slip flow with thermal radiation over a porous rotating infinite disk,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 4, pp. 1917–1930, 2011.

A. A. Afify, “Similarity solution in MHD: effects of thermal diffusion and diffusion thermo on free convective heat and mass transfer over a stretching surface considering suction or injection,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 2202–2214, 2009.

M. A. A. Hamad, M. J. Uddin, and A. I. M. Ismail, “Investigation of combined heat and mass transfer by Lie group analysis with variable diffusivity taking into account hydrodynamic slip and thermal convective boundary conditions,” International Journal of Heat and Mass Transfer, vol. 55, no. 4, pp. 1355–1362, 2012.

M. Ali and F. Al-Yousef, “Laminar mixed convection from a continuously moving vertical surface with suction or injection,” Heat and Mass Transfer, vol. 33, no. 4, pp. 301–306, 1998.

M. Ali and F. Al-Yousef, “Laminar mixed convection boundary layers induced by a linearly stretching permeable surface,” International Journal of Heat and Mass Transfer, vol. 45, no. 21, pp. 4241–4250, 2002.

M. E. Ali, “The effect of variable viscosity on mixed convection heat transfer along a vertical moving surface,” International Journal of Thermal Sciences, vol. 45, no. 1, pp. 60–69, 2006.

S. S. Das, A. Satapathy, J. K. Das, and J. P. Panda, “Mass transfer effects on MHD flow and heat transfer past a vertical porous plate through a porous medium under oscillatory suction and heat source,” International Journal of Heat and Mass Transfer, vol. 52, no. 25-26, pp. 5962–5969, 2009.

C.-H. Chen, “Combined heat and mass transfer in MHD free convection from a vertical surface with Ohmic heating and viscous dissipation,” International Journal of Engineering Science, vol. 42, no. 7, pp. 699–713, 2004.

N. F. M. Noor, S. Abbasbandy, and I. Hashim, “Heat and mass transfer of thermophoretic MHD flow over an inclined radiate isothermal permeable surface in the presence of heat source/sink,” International Journal of Heat and Mass Transfer, vol. 55, no. 7-8, pp. 2122–2128, 2012.

C. R. A. Abreu, M. F. Alfradique, and A. S. Telles, “Boundary layer flows with dufour and soret effects: I: forced and natural convection,” Chemical Engineering Science, vol. 61, no. 13, pp. 4282–4289, 2006.

D. Pal and B. Talukdar, “Influence of fluctuating thermal and mass diffusion on unsteady MHD buoyancy-driven convection past a vertical surface with chemical reaction and Soret effects,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 4, pp. 1597–1614, 2012.

J. A. Gbadeyan, A. S. Idowu, A. W. Ogunsola, O. O. Agboola, and P. O. Olanrewaju, “Heat and mass transfer for Soret and Dufours effect on mixed convection boundary layer flow over a stretching vertical surface in a porous medium filled with a viscoelastic fluid in the presence of magnetic field,” Global Journal of Science Frontier Research, vol. 11, pp. 97–114, 2011.

V. R. Prasad, B. Vasu, O. A. B{\'e}g, and R. D. Parshad, “Thermal radiation effects on magnetohydrodynamic free convection heat and mass transfer from a sphere in a variable porosity regime,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 2, pp. 654–671, 2012.

D. Pal and H. Mondal, “MHD non-Darcian mixed convection heat and mass transfer over a non-linear stretching sheet with Soret-Dufour effects and chemical reaction,” International Communications in Heat and Mass Transfer, vol. 38, no. 4, pp. 463–467, 2011.

D. Pal and H. Mondal, “Effects of Soret Dufour, chemical reaction and thermal radiation on MHD non-Darcy unsteady mixed convective heat and mass transfer over a stretching sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 4, pp. 1942–1958, 2011.

D. Pal and H. Mondal, “MHD non-Darcy mixed convective diffusion of species over a stretching sheet embedded in a porous medium with non-uniform heat source/sink, variable viscosity and Soret effect,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 2, pp. 672–684, 2012.

M. A. Mansour, N. F. El-Anssary, and A. M. Aly, “Effects of chemical reaction and thermal stratification on MHD free convective heat and mass transfer over a vertical stretching surface embedded in a porous media considering Soret and Dufour numbers,” Chemical Engineering Journal, vol. 145, no. 2, pp. 340–345, 2008.

O. A. Bég, V. R. Prasad, B. Vasu, N. B. Reddy, Q. Li, and R. Bhargava, “Free convection heat and mass transfer from an isothermal sphere to a micropolar regime with Soret/Dufour effects,” International Journal of Heat and Mass Transfer, vol. 54, no. 1–3, pp. 9–18, 2011.

M. S. Alam, M. M. Rahman, and M. A. Sattar, “Effects of variable suction and thermophoresis on steady MHD combined free-forced convective heat and mass transfer flow over a semi-infinite permeable inclined plate in the presence of thermal radiation,” International Journal of Thermal Sciences, vol. 47, no. 6, pp. 758–765, 2008.

B.-C. Tai and M.-I. Char, “Soret and Dufour effects on free convection flow of non-Newtonian fluids along a vertical plate embedded in a porous medium with thermal radiation,” International Communications in Heat and Mass Transfer, vol. 37, no. 5, pp. 480–483, 2010.

A. Mahdy, “MHD non-Darcian free convection from a vertical wavy surface embedded in porous media in the presence of Soret and Dufour effect,” International Communications in Heat and Mass Transfer, vol. 36, no. 10, pp. 1067–1074, 2009.

A. Mahdy, “Soret and Dufour effect on double diffusion mixed convection from a vertical surface in a porous medium saturated with a non-Newtonian fluid,” Journal of Non-Newtonian Fluid Mechanics, vol. 165, no. 11-12, pp. 568–575, 2010.

D. Pal and S. C. Sewli, “Mixed convection magnetohydrodynamic heat and mass transfer past a stretching surface in a micropolar fluid-saturated porous medium under the influence of Ohmic heating, Soret and Dufour effects,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 3, pp. 1329–1346, 2011.

R. Tsai and J. S. Huang, “Numerical study of Soret and Dufour effects on heat and mass transfer from natural convection flow over a vertical porous medium with variable wall heat fluxes,” Computational Materials Science, vol. 47, no. 1, pp. 23–30, 2009.

S. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, Chapman & Hall/CRC Press, 2004.

S. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004.

M. M. Rashidi, T. Hayat, E. Erfani, S. A. M. Pour, and A. A. Hendi, “Simultaneous effects of partial slip and thermal-diffusion and diffusion-thermo on steady MHD convective flow due to a rotating disk,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 11, pp. 4303–4317, 2011.

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 and S. A. M. Pour, “Analytic approximate solutions for unsteady boundary-layer flow and heat transfer due to a stretching sheet by homotopy analysis method,” Nonlinear Analysis: Modelling and Control, vol. 15, no. 1, pp. 83–95, 2010.

Z. Abbas, Y. Wang, T. Hayat, and M. Oberlack, “Mixed convection in the stagnation-point flow of a Maxwell fluid towards a vertical stretching surface,” Nonlinear Analysis: Real World Applications, vol. 11, no. 4, pp. 3218–3228, 2010.

S. Dinarvand, A. Doosthoseini, E. Doosthoseini, and M. M. Rashidi, “Series solutions for unsteady laminar MHD flow near forward stagnation point of an impulsively rotating and translating sphere in presence of buoyancy forces,” Nonlinear Analysis: Real World Applications, vol. 11, no. 2, pp. 1159–1169, 2010.

T. Hayat, M. Nawaz, S. Asghar, and S. Mesloub, “Thermal-diffusion and diffusion-thermo effects on axisymmetric flow of a second grade fluid,” International Journal of Heat and Mass Transfer, vol. 54, no. 13-14, pp. 3031–3041, 2011.

Z. Ziabakhsh, G. Domairry, and H. R. Ghazizadeh, “Analytical solution of the stagnation-point flow in a porous medium by using the homotopy analysis method,” Journal of the Taiwan Institute of Chemical Engineers, vol. 40, no. 1, pp. 91–97, 2009.

P. D. Weidman and M. E. Ali, “Aligned and nonaligned radial stagnation flow on a stretching cylinder,” European Journal of Mechanics, B/Fluids, vol. 30, no. 1, pp. 120–128, 2011.

M. M. Rashidi, M. Ashraf, B. Rostami, M. T. Rastegari, and S. Bashir, “Mixed convection boundary-layer flow of a micro polar fluid towards a heated shrinking sheet by homotopy analysis method,” Thermal Science, 2013.

M. M. Rashidi, E. Momoniat, and B. Rostami, “Analytic approximate solutions for MHD boundary-layer viscoelastic fluid flow over continuously moving stretching surface by homotopy analysis method with two auxiliary parameters,” Journal of Applied Mathematics, vol. 2012, Article ID 780415, 19 pages, 2012.

C.-H. Chen, “On the analytic solution of MHD flow and heat transfer for two types of viscoelastic fluid over a stretching sheet with energy dissipation, internal heat source and thermal radiation,” International Journal of Heat and Mass Transfer, vol. 53, no. 19-20, pp. 4264–4273, 2010.

T. Hayat, M. Mustafa, and I. Pop, “Heat and mass transfer for Soret and Dufour's effect on mixed convection boundary layer flow over a stretching vertical surface in a porous medium filled with a viscoelastic fluid,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 5, pp. 1183–1196, 2010.

A. Ishak, R. Nazar, and I. Pop, “MHD mixed convection boundary layer flow towards a stretching vertical surface with constant wall temperature,” International Journal of Heat and Mass Transfer, vol. 53, no. 23-24, pp. 5330–5334, 2010.

O. A. Bég, A. Bakier, R. Prasad, and S. K. Ghosh, “Numerical modelling of non-similar mixed convection heat and species transfer along an inclined solar energy collector surface with cross diffusion effects,” World Journal of Mechanics, vol. 1, pp. 185–196, 2011.