Research Article Effect of Suction/Injection on Unsteady Hydromagnetic Convective Flow of Reactive Viscous Fluid between Vertical Porous Plates with Thermal Diffusion I. J. Uwanta and M. M. Hamza Department of Mathematics, Usmanu Danfodiyo University, PMB 2346, Sokoto, Nigeria Correspondence should be addressed to M. M. Hamza; [email protected] Received 11 April 2014; Revised 23 June 2014; Accepted 7 July 2014; Published 4 November 2014 Academic Editor: Majid Ghassemi Copyright Β© 2014 I. J. Uwanta and M. M. Hamza. 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. An investigation is performed to study the effect of suction/injection on unsteady hydromagnetic natural convection flow of viscous reactive fluid between two vertical porous plates in the presence of thermal diffusion. The partial differential equations governing the flow have been solved numerically using semi-implicit finite-difference scheme. For steady case, analytical solutions have been derived using perturbation series method. Suction/injection is used to control the fluid flow in the channel, and an exothermic chemical reaction of Arrhenius kinetic is considered. Numerical results are presented graphically and discussed quantitatively with respect to various parameters embedded in the problem.

1. Introduction Suction or injection on the boundary layer control played significant role in the field of aerodynamics and space sciences. Shojaefard et al. [1] used suction/injection to control fluid flow on the surface of subsonic aircraft. By controlling the flow as such, fuel consumption might be decreased by 30%, a considerable reduction in pollutant emission is achieved, and operating costs of commercial airplanes are reduced by at least 8%; see the study by Braslow in [2]. In mass transfer cooling, suction or injection of a fluid through the bounding surface can significantly change the flow field and, as a result, affect the heat transfer rate from the plate; see the study by Ishak et al. in [3]. Many interests have been built in the study of flow of heat and mass transfer with suction or injection because of its extensive engineering applications. In the area of steady flow of viscous incompressible fluid over infinite porous plates subject to suction or injection, various aspects of the problem have been investigated by many authors. To be more specific, Griffith and Meredith [4] investigated the steady flow of an incompressible viscous fluid over an infinite porous flat plate subject to uniform suction. Jena and Mathur [5] studied free convection in the laminar boundary layer flow of a thermomicropolar fluid over a vertical flat plate subject

to uniform suction or injection. Boundary layer controls by suction or injection in the flow of incompressible fluid over an infinite porous wedge are to be found in the studies by Devi and Kandasamy [6] and Kandasamy et al. [7, 8]. Layek et al. [9], Shateyi [10], and Cortell [11] have analyzed the stretching sheet problem with suction or injection. Attia [12] reported the unsteady flow due to a rotating disk with uniform suction or injection. Al-Sanea [13] investigated mixed convection heat transfer along a continuously moving heated vertical plate with suction or injection. Unsteady free convection and mass transfer flow over an infinite vertical porous plate considering suction or injection are to be found in the study by Takhar et al. in [14]. Recently, Cortell [15] studied the effects of suction, viscous dissipation, and thermal radiation on flow and heat transfer of a power-law fluid past an infinite porous plate. Effect of suction and injection on unsteady free convection Couette flow and heat transfer of reactive viscous fluid in vertical porous plate is to be found in the study by Jha et al. in [16]. The investigation of the flow of an electrically conducting fluid in a porous channel in the presence of a transverse magnetic field is important because of its widespread engineering and industrial applications such as MHD marine propulsion, electronic packages, microelectronic devices,

2 thermal insulation, petroleum reservoirs, MHD stirring of molten metal, exothermic reaction in packaged reactors, and magnetic-levitation casting. On the other hand, Soret or thermal diffusion is important where more than one chemical species are present under very large temperature gradients, such as chemical reactions and isotope separation, and in mixtures of gases with very light molecular weight such as hydrogen or helium and of medium molecular weight such as nitrogen or air. Because of the applications of MHD and Soret, many authors investigated their effects on natural convection heat and mass transfer flow. Postelnicu [17] reported the influence of magnetic field on heat and mass transfer by natural convection from vertical surfaces in porous media considering Soret and Dufour effects. Osalusi et al. [18] determined numerically the effects of Soret and Dufour on heat and mass transfer of a steady MHD convective and slip flow due to a rotating disk with viscous dissipation and ohmic heating. Recently, Turkyilmazoglu and Pop [19] reported the effects of Soret and heat source on unsteady radiative MHD free convection flow from an impulsively infinite vertical plate. In nutshell, there have been considerable published works dealing with steady flow with Soret and Dufour effect; some of them are the works of Alam et al. [20], Chamkha and Ben-Nakhi [21], Tsai and Huang [22], Tak et al. [23], and Magyari and Postelnicu [24]. Steady flows with chemical reaction considering Soret and Dufour effect are to be found in the studies by Mansour et al. [25], Beg et al. [26], El-Kabeir et al. [27], and Gangadhar [28]. Unsteady fluid flow problems in the presence of Soret and Dufour effects with chemical reaction can be found in the studies by Bhargava et al. [29] and Pal and Mondal [30]. Nandkeolyar et al. [31] investigated numerically and analytically the effect of suction/injection on unsteady hydromagnetic heat and mass transfer flow of a radiating and chemically reactive fluid past a flat porous plate with ramped wall temperature. Numerical investigation of buoyancy effects on hydromagnetic unsteady flow through a porous channel considering suction and injection is to be found in the study by Makinde and Chinyoka [32]. The aim of the present analysis is to study the effect of suction/injection on unsteady hydromagnetic convective flow of viscous reactive fluid between two infinite vertical parallel porous plates in the presence of transverse magnetic field and thermal diffusion. In this paper, an exothermic chemical reaction of Arrhenius kinetics is employed and suction/injection is used to control fluid flow in the channel.

2. Governing Equations Consider the transient natural convection and mass transfer flow of viscous reactive, incompressible, and electrically conducting fluid between infinite vertical parallel porous plates under the influence of a transversely magnetic field of strength π΅0 ; see Figure 1. The magnetic Reynolds number is assumed to be small so that the induced magnetic field and the hall effect of MHD are negligible. At time π‘σΈ β€ 0, both the fluid and the plates are at rest and at the same temperature and concentration π0σΈ and πΆ0σΈ , respectively. At π‘σΈ > 0, the temperature and concentration of the plate π¦σΈ = 0 are raised to ππσΈ and πΆπσΈ and thereafter remain constant

