Journal of Fluids

Volume 2013 (2013), Article ID 297493, 14 pages

http://dx.doi.org/10.1155/2013/297493

## Free Convection Heat and Mass Transfer MHD Flow in a Vertical Porous Channel in the Presence of Chemical Reaction

^{1}Department of Mathematics, Trident Academy of Technology, Infocity, Bhubaneswar, Odisha 751024, India^{2}Department of Mathematics, S.O.A. University, Bhubaneswar, Odisha 751030, India^{3}Department of Mathematics, Christ College, Cuttack, Odisha, India

Received 20 May 2013; Revised 18 September 2013; Accepted 23 September 2013

Academic Editor: Hideki Tsuge

Copyright © 2013 R. N. Barik et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The objective of the present study is to examine the fully developed free convective MHD flow of an electrically conducting viscous incompressible fluid in a vertical porous channel under influence of asymmetric wall temperature and concentration in the presence of chemical reaction. The heat and mass transfer coupled with diffusion-thermo effect renders the present analysis interesting and curious. The analytical solution by Laplace transform technique of partial differential equations is used to obtain the expressions for the velocity, temperature, and concentration. It is observed that under the influence of dominating mass diffusivity over thermal diffusivity with stronger Lorentz force the velocity is reduced at all points Further, low rate of thermal diffusion delays the attainment of free stream state. Flow of aqueous solution in the presence of heavier species is prone to back flow.

#### 1. Introduction

In many transport processes and industrial applications, transfer of heat and mass occurs simultaneously as a result of combined buoyancy effects of thermal diffusion and diffusion of chemical species. Unsteady natural convection of heat and mass transfer is of great importance in designing control systems for modern free convection heat exchangers. More recently, Jha and Ajibade [1] have studied the heat and mass transfer aspect of the flow of a viscous incompressible fluid in a vertical channel considering the Dufour effect.

Soundalgekar and Akolkar [2] studied the effect of mass transfer and free convection currents on the flow past an impulsively started infinite vertical plate and observed that the presence of foreign gasses in the flow domain leads to reduce the shear stress and rate of mass transfer significantly.

In nature, flow occurs due to density differences caused by temperature as well as chemical composition gradients. Therefore, it warrants the simultaneous consideration of temperature difference as well as concentration difference when heat and mass transfer occurs simultaneously. It has been found that an energy flux can be created not only by temperature gradients but by composition gradients also. This is called Dufour effect. If, on the other hand, mass fluxes are created by temperature gradients, it is called the Soret effect.

The Soret and the Dufour effects have been found to be useful as the Soret effect is utilized for isotope separation and, in a mixture of gases of light and medium molecular weight, the Dufour effect is found to be of considerable order of magnitude such that it cannot be neglected. In view of importance of the diffusion-thermo effect Kafoussias and Williams [3] have studied the effects of thermal diffusion and diffusion-thermo on mixed free and forced convective and mass transfer boundary layer flow with temperature dependent viscosity. Jha and Singh [4], Kafoussias [5], and Alam et al. [6, 7] have contributed significantly to this field of study. Recently, Beg et al. [8] studied chemically reacting mixed convective heat and mass transfer along inclined and vertical plates considering Soret and Dufour effects. More recently, Dursunkaya and Worek [9] have studied the diffusion-thermo and thermal diffusion effects in transient and steady natural convections from a vertical surface.

MHD flow with thermal diffusion and chemical reaction finds numerous applications in various areas such as thermo nuclear fusion, liquid metal cooling of nuclear reactions, and electromagnetic casting of metals. Abreu et al. [10] have studied the boundary layer flows with Dufour and Soret effects. Osalusi et al. [11] have worked on mixed and free convective heat and mass transfer of an electrically conducting fluid considering Dufour and Soret effects. Anghel et al. [12] and Postelnicu [13] obtained numerical solutions of chemical reacting mixed convective heat and mass transfer along inclined and vertical plates with the Soret and the Dufour effects and concluded that skin friction increases with a positive increase in the concentration-to-thermal-buoyancy ratio parameter.

