Journal of Applied Mathematics

Volume 2013 (2013), Article ID 584534, 11 pages

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

Research Article

Thermal Diffusion and Diffusion Thermo Effects on the Viscous Fluid Flow with Heat and Mass Transfer through Porous Medium over a Shrinking Sheet

Nabil Eldabe¹ and Mahmoud Abu Zeid²

¹Department of Mathematics, Faculty of Education, Ain Shams University, Roxy, Heliopolis, Cairo, Egypt

²Basic Science Department, Higher Technological Institute, 10th of Ramadan City, Egypt

Received 9 February 2013; Revised 16 April 2013; Accepted 18 April 2013

Academic Editor: Magdy A. Ezzat

Copyright © 2013 Nabil Eldabe and Mahmoud Abu Zeid. 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 motion of the viscous, incompressible fluid through a porous medium with heat and mass transfer over a shrinking sheet is investigated. The cross-diffusion effect between temperature and concentration is considered. This phenomenon is modulated mathematically by a set of partial differential equations which govern the continuity, momentum, heat, and mass. These equations are transformed to a set of ordinary differential equations by using similarity solutions. The analytical solutions of these equations are obtained. The velocity, temperature, and concentration of the fluid as well as the heat and mass transfer with shear stress at the sheet are obtained as a function of the physical parameters of the problem. The effects of Prandtl number, mass transfer parameter, the wall shrinking parameter, the permeability parameter, and Dufour and Soret numbers on temperature and concentration are studied. Also, the effects of mass transfer parameter, permeability parameter, and shrinking strength on the velocity and shear stress are discussed. These effects are illustrated graphically through a set of figures.

1. Introduction

The study of two-dimensional boundary layer flow of an incompressible fluid over a stretching or shrinking surface has assumed significance in recent years because of its extensive applications in engineering disciplines. Some applications are glass fiber production, hot rolling and wire drawing, the aerodynamic extrusion of plastic sheets, paper production, glass blowing, and drawing plastic films. The flow over a continuously stretching surface with a constant speed was initiated by Sakiadis [1, 2] and had been extended to the different physical situations [3–5]. The effect of large suction on the MHD forced and free convection flow past a vertical porous plate was studied by Ferdows et al. [6]. Lin et al. [7] studied the flow characteristics of an electrically conducting second-order fluid over a stretching sheet with suction velocity in the presence of a transverse magnetic field. P. S. Gupta and A. S. Gupta [8] studied the heat and mass transfer in a stretching surface with suction or injection. Chen and Char [9] studied the effects of variable surface temperature and variable surface heat flux on the heat transfer characteristics of a linearly stretching sheet. The MHD effect on a vertical stretching surface with suction and blowing was studied by Gorla et al. [10]. Yao et al. [11] investigated heat transfer problem with a convective boundary condition for a viscous and incompressible fluid over a permeable (with mass flux) stretching/shrinking sheet in a quiescent fluid. Fang et al. [12] investigate the behavior of the steady boundary layer flow and heat transfer of a viscous and incompressible fluid towards a permeable stretching/shrinking sheet in a quiescent fluid.

Recently, the flow over a shrinking sheet has attracted the attention of many researchers. Fang [13], Fang et al. [14], Fang and Zhang [15], and Fang et al. [16] studied problems of boundary layer flow over shrinking sheets with mass transfer. The flow induced by a shrinking sheet shows physical phenomena distinct from the forward stretching flow. A steady flow over a shrinking sheet is not possible because the generated vorticity is not confined within the boundary layer. So, to overcome this difficulty the flow needs a certain amount of external opposite force at the sheet. The most suitable external force is the suction at the sheet. Miklavčič and Wang [17] investigated the shrinking flow where velocity of boundary layer moves toward a fixed point. They found an exact solution of the Navier-Stokes equations and reported that an adequate suction is necessary to maintain the flow over a shrinking sheet. This phenomenon can be found, for example, on a rising and shrinking balloon [18].

