Research Article  Open Access
On the Solution of DoubleDiffusive Convective Flow due to a Cone by a Linearization Method
Abstract
The paper details the use of a nonperturbation successive linearization method to solve the coupled nonlinear boundary value problem due to doublediffusive convection from an inverted cone. Diffusionthermo and thermaldiffusion effects have been taken into account. The governing partial differential equations are transformed into ordinary differential equations using a suitable similarity transformation. The SLM is based on successively linearizing the governing nonlinear boundary layer equations and solving the resulting higherorder deformation equations using spectral methods. The results are compared with the limited cases from previous studies and results obtained using the Matlab inbuilt bvp4c numerical algorithm and a shooting technique that uses RungeKuttaFehlberg (RKF45) and NewtonRaphson schemes. These comparisons reveal the robustness and validate the usage of the linearisation method technique. The results show that the nonperturbation technique in combination with the Chebyshev spectral collocation method is an efficient numerical algorithm with assured convergence that serves as an alternative to numerical methods for solving nonlinear boundary value problems.
1. Introduction
The convection driven by two different density gradients with differing rates of diffusion is widely known to as “doublediffusive convection’’ and is an important fluid dynamics phenomenon (see Mojtabi and CharrierMojtabi [1]). The study of doublediffusive convection has attracted attention of many researchers during the recent past due to its occurrence in nature and industry. Oceanography is the root of doublediffusive convection in natural settings. The existence of heat and salt concentrations at different gradients and the fact that they diffuse at different rates lead to spectacular doublediffusive instabilities known as “saltfingers’’ (see Stern [2, 3]). The formation of saltfingers can also be observed in laboratory settings. Doublediffusive convection occurs in the sun where temperature and helium diffusions take place at different rates. Convection in magma chambers and seawind formations are among other manifestations of doublediffusive convection in nature. Migration of moisture through air contained in fibrous insulations, grain storage systems, the dispersion of contaminants through watersaturated soil, crystal growth, the underground disposal of nuclear wastes, the formation of microstructures during the cooling of molten metals, and fluid flows around shrouded heatdissipation fins are among other industrial applications of doublediffusive convection.
The inherent instabilities due to doublediffusive convection have been investigated by, among others, Nield [4], Baines and Gill [5], Guo et al. [6], Khanafer and Vafai [7], Sunil et al. [8], and Gaikwad et al. [9]. Doublediffusive convection due to horizontal, inclined, and vertical surfaces embedded in a porous medium has been studied by, among others, Cheng [10, 11], Nield and Bejan [12], and Ingham and Pop [13]. Chamkha [14] investigated the coupled heat and mass transfer by natural convection of Newtonian fluids about a truncated cone in the presence of magnetic field and radiation effects. Yih [15] examined the effect of radiation in convective flow over a cone.
Though heat and mass transfer happens simultaneously in a moving fluid, the relations between the fluxes and the driving potentials are generally complicated. It should be noted that the energy flux can be generated by both temperature and composition gradients. The energy flux caused by a composition gradient gives rise to the Dufour or diffusionthermo effect. Mass fluxes created by temperature gradient lead to the Soret or thermaldiffusion effect. These effects are in collective known as crossdiffusion effects. The crossdiffusion effect has been extensively studied in gases, while the Soret effect has been studied both theoretically and experimentally in liquids, see Mortimer and Eyring [16]. They used an elementary transition state approach to obtain a simple model for Soret and Dufour effects in thermodynamically ideal mixtures of substances with molecules of nearly equal size. In their model, the flow of heat in the Dufour effect was identified as the transport of the enthalpy change of activation as molecules diffuse. The results were found to fit the Onsager reciprocal relationship, Onsager [17].
In general, the crossdiffusion effects are small compared to the effects described by Fourier and Fick’s laws (Mojtabi and CharrierMojtabi [1]) and can therefore be neglected in many heat and masstransfer processes. However, it has been shown in a number of studies that there are exceptions in areas such as in geosciences where crossdiffusion effects are significant and cannot be ignored, see for instance Kafoussias and Williams [18], Awad et al. [19], and the references therein. With this view point, many investigators included crossdiffusion effects in the study of doublediffusive convection in fluid flows involving bodies of various geometries. Alam et al. [20] investigated the Dufour and Soret effects on steady combined freeforced convective and mass transfer flow past a semiinfinite vertical flat plate of hydrogenair mixtures. They used the fourthorder RungeKutta method to solve the governing equations of motion. Their study showed that the Dufour and Soret effects should not be neglected. Shateyi et al. [21] investigated the effects of diffusionthermo and thermaldiffusion on MHD fluid flow over a permeable vertical plate in the presence of radiation and hall current. Awad and Sibanda [22] used the homotopy analysis method to study heat and mass transfer in a micropolar fluid subject to Dufour and Soret effects.
Most boundary value problems in fluid mechanics are solved numerically using either the shooting method or the implicit finite difference scheme in combination with a linearization technique. These methods have their associated difficulties and failures in handling situations where solutions either vary sharply over a domain or problems that exhibit multiple solutions. These limitations necessitate the development of computationally improved semianalytical methods for solving strongly nonlinear problems. There are many different semianalytical methods to solve nonlinear boundary value problems, among them, the variational iteration method, the homotopy perturbation method [23–25], the Adomian decomposition method [26, 27], homotopy analysis method [28], and the spectralhomotopy analysis methods [29, 30]. These iterative methods may sometimes fail to converge or give slow convergence for strongly nonlinear problems or problems involving large parameters. Yildirim [31] applied He’s homotopy perturbation method to solve the Cauchy reactiondiffusion problem and compared his results with analytical solutions in certain test cases. Yildirim and Pinar [32] obtained periodic solutions of nonlinear reactiondiffusion equations arising in mathematical biology using the expfunction method. Yildirim and Sezer [33] found analytical solutions of linear and nonlinear spacetime fractional reactiondiffusion equations (STFRDE) on a finite domain using the homotopy perturbation method (HPM). Yildirim et al. [34] presented approximate analytical solutions of the biochemical reaction model by the multistep differential transform method (MsDTM) and validated the results by comparing with the fourthorder RungeKutta method.
Ganji et al. [35] solved the nonlinear JefferyHamel flow problem using two semianalytical methods, the variational iteration method (VIM) and the homotopy perturbation method. Ghafoori et al. [36] solved the equation for a nonlinear oscillator using the differential transform method (DTM). They compared DTM solutions with those obtained using the variational iteration method and the homotopy perturbation method. Joneidi et al. [37] used three analytical methods, the homotopy analysis method (HAM), homotopy perturbation method, and the differential transform method, to solve the JefferyHamel flow problem. Babaelahi et al. [38] studied the heat transfer characteristics in an incompressible electrically conducting viscoelastic boundary layer fluid flow over a linear stretching sheet. They solved the flow equations using the optimal homotopy asymptotic method (OHAM) and validated their results by comparing the OHAM solutions with RungeKutta solutions.
In this study, we use a nonperturbation, semianalytic successive linearization method (see Makukula et al. [39, 40]) to investigate doublediffusive convection from a cone in a viscous incompressible fluid subject to crossdiffusion effects. The study is an extension of the work by Ece [41] to include mass transfer and crossdiffusion effects. The linearization method iteratively linearizes the nonlinear equations to give a system of higherorder deformation equations that are then solved using the Chebyshev spectral collocation method.
2. Mathematical Formulation
Consider a vertical downpointing cone with halfangle immersed in a viscous incompressible liquid. The axis is along the surface of the cone, and the axis coincides with the outward normal to the surface of the cone. The origin is at the vertex of the cone, see Figure 1. The surface of the cone is subject to a linearly varying temperature (>) where is the ambient temperature.
Following the usual boundary layer and Boussinesq approximations, the basic equations governing the steady state dynamics of a viscous incompressible liquid are given by where and are the velocity components in the and directions, respectively, is the local radius of the cone, is the kinematic viscosity, is the density, is the acceleration due to gravity, is the coefficient of thermal expansion, is the coefficient of solutal expansion, is the temperature, is the concentration, is the thermal diffusivity, is the species diffusivity, and , are crossdiffusion coefficients.
The boundary conditions for (2.1) have the form Here, the subscripts and refer to the surface and ambient conditions, respectively, and are positive constants, and is a characteristic length.
We introduce the dimensionless variables where the reference velocity and Grashof number are defined, respectively, as On using the variables (2.3), the boundarylayer equations (2.1) reduce to The nondimensional parameters appearing in (2.5)–(2.8) are the buoyancy ratio , the Prandtl number , the Schmidt number , Dufour number D_{f}, and Soret number defined, respectively, as Assuming and , the boundary conditions (2.2) can be written as We now introduce the stream function such that so that the continuity equation (2.5) is satisfied identically. The boundary layer equations (2.6)–(2.8) can be written in terms of the stream function as The boundary conditions (2.10) in terms of the stream function are We further introduce the following similarity variables Using (2.14), (2.12) along with boundary conditions (2.13) reduces to the following twopoint boundary value problem The primes in (2.15)–(2.18) denote differentiation with respect to .
3. Successive Linearization Method
The successive linearization method (see Makukula et al. [39, 40]) is used to solve the boundary value problem (2.15)–(2.18). We assume that the functions , , and may be expanded in series form as where , and are unknown functions and , and , are approximations that are obtained by recursively solving the linear part of the equation that results from substituting (3.1) in (2.15)–(2.17). Substituting (3.1) in the governing (2.15)–(2.17), we obtain where the coefficient parameters (), , (), and () are defined as
Starting from the initial approximations which are chosen to satisfy the boundary conditions (2.18), the subsequent solutions , , , are obtained by successively solving the linearized form of (3.2) given below subject to the boundary conditions Once each solution , , and () has been found from iteratively solving (3.5) for each , the functions , , and are obtained as series where is the order of SLM approximation. Equations (3.5) are integrated using the Chebyshev spectral collocation method [42–44]. The unknown functions are defined by the Chebyshev interpolating polynomials with the GaussLobatto points defined as where is the number of collocation points used. The physical region is transformed into the domain using the domain truncation technique in which the problem is solved on the interval instead of . This leads to the mapping where is the known number used to invoke the boundary condition at infinity. The unknown functions , , and are approximated at the collocation points by where is the th Chebyshev polynomial defined as The derivatives of the variables at the collocation points are represented as where is the order of differentiation and where is the Chebyshev spectral differentiation matrix (see, [42–44]). Substituting (3.8)–(3.12) in (3.5)–(3.6) leads to the matrix equation in which is a square matrix of order and , are column vectors defined by with
In the above definitions, (), , (), , and are diagonal matrices of order and is the identity matrix of order . Finally, the solution of the problem is obtained as Thus, starting with the initial solutions , and , a sequence of approximations , , and , are obtained until (3.7) holds. The convergence of this iteration process depends on the parameter values, that is, for small parameter values, the iterates converge faster as compared to large parameter values.
4. Skin Friction, Heat and Mass Transfer Coefficients
The parameters of engineering interest in heat and mass transport problems are the skin friction coefficient , the Nusselt number and the Sherwood number . These parameters characterize the surface drag, the wall heat and mass transfer rates, respectively.
The shearing stress at the surface of the cone is defined as where is the coefficient of viscosity. The skin friction coefficient at the surface of the cone is defined as Using (4.1) in (4.2), we obtain the following relation The heat transfer rate at the surface of the cone is defined as where is the thermal conductivity of the fluid. The Nusselt number is defined as Using (4.4) in (4.5), the dimensionless wall heat transfer rate is obtained as follows: The mass flux at the surface of the cone is defined as and the Sherwood is defined as Using (4.7) in (4.8), the dimensionless wall mass transfer rate is obtain as
5. Results and Discussion
The successive linearization method (SLM) has been applied to solve the nonlinear coupled boundary value problem arising due to doublediffusive convection from a vertical cone immersed in a viscous liquid. Crossdiffusion effects are taken into consideration. The parameters controlling the flow dynamics are the Prandtl number , Schmidt number , buoyancy ratio , Dufour number D_{f,} and the Soret number . We, however, do not discuss the effects of parameters such as the Prandtl and Schmidt numbers whose significance has been widely studied in the literature on doublediffusive convection in viscous liquids. We have thus fixed and instead focus attention on results pertaining to the other three important parameters. In addition, we restrict ourselves to parameter values in the interval , . To highlight the effect of buoyancy, for aiding buoyancy condition, we take while, for opposing buoyancy, .
We first establish the robustness and accuracy of the successive linearization method (SLM) by comparing the SLM results with those obtained numerically and previous related studies in the literature. The Matlab inbuilt bvp4c routine and the shooting technique with RungeKuttaFehlberg (RKF45) and NewtonRaphson schemes are used to obtain the numerical solutions.
Tables 1 and 2 show the results of , , and for different parameter values. Table 1 gives the comparison of the SLM results in the absence of crossdiffusion (i.e., ) with those presented by Ece [41]. The SLM solutions are found to be in excellent agreement with those of Ece [41] indicating the accuracy of the linearisation method.