International Scholarly Research Notices xσ³° B0

yσ³°

g

Cσ³° = Cπσ³°

Cσ³° = C0σ³°

Tσ³° = Tπσ³°

Tσ³° = T0σ³°

uσ³° = 0 Suction

Injection

yσ³° = 0

yσ³° = H

Figure 1: Schematic diagram of the problem.

and those of π¦σΈ = π» are lowered to π0σΈ and πΆ0σΈ , where ππσΈ > π0σΈ and πΆπσΈ > πΆ0σΈ . It is assumed that the flow is subjected to suction of the fluid from one porous plate and at the same rate fluid is being injected through the other porous plate. We chose a Cartesian coordinate system with the π₯σΈ axis along the upward direction and the π¦σΈ axis normal to it. The physical properties are assumed to be constant excluding density in the buoyancy term. The fluid is assumed to be Newtonian and obeys the Boussinesqβs approximation. Under the previous assumptions, the momentum, energy, and concentration equations in the dimensional form are the following: ππ’σΈ π2 π’ σΈ ππ’σΈ β V = ] 0 ππ‘σΈ ππ¦σΈ ππ¦σΈ 2 + ππ½ (πσΈ β π0σΈ ) + ππ½β (πσΈ β π0σΈ ) β

ππ΅02 σΈ π’, π

ππσΈ ππσΈ π π2 πσΈ ππΆ0 π΄ (βπΈ/π πσΈ ) β V0 σΈ = + π , σΈ ππ‘ ππ¦ ππΆπ ππ¦σΈ 2 ππΆπ ππΆσΈ π2 πΆσΈ π· π π2 πσΈ ππΆσΈ β V0 σΈ = π·π σΈ 2 + π π σΈ 2 . σΈ ππ‘ ππ¦ ππ¦ ππ ππ¦

(1)

The initial and boundary conditions for the present problem are the following: π‘σΈ β€ 0 : π’σΈ = 0,

πσΈ σ³¨β π0σΈ ,

π‘σΈ > 0 : π’σΈ = 0, π’σΈ = 0,

πσΈ = ππσΈ ,

πσΈ = π0σΈ ,

πΆσΈ σ³¨β πΆ0σΈ , πΆσΈ = πΆπσΈ

πΆσΈ = πΆ0σΈ

0 β€ π¦σΈ β€ π», at π¦σΈ = 0,

as π¦σΈ σ³¨β π», (2)

where π is the conductivity of the fluid, π΅0 is the electromagnetic induction, π½ is the coefficient of thermal expansion, π½β is the coefficient of concentration expansion, π is the heat of reaction, π΄ is the rate constant, πΈ is the activation energy, π is the universal gas constant, ] is the kinematic viscosity, πΆ0σΈ is the initial concentration of the reactant species, π is the gravitational force, πΆπ is the specific heat at constant pressure, π is the thermal conductivity of the fluid, π is the density of

International Scholarly Research Notices

3

the fluid, π·π is the coefficient of mass diffusivity, ππ is the mean fluid temperature, and ππ is the thermal diffusion ratio. In order to solve (1) to (2), we employ the following dimensionless parameters:

π¦= Pr = π=

π¦σΈ , π»

π‘=

]ππΆπ π

,

π‘σΈ ] , π»2

Gc =

πΈ (πσΈ β π0 ) π π02

,

π’=

π=

π=

π π0 , πΈ

πΎ=

]0 π» , ]

ππ =

πΈ (ππ β π0 ) , π π02

π=

Gr =

ππ½π π02 π»2 , πΈ]]0

ππΆ0 π΄πΈπ»2 (βπΈ/π π0 ) π , π π02 πΆ=

πΈ (πΆσΈ β πΆ0 ) π πΆ02 πΆπ =

,

πΈ (πΆπ β πΆ0 ) , π πΆ02

ππ΅02 π»2 . ]π

ππ ππ π2 π + ππ(π/(1+ππ)) , βπΎ )= ππ‘ ππ¦ ππ¦2

(4)

π’ = 0,

π = 0,

πΆ = πΆπ

πΆ = 0,

at π¦ = 0,

π = 1,

πΆ = 1,

at π¦ = 0,

π’ = 0,

π = 0,

πΆ = 0,

at π¦ = 1.

(7)

In order to construct an approximate solution to (6) subject to (7), we employed a regular perturbation method by taking a power series expansion in the Frank-Kamenetskii parameter π as follows: π’ = π’0 + ππ’1 , (8)

π’ = π»1 πβ1 π¦ + π»2 πββ2 π¦ + π»3

+ π [π»6 πβ1 π¦ + π»7 πββ2 π¦ + π»8 π¦ + π»9 + (π»10 + π»11 ) πβπΎPrπ¦ + π»12 π¦πβπΎPrπ¦ + π»13 πβ2πΎPrπ¦ + π»14 πβπΎScπ¦ ] , π = 1 + π΅ (πβπΎPrπ¦ β 1)

0 β€ π¦ β€ 1, π‘ β€ 0,

π = ππ ,

π’ = 0,

+ π»4 πβπΎPrπ¦ + π»5 πβπΎScπ¦

The initial and boundary conditions in dimensionless form are the following:

π‘ > 0 : π’ = 0,

ππΆ π2 πΆ π2 π + Sr . = ππ¦ ππ¦2 ππ¦2

Substituting (8) into (6), the solutions of the governing equations are obtained as follows:

ππΆ π2 π ππΆ π2 πΆ + Sr 2 . βπΎ )= 2 ππ‘ ππ¦ ππ¦ ππ¦

π = 0,

(6)

π = π0 + ππ1 .

ππ’ ππ’ π2 π’ + Grπ + GcπΆ β ππ’, βπΎ = ππ‘ ππ¦ ππ¦2

π’ = 0,

ππ π2 π + ππ(π/(1+ππ)) , = ππ¦ ππ¦2

πΆ = πΆ0 + ππΆ1 ,

Using (3), (1) to (2) can take the following form:

Sc (

βπΎSc

ππ’ π2 π’ = + Grπ + GcπΆ β ππ’, ππ¦ ππ¦2

The boundary conditions are the following:

(3)

Pr (

βπΎ βπΎPr

ππ½β π πΆ02 π»2 , πΈ]]0

π π2 Sr = π 02 , ππ πΆ0

] , Sc = π·π

π’σΈ , ]0