The heat and mass transfer that simultaneously affecting each other will cause the cross-diffusion effect. The heat transfer caused by concentration gradient is called the diffusion-thermo or Dufour effect. On the other hand mass transfer caused by temperature gradient is called Soret or thermal-diffusion effect. Alam and Rahman [19] investigated the Dufour and Soret effects on mixed convection flow past a vertical porous flat plate with variable suction. El-Arabawy [20] investigated the heat and mass transfer by natural convection from vertical surface embedded in a fluid-saturated porous media considering Soret and Dufour effects with variable surface temperature and constant concentration.

Motivated by the above-mentioned investigations, we investigate in this present paper the behavior of flow of a viscous and incompressible fluid through a porous medium over a continuously shrinking sheet with heat and mass transfer in a stationary fluid. The solutions for the momentum, heat, and mass transfer equations are solved analytically, and the solutions are presented in terms of algebraically decaying function. The effects of various physical parameters like the Prandtl number, the mass transfer parameter, the wall shrinking parameter, the permeability parameter, Dufour number, and Soret number on the temperature and concentration profiles are obtained. Also, the effect of mass suction, permeability of the porous medium, and shrinking strength on the velocity and shear stress are studied.

2. Basic Equations and Exact Solutions

Consider the steady two-dimensional flow of a viscous and incompressible fluid through a porous medium over a continuously shrinking sheet with heat and mass transfer in a stationary fluid. It is assumed that the velocity of the shrinking sheet is , where is a real number. It is also assumed that constant mass transfer velocity is . In the present study, the effective viscosity is considered to be identical to the dynamic viscosity . The ambient fluid temperature and concentration are and , and these at the sheet surface are and . The -axis is along the shrinking surface and -axis is perpendicular to it (see Figure 1).

Figure 1: The physical configuration of the problem.

With these assumptions the governing equations can be written as

where and are the velocity components along the - and -axes, respectively, is the pressure, is the density, is the fluid kinematic viscosity, is the permeability of porous medium, is the fluid thermal diffusivity, is the coefficient of mass diffusivity, is the specific heat at constant pressure, is the mean fluid temperature, is the thermal diffusion ratio, is the concentration susceptibility, is the fluid temperature, and is the concentration of the fluid. The appropriate boundary conditions are

We assume that (1) to (5) subject to the boundary conditions (6) admit the similarity solutions

where prime denotes differentiation with respect to , and is a positive constant. Also we denote as a reference velocity for this problem and . Using (3) and the boundary conditions (6), we obtain the following expression for the pressure :

where is the stagnation pressure. Substituting (7) and (8) into (2), (4), and (5), we obtain the following nonlinear ordinary differential equations:

subject to the boundary conditions

where is the permeability parameter, is the Prandtl number of the fluid with , is the Dufour number, is the Schmidt number, is the Soret number, and is the mass transfer parameter showing the strength of the mass transfer at the sheet. Also is for wall stretching and is for shrinking, respectively.

A physical quantity of interest is the skin friction coefficient , which is defined as

where is the shear stress and is given by

Using the similarity variables (7), we obtain

where is the local Reynolds number.

The momentum equation (9) together with the boundary conditions (12) yields an algebraically decaying solution as follows:

The velocity fields become

It is seen that the algebraically decaying function exists only for a shrinking sheet with . The mass suction at the sheet is

There exists another exact solution for (9) with the boundary conditions (12) as

with

Eliminating from (10) and (11), we obtain

Differentiating (20) with respect to :

Eliminating between (10) and (21), we obtain the Euler-Cauchy equation:

Or,

The solution can be found as follows:

where

In order to match the BC at , it is required that and , and

Integrating (20) with respect to :

where is a constant.

Substituting from (26) into (27), we obtain

where

In order to match the boundary conditions at , it is required that and , and , so (28) becomes

Applying the boundary conditions at the wall yields

Solving (31), we obtain

Then the temperature and concentration solutions become

The heat and mass transfer rates at the wall are related to the temperature and concentration gradients at the wall are, respectively,

where is the thermal conductivity of the fluid.

The heat and mass fluxes at the wall read

