# Non-Darcy mixed convective heat and mass transfer flow of a viscous electrically conducting fluid through a porous medium in a circular annulus in the presence of temperature gradient heat sources with Soret and Dufour effects--a finite element study.

IntroductionConvection flows driven by temperature and concentration differences have been studied extensively in the past and various extensions of the problems have been reported in the literature with both temperature and concentration interacting simultaneously, the convection have become quite complex. Bejan and Khair [3] have investigated the vertical free convection flow embedded in a porous medium resulting from the combined heat and mass transfer. Jang and Chang [8] have used an implicit finite difference method to study the buoyancy induced inclined boundary layer in a porous medium resulting from the combined heat and mass buoyancy effects.

Heat transfers in thermal insulation within vertical cylindrical annuli provide us insight into the mechanism of energy transport and enable engineers to use insulation more efficiently. In particular design engineers require relationship between heat transfer, geometry and boundary conditions which can be utilized cost-benefit analysis to determine the amount of insulation that will yield the maximum investment. An understanding of convective heat transfer in porous annuli is essential for its applications in packed-bed catalytic reactors, Geophysics, thermal insulation, design of regenerative heat exchangers, Geological disposal of high-level nuclear waste, petroleum resources and many other uses. Free convection in a vertical porous annulus has been extensively studied by Prasad [11], Prasad and Kulacki [12] and Prasad et al., [13] both theoretically and experimentally. Convection through annular regions under steady conditions has also been discussed with the two cylindrical surfaces kept at different temperatures [12]. This work has been extended in temperature dependent convection flow [9] as well as convection flow through horizontal porous channel whose inner surface is maintained at constant temperature, while the other surface is maintained at circumferentially varying sinusoidal temperature [17].

The applications of electromagnetic fields in controlling the heat transfer as in aerodynamic heating leads to the study of Magneto hydrodynamic heat transfer. This MHD heat transfer has gained significance owing to recent advancement of space technology. The MHD heat transfer can be divided into two sections. One contains problems in which the heating is an incidental by product of the electromagnetic fields as in MHD generators and pumps etc. and the second consist of problems in which the primary use of electromagnetic fields is to control the heat transfer(5). With fuel crisis deepening all over the world, there is a great concern to utilize the enormous power beneath the earth's crust in the geothermal region (10). Liquid in the geothermal region is an electrically conducting liquid because of high temperature and that they undergo the influence of magnetic field.

In many industrial applications of transient free convection flow problems, there occurs a heat source or a sink which is either a constant or temperature gradient or temperature dependent heat source. This heat source occurs in the form of a coil or a battery. Gokhale and Behnaz-Farman [7] analyzed Transient free convection flow of a n incompressible fluid past an isothermal plate with temperature gradient dependent heat sources. Implicit finite difference scheme which is unconditionally stable has been used to solve the governing partial differential equations of the flow. Transient temperature and velocity profiles are plotted to show the effect of heat source. Muthukumara swamy [9] et al has analyzed the radiation effect on moving vertical plate with variable temperature and mass diffusion. Sreevani [16] has analyzed the Soret effect on convective heat and mass transfer flow of a viscous fluid in a cylindrical annulus with heat generating sources. Sivaiah [15] has discussed the convective heat and mass transfer flow in a circular duct with Soret effect.

Literature suggests that the effect of viscous dissipation on heat transfer has been studied for different geometries. Brinkman [5] have studied the viscous dissipation effect on natural convection in horizontal cylinder embedded in porous medium. Their study showed that the viscous dissipation effect might not be neglected. Saffman [14] have studied the viscous dissipation effect on natural convection in a porous cavity and found that the heat transfer rate at hot surface decreases with increase of viscous dissipation parameter. Thermal radiation plays a significant role in the overall surface hear transfer where convective heat transfer is small. Verschoor et al [18] have studied the effect of viscous dissipation and radiation on unsteady magneto hydro dynamic free convection flow fast vertical plate in porous medium. They found that the temperature profile increases when viscous dissipation increases. A good amount of work has been done to understand natural convection in porous cavity. In spite of endeavor efforts to study heat transfer in porous cavity, the combined effect of viscous dissipation and radiation on porous medium filled inside a square cavity has not received attention.