ππ’/ππ‘ = 0, ππ/ππ‘ = 0, and ππΆ/ππ‘ = 0 into (4) to (5) and by taking ππ = 1 and πΆπ = 1 at the boundary, we get

(5)

as π¦ = 1.

3. Analytical Solutions The analytical solutions have played an important role in validating and exploring computer routines of complicated problems. They are also used to inspect the internal consistency of mathematical models and of the approximations adopted by Jha et al. in [16]. Therefore, we reduce the governing equations of this problem due to its nonlinearity into a form that can be solved analytically. By setting

(9)

+ π [π·1 + π·2 πβπΎPrπ¦ + π·3 π¦ + π·4 π¦πβπΎPrπ¦ + π·5 πβ2πΎPrπ¦ ] , πΆ = πΈ1 + πΈ2 πβπΎScπ¦ + πΈ3 πβπΎPrπ¦ + π [πΈ4 + πΈ5 πβπΎScπ¦ + (πΈ6 + πΈ7 ) πβπΎPrπ¦ + πΈ8 π¦πβπΎPrπ¦ + πΈ9 πβ2πΎPrπ¦ ] . Using (9), we write the steady-state skin friction, rate of heat transfer, and rate of mass transfer on the boundaries as follows.

4

International Scholarly Research Notices

Steady-state skin frictions on the boundary plates are the following: ππ’ σ΅¨σ΅¨σ΅¨σ΅¨ = π»1 β1 β π»2 β2 β πΎPrπ»4 β πΎScπ»5 σ΅¨ ππ¦ σ΅¨σ΅¨σ΅¨π¦=0 + π [π»6 β1 β π»7 β2 + π»8 β (π»10 + π»11 + 2π»13 ) + π»12 β π»14 ScπΎ] , ππ’ σ΅¨σ΅¨σ΅¨σ΅¨ = π»1 β1 πβ1 β π»2 β2 πββ2 β πΎPrπ»4 πβπΎPr σ΅¨ ππ¦ σ΅¨σ΅¨σ΅¨π¦=1

(π) = π2 π’πβ1 + (1 β 2π2 β πΎπ3 β πΞπ‘) π’π(π)

β πΎScπ»5 π

ββ2

+ π [π»6 β1 π β π»7 β2 π

(π) + (π2 + πΎπ3 ) π’π+1 + Ξπ‘Grπππ + Ξπ‘GcπΆππ,

+ π»8 βπΎPr

+ (π»12 + π»12 πΎPr β π»11 πΎPr β π»10 πΎPr) π

(10) The steady-state rate of heat transfer on the boundary plates is ππ σ΅¨σ΅¨σ΅¨σ΅¨ = β (βπ΅πΎPr + π [βπ·2 πΎPr + π·3 + π·4 β 2π·5 πΎPr]) , β σ΅¨ ππ¦ σ΅¨σ΅¨σ΅¨π¦=0 ππ σ΅¨σ΅¨σ΅¨σ΅¨ = βπ΅πΎPrπβπΎPr σ΅¨ ππ¦ σ΅¨σ΅¨σ΅¨π¦=1 + π [π·3 + (π·4 β π·4 πΎPr β π·2 πΎPr) πβπΎPr β2πΎPrπβ2πΎPr ] . (11) The steady-state rate of mass transfer on the boundary plates is given as follows: ππΆ σ΅¨σ΅¨σ΅¨σ΅¨ σ΅¨ ππ¦ σ΅¨σ΅¨σ΅¨π¦=0 = β (βπΎScπΈ2 β πΎPrπΈ3 +π [βπΈ5 πΎSc β πΎPr (πΈ6 + πΈ7 + 2πΈ9 ) + πΈ8 ]) , ππΆ σ΅¨σ΅¨σ΅¨σ΅¨ σ΅¨ ππ¦ σ΅¨σ΅¨σ΅¨π¦=1

(π+1) (π+1) + (Sc + 2π1 ) πΆπ(π+1) β π1 πΆπ+1 β π1 πΆπβ1 (π) = π2 πΆπβ1 + (Sc β 2π2 β ScπΎπ3 ) πΆπ(π)

β 2π»13 πΎPrπβ2πΎPr β π»14 ScπΎπβπΎSc ] .

β

The complete forms of (4) are solved numerically using semi-implicit finite-difference method given in the study by Makinde and Chinyoka in [32]. We used forward difference formulas for all time derivatives and approximate the spatial derivatives with central difference formula. The semi-implicit finite-difference equations corresponding to (4) are as follows: (π+1) (π+1) + (1 + 2π1 ) π’π(π+1) β π1 π’π+1 β π1 π’πβ1

βπΎSc β1

4. Numerical Solutions

(12)

= βπΎScπΈ2 πβπΎSc β πΎPrπΈ3 πβπΎPr + π [βπΎScπΈ5 πβπΎSc + (πΈ8 β πΈ8 πΎPr β πΈ6 πΎPr β πΈ7 πΎPr) πβπΎPr β2πΎPrπΈ9 πβ2πΎPr ] . The constants π΅, π»1 , π»2 , π»3 , π»4 , π»5 , π»6 , π»7 , π»8 , π»9 , π»10 , π»11 , π»12 , π»13 , π»14 , π·1 , π·2 , π·3 , π·4 , π·5 , πΈ1 , πΈ2 , πΈ3 , πΈ4 , πΈ5 , πΈ6 , πΈ7 , πΈ8 , πΈ9 , β1 , and β2 are defined in the appendix section.

(13)

(π) π π + (π2 + ScπΎπ3 ) πΆπ+1 + π4 (ππβ1 β 2πππ + ππ+1 ), (π+1) (π+1) + (Pr + 2π1 ) ππ(π+1) β π1 ππ+1 β π1 ππβ1 (π) = π2 ππβ1 + (Pr β 2π2 β πΎPrπ3 ) ππ(π) (π) + (π2 + πΎPrπ3 ) ππ+1 + πΞπ‘ exp (

ππ(π) 1 + πππ(π)

),

where π1 = πΞπ‘/Ξπ¦2 , π2 = (1 β π)Ξπ‘/Ξπ¦2 , π3 = Ξπ‘/Ξπ¦, π4 = SrΞπ‘/Ξπ¦2 , and 0 β€ π β€ 1. We chose π = 1; the detailed reasons to this particular selection are documented in [32]. Also the analytical solutions displayed in the previous section are used as a check on the accuracy and effectiveness of the numerical scheme. Again, in order to reconfirm the accuracy of the scheme, the numerical results for velocity, concentration, and temperature are compared with the analytical solutions. It has been found that the numerical values of the velocity, concentration, and temperature fields calculated from the expressions (9) have matched very well with the numerical obtained from the expressions (13) at the steady-state time. See Figure 2 for the graph of the numerical solutions at steady state and steady-state analytical solutions for velocity, concentration, and temperature fields.