The analysis of natural convection heat and mass transfer near a moving vertical plate has received much attention in recent times due to its wide application in engineering and technological processes. There are applications of interest in which combined heat and mass transfer by natural convection occurs between a moving material and ambient medium, such as the design and operation of chemical processing equipment, design of heat exchangers, transpiration cooling of a surface, chemical vapour deposition of solid layer, nuclear reactor, and many manufacturing processes like hot rolling, hot extrusion, wire drawing, continuous casting, and fiber drawing. Gebhart and Pera [14] studied the effects of mass transfer on a steady free convection flow past a semi-infinite vertical plate by the similarity method and it was assumed that the concentration level of the diffusing species in fluid medium was very low. This assumption enabled them to neglect the diffusion-thermo and thermo-diffusion effects as well as the interfacial velocity at the wall due to species diffusion. Following this assumption, Das et al. [15] investigated the effects of simultaneous heat and mass transfer on free convection flow past an infinite vertical plate under different physical situations.

Recently, the unsteady MHD heat and mass transfer free convection flow of a polar fluid past a vertical moving porous plate in a porous medium with heat generation and thermal diffusion has been studied by Saxena and Dubey [16]. Raveendra Babu et al. [17] studied diffusion-thermo and radiation effects on MHD free convective heat and mass transfer flow past an infinite vertical plate in the presence of a chemical reaction of first order. Sudhakar et al. [18] discussed chemical reaction effect on an unsteady MHD free convection flow past an infinite vertical accelerated plate with constant heat flux, thermal diffusion, and diffusion thermo.

In the present paper, the flow of an electrically conducting viscous incompressible fluid in a vertical porous channel formed by two vertical parallel porous plates in the presence of a transverse magnetic field is studied. Further, the mass transfer phenomena considered in this problem is associated with chemically reacting species. The objective of the present study is to extend the work of Jha and Ajibade [1] by incorporating the permeability and the magnetic field effect as well as chemical reaction on the flow, heat and mass transfer phenomena.

#### 2. Mathematical Formulation

Let us consider the free convective heat and mass transfer MHD flow of a viscous incompressible fluid in a vertical porous channel formed by two infinite vertical parallel porous plates in the presence of chemical reaction. A magnetic field of uniform strength is applied transversely to the plate. The induced magnetic field is neglected as the magnetic Reynolds number of the flow is taken to be very small. The convection current is induced due to both the temperature and concentration differences. The flow is assumed to be in the direction which is taken to be vertically upward along the channel walls and *axis* is taken to be normal to the plates that are distance apart.

Under the usual Boussinesq’s approximations, the governing equations for flow are given by where , , , , , , , , , , , and are kinematic viscosity, acceleration due to gravity, coefficient of thermal expansion, coefficient of mass expansion, electrical conductivity of the fluid, uniform magnetic field, density of the fluid, chemical molecular diffusivity, the dimensional chemical reaction parameter in the diffusion equation, thermal diffusivity, the dimensional coefficient of the diffusion-thermo effect, and the permeability of the medium, respectively.

Initial and boundary conditions of the problem are We now introduce the following dimensionless quantities:

Here , , , , , , and are the Prandtl number, the Schmidt number, the Sustentation parameter, the Dufour number, the magnetic field parameter, the chemical reaction parameter, and the Porosity parameter, respectively.

The nondimensional form of (1) is given by Subject to boundary conditions,

#### 3. Solution of the Problem

Applying Laplace transform to (5) with boundary condition (6), we get

The boundary condition (6) becomes Applying inverse Laplace transform to (7) and (8) subject to boundary condition (9), following Debanath and Bhatta [19], and using the formula The solutions of (5) subject to the boundary conditions (6) are obtained by Laplace transform technique.

*Case 1 (, ). *The solutions of (5) subject to boundary conditions (6) are given by
where and .

The functions and erfc are given in Appendix A. The constants , , , , , and are given in Appendix B.