The Soret and Dufour effects have garnered considerable interest in both Newtonian and non-Newtonian convective heat and mass transfer. Such effects are significant when density differences exist in the flow regime. Soret and Dufour effects are important for intermediate molecular weight gases in coupled heat and mass transfer in binary systems, often encountered in chemical process engineering and also in high-speed aerodynamics. Soret and Dufour effects are also critical in various porous flow regimes occurring in chemical and geophysical systems. There are few studies about the Soret and Dufour effects in a Darcy or non-Darcy porous medium. Anghel et al [1] has examined the composite Soret and Dufour effects on free convective heat and mass transfer in a Darcian porous medium with Soret and Dufour effects. Very recently, Barletta. A, Lazzari.S, and others [2] have studied on Mixed convection with heating effects in a vertical porous annulus with a radially varying magnetic field. Emmunuel Osalusi, Jonathan Side, Robert Harris [6] have discussed 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.

The Weighted residual method is of the generalization of the Ritz-variational method where in we seek an approximate solution in the form of linear combination of suitable approximation functions. The parameters in the linear combination as determined by setting integral of a weighted residual of the approximation over the domain zero. A comprehensive description of a weighted residual method. In many situations the Galerkin method which is one of the important weighted residual methods is equivalent to the Ritz method for solving variational problems. The finite element method is piece--wise applications of weighted residual method in which the Ritz-Galerkin type methods are employed over each element of the domain. The finite element method was initially developed as an adhoc engineering procedure for constructing matrix solutions to stress and displacement calculations in structural analysis. Very few fluid dynamic problems can be expressed in a variational form. Consequently most of the finite element applications in fluid dynamics have used in Galerkin finite element formulation. The Galerkin finite element method has two important futures. Firstly the approximate solution is written directly as a linear combination of approximating functions in terms of the nodal unknowns. Secondly the approximating functions or the shape functions are chosen exclusively from low order piecewise polynomials restricted to contiguous elements.

In this paper we discuss the mixed convective viscous dissipative flow through a porous medium in a circular cylindrical annulus with Thermal-Diffusion and Diffusion-Thermo effects in the presence of temperature gradient heat source, where the inner wall is maintained constant temperature while the outer wall is maintained constant heat flux and the concentration is constant on the both walls. The coupled momentum, energy and diffusion equations are solved by using finite element analysis with quadratic interpolation polynomials. The effect of temperature gradient heat sources on the flow and heat transfer characteristics are analyzed. The stress, rate of heat transfer and the rate of mass transfer are discussed numerically for different variations of the governing parameters.

Formulation of the Problem

We consider free and force convective flow of a viscous, electrically conducting fluid through a porous medium in a circular cylindrical annulus with Thermal-Diffusion and Diffusion-Thermo effects in the presence of temperature gradient heat source, whose inner wall is maintained at a constant temperature and the outer wall is maintained constant heat flux. Also the concentration is constant on the both walls. A uniform radial magnetic field is applied on the flow. The flow, temperature and concentration in the fluid are assumed to be fully developed. Both the fluid and porous region have constant physical properties and the flow is a mixed convection flow taking place under thermal and molecular buoyancies and uniform axial pressure gradient. The boussenissque approximation is invoked so that the density variation is confined to the thermal and molecular buoyancy forces. The Brinkman-Forchhimer-Extended Darcy model which accounts for the inertia and boundary effects has been used for the momentum equation in the porous region. In the momentum, energy and diffusion are coupled and non-linear. Also the flow in is unidirectional along the axial cylindrical annulus. Making use of the above assumptions the governing equations are Equation of linear momentum

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

Equation of Energy

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

Equation of diffusion

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

Equation of state

[rho] - [[rho].sub.0] = -[beta][[rho].sub.0](T - [T.sub.0]) - [[beta].sup.*][[rho].sub.0](C - [C.sub.0]) (4)

Where u is the axial velocity in the porous region, T & C are the temperature and concentrations of the fluid, k is the permeability of porous medium, F is a function that depends on Reynolds number and the microstructure of the porous medium and [D.sub.1] is the Molecular diffusivity, [D.sub.m] is the coefficient of mass diffusitivity, [T.sub.m] is the mean fluid temperature, [K.sub.t] is the thermal diffusion, [C.sub.s] is the concentration susceptibility, [C.sub.p] is the specific heat, [rho] is density, g is gravity, [beta] is the coefficient of thermal expansion, [[beta].sup.*] is the coefficient of volume expansion. [sigma] is the electrical conductivity, [[mu].sub.e] is the magnetic permeability.

u = 0, T = [T.sub.i], C = [C.sub.i] at r = a (6)

u = 0, [partial derivative]T/[partial derivative]r = [Q.sub.1], C = Co at r = a + s (7)

The axial temperature gradient [partial derivative]T/[partial derivative]r and concentration gradient [partial derivative]C/[partial derivative]z are assumed to be constant say A and B respectively.

Using equations (2.5) and (2.6) equation (2.2) reduces to

