ARTICLES
Dufour and Soret effects on unsteady MHD convective heat and mass transfer flow due to a rotating disk
Kh. A. Maleque
Associate
Professor, Department of Mathematics, American International
University-Bangladesh, House-23, Road-17, Kemal Ataturk Avenue,
Banani, Dhaka-1213, Bangladesh.
maleque@aiub.edu, khmaleque@yahoo.com
Abstract - An unsteady three dimensional MHD convective heat and mass transfer flow in an incompressible fluid due to a rotating disk is studied by taking into account the Dufour and Soret effects. The system of axial symmetric non-linear partial differential equations governing the unsteady flow, heat and mass transfer is written in cylindrical polar coordinates and reduced to nonlinear ordinary differential equations by similarity transformations. The resulting system of ordinary nonlinear differential equations is then solved numerically by a shooting method using Range-Kutta six order integration scheme. The flow, temperature and concentration fields are affected by the magnetic interaction parameter M, Rotational parameter R, Soret Number S0 and Dufour number Df respectively. The results of the numerical solution are presented graphically in the form of velocity, temperature and concentration profiles. The results for the wall flow, temperature and concentration gradients obtained are presented in tabular form for various values of the parameters M, R, S0 and Df.
Keywords - Dufour and Soret Effects; MHD Flow; Rotating Disk.
I. INTRODUCTION
The solution of the Navier-Stokes equations is furnished by the flow around a flat disk which rotates about an axis perpendicular to its plane with uniform angular velocity O, in a fluid otherwise at rest. The layer near the disk is carried by it through friction and is thrown outwards owing to the action of centrifugal forces. This is compensated by particles which flow in an axial direction towards the disk to be in turn carried and ejected centrifugally. Thus the case is seen to be one of fully three dimensional flows. The problem of flow over a rotating disk has a long history, with von-Karman (1921) originally describing similarity transformations that enable the Navier-Stokes equations for isothermal, impermeable rotating disk to be reduced to a system of coupled ordinary differential equations. Later, Cochran (1934) obtained more accurate results by patching two series expansions. It is found that the disk acts like a centrifugal fan, the fluid near the disk being thrown radially outwards. This in turn impulses an axial flow towards the disk to maintain continuity. Benton (1965) improved Cochran's solution and extended the hydrodynamics problem to flow starting impulsively from rest. Subsequently comprehensive studies have been carried out on rotating disk flow by many researchers, some of them are Wagner (1948), Sparrow and Gregg (1959), Kuiken (1971), Hassan and Attia (1997), Attia (1998) and Kelson and Desseaux (2000). Maleque and Sattar (2003) considered an implicit finite difference solution for an unsteady convective incompressible MHD flow due to an infinite rotating disk, taking into account heat absorption effect. Maleque and Sattar (2005a and 2005b) studied the effects of variable properties and Hall current on steady MHD laminar flow of a compressible electrically conducting fluid on a porous rotating disk. Sajid and Hayat (2008) studied the influence of thermal radiation on the boundary layer flow due to an exponentially stretching sheet. Recently, Hayat et al. (2008c) investigated the influence of heat transfer in an MHD second grade fluid film over an unsteady stretching sheet. More recently, Hayat et al. (2008d) studied hall and heat transfer effects on the steady flow of generalized Burgers' fluid induced by a sudden pull of eccentric rotating disks.
The effects of thermal diffusion on MHD convection of mass transfer flows have been considered by many investigators due to its important role particularly in isotope separation and in mixtures between gases with very light molecular weight (H2, He) and medium molecular weight (N2, air) (Eckert and Drake, 1972). Considering these aspects, many papers have been published on the effects on thermal diffusion on MHD convective and mass transfer flow. Some of them are Jha and Singh (1990), Kafousias (1992) and Alam and Sattar (1998). Sattar and Alam (2001) obtained an analytical solution on the free convection and mass transfer flow with thermal diffusion. Maleque and Sattar (2002) obtained the numerical solution of MHD free-convective and mass transfer flow over an infinite vertical porous plate with thermal diffusion effects. Recently, Hayat et al. (2008a) considered mixed convection flow of a micropolar fluid over a non-linearly stretching sheet. More recently Hayat et al. (2008b) studied heat and mass transfer analysis on the flow of a second grade fluid in the presence of chemical reaction. Considering the importance of thermal and mass diffusions of MHD convective heat and mass transfer flow due to rotating disk, in the present study a numerical solution is obtained for an unsteady three dimensional MHD convective heat and mass transfer flow in an incompressible fluid due to a rotating disk with thermal and mass diffusions. The governing partial differential equations of the MHD convective mass transfer flow are reduced to nonlinear ordinary differential equations by introducing suitable similarity transformations. The nonlinear similarity equations are then solved numerically by Nachtshiem-Swigert (1965) iteration technique. The results of the numerical solution are then presented graphically in the form of velocity, temperature and concentration profiles. The corresponding skin friction coefficients, the Nusselt number and the Sherwood number are also calculated and displayed in tables showing the effects of various parameters on them.
II. GOVERNING EQUATIONS
Unsteady, axially-symmetric, incompressible flow of a homogeneous electrically conducting fluid with heat and mass transfer due to an infinite rotating disk has been considered, assuming that the fluid is infinite in extent in the positive z-direction.
The disk
rotates with constant angular velocity O and is placed at
z=0, where z is the vertical axis in the cylindrical
coordinates system with r and
as the radial and tangential axes respectively. The components of
the flow velocity are (u,v,w) in the
directions of increasing (r,
,z) respectively, the pressure is
P and the density of the fluid is ?. T and
C are the fluid temperature and concentration respectively
and the surface of the rotating disk is maintained at a uniform
temperature cw and uniform concentration
Tw. Far away from the surface, the free
stream is kept at a constant temperature
T∞, at a constant concentration
C∞ and at a constant pressure
P∞.
The fluid is assumed to be Newtonian, viscous and electrically conducting. The external uniform magnetic field is applied perpendicular to the surface of the disk and has a constant magnetic flux density B0 which is assumed unchanged by taking small magnetic Reynolds number (Rem<<1). In addition, no electric field is assumed to exist and the Hall effect is negligible.
The physical model and geometrical coordinates are shown in Fig. 1. Due to axial symmetry, incompressible MHD laminar flow of a homogeneous fluid satisfies the following continuity, Navier-Stokes, energy and mass equations:
![]() |
(1) |
![]() |
(2) |
![]() |
(3) |
![]() |
(4) |
![]() |
(5) |
![]() |
(6) |
here, ? is the thermal conductivity of heat, S is the electrical conductivity, cp is the specific heat with constant pressure, Dm is the molecular diffusivity, Tm is the mean fluid temperature, kT is the thermal diffusion ratio and cs is the concentration susceptibility.
Fig.1 The flow
configuration and the coordinate system.
Appropriate boundary conditions are given by
![]() |
(7) |
III. MATHEMATICAL FORMULATIONS
To solve the governing Eqs. (1-7), they are converted to some suitable form by introducing the following dimensionless quantities:
![]() |
(8) |
where ? is a scale factor and is a function of time as
![]() |
(9) |
![]() |
(10) |
![]() |
(11) |
![]() |
(12) |
is the rotational parameter.
is the magnetic interaction parameter which
represents the ratio of the magnetic force to the fluid inertia
force,
is the Prandtl number,
is the Schmidt number and
is the Soret number and
is the Dufour number, where, ?T =
Tw - T∞ and
?C = Cw -
C∞.
The Eqs.
(10)- (12) depend on a single similarity variable ?
except a time dependent term where
time t appears explicitly. Thus the similarity condition
requires that
must be a constant
quantity. We therefore assume that
![]() |
(13) |
Integrating (13), we obtain
![]() |
(14) |
where the constant of integration L is determined through the condition that ?=L when t=0.. Here C=0 implies that ?=L represents the length scale for steady flow and C≠0 that is, ? represents the length scale for unsteady flow. Let us now consider a class of solutions for which C=2 and hence the length scale ? from Eq. (14) becomes
![]() |
(15) |
which exactly corresponds to the usual scaling factor considered for various unsteady boundary layer flows Schlichting (1968). Since ? is a scaling factor as well as a similarity parameter, any other values of C in Eq. (14) would not change the nature of the solution except that the scale would be different. Finally, introducing Eq. (15) in Eqs (9-12) respectively, we have the following dimensionless ordinary non-linear differential equations
![]() |
(16) |
![]() |
(17) |
![]() |
(18) |
![]() |
(19) |
The boundary conditions (7) then become
![]() |
(20) |
In all the above equations primes denote the differentiation with respect to ?. Equation (16) to (19) are solved numerically under the boundary conditions (20) using Nachtsheim-Swigert iteration technique.
IV. NUMERICAL SOLUTIONS
Equations (15)-(18) are solved numerically under the boundary conditions (19) using Nachtsheim and Swigert (1965) iteration technique. Thus adopting this numerical technique, a computer program was set up for the solutions of the basic non-linear differential equations of our problem where the integration technique was adopted as the six ordered Range-Kutta method of integration. The results of this integration are then displayed graphically in the form of velocity, temperature and concentration profiles in Figs. (2)-(10). In the process of integration, the local skin-friction coefficients, the local rates of heat and mass transfer to the surface, which are of chief physical interest are also calculated out. From these solutions we may obtain derivatives of the G and F' functions at the disk surface (?=0) which are required for the calculation of the tangential shear stress τt and the radial shear stress τr. To find the tangential shear stress τt and the radial shear stress τr, we apply the Newtonian formulae
![]() |
(20) |
and
![]() |
(21) |
The heat flux (qw) and the mass flux (Mw) at the surface are given by
![]() |
(22) |
![]() |
(23) |
Fig. 2. Effects
of M on the radial velocity profiles for
Pr=0.71, R=0.5, Sc=0.6,
Df=0.2 and So=0.5.
Fig. 3. Effects
of M on the tangential velocity profiles for
Pr=0.71, R=0.5, Sc=0.6,
Df=0.2 and So=0.5.
Fig. 4. Effects
of M on the axial velocity profiles for
Pr=0.71, R=0.5,
Sc=0.6, Df=0.2
and So=0.5.
Fig. 5. Effects
of R on the radial velocity profiles for
Pr=0.71, M=0.5,
Sc=0.6, Df=0.2
and So=0.5.
Fig. 6. Effects
of R on the tangential velocity profiles for
Pr=0.71, M=0.5,
Sc=0.6, Df=0.2
and So=0.5.
Fig. 7. Effects
of R on the axial velocity profiles for
Pr=0.71, M=0.5,
Sc=0.6, Df=0.2
and So=0.5.
Fig. 8. Effects
of Df on the temperature profiles for
Pr=0.71, M=0.5,
Sc=0.6, R=0.5 and
So=0.5.
Fig. 9. Effects
of So on the concentration profiles for
Pr=0.71, M=0.5,
Sc=0.6, R=0.5 and
Df =0.2.
These above coefficients are then obtained from the procedure of the numerical computations and are sorted in Table 1. In the above equations prime denotes derivative with respect to similarity variable ?.
Table 1(a): Numerical values of the
radial and tangential skin frictions and the rate of heat and mass
transfer coefficients obtained for Pr=0.71,
Sc=0.6, M=0.5, Df=0.2
and So=0.5.
Table 1(b): Numerical values of the
radial and tangential skin frictions obtained for
Pr=0.71, Sc=0.6, R=0.5,
Df=0.2 and So=0.5.
Table 1(c): Numerical values of the
rate of heat and mass transfer coefficients obtained for
Pr=0.71, Sc=0.6 and
R=0.5.
V. RESULTS AND DISCUSSIONS
Considered here involves a number of parameters on the basis of which a wide range of numerical results have been derived. Of this results a small section are presented here for brevity. The numerical results of the velocity, temperature and concentration profiles are shown in Figs.2-4.
The parameters entering into the fluid flow are Magnetic parameter M, Rotational parameter R, Dufour number Df and Soret number So. It is, therefore, pertinent to inquire the effects of variation of each of them when the others are kept constant. The numerical results are thus presented for different values of M, R, Df and So for Prandlt number Pr=0.71 and Schmidt number Sc=0.6 are considered.
The representative velocity, temperature and concentration profiles and the values of the physically important parameters, i.e. the local shear stress, the local rates of heat and mass transfer, are illustrated for uniform wall temperature and species concentration in Figs. 2-9 and in Table 1(a, b, c). In Figs. 2-4, the radial, tangential and axial velocity profiles are shown respectively for different values M (0, 1 and 2). From these figures it is observed that the increase in the magnetic field leads to the decrease in the velocity field indicating that the magnetic field retards the velocity field. In Fig. 5 - Fig. 7, we have plotted the radial, tangential and axial velocity for air (Pr=0.71) at 20ºC and at one atmospheric pressure showing the effect of variation of Rotational parameter R. It is thus apparent that the rotational parameter R has increasing effects on radial and axial velocity profiles but decreasing effects on tangential velocity profiles.
The parameter Df (Dufour number) does not enter directly into the momentum and mass equations. Thus the effect of Dufour number on velocity and mass profiles is not apparent. Figure 8 shows the variation of temperature profiles for different values of Df.
The parameter Df has marked effects on the temperature profiles. It is observed that the temperature profiles increase with the increasing values of Df. It is also observed from this figure that when Df = 1.0, that is, when the ratio between temperature and concentration gradient is very small the temperature profile shows its usual trend of gradual decay. As Dufour number Df becomes large the profiles overshoot the uniform temperature close to the boundary.
The parameter So (Soret number) does not enter directly into the momentum and energy equations. Thus the effect of Soret number on velocity and temperature profiles is not apparent. Figure 9 shows the variation of concentration profiles for different values of So. The parameter So has marked effects on the concentration profiles. It is observed that the concentration profiles increase with the increasing values of So. It is also observed from this figure that when So=1.0, that is, when the ratio between concentration and temperature gradient is very small the concentration profile shows its usual trend of gradual decay. As Soret number So becomes large the profiles overshoot the uniform concentration close to the boundary.
The radial and tangential skin-frictions, the heat transfer and the mass transfer coefficients are tabulated in Table 1 for various values of M, R, Df and So. We observe from Table 1(a) that the radial skin-friction, heat and mass transfer coefficients increase while the tangential skin-friction decreases with the increase in the rotational parameter R. It has been found from Table 1(b) that the radial skin-friction and tangential skin friction decrease with the increase in the magnetic parameter M. The Table 1(c) shows that the heat transfer coefficient decreases with the increasing values of the Dufour number Df . It is also observed from Table 1(c) that as thermal diffusion (S0) intensifies mass transfer coefficient decreases.
In order to highlight the validity of the numerical computations adapted in the present investigation, some of our results have been compared with those of Sparrow and Gregg (1959) in Table 2. The comparison shows excellent agreement, hence an encouragement for the use of the present numerical computations.
Table 2: Comparison between the
(analytical) results of Sparrow and Gregg (1959) and present work.
Heat transfer coefficients -T'(0) for
Sc=0.6, Df=0,
So=0, M=C=0 and R=1.
VI. CONCLUSIONS
In this paper, an unsteady three dimensional MHD convective heat and mass transfer flow in an incompressible fluid due to an infinite rotating disk is studied. The Nachtsheim and Swigert (1965) iteration technique based on sixth-order Range-Kutta and shooting methods has been employed to complete the integration of the resulting solutions.
The following conclusions can be drawn as a result of the computations:
REFERENCES
1. Alam, M.M. and M.A. Sattar, "Unsteady MHD Free Convection and Mass Transfer Flow in a Rotating System with Thermal Diffusion", Journal of Energy Heat and Mass Transfer, 20, 77-87 (1998).
2. Attia, H.A., "Unsteady MHD flow near a rotating porous disk with uniform suction or injection", Fluid Dynamics Research, 23, 283-290 (1998).
3. Benton, E.R., "On the flow due to a rotating disc", J. Fluid Mech. Part 4, 24, 781-800 (1966).
4. Cochran, W.G., "The flow due to a rotating disc", Proc. Cambridge Phil. Soci., 30, 365-375 (1934).
5. Eckert, E.R.G. and R.M. Drake, Analysis of Heat and Mass Transfer, McGraw-Hill Book Co., NewYork (1972).
6. Hassan, A.L.A. and H.A. Attia, "Flow due to a rotating disk with Hall effect", Physics Letters A, 228, 246-290 (1997).
7. Hayat, T., S. Saif and Z. Abbas, "Mixed convection flow of a micropolar fluid over a non-linearly stretching sheet", Physics letters A, 372,637-647 (2008a).
8. Hayat, T., Z. Abbas and M. Sajid, "Heat and mass transfer analysis on the flow of a second grade fluid in the presence of chemical reaction", Physics letters A, 372, 2400-2408 (2008b).
9. Hayat, T., S. Saif and Z. Abbas, "The influence of heat transfer in an MHD second grade fluid film over an unsteady stretching sheet", Physics letters A, 372, 5037-5045 (2008c).
10. Hayat, T., K. Maqbool and M. Khan, "Hall and heat transfer effects on the steady flow of generalized Burgers' fluid induced by a sudden pull of eccentric rotating disks", Nonlinear Dynamics, 51, 267-276 (2008d).
11. Jha, B.K. and A.K. Singh, "Soret Effects Free Convection and Mass Transfer Flow over an Infinite Vertical Moving Plate", Astrophysics and Space Science, 173, 251-255 (1990).
12. Kafoussias, N.G., "MHD Thermal-Diffusion Effects on Free-Convective and Mass Transfer Flow over an Infinite Vertical Moving Plate", Astrophysics and Space Science, 192, 11-19 (1992).
13. Kelson, N. and A. Desseaux, "Notes on porous rotating disk flow", ANZIAM J., 42(E), C837-C855 (2000).
14. Kuiken, H.K., "The effect of normal blowing on the flow near a rotating disk of infinite extent", J. Fluid Mech., 47, 789-798 (1971).
15. Maleque, Kh.A. and M.A. Sattar, "Similarity Solution of MHD Free-Convective and Mass Transfer Flow Over an Vertical Porous Plate with Thermal Diffusion Effects", The Aiub Journal of Science and Engineering (AJSE), 1, 44-55 (2002).
16. Maleque, Kh.A. and M.A. Sattar, "Transient Convective Flow Due to a Rotating Disc with Magnetic Field and Heat Absorption Effects", Journal of Energy, Heat and Mass Transfer, 25, 279-291 (2003).
17. Maleque, Kh.A. and M.A. Sattar, "The effects of variable properties and Hall current on steady MHD compressible laminar convective fluid flow due to a porous rotating disc", International Journal of Heat and Mass Transfer, 48, 4963-4972 (2005a).
18. Maleque, Kh.A. and M. A. Sattar, "The effects of variable properties on steady laminar convective flow due to a porous rotating disc", ASME Journal of Heat Transfer, Technical Brief, 127, 1406-1409 (2005b).
19. Nachtsheim, P.R. and P. Swigert, Satisfaction of asymptotic boundary conditions in numerical solution of system of nonlinear of boundary layer type, NASA TN-D3004 (1965).
20. Sajid, M. and T. Hayat, "Influence of thermal radiation on the boundary layer flow due to an exponentially stretching sheet", International Communications in Heat and Mass Transfer, 35, 347-356 (2008).
21. Sattar, M.A., and M.M. Alam, "Analytical Solution of the Free Convection and Mass Transfer Flow with Thermal Diffusion", Dhaka University Journal of Science, 49, 95-104 (2001).
22. Schlichting, H., Boundary Layer Theory, McGraw Hill Book Co., New York (1968).
23. Sparrow, E.M. and J.L. Gregg, "Heat transfer from a rotating disk to fluids of any Prandtl number", ASME J. Heat Transfer, C., 8, 249-251 (1959).
24. von Karman, T., "Ueber laminare und turbulente reibung", ZAMM, 1, 233-255 (1921).
25. Wagner, C., "Heat transfer from a rotating disk to ambient air", Journal of Applied Physics, 19, 837-839 (1948).
Received:
July 22, 2008.
Accepted: June 10, 2009.
Recommended by Subject Editor Walter Ambrosini.