Advances in Mathematical Physics
Volume 2016 (2016), Article ID 7036728, 14 pages
http://dx.doi.org/10.1155/2016/7036728

Research Article

Homotopy Analysis of the Radiation Effect on MHD Flow with Heat and Mass Transfer due to a Point Sink

C. N. Guled and B. B. Singh

Department of Mathematics, Dr. Babasaheb Ambedkar Technological University, Lonere, Raigad District 402103, India

Received 4 June 2016; Revised 8 September 2016; Accepted 27 September 2016

Academic Editor: Alkesh Punjabi

Copyright © 2016 C. N. Guled and B. B. Singh. 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.

An analytical solution of the magnetohydrodynamic, steady, and incompressible laminar boundary layer flow in the presence of heat and mass transfer as well as magnetic field on a cone due to a point sink by using the homotopy analysis method (HAM) has been studied under the radiative fluid properties. The HAM produces an analytical solution of the governing self-similar nonlinear two-point boundary layer equations. The effects of the suction/injection, magnetic, and radiation parameters over the obtained solution have been discussed. The effects of Prandtl number on temperature and Schmidt number on concentration profiles have also been studied. It has been observed that the temperature profiles exhibit an increasing trend with radiation in case of injection while an opposite trend is observed in case of suction. The results obtained in the present study have also been compared numerically as well as graphically with the corresponding results obtained by using other methods. An excellent agreement has been found between them. The analytical solution obtained by the HAM is very near to the exact solution for a properly selected initial guess, auxiliary, and convergence control parameters and for higher orders of deformations.

1. Introduction

In many areas of applied engineering and industry, we often come across the boundary layer flow of a steady incompressible laminar fluid flow in the presence of mass transfer and applied magnetic field. This kind of study is of interest for the fields of different branches of innovation, for example, in vortex chambers, magnetohydrodynamic (MHD) power generators, atomic reactors, and geophysical liquid flow. The examination of the boundary layer flow of an electrically conducting fluid on a cone because of a point sink with an applied magnetic field is significant in the investigation of conical nozzle or diffuser-flow problems and it was first of all concentrated on by Choi and Wilhelm [1]. Prior to this problem, Rosenhead [2] studied the same problem in the absence of magnetic field, mass flux diffusion, and heat transfer. Ackerberg [3] presented the series solution for the converging motion of the viscous fluid inside a cone. Takhar et al. [4] extended the same problem for electrically conducting fluid and discussed the heat and mass transfer effects. Eswara et al. [5] investigated the problem for the transient case. Eswara and Bommaiah [6] revisited the problem by taking into account temperature dependent viscosity. Turkyilmazoglu [7] considered the Falkner-Skan flows past stretching boundaries when the momentum and thermal slip boundary conditions are allowed at the boundary. Turkyilmazoglu extended the flow model set-up in a moving convergent channel by Magyari [8] by taking into consideration the momentum slip condition at the wall and found exact analytical solutions for the converging channel, for example, wedge nozzle.

In the context of space technology and processes involving temperatures, the effects of radiation are of vital significance. Recent developments in hypersonic flights, missile reentry, rocket combustion chambers, power plants for interplanetary flight, and gas cooled thermal reactors have focused their attention on thermal radiation as a mode of energy transfer. As a consequence of this, Vyas and Rai [9] made an elaborate analysis of the radiative flow inside a circular cone due to a point sink at the vertex of the cone. It is here worth mentioning that, unlike convection/conduction, the radiative heat transfer mechanism is rather more complex. However, some reasonable approximations have been found satisfactory to make the radiative systems solvable. The works of Sparrow and Cess [10] and Howell [11] describe the essentials of the radiative heat transfer. Many other pertinent radiative heat transfer studies for different configurations have also been reported by authors like Plumb et al. [12], Hossain and Takhar [13], Raptis [14], Sedeek and Salem [15], Al-Odat et al. [16], Prasad et al. [17], Mukhopadhyay [18], Vyas and Srivastava [19], Vyas and Ranjan [20], Chauhan and Kumar [21], Baoku et al. [22], Babu et al. [23], and so forth.