5. Results and Discussion The numerical results are obtained by solving (13) using the method described in the previous section for various values of physical parameters to describe the physics of the problem. The nondimensional parameters that govern the flow are the Prandtl number (Pr), which is inversely proportional to the thermal diffusivity of the working fluid, the Frank-Kamenetskii parameter (π), the Soret number (Sr), the magnetic parameter (π), the thermal Grashof number (Gr), the solutal Grashof number (Gc), the nondimensional time (π‘), the Schmidt number (Sc), which is inversely proportional to the mass diffusivity of the working fluid, and

Velocity/temperature/concentration

International Scholarly Research Notices

5

1.2 1 0.8 0.6 0.4 0.2 0

0

0.1

0.2

0.3

0.4

0.5 y

0.6

0.7

0.8

0.9

1

Num. soln. for velocity at steady state Steady-state analy. soln. for velocity Num. soln. for temp. at steady state Steady-state analy. soln. for temp. Num. soln. for concentration at steady state Steady-state analy. soln. for concentration

Figure 2: Unsteady and steady-state solutions for velocity, concentration, and temperature profiles.

suction/injection parameter (πΎ), which were simultaneously applied each to opposite porous plates of the channel at the same rate. For the purpose of discussion, some numerical calculations are carried out for dimensionless velocity (π’), temperature (π), concentration (πΆ), skin friction, rate of heat transfer in terms of Nusselt number, and the rate of mass transfer in terms of Sherwood number. Unless otherwise stated, the values π = 0.1, Gr = 0.1, π = 1, Gc = 0.1, Sr = 0.1, ππ = 1, πΆπ = 1, Pr = 0.71, Sc = 0.62, π‘ = 0.1, πΎ = 0.5, and π = 0.01 are used for the investigation. Results obtained are displayed graphically for velocity, temperature, concentration, skin friction, Nusselt number, and Sherwood number for various flow parameters. Figure 3 shows the effects of the Frank-Kamenetskii parameter (π) and suction/injection (πΎ) on the temperature profiles. From Figure 3(a), it is observed that temperature of the fluid increases with increasing values of π in case of suction and injection, respectively. This is physically true since an increase in π leads to significant increases in the reaction and viscous source terms and hence significantly increases the fluid temperature. It is evident from Figure 3(a) that temperature of the fluid is greater in case of injection than suction. Figure 3(b) represents the influence of suction/injection parameter on the temperature field. From Figure 3(b), it is seen that temperature decreases due to suction but increases due to injection. In case of suction, the fluid at ambient conditions is brought closer to the surface and reduces the thermal boundary layer thickness. The same principle operates but in reverse direction in case of injection. Figures 4(a) and 4(b) illustrated the effects of the FrankKamenetskii parameter (π) and suction/injection parameter (πΎ) on velocity distribution, respectively. Figure 4(a) revealed that increasing π accelerates the velocity of the fluid in case of suction and injection. It is evident from this figure that the velocity is higher in case of injection than suction. From Figure 4(b), it is seen that velocity of the fluid decelerates due

to suction while it accelerates due to blowing. The physical explanation for such a behavior is that while stronger blowing is provided, the heated fluid is pushed farther from the wall where the buoyancy forces can act to accelerate the flow with less influence of the viscosity. This effect acts to increase the shear by increasing the maximum velocity within the boundary layer. The same principle operates but in reverse direction in case of suction. It is also noticed that, in case of suction, velocity of the fluid moves away from the channel centerline towards the plate (π¦ = 0) and, in case of injection, the maximum velocities are shifted towards the right porous plate (π¦ = 1). In Figure 5, we have presented the response of the fluid velocity to variations in the Soret number (Sr) and magnetic parameter (π) in the presence of suction and injection parameter. Figure 5(a) shows the effect of the Soret number on the velocity with constant suction and injection. From this figure, it is noted that Soret number accelerates the fluid velocity in the presence of suction and injection. In Figure 5(a), it is seen that, in the presence of suction, the velocity of the fluid moves toward the left porous plate and in case of injection the maximum velocities move towards the right porous plate. Figure 5(b) revealed that increasing the strength of magnetic parameter is to decrease velocity profiles. This is due to the fact that transverse magnetic field produces a resistivity force (Lorentz force) similar to the drag force which retards the velocity. It is seen that the velocity of the fluid is greater in case of injection than suction. The effects of the Soret number (Sr) and suction/injection parameter (πΎ) on the concentration distribution are shown graphically in Figures 6(a) and 6(b), respectively. From Figure 6(a), it can be noticed that concentration of the fluid increases with increasing values of Sr in case of suction and injection. Further, Figure 6(a) reveals that the concentration of the fluid is high in the vicinity of the wall (π¦ = 0), where suction takes place more than at the wall (π¦ = 1) where injection takes place. From Figure 6(b), it can be noticed that concentration of the fluid decreases due to suction but increases due to injection. The influence of the thermal Grashof number (Gr) and solutal Grashof number (Gc) is illustrated in Figures 7(a) and 7(b), respectively. These plots of Figures 7(a) and 7(b) indicate that the momentum boundary layer thickness increases with increasing values of Gr and Gc. It is further noticed from these figures that velocity of the fluid is greater in case of injection than suction. The wall shear stress dependence on reaction parameter π is illustrated in Figure 8 for varying values of the nondimensional time when suction and injection take place. Figures 8(a) and 8(b) represent the wall shear stress at the walls π¦ = 0 and π¦ = 1, respectively. From Figure 8, it is observed that skin friction increases as time and π increase until a steady-state condition is reached in case of suction and injection. Figure 8(a) reveals that skin friction is greater in case of suction than injection, but opposite trend is noticed in Figure 8(b). The rate of heat transfer in terms of Nusselt number dependence on the reaction parameter π is displayed in Figure 9 for varying values of time in the presence of suction and injection. Figures 9(a) and 9(b) represent the Nusselt number at π¦ = 0 and π¦ = 1 for different values

International Scholarly Research Notices 1.4

1.4

1.2

1.2

1

