Entropy Generation in MHD Porous Channel Flow Under Constant Pressure Gradient

Transcription

1 Applied Mathematics and Physics, 013, Vol. 1, No. 3, Available online at Science and Education Publishing DOI: /amp Entropy Generation in MHD Porous Channel Flow Under Constant Pressure Gradient Sanatan Das 1,*, abindra Nath Jana 1 Department of Mathematics, University of Gour Banga, Malda, India Department of Applied Mathematics, Vidyasagar University, Midnapore, India *Corresponding author: jana61171@yahoo.co.in eceived September 19, 013; evised October 04, 013; Accepted October 09, 013 Abstract Effects of magnetic field and suction/injection on the entropy generation in a flow of a viscous incompressible electrically conducting fluid between two infinite horizontal parallel porous plates under a constant pressure gradient have been studied. An exact solution of governing equation has been obtained in closed form. The entropy generation number and the Bejan number are also obtained. The influences of each of the governing parameters on velocity, temperature, entropy generation and Bejan number are discussed with the aid of graphs. It is observed that the magnetic field tends to increase the entropy generation within the channel. Keywords: MHD flow, porous channel, eynolds number, entropy generation number, Bejan number Cite This Article: Sanatan Das, and abindra Nath Jana, Entropy Generation in MHD Porous Channel Flow Under Constant Pressure Gradient. Applied Mathematics and Physics 1, no. 3 (013): doi: /amp Introduction Entropy generation minimization studies are vital for ensuring optimal thermal systems in contemporary idustrial and technological fields like geothermal systems, electronic cooling, heat exchangers to name a few. All thermal systems confront with entropy generation. Entropy generation is squarely associated with thermodynamic irreversibility. Analysis of the flow of an electrically conducting fluid in a porous channel in the presence of a transverse magnetic field is of special technical significance because of its widespread engineering and industrial applications such as in geothermal reservoirs, nuclear reactor cooling, MHD marine propulsion, electronic packages, micro electronic devices, thermal insulation, and petroleum reservoirs. Some other quite promising applications are in the field of metallurgy such as MHD stirring of molten metal and magnetic-levitation casting. This type of problem also arises in electronic packages and microelectronic devices during their operation. The contemporary trend in the field of heat transfer and thermal design is the second law analysis and its design-related concept of entropy generation minimization. The foundation of knowledge of entropy production goes back to Clausius and Kelvins studies on the irreversible aspects of the second law of thermodynamics. Since then the theories based on these foundations have rapidly developed. However, the entropy production resulting from combined effects of velocity and temperature gradients has remained untreated by classical thermodynamics, which motivates many researchers to conduct analyses of fundamental and applied engineering problems based on second law analysis. Entropy generation is associated with thermodynamic irreversibility, which is common in all types of heat transfer processes. Moreover, in thermodynamical analysis of flow and heat transfer processes, one thing of core interest is to improve the thermal systems to avoid the energy losses and fully utilize the energy resources. The magnetohydrodynamic channel flow with heat transfer has attracted the attention of many researchers due to its numerous engineering and industrial applications. Such flows finds applications in thermofluid transport modeling in magnetic geosystems, meteorology, turbo machinery, solidification process in metallurgy and in some astrophysical problems. One of the methods used for predicting the performance of the engineering processes is the second law analysis. The second law of thermodynamics is applied to investigate the irreversibilities in terms of the entropy generation rate. Since entropy generation is the measure of the destruction of the available work of the system, the determination of the factors responsible for the entropy generation is also important in upgrading the system performances. The method is introduced by Bejan [1,]. The entropy generation is encountered in many of energyrelated applications, such as solar power collectors, geothermal energy systems and the cooling of modern electronic systems. Efficient utilization of energy is the primary objective in the design of any thermodynamic system. This can be achieved by minimizing entropy generation in processes. The theoretical method of entropy generation has been used in the specialized literature to treat external and internal irreversibilities. The irreversibility phenomena, which are expressed by entropy generation in a given system, are related to heat and mass

2 Applied Mathematics and Physics 79 transfers, viscous dissipation, magnetic field etc. Several researchers have discussed the irreversibility in a system under various flow configurations [3-14]. They showed that the pertinent flow parameters might be chosen in order to minimize entropy generation inside the system. Salas et al. [15] analytically showed a way of applying entropy generation analysis for modelling and optimization of magnetohydrodynamic induction devices. Salas et al. [15] restricted their analysis to only Hartmann model flow in a channel. Thermodynamics analysis of mixed convection in a channel with transverse hydromagnetic effect has been investigated by Mahmud et al. [16]. Flow, thermal and entropy generation characteristic inside a porous channel with viscous dissipation have been investigated by Mahmud [17]. The heat transfer and entropy generation during compressible fluid flow in a channel partially filled with porous medium have been analyzed by Chauhan and Kumar [18]. Tasnim et al. [19] have studied the entropy generation in a porous channel with hydromagnetic effects. Eegunjobi and Makinde [0] have studied the combined effect of buoyancy force and Navier slip on entropy generation in a vertical porous channel. The second law analysis of laminar flow in a channel filled with saturated porous media has been studied by Makinde and Osalusi [1]. Makinde and Maserumule [] has presented the thermal criticality and entropy analysis for a variable viscosity Couette flow. Makinde and Osalusi [3] have investigated the entropy generation in a liquid film falling along an incline porous heated plate. Cimpean and Pop [4,5] have presented the parametric analysis of entropy generation in a channel. The effect of an external oriented magnetic field on entropy generation in natural convection has been investigated by Jery et al. [6]. Dwivedi et al. [7] have made an analysis on the incompressible viscous laminar flow through a channel filled with porous media. ecently, Makinde1 and Chinyoka [8] have discussed the numerical investigation of buoyancy effects on hydromagnetic unsteady flow through a porous channel with suction/injection. In the present paper, we study the entropy generation in MHD flow in a porous channel under the constant pressure gradient. Both the channel walls are maintained at constant but different temperatures. A parametric study is carried out to see how the pertinent parameters of the problem affect the flow field, temperature field and the entropy generation.. Mathematical Formulation and its Solution Consider the viscous incompressible electrically conducting fluid bounded by two infinite horizontal parallel porous plates separated by a distance h. Choose a cartesian co-ordinate system with x -axis along the centerline of the channel, the y -axis is perpendicular to the planes of the plates y = ± h. A uniform transverse magnetic field B 0 is applied perpendicular to the channel plates. The plates are maintained at constant temperatures T 0 and T h respectively. Since the magnetic eynolds number is very small for most fluid used in industrial applications, we assume that the induced magnetic field is negligible. Since the plates are infinitely long, all physical variables, except pressure, depend on y only. The equation of continuity then gives v = v0 everywhere in the fluid where v 0 is the suction velocity at the plates. Figure 1. Geometry of the problem The Ohms law neglecting Hall currents, the ion-slip and thermoelectric effects as well as the electron pressure gradient for a conducting fluid is J = σ[ E+ q B], (1) where E is the the electric field vector, J the current density vector, q the velocity vector, σ the electrical conductivity of the fluid. The equations of motion along x -direction is du 1 p d u B0 v0 = + ν + J z dy ρ x dy ρ, () where u is the fluid velocity in the x -direction, ρ is the fluid density, ν the kinematic viscosity and p is the fluid pressure. The energy equation is dt k d T µ du J v0 = + λ, dy ρ c + (3) p dy ρ c p dy σ where µ is the coefficient of viscosity, k the thermal conductivity, c p the specific heat at constant pressure, the index λ in the equation (3) is set equal to 0 for excluding Jule dissipation and 1 for including Jule dissipation. The boundary conditions are u = 0, T = T at y = h, 0 u = 0, T = Th at y = h, (4) where T is the fluid temperature, T h the temperature at upper plate and T 0 the temperature at the lower plate. de For steady motion E = 0 which yields x 0 dy = de and z 0 dy = and hence E x = constant and E z = constant. On taking Ex = 0 = Ez, we have

3 80 Applied Mathematics and Physics Jx = 0, Jz = σ ub0. (5) On the use of (5), equations () and (3) become du 1 p d u σ B0 v0 = + ν u, (6) dy ρ x dy ρ dt k d T µ du v0 = + B + λσ 0u. dy ρcp y ρcp dy Introducing the non-dimensional variables y hu T T0 η =, u1 =, θ =, h ν Th T0 equations (6) and (7) become where (7) (8) d u1 du1 M u 1 P, dη + dη = (9) dθ d θ du 1 Pe = + Br λ M u + 1, dη dη dη µν Br = kt ( h T0 ) h (10) σ B h M = is the magnetic parameter, ρν the Peclet number, 3 h p P = x the Brinkmann number, v0 hρcp Pe = k v0 h = the eynolds number and ν the non-dimensional pressure gradient. ρν The boundary conditions for u 1 ( η ) and θη ( ) are u1 = 0, θ = 0 atη = 1 and u1 = 0, θ = 1 at η = 1. (11) The solution of the equation (9) subject to the boundary conditions (11) is P η u ( ) 1( η) 1 a 1cosh rη asinh rη e = +, (1) M where 1 cosh sinh r = + M, a 1 = a nd a =.(13) 4 cosh r sinh r The solution given by the equation (1) exists for both <0(corresponding to v 0 <0 for the blowing at the plates) and >0(corresponding to v 0 >0 for the suction at the plates). On the use of (1), the equation (10) becomes d θ dθ + Pe dη d η PBr {( a1 a 4 ) λ M ( a3 a4) } M { a1 a λm ( a3 a4) } cosh r PBr η = 4( a ) 4 1cosh rη+ asinh rη e M η + ( aa 1 + λm aa 3 4)sinh rη e. η (14) Solution of the equation (14) subject to boundary condition (11) is given by Peη ( ) = c1 + ce θη where PBr η η e ( X 4 3cosh rη + X4sinh rη) M Pe η + e ( X7 cosh rη + X8 sinh rη) + X9, 1 1 a3 = a1 ra, a4 = a ra1, (15) a5= a1+ a+ λm ( a+ a3), a6= ( a1a+ λm a3a4), X1 = r r + Pe, X = r + r + Pe, X3 = a1 + + a, X1 X X1 X X4 = a1 + a +, X1 X X1 X ( )( ) ( )( ) X5 = r r + Pe, X = r+ r+ Pe, X7 = a5 + + a6, X5 X6 X5 X X8 = a5 + a6 +, X5 X6 X5 X6 a1 a + λ M ( a3 a4) X9 =, ( Pe) Pr Br A= e ( X 4 3cosh r+ X4sinh r) M Pe + e ( X7cosh r+ X8sinh r) + X 9, Pr Br B = e ( X 4 3cosh r+ X4sinh r) M Pe + e ( X7cosh r+ X8sinh r) + X 9 Pe Pe (1 + B) e Ae A B c1 =, c =. sinh Pe sinh Pe 3. esults and Discussion (16) To study the effects of magnetic field and eynolds number on the velocity field we have presented the nondimensional velocity u 1 against η in Figure and Figure 3 for several values of magnetic parameter M and eynolds number when P = 1. It is seen from Figure that the fluid velocity u 1 decreases with an increase in

4 Applied Mathematics and Physics 81 magnetic parameter M. This can be attributed to the presence of Lorentz force acting as a resistance to the flow as expected. Figure 3 shows that the fluid velocity u 1 increases near the lower plate and decreases near the upper plate with an increase in eynolds number e. We have plotted the temperature distribution θ against η in Figure 4, Figure 5 and Figure 6 for several values of magnetic parameter M, Brinkmann Br and Peclet number Pe. It is seen from Figure 4, Figure 5 and Figure 6 that the fluid temperature θ increases with an increase in either magnetic parameter M or Brinkmann Br or Peclet number Pe. As M increases due to increasing magnetic field intensity, the fluid temperature increases within the channel. So, an increase in the fluid temperature is due to the presence of Ohmic heating (or Lorentz heating) which serves as additional heat source to the flow system. The terms linked to the Brinkman number act as strong heat sources in the energy equation. Increases in the Brinkman number hence significantly also increase the fluid temperature shown in Figure 5. Figure. Velocity profiles for different M when = 0.5 Figure 3. Velocity profiles for different when M = 5

5 8 Applied Mathematics and Physics Figure 4. Temperature profiles for different M when = 0.5 and Pe = Figure 5. Temperature profiles for different Br when M = 5 and Pe =

6 Applied Mathematics and Physics 83 Figure 6. Temperature profiles for different Pe when M = 5 and = 0.5 The rate of heat transfer at the plates η = 1 and η = 1 can be obtained from (15) as θ ( 1) = c Pee θ Pe ( X 3 rx 4)cosh r + e PBr Pe ( X 4 rx3)sinh r (17) 4 M ( X7 rx8)cosh r e ( X8 rx7)sinh r X 9 Pe (1) = c Pee ( X 3 rx 4) cosh r + e PBr Pe ( X 4 rx 3) sinh r (18), 4 M ( X 7 rx8) cosh r e + ( X8 rx 7) sinh r X 9 where X 3, X 4, X 7, X 8 and c are given by (16). The numerical values of the rate of heat transfers θ ( 1) and θ (1) are entered in the Table 1 and Table for several values of M,, Br and Pe. It is seen from the Table that the rate of heat transfer θ ( 1) at the plate η = 1 increases whereas the rate of heat transfer θ (1) at the plate η = 1 decreases with an increase in either magnetic parameter M or eynolds number. Table shows that the rate of heat transfer θ ( 1) increases with an increase in either Peclet number Pe or Brinkman number Br. The rate of heat transfer θ (1) increases with an increase in Brinkman number Br while it decreases with an increase in Peclet number Pe. On the other hand, θ (1) < 0 means the heat transfer take places from fluid to the upper plate, because when there is significant heat generation in the fluid due to viscous and Ohmic dissipations, the temperature of the fluid exceeds the plate temperature which causes heat flow from fluid to the plate. Table 1. ate of heat transfer at the plates η = 1 and η = 1 when Pe = 0. and Br = 1.4 θ ( 1) θ (1) \ M

7 84 Applied Mathematics and Physics Table. ate of heat transfer at the plates η = 0 and η = 1 with = 0.5 and M = 5 θ ( 1) θ (1) Br \ Pe N where Φ= is the irreversibility ratio. Heat transfer N1 4. Entropy Generation Second law analysis in terms of entropy generation rate is a useful tool for predicting the performance of the engineering processes by investigating the irreversibility arising during the processes. According to Wood [9], the local entropy generation rate is defined as k dt µ du σ B E 0 G = u. + + λ (19) T dy T 0 0 dy T0 The first term in equation (0) is the irreversibility due to heat transfer, the second term is the entropy generation due to viscous dissipation and third term is due to magnetic field. The dimensionless entropy generation number may be defined by the following relationship: N S 0 T h E = kt ( T ) h G 0. (0) In terms of the dimensionless velocity and temperature, the entropy generation number becomes dθ Br du N 1 S = + + dη Ω dη λ M u1, (1) µ e u0 where Br = k ( Th T0 ) is the Brinkmann number and Tw T0 Ω= T0 is the non-dimensional temperature difference. The entropy generation number N S can be written as a summation of the entropy generation due to heat transfer denoted by N 1 and the entropy generation due to fluid friction with magnetic field denoted by N given as dθ Br du 1 N1 =, N = + λ M u1. dη Ω dη () The entropy in a system is associated with the presence of irreversibility. We have to notice that the contribution of the heat transfer entropy generation to the overall entropy generation rate is needed in many engineering applications. The Bejan number Be is an alternative irreversibility distribution parameter and it represents the ratio between the heat transfer irreversibility N 1 and the total irreversibility due to heat transfer and fluid friction with magnetic field N S. It is defined by N1 1 Be = =, (3) NS 1+Φ dominates for 0 Φ <1 and fluid friction with magnetic effects dominates when Φ >1. The contribution of both heat transfer and fluid friction to entropy generation are equal when Φ =1. The Bejan number Be takes the values between 0 and 1 (see Cimpean et al. [30]). The value of Be = 1 is the limit at which the heat transfer irreversibility dominates, Be = 0 is the opposite limit at which the irreversibility is dominated by the combined effects of fluid friction and magnetic field and Be = 0.5 is the case in which the heat transfer and fluid friction with magnetic field entropy production rates are equal. Further, the behavior of the Bejan number is studied for the optimum values of the parameters at which the entropy generation takes its minimum. The influences of the different governing parameters on entropy generation within the channel are presented in Figures It is seen from Figure 7 that the entropy generation number N S decreases with an increase in magnetic parameter M. The effect of magnetic parameter M on the entropy generation number is displayed in Figure 7. This figure shows that the entropy generation is slightly increases with magnetic parameter M. The magnetic parameter M is not too much dominating on entropy generation. A large variation of M causes a small variation in the rate of entropy generation. Figure 8 and Figure 9 show that the entropy generation number N S increases near lower plate with an increase in either eynolds number or Peclet number Pe. As expected, since the fluid velocity and temperature increase near the lower plate. An increase in Peclet number Pe has no effect on the entropy generation in the central region of the channel as well as near the upper plate. Figure 10 reveals that the entropy generation number N S increases with an increase in group parameter BrΩ. It is seen from Figure 11 that the Bejan number Be decreases near the central region of channel with an increase in magnetic parameter M. Figure 1 shows that the Bejan number Be decreases near the central region of channel with an increase in eynolds number. Figure 13 and Figure 14 shows that the Bejan number Be decreases near the central region of channel with an increase in either Peclet number Pe or group parameter BrΩ. The Bejan number Be is unaffected at the central point of the channel width for Peclet number Pe and group parameter BrΩ. The group parameter is an important dimensionless number for irreversibility analysis. It determines the relative importance of viscous effects to that of temperature gradient entropy generation.

8 Applied Mathematics and Physics 85 Figure 7. N S for different M when Pe =, BrΩ = 1 and = 0.5 Figure 8. N S for different when Pe = and M = 5, BrΩ = 1

9 86 Applied Mathematics and Physics Figure 9. N S for different Pe when M = 5, Pe =, BrΩ = 1 and = 0.5 Figure 10. N S for different BrΩ when M = 5, Pe =, BrΩ = 1 and = 0.5

10 Applied Mathematics and Physics 87 Figure 11. Bejan number for different M when Pe =, BrΩ = 1 and = 0.5 Figure 1. Bejan number for different when M = 5, Pe = and BrΩ = 1

11 88 Applied Mathematics and Physics Figure 13. Bejan number for different Pe when M = 5, BrΩ = 1 and = 0.5 Figure 14. Bejan number for different BrΩ, M = 5, Pe = and = 0.5

12 Applied Mathematics and Physics Conclusion We investigate the entropy generation in an MHD flow between two infinite parallel porous plates. The analytical results obtained for the velocity and temperature profiles are used in order to obtain the entropy generation production. The non-dimensional entropy generation number and the Bejan number are calculated for the problem involved. It is found that, the entropy generation decreases with an increase in magnetic parameter. The entropy generation increases with an increase in either eynolds number or Peclet number or group parameter. The rate of heat transfer at the lower plate increases whereas the rate of heat transfer at the upper plate decreases with an increase in either magnetic parameter or eynolds number. It is important to note that the fluid velocity, the fluid temperature as well as entropy generation are significantly influenced by Jule dissipation. Acknowledgement The authors would like to express thanks to the anonymous referees for their valuable suggestions. eferences [1] Bejan, A., Second law analysis in heat transfer. Energy Int. J. 5: 71-73, [] Bejan, A., Entropy Generation Through Heat and Fluid Flow, Wiley, Canada, [3] Bejan, A. Second-law analysis in heat transfer and thermal design. Adv. Heat Transf. 15: 1-58, 198. [4] Bejan, A., Entropy Generation Minimization, CC Press: New York, USA, [5] Bejan, A., A study of entropy generation in fundamental convective heat transfer. J. Heat Transfer. 101: , [6] Bejan, A., Tsatsaronis, G. and Moran, M. Thermal Design and Optimization, Wiley: New York, USA, [7] Arpaci, V.S. and Selamet, A. Entropy production in flames. Combust. Flame.73: 54-59, [8] Arpaci, V.S. and Selamet, A. Entropy production in boundary layers. J. Thermophys. Heat Tr. 4: , [9] Arpaci, V.S. adiative entropy production-heat lost to entropy. Adv. Heat Transfer. 1: 39-76, [10] Arpaci, V.S. Thermal deformation: From thermodynamics to heat transfer. J. Heat Transfer. 13: 81-86, 001. [11] Arpaci, V.S. and Esmaeeli, A., adiative deformation. J. Appl. Phys. 87: , 000. [1] Magherbi, M., Abbassi, H. and Ben Brahim, A. Entropy generation at the onset of natural convection. Int. J. Heat Mass Transfer 46: , 003. [13] Magherbi, M., Abbassi, H., Hidouri N. and Ben Brahim, A. Second law analysis in convective heat and mass transfer. Entropy. 8: [14] Abbassi, H., Magherbi, M. and Ben Brahim, A. Entropy generation in Poiseuille-Benard channel flow. Int. J. Therm. Sci. 4: , 003. [15] Salas, S., Cuevas, S. and Haro, M.L., Entropy generation analysis of magnetohydrodynamic induction devices. J. Phys. D: Appl. Phys. 3: , [16] Mahmud, S., Tasnim, S.H. and Mamun, H. A. A., Thermodynamics analysis of mixed convection in a channel with transverse hydromagnetic effect. Int. J Therm. Sci. 4: , 003. [17] Mahmud, S. and Fraser,.A., Flow, thermal and entropy generation characteristic inside a porous channel with viscous dissipation. Int. J. Therm. Sci. 44: 1-3, 005. [18] Chauhan, D.S. and Kumar, V., Heat transfer and entropy generation during compressible fluid flow in a channel partially filled with porous medium. Int. J. Energ. Tech. 3: 1-10, 011. [19] Tasnim, S.H., Mahmud, S. and Mamun, M.A.H., Entropy generation in a porous channel with hydromagnetic effects. Exergy. : , 00. [0] Eegunjobi, A.S. and Makinde, O.D., Combined effect of buoyancy force and navier slip on entropy generation in a vertical porous channel. Entropy. 14, , 01. [1] Makinde, O.D. and Osalusi, E., Second law analysis of laminar flow in a channel filled with saturated porous media. Entropy. 7 (): , 005. [] Makinde, O.D. and Maserumule,.L., Thermal criticality and entropy analysis for a variable viscosity Couette flow. Phys. Scr. 78: 1-6, 008. [3] Makinde, O.D. and Osalusi, E., Entropy generation in a liquid film falling along an incline porous heated plate. Mech. es. Commun. 33: , 006. [4] Cimpean, D. and Pop I., Parametric analysis of entropy generation in a channel filled with a porous medium. ecent esearches in Applied and Computational Mathematics, WSEAS ICACM, 011, [5] Cimpean, D. and Pop I., A study of entropy generation minimization in an inclined channel. 6(): 31-40, 011. [6] Jery, A. E., Hidouri, N., Magherbi, M. and Ben Brahim, A., Effect of an external oriented magnetic field on entropy generation in natural convection. Entropy. 1: , 010. [7] Dwivedi,., Singh, S. P. and Singh, B. B., Analysis of incompressible viscous laminar flow through a channel filled with porous media. Int. J. Stability and Fluid Mechanics. 1(1): , 010. [8] Makinde, O. D. and Chinyoka, T., Numerical investigation of buoyancy effects on hydromagnetic unsteady flow through a porous channel with suction/injection. J. Mechanical Science and Technology. 7 (5): , 013. [9] Woods, L. C., Thermodynamics of Fluid Systems, Oxford University Press, Oxford, UK, [30] Cimpean, D., Lungu, N. and Pop, I., A problem of entropy generation in a channel filled with a porous medium. Creative Math and Inf. 17: , 008.