## Services on Demand

## Article

## Indicators

- Cited by SciELO

## Related links

- Similars in SciELO

## Share

## Latin American applied research

##
*On-line version* ISSN 1851-8796

### Lat. Am. appl. res. vol.44 no.1 Bahía Blanca Jan. 2014

**A binary chemical reaction on unsteady free convective boundary layer heat and mass transfer flow with arrhenius activation energy and heat generation/absorption**

**Kh. A. Maleque**

*Professor of Mathematics, American International University-Bangladesh, House-23, 17, Kamal Ataturk avenue, Banani, Dhaka-1213, Bangladesh, email: maleque@aiub.edu, malequekh@gmail.com*

*Abstract* — We investigate a local similarity solution of an unsteady natural convection heat and mass transfer boundary layer incompressible fluid flow past a moving vertical porous plate in the presence of the heat absorption and generation. The effects of chemical reaction rate which is function of temperature and Arrhenius activation energy on the velocity, temperature and concentration are also studied in this paper. The governing partial differential equations are reduced to ordinary differential equations by introducing local similarity transformation (Maleque, 2010a). Numerical solutions to the reduced non-linear similarity equations are then obtained by adopting Runge-Kutta and shooting methods using the Nachtsheim- Swigert iteration technique. The results of the numerical solution are then presented graphically in the form of velocity, temperature and concentration profiles. The corresponding skin friction coefficient, the Nusselt number and the Sherwood number are also calculated and displayed in table showing the effects of various parameters on them.

*Keywords* — Heat and Mass Transfer; Chemical Reaction; Activation Energy; Porous Plate; Heat Generation/Absorption.

**I. INTRODUCTION**

In the past many researchers have published their papers on combined heat and mass transfer laminar free convective boundary layer flows for plate with various orientations. More information on this subject can be found in the monographs by Blums *et al.* (1987) and Gebhart *et al. *(1988) and in the papers by Khair and Bejan (1985), Mongruel *et al.* (1996), Lin and Wu (1997), Hossain and Rees(1999), Sattar and Maleque (2000), Maleque and Sattar (2002), Chaudhary and Sharma (2006), Ezzat *et al*. (2011), Maleque (2009, 2010b, 2011, 2013a).

In free convection boundary layer flows with simultaneous heat mass transfer, one important criteria that is generally not encountered is the species chemical reactions with finite Arrhenius activation energy. The Arrhenius law is usually of the form (Tencer *et al., *2004).

(1) |

where *K* is the rate constant of chemical reaction and *B* is the pre-exponential factor simply prefactor (constant), is based on the fact that increasing the temperature frequently causes a marked increase in the rate of reactions. *E _{a}* is the activation energy and

*k=*8.61×10

^{-5}eV/°K is the Boltzmann constant. In areas such as geothermal or oil reservoir engineering, the above phenomenon is usually applicable. Apart from experimental works in these areas, it is also important to make some theoretical efforts to predict the effects of the activation energy in flows mentioned above. But in this regard very few theoretical works are available in literature. The reason is that the chemical reaction processes involved in the system are quite complex and generally the mass transfer equation that is required for all the reactions involved also become complex. Theoretically, such an equation is rather impossible to tackle. Form chemical kinetic viewpoint this is a very difficult problem, but if the reaction is restricted to binary type a lot of progress can be made. The thermomechanical balance equations for a mixture of general materials were first formulated by Truesdell (1957). Thereafter Mills (1966) and Beevers and Craine (1982) have obtained some exact solutions for the boundary layer flow of a binary mixture of incompressible Newtonian fluids. Several problems relating to the mechanics of oil and water emulsions, particularly with regard to applications in lubrication practice, have been considered within the context of a binary mixture theory by Al- Sharif

*et al.*(1993) and Wang