[[rho].sub.0][C.sub.p]uA = [lambda]([T.sub.rr] + [1/r] [T.sub.r]) + [[D.sub.m][K.sub.t]/[C.sub.s][C.sub.P]] ([C.sub.rr] + [1/r] [C.sub.r]) (8)

[[rho].sub.0][C.sub.P]uB = [lambda]([C.sub.rr] + [1/r] [C.sub.r]) + [[D.sub.m][K.sub.t]/[C.sub.s][C.sub.P]] ([T.sub.rr] + [1/r] [T.sub.r]) (9)

We now define the following non-dimensional variables

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)

Introducing these non-dimensional variables, the governing equations in the non-dimensional form are (on removing the stars)

[[d.sup.2]u/d[r.sup.2] + [1/r][du/dr]] = P + [delta]([D.sup.-1] + [[M.sup.2]/[r.sup.2]])u + [[delta].sup.2][LAMBDA] [u.sup.2] - [delta]G[theta] (11)

[[d.sup.2][theta]/d[r.sup.2] + [1/r][d[theta]/dr] + [[alpha]/r] [d[theta]/dr]] = [P.sub.r][N.sub.T]u + Du[N.sub.t] ([d.sup.2]C/d[r.sup.2] + [1/r][d[theta]/dr]) + [P.sub.r][E.sub.c] [(du/dr).sup.2] (12)

[d[r.sup.2]C/d[r.sup.2] + [1/r][dC/dr]] = Sc[N.sub.c]u + ScSr ([d.sup.2][theta]/d[r.sup.2] + [1/r][d[theta]/dr]) (13)

where

[LAMBDA] = F[D.sup.-1] (Forchheimer number)

[P.sub.r] = [mu]/[C.sub.p]/[lambda] (Prandtl number)

G = g[beta]([T.sub.1] - [T.sub.0])[a.sup.3]/[v.sup.2] (Grashof number)

[D.sup.-1] = [a.sup.2]/k (Inverse Darcy parameter)

[N.sub.t] = Aa/[T.sub.1] - [T.sub.0] (Temperature gradient)

[N.sub.c] = Ba/[C.sub.1] - [C.sub.0] (Non-dimensional concentration gradient)

Du = ([D.sub.m][K.sub.t][DELTA][ca.sup.2]/[C.sub.s][C.sub.p][DELTA]T][lambda]) (Dufour Number)

Sc = v/[D.sub.1] (Schmidt number)

Sr = ([D.sub.m][K.sub.t][DLETA]T/v[T.sub.m][DELTA]C) (Soret number)

[alpha] = Q[L.sup.2]/[lambda][C.sub.p] (Heat source parameter)

[M.sup.2] = [sigma][[mu].sup.2.sub.e][H.sup.2.sub.0][a.sup.2]/[[gamma].sup.2] (Hartman number)

[E.sub.c] = [v.sup.2]/[a.sup.2](Aa)[c.sub.p] (Eckert number)

With the corresponding boundary conditions as:

u = 0, [theta] = 0, C = 1 at r = 1 (14)

u = 0, [partial derivative][theta]/[partial derivative]r = [Q.sub.1], C = 0 at r = 1 + s (15)

Numerical Analysis

The finite element method has been implemented to obtain numerical solutions of equations (11) to (13) under boundary conditions (14) and (15). This technique is extremely efficient and allows robust solutions of complex coupled, nonlinear multiple degree differential equation systems. The fundamental steps comprising the method are now summarized:

Phase 1] Discretization of the domain into elements

Phase 2] Derivation of element equations

Phase 3] Assembly of Element Equations

Phase 4] Imposition of boundary conditions

Phase 5] Solution of assembled equations

The shear stress are evaluated on the cylinder using the formula

[tau] = [(du/dr).sub.r=1,1+s]

The rate of heat transfer (Nusselt number) are evaluated on the cylinder using the formula

Nu = [(d[theta]/dr).sub.r=1]

The rate of mass transfer (Sherwood Number) is evaluated using the formula

Sh = -[(dC/dr).sub.r=1,1+s]

Discussion of the Numerical Results

In this analysis we investigate thermo-Diffusion and Diffusion-Thermo effects on convective heat and mass transfer flow of a viscous conducting fluid through a porous medium in circular annulus in the presence of temperature gradient dependent heat source with viscous dissipation. The inner cylinder is maintained at constant temperature and the outer wall is maintained constant heat flux while the concentration is maintained constant on both the cylinders. The axial flow is in vertically downword direction, u > 0 indicates a reversal flow. The velocity, temperature and concentration distributions are shown in figures 1-27 for different values of the parameters G, [D.sup.-1], M, Sc, Sr, Du, N, [alpha] and Ec.