3. Results and Discussion

The nondimensional temperature profiles for a shrinking sheet are shown in Figures 2(a)–2(f). The effect of the Prandtl number on the temperature profile is shown in Figure 2(a). It is observed that the temperature increases with decreasing values of the Prandtl number. With the increase of Prandtl number, the temperature drops faster and the boundary layer becomes thinner. Figure 2(b) illustrates the influence of the mass transfer parameter (mass suction) for a shrinking sheet problem on the temperature profile. It is found that the temperature increases with the decrease of suction parameter. With the increase in mass suction, the temperature drops faster and the boundary layer becomes thinner. Figure 2(c) represents the temperature profile for different values of shrinking sheet strength. It is noticed from the figure that the temperature decreases with decreasing the magnitude of shrinking sheet strength. The boundary layer becomes thicker for a high magnitude of shrinking strength. Figure 2(d) represents the temperature profile for different values of Permeability parameter . From the graph we observe that the increase of Permeability parameter leads to the increase of the temperature profile. The temperature distribution for different values of is plotted in Figure 2(e). It is observed from this figure that the temperature distribution is more for higher values of which increases with an increase in the value of . Higher values of gives rise to the formation of peak near the sheet. Figure 2(f) shows the variation of concentration profiles for various values of Soret number . We notice that for , the temperature distribution increases with increasing .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 2: (a) Temperature profiles for a shrinking sheet problem under different values of Prandtl number. (b) Effects of the mass transfer parameter (mass suction) on the temperature profiles of a shrinking sheet. (c) Temperature profiles for different values of shrinking sheet strength with mass suction. (d) Effects of the permeability parameter on the temperature profiles of a shrinking sheet with mass suction. (e) Effects of Dufour number on the temperature profiles of a shrinking sheet with mass suction. (f) Effects of Soret number on the temperature profiles of a shrinking sheet with mass suction.

Figure 3(a) displays the effect of Prandtl number on the concentration profile under certain shrinking strength, and Permeability parameter. It can be seen that for , the concentration increases when the Prandtl number decreases. The concentration boundary layer thickness reduces with the increase of the Prandtl number. The effect of the mass transfer parameter (mass suction) on the dimensionless concentration profile for sheet shrinking problem is shown in Figure 3(b). We observe that the concentration increases with the decrease of mass suction. The concentration boundary layer thickness decreases with the increase of mass transfer parameter (mass suction). Figure 3(c) represents the concentration profile for different values of shrinking sheet strength. From the graph we observe that the concentration decreases with decreasing the magnitude of shrinking sheet strength and the concentration boundary layer thickness becomes thinner. Figure 3(d) displays the concentration profile for different values of Permeability parameter. We notice from the figure that the effect of increasing values of Permeability parameter is to increase the concentration profile. Figure 3(e) shows the concentration profiles for different values of Dufour number . It is observed that for , the concentration distribution increases with increasing the Dufour number . The concentration distribution for different values of Soret number is plotted in Figure 3(f). It is noticed from this figure that the temperature distribution is more for higher values of which increases with an increase in the value of and peaks appear near the sheet for higher values of .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3: (a) Concentration profiles for a shrinking sheet under different values of Prandtl number. (b) Effects of the mass transfer parameter (mass suction) on the concentration profiles of a shrinking sheet. (c) Concentration profiles for different values of shrinking sheet strength with mass suction. (d) Effects of the permeability parameter on the concentration profiles of a shrinking sheet with mass suction. (e) Effects of Dufour number on the concentration profiles of a shrinking sheet with mass suction. (f) Effects of Soret number on the concentration profiles of a shrinking sheet with mass suction.

The effect of mass suction parameter on the velocity profile for shrinking sheet is displayed in Figure 4(a) which shows that the velocity decreases with the increase of the mass suction parameter. The boundary layer is sucked closer to the wall for smaller mass suction parameters. Figure 4(b) displays the velocity profile for various values of shrinking sheet parameter. The velocity decreases and drops faster with high magnitude of sheet shrinking parameter . The velocity at the sheet decreases with increasing the magnitude of shrinking sheet parameter . There exists cross-over points among the velocity profiles for different values of . Figure 4(c) shows velocity profile for different values of the permeability parameter for shrinking sheet. We remark that as permeability parameter increases the velocity decreases.