*et al*. (1993). A simple model involving binary reaction was studied by Bestman (1990). He considered the motion through the plate to be large which enabled him to obtain analytical solutions (subject to same restrictions) for various values of activation energy by employing the perturbation technique proposed by Singh and Diskshit (1988). Bestman (1991) and Alabraba

*et al*. (1992) took into account the effect of the Arrhenius activation energy under the different physical conditions. Kandasamy

*et al.*(2005) studied the combined effects of chemical reaction, heat and mass transfer along a wedge with heat source and concentration in the presence of suction or injection. Their result shows that the flow field is influenced appreciably by chemical reaction, heat source and suction or injection at the wall of the wedge. Recently Makinde

*et al.*(2011) and Makinde and Olanrewaju (2011) studied the problems of unsteady convection with chemical reaction and radiative heat transfer past a flat porous plate moving through a binary mixture in an optically thin environment. More recently Maleque (2013b) investigated an unsteady natural convection boundary layer flow with mass transfer and a binary chemical reaction and also studied the effects of exothermic chemical reaction. In the present paper, we investigate a numerical solution of unsteady natural convection heat and mass transfer boundary layer incompressible fluid flow past a uniform moving flat porous plate taking into account the effect of Arrhenius activation energy and heat generation /absorption. The rate of chemical reaction is also considered the function of temperature. This problem is an extension work studied by Maleque (2013c).

**II. GOVERNING EQUATIOS**

We consider the boundary wall to be of infinite extent so that all quantities are homogeneous in *x *and hence all derivatives with respect to *x* are neglected. Thus the governing equations are

(2) | |

(3) | |

(4) | |

(5) |

The boundary conditions of above system are

(6) |

where (*u*,*v*) is the velocity vector, *T *is the temperature, *C* is the concentration of the fluid, *υ* is the kinemetic coefficient of viscosity, *β* and *β*^{*} are the coefficients of volume expansions for temperature and concentration respectively, *κ* is the heat diffusivity coefficient, *c _{p}* is the specific heat at constant pressure,

*D*is the coefficient of mass diffusivity,

_{m}*k*

_{r}^{2 }is the chemical reaction rate constant and (

*T*-

*T*

_{∞})

*exp[-*

^{w}*E*/

_{a}*k*(

*T*-

*T*

_{∞})] is the Arrhenius function where

*w*is a unit less constant exponent fitted rate constants typically lie in the range -1<