The objective of the present paper is to extend the work of Takhar et al. [4] by taking into account the radiative properties of the fluid at the wall. The flow problem which is governed by nonlinear equations with two-point boundary conditions has been solved by using HAM. Using the recursive method derived by Liao [24, 25], rigorous recursive formulae have been developed. Symbolic computation software and high performance computers have been used to derive the analytic solutions. The flow characteristics have been analyzed, and the results have been compared with those of [4] by setting the radiation parameter as zero.

In 1992, Liao [26] investigated the homotopy analysis method (HAM). The strength of HAM is that it leads to convergent analytic series solutions of strongly nonlinear problems faster than any other existing methods, independent of small or large physical parameter/s involved in the problem [27]. This behaviour of HAM makes it a superior technique to the conventional perturbation methods. The methods such as Adomian decomposition method [28–30], -expansion method [31], and Lyapunov artificial small parameter method [32] may not be valid for strong nonlinear problems due to the divergent nature of their obtained solution series. Liao [33], indeed, showed that HAM is the general case and the Adomian decomposition method, -expansion method and Lyapunov artificial small parameter method are the special cases of HAM. Moreover, He’s homotopy perturbation method (HPM) [34, 35] is also a special case of the HAM (cf. Liao [36]). Actually, for some auxiliary linear operators, the traditional HPM turns out to be the Taylor series expansion (cf. Turkyilmazoglu [37]).

The HAM has been applied to an extensive variety of nonlinear problems in science and engineering ever since it was first introduced in 1992. The problems of viscous flows of non-Newtonian fluids that have been mainly tackled by Hayat and his coworkers [38–40] and problems of heat transfer [41, 42] are some of the examples of applications of HAM. Even a much wider range of applications of HAM can be found in [27].

2. Governing Equations

We consider here the boundary layer flow of a fluid which is electrically conducting. The flow is assumed to be steady, laminar, incompressible, and axisymmetric in a circular cone having three-dimensional sink at the vertex (Figure 1).

Figure 1: Schematic diagram of the flow model.

Here is the magnetic field applied in -direction and in comparison with the fluid; it is fixed. With a specific end goal to disregard the induced magnetic field as compared to the applied magnetic field, the magnetic Reynolds number is considered to be small. The temperature and concentration of the wall and the free stream are preserved at constant level. Importance has been given for studying the effects of suction/injection, magnetic, and radiation parameters and hence they have been included in the analysis. On the other hand, the dissipation terms and the Hall effect have been neglected. It is also considered that the injected gas and the boundary layer gas both possess the same physical properties and are perfect gases. Also, the static temperature of injected gas is equal to the wall temperature. The basic boundary layer equations for the MHD flow due to point sink, including the radiative effect of the fluid at the wall, are (cf. Takhar et al. [4] and Vyas and Rai [9])

where

Here the subscripts and denote derivatives with respect to and , respectively, and represents the static pressure.

The boundary conditions are given by

We now use the following similarity transformations in (1)–(6):

Here is the distance along the cone from the vertex and is the corresponding velocity component along direction, whereas is the distance perpendicular to the cone; is velocity component along direction; is radius of the cone given by , where is the semivertical angle of the cone (see Figure 1); is the dimensional stream function and represents the corresponding dimensionless stream function; is the dimensional concentration and is the corresponding dimensionless concentration; is the dimensional temperature while represents dimensionless temperature; , , and are the electrical conductivity, density, and kinematic viscosity, respectively; represents similarity variable; stands for magnetic field; and represent the binary diffusion coefficient and thermal diffusivity, respectively; is the strength of point sink; is the magnetic parameter; represents the inviscid flow velocity; denotes the mass transfer parameter; is the specific heat at constant pressure; is the radiative heat flux; is coefficient of thermal conductivity; is radiation parameter; is the Stefan-Boltzmann constant, and is the mean absorption coefficient.

By using the Rosseland approximation (Brewster [43]), the radiative heat flux is given by

Here it should be noted that by using the Rosseland approximation, the present analysis is limited to optically thick fluids. If the temperature differences within the flow are sufficiently small, then (8) can be linearized by expanding into the Taylor series about , which after neglecting the higher order terms takes the form

In view of the result given in (9) and the similarity transformations given in (7), the continuity equation (1) is satisfied identically, and (2)–(4) get reduced to self-similar equations given by

Also, the boundary conditions (6) get reduced to