1

0.8

0.8

Temperature

Temperature

6

0.6

0.6

0.4

0.4

0.2

0.2

0

0

0.2

0.4

0.6

0.8

0

1

y

0

0.2

0.4

0.6

0.8

1

y

π = 0.1, πΎ = 5 π = 0.5 π=1 π=2 π=3

π π π π π

= = = = =

β0.1, πΎ = β5 β0.5 β1 β2 β3

πΎ πΎ πΎ πΎ πΎ

= = = = =

4 6 8 10 0

(a)

πΎ = β4 πΎ = β6 πΎ = β8 πΎ = β10

(b)

Figure 3: Effects of Frank-Kamenetskii parameter (π) and suction/injection (πΎ) on temperature. Γ10β3 9

0.01

8

0.009

7

0.008 0.007 0.006

5

Velocity

Velocity

6

4

0.005 0.004

3

0.003

2

0.002

1 0

0.001 0

0.2

0.4

0.6

0.8

y π = 0.1, πΎ = 5 π = 0.5 π=1 π=2 π=3

0

0

0.2

0.4

0.6

0.8

y π π π π π

(a)

1

= = = = =

β0.1, πΎ = β5 β0.5 β1 β2 β3

πΎ πΎ πΎ πΎ πΎ

= = = = =

4 6 8 10 0

πΎ = β4 πΎ = β6 πΎ = β8 πΎ = β10

(b)

Figure 4: Effects of Frank-Kamenetskii parameter (π) and suction/injection (πΎ) on velocity.

1

International Scholarly Research Notices

7 Γ10β3 8

0.01 0.009

7

0.008 6 0.007 5 Velocity

Velocity

0.006 0.005 0.004

4 3

0.003 2 0.002 1

0.001 0

0

0.2

0.4

0.6

0.8

0

1

0

0.2

0.4

y Sr Sr Sr Sr Sr

Sr = β0.1, πΎ = β5 Sr = β0.5 Sr = β1 Sr = β2 Sr = β3

= 0.1, πΎ = 5 = 0.5 =1 =2 =3

0.6

0.8

1

y M = β1, πΎ = β5 M = β5 M = β10 M = β15 M = β20

M = 1, πΎ = 5 M=5 M = 10 M = 15 M = 20

(a)

(b)

Figure 5: Effects of Soret number (Sr) and magnetic parameter (π) on velocity. 1.5

1 0.9 0.8 0.7 Concentration

Concentration

1

0.5

0.6 0.5 0.4 0.3 0.2 0.1

0

0

0.2

0.4

0.6

0.8

y Sr Sr Sr Sr Sr

= = = = =

0.1, πΎ = 5 0.5 1 2 3

(a)

1

0

0

0.2

0.4

0.6 y

Sr = β0.1, πΎ = β5 Sr = β0.5 Sr = β1 Sr = β2 Sr = β3

πΎ πΎ πΎ πΎ πΎ

= = = = =

4 6 8 10 0

πΎ = β4 πΎ = β6 πΎ = β8 πΎ = β10

(b)

Figure 6: Effects of Soret number (Sr) and suction/injection (πΎ) on concentration.

0.8

1

8

International Scholarly Research Notices

0.035

0.06

0.03

0.05

0.025

Velocity

Velocity

0.04

0.03

0.02 0.015

0.02 0.01 0.01

0

0.005

0.2

0

0.4

0.6

0.8

0

1

0

0.2

0.4

0.6

y Gr Gr Gr Gr

= 0.1, πΎ = 5 = 0.5 =1 = 1.5

0.8

1

y Gr = β0.1, πΎ = β5 Gr = β0.5 Gr = β1 Gr = β1.5

Gc = β0.2, πΎ = β5 Gc = β0.4 Gc = β0.6 Gc = β0.8

Gc = 0.2, πΎ = 5 Gc = 0.4 Gc = 0.6 Gc = 0.8

(a)

(b)

Figure 7: Effects of thermal Grashof number (Gr) and solutal Grashof number (Gc) on velocity.

0.018

0.03

0.016 0.014 0.025 Skin friction

Skin friction

0.012 0.01 0.008

0.02

0.006 0.004 0.015 0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.002 0.1

0.2

0.3

0.4

0.5

t = 0.1 suction t = 0.2 t = 0.3 Steady state

(a)

0.6

0.7

0.8

0.9

π

π t = β0.1 injection t = β0.2 t = β0.3 Steady state

t = 0.1 suction t = 0.2 t = 0.3 Steady state

t = β0.1 injection t = β0.2 t = β0.3 Steady state

(b)

Figure 8: Variation of skin friction with π and time (π‘) at π¦ = 0 and π¦ = 1.

1

International Scholarly Research Notices

9

0.7

0.9

0.6

0.8 0.7

0.4

Nusselt number

Nusselt number

0.5

0.3 0.2

0.5 0.4

0.1

0.3

0

0.2

β0.1 β0.2 0.1

0.6

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.1 0.1

0.2

0.3

0.4

0.5

t = 0.1 suction t = 0.2 t = 0.3 Steady state

0.6

0.7

0.8

0.9

1

π

π t = β0.1 injection t = β0.2 t = β0.3 Steady state

t = 0.1 suction t = 0.2 t = 0.3 Steady state

(a)

t = β0.1 injection t = β0.2 t = β0.3 Steady state

(b)

Figure 9: Variation of Nusselt number with π and time (π‘) at π¦ = 0 and π¦ = 1.

of π, respectively. From Figure 9(a), it is seen that the rate of heat transfer decreases by increasing time and π until a steady-state condition is attained in case of suction and injection. It is also evident from Figure 9(a) that Nusselt number decreases more in case of injection than suction. Figure 9(b) reveals that rate of heat transfer at the plate π¦ = 1 increases as time and π increase until a steady-state value is achieved. From this figure, it is seen that Nusselt number is higher in case of injection than suction. The skin friction on suction/injection parameter for varying values of the Soret number is illustrated in Figure 10. Figures 10(a) and 10(b) represent the skin friction at the walls π¦ = 0 and π¦ = 1, respectively. From Figure 10(a), it is noticed that skin friction increases as Sr increases in case of suction while it decreases due to injection. Figure 10(b) reveals that skin friction increases with increasing values of Sr in case of injection while it decreases due to suction. Figure 11 shows the influence of Sr on the Sherwood number. Figures 11(a) and 11(b) represent the Sherwood number at the plates π¦ = 0 and π¦ = 1. Sherwood number decreases as Sr increases in case of suction and injection; see Figure 11(a). Figure 11(b) reveals that Sherwood number increases as Sr increases in case of injection while it decreases as Sr increases due to suction. Figure 12 shows the wall shear stress dependence on magnetic parameter π at the plates π¦ = 0 and π¦ = 1 for varying values of nondimensional time in case of suction and injection. Figures 12(a) and 12(b) represent the skin friction at the plates π¦ = 0 and π¦ = 1, respectively. Figure 12 reflected that skin friction increases with increasing time