*w*<1 (http://www.iupac.org/goldbook/mo3963.pdf).

**III. MATHEMATICAL FORMULATIONS**

In order to solve the governing Eqs. (2) to (5) under the boundary conditions (6), we adopt the well defined similarity technique to obtain the similarity solutions. For this purpose the following non-dimensional variables are now introduced:

(7) |

From the equation of continuity (3) we have

(8) |

where *v*_{0} is the dimensionless suction/injection velocity at the plate, *v*_{0}>0 corresponds to suction and *v*_{0}<0 corresponds to injection.

Introducing the dimensionless quantities from Eq. (8) and *v* from Eq. (9) in Eqs. (4), (5) and (6), we finally obtain the nonlinear ordinary differential equations as

(9) | |

(10) | |

(11) |

Here, Grashof number *G _{r}=gβ*(

*T*-

_{w}*T*

_{∞})δ

^{2}/(

*υU*

_{0}), Modified (Solutal) Grashof number

*G*(

_{m}=gβ^{*}*C*-

_{w}*C*

_{∞})

*δ*

^{2}/(

*υU*

_{0}), Prandtl Number

*P*, Schmidt number

_{r}=ρυc_{p}/κ*S*, The chemical reaction rate constant λ

_{c}=υ/D_{M}^{2}=

*k*

_{r}^{2}

*δ*

^{2}(

*T*-

_{w}*T*

_{∞})

*/*

^{w}*υ*, dimensionless activation energy

*E*=

*E*/[

_{a}*k*(

*T*-

_{w}*T*

_{∞})]. The dimensionless heat generation / absorption coefficient

*γ*=

*δ*

^{2}

*Q*

_{0}/(

*υρc*). The Eqs. (9) to (11) are similar except for the term

_{p}*δδ '*/

*υ*, where time

*t*appears explicitly. Thus the similarity condition requires that

*δδ*'/

*υ*must be a constant quantity. Hence following Maleque (2010a) one can try a class of solutions of the Eqs. (10) to (12) by assuming that

(12) |

From Eq. (10) we have

(13) |

where the constant of integration *L* is determined through the condition that *δ*=*L *when *t*=0 . Here *A*=0 implies that *δ*=*L *represents the length scale for steady flow and *A*≠0 that is, *δ* represents the length scale for unsteady flow. Let us now consider a class of solutions for which *A*=2 and hence the length scale *δ* from Eq. (12) becomes

(14) |

which exactly corresponds to the usual scaling factor considered for various unsteady boundary layer flow (Schlichting, 1968). Since *δ* is a scaling factor as well as similarity parameter any other values of *A* in Eq. (14) would not be change the nature of the solution except that the scale would be different. Finally introducing Eq. (14) in Eqs. (9) to (11) respectively, we have the following dimensionless nonlinear ordinary differential equations

(15) | |

(16) | |

(17) |

The boundary conditions Eq. (6) then becomes

(18) |

In all over equations primes denote the differentiation with respect to *η*. Equations (15) to (17) are solved numerically under the boundary conditions (18) using Nachtsheim-Swigert iteration technique.

**IV. NUMERICAL SOLUTIONS**

Equations (15) -(17) are solved numerically under the boundary conditions (18) using Nachtsheim and Swigert (1965) iteration technique. Thus adopting this numerical technique, a computer program was set up for the solutions of the basic non-linear differential equations of our problem where the integration technique was adopted as the six ordered Range-Kutta method of integration. The results of this integration are then displayed graphically in the form of velocity, temperature and concentration profiles in Figs. (1)-(15). In the process of integration, the local skin-friction coefficient, the local rates of heat and mass transfer to the surface, which are of chief physical interest are also calculated out. The equations defining the wall skin-friction, Nusselt number and Sherwood number are

(19) | |

(20) |

and

(21) |

Here, the Reynolds number *R _{e}*=

*ρU*

_{0}

*δ/μ.*

Fig. 1: Effect of *G _{r}* on the velocity profiles.

Fig. 2: Effect of *G _{m}* on the velocity profiles.

Fig. 3: Effects of *λ* on the concentration profiles.

Fig. 4: Effects of *λ* on the temperature profiles.

Fig. 5: Effects of *λ* on the velocity profiles.

Fig. 6: Effects of *E* on the concentration profiles.

Fig. 7: Effects of *E* on the temperature profiles.

Fig. 8: Effects of *E* on the velocity profiles.

Fig. 9: Effects of *v*_{0} on the velocity profiles.

Fig. 10: Effects of *v*_{0} on the temperature profiles.

Fig. 11: Effects of *v*_{0} on the concentration profiles.

Fig. 12: Effects of *γ* on the velocity profiles.

Fig. 13: Effects of *γ* on the temperature profiles.

Fig. 14: Effects of *γ* on the concentration profiles.

Fig. 15: Comparison of our calculated velocity profile and the velocity profile of Bestman (1990).

These above coefficients are then obtained from the procedure of the numerical computations and are sorted in Table 1.

Table 1: Effects of *G _{r}*,

*λ*,

*w*,

*E*and

*v*

_{0}on the skin friction coefficient, Nusselt number and the Sherwood number for

*G*=10,

_{r}*G*=1,

_{m}*P*=0.7,

_{r}*v*

_{0}=0.6,

*λ*=5,

*E*=1 and

*γ*=0.

**V. RESULTS AND DISCUSSIONS**

The parameters entering into the fluid flow are Grashof number *G _{r}*, Solutal (modified) Grashof number

*G*, suction parameter

_{m}*v*

_{0}, Prandtl number

*P*, the non dimensional chemical reaction rate constant λ

_{r}^{2}, Schmidt number

*S*and the non-dimensional activation energy

_{c}_{ }*E*.

It is, therefore, pertinent to inquire the effects of variation of each of them when the others are kept constant. The numerical results are thus presented in the form of velocity profiles, temperature profiles and concentration profiles in Figs.1-14 for the different values of *G _{r}*,

*G*, λ,

_{m }*E*,

*v*

_{0}and γ. The value of

*G*is taken to be both positive and negative, since these values represent respectively cooling and heating of the plate. The values of

_{r}*G*is taken to be large (

_{r }*G*=10), since the value corresponds to a cooling problem that is generally encountered in nuclear engineering in connection with the cooling of reactors. In air (

_{r }*P*= 0.71) the diluting chemical species of most common interest have Schmidt number in the range from 0.6 to .75. Therefore, the Schmidt number

_{r}*S*=0.6 is considered. In particulars, 0.6 corresponds to water vapor that represents a diffusing chemical species of most common interest in air. The values of the suction parameter

_{c}*v*

_{0}are taken to be large. Apart from the above Figures and Tables, the representative velocity, temperature and concentration profiles and the values of the physically important parameters, i.e. the local shear stress, the local rates of heat and mass transfers, are illustrated for uniform wall temperature and species concentration in Figs. 1-14 and in Table 1.

The velocity profiles generated due to impulsive motion of the plate are plotted in Fig.1 for both cooling (*G _{r}*>0) and heating (

*G*<0) of the plate keeping other parameters fixed

_{r}*G*=1.0,

_{m}*P*=0.71,

_{r}*S*=0.6,

_{c}*λ*=5,

*v*

_{0}=3.0,

*E*=1.0 and

*w*=1. In Fig.1, velocity profiles are shown for different values of

*G*. We observe that velocity increases with increasing values of

_{r}*G*for the cooling of the plate. From this figure it is also observed that the negative increase in the Grashof number leads to the decrease in the velocity field. That is, for heating of the plate (Fig.1), the effects of the Grashof number

_{r }*G*on the velocity field have also opposite effects, as compared to the cooling of the plate. Solutal Grashof number

_{r}*G*>0 corresponds that the chemical species concentration in the free stream region is less than the concentration at the boundary surface. Figure 2 presents the effects of Solutal Grashof number

_{m}*G*on the velocity profiles. It is observed that the velocity profile increases with the increasing values of Solutal Grashof number

_{m }*G*

_{m}. Considering chemical reaction rate constant λ^{2} is always positive. Figures 3-5 represent the effect of chemical rate constant λ on the concentration, the temperature and the velocity profiles respectively. We observe from the last part of Eq. (18) that λ^{2}exp(-*E*/θ) increases with the increasing values of λ. We also observe from this equation that increase in λ^{2}exp(-*E*/θ) means increase in λ leads to the decrease in the concentration profiles. This is in great agreement with Fig.3. It is observed from the Eq. (1) that increasing temperature frequently causes a marked increase in the rate of reactions is shown in Fig. 4. The parameter λ does not enter directly into the momentum equation but its influence comes through the mass equation. Figure5 shows the variation of the velocity profiles for different values of λ. The velocity profile increases with the increasing values of λ.

Activation energy may be defined as the minimum energy required starting a chemical reaction. The activation energy of a reaction is usually denoted by *Ea*, and given in units of kilojoules per mole. Effects of Activation Energy (*E*) on the concentration, the temperature and velocity profiles are shown in fig.6 to Fig.8 respectively. From Eq. (1) we observe that chemical reaction rate (*K*) decreases with the increasing values of activation energy (*Ea*). We also observe from Eq. (18) that increase in activation energy (*E*) leads to decrease λ^{2}exp(-*E*/θ) as well as to increase the concentration profiles shown in Fig. 6. Equation (1) suggests that the activation energy is dependent on temperature. For fixed value of chemical reaction rate constant the temperature profile increases with the increasing values of activation energy shown in Fig.7. The parameter *E* does not enter directly into the momentum equation but its influence comes through the mass equation. Figure 8 shows the variation of the velocity profiles for different values of *E*. From this figure it has been observed that the velocity profile increases with the increasing values of *E.*

The effects of suction and injection (*v*_{0}) for λ= 5, *E*=1, *G _{r}*= 5,

*G*=1,

_{m}*S*=0.6 and

_{c}*P*=0.71 on the velocity profiles, temperature profiles and concentration profiles are shown respectively in Fig. 9 to Fig. 11. For strong suction (

_{r}*v*

_{0}>0), the velocity, the temperature and the concentration profiles decay rapidly away from the surface. The fact that suction stabilizes the boundary layer is also apparent from these figures. As for the injection (

*v*

_{0}<0) , from Fig. 9 to Fig. 11 it is observed that the boundary layer is increasingly blown away from the plate to form an interlayer between the injection and the outer flow regions.

The positive value and the negative value of γ represent the heat generation and heat absorption respectively. The dimensionless heat generation/ absorption coefficient γ does not enter directly into the equation of continuity, Navier Stokes equations and concentration equation, its influence comes through the energy equation. In the presence of the buoyancy force, it appears that the heat generation/ absorption coefficient γ has marked effect on the boundary layer and the velocity profiles of the fluid is shown in Fig. 12. From this figure it has been observed that the effect of heat generation coefficient is to expand the boundary layer thickness and opposite effect is found for the effect of heat absorption. Thus the dimensionless heat generation/absorption coefficient has same effects on skin-friction coefficients shown in Table 1. The effects of heat absorption coefficient (γ<0) on temperature profiles and the Nusselt number are shown in Fig. 13 and Table 1. In general, the effect of heat absorption coefficient (γ<0) is to reduce the thickness of the thermal boundary layer. The temperature profiles decrease as the heat absorption coefficient γ (<0) increases as shown in Fig. 13. On the other hand, the heat generation coefficient γ (>0) leads to increase the temperature profiles. Thermal boundary layer thickness expands for increasing values of the heat generation coefficient γ (>0) as shown in Fig. 13

Finally the effects of heat generation/absorption coefficient (γ) on the dimensionless Nusselt number *N _{u}* is shown in Table 1. The Nusselt number increases with the increasing values of the heat absorption coefficient (γ<0). This is due to the fact that the heat absorption increases, the disk tends to supply more energy through its wall to maintain the constant temperature of the wall, and hence the Nusselt number increases. Opposite effects are found for the increasing values of the heat generation coefficient (γ> 0).

The dimensionless heat generation/absorption coefficient *γ* does not enter directly into the concentration equation, its influence comes through the energy equation. But in the presence of the activation energy, it appears that the heat generation/absorption coefficient *γ* has negligible effect on the concentration profiles of the fluid. The concentration profiles decrease for increasing values of *γ* as shown in Fig. 14. Also it has been found from Table 1 that Sherwood number *S _{h} *increases with the increasing values of

*γ.*

The skin-friction, the heat transfer and the mass transfer coefficients are tabulated in Table 1 for different values of *G _{r}*, λ,

*E*,

*w*,

*v*

_{0}and γ. We observe from Table 1 that the skin-friction coefficient (-

*f*′(0)) decreases with increase in the Grashop number

*G*. It is also observed from Table-1 that the rate of mass transfer coefficient increases with the increasing values of

_{r}*λ*and the mass transfer coefficient decreases with the increasing values of

*E*and

*w*. It is observed from this table that this is an increasing effect of suction parameter

*v*

_{0}on the skin friction coefficient, heat and mass transfer coefficients. Finally the effects of heat generation/ absorption coefficient (

*γ*) on the dimensionless Nusselt number

*N*is shown in Table 1. The Nusselt number increases with the increasing values of the heat absorption coefficient (

_{u}*γ*<0). This is due to the fact that the heat absorption increases, the disk tends to supply more energy through its wall to maintain the constant temperature of the wall, and hence the Nusselt number increases. Opposite effects are found for the increasing values of the heat generation coefficient (γ>0).

Comparison is made between our numerical result and the analytical result of Bestman (1990). Bestman (1990) studied steady natural convective boundary layer flow with large suction. He solved his problem analytically by employing the perturbation technique proposed by Singh and Diskshit (1988). In our present work, taking *A *= 0 in Eq. (13) for considering steady flow. The values of the suction parameter *v*_{0}^{2} = 10 is taken to see the effects of large suction. *G _{r}=G_{m}=*1.0, λ=5,

*E*=5 and

*w*=1 are also chosen with a view to compare our numerical results with the analytical results of Bestman (1990). The comparison of velocity profiles as seen in Fig. 15 highlight the validity of the numerical computations adapted in the present investigation.

**VI. CONCLUSIONS**

In this paper, we investigate the effects of chemical reaction rate and Arhenious activation energy on an unsteady natural convection heat and mass transfer boundary layer flow past a flat porous plate. The Nachtsheim and Swigert (1965) iteration technique based on sixth-order Range-Kutta and Shooting method has been employed to complete the integration of the resulting solutions.

The following conclusions can be drawn as a result of the computations:

- Velocity increases with increasing values of
*G*for the cooling of the plate and the negative increase in the Grashof number leads to the decrease in the velocity field. That is, for heating of the plate the effects of the Grashof number_{r}_{ }*G*on the velocity field have also opposite effects, as compared to the cooling of the plate._{r } - Solutal Grashof number
*G*>0 corresponds that the chemical species concentration in the free stream region is less than the concentration at the boundary surface. It is observed that the velocity profile increases with the increasing values of Solutal Grashof number_{m}*G*._{m} - Increase in λ leads to the decrease in the concentration profiles. It is observed from the Eq. (1) that increasing temperature frequently causes a marked increase in the rate of reactions
*λ*. The velocity profile slightly increases with the increasing values of*λ*. - Increase in activation energy (
*E*) leads to increase the concentration, temperature and velocity profiles. - For strong suction (
*v*_{0}>0), the velocity, the temperature and the concentration profiles decay rapidly away from the surface. As for the injection (*v*_{0}<0), it is observed that the boundary layer is increasingly blown away from the plate to form an interlayer between the injection and the outer flow regions. - The Nusselt number increases with the increasing values of the heat absorption coefficient (γ<0). This is due to the fact that the heat absorption increases, the disk tends to supply more energy through its wall to maintain the constant temperature of the wall, and hence the Nusselt number increases.

**NOMENCLATURE**

(*u*, *v*): components of velocity field.

*c _{p} *: specific heat at constant pressure

*N _{u} *: Nusselt number

*P _{r}*: Prandtl number

*Q*_{0}: Heat generation/absorption coefficient

*S _{h} *: Sherwood number

*T* : Temperature of the flow field

*T _{∞}* : Temperature of the fluid at infinity

*T _{w}*: Temperature at the plate

*f *: Dimensionless similarity functions

*E _{a}*: The activation energy

*k* : The Boltzmann constant

*S _{c}*: Schmidt number

*E : *The non-dimensional activation energy

*G _{r}*: Grashof number

*G _{m}* : Modified (Solutal) Grashof number

**GREEK SYMBOLS**

*θ* : Dimensionless temperature

: Dimensionless concentration

*η* : Dimensionless similarity variable

*δ* : Scale factor

*υ* : Kinematic viscosity

*ρ* : Density of the fluid

*κ* : Thermal conductivity,

*τ* : Shear stress

*μ* : Fluid viscosity

*λ*^{2 }: Non-dimensional chemical reaction rate constant

*β* and *β*^{*} : the coefficients of volume expansions for temperature and concentration respectively,

γ : The dimensionless heat generation/ absorption coefficient

**REFERENCES**

1. Alabraba, M.A., A.R. Bestman and A. Ogmu, "Laminar convection in binary mixture of hydromagnetic flow with radiative heat transfer," *Astrophysics and Space Science*, **195**, 431-439 (1992). [ Links ]

2. Al-Sharif, A., K. Chamniprasart, K.R. Rajagopal and A.Z. Szeri, "Lubrication with binary mixtures: Liquid-liquid emulsion," *J. Tribol*., **115**, 46-55 (1993). [ Links ]

3. Beevers, C.E. and R.E. Craine, "On the determination of response function for a binary mixture of incompressible Newtonian fluid," *Int. J. Engg. Sci.*, **20**, 737-745 (1982). [ Links ]

4. Bestman, A.R., "Natural convection boundary layerwith and mass transfer in a porous medium," *International Journal of Energy Research*, **14**, 389-396 (1990). [ Links ]

5. Bestman, A.R., "Radiative heat transfer to flow of a combustible mixture in a vertical pipe," *International Journal of Energy Research*, **15**, 179-184 (1991). [ Links ]

6. Blums, E., Y.A. Mikhailov and R. Ozols, *Heat and Mass Transfer in MHD Flow*, World Scientific Publishing Co. Pts. Ltd. (1987). [ Links ]

7. Chaudhary, R.C., and B.K. Sharma, "Combined heat and mass transfer by laminar mixed convection flow from a vertical surface with induced magnetic field," *J. Appl. Phys.*, **99**, 034901 (2006). [ Links ]

8. Ezzat, M., A.A. El-Bary and S. Ezzat, "Combined heat and mass transfer for unsteady Mhd flow of perfect conducting micro polar fluid with thermal relaxation," *Energy Conversion and Management*, **52**, 934-945 (2011). [ Links ]

9. Gebhart, B., Y. Jaluria, R.L. Mahajan and B. Sammakia, *Buoancy-induce flow and transport*, Washington, Hemisphere (1988). [ Links ]

10. Hossain, M.A. and A.S. Rees, "Combined heat and mass transfer in natural convection flow from a vertical wavy surface," *Acta Mechanica*, **136**, 133-141 (1999). [ Links ]

11. Kandasamy, R., K. Periasami and K.K. Sivagnana Probhu, "Effects of chemical reaction, heat and mass transfer along a wedge with heat source and concentration in the presence of suction and injection," *International Journal of Heat and Mass Transfer*, **48**, 1388-1394 (2005). [ Links ]

12. Khair, K.R. and A. Bejan, "Mass transfer to natural convection boundary flow driven by heat transfer," *Int. J. Heat Mass Transfer*, **30**, 369-376 (1985). [ Links ]

13. Lin, H.T. and C.M. Wu, "Combined heat and mass transfer by laminar natural convection from a vertical plate with uniform heat flux and concentrations," *Heat Mass Transfer*, **32**, 293-299 (1997). [ Links ]

14. Maleque, Kh.A, "Magnetohydrodynamic convective heat and mass transfer due to a rotating disk with thermal diffusion effect," *ASME Journal of Heat Transfer*, **131**, 082001 (2009). [ Links ]

15. Maleque, Kh.A.,"Dufour and Soret effects on unsteady MHD convective heat and mass transfer flow due to a rotating disk," *Latin American Applied Research*, **40**, 105-111 (2010a). [ Links ]

16. Maleque, Kh.A., "Effects of combined temperature-and depth-dependent viscosity and hall current on an unsteady Mhd laminar convective flow due to a rotating disk," *Chemical Engineering Communications*, **197**, 506-521 (2010b). [ Links ]

17. Maleque, Kh.A.,"Soret effect on steady Mhd convective heat and mass transfer incompressible fluid flow due to a porous rotating disk," *The Aiub Journal of Science and Engineering(AJSE),* **10**, 105-116 (2011). [ Links ]

18. Maleque, Kh.A., "Effects of Exothermic/Endothermic Chemical Reactions with Arrhenius Activation Energy on MHD Free Convection and Mass Transfer Flow in Presence of Thermal Radiation," J*ournal of Thermodynamics*, **2013**, ID 692516 (2013a). [ Links ]

19. Maleque, Kh.A., "Unsteady Natural Convection Boundary Layer Flow with Mass Transfer and a Binary Chemical Reaction," *British Journal of Applied Science and Technology,* **3**, 131-149 (2013b). [ Links ]

20. Maleque, Kh.A., "Unsteady natural convection boundary layer heat and mass transfer flow with exothermic chemical reactions," *Journal of Pure and Applied Mathematics: Advances and applications*, **9**, 17-41 (2013c). [ Links ]

21. Maleque, Kh.A. and M.A. Sattar, "Similarity solution of Mhd free convective and mass transfer flow over an infinite vertical porous plate," *The Aiub J. Science and Engineering*, **1**, 44-55 (2002). [ Links ]

22. Makinde, O.D., P.O. Olanrewaju and W.M. Charles, "Unsteady convection with chemical reaction and radiative heat transfer past a flat porous plate moving through a binary mixture," *Africka Matematika AM,* **21**, 1-17 (2011). [ Links ]

23. Makinde, O.D. and P.O. Olanrewaju, "Unsteady mixed convection with Soret and Dufour effects past a porous plate moving through a binary mixture of chemically reacting fluid,"* Chemical Engineering Communications*, **198**, 920-938 (2011). [ Links ]

24. Mills, N, "Incompressible mixtures of Newtonian fluids,"* Int. J. Engg. Sci*., **4**, 97-112 (1966). [ Links ]

25. Mongruel, A., M. Cloitre and C. Allain, "Scaling of boundary layer flows driven by double-diffusive convection," *Int. J. Heat Mass Transfer*, **39**, 3899-3910 (1996). [ Links ]

26. Nachtsheim, P.R. and P. Swigert, *Satisfaction of asymptotic boundary conditions in numerical solution of system of nonlinear of boundary layer type,* NASA TN-D3004 (1965). [ Links ]

27. Sattar, M.A. and Kh.A. Maleque, "Numerical solution of Mhd free convection and mass transfer flow with hall current, viscous dissipation, and joule heating," *J. Energy, Heat and Mass Transfer*, **22**, 237-244 (2000). [ Links ]

28. Schlichting, H., *Boundary Layer Theory*, McGraw Hill Book Co., New York (1968). [ Links ]

29. Singh, A.K and C.K. Dikshit, "Hydromagnetic flow past a continuously moving semi-infinite plate for large suction," *Astrophysics and Space Science*, **14**, 219-256 (1988). [ Links ]

30. Tencer, M., J.S. Moss and T. Zapach "Arrhenius Average Temperature: The Effective Temperature for Non-Fatigue Wearout and Long Term Reliability in Variable Thermal Conditions and Climates," *IEEE Transactions on Components and Packaging Technologies*, **27**, 602-607 (2004). [ Links ]

31. Truesdell, C, "Sulle basi della thermomeccanica," *Rend. Lince*, **22**, 33-38, 158-166 (1957). [ Links ]

32. Wang, S.R., A. Al-Sharif, K.R. Rajagopal and A.Z. Szeri, "Lubrication with binary mixtures: Liquid-liquid emulsion an EHL conjunction," *J. Tribol.*, 115, 515-522 (1993). [ Links ]

**Received: November 6, 2012 Accepted: July 4, 2013 Recommended by Subject Editor: Walter Ambrosini**