where and are Schmidt and Prandtl numbers, respectively; the subscripts denote conditions at the wall and denote conditions in the free stream; and prime denotes derivative with respect to .

As pointed out by Rosenhead [2], it is to be noted that the abovementioned boundary layer approximation is not valid in the immediate neighbourhood of the hole (cf. Rosenhead [2], pp. 428). Also for our mathematical analysis, it is remarked that the mass transfer parameter is treated as constant (Ref. Takhar et al. [4]). Also, the magnetic parameter can be treated locally as a constant for a fixed as it was first considered by Takhar et al. [44] and then by Takhar and Nath [45]. Also in a sink flow, is referred to as suction and as injection (cf. Takhar et al. [4] and Schlichting and Gersten [46], pp. 294–298).

3. Homotopy Analysis

3.1. Solution for Skin Friction ((10) and (13))

In order to find the analytic solution of (10) along with the boundary conditions given by (13), we first select the linear operator as

and we choose as an embedding parameter. We, now, construct the following zeroth-order deformation equation:

with boundary conditions

where the prime denotes the partial derivative with respect to and

When , we have

and when , we have

Hence, as varies from to , varies from initial solution to the exact solution .

Here we choose , the initial guess, such that it satisfies and the boundary conditions (13). We select

We here assume that the -order deformation derivative given by

exists. By using (20) and Taylor’s formula, we have

We here assume that both and are properly chosen in such a way that the series (25) is convergent at . From (21) and (23) at , we find the following relationship between known initial solution and the unknown solution :

where we define

In order to find the -order deformation equation, we first differentiate (17) and (18) times with respect to and then we set , and finally we divide it by , to obtain

with the corresponding boundary conditions

where prime denotes the partial derivative with respect to .

Using (23) and (30), we can first calculate and then by solving linear equation (28) with boundary conditions (29), we can find . In a similar manner, we can calculate by using (31) and then find and so on. In this way, we solve first few deformation equations by using symbolic software MATHEMATICA and observe that can be explicitly expressed as

where

The first order of approximation of (10) and (13), calculated by using MATHEMATICA, is

After knowing the structure (32) of , we deduce a recurrence formula for the coefficients of rigorously. Kindly refer to the Appendix for details. Hence, we obtain an analytic solution explicitly for the MHD flow equations with heat and mass transfer due to a point sink, governed by (10) and (13), as

The solution (35) in the form of infinite series contains two parameters and .

3.2. Convergence of the Analytic Solution

As suggested by Liao [33], the auxiliary parameter plays a significant role in controlling the convergence and the rate of approximation for the HAM. It is also to be noted that the HAM provides a great deal of flexibility and freedom for choosing appropriate values of and so as to ensure the convergence of the solution, obtained in the form of infinite series, to . To choose , Liao introduced the concept of -curve that gives an admissible range, called convergence region, for the selection of the suitable values of . It is here to be emphasized that there exists a value of parameter for every belonging to the convergence region which is the most suitable value in the sense that it guarantees the fastest convergence of the given series.

In the present analysis, the values of for each particular solution of the distinct magnetic parameter and mass transfer parameter have been found by way of selecting suitable values of and with the help of -curves given in Figures 2–4. These values have been found after appropriate orders of approximations. These values agree well with the corresponding numerical values of Takhar et al. [4] who tackled the problem by shooting method in conjunction with Runge-Kutta fourth order method. Figure 7 is also an indicator of the fact that obtained from the series given in (35) is also in good agreement with Takhar et al.’s numerical solution.

Figure 2: -curve for , and .

Figure 3: -curve for , and .

Figure 4: -curve for , and .

It is obvious that if (35) converges, its second-order derivative with respect to at , say,

must converge. By using (A.5), we have its -order approximation as

In this context, it is worth mentioning that Turkyilmazoglu [47] has also proposed a new and novel way of finding the optimum value of convergence control parameter to ensure the convergence of the HAM series in the fastest manner. The proposed method constitutes an alternative to the classical -level curves method (refer to Liao [33]) and the squared residual approach proposed by Liao [48] for the determination of optimal value of the convergence control parameter.

3.3. Solutions for Temperature and Concentration ((11)–(15))

The heat and concentration equations can now easily be solved, for the velocity equation (10) with boundary conditions (13) already possesses a suitable and accurate solution. It is very much precise to obtain the solution of (11) by using boundary conditions (14), as

