Research Article - (2017) Volume 7, Issue 2
In this paper, entropy generation of double-diffusive natural convection in a 2D dimensional enclosure with magnetic and Dufour effects has been numerically performed. Dirichlet boundary conditions for temperature and solute concentration are applied to the two vertical walls of the enclosure; wheras the two horizontal walls are adiabatic and insulated. The governing equations of continuity, momentum, energy and concentration are numerically solved by using a Control Volume Fined Elements Method, CVFEM of Patankar. The governing parameters of the problem are the thermal Grashof number (GrT), the buoyancy ratio (N), the Hartmann number (Ha), the Dufour parameter (Du) and the Prandtl number (Pr). The obtained results were presented graphically via the velocity field components, temperature and concentration distributions, entropy genertion rate behaviour and by isotherms, streamlines and isentropic lines maps. The average Nusselt and Sherwood numbers are also derived and discussed numerically. The investigated results showed that the flow field and then entropy generation are notably influenced by the considering parameters.
Keywords: MHD double-diffusive convection; Square cavity; Entropy generation; Dufour effect; CVFEM.
B | Magnetic Field (T) |
C | Dimensionless Concentration |
c | Concentration (mol·m^{−3}) |
C_{p} | Specific Heat (J∙Kg∙K^{−1}) |
D | Molecular Diffusivity (m^{2}∙s^{-1}) |
D_{F} | Dufour coefficient (m^{2}∙s^{−1}∙K^{-1}) |
Du | Dufour parameter |
E | Electric Field (V·m^{−1}) |
F | Electromagnetic Force (N) |
g | Gravitational Acceleration (m·s^{−2}) |
Gr_{T } | Thermal Grashof Number |
Grc | Solutal Grashof Number |
Ha | Hartmann Number |
J | Electric Current Density (A·m^{−2}) |
J_{k } | Diffusion Flux (k=u, v, T) |
L | Length of the Cavity (m) |
Le | Lewis number |
N | Buoyancy Ratio |
Nu | Nusselt Number |
P | Dimensionless Pressure |
P | Pressure (N·m^{−2}) |
Pr | Prandtl Number |
Ra_{T } | Thermal Rayleigh Number, RaT=GrT.Pr |
Sc | Schmidt Number |
Sh | Sherwood Number |
S | Entropy Generation Rate (J·s^{−1}·K^{−1}·m^{−3}) |
T | Dimensionless Temperature |
t, τ | Time (s) |
u, v | Velocity Components in X And Y Directions (m•s^{−1}) |
U, V | Dimensionless Velocity Components |
U_{d } | Dimensionless Velocity Vector |
α_{T } | Thermal Diffusivity (m^{2}•s^{−1}) |
ø | Magnetic Inclination Angle (°) |
β_{T } | Thermal Expansion Coefficient (K^{−1}) |
βc | Solutal Expansion Coefficient (m^{3}•mol^{-1}) |
λ | Thermal Conductivity (J·m^{−1}·s^{−1}·K^{−1}) |
λi | Irreversibility Coefficient (i=1,2,3,4) |
μ | Dynamic Viscosity (kg·m^{−1}·s^{−1}) |
μsp | Chemical Potential Of Species (J•kg ^{−1}) |
W | Thermodynamic Parameter W=[(C_{0}. θ_{0}/CP].(∂µ_{sp}/∂C)θ_{0},C_{0} (K^{2}) |
x, y | Cartesian Coordinates (m) |
X, Y | Dimensionless Cartesian Coordinates |
Greek letters | |
ν | Kinematic Viscosity (m^{2}·s^{−1}) |
ρ | Mass Density (kg·m^{−3}) |
σ | Dimensionless Entropy Generation |
σ_{e} | Electric Conductivity (Ω^{−1}·m^{−1}) |
θ | Temperature (K) |
θ_{0} | Mean Temperature (K) |
Subscripts | |
a | Dimensionless |
Mag | Magnetic |
c | Cold Wall |
h | Hot Wall |
In the past few years, the problems of double-diffusive convective and heat transfer flows through a cavity under the influence of a magnetic field have been attracted the attention of a number of researchers because of their several applications in many engineering problems, such as in thermal energy storage systems, industrial magnetic devices, petroleum industry, polymer technology and glass fiber production. The interest for analysing the double-diffusive convective flows induced by the cross thermo diffusion effects has surged in view of its importance in many engineering problems such as in insulation, the underground disposal of nuclear wastes and drying processes.
Trevisan and Bejan [1] reported a fundamental numerically and analytically study of scale analysis relative to heat and mass transfer within a vertical slot filled with a porous medium submitted to horizontal combined and pure temperature and concentration gradients. The more observed effect was due to heat transfer and then to mass transfer. They obtained that the magnetic effect is more pronounced than friction one. The magnetic field parameter suppresses the flow in the cavity and this lead to a decrease of entropy generation. At local level, results show that entropy generation lines are localized on lower heated and upper cooled regions of the active walls.
Other related numerical heat and mass transfer studies of a steady free convection in MHD micropolar electrically conducting fluid flow on a vertical plate in the presence of Soret and Dufour effects under the influence of an applied uniform magnetic field were considered by Ali and Alam [2], they observed that the magnetic field exerts a retarding force on the free convection flow, the temperature profile increases as Prandtl and Dufour parameter increases and the concentration increases for increasing values of Shmidt and Soret parameter. Subhakar et al. [3] considered the effect of linear thermal stratification in stable stationary ambient fluid on steady MHD convective flow of a viscous incompressible electrically conducting fluid along a moving, non-isothermal vertical plate in the presence of mass transfer, Soret and Dufour effects and heat generation or absorption, they have found that the local skin-friction coefficient, local heat and mass transfer rates at the plate increase with an increase in the buoyancy forces or Soret or Dufour numbers or heat generation/absorption parameter. The local skin-friction coefficient and the local Nusselt number decreases with increasing the Dufour number (decreasing Soret number) while the local Sherwood number increases. Sharma et al. [4] analyzed The Soret and Dufour effects on unsteady MHD mixed convection flow past an infinite radiative vertical porous plate embedded in a porous medium in the presence of chemical reaction, they have obtained that An increase in radiation parameter, the velocity and temperature profiles increase, with increasing heat source parameter, the velocity and temperature profiles are increasing. The velocity and concentration profiles decrease as an increase in chemical reaction parameter. For fluids with medium molecular weight (H_{2}, air), Dufour and Soret effects should not be neglected. Bouabid et al. [5] have studied the contributions of thermal, diffusive, friction and magnetic terms on entropy generation are investigated. Hidouri et al. [6] have reported the influence of Soret and Dufour effects on entropy generation dealing with transient doublediffusive convection of a binary gas mixture for the special case of opposing buoyancy forces with equal intensity. It was found that for moderate thermal Grashof number, Soret and Dufour parameters induced a slight increase of entropy generation, but for relatively higher thermal Grashof number, oscillatory behavior of entropy generation was obtained.
Magherbi et al. [7] have determined numerically the entropy generation due to heat transfer, mass transfer and fluid friction in steady state for laminar double-diffusive convection, in an inclined enclosure with heat and mass diffusive walls.
The more observed effect was due to heat transfer and then to mass transfer. They obtained that the magnetic effect is more pronounced than friction one. The magnetic field parameter suppresses the flow in the cavity and this lead to a decrease of entropy generation. At local level, results show that entropy generation lines are localized on lower heated and upper cooled regions of the active walls.
Recently, Malashetty and Gaikwad [8] have analyzed the effect of cross diffusion namely Soret coefficient and Dufour coefficient on the double diffusive convection in an unbounded vertically stratified two component system with compensating horizontal thermal and solute gradients. They obtained the effect of the Soret coefficient, S is more significant compared to the Dufour effect. The fluid particle moves with greater velocity and high heat transfer rate for the high value of Dufour coefficient and the velocity is less and high mass transfer rate for the presence of Soret coefficient when the left active concentration wall C_{c} and right wall C_{h}. The opposite behavior was observed for the inversed configuration.
Ahammad and Mollah [9] have dealt numerically with MHD free convection flow and mass transfer over a stretching sheet including the Dufour and Soret effects with a magnetic field, the have concluded that for the combined effects of Dufour inreasing and Soret decreasing, the temperature profiles increase while velocity and concentration profiles decrease. The local skin-friction coefficient and the local Nusselt number decreases with increasing the Dufour number (decreasing Soret number) while the local Sherwood number increases. Despite all these studies, the analysis of thermodynamics second law with cross thermodiffusion Dufour and Soret effects at MHD mixed convection, for an electrically conducting fluid have received few attentions. Hence, the aim of this study is to investigate the entropy generation on steady magneto double-diffusive convection in a square confined cavity filled with an incompressible Newtonian binary gas mixture fluid under Dufour’s effect in the presence of magnetic field. It is assumed that the enclosure is heated and concentrated at their vertical walls in the presence of a transverse magnetic field.
Mathematical formulation
Double-diffusive natural convective flow inside a vertical square cavity of length, L, filled with a binary perfect gas mixture composed of air and a pollutant species in the presence of an oriented magnetic field is considered as shown in Figure 1. Different temperature and concentrations are imposed between the left (θ_{h}, c_{h}) and right vertical walls (θ_{c}, c_{c}), where the θ_{h}>θ_{c} and c_{h}>c_{c}. Adiabatic and impermeable boundary conditions are imposed on the remaining walls. The gravity acts in the downward direction. It is also assumed that the enclosure is permeated by a uniform magnetic field, B=B_{x}e_{x}+B_{y}e_{y} of constant magnitude B_{0}, where B_{x} and B_{y} are space-independent, ex and ey are unit vectors in Cartesian coordinate system. U=ue_{x}+ve_{y} is the field velocity. The orientation of the magnetic field forms an angle ø with horizontal axis. The electric current, J, and the electromagnetic force, F, are defined by J=σe (V×B) and F=σ_{e} (V × B) × B, respectively. The solution of the governing conservation equations of mass, momentum, energy and the chemical species describing the phenomenon is obtained by considering the following assumptions:
- The flow in the cavity is laminar and two-dimensional.
- The Soret effect is neglected.
The binary fluid is assumed to be Newtonian, incompressible and to satisfy the Boussinesq approximation.
All physical properties of the fluid are assumed to be constant, except its density variations upon temperature and concentration which is described by the following state equation:
ρ(θ,c)=ρ_{0}[1−β_{T}(θ−θ_{0})−β_{c}(c−c_{0})] → (1)
ρ_{0}, θ_{0}, and c_{0} are the mass density, the temperature and the concentration at reference conditions, respectively. and are the thermal and the solute expansion coefficients, respectively.
Under these above assumptions, the governing equations for the continuity, momentum, and energy and species diffusion in laminar incompressible flow can be written as follows:
The appropriate initial and boundary conditions are:
For the whole space, at: t=0.
u=v=0, P=0, (θ=θ_{0}) and (c=c_{0}) for 0 ≤ x ≤ L and 0 ≤ y ≤ L → (8.1)
At t}0:
u=v=0, (θ=θ_{h}) and (c=c_{h}) for x=0 and 0 ≤ y ≤ L → (8.3)
u=v=0, (θ=θ_{c}) and (c=c_{c}) for x=L and 0 ≤ y ≤ L → (8.4)
The set of the dimensionless governing equations using dimensionless variables defined below is, given by:
divUd=0 → (9)
In this study, the height of the enclosure L is taken as a reference length for the dimensionless variables. Thus, the dimensionless above equations are established using the following dimensionless variables:
The initial and boundary conditions in dimensionless form that solve the stated Equations system (9-13) are:
For the whole domain, at τ=0:
U=V=0, P=0 (T=0.5−X) and (C=0.5−X) for 0 ≤ X ≤ 1 and 0 ≤ Y≤ 1 → (15.1)
As τ>0
Conditions | Velocity | Temperature | Concentration | |
Bottom and top wall | ||||
Hot wall | X=0 0 ≤ Y ≤ L | U=V=0 | T=0.5 | C=0.5 → (15.3) |
Cold wall | X=1 0 ≤ Y ≤ L | U=V=0 | T=-0.5 | C=-0.5 → (15.3) |
The average heat and mass transfer along the heated wall are given in dimensionless form by Nusselt and Sherwood numbers, respectively. They are defined as follow:
Second law formulation
Entropy generation of a system is the sum of Irreversibilities (degraded energy) that are produced due to different phenomena inside the system. The volumetric local entropy generation in double-diffusive convection through a square medium with the effect of magnetic field can be calculated by the following equation [10]. The first term of the right hand side of Equation (18) is irreversibility due to heat transfer, the second is due to mass transfer, the third to fluid friction and the fourth is relative to the magnetic field.
Total dimensionless entropy generation is obtained by a numerical integration of the dimensionless local entropy generation over the entire volume of the cavity, Ω. It is given by:
Numerical method
The dimensionless governing Equations (9) to (13) subjected to the boundary conditions (15) are integrated numerically using the control volume finite element method CVFEM [11,12]. The above equations are resolved by applying the SIMPLE algorithm of [11]. The SIMPLER algorithm and the SIMPLEC approximation of Van Doormal and Raithby are used in conjunction with an Alternating Direction Implicit (ADI) scheme for performing the time evolution. For more details, the used numerical code written in FORTRAN language and discussions about these functions, reader can see [13-16].
The power law scheme proposed by Patankar [11] is used for the convection terms formulation of the energy and mass transfer equations.
The resulting algebraic equations were solved by line-by-line using the Tri-Diagonal Matrix Algorithm iteration. We assumed linear and exponential variations when the dependent variable χ is calculated in the diffusive and the convective terms of the conservation equations, respectively. χ is the general dependent variable which can stands for u, v, T, or C.
The iteration process is terminated under the following condition.
The local volumetric entropy generation is calculated after calculating the temperature, concentration and stream function in each control volume in the mesh. Finally the entropy generation number for the cavity volume is calculated using Equation (9) via the numerical integration.
Results shown in the following section have been calculated by taking into account the initial and the boundary conditions given above. Table 1 depicts the used grids (i.e., sizes), the dimensionless times as well as the time steps that are found sufficient to reach the steady state situation and the imposed convergence criteria. Error percentage is given by:
Grid size | Nu | Er (%) | Ra=Gr × Pr |
---|---|---|---|
21×21 | 1.1024 | -- | 10^{3} |
31×31 | 1.0992 | 0.29 | |
41×41 | 1.0979 | 0.408 | |
31×31 | 2.3201 | -- | 10^{4} |
41×41 | 2.2956 | 1.055 | |
51×51 | 2.2832 | 1.59 | |
41×41 | 4.7097 | -- | 10^{5} |
51×51 | 4.6641 | 0.968 | |
61×61 | 4.6593 | 1.07 |
Table 1: Grid size independence study (Pr=0.71, N=Ha=0).
Er(%)=(│Nu(x+10, x+10) −Nu (x,x)│/ Nu(x,x)) × 100
It is found that 51 × 51 grid size is an appropriate size for the most calculations in the present problem, x represents the Grid size. Comparison of the calculated average Nusselt number using the present numerical method for different thermal Grashof numbers with literature is presented by Davis and Nithyadevi [16,17] (see Table 2). For various Prandtl number values, Table 3 illustrates a comparison of average Nusselt number results with those published by Pesso et al. [18] at fixed thermal Grashof number, Gr_{T}=10^{5}. Isotherms and streamlines are plotted in Figure 2 for different thermal Grashof number and compared to results obtained by Abbassi et al. [13] For thermosolutal natural convection in a square cavity, the present simulation is validated against numerical results of Sezai and Teamah [18,19] for different Lewis numbers and Buoyancy ratios. Figures 3 and 4 illustrate a validation by; avearage Nusselt and Sherwood numbers versus Lewis number and versus Buoyancy ratio variations and by isotherms and isoconcentrations contours at fixed Rayleigh number, Ra=10^{5}, Prandtl number, Pr=10 and Buoyancy ratio, N=-0.2 [20-23].
Ra | 10^{3} | 10^{4} | 10^{5} | |
Present study | 1.099 | 2.295 | 4.664 | |
Davis | 1.118 | 2.243 | 4.519 | |
Nithyadevi et al. | 1.123 | 2.304 | 4.899 |
Table 2: Average Nusselt number for different values of Rayleigh number (Pr=0.71, N=Ha=0).
Pr | 0.071 | 0.71 | 1 | 7.1 |
Present study | 3.948 | 4.664 | 4.745 | 4.863 |
Pesso et al. | 3.813 | 4.521 | 4.685 | 4.732 |
Relative error (Er(%)) | 3.54 | 3.16 | 1.28 | 2.77 |
Table 4: Average Nusselt number for different values of Prandtl numbers (Gr_{T}=10^{5}, N=Ha=0).
Relatively good agreement is obtained. We are, therefore, confident that our results are accurate.
The present work focused the study of the fluid flow structure and entropy generation in double-diffusive natural convection through a square cavity submitted to an oriented magnetic field with Dufour effect and saturated by a binary perfect gas mixture of air and a pollutant species characterized by Pr=0.71 and Le=1.2.
Numerical computations of double-diffusive natural convection in a two-dimensional vertical square cavity in the presence of Magnetic and Dufour effects have been carried out for different values of N, Ha, Du, and for fixed values of Pr=0.71, which corresponds to air at 20°C and 1 atmospheric pressure, Ø=0°, Le=1.2 and Ra_{T}=10^{5}. With the abovementioned flow parameters, the results are displayed in Figures 5-10, for the velocity, temperature and concentration profiles and for entropy generation rate behaviour. The rate of heat and mass transfer in the enclosure is measured in terms of the average Nusselt number and average Sherwood number.
The Dufour effect consists on concentration grdients that drive heat currents or change the temperature field on convection in binary heated fluid layers. The Dufour effect changes the temperarure field equation “diagonally” via the term (−Δ^{2}T) which enhances heat diffusion and “off-diagonally” via the term (−Δ^{2}C) which reflects concentrationinduced temperature variations.
The effects of Dufour parameter on the velocity field for N=-1, Ha=0 and Ra_{T}=10^{5} are shown in Figure 5. It exhibits a comparison of center line velocity profiles along x and y directions at the horizontal and vertical center lines in the cavity for different values of Du. the velocity is substantially zero in the absence of the Dufour effect (Du=0) indicating a state of stable equilibrium. Increasing Du leads to an increase in the fluid velocity field by means of an increase in the temperature boundary layer thickness due to augmentation of the fluid temperature which is proportional to Dufour effect, as a result the birth of a convective flow following the disappearance of the diffusive flow. The increase in velocity with Dufour parameter indicates that the buoyant convection can be enhanced by an increase in Du. these results may be explained by considering the temperature trend as depicted in Figure 6a. The influence of Dufour effect on the concentration field is illustrated in Figure 6b.
Figure 7 displays the variations of entropy generation rate (S) on the buoyancy ratio (N) at two Hartmann number values (Ha=0 and Ha=25) for different Dufour parameter (Du). It can be noticed that the amplitude of the entropy generation increases (respectively decreases) with the Dufour parameter in the case of cooperating buoyancy forces (respectively opposing case) whereas it decreases with Hartmann number increasing. The entropy generation takes a minimum value obtained at N=−1. Dufour effect causes a significant increase in irreversibility due to friction. Therefore the movement of the fluid within the cavity is high, resulting from an increase in the fluid velocity field caused by the buoyancy forces. Increasing Hartmann number leads to reduce the flow circulation strength in the cavity and this leads to decrease the rate of heat and mass transfer as well as the rate of entropy generation. Enhancing the magnetic effect leads to a decrease in the velocity field as can be seen in Figure 8, the centreline velocities are signicantly reduced due to a development of MHD Hartmann boundary layer. The magnetic field presents a damping effect on the velocity by creating a drag force that tends to suppress the fluid flow.
The effects of the Dufour parameter (Du) on the average Nusselt and Sherwood numbers and entropy generation rate at fixed buoyancy ratio (N=-1) are presented in Figure 9 which shows that for absence of magnetic field (Ha=0), when the Dufour parameter increases, they increase. Additionally, for moderate magnetic effect (Ha=25), the average Nusselt and Sherwood numbers slightly decrease until Du=0.2, then they increase as Dufour parameter increases. For high values of Hartmann number (Ha ≥ 100), the average Nusselt and Sherwood numbers decrease until Du=0.2, then they are no influenced by Dufour effect and remain constant approching unity. The entropy generation ratio shows a monotonic behaviour with Dufour increase for low magnetic effect, exponetial trend for moderate magnetic field and no significant effect at high values of Hartmann number.
Figure 10 shows the variations of average Nusselt (Nu) and Sherwood (Sh) numbers and entropy generation rate (S) with buoyancy ratio (N) for different Hartmann number (Ha) at fixed value of Dufour parameter (Du=0.2). Increasing the magnetic field leads to increase the Lorentz forces which tends to suppress the fluid flow leading to a decrease of Nu and Sh rapidly towards unity, therefore the entropy generation decreases but it increases with the absolute value of buoyancy ratio, N. For all the values of Ha, minimum values of Nu, Sh and S are obtained at N=-1 which means that the buoyancy forces generated due to temperature and concentration gradients are acting equally in opposite directions. This case corresponds to stagnant fluid where the heat and mass transfers are in pure diffusion processes.
Figure 11 show the effect of Hartmann number at fixed Dufour parameter, RaT=105 and N=-1 on isotherms, streamlines and isentropic lines. As magnetic effect increases, it is observed that the vortex strength increases, the streamline behaviour is very pronounced at high Hartmann numbers (Ha=100), it can be seen that the streamline is elongated and two secondary vortices appear inside the cavity. These influences on the vortex strength are proportional to the Dufour parameter increasing. The flow is driven by the buoyancy force due to the effect of both temperature and concentration variations and draged by Lorentz force which has a retarding effect to the fluid flow.
For low Hartmann number situations, most of the buoyancy driver is due to the effect of temperature variations. The effects of Hartmann number are clearly seen on the isotherms at moderate and high Dufour parameter, it could be seen that for low Hartmann number, the isotherms are nearly parallel near the bottom left wall and near the upper right corner as a result of the fluid rising due to the influence of the buoyancy force. Further increase in Ha eventually leads to a linear temperature variation between the hot and cold regions, almost like heat conduction in a solid. The isotherm lines become parallel corresponding to pure conduction regime. An increase in Dufour parameter enhance entropy generation in the enlosure where its maximum is concentrated near the hot and cold walls but with increasing magnetic effect, isentropic lines prolfie tends to cover the whole domain. This is due to the anisotropic nature of MHD flow.
The influence of magnetic and Dufour effects on entropy generation inside a square cavity filled with a binary perfect gas mixture at doublediffusive convection is accounted. Various Hartmann numbers, Buoyancy ratios and Dufour parameters have been considered for the flow, temperature and concentration fields as well as the entropy generation rate. The results of the numerical analysis lead to the following conclusions:
• Increasing Du leads to an increase in the fluid velocity field by means of an increase in the temperature boundary layer thickness,
• The amplitude of the entropy generation increases (respectively decreases) with the Dufour parameter in the case of cooperating buoyancy forces (respectively opposing case),
• The magnetic field presents a damping effect on the velocity by creating a drag force that tends to suppress the fluid flow,
• The entropy generation ratio shows a monotonic behaviour with Dufour increase for low magnetic effect, exponetial trend for moderate magnetic field and no significant effect at high values of Hartmann number,
• For all the values of Ha, minimum values of Nu, Sh and S are obtained at N=−1,
• The structure of the fluid streamlines, isotherms and isoentropic lines within the enclosure is found to significantly depend upon the Ha and Du.
The authors are grateful to a referee and an editor for great suggestions and help in polishing the manuscript.