until a steady-state condition is attained in case of suction and injection. Figure 12(a) revealed that higher values of π reduce the skin friction in case of suction and injection. Reduction is noticed in case of injection more than suction. From Figure 12(b), it is seen that as π increases the skin friction decreases in the presence of suction and injection. Further, it is observed from Figure 12(b) that reduction is seen in the presence of suction more than injection when π is large. The wall shear stress dependence on Gr for varying values of time at the porous plates π¦ = 0 and π¦ = 1 when suction and injection take place is illustrated in Figure 13. Figures 13(a) and 13(b) represent the skin friction at the walls π¦ = 0 and π¦ = 1, respectively. Figure 14 shows the wall shear stress dependence on Gc for varying values of time when suction and injection take place. Figures 14(a) and 14(b) represent the skin friction at the porous plates π¦ = 0 and π¦ = 1, respectively. Figures 13 and 14 revealed that skin friction increases by increasing nondimensional time until a steady-state value is achieved. In both Figures 13(a) and 14(a), skin friction when suction takes place is greater than when injection takes place, but reverse effect is observed in Figures 13(b) and 14(b).

6. Conclusion In the present study, the effect of suction/injection on transient hydromagnetic convective flow of viscous reactive fluid between vertical porous plates in the presence of transverse magnetic field and thermal diffusion is investigated. It is

10

International Scholarly Research Notices Γ10β3 7

0.019 0.0185

6

0.018

Skin friction

Skin friction

5 0.0175 0.017 0.0165

4

3

0.016 2

0.0155 0.015 0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1 0.1

1

0.2

0.3

0.4

0.5

Sr Sr Sr Sr

0.6

0.7

0.8

0.9

1

πΎ

πΎ Sr = 0.1 injection Sr = 0.5 Sr = 1 Sr = 1.5

= 0.1 suction = 0.5 =1 = 1.5

Sr Sr Sr Sr

Sr = 0.1 injection Sr = 0.5 Sr = 1 Sr = 1.5

= 0.1 suction = 0.5 =1 = 1.5

(a)

(b)

0.7

0.8

0.6

0.7

0.5

0.6 Sherwood number

Sherwood number

Figure 10: Variation of skin friction with Sr and suction/injection (πΎ) at π¦ = 0 and π¦ = 1.

0.4 0.3

0.5 0.4

0.2

0.3

0.1

0.2

0 0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.1 0.1

0.2

0.3

0.4

0.5

Sr Sr Sr Sr

= 0.1 suction = 0.5 =1 = 1.5

(a)

0.6

0.7

0.8

0.9

πΎ

πΎ Sr = 0.1 injection Sr = 0.5 Sr = 1 Sr = 1.5

Sr Sr Sr Sr

= 0.1 suction = 0.5 =1 = 1.5

Sr = 0.1 injection Sr = 0.5 Sr = 1 Sr = 1.5

(b)

Figure 11: Variation of Sherwood number with Sr and suction/injection (πΎ) at π¦ = 0 and π¦ = 1.

1

11

0.028

Γ10β3 16

0.026

14

0.024

12

0.022

Skin friction

Skin friction

International Scholarly Research Notices

0.02

10 8

0.018

6

0.016

4

0.014

1

1.5

2

2.5

3

3.5

4

4.5

2

5

1

1.5

2

2.5

3

3.5

4

4.5

5

M

M t = 0.1 injection t = 0.2 t = 0.3 Steady state

t = 0.1 suction t = 0.2 t = 0.3 Steady state

t = 0.1 injection t = 0.2 t = 0.3 Steady state

t = 0.1 suction t = 0.2 t = 0.3 Steady state

(a)

(b)

Figure 12: Variation of skin friction with magnetic parameter (π) at π¦ = 0 and π¦ = 1.

0.16

0.09

0.14

0.08 0.07 0.06

0.1 Skin friction

Skin friction

0.12

0.08 0.06

0.04 0.03

0.04

0.02

0.02 0 0.1

0.05

0.01

0.2

0.3

0.4

0.5 0.6 Gr

t = 0.1 suction t = 0.2 t = 0.3 Steady state

(a)

0.7

0.8

0.9

t = 0.1 injection t = 0.2 t = 0.3 Steady state

1

0 0.1

0.2

0.3

0.4

0.5 0.6 Gr

t = 0.1 suction t = 0.2 t = 0.3 Steady state

0.7

0.8

0.9

t = 0.1 injection t = 0.2 t = 0.3 Steady state

(b)

Figure 13: Variation of skin friction with thermal Grashof number (Gr) at π¦ = 0 and π¦ = 1.

1

12

International Scholarly Research Notices 0.11

0.06

0.1 0.05

0.09

0.04

0.07

Skin friction

Skin friction

0.08

0.06 0.05

0.03

0.02

0.04 0.03

0.01

0.02 0.01 0.1

0.2

0.3

0.4 Gc

0.5

0.6

0.7

0 0.1

t = 0.1 injection t = 0.2 t = 0.3 Steady state

t = 0.1 suction t = 0.2 t = 0.3 Steady state

0.2

0.3

0.4 Gc

0.6

0.7

t = 0.1 injection t = 0.2 t = 0.3 Steady state

t = 0.1 suction t = 0.2 t = 0.3 Steady state

(a)

0.5

(b)

Figure 14: Variation of skin friction with solutal Grashof number (Gc) at π¦ = 0 and π¦ = 1.

found that suction/injection, thermal diffusion, reaction consumption, and thermal and solutal buoyancy play an important role in controlling the transport phenomena. Formation of the minimum flow occurs near the wall where suction takes place, while the maximum flow forms near the wall where injection takes place except in the concentration distribution for varying values of thermal diffusion. It is hoped that the present work may be useful in engineering applications where the formation of boundary layer is to be delayed or enhanced.

Consider the following: β1 , πβπΎPr β 1