So, the heat transfer parameter is given by

where . Similarly, the corresponding solution of (12), along with the conditions (15), can be obtained from (38)-(39) simply by replacing by and by .

4. Results and Discussion

By using the recurrence formula (37), the numerical values of the skin-friction parameter for different values of magnetic parameter and mass transfer parameter (for suction and injection both) have been calculated and are given in Tables 1 and 2. From Table 1, it is obvious that the skin-friction parameter increases along with the increasing values of the mass suction parameter and magnetic parameter . From Table 2, it is clear that the skin-friction parameter increases along with the increasing values of the magnetic parameter but decreases with the increasing values of the mass injection parameter . The reason for such a behaviour is that both the suction and the magnetic parameters reduce the thickness of the momentum boundary layer which results in an increase in skin friction. The effect of injection is just opposite.

Table 1: Values of for different values of and

Present
Takhar et al. [4]

0
−2
3.5211
3.5182

−1
2.8517
2.8772

0.5
−2
3.6172
3.6162

−1
2.9554
3.0231

1
−2
3.7098
3.7124

−1
3.0542
3.1121

Table 2: Values of for different values of and

Present
Takhar et al. [4]
Rosenhead [2]

0
0
2.2721
2.2728
2.273

1
1.7861
1.7505

2
1.4167
1.4121

0.5
0
2.3827
2.392

1
1.9117
1.973

2
1.5232
1.5529

1
0
2.4604
2.4552

1
2.0158
2.0825

2
1.6252
1.6345

The numerical values of the skin-friction parameter for the present case and for the case of Takhar et al. [4] have been given in the Tables 1 and 2 for different values of mass suction/injection and magnetic parameters. These results are in excellent agreement with each other. The skin-friction parameter and the mass flux diffusion parameter () have been plotted in Figures 5 and 6, respectively, for different values of the magnetic and mass transfer parameters. The graphs thus obtained exhibit excellent agreement with the corresponding graphical results obtained by [4] who visited the problem in the absence of radiation parameter ().

Figure 5: Variation of skin friction with for (cf. [4]).

Figure 6: Variation of heat transfer and mass diffusion parameters with for and (cf. with [4]).

Figure 7: Variation of velocity profiles with for .

From Figure 6, it is evident that the heat transfer parameter () decreases for increasing values of the magnetic parameter but increases with increasing values of mass suction (). On the other hand, the parameter decreases with the increasing values of both the magnetic and mass injection parameters. As a result, the suction () reduces the thermal boundary layer, whereas the injection () and the magnetic parameters increase them.

It is here to be noted that the mass flux diffusion parameter () is similar to the heat transfer parameter () in the absence of radiative fluid properties (i.e., for ). So, the parameter will exhibit similar behaviour as shown in Figure 6.

From Figure 7, it is evident that the velocity profiles exhibit an increasing trend with increasing values of the magnetic parameter () in both the cases of suction and injection.

From Figure 7, it is also clear that the velocity profiles due to suction () are steeper than those due to injection ().

As the velocity increases with an increase in magnetic parameter (), the thickness of the momentum boundary layer also decreases. This happens due to Lorentz’s force arising from the interaction of the magnetic and electric fields during the motion of the electrically conducting fluid.

In Figure 8, the nature of temperature profiles with the Prandtl number () for a fixed value of magnetic parameter () has been studied in the absence of the radiative fluid properties (i.e., ). From Figure 8, it is clear that the Prandtl number and hence the Schmidt number , respectively, have significant effects on temperature and concentration profiles. Both and , respectively, increase the temperature and concentration profiles in case of injection. An opposite trend is observed in case of suction.

Figure 8: Variation of temperature profiles for , and .

In Figure 9, the effects of suction and injection parameters on temperature profiles have been studied for fixed values of the magnetic parameter (), radiation parameter (), and Prandtl number (). In both the cases of suction and injection, the temperature profiles are showing a decreasing trend along with the similarity variable . However, the decreasing trend is more effective in case of suction as compared to that of injection.

Figure 9: Variation of temperature profiles with for .

Figures 10 and 11 show the variation of temperature profiles with Prandtl number for fixed values of the magnetic and radiation parameters. Figure 10 shows that the temperature profiles exhibit an increasing trend with increasing in case of injection. But an opposite trend is observed in case of suction, as is clear from the Figure 11.