The rate of mass transfer (), the rate of heat transfer (), and the skin function () at the walls of the channel are obtained as follows:

*Sherwood Number (**).*
Consider

*Nusselt Number (**).*
Consider

*Skin Friction (**).*
Consider
where the function is given in Appendix A.

The mass flux (volumetric flow rate) for the problem is given by This is computed by using the trapezoidal rule for numerical integration.

*Case 2 (). *The solutions of (5) subject to boundary conditions (6) are given by
where and are given in Appendix B.

In this case, the rate of mass transfer (), the rate of heat transfer (), and the skin frictions () on the walls of the channel are obtained as follows:

*Sherwood Number.*
Consider

*Nusselt Number.*
Consider

*Skin Friction.*
Consider

*Steady State.* Setting , and in (5), the steady state of the problem is obtained as

The solutions of (21) subject to boundary conditions (6) are given by where is given in Appendix B.

In this case, the Sherwood number (), the Nusselt number (), and the skin friction () are given by the following:

*Sherwood Number (**)*. Consider

*Nusselt Number (**)*. Consider

*Skin Friction (**).*
Consider

*Mass Flux.*
Consider

#### 4. Results and Discussion

Computations have been carried out by assigning the values to the pertinent parameters characterising the fluids of practical interest. The flow phenomenon is characterized by magnetic parameter (), chemical reaction parameter (), sustentation parameter (), Dufour number (), Schmidt number (), Prandtl number (), and porosity parameter (). The effects of various parameters on velocity, temperature, and concentration profiles are shown graphically and in tabulated form for Sherwood number (), Nusselt number (), and skin friction ().

The particular case without magnetic field and porosity (, ) can be obtained from (8) which is in good agreement with the work reported earlier [1].

Another interesting point to note is that diffusion in aqueous solution, that is, for higher values of gives rise to oscillations in the velocity as well as temperature field in the presence of chemical reaction parameter ().

Figure 1 depicts the velocity profile for various values of pertinent parameters characterising flow fields. The common feature of the profiles is parabolic. This is evident from curves II and X (, Curve II; the steady state curve X) in case of steady state as well as ; the velocity profiles remain positive throughout the flow field. This is evident from curves II and X (, Curve II; the steady state curve X). Thus it is concluded that the time span plays a vital role for engendering back flow (Curve I and Curve II).

Further, magnetic parameter (), chemical reaction parameter (), sustentation parameter (), and porosity parameter () decrease the velocity at all points. In all other cases such as , , and , the reverse effect is observed. If the mass diffusivity becomes greater than the thermal diffusivity (i.e., , the sustentation parameter >1.0) then the velocity decreases. Moreover, with the stronger magnetic field Lorentz force also reduces the velocity field which is in conformity with the result of Cramer and Pai [20].

An increase in leads to decrease in the velocity which suggests that low rate of thermal diffusion leads to increase in the velocity boundary layer thickness.

Figure 2 depicts the velocity distribution in the absence of chemical reaction. It is noteworthy to observe that the back flow occurs for higher values of ; that is, and and the permeability of the medium representing the diffusion in aqueous solution in the presence of constant magnetic field and Dufour effect. Further, it is to mention that for and , no velocity profile could be graphed/presented. Further, it is to note that heavier species that is, with increasing lead to flow reversal. It is most interesting to note that in the presence of destructive reaction that is, , higher value of leads to oscillatory flow and, in case of diffusion in the aqueous solution, it leads to back flow. Thus, this suggests that the back flow is prevented by diffusing lighter species which is in conformity with the results reported earlier by Pop [21], Hossain and Mohammad [22], and Rath et al. [23].

It is also remarked that when varies from 0 to 5 the back flow is prevented in the absence of chemical reaction. But in the presence of heavier species the sharp fall of velocity is well marked.

Figure 3 exhibits the steady state velocity profiles indicating no back flow irrespective of higher or lower value of mass transfer coefficient and Dufour effect.