π·3 = β π·5 =

π=

πβ1 , 2

(2 β π΅ β π(1 β π΅)2 ) πΎPr 2

π΅π , πΎ2 Pr2

π·2 = β

πΈ6 = β

Srπ·4 πΎ2 Pr2 , β ScπΎ2 Pr

πΈ7 = β

πΈ8 = β πΈ9 = β

β1 =

Appendix

π΅ =

SrπΎ2 Pr2 π΅ , πΈ3 = β 2 2 πΎ Pr β ScπΎ2 Pr

π·4 = β

(2 (1 β π΅) π΅π β π΅) , πΎPr

π·1 = βπ·2 β π·5 ,

π·5 (πβ2πΎPr β 1) (πβπΎPr β 1)

πΈ1 = 1 β πΈ2 β πΈ3 ,

,

β

πΈ2 = β

1 (πβπΎSc

β 1)

β

πΈ3 (πβ2πΎPr β 1) (πβπΎSc

β 1)

,

4Srπ·5 πΎ2 Pr2 , 4πΎ2 Pr2 β 2ScπΎ2 Pr

βπΎ + βπΎ2 + 4π 2

πΈ8 (πΎSc β 2πΎPr) , πΎ2 Pr2 β ScπΎ2 Pr

πΈ4 = βπΈ5 β πΈ6 β πΈ7 β πΈ9 ,

β2 =

,

πΎ2 Pr2 β ScπΎ2 Pr

πΎ + βπΎ2 + 4π 2

,

β (Grπ΅ + GcπΈ3 ) , πΎ2 Pr2 β πΎ2 Pr β π

π»3 =

Gr (1 β π΅) + GcπΈ1 , π

π»4 =

π»5 =

βGcπΈ2 , πΎ2 Sc2 β πΎ2 Sc β π

π»1 = βπ»2 β π»3 β π»4 β π»5 ,

π»2 = π»8 =

π· π·4 πβπΎPr β βπΎPr 3 , βπΎPr (π β 1) (π β 1)

πΎ2 Pr2

Sr (π·2 πΎ2 Pr2 β 2πΎPrπ·4 )

π»3 (πβ1 β 1) (πββ2 β πβ1 ) Grπ·3 , π

+

π»4 (πβ1 β πβπΎPr )

π»9 =

(πββ2 β πβ1 )

+

π»5 (πβ1 β πβπΎSc ) (πββ2 β πβ1 )

Grπ·1 + GcπΈ4 + πΎπ»8 , π

π»10 =

β (Grπ·2 + GcπΈ6 + GcπΈ7 ) , πΎ2 Pr2 β πΎ2 Pr β π

π»11 =

π»12 (πΎ + 2πΎPr) , 2 πΎ Pr2 β πΎ2 Pr β π

,

,

International Scholarly Research Notices π»12 = π»14 =

β (Grπ·4 + GcπΈ8 ) , πΎ2 Pr2 β πΎ2 Pr β π πΎ2 Sc2

π»13 =

13

β (Grπ·5 + GcπΈ9 ) , 4πΎ2 Pr2 β 2πΎ2 Pr β π

βGcπΈ5 , β πΎ2 Sc β π

π»6 = βπ»7 β π»9 β π»10 β π»11 β π»13 β π»14 , π»7 =

π»9 (πβ1 β 1) βπ»8 + (πββ2 β πβ1 ) (πββ2 β πβ1 ) + β +

π»10 (πβ1 β πβπΎPr ) (πββ2 β πβ1 )

+

π»11 (πβ1 β πβπΎPr ) (πββ2 β πβ1 )

π»12 πβπΎPr (πββ2 β πβ1 ) π»13 (πβ1 β πβ2πΎPr ) (πββ2 β πβ1 )

+

π»14 (πβ1 β πβπΎSc ) (πββ2 β πβ1 )

. (A.1)

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

Acknowledgment The author M. M. Hamza is thankful to Usmanu Danfodiyo University, Sokoto, Nigeria, for financial assistance.

References [1] M. H. Shojaefard, A. R. Noorpoor, A. Avanesians, and M. Ghaffapour, βNumerical investigation of flow control by suction and injection on a subsonic airfoil,β The American Journal of Applied Sciences, vol. 20, pp. 1474β1480, 2005. [2] A. I Braslow, A History of Suction Type Laminar Flow Control with Emphasis on Flight Research, American Institute of Aeronautics and Astronautics, Washington, Wash, USA, 1999. [3] A. Ishak, J. H. Merkin, R. Nazar, and I. Pop, βMixed convection boundary layer flow over a permeable vertical surface with prescribed wall heat flux,β ZAMP: Zeitschrift fΒ¨ur angewandte Mathematik und Physik, vol. 59, no. 1, pp. 100β123, 2008. [4] A. A. Griffith and F. W. Meredith, βThe possible improvement in Aircraft performance due to use of boundary layer suction,β Tech. Rep. 2315, Aeronautical Research Council, London, UK, 1936. [5] S. K. Jena and M. N. Mathur, βFree convection in the laminar boundary layer flow of a thermomicropolar fluid past a vertical flat plate with suction/injection,β Acta Mechanica, vol. 42, no. 3-4, pp. 227β238, 1982. [6] S. P. A. Devi and R. Kandasamy, βEffects of chemical reaction, heat and mass transfer on non-linear MHD laminar boundarylayer flow over a wedge with suction or injection,β International Communications in Heat and Mass Transfer, vol. 29, no. 5, pp. 707β716, 2002.