The variation of u with Grashof number G shows that the axial flow enhances with increase in G and the region of reversal flow enhances with increase in G (fig.1). With respect to the variation of u with [D.sup.-1] we find that lesser the permeability of porous medium smaller the magnitude of u and for further lowering of the permeability larger th-1e magnitude of u in the entire flow region, the region of reversal flow shrinks with [D.sup.-1] [less than or equal to] 2 X [10.sup.3] and enhances with higher [D.sup.-1] [greater than or equal to] 3 X [10.sup.3] (fig.2). From fig.3 we find that higher the Lorentz force larger the velocity in the flow region. Also the region of reversal flow enhances with increase in M. Fig.4 represents the variation of u with Sc. We notice that lesser the molecular diffusitivity smaller [absolute value of u] and for further lowering of molecular diffusitivity it experiences a depreciation in the entire flow region and it attains maximum at r = 1.5.

[FIGURE 1 OMITTED]

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

The variation of u with Soret parameter Sr shows that the velocity experiences an enhancement with increase in Sr [less than or equal to] 0.8 and for further increase in Sr [greater than or equal to] 1 it depreciates in its magnitude (fig.5). From fig.6 we observe that the region of reversal velocity enlarges with increase in Du and [absolute value of u] enhances with Du. The variation of u with N shows that when the molecular buoyancy force dominates over the thermal buoyancy force the actual axial velocity experiences a depreciation when the buoyancy forces act in the same direction while for the forces acting in the opposite directions it experiences an enhancement in the flow region (fig.7). The influence of heat source parameter [alpha] on u is shown in fig.8.

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

[FIGURE 7 OMITTED]

[FIGURE 8 OMITTED]

An increase in [alpha] < 0 enhances the actual axial velocity u while it depreciates with [alpha] > 0. This shows that in the presence of the temperature heat source depreciates the velocity in the flow region with maximum in the mid region. The influence of dissipative effect on u is shown in fig.9. We conclude that the axial velocity u experiences an enhancement with Ec.

[FIGURE 9 OMITTED]

The non-dimensional temperature ([theta]) is shown in fig 10-18 for different values of the parameters. It is found that the non-dimensional temperature gradually increases from its prescribed value 0 on r = 1 to attain its prescribed value 1 at r = 2. An increase in G enhances the temperature (fig.10). The variation of [theta] with [D.sup.-1] shows that lesser the permeability of porous medium larger the actual temperature in the flow region (fig.11). From fig.12 we find that lesser the Lorentz force larger the temperature. With respect to Sc, we notice that lesser the molecular diffusitivity smaller the temperature in the flow region (fig.13).

[FIGURE 10 OMITTED]

[FIGURE 11 OMITTED]

[FIGURE 12 OMITTED]

[FIGURE 13 OMITTED]

An increase in Soret parameter Sr results in a depreciation in the actual temperature in the region (fig.14). The variation of [theta] with Dufour parameter Du shows that the actual temperature enhances gradually with increase in Du (fig.15). When the molecular buoyancy force dominates over the thermal buoyancy force the actual temperature decreases irrespective of the directions of the buoyancy forces (fig.16). The influence of temperature gradient heat source parameter [alpha] on [theta] is shown in fig.17. It found that the temperature is negative for [alpha] = 0, [alpha] < 0 and positive for [alpha] > 0. The actual temperature experiences a depreciation with increase in the strength of the heat sources. The variation of [theta] with Eckert number Ec is shown in fig.18. We find that the actual temperature enhances with increase in Ec.

[FIGURE 14 OMITTED]

[FIGURE 15 OMITTED]

[FIGURE 16 OMITTED]

[FIGURE 17 OMITTED]

[FIGURE 18 OMITTED]

The non-dimensional concentration ([phi]) is shown in fig 19-27 for different values of the parameters G, [D.sup.-1], M, Sc, Sr, Du, N, [alpha] and Ec. It is found that the non-dimensional concentration gradually increases from its prescribed value 0 on r = 1 and attain its prescribed value 1 at r = 2. Fig.19. shows the variation of [phi] with G. It is noticed that the concentration depreciates with increase in the Grashof number G. The variation of [phi] with [D.sup.-1] shows that lesser the permeability of porous medium higher the actual concentration in the flow region and for further lowering of the permeability smaller the actual concentration in the flow region (fig.20). From fig.21 we find that higher the Lorentz force larger the concentration in the flow region. Lesser the molecular diffusitivity smaller the concentration in the flow field (fig.22). An increase in the Soret parameter Sr enhances the concentration everywhere in the flow region (fig.23). The variation of [phi] with Dufour parameter Du shows that the concentration experiences a marginal depreciation in the flow region (fig.24). The variation of [phi] with N shows that when the molecular buoyancy force dominates over the thermal buoyancy force the actual concentration experiences an enhancement when the buoyancy forces act in the same direction while for the forces acting in the opposite directions it experiences a depreciation in the flow region (fig.25). From fig.26, we observe that [phi] experiences a marginal enhancement with [alpha] > 0 and a depreciation with [alpha] < 0. The inclusion of the dissipation in the flow enhances the concentration in the flow region (fig.27).