Figure 10: Variation of temperature profiles with for and .

Figure 11: Variation of temperature profiles with for and .

In Figures 12 and 13, the variation of the temperature profiles along with the radiation parameter () has been shown for fixed values of magnetic parameter () and Prandtl number () in cases of injection and suction both. From Figure 12, it is obvious that the temperature profiles exhibit an increasing trend with increasing values of in case of injection. But an opposite trend is observed in case of suction (see Figure 13).

Figure 12: Variation of temperature profiles with for , and .

Figure 13: Variation of temperature profiles with for , and .

5. Concluding Remarks

(1)Our results for have been found in good agreement with those obtained by Takhar et al. [4] which verifies the great potential and validity of the HAM.(2)The skin-friction increases with increasing magnetic field. The skin friction is greater for suction parameter as compared to injection parameter .(3)The skin-friction (), heat transfer parameter (), and mass flux diffusion parameter () decrease by injection and increase by suction .(4)The effect of magnetic field on the skin friction, heat transfer, and mass flux diffusion is less than the effect of mass transfer.(5)The temperature and concentration profiles get affected by Prandtl number and Schmidt number, respectively.(6)To get better results, that is, to get better approximations, HAM can offer us large flexibility and great freedom to choose better auxiliary linear operator (), nonzero auxiliary parameters (), and spatial-scale parameter () for satisfying the rule for solution expression and initial approximations.(7)With the help of high-speed computers and symbolic computation software like MATHEMATICA, Maple, and so forth, the HAM might become more powerful and perfect analytic tool to solve rigorous non-linear problems in science and engineering.(8)HAM is much superior and accurate method as compared to the shooting method used by Takhar et al. [4] to solve the present problem, for the shooting method gives numerical solution to the problem and HAM gives an analytic solution very near to the exact solution of the problem. So, the results obtained in the present study by using HAM are more credible than those of Takhar et al. [4]. The essence of the shooting method to solve a boundary value problem (BVP) is to convert it into a system of initial value problems where systematic guesses are made for unknown quantities till the end conditions are satisfied. As a consequence of this, this method involves a large number of repetitions if the initial value of the step size (, say) is far away from the optimal value. So, for the solution of every boundary value problem, there is an optimum step size. For large , the accuracy in the results diminishes, and if is too small, the round-off errors would dominate and reduce the accuracy of results.

Appendix

Derivation of the Coefficients in (32)

(i) It is to be noted that the real valued function is defined by (33) and the initial guess defined in (23) has the same structure as (32).

(ii) It can be proved that too has the same structure as (32), if we consider that the solutions have the same structure as (32).

We redefine the real function by introducing for the sake of ease of calculations as (cf. Liao [25])

where

Now we compute , and as follows:

By using (A.1)–(A.3), we have

where

Now we build recurrence formulae for in order to find . When , we get by using (32), (A.2), and (A.5)

Hence,

where

Here, we point out that for all . Hence, we have from (A.9) and (A.10) for

Similarly, we have

where

Now, (A.12) leads to

But it is observed that for

So, we have

Finally, putting all these in (31), becomes

where, for ,

and, for ,

Hence, by substituting (A.16) into (28), we have the -order deformation equation

That is,

where and is given by (16).

Hence,

where , and are the integral constants:

Now applying boundary conditions we get

This shows that has the same structure as (32) and hence the related coefficients can be obtained by the following recurrence relations:

(iii) In (i), we remarked that the initial guess approximation has the same structure as (32). In (ii), we derived all the recurrence formulae (A.1)–(A.24) and also we have proved that too should have the same structure as (32), if we consider that the first solutions have the same structure as (32). Hence, from (i) and (ii), all have the same structure as (32). Along these lines, we can figure all coefficients in a steady progression by utilizing the initial four known coefficients:

which are dictated by the initial guess estimation characterized by (13).

Competing Interests

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

References

S. H. Choi and H. E. Wilhelm, “Self‐similar magnetohydrodynamic diffuser flows with induced magnetic fields,” Physics of Fluids, vol. 20, no. 11, pp. 1821–1826, 1977.

L. Rosenhead, Laminar Boundary Layers, Clarendon Press, Oxford, UK, 1963.