(a)

(b)

(c)

Figure 4: (a) Velocity profiles for a shrinking sheet under different values of mass suction parameter. (b) Velocity profiles under different values of shrinking sheet strength with mass suction. (c) Velocity profiles for a shrinking sheet under different values of permeability parameter.

Shear stress results for shrinking sheet are seen in Figure 5(a) under the effect of mass suction. The shear stress drops faster as the mass suction becomes smaller. The increasing effect of shear stress is found in the interval with decreasing . The wall shear stress decreases with increasing mass suction. There exists cross-over points among the shear profiles for different values of . In Figure 5(b) the shear stress increases with increasing magnitude of shrinking sheet parameter. The shear drops faster when magnitude of shrinking sheet parameter increases and the wall shear stress decreases with decreasing magnitude of shrinking parameter. Figure 5(c) illustrates the shear stress profile for a shrinking sheet with mass suction. The wall shear stress increases permeability parameter. The increasing effect of shear stress is found in the interval with increasing .

(a)

(b)

(c)

Figure 5: (a) Shear stress profiles for a shrinking sheet under different values of mass suction parameter. (b) Shear stress profiles under different values of shrinking sheet strength with mass suction. (c) Shear stress profiles for a shrinking sheet under different values of permeability parameter.

Figures 6(a) and 6(b) and Figures 7(a) and 7(b) show the rate of heat and mass transfer at the sheet surface versus mass transfer parameter and sheet shrinking parameter , respectively, with different values of Prandtl number . We observe that increase with increasing the values of whereas reverse effect is observed on .

(a)

(b)

Figure 6: (a) Wall heat flux, with for various values of . (b) Wall heat flux, with for various values of .

(a)

(b)

Figure 7: (a) Wall mass flux, with for various values of . (b) Wall mass flux, with for various values of .

Figures 8(a) and 8(b) depict the effect of sheet shrinking parameter on wall shear stress versus mass suction and permeability parameter , respectively. We notice that increase with increasing the magnitude of shrinking sheet parameter .

(a)

(b)

Figure 8: (a) Wall shear stress, with for various values of . (b) Wall shear stress, with for various values of .

4. Conclusion

In this work, the solutions for the momentum, heat, and mass transfer equations are obtained in a closed analytical form and the solutions presented in terms of algebraically decaying function. The effects of the Prandtl number, the mass transfer parameter, the wall shrinking parameter, and the Permeability parameter on the temperature and concentration profiles were analyzed. Also, the effects of mass suction, Permeability of the porous medium, and shrinking strength on the velocity and shear stress are studied. We found some of the important results from the graphical representation which are listed below.

(1)The effect of Prandtl number is to reduce the thickness of both the temperature and the concentration boundary layer.(2)The effect of mass transfer parameter (mass suction) is seen to decrease both of the temperature and the concentration distributions in the flow region.(3)Increasing magnitude of shrinking sheet strength is to increase both of the temperature and the concentration distributions in the flow region (4)The effect of Permeability parameter is seen to increase both of the temperature and the concentration distributions in the flow region.(5)Dufour effect is to increase both temperature and concentration distributions with formation of peak for higher values of Dufour parameter in the thermal boundary layer.(6)Soret effect is to increase both temperature and concentration distributions with formation of peak for higher values of Soret parameter in the concentration boundary layer.(7)The velocity decreases with the increase of the mass suction parameter , the magnitude of shrinking sheet strength , and the Permeability parameter .(8)The shear stress increases with the increase of the mass suction parameter , magnitude of shrinking sheet strength parameter , and Permeability parameter .(9) versus , increases with increasing the values of whereas reverse effect is observed on . Also, versus increases with increasing the values of whereas reverse effect is observed on . (10) increases with increasing the magnitude of shrinking sheet parameter .

References