On careful observation, it is revealed that for higher values of Dufour number (curve VI), sustentation parameter (Curve VIII), magnetic parameter (Curve V), and porosity parameter (Curve IV) steady state velocity is effected significantly whereas variation in reaction parameter and Schmidt number produces no change. Thus, it is important to note that mass transfer phenomena with chemical reaction does not affect the steady flow significantly but increase in , , , and , decrease the velocity, increase the velocity, increase the velocity, and decrease the velocity, respectively. It may be inferred that magnetic parameter, porosity parameter, and Dufour effect produce the same effect irrespective of steady or unsteady state except for sustentation parameter.

Figures 4, 5, and 6 present the temperature distribution. It is seen that for an increase in destructive reaction parameter () temperature increases whereas increase in , , and reduces the temperature at all points for (water) when . It is noted that, the slow rate of diffusion leaves the heat energy spread in the fluid mass and it is enhanced for heavier species.

It is observed in case of (air) that the temperature falls sharply at all points with almost linear distribution regardless of other parameters. Thus, it may be pointed out that thermal diffusivity property of the fluid plays a vital role in controlling the temperature distribution and hence contributing to thermal boundary layer thickness.

Figure 5 also exhibits the temperature distribution when and which means that thermal diffusivity and mass diffusivity are at par. This contributes to the uniform variation/decrease in temperature. Increase in reaction parameter gives rise to increase in temperature whereas allowing for a large span of time reduces it.

Figure 6 presents the nonlinearity distribution of temperature due to high value of , that is, for large Dufour effect. Thus, sharp rise and fall of temperature is due to low diffusivity and high Dufour effect.

Figures 7, 8 and 9 exhibit the concentration variation in the flow domain. Sharp fall of concentration is indicated in Figures 7 and 9 along with mass absorption near the plate for very high value of (i.e., ). This means that heavier species give rise to sharp fall of concentration accelerating the process of mass diffusion whereas in case of lighter species the variation is smooth.

Skin friction plays an important role in flow characteristics. It measures the frictional forces encountered at the solid surfaces due to the motion of the fluid. One striking feature of the entries of the skin friction in Table 1, in case of unsteady motion, is that the skin frictions bear same sign in case of lower and upper plate. Further, on careful observation it is revealed that at the lower plate all are negative for except for the case when . Thus, it may be concluded that time span plays an important role to modify the frictional drag due to shear stress at the plates.

Moreover, from (18)-(19) it is clear that the Prandtl number () has no role to play to affect the velocity, temperature, and concentration fields in the steady flow and hence the skin friction. In case of unsteady flow, as increases, the shearing stress increases at both the plates.

Considering Tables 1 and 2 for both steady and unsteady cases, it is concluded that an increase in magnetic parameter, porosity parameter, and sustentation parameter leads to reduce the magnitude of frictional drag at both plates. Moreover, an increase in , , and leads to enhance the frictional drag at the plates for steady and unsteady flow.

Physically, we may interpret the above results as follows.

Lorentz force and sustentation parameter, that is, ratios of buoyancy effects due to temperature and concentration differences reduce the frictional drag which agrees well with the established result reported earlier. Moreover, presence of heavier species and exothermic reaction and Dufour effect enhance the magnitude of the shearing stress at both plates.

From Table 3 it is seen that an increase in , , and leads to decrease the Nusselt number at both plates. One interesting point is that rate of heat transfer remains negative for all the parameters at both the plates for aqueous solution but in case of air, that is, , Nusselt number is positive at the lower plate and negative at the upper plate.

Table 4 presents the variation of mass transfer for various values of , , and . It is seen that an increase in chemical reaction parameter increases the mass transfer at the lower plate and decreases it at the other. Further, it is to note that Sherwood number at both the plates increases as increases, that is, for heavier species mass transfer increases at the plates. The effect of increase in time span is to reduce the rate of mass transfer at both plates.