R. C. Ackerberg, “The viscous incompressible flow inside a cone,” Journal of Fluid Mechanics, vol. 21, pp. 47–81, 1965.

H. S. Takhar, C. D. Surma Devi, and G. Nath, “MHD flow with heat and mass transfer due to a point sink,” Indian Journal of Pure and Applied Mathematics, vol. 17, no. 10, pp. 1242–1247, 1986.

A. T. Eswara, S. Roy, and G. Nath, “Unsteady MHD forced flow due to a point sink,” Acta Mechanica, vol. 145, no. 1–4, pp. 159–172, 2000.

A. T. Eswara and B. C. Bommaiah, “The effect of variable viscosity on laminar flow due to a point sink,” Indian Journal of Pure and Applied Mathematics, vol. 35, no. 6, pp. 811–815, 2004.

M. Turkyilmazoglu, “Exact multiple solutions for the slip flow and heat transfer in a converging channel,” Journal of Heat Transfer, vol. 137, no. 10, Article ID 101301, 2015.

E. Magyari, “Backward boundary layer heat transfer in a converging channel,” Japan Society of Fluid Mechanics, vol. 39, no. 6, pp. 493–504, 2007.

P. Vyas and A. Rai, “Radiative variable fluid properties flow due to a point sink inside a cone filled with porous medium,” Applied Mathematical Sciences, vol. 6, no. 87, pp. 4307–4317, 2012.

E. M. Sparrow and R. D. Cess, Radiation Heat Transfer, Brooke/Cole, Belmont, Calif, USA, 1970.

J. R. Howell, “Radiative transfer in porous media,” in Handbook of Porous Media, K. Vafai, Ed., pp. 663–698, CRC Press, New York, NY, USA, 2000.

O. A. Plumb, J. S. Huenefeld, and E. J. Eschbach, “The effect of crossbow and radiation on natural convection from vertical heated surfaces in saturated porous media,” in Proceedings of the AIAA 16th Thermophysics Conference, Polo Alto, Calif, USA, 1981.

M. A. Hossain and H. S. Takhar, “Radiation effect on mixed convection along a vertical plate with uniform surface temperature,” Heat and Mass Transfer, vol. 31, no. 4, pp. 243–248, 1996.

A. Raptis, “Radiation and free convection flow through a porous medium,” International Communications in Heat and Mass Transfer, vol. 25, no. 2, pp. 289–295, 1998.

M. A. Seddeek and A. M. Salem, “Laminar mixed convection adjacent to vertical continuously stretching sheets with variable viscosity and variable thermal diffusivity,” Heat and Mass Transfer, vol. 41, no. 12, pp. 1048–1055, 2005.

M. Q. Al-Odat, F. M. S. Al-Hussien, and R. A. Damseh, “Influence of radiation on mixed convection over a wedge in non-Darcy porous medium,” Forschung im Ingenieurwesen, vol. 69, no. 4, pp. 209–215, 2005.

V. R. Prasad, N. B. Reddy, and R. Muthucumaraswamy, “Radiation and mass transfer effects on two-dimensional flow past an impulsively started infinite vertical plate,” International Journal of Thermal Sciences, vol. 46, no. 12, pp. 1251–1258, 2007.

S. Mukhopadhyay, “Effects of radiation and variable fluid viscosity on flow and heat transfer along a symmetric wedge,” Journal of Applied Fluid Mechanics, vol. 2, no. 2, pp. 29–34, 2009.

P. Vyas and N. Srivastava, “Radiative MHD flow over a non-isothermal stretching sheet in a porous medium,” Applied Mathematical Sciences, vol. 4, no. 50, pp. 2475–2484, 2010.

P. Vyas and A. Ranjan, “Dissipative MHD boundary-layer flow in a porous medium over a sheet stretching nonlinearly in the presence of radiation,” Applied Mathematical Sciences, vol. 4, no. 61–64, pp. 3133–3142, 2010.

D. S. Chauhan and V. Kumar, “Radiation effects on mixed convection flow and viscous heating in a vertical channel partially filled with a porous medium,” Tamkang Journal of Science and Engineering, vol. 14, no. 2, pp. 97–106, 2011.