B. C. Sakiadis, “Boundary-layer behavior on continuous solid surface: I. Boundary-layer equations for two-dimensional and axisymmetric flow,” AIChE Journal, vol. 7, pp. 26–28, 1961.

B. C. Sakiadis, “Boundary-layer behavior on continuous solid surface: II. Boundary-layer equations for two-dimensional and axisymmetric flow,” AIChE Journal, vol. 7, pp. 221–225, 1961.

L. J. Crane, “Flow past a stretching plate,” Journal of Applied Mathematics and Physics, vol. 21, no. 4, pp. 645–647, 1970.

P. S. Gupta and A. S. Gupta, “Heat and mass transfer on a stretching sheet with suction or blowing,” The Canadian Journal of Chemical Engineering, vol. 55, pp. 744–746, 1977.

C. K. Chen and M. I. Char, “Heat transfer of a continuous, stretching surface with suction or blowing,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 568–580, 1988.

M. Ferdows, M. Ota, A. Sattar, and M. Alam, “Similarity solutions for MHD through vertical porous plate with suction,” Journal of Computational and Applied Mathematics, vol. 6, pp. 15–25, 2005.

J. R. Lin, L. J. Liang, and R. D. Chien, “Magneto-hydrodynamic flow of a second order fluid over a stretching sheet with suction,” Journal of the Chinese Institute of Engineers, vol. 30, no. 1, pp. 183–188, 2007.

P. S. Gupta and A. S. Gupta, “Heat and mass transfer on a stretching sheet with suction or blowing,” The Canadian Journal of Chemical Engineering, vol. 55, pp. 744–746, 1977.

C. K. Chen and M. I. Char, “Heat transfer of a continuous, stretching surface with suction or blowing,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 568–580, 1988.

R. S. R. Gorla, D. E. Abboud, and A. Sarmah, “Magnetohydrodynamic flow over a vertical stretching surface with suction and blowing,” Heat and Mass Transfer, vol. 34, no. 2-3, pp. 121–125, 1998.

S. Yao, T. Fang, and Y. Zhong, “Heat transfer of a generalized stretching/shrinking wall problem with convective boundary conditions,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 2, pp. 752–760, 2011.

T. Fang, S. Yao, and I. Pop, “Flow and heat transfer over a generalized stretching/shrinking wall problem-Exact solutions of the Navier-Stokes equations,” International Journal of Non-Linear Mechanics, vol. 46, pp. 1116–1127, 2011.

T. Fang, “Boundary layer flow over a shrinking sheet with power-law velocity,” International Journal of Heat and Mass Transfer, vol. 51, no. 25-26, pp. 5838–5843, 2008.

T. Fang, W. Liang, and C. F. Lee, “A new solution branch for the Blasius equation— shrinking sheet problem,” Computers & Mathematics with Applications, vol. 56, no. 12, pp. 3088–3095, 2008.

T. Fang and J. Zhang, “Closed-form exact solutions of MHD viscous flow over a shrinking sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 7, pp. 2853–2857, 2009.

T. Fang, J. Zhang, and S. Yao, “Slip MHD viscous flow over a stretching sheet—an exact solution,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 11, pp. 3731–3737, 2009.

M. Miklavčič and C. Y. Wang, “Viscous flow due to a shrinking sheet,” Quarterly of Applied Mathematics, vol. 64, no. 2, pp. 283–290, 2006.

C. Y. Wang, “Stagnation flow towards a shrinking sheet,” International Journal of Non-Linear Mechanics, vol. 43, no. 5, pp. 377–382, 2008.

M. S. Alam and M. M. Rahman, “Dufour and Soret effects on mixed convection flow past a vertical porous flat plate with variable suction,” Nonlinear Analysis. Modelling and Control, vol. 11, no. 1, pp. 3–12, 2006.

H. A. M. El-Arabawy, “Soret and dufour effects on natural convection flow past a vertical surface in a porous medium with variable surface temperature,” Journal of Mathematics and Statistics, vol. 5, no. 3, pp. 190–198, 2009.