Tables 5 and 6 show the effects of various parameters on mass flux both for unsteady and steady cases, respectively. From Table 5, it is observed that mass flux increases with the increase in the values in , , , and but the reverse effect is observed in case of , , , and when the motion is unsteady. Hence, it is concluded that in case of unsteady motion more flux fluid was experienced with heavier species having low thermal diffusivity and increasing Dufour effect in the presence of constructive chemical reaction.

Now, if we analyse the steady case,we observe that more flux is measured with higher value of and whereas the reverse effect is observed in case of and . In case of and there is no significant effect on mass flux.

Comparing both cases steady and unsteady motion, it is observed that an increase in and gives rise to higher flux whereas magnetic field yields less amount of flux.

#### 5. Conclusion

(i)Under the influence of dominating mass diffusivity over thermal diffusivity with stronger Lorentz force the velocity is reduced at all points of the channel.(ii)Low rate of thermal diffusion leads to increase the thickness of boundary layer.(iii)Diffusion of heavier species leads to flow reversal. It is well marked in case of aqueous solution. Flow of aqueous solution in the presence of heavier species is prone to back flow.(iv)Destructive reaction in the presence of heavier species leads to oscillatory flow.(v)In case of steady state, no back flow occurs irrespective of high or low value of mass transfer coefficient and Dufour effects.(vi)Sharp rise and fall of temperature is the outcome of low diffusivity and high Dufour effects.(vii)Dufour effect and chemical reaction rate in the presence of heavier species enhance the frictional drag.(viii)Heat transfer bears same sign at both the plates in the presence of aqueous solution but in case of air it is of opposite sign.(ix)Mass transfer increases at both plates due to the presence of heavier species but it reduces with increasing time span.(x)In case of unsteady motion more flux of fluid is experienced due to heavier species with low thermal diffusivity. Dufour effect enhances the flux both in steady and unsteady motion.

#### Appendices

#### A.

The complementary error function is given by

#### B.

#### References

- B. K. Jha and A. O. Ajibade, “Free convection heat and mass transfer flow in a vertical channel with the Dufour effect,”
*Journal of Process Mechanical Engineering*, vol. 224, no. 2, pp. 91–101, 2010. View at Publisher · View at Google Scholar · View at Scopus - V. M. Soundalgekar and S. P. Akolkar, “Effects of free convection currents and mass transfer on flow past a vertical oscillating plate,”
*Astrophysics and Space Science*, vol. 89, no. 2, pp. 241–254, 1983. View at Publisher · View at Google Scholar · View at Scopus - N. G. Kafoussias and E. W. Williams, “Thermal-diffusion and diffusion-thermo effects on mixed free-forced 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 Google Scholar · View at Scopus - B. K. Jha and A. K. Singh, “Soret effects on free-convection and mass transfer flow in the stokes problem for a infinite vertical plate,”
*Astrophysics and Space Science*, vol. 173, no. 2, pp. 251–255, 1990. View at Publisher · View at Google Scholar · View at Scopus - N. G. Kafoussias, “MHD thermal-diffusion effects on free-convective and mass-transfer flow over an infinite vertical moving plate,”
*Astrophysics and Space Science*, vol. 192, no. 1, pp. 11–19, 1992. View at Publisher · View at Google Scholar · View at Scopus - M. S. Alam, M. M. Rahman, and M. A. Samad, “Dufour and Soret effects on unsteady MHD free convection and mass transfer flow past a vertical porous plate in a porous medium,”
*Nonlinear Analysis. Modelling and Control*, vol. 11, no. 3, pp. 217–226, 2006. View at Google Scholar - M. S. Alam, M. M. Rahman, M. Ferdous, K. Maino, E. Mureithi, and A. Postelnicu, “Diffusion-thermo and thermal-diffusion effects on free convective heat and mass transfer flow in a porous medium with time dependent temperature and concentration,”
*International Journal of Applied Engineering Research*, vol. 2, no. 1, pp. 81–96, 2007. View at Google Scholar - O. A. Beg, T. A. Beg, A. Y. Bakier, and V. R. Prasad, “Chemically-reacting mixed convective heat and mass transfer along inclined and vertical plates with soret and Dufour effects: numerical solution,”
*International Journal of Applied Mathematics and Mechanics*, vol. 5, no. 2, pp. 39–57, 2009. View at Google Scholar - Z. Dursunkaya and W. M. Worek, “Diffusion-thermo and thermal-diffusion effects in transient and steady natural convection from a vertical surface,”
*International Journal of Heat and Mass Transfer*, vol. 35, no. 8, pp. 2060–2065, 1992. View at Google Scholar · View at Scopus - C. R. A. Abreu, M. F. Alfradique, and A. S. Telles, “Boundary layer flows with Dufour and Soret effects: I: forced and natural convection,”
*Chemical Engineering Science*, vol. 61, no. 13, pp. 4282–4289, 2006. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - M. Anghel, H. S. Takhar, and I. Pop, “Dufour and Soret effects on free convection boundary layer over a vertical surface embedded in a porous medium,”
*Studia Universitatis Babes-Bolyai. Matematica*, vol. 11, no. 4, pp. 11–21, 2000. View at Google Scholar - 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. View at Publisher · View at Google Scholar · View at Scopus - B. Gebhart and L. Pera, “The nature of vertical natural convection flows resulting from the combined buoyancy effects of thermal and mass diffusion,”
*International Journal of Heat and Mass Transfer*, vol. 14, no. 12, pp. 2025–2050, 1971. View at Google Scholar · View at Scopus - U. N. Das, R. Deka, and V. M. Soundalgekar, “Effects of mass transfer on flow past an impulsively started infinite vertical plate with constant heat flux and chemical reaction,”
*Forschung im Ingenieurwesen*, vol. 60, no. 10, pp. 284–287, 1994. View at Google Scholar · View at Scopus - S. S. Saxena and G. K. Dubey, “Unsteady MHD heat and mass transfer free convection flow of a polar fluid past a vertical moving porous plate in a porous medium with heat generation and thermal diffusion,”
*Advances in Applied Science Research*, vol. 2, no. 4, pp. 259–278, 2011. View at Google Scholar - K. Raveendra Babu, A. G. Vijaya Kumar, and S. V. K. Varma, “Diffusion-thermo and radiation effects on MHD free convective heat and mass transfer flow past an infinite vertical plate in the presence of a chemical reaction of first order,”
*Advances in Applied Science Research*, vol. 3, no. 4, pp. 2446–2462, 2012. View at Google Scholar - K. Sudhakar, R. Srinibasa Raju, and M. Rangamma, “Chemical reaction effect on an unsteady MHD free convection flow past an infinite vertical accelerated plate with constant heat flux, thermal diffusion and diffusion thermo,”
*International Journal of Modern Engineering Research*, vol. 2, no. 5, pp. 3329–3339, 2012. View at Google Scholar - L. Debanath and D. Bhatta,
*Integral Transforms and Their Applications, Applications of Laplace Transforms*, Chapman and Hall/CRC, Taylor and Francies Group, London, UK, 2007. - K. R. Cramer and S. I. Pai,
*Magnetofluid Dynamics for Engineers and Applied Physics*, McGraw-Hill, New York, NY, USA, 1973. - I. Pop, “Effect of Hall current on hydromagnetic flow near an accelerated plate,”
*International Journal of Physics and Mathematical Sciences*, vol. 5, p. 375, 1971. View at Google Scholar - M. A. Hossain and K. Mohammad, “Effect of hall current on hydromagnetic free convection flow near an accelerated porous plate,”
*Japanese Journal of Applied Physics*, vol. 27, no. 8, pp. 1531–1535, 1988. View at Google Scholar · View at Scopus - P. K. Rath, G. C. Dash, and A. K. Patra, “Effect of Hall current and chemical reaction on MHD flow along an exponentially accelerated porous flat plate with internal heat absorption/generation,”
*Proceedings of the National Academy of Sciences India A*, vol. 80, no. 4, pp. 295–308, 2010. View at Google Scholar · View at Scopus