Table 2 highlights both the accuracy and the accelerated convergence of the SLM for different values of D_{f} and . The linearisation method converges to the numerical solutions at the fourthorder SLM for all values of D_{f} and . However, for larger values, convergence may require extra terms in the SLM solution series. It is evident that the SLM results are highly accurate as they match with those obtained by the bvp4c and the shooting technique up to the sixth significant digit.
It is to be noted from Table 2 that simultaneously increasing D_{f} and decreasing lead to initial decreases in the skinfriction coefficient up to and then start increasing. The heat transfer coefficient shows monotonic decrease, while the mass transfer coefficient exhibits the opposite change when subjected to simultaneous increase in D_{f} and decrease in .
To gain some insight into the dynamics of the problem, the temperature and concentration distributions are shown graphically in Figures 2–6. The Nusselt number and Sherwood number which highlight the heat and mass transfer are shown in Figures 7 and 8, as functions of for different values of D_{f} in the aiding and opposing buoyancy cases.
(a)
(b)
The variation of temperature and concentration profiles subject to a simultaneous increase in the crossdiffusion parameters D_{f} and is shown in Figure 2. We observe enhanced heat and mass transfer in the presence of the crossdiffusion effect as compared to the case (no crossdiffusion). Increasing the crossdiffusion parameters increases both the thermal and species boundary layer thickness in both the aiding and opposing buoyancy situations. Hence, the crossdiffusion effect plays an important role in enhancing heat and mass transfer in doublediffusion convection processes.
Figure 3 shows the effect of the Dufour number on the temperature distributions. The energy flux created by the concentration gradient gives rise to the Dufour effect or diffusionthermo effect and due to the increase in the energy flux created by concentration gradients, the temperature in the boundary layer increases significantly. The Dufour effect thus serves to thicken the thermal boundary layer. This trend is true for both aiding and opposing buoyancy scenarios.
Due to the coupling between the momentum, energy, and species balance equations, the Dufour parameter has an effect on the concentration boundary layer as well. This is shown in Figure 4 where it is evident that D_{f} reduces the concentration in the boundary layer in both the cases of aiding and opposing buoyancy.
The effect of the Soret number on the temperature distribution is shown in Figure 5. The Soret parameter has a mixed effect on profiles. In the case of opposing buoyancy, increasing Soret parameter results in the thickening of the thermal boundary layer, while, in the aiding buoyancy case, the effect of is exactly the opposite.
Figure 6 shows the effect of the Soret number on the species distribution in aiding and opposing buoyancy cases. The mass flux created by the temperature gradient gives rise to Soret or thermaldiffusion or thermophoresis effect. The thermophoretic force developed due to temperature gradients drives solute particles into the boundary layer region thereby increasing the concentration boundary layer as can be seen from Figure 6. The increase in concentration boundary with is observed in both aiding and opposing buoyancy cases.
Figure 7 shows the Nusselt number as a function of for different values of D_{f} in aiding and opposing buoyancy conditions. In the opposing buoyancy situation, decreases with for the case of pure thermophoresis () and increases in the crossdiffusion case (). In the aiding buoyancy situation, increases monotonically with for both and . The Dufour number reduces the heat transfer coefficient in both aiding and opposing flow situations. Further, we observe enhanced heat transfer in the case of aiding buoyancy () as compared to the opposing buoyancy () case.
Figure 8 shows the mass transfer coefficient as a function of for different values of D_{f} in aiding and opposing buoyancy conditions. In both aiding and opposing buoyancy situations, is a decreasing function of and an increasing function of D_{f}. There is also an increased mass transfer in the case of aiding buoyancy () as compared to the opposing buoyancy () case.
6. Conclusions
The problem of doublediffusive convection from a vertical cone was solved using a successive linearization algorithm in combination with a Chebyshev spectral collocation method. A comparison with results in the literature and numerical approximations showed that the SLM is highly accurate with assured and accelerated convergence rate thus confirming the SLM as an alternative semianalytic technique for solving nonlinear boundary value problems with a strong coupling. We found that the Dufour parameter reduces the heat transfer coefficient while increasing the mass transfer rate. In general, the effect of the Soret parameter is to increase the heat transfer coefficient and to reduce the mass transfer coefficient. Aiding buoyancy enhances heat and mass transfer compared to the opposing buoyancy condition.
Nomenclature
:  Concentration 
:  Dimensionless concentration 
:  Local skin friction coefficient 
:  Concentration difference, 
:  Boundary layer stream function 
:  Solutal diffusivity 
D_{f}:  Dufour number 
:  Acceleration due to gravity 
Grashof number  
:  Mass flux 
:  Thermal conductivity 
, :  Crossdiffusion coefficients 
:  Characteristic length 
:  Order of successive linearization method 
:  Number of collocation points 
Local Nusselt number  
Prandtl number  
:  Heat flux 
:  Local radius of the cone, 
:  Dimensionless local radius of the cone, 
Schmidt number  
Local Sherwood number  
Soret number  
:  Temperature 
:  Dimensionless temperature 
:  Temperature difference, 
:  Reference velocity 
, :  Velocity component in the directions 
, :  Dimensionless velocity component in the directions 
, :  Coordinate measured along the surface and normal to it 
, :  Dimensionless coordinates. 
:  Thermal diffusivity of the fluid 
:  Coefficient of thermal expansion of the fluid 
:  Coefficient of solutal expansion 
:  Vertex half angle of the cone 
:  Buoyancy ratio 
:  Coefficient of viscosity 
:  Coefficient of kinematic viscosity, 
:  Boundary layer temperature 
:  Density of the fluid 
:  Dimensionless stream function 
:  Boundary layer concentration 
:  Collocation point 
:  Shearing stress. 
:  Quantities at the surface of the cone 
:  Quantities far away from the surface of the cone. 
Acknowledgments
The authors wish to thank University of KwaZuluNatal and the National Research Foundation (NRF) for financial support.
References
 A. Mojtabi and M. C. CharrierMojtabi, “Double diffusive convection in porous media,” in Handbook of Porous Media, pp. 269–320, 2nd edition, 2005. View at: Google Scholar
 M. E. Stern, “The “salt fountain” and thermohaline convection,” Tellus, vol. 12, pp. 172–175, 1960. View at: Publisher Site  Google Scholar
 M. E. Stern, “Collective instability of salt fingers,” Journal of Fluid Mechanics, vol. 35, no. 2, pp. 209–218, 1969. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. A. Nield, “Onset of thermohaline convection in a porous medium,” Water Resources Research, vol. 5, pp. 553–560, 1968. View at: Google Scholar
 P. G. Baines and A. E. Gill, “On thermohaline convection with linear gradients,” Journal of Fluid Mechanics, vol. 37, no. 2, pp. 289–306, 1969. View at: Publisher Site  Google Scholar
 J. Guo, Y. Qin, and P. N. Kaloni, “Nonlinear stability problem of a rotating doubly diffusive fluid layer,” International Journal of Engineering Science, vol. 32, no. 8, pp. 1207–1219, 1994. View at: Publisher Site  Google Scholar
 K. Khanafer and K. Vafai, “Doublediffusive mixed convection in a liddriven enclosure filled with a fluidsaturated porous medium,” Numerical Heat Transfer, Part A, vol. 42, no. 5, pp. 465–486, 2002. View at: Publisher Site  Google Scholar
 Sunil, A. Sharma, and R. C. Sharma, “Effect of dust particles on ferrofluid heated and soluted from below,” International Journal of Thermal Sciences, vol. 45, no. 4, pp. 347–358, 2006. View at: Publisher Site  Google Scholar
 S. N. Gaikwad, M. S. Malashetty, and K. R. Prasad, “An analytical study of linear and nonlinear double diffusive convection in a fluid saturated anisotropic porous layer with Soret effect,” Applied Mathematical Modelling, vol. 33, no. 9, pp. 3617–3635, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C. Ping, “Similarity solutions for mixed convection from horizontal impermeable surfaces in saturated porous media,” International Journal of Heat and Mass Transfer, vol. 20, no. 9, pp. 893–898, 1977. View at: Publisher Site  Google Scholar
 P. Cheng, “Natural convection in a porous medium: external flow,” in Proceedings of the NATO Advanced Study in Natural Convection, Ezmir, Turkey, 1985. View at: Google Scholar
 D. A. Nield and A. Bejan, Convection in Porous Media, Springer, Berlin, Germany, 3rd edition, 2010.
 D. B. Ingham and I. Pop, Transport Phenomenon in Porous Media, Elsevier, Oxford, UK, 2nd edition, 2002.
 A. J. Chamkha, “Coupled heat and mass transfer by natural convection about a truncated cone in the presence of magnetic field and radiation effects,” Numerical Heat Transfer, Part A, vol. 39, no. 5, pp. 511–530, 2001. View at: Google Scholar
 K. A. Yih, “Coupled heat and mass transfer by free convection over a truncated cone in porous media: VWT/VWC or VHF/VMF,” Acta Mechanica, vol. 137, no. 1, pp. 83–97, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. G. Mortimer and H. Eyring, “Elementary transition state theory of the Soret and Dufour effects,” Proceedings of the National Academy of Sciences of the United States of America, vol. 77, no. 4, pp. 1728–1731, 1980. View at: Publisher Site  Google Scholar
 L. Onsager, “Reciprocal relations in irreversible processes. I,” Physical Review, vol. 37, no. 4, pp. 405–426, 1931. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 N. G. Kafoussias and E. W. Williams, “Thermaldiffusion and diffusionthermo effects on mixed freeforced convective and mass transfer boundary layer flow with temperature dependent viscosity,” International Journal of Engineering Science, vol. 33, no. 9, pp. 1369–1384, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. G. Awad, P. Sibanda, and S. S. Motsa, “On the linear stability analysis of a Maxwell fluid with doublediffusive convection,” Applied Mathematical Modelling, vol. 34, no. 11, pp. 3509–3517, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. S. Alam, M. M. Rahman, M. A. Maleque, and M. Ferdows, “Dufour and Soret effects on steady MHD combined freeforced convective and mass transfer flow past a semiinfinite vertical plate,” Thammasat International Journal of Science and Technology, vol. 11, pp. 1–12, 2006. View at: Google Scholar
 S. Shateyi, S. S. Motsa, and P. Sibanda, “The effects of thermal radiation, hall currents, soret, and dufour on MHD flow by mixed convection over a vertical surface in porous media,” Mathematical Problems in Engineering, vol. 2010, Article ID 627475, 20 pages, 2010. View at: Publisher Site  Google Scholar
 F. G. Awad and P. Sibanda, “Dufour and Soret effects on heat and mass transfer in a micropolar fluid in a horizontal channel,” WSEAS Transaction on Heat and Mass Trasfer, vol. 5, pp. 165–177, 2010. View at: Google Scholar
 J. H. He and X. H. Wu, “Variational iteration method: new development and applications,” Computers and Mathematics with Applications, vol. 54, pp. 115–123, 2007. View at: Google Scholar
 J. H. He, “Coupling method of a homotopy technique and a perturbation technique for nonlinear problems,” International Journal of NonLinear Mechanics, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site  Google Scholar
 J. H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 12, pp. 87–88, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. Adomian, “Nonlinear stochastic differential equations,” Journal of Mathematical Analysis and Applications, vol. 55, no. 2, pp. 441–452, 1976. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Liao, Beyond Perturbation: Introduction to Homotopy Analysis Method, vol. 2 of CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2003.
 S. S. Motsa, P. Sibanda, and S. Shateyi, “A new spectralhomotopy analysis method for solving a nonlinear second order BVP,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 9, pp. 2293–2302, 2010. View at: Publisher Site  Google Scholar
 S. S. Motsa, P. Sibanda, F. G. Awad, and S. Shateyi, “A new spectralhomotopy analysis method for the MHD JefferyHamel problem,” Computers and Fluids, vol. 39, no. 7, pp. 1219–1225, 2010. View at: Publisher Site  Google Scholar
 A. Yildirim, “Application of He's homotopy perturbation method for solving the Cauchy reactiondiffusion problem,” Computers and Mathematics with Applications, vol. 57, no. 4, pp. 612–618, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Yıldırım and Z. Pınar, “Application of the expfunction method for solving nonlinear reactiondiffusion equations arising in mathematical biology,” Computers & Mathematics with Applications, vol. 60, no. 7, pp. 1873–1880, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Yildirim and S. A. Sezer, “Analytical solution of linear and nonlinear spacetime fractional reactiondiffusion equations,” International Journal of Chemical Reactor Engineering, vol. 8, article A110, 2010. View at: Google Scholar
 A. Yildirim, A. Gökdogan, and M. Merdan, “Numerical approximations to solution of biochemical reaction model,” International Journal of Chemical Reactor Engineering, vol. 9, article A92, 2011. View at: Google Scholar
 Z. Z. Ganji, D. D. Ganji, and M. Esmaeilpour, “Study on nonlinear JefferyHamel flow by He's semianalytical methods and comparison with numerical results,” Computers and Mathematics with Applications, vol. 58, no. 1112, pp. 2107–2116, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Ghafoori, M. Motevalli, M. G. Nejad, F. Shakeri, D. D. Ganji, and M. Jalaal, “Efficiency of differential transformation method for nonlinear oscillation: comparison with HPM and VIM,” Current Applied Physics, vol. 11, no. 4, pp. 965–971, 2011. View at: Publisher Site  Google Scholar
 A. A. Joneidi, G. Domairry, and M. Babaelahi, “Three analytical methods applied to JefferyHamel flow,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 11, pp. 3423–3434, 2010. View at: Publisher Site  Google Scholar
 M. Babaelahi, G. Domairry, and A. A. Joneidi, “Viscoelastic MHD flow boundary layer over a stretching surface with viscous and ohmic dissipations,” Meccanica, vol. 45, no. 6, pp. 817–827, 2010. View at: Publisher Site  Google Scholar
 Z. Makukula, S. Motsa, and P. Sibanda, “On a new solution for the viscoelastic squeezing flow between two parallel plates,” Journal of Advanced Research in Applied Mathematics, vol. 2, no. 4, pp. 31–38, 2010. View at: Publisher Site  Google Scholar
 P. Sibanda, Z. G. Makukula, and S. S. Motsa, “A novel numerical technique for twodimensional laminar flow between two moving porous walls,” Mathematical Problems in Engineering, vol. 2010, Article ID 528956, 15 pages, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. C. Ece, “Free convection flow about a cone under mixed thermal boundary conditions and a magnetic field,” Applied Mathematical Modelling, vol. 29, no. 11, pp. 1121–1134, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer Series in Computational Physics, Springer, Berlin, Germany, 1988.
 W. S. Don and A. Solomonoff, “Accuracy and speed in computing the Chebyshev collocation derivative,” SIAM Journal on Scientific Computing, vol. 16, no. 6, pp. 1253–1268, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L. N. Trefethen, Spectral Methods in MATLAB, vol. 10 of Software, Environments, and Tools, SIAM, Philadelphia, Pa, USA, 2000.
Copyright
Copyright © 2012 Mahesha Narayana and Precious Sibanda. 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.