[7] R. Kandasamy, K. Periasamy, and K. K. S. Prabhu, βEffects of chemical reaction, heat and mass transfer along a wedge with heat source and concentration in the presence of suction or injection,β International Journal of Heat and Mass Transfer, vol. 48, no. 7, pp. 1388β1394, 2005. [8] R. Kandasamy, A. W. B. Raj, and A. B. Khamis, βEffects of chemical reaction, heat and mass transfer on boundary layer flow over a porous wedge with heat radiation in the presence of suction or injection,β Theoretical and Applied Mechanics, vol. 33, no. 2, pp. 123β148, 2006. [9] G. C. Layek, S. Mukhopadhyay, and S. A. Samad, βHeat and mass transfer analysis for boundary layer stagnation-point flow towards a heated porous stretching sheet with heat absorption/generation and suction/blowing,β International Communications in Heat and Mass Transfer, vol. 34, no. 3, pp. 347β356, 2007. [10] S. Shateyi, βThermal radiation and buoyancy effects on heat and mass transfer over a semi-infinite stretching surface with suction and blowing,β Journal of Applied Mathematics, vol. 2008, Article ID 414830, 12 pages, 2008. [11] R. Cortell, βFlow and heat transfer of a fluid through a porous medium over a stretching surface with internal heat generation/absorption and suction/blowing,β Fluid Dynamics Research, vol. 37, no. 4, pp. 231β245, 2005. [12] H. A. Attia, βOn the effectiveness of uniform suction-injection on the unsteady flow due to a rotating disk with heat transfer,β International Communications in Heat and Mass Transfer, vol. 29, no. 5, pp. 653β661, 2002. [13] S. A. Al-Sanea, βMixed convection heat transfer along a continuously moving heated vertical plate with suction or injection,β International Journal of Heat and Mass Transfer, vol. 47, no. 6-7, pp. 1445β1465, 2004. [14] H. S. Takhar, S. Roy, and G. Nath, βUnsteady free convection flow over an infinite vertical porous plate due to the combined effects of thermal and mass diffusion, magnetic field and Hall currents,β Heat and Mass Transfer, vol. 39, no. 10, pp. 825β834, 2003. [15] R. Cortell, βSuction, viscous dissipation and thermal radiation effects on the flow and heat transfer of a power-law fluid past an infinite porous plate,β Chemical Engineering Research and Design, vol. 89, no. 1, pp. 85β93, 2011. [16] B. K. Jha, A. K. Samaila, and A. O. Ajibade, βUnsteady/steady natural convection flow of reactive viscous fluid in a vertical annulus,β International Journal of Fluid Mechanics Research, vol. 39, no. 4, pp. 301β311, 2012. [17] A. Postelnicu, βInfluence of a magnetic field on heat and mass transfer by natural convection from vertical surfaces in porous media considering Soret and Dufour effects,β International Journal of Heat and Mass Transfer, vol. 47, no. 6-7, pp. 1467β1472, 2004. [18] E. Osalusi, J. Side, and R. Harris, βThermal-diffusion and diffusion-thermo effects on combined heat and mass transfer of a steady MHD convective and slip flow due to a rotating disk with viscous dissipation and Ohmic heating,β International Communications in Heat and Mass Transfer, vol. 35, no. 8, pp. 908β915, 2008. [19] M. Turkyilmazoglu and I. Pop, βSoret and heat source effects on the unsteady radiative MHD free convection flow from an impulsively started infinite vertical plate,β International Journal of Heat and Mass Transfer, vol. 55, no. 25-26, pp. 7635β7644, 2012.

14 [20] M. S. Alam, M. M. Rahman, A. Maleque, and M. Ferdows, βDufour and Soret effects on steady MHD combined freeforced convective and mass transfer flow past a semi-infinite vertical plate,β The Thammasat International Journal of Science and Technology, vol. 11, pp. 1β12, 2006. [21] A. J. Chamkha and A. Ben-Nakhi, βMHD mixed convectionradiation interaction along a permeable surface immersed in a porous medium in the presence of Soret and Dufourβs effects,β Heat and Mass Transfer, vol. 44, no. 7, pp. 845β856, 2008. [22] R. Tsai and J. S. Huang, βNumerical study of Soret and Dufour effects on heat and mass transfer from natural convection flow over a vertical porous medium with variable wall heat fluxes,β Computational Materials Science, vol. 47, no. 1, pp. 23β30, 2009. [23] S. S. Tak, R. Mathur, R. K. Gehlot, and A. Khan, βMHD free Convection-Radiation interaction along a vertical surface embedded in darcian porous medium in presence of soret and Dufourβs effects,β Thermal Science, vol. 14, no. 1, pp. 137β145, 2010. [24] E. Magyari and A. Postelnicu, βDouble-diffusive natural convection flows with thermosolutal symmetry in porous media in the presence of the Soret-Dufour effects,β Transport in Porous Media, vol. 88, no. 1, pp. 149β167, 2011. [25] M. A. Mansour, N. F. El-Anssary, and A. M. Aly, βEffects of chemical reaction and thermal stratification on MHD free convective heat and mass transfer over a vertical stretching surface embedded in a porous media considering Soret and Dufour numbers,β Chemical Engineering Journal, vol. 145, no. 2, pp. 340β345, 2008. [26] O. A. Beg, T. A. Beg, A. Y. Bakier, and V. R. Prasad, βChemicallyreacting mixed convective heat and mass transfer along inclined and vertical plates with soret and dufour effects: numerical solutions,β International Journal of Applied Mathematics and Mechanics, vol. 5, pp. 39β57, 2009. [27] S. M. M. El-Kabeir, A. Chamkha, and A. M. Rashad, βHeat and mass transfer by MHD stagnation-point flow of a powerlaw fluid towards a stretching surface with radiation, chemical reaction and Soret and Dufour effects,β International Journal of Chemical Reactor Engineering, vol. 8, pp. 1β18, 2010. [28] K. Gangadhar, βSoret and Dufour effects on hydro magnetic heat and mass transfer over a vertical plate with a convective surface boundary condition and chemical reaction,β Journal of Applied Fluid Mechanics, vol. 6, no. 1, pp. 95β105, 2013. [29] R. Bhargava, R. Sharma, and O. A. Beg, βOscillatory chemicallyreacting MHD free convection heat and mass transfer in a porous medium with Soret and Dufour effects: finite element modeling,β International Journal of Applied Mathematics and Mechanics, vol. 5, pp. 15β37, 2009. [30] D. Pal and H. Mondal, βEffects of Soret Dufour, chemical reaction and thermal radiation on MHD non-Darcy unsteady mixed convective heat and mass transfer over a stretching sheet,β Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 4, pp. 1942β1958, 2011. [31] R. Nandkeolyar, M. Das, and P. Sibanda, βUnsteady hydromagnetic heat and mass transfer flow of a heat radiating and chemically reactive fluid past a flat porous plate with ramped wall temperature,β Mathematical Problems in Engineering, vol. 2013, Article ID 381806, 12 pages, 2013. [32] O. D. Makinde and T. Chinyoka, βNumerical investigation of buoyancy effects on hydromagnetic unsteady flow through a porous channel with suction/injection,β Journal of Mechanical Science and Technology, vol. 27, no. 5, pp. 1557β1568, 2013.

International Scholarly Research Notices