I. G. Baoku, C. Israel-Cookey, and B. I. Olajuwon, “Influence of thermal radiation on a transient MHD couette flow through a porous medium,” Journal of Applied Fluid Mechanics, vol. 5, no. 1, pp. 81–87, 2012.

P. R. Babu, J. A. Rao, and S. Sheri, “Radiation effect on Mhd heat and mass transfer flow over a shrinking sheet with mass suction,” Journal of Applied Fluid Mechanics, vol. 7, no. 4, pp. 641–650, 2014.

S.-J. Liao, “An explicit, totally analytic approximate solution for Blasius' viscous flow problems,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 759–778, 1999.

S.-J. Liao, “A uniformly valid analytic solution of two-dimensional viscous flow over a semi-infinite flat plate,” Journal of Fluid Mechanics, vol. 385, pp. 101–128, 1999.

S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, 1992.

S. J. Liao, “Notes on the homotopy analysis method: some definitions and theorems,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 983–997, 2009.

G. Adomian, “Nonlinear stochastic differential equations,” Journal of Mathematical Analysis and Applications, vol. 55, no. 2, pp. 441–452, 1976.

G. Adomian, “A review of the decomposition method and some recent results for nonlinear equations,” Computers & Mathematics with Applications, vol. 21, no. 5, pp. 101–127, 1991.

R. Rach, “On the Adomian (decomposition) method and comparisons with Picard's method,” Journal of Mathematical Analysis and Applications, vol. 128, no. 2, pp. 480–483, 1987.

V. M. Soundalgekar, M. Singh, and H. S. Takhar, “MHD free convection past a semi-infinite vertical plate with suction and injection,” Nonlinear Analysis: Theory, Methods & Applications, vol. 7, no. 9, pp. 941–944, 1983.

A. M. Lyapunov, General Problem on Stability of Motion, Taylor & Francis, London, UK, 1992.

S. Liao, Beyond Perturbation: Introduction to Homotopy Analysis Method, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2003.

J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.

J.-H. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.

S. J. Liao, “Comparison between the homotopy analysis method and homotopy perturbation method,” Applied Mathematics and Computation, vol. 169, no. 2, pp. 1186–1194, 2005.

M. Turkyilmazoglu, “Is homotopy perturbation method the traditional Taylor series expansion,” Hacettepe Journal of Mathematics and Statistics, vol. 44, no. 3, pp. 651–657, 2015.

T. Hayat and M. Sajid, “Analytic solution for axisymmetric flow and heat transfer of a second grade fluid past a stretching sheet,” International Journal of Heat and Mass Transfer, vol. 50, no. 1-2, pp. 75–84, 2007.

T. Hayat, Z. Abbas, M. Sajid, and S. Asghar, “The influence of thermal radiation on MHD flow of a second grade fluid,” International Journal of Heat and Mass Transfer, vol. 50, no. 5-6, pp. 931–941, 2007.

T. Hayat and M. Sajid, “Homotopy analysis of MHD boundary layer flow of an upper-convected Maxwell fluid,” International Journal of Engineering Science, vol. 45, no. 2–8, pp. 393–401, 2007.

S. Abbasbandy, “The application of homotopy analysis method to nonlinear equations arising in heat transfer,” Physics Letters. A, vol. 360, no. 1, pp. 109–113, 2006.

S. Abbasbandy, “Homotopy analysis method for heat radiation equations,” International Communications in Heat and Mass Transfer, vol. 34, no. 3, pp. 380–387, 2007.

M. Q. Brewster, Thermal Radiative Transfer and Properties, John Wiley & Sons, 1992.

H. S. Takhar, “Hydromagnetic free convection from a at plate,” Indian Journal of Physics, vol. 45, pp. 289–311, 1971.

H. S. Takhar and G. Nath, “Similarity solution of unsteady boundary layer equations with a magnetic field,” Meccanica, vol. 32, no. 2, pp. 157–163, 1997.

H. Schlichting and K. Gersten, Boundary Layer Theory, Springer, Berlin, Germany, 2000.

M. Turkyilmazoglu, “An effective approach for evaluation of the optimal convergence control parameter in the homotopy analysis method,” Filomat, vol. 30, no. 6, pp. 1633–1650, 2016.

S. Liao, Advances in the Homotopy Analysis Method, World Scientific, Singapore, 2014.