[FIGURE 19 OMITTED]

[FIGURE 20 OMITTED]

[FIGURE 21 OMITTED]

[FIGURE 22 OMITTED]

[FIGURE 23 OMITTED]

[FIGURE 24 OMITTED]

[FIGURE 25 OMITTED]

[FIGURE 26 OMITTED]

[FIGURE 27 OMITTED]

References

[1] Angel M, Takhar HS, and Pop I (2000): Dofour and Soret effects on free convection boundary layer over a vertical surface embedded in a porous medium. Studia universities-Bolyai, Mathematica XLV, pp. 11-21.

[2] Barletta.A, Lazzari.S,(2008): Mixed convection with heating effects in a vertical porous annulus with a radially varying magnetic field.

[3] Bejan, A and Khair, K.R: Heat and mass transfer by natural convection in a pourous medium., Int.J. Heat and mass transfer, V.28.PP.908-818(1985).

[4] Brinkman H.C (1948): A Calculation of the viscous force eternal by a flowing fluid on a dense swarm of particles. Appl.Science Research, Ala, p81.

[5] Chandra sekhar. S: Hydrodynamic and Hydro magnetic stability, Clarandon press, oxford (1961).

[6] Emmunuel Osalusi, Jonathan Side, Robert 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: Int.Communications in heat and mass transfer,Vol.35,PP.908-915(2008).

[7] Gokhale, M.Y, Behnaz-Farnam: "Transient free convection flow on an isothermal plate with temperature dependent heat sources", International review of Pure and Applied Mathematics, Vol.3, No.1, pp.129-136(2007).

[8] Jang J.Y and Chang W.J: The flow and vertex instability of horizontal natural convection in a porous medium resulting from combined heat and mass buoyancy effects. Int. Heat and Mass Transfer, V. 31, p.769-777 (1987).

[9] Muthukumara Swamy, Maheswari, J, Pandurangan, J: "Study of MHD and Radiation effects on moving vertical plate with variable temperature and mass diffusion", International Review of Pure and Applied Mathematics, Vol.3, No.1, pp.95-103(2007).

[10] Nanda, R.S and Mohan, M: Proc. Ind Acad.Sci., V.876 A,No.5,P.147(1978).

[11] Prasad, V: Natural convection in a vertical porous annulus, Int. J. Heat and Mass Transfer, V. 27, p 207-219 (1984).

[12] Prasad, V and Kulacki, F. A: Natural convection in porous media bounded by short concentric vertical cylinders, ASME J. Heat Transfer, V. 107, p 147-154 (1985).

[13] Prasad, V, Kulacki, F.A and Keyhani, M: Natural convection in porous media, J. fluid mech., V. 150, p 89-119 (1985).

[14] Saffman .P.G (1971): on the boundary conditions at the free surface of a porous medium. Stud.Appl.Maths, V.2, p 93.

[15] Sivaiah (2004): Thermo-diffusion effect on convective heat and mass transfer flow through a porous medium in ducts. Ph.D thesis, S.K. University. Anantapur.

[16] Sreevani, M: Mixed convection heat and mass transfer through a porous medium in channels with dissipative effects, Ph.D thesis, S, K.University, Anantpur, India(2003).

[17] Vasseur, P, Nguyen, T.H., Robillard and Thi. V.K.T: Int. Heat Mass Transfer, V. 27, p 337 (1984).

[18] Verschoor et al: Int. Heat and Mass Transfer, V. 31(1992).

P. Sudarsan Reddy (1), G. Srinivas (2), P. Sreenivasa Rao (3), DRV. Prasada Rao (4)

(1) Asst. Professor, Dept of Mathematics, RGM Engg. College, Nandyal, JNTU Anantapur, India Email: suda1983@gmail.com

(2) Assoc. Professor, Dept of Mathematics, Vignan Engg. College, JNTU Hyderabad, India

(3) Assoc. professor, Dept. of Physics, Jyothismathi Engg.College, Karimnagar, JNTU Kakinada, India.

(4) Professor, Dept of Mathematics, SK University, Anantapur. India