Abstract

The heat and mass transfer characteristics of natural convection about a vertical surface embedded in a saturated porous medium subject to variable viscosity are numerically analyzed, by taking into account the diffusion-thermo (Dufour) and thermal-diffusion (Soret) effects. The governing equations of continuity, momentum, energy, and concentrations are transformed into nonlinear ordinary differential equations, using similarity transformations, and then solved by using Runge-Kutta-Gill method along with shooting technique. The parameters of the problem are variable viscosity, buoyancy ratio, Lewis number, Prandtl number, Dufour effect, Soret effect, and Schmidt number. The velocity, temperature, and concentration distributions are presented graphically. The Nusselt number and Sherwood number are also derived and discussed numerically.

1. Introduction

Coupled heat and mass transfer by natural convection in a fluid-saturated porous medium has received great attention during the last decades due to the importance of this process which occurs in many engineering, geophysical, and natural systems of practical interest such as geothermal energy utilization, thermal energy storage, and recoverable systems and petroleum reservoirs. A comprehensive account of the available information in this field is provided in recent books by Pop and Ingham [1], Ingham and Pop [2], and Vafai [3].

The previous studies, dealing with the transport phenomena of momentum and heat transfer, have dealt with one-component phases which posses a natural tendency to reach equilibrium conditions. However, there are activities, especially in industrial and chemical engineering processes, where a system contains two or more components whose concentrations vary from point to point. In such a system, there is a natural tendency for mass to be transferred, minimizing the concentration differences within the system and the transport of one constituent, from a region of higher concentration to that of a lower concentration. This is called mass transfer.

When heat and mass transfer occurs simultaneously between the fluxes, the driving potential is of more intricate nature, as energy flux can be generated not only by temperature gradients but by composition gradients as well. The energy flux caused by a composition gradient is called the Dufour or diffusion-thermo effect. Temperature gradients can also create mass fluxes, and this is the Soret or thermal-diffusion effect. Generally, the thermal-diffusion and the diffusion-thermo effects are of smaller-order magnitude than the effects prescribed by Fourier’s or Fick’s laws and are often neglected in heat and mass transfer processes. There are, however, exceptions. The Soret effect, for instance, has been utilized for isotope separation, and in mixtures between gases with very light molecular weight (H2, He). For medium molecular weight (N2, air), the Dufour effect was found to be of a considerable magnitude such that it cannot be neglected (Eckert and Drake) [4].

Kassoy and Zebib [5] studied the effect of variable viscosity on the onset of convection in porous medium. Cheng and Minkowycz [6] studied the effect of free convection about a vertical plate embedded in a porous medium with application to heat transfer from a dike. Bejan and Khair [7] studied the buoyancy-induced heat and mass transfer from a vertical plate embedded in a saturated porous medium. Lai and Kulacki [8] studied the coupled heat and mass transfer by natural convection from vertical surface in a porous medium. The same authors [9] also studied the effect of variable viscosity on convection heat transfer along a vertical surface in a saturated porous medium. Elbashbeshy and Ibrahim [10] investigated the effect of steady free convection flow with variable viscosity and thermal diffusivity along a vertical plate. Kafoussias and Williams [11] studied the thermal-diffusion and diffusion-thermo effects on the mixed free-forced convective and mass transfer steady laminar boundary layer flow, over a vertical plate, with temperature-dependent viscosity. Yih [12] analyzed the coupled heat and mass transfer in mixed convection about a wedge for variable wall temperature and concentration. Jumah and Mujumdar [13] studied the coupled heat and mass transfer for non-Newtonian fluids. Anghel et al. [14] investigated the Dufour and Soret effects on free convection boundary layer over a vertical surface embedded in a porous medium. Kumari [15] analyzed the effect of variable viscosity on free and mixed convection boundary layer flow from a horizontal surface in a saturated porous medium. Postelnicu et al. [16] investigated the effect of variable viscosity on forced convection over a horizontal flat plate in a porous medium with internal heat generation. Seddeek [17], Seddeek and Salem [18] studied the effects of chemical reaction, variable viscosity, and thermal diffusivity on mixed convection heat and mass transfer through porous media. Ali [19] studied the effect of variable viscosity on mixed convection along a vertical plate. Alam et al. [20] analyzed the study of the combined free-forced convection and mass transfer flow past a vertical porous plate in a porous medium with heat generation and thermal diffusion. Pantokratoras [21] analyzed the effect of variable viscosity with constant wall temperature. Partha et al. [22] looked for the effect of double dispersion, thermal-diffusion, and diffusion-thermo effects in free convection heat and mass transfer in a non-Darcy electrically conducting fluid saturating a porous medium. Alam and Rahman [23] studied the Dufour and Soret effects on mixed convective flow past a vertical porous plate with variable suction. Seddeek et al. [24] studied the effects of chemical reaction and variable viscosity on hydromagnetic mixed convection heat and mass transfer through porous media. Another contribution to the theme of Dufour and Soret effects can be found in the paper by Afify [25], where there is a non-Darcy free convection past a vertical surface with temperature viscosity. Narayana and Murthy [26] studied the Soret and Dufour effects in a doubly stratified Darcy porous medium. Postelnicu [27] studied the influence of chemical reactions along with Soret and Dufour’s effect in the absence of magnetic field on free convection. Singh and Chandarki [28] used integral treatment to obtain the expressions for Nusselt number and Sherwood number. El-Arabawy [29] studied the Soret and Dufour effects in a vertical plate with variable surface temperature. Postelnicu [30] analyzed the effect of Soret and Dufour on heat and mass transfer at a stagnation point. Tak et al. [31] investigated the MHD free convection-radiation in the presence of Soret and Dufour. Vempati and Laxmi-Narayana-Gari [32] studied the Soret and Dufour on MHD with thermal radiation. Recently, Cheng [33] studied the Soret and Dufour on heat and mass transfer on a vertical truncated cone with variable wall temperature and concentration. The aim of this paper is to study the effect of Soret and Dufour on heat and mass transfer by natural convection from a vertical plate with variable viscosity.

2. Analysis

Consider a vertical surface embedded in a saturated porous medium. The properties of the fluid and porous medium are isotropic, and the viscosity of the fluid is assumed to be an inverse linear function of temperature. Using Boussinesq and boundary layer approximations, the governing equations of continuity, momentum, energy, and concentration are given byπœ•π‘’+πœ•π‘₯πœ•π‘£πœ…πœ•π‘¦=0,(2.1)𝑒=βˆ’πœ‡ξ‚΅πœ•π‘ξ‚Άπœ…πœ•π‘₯+πœŒπ‘”,(2.2)𝑣=βˆ’πœ‡ξ‚΅πœ•π‘ξ‚Άπ‘’πœ•π‘¦,(2.3)πœ•π‘‡πœ•π‘₯+π‘£πœ•π‘‡πœ•πœ•π‘¦=𝛼2π‘‡πœ•π‘¦2+π·π‘šπ‘˜π‘‡π‘π‘ π‘π‘πœ•2π‘πœ•π‘¦2𝑒,(2.4)πœ•π‘πœ•π‘₯+π‘£πœ•π‘πœ•π‘¦=π·π‘šπœ•2π‘πœ•π‘¦2+π·π‘šπ‘˜π‘‡π‘‡π‘šπœ•2π‘‡πœ•π‘¦2,(2.5)𝜌=πœŒβˆžξ€½ξ€·1βˆ’π›½π‘‡βˆ’π‘‡βˆžξ€Έβˆ’π›½βˆ—ξ€·π‘βˆ’π‘βˆžξ€Έξ€Ύ.(2.6) The viscosity of the fluid is assumed to be an inverse linear function of temperature, and it can be expressed as1πœ‡=1πœ‡βˆžξ€½ξ€·1+π›Ύπ‘‡βˆ’π‘‡βˆžξ€Έξ€Ύ,(2.7) which is reasonable for liquids such as water and oil. Here, 𝛾 is a constant. The boundary conditions are𝑦=0,𝑣=0,𝑇=𝑇𝑀,𝑐=𝑐𝑀,π‘¦βŸΆβˆž,𝑒=0,𝑇=π‘‡βˆž,𝑐=π‘βˆž.(2.8)

3. Method of Solution

Introducing the stream function Ξ¨(π‘₯,𝑦) such that𝑒=πœ•πœ“πœ•π‘¦,𝑣=βˆ’πœ•πœ“πœ•π‘₯,(3.1) where ξ€·πœ“=𝛼𝑓Raπ‘₯ξ€Έ1/2,𝑦(3.2)πœ‚=π‘₯ξ€·Raπ‘₯ξ€Έ1/2,(3.3)Raπ‘₯=ξ‚»π‘˜π‘”π›½Ξ”π‘‡π‘₯ξ‚Όπœˆπ›Ό,istheRayleighnumber(3.4)

Define πœƒ=π‘‡βˆ’π‘‡βˆžπ‘‡π‘€βˆ’π‘‡βˆž,(3.5)πœ™=π‘βˆ’π‘βˆžπ‘π‘€βˆ’π‘βˆžπ›½,(3.6)𝑁=βˆ—ξ€·π‘π‘€βˆ’π‘βˆžξ€Έπ›½ξ€·π‘‡π‘€βˆ’π‘‡βˆžξ€Έ.(3.7) Substituting these transformations (3.1) to (3.7) to Equations (2.2) to (2.5) along with (2.6) and (2.7), the resulting nonlinear ordinary differential equations areπ‘“β€²ξ…ž=π‘“β€²πœƒβ€²πœƒβˆ’πœƒπ‘Ÿβˆ’ξ‚΅πœƒβˆ’πœƒπ‘Ÿπœƒπ‘Ÿξ‚Άπœƒ(πœƒβ€²+πœ™β€²π‘),(3.8)ξ…žξ…ž+PrDuπœ™β€²ξ…ž+12πœ™π‘“πœƒβ€²=0,(3.9)ξ…žξ…ž+ScSrπœƒβ€²ξ…ž+Le2π‘“πœ™ξ…ž=0,(3.10) together with the boundary conditionsπœ‚=0,𝑓=0,πœƒ=1,πœ™=1,πœ‚βŸΆβˆž,𝑓′=0,πœƒ=0,πœ™=0,(3.11) where Pr=𝜈/𝛼 is the Prandtl number, Du=π·π‘šπ‘˜π‘‡(π‘π‘€βˆ’π‘βˆž)/π‘π‘ π‘π‘πœˆ(π‘‡π‘€βˆ’π‘‡βˆž) is the Dufour number, Sc=𝜈/π·π‘š is the Schmidt number, Sr=π·π‘šπ‘˜π‘‡(π‘‡π‘€βˆ’π‘‡βˆž)/π‘‡π‘šπœˆ(π‘π‘€βˆ’π‘βˆž) is the Soret number, Le=𝛼/π·π‘š is the Lewis number, and πœƒπ‘Ÿ=βˆ’1/𝛾(π‘‡π‘€βˆ’π‘‡βˆž) is the parameter characterizing the influence of viscosity. For a given temperature differential, large values of πœƒπ‘Ÿ imply that either 𝛾 or (π‘‡π‘€βˆ’π‘‡βˆž) is small. In this case, the effect of variable viscosity can be neglected.

The effect of variable viscosity is important if πœƒπ‘Ÿ is small. Since the viscosity of liquids decreases with increasing temperature while it increases for gases, πœƒπ‘Ÿ is negative for liquids and positive for gases. The concept of this parameter πœƒπ‘Ÿ was first introduced by Ling and Dybbs [34] in their study of forced convection flow in porous media. The parameter 𝑁 measures the relative importance of mass and thermal diffusion in the buoyancy-driven flow. It is clear that 𝑁 is zero for thermal-driven flow, infinite for mass-driven flow, positive for aiding flow, and negative for opposing flow.

4. Numerical Analysis and Discussion

Equations (3.8), (3.9), and (3.10) are integrated numerically by using Runge-Kutta-Gill method along with shooting technique. The parameters involved in this problem are πœƒπ‘Ÿ: the variable viscosity, Le: Lewis number, 𝑁: the buoyancy parameter, Pr: Prandtl number, Du: Dufour number, Sr: Soret number, and Sc: Schmidt number. To observe the effect of variable viscosity on heat and mass transfer, we have plotted the velocity function 𝑓′, temperature function πœƒ, and the concentration function Ξ¦ against πœ‚ for various values of πœƒπ‘Ÿ, Le, and 𝑁. The value of Prandtl number Pr is taken equal to 0.71 which corresponds to air. The values of the Dufour number and Soret number are taken in such a way that their product is constant according to their definition provided that the mean temperature is kept constant as well. The parameter πœƒπ‘Ÿ is used to represent the effect of variable viscosity. The case πœƒπ‘Ÿ<0 corresponds to the case of liquids, and πœƒπ‘Ÿ>0 corresponds to the case of gases.

The influence of Dufour number Du and Soret number Sr on velocity, temperature, and concentration profiles is shown in Figures 1, 2, and 3, respectively, for πœƒπ‘Ÿ=5, 𝑁=1, Le=0.1, Pr=.71, and Sc=0.1. Figure 1 shows that as Dufour number increases, the velocity decreases slightly for decrease in Soret number. From Figure 2, it is evident that as Dufour number increases, the temperature increases for decrease in Soret number. Figure 3 shows that as Dufour number increases, the concentration decreases for decrease in Soret number.

The effect of variable viscosity πœƒπ‘Ÿ on velocity, temperature, and concentration is shown in Figures 4, 5, and 6, respectively, for 𝑁=1, Le=0.1, Pr=.71, Sc=0.1, Du=0.03, and Sr=2. From Figure 4, it is realized that the velocity increases near the plate and decreases away from the plate as πœƒπ‘Ÿβ†’0 in the case of liquids (πœƒπ‘Ÿ<0) and decreases near the plate and increases away from the plate as πœƒπ‘Ÿβ†’0 in the case of gases (πœƒπ‘Ÿ>0). From Figures 5 and 6, it is evident that the temperature and concentration increase as πœƒπ‘Ÿβ†’0 for πœƒπ‘Ÿ>0 (i.e., for gases) and decrease as πœƒπ‘Ÿβ†’0 for πœƒπ‘Ÿ<0 (i.e., for liquids).

The effect of buoyancy ratio 𝑁 on velocity, temperature, and concentration is shown in Figures 7, 8, and 9, respectively, for πœƒπ‘Ÿ=5, Le=0.1, and Pr=.71 and Sc=0.1, Du=0.03, and Sr=2. From Figure 7, it is observed that the velocity increases near the plate and decreases away from the plate as 𝑁 increases. From Figures 8 and 9, it is realized that the temperature and concentration decrease as 𝑁 increases.

The effect of Lewis number Le on temperature and concentration is shown in Figures 10 and 11, respectively, for 𝑁=1, Pr=.71, Sc=0.1, Du=0.03, Sr=2, and for different values of πœƒπ‘Ÿ and Le. From Figure 10, it is observed that as Le increases, the heat transfer decreases for both πœƒπ‘Ÿ>0 and πœƒπ‘Ÿ<0. From Figure 11, it is realized that as Le increases, the mass transfer increases for both πœƒπ‘Ÿ>0 and πœƒπ‘Ÿ<0.

The effect of Prandtl number Pr on velocity, temperature, and concentration is shown in Figures 12, 13, and 14, respectively, for 𝑁=1, Sc=0.1, Du=0.03, Sr=2, Le=0.1, and πœƒπ‘Ÿ=5. From Figures 12 and 13, it is evident that as the Prandtl number increases, the velocity and temperature increase. Figure 14 shows that as the prandtl number increases, the concentrations decrease.

The parameters of engineering interests for the present problem are the local Nusselt number Nuπ‘₯ and Sherwood number Shπ‘₯, which are given byNuπ‘₯ξ€·Raπ‘₯ξ€Έ1/2=βˆ’πœƒξ…ž,(0)Shπ‘₯ξ€·Raπ‘₯ξ€Έ1/2=βˆ’πœ™ξ…ž.(0)(4.1)

The values of Nusselt number and Sherwood number for different values of variable viscosity πœƒπ‘Ÿ and the buoyancy ratio 𝑁 are presented in Table 1 for Le=0.1, Sc=0.1, Du=0.03, Sr=2, and Pr=.71. It is evident that as πœƒπ‘Ÿβ†’0 for πœƒπ‘Ÿ>0, the Nusselt number decreases and the Sherwood number increases. It is also realized that as πœƒπ‘Ÿβ†’0 for πœƒπ‘Ÿ<0, the Nusselt number increases and the Sherwood number decreases for other parameters fixed.

5. Conclusion

The Dufour and Soret effect on free convective heat and mass transfer flow past a semi-infinite vertical plate under the influence of variable viscosity has been studied. Using usual similarity transformations, the governing equations have been transformed into nonlinear ordinary differential equations. The similarity solutions are obtained numerically by applying Runge-Kutta-Gill method along with shooting technique. The effects of the variable viscosity parameter πœƒπ‘Ÿ, the Lewis number Le, the buoyancy ratio 𝑁, the Dufour number Du, the Soret number Sr, the prandtl number Pr, and the Schmidt number Sc on the velocity, temperature, and concentration profiles are examined. From the present study, we see that the thermal and species concentration boundary layer thickness increases for gases and decreases for liquids.

Nomenclature

𝑐:Concentration at any point in the flow field
𝑐𝑝:Specific heat at constant pressure
𝑐𝑠:Concentration susceptibility
𝑐𝑀:Concentration at the wall
π‘βˆž:Concentration at the free stream
π·π‘š:Mass diffusivity
Du:Dufour number
𝑓:Dimensionless velocity function
𝑔:Acceleration due to gravity
π‘˜:Permeability
π‘˜π‘‡:Thermal diffusion ratio
Le:Lewis number [Le=𝛼/π·π‘š]
𝑁:Buoyancy ratio [𝑁=π›½βˆ—(π‘π‘€βˆ’π‘βˆž)/𝛽(π‘‡π‘€βˆ’π‘‡βˆž)]
Nuπ‘₯:Nusselt number [Nuπ‘₯=βˆ’π‘₯(πœ•π‘‡/πœ•π‘¦)𝑦=0/(π‘‡π‘€βˆ’π‘‡βˆž)]
𝑝:Pressure
Pr:Prandtl number
Raπ‘₯:Rayleigh number Raπ‘₯=[π‘˜π‘”π›½Ξ”π‘‡π‘₯/πœˆπ›Ό]
Sc:Schmidt number
Shπ‘₯:Sherwood number [Shπ‘₯=π‘šπ‘₯/𝐷(π‘π‘€βˆ’π‘βˆž)]
Sr:Soret number
𝑇:Temperature of the fluid
𝑇𝑀:Temperature of the plate
π‘‡βˆž:Temperature of the fluid far from the plate
𝑒,𝑣:Velocity components in π‘₯ and 𝑦 direction
π‘₯,𝑦:Coordinate system.
Greek Letters
𝛼:Thermal diffusivity
𝛽:Coefficient of thermal expansion
π›½βˆ—:Concentration expansion coefficient
𝛾:Constant defined in(2.7)
πœ‚:Dimensionless similarity variable
πœƒ:Dimensionless temperature
πœƒπ‘Ÿ:βˆ’1/𝛾(π‘‡π‘€βˆ’π‘‡βˆž)
πœ‡:Viscosity [pas]
𝜈:Kinematic viscosity
𝜌:Density
Ξ¦:Dimensionless concentration
πœ“:Dimensionless stream function.

Acknowledgments

The authors tender their heartfelt thanks to Dr. T. Govindarajulu the former Professor and Head of the Department of Mathematics, Anna University, Chennai, for his generous help. The authors wish to thank the Director and the Principal of IRTT for their kind support. The authors also wish to thank Dr. K. Thyagarajah, The Principal, K. S. Rangasamy College of Technology, for his help and guidance to do this work.