# Entropy Generation Analysis in a Variable Viscosity MHD Channel Flow with Permeable Walls and Convective Heating

Authors: A. S. Eegunjobi; O. D. Makinde
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

This paper examines the effects of the thermodynamic second law on steady flow of an incompressible variable viscosity electrically conducting fluid in a channel with permeable walls and convective surface boundary conditions. The nonlinear model governing equations are solved numerically using shooting quadrature. Numerical results of the velocity and temperature profiles are utilised to compute the entropy generation number and the Bejan number. The results revealed that entropy generation minimization can be achieved by appropriate combination of the regulated values of thermophysical parameters controlling the flow systems.


---

## Body

## 1. Introduction

Hydromagnetic channel flows have attracted the attention of many researchers due to their numerous engineering and industrial applications. Such flow can be found in magnetohydrodynamic (MHD) generator, geothermal reservoirs, cooling of nuclear reactors, petroleum reservoirs, accelerators, pumps, flow meter, astrophysics, metallurgy, crystal growth, magnetic filtration and separation, jet printers, and microfluidic devices [1]. Several researchers have discussed MHD fluid flow under various physical situations [2–4]. Lehnert [5] presented a theoretical investigation on the behavior of electrical conducting liquid under magnetic field. Makinde and Mhone [6] investigated the combined effect of transverse magnetic field and radiative heat transfer on unsteady flow of a conducting optically thin fluid through a channel filled with saturated porous medium and nonuniform wall temperature. Seth et al. [7] studied unsteady MHD Couette flow of a viscous incompressible electrically conducting fluid between two parallel porous plates in the presence of a transverse magnetic field. Agarwal [8] analyzed the effect of magnetic field on generalized Couette flow. The combined effects of variable viscosity and electrical conductivity on MHD generalized Couette flow and heat transfer were numerically investigated by Makinde and Onyejekwe [9].Meanwhile, most industrial and engineering flow processes and thermal systems are unable to work at optimal level due to entropy production. Therefore, it is imperative to determine the factors that contributed to entropy generation in order to minimize their effects and maximize the flow system efficiency. The analysis of entropy generation minimization in a thermal system was pioneered by Bejan [10]. Thereafter, several researchers have theoretically studied entropy generation in thermal and flow systems under many physical situations [11–13]. Sahin and Ben-Mansour [14] reported a numerical solution of the entropy generation in a circular pipe. Hooman [15] studied the effects of different thermal boundary conditions on entropy generation in a microscale forced convection with velocity slip. The effect of Navier slip on entropy generation in a porous channel with suction/injection was investigated by Eegunjobi and Makinde [16]. Recently, Makinde and Eegunjobi [17] reported a numerical solution for the effects of convective heating on entropy production in a channel with permeable walls.In this present study, the recent work of Makinde and Eegunjobi [17] is extended to include the combined effects of variable viscosity and asymmetric convective boundary conditions on the entropy generation rate in MHD porous channel flow. In the following sections, the model problem is formulated, analyzed, and numerically solved. Pertinent results are presented graphically and discussed quantitatively, with respect to various thermophysical parameters controlling the flow system.


## 2. Mathematical Model

We consider a steady, incompressible flow of an electrically conducting variable viscosity fluid between two fixed permeable parallel infinite plates of widthh. The flow is fully developed and the edge effects are disregarded. A constant magnetic field of strength B0 is imposed transversely in the y-direction. In addition, both the electric field and Hall effect are not present (see Seth et al. [7], Turkyilmazoglu [18]). The applied magnetic field is assumed to be strong enough so that the induced magnetic field due to the fluid motion is weak and can be neglected. It is assumed that the lower permeable plate, where fluid injection occurs, is convectively heated; while at the upper permeable plate both fluid suction and convective heat loss take place as shown in Figure 1.Figure 1
Schematic diagram of the problem.Under these assumptions, the governing equations for the momentum and energy balance in one dimension can be written as follows [7–9, 16, 17]:
(1)Vdudy=-1ρdPdx+1ρddy(μ-(T)dudy)-σB02uρ,(2)VdTdy=αd2Tdy2+μ-(T)ρcP(dudy)2+σB02u2ρcp.
The boundary conditions are
(3)u(0)=0,-kdTdy(0)=γ0(Tf-T(0)),u(h)=0,-kdTdy(h)=γ1(T(h)-T∞),
where (x,y) is the axial and normal coordinates, u is the velocity of the fluid, P is the fluid pressure, V is the uniform suction/injection velocity at the channel walls, γ0 is the heat transfer coefficient at the lower plate, γ1 is the heat transfer coefficient at the upper plate, α is the thermal diffusivity, ρ is the fluid density, σ is the fluid electrical conductivity, k is the thermal conductivity coefficient, cp is the specific heat at constant pressure, Tf is the temperature, of the hot fluid at the lower permeable plate, T is the channel fluid temperature, and T∞ is the ambient temperature above the upper plate. The temperature dependent viscosity μ- can be expressed as [9]
(4)μ-(T)=μ0e-m(T-T∞),
where m is a viscosity variation parameter and μ0 is the fluid dynamic viscosity at the ambient temperature. We introduce the following nondimensional quantities:
(5)η=yh,α=kρcp,w=uV,θ=T-T∞Tf-T∞,X=xh,P-=phμ0V,G=-∂P-∂X,μ=μ-μ0,υ=μ0ρ.
Substituting (5) into (1)–(4), we obtain
(6)d2wdη2-εdθdηdwdη-eεθ(Redwdη+Haw-G)=0,(7)d2θdη2-RePrdθdη+EcPre-εθ(dwdη)2+EcPrHaw2=0,
with the boundary conditions
(8)w(0)=0,dθdη(0)=Bi0(θ(0)-1),w(1)=0,dθdη(1)=-Bi1θ(1),
where G is the pressure gradient parameter,Re=Vh/υ (Reynolds number),Pr=υ/α (Prandtl number),Ec=V2/cp(Tf-T∞) (Eckert number),Ha=σB02h2/μ0 (magnetic field parameter or square of Hartmann number),ε=m(Tf-T∞) (variable viscosity parameter),Bi0=γ0h/k (lower plate Biot number),Bi1=γ1h/k (upper plate Biot number).It is important to note that ε=0 corresponds to the case of constant viscosity conducting fluid. The exact solution of (6) for the fluid velocity is possible under this constant viscosity scenario and we obtain
(9)w(η)=GHa[eαη(eβ-1)-eβη(eα-1)eα-eβ+1],
where α=(Re+Re2+4Ha)/2 and β=(Re-Re2+4Ha)/2. Moreover, the coupled nonlinear boundary value problem represented by (6)-(7) together with their boundary conditions in (8) has been solved numerically using an efficient fourth-order Runge–Kutta method along with a shooting technique [19].


## 3. Entropy Analysis

In many engineering and industrial processes, entropy production destroys the available energy in the system. It is therefore imperative to determine the rate of entropy generation in a system, in order to optimize energy in the system for efficient operation in the system. The convection process in a channel is inherently irreversible and this causes continuous entropy generation. Wood [11] gave the local volumetric rate of entropy generation for a viscous incompressible conducting fluid in the presence of magnetic field as follows:
(10)EG=kT∞2(dTdy)2+μT∞(dudy)2+σB02T∞u2.
In (10), the first term represents irreversibility due to heat transfer; the second term is entropy generation due to viscous dissipation, while the third term is local entropy generation due to the effect of the magnetic field (Joule heating or Ohmic heating). Using (5), the dimensionless form of local entropy generation rate in (10) is given as follows:
(11)Ns=T∞2h2EGk(Tf-T∞)2=(dθdη)2+BrΩ[e-εθ(dwdη)2+Haw2],
where Ω=(Tf-T∞)/T∞ is the temperature difference parameter and Br=EcPr is the Brinkmann number. The Bejan number (Be) is defined as
(12)Be=N1Ns=11+Φ,
where Ns=N1+N2,  N1=(dθ/dη)2 (heat transfer irreversibility due to heat transfer), N2=(Br/Ω)[e-εθ(dw/dη)2+Haw2] (fluid friction and magnetic field irreversibility), Φ=N2/N1(irreversibility ratio).The Bejan number (Be) as shown in (12) has a range of 0≤Be≤1. If Be=0, then the irreversibility is dominated by the combined effects of fluid friction and magnetic fields, but if Be=1, then the irreversibility due to heat transfer dominates the flow system by the virtue of finite temperature differences.


## 4. Results and Discussion

The numerical results for the fluid velocity, temperature, entropy generation rate, and Bejan number distributions are reported in Table1 and Figures 2–30. Representative values of various parameters are utilized and the Prandtl number Pr is assumed to range from 0.71 (Air) ≤ Pr ≤ 7.1 (water). In order to validate the accuracy of our numerical procedure, we compare a special case of our result (ε=0) with the exact solution for the velocity profile in (9). The results displayed in Table 1 show perfect agreement and attest the correctness of our results.Table 1
Comparison between the exact and numerical solution of velocity profile forG=1, Re=1, Ha=1, and ε=0.
η
Exact solutionw(η)
Numerical solutionw(η)
0
0
0
0.1
0.035822
0.035822
0.2
0.065264
0.065264
0.3
0.087963
0.087963
0.4
0.103449
0.103449
0.5
0.111127
0.111279
0.6
0.110257
0.110257
0.7
0.099930
0.099930
0.8
0.079042
0.079042
0.9
0.046264
0.046264
1.0
0
0Figure 2
Velocity profiles increasing Ha.Figure 3
Velocity profiles increasingε.Figure 4
Velocity profiles with increasing Re.Figure 5
Velocity profiles with increasing Ec.Figure 6
Velocity profiles increasing Pr.Figure 7
Velocity profiles with increasing Bi0.Figure 8
Velocity profiles with increasing Bi1.Figure 9
Temperature with increasing Ha.Figure 10
Temperature with increasingε.Figure 11
Temperature with increasing Re.Figure 12
Temperature with increasing Ec.Figure 13
Temperature with increasing Pr.Figure 14
Temperature with increasing Bi0.Figure 15
Temperature with increasing Bi1.Figure 16
Ns with increasing Ha.Figure 17
Ns with increasing Re.Figure 18
Ns with increasingε.Figure 19
Ns with increasing BrΩ−1.Figure 20
Ns with increasingG.Figure 21
Ns with increasing Bi0.Figure 22
Ns with increasing Bi1.Figure 23
Bejan number with increasing Ha.Figure 24
Bejan number with increasing Re.Figure 25
Bejan number with increasingε.Figure 26
Bejan number with increasing BrΩ−1.Figure 27
Bejan number with increasing Bi1.Figure 28
Bejan number with increasing Bi0.Figure 29
Bejan number with increasingG.Figure 30
Bejan number with increasing Pr.


### 4.1. Effects of Parameter Variations on Velocity Profiles

The effects of variation in key parameters on the velocity profiles are shown in Figures2–8. Generally, the velocity profiles are parabolic in geometries with zero values at the channel walls due to no slip condition and attain their maximum value within channel. In Figure 2, it is observed that the fluid velocity decreases with increasing magnetic field intensity (Ha). This can be attributed to the presence of Lorentz force acting as a resistance to the flow as expected and is in perfect agreement with earlier results as reported in the literature [3–9]. As the fluid viscosity decreases as shown in Figure 3 with increasing values of ε, the velocity profiles increase. The fluid velocity decreases and skews towards the upper plate as Reynolds number (Re) increases due to increasing injection at the lower plate and increasing suction at the upper plate as shown in Figure 4. Figures 5 and 6 show that the fluid velocity increases with increase in the values of Eckert number (Ec) and Prandtl number (Pr). As Ec increases, the velocity gradient increases as a result of a decrease in the fluid viscosity, consequently, the fluid velocity increases. Figures 7 and 8 show that the fluid velocity increases with increasing convective heating (Bi0) at the lower plate and decreases with increasing convective cooling (Bi1) at the upper plate. This is expected, since the fluid becomes lighter and flows faster with increasing temperature due to convective heating.


### 4.2. Effects of Parameter Variations on Temperature Profiles

Figures9–15 demonstrate the effects of various parameters on the temperature profiles. The imposed thermal boundary conditions ensure that the fluid temperature at the lower plate is highest due to convective heating and decreases gradually to its lowest value at the upper plate due to convective heat loss to the ambient. Figure 9 shows the influence of magnetic field (Ha) on the flow field. As Ha increases due to increasing magnetic field intensity, the fluid temperature decreases within the channel. This decrease in the fluid temperature may be attributed to the combined effects of fluid suction and convective heat loss, despite the presence of Ohmic heating (or Lorentz heating) which serves as additional heat source to the flow system. The effects of increasing ε, Re, and Ec are shown in the Figures 10, 11, and 12. The rise in the fluid temperature is observed with increasing values of these parameters. This may be attributed to the facts that as ε, Re, and Ec increase, the fluid viscosity becomes lighter and viscous heating increases due to increasing convective heating at the lower plate increases leading to a rise in the fluid temperature. As the Prandtl number increases from Pr=0.71 (Air) to Pr=7.1 (water) the fluid temperature decreases as illustrated in Figure 13. In Figure 14, a rise in the fluid temperature is observed with increasing convective heating at the lower plate as expected. Figure 15 shows the effect of increasing Bi1 on the temperature. As expected, the fluid temperature decreases due to increasing convective heat loss at the upper plate.


### 4.3. Effects of Parameter Variations on Entropy Generation Rate

The effects of key parameters variation on entropy generation rate (Ns) are shown in Figures16–22. Generally, the entropy production is more pronounced at the permeable channel walls and decreases towards the channel centerline region. Figure 16 reveals the effect of increasing magnetic field intensity (Ha) on entropy generation rate (Ns). As Ha increases, the entropy generation decreases at the walls and increases at the centerline region of the channel. Meanwhile, it is interesting to note that two points exist, that is, η=0.3 and η=0.7 within the flow field where the entropy production is not affected by increasing Ha. In Figure 17, it is observed that entropy production is enhanced with increasing suction (Re) at the upper wall region, while a decrease in the entropy generation occurs at the lower wall with increasing fluid injection. The effects of ε, BrΩ−1 and G on the entropy generation rate are shown in Figures 18–20. As the fluid viscosity decreases with increasing values of ε, the entropy production increases at both walls and decreases towards the channel centerline as indicated in Figure 18. However, a point exists at η=0.5, where entropy production is virtually zero. This may be attributed to the presence of zero velocity gradients at the region. A similar trend of entropy production is observed with increasing values of group parameter (BrΩ−1) and constant pressure gradient (G) as demonstrated in Figures 19 and 20. Figures 21 and 22 show that the entropy generation rate increases with combined increase in convective heating at the lower wall and convective cooling at the upper wall that is, as Bi0 and Bi1 increase.


### 4.4. Effects of Parameter Variations on Bejan Number

Figures23–30 illustrate the effects of different values of key parameters on Bejan number (Be). Generally, the Bejan number is highest along the channel centerline region with irreversibility due to heat transfer dominating the flow, while near the channel walls the fluid friction and magnetic field irreversibility dominate. Moreover, it is evident in Figure 23 that an increase in Ha results in a decrease of Be along the channel centerline. As Re increases (see Figure 24), the Bejan number decreases near the lower wall due to injection and increases toward the upper wall due to suction. Figures 25 and 26 show a general decrease in Bejan number with increasing parameter values of ε and BrΩ−1 due to a decrease in fluid viscosity and an increase in viscous dissipation irreversibility. In Figures 27 and 28, an increase in the dominant influence of heat transfer irreversibility is observed as the parameter values of Bi0 and Bi1 increase, consequently, the Bejan number increases. Hence, the convective thermal boundary conditions enhance the dominant effects of heat transfer irreversibility on the flow system. An increase in the pressure gradient parameter causes a decrease in the Bejan number within the channel leading to an increase in the irreversibility due to fluid friction as shown in Figure 29. Figure 30 shows that the Bejan number slightly decreases at the lower plate and slightly increases at the upper plate with increasing Prandtl number Pr.


## 4.1. Effects of Parameter Variations on Velocity Profiles

The effects of variation in key parameters on the velocity profiles are shown in Figures2–8. Generally, the velocity profiles are parabolic in geometries with zero values at the channel walls due to no slip condition and attain their maximum value within channel. In Figure 2, it is observed that the fluid velocity decreases with increasing magnetic field intensity (Ha). This can be attributed to the presence of Lorentz force acting as a resistance to the flow as expected and is in perfect agreement with earlier results as reported in the literature [3–9]. As the fluid viscosity decreases as shown in Figure 3 with increasing values of ε, the velocity profiles increase. The fluid velocity decreases and skews towards the upper plate as Reynolds number (Re) increases due to increasing injection at the lower plate and increasing suction at the upper plate as shown in Figure 4. Figures 5 and 6 show that the fluid velocity increases with increase in the values of Eckert number (Ec) and Prandtl number (Pr). As Ec increases, the velocity gradient increases as a result of a decrease in the fluid viscosity, consequently, the fluid velocity increases. Figures 7 and 8 show that the fluid velocity increases with increasing convective heating (Bi0) at the lower plate and decreases with increasing convective cooling (Bi1) at the upper plate. This is expected, since the fluid becomes lighter and flows faster with increasing temperature due to convective heating.


## 4.2. Effects of Parameter Variations on Temperature Profiles

Figures9–15 demonstrate the effects of various parameters on the temperature profiles. The imposed thermal boundary conditions ensure that the fluid temperature at the lower plate is highest due to convective heating and decreases gradually to its lowest value at the upper plate due to convective heat loss to the ambient. Figure 9 shows the influence of magnetic field (Ha) on the flow field. As Ha increases due to increasing magnetic field intensity, the fluid temperature decreases within the channel. This decrease in the fluid temperature may be attributed to the combined effects of fluid suction and convective heat loss, despite the presence of Ohmic heating (or Lorentz heating) which serves as additional heat source to the flow system. The effects of increasing ε, Re, and Ec are shown in the Figures 10, 11, and 12. The rise in the fluid temperature is observed with increasing values of these parameters. This may be attributed to the facts that as ε, Re, and Ec increase, the fluid viscosity becomes lighter and viscous heating increases due to increasing convective heating at the lower plate increases leading to a rise in the fluid temperature. As the Prandtl number increases from Pr=0.71 (Air) to Pr=7.1 (water) the fluid temperature decreases as illustrated in Figure 13. In Figure 14, a rise in the fluid temperature is observed with increasing convective heating at the lower plate as expected. Figure 15 shows the effect of increasing Bi1 on the temperature. As expected, the fluid temperature decreases due to increasing convective heat loss at the upper plate.


## 4.3. Effects of Parameter Variations on Entropy Generation Rate

The effects of key parameters variation on entropy generation rate (Ns) are shown in Figures16–22. Generally, the entropy production is more pronounced at the permeable channel walls and decreases towards the channel centerline region. Figure 16 reveals the effect of increasing magnetic field intensity (Ha) on entropy generation rate (Ns). As Ha increases, the entropy generation decreases at the walls and increases at the centerline region of the channel. Meanwhile, it is interesting to note that two points exist, that is, η=0.3 and η=0.7 within the flow field where the entropy production is not affected by increasing Ha. In Figure 17, it is observed that entropy production is enhanced with increasing suction (Re) at the upper wall region, while a decrease in the entropy generation occurs at the lower wall with increasing fluid injection. The effects of ε, BrΩ−1 and G on the entropy generation rate are shown in Figures 18–20. As the fluid viscosity decreases with increasing values of ε, the entropy production increases at both walls and decreases towards the channel centerline as indicated in Figure 18. However, a point exists at η=0.5, where entropy production is virtually zero. This may be attributed to the presence of zero velocity gradients at the region. A similar trend of entropy production is observed with increasing values of group parameter (BrΩ−1) and constant pressure gradient (G) as demonstrated in Figures 19 and 20. Figures 21 and 22 show that the entropy generation rate increases with combined increase in convective heating at the lower wall and convective cooling at the upper wall that is, as Bi0 and Bi1 increase.


## 4.4. Effects of Parameter Variations on Bejan Number

Figures23–30 illustrate the effects of different values of key parameters on Bejan number (Be). Generally, the Bejan number is highest along the channel centerline region with irreversibility due to heat transfer dominating the flow, while near the channel walls the fluid friction and magnetic field irreversibility dominate. Moreover, it is evident in Figure 23 that an increase in Ha results in a decrease of Be along the channel centerline. As Re increases (see Figure 24), the Bejan number decreases near the lower wall due to injection and increases toward the upper wall due to suction. Figures 25 and 26 show a general decrease in Bejan number with increasing parameter values of ε and BrΩ−1 due to a decrease in fluid viscosity and an increase in viscous dissipation irreversibility. In Figures 27 and 28, an increase in the dominant influence of heat transfer irreversibility is observed as the parameter values of Bi0 and Bi1 increase, consequently, the Bejan number increases. Hence, the convective thermal boundary conditions enhance the dominant effects of heat transfer irreversibility on the flow system. An increase in the pressure gradient parameter causes a decrease in the Bejan number within the channel leading to an increase in the irreversibility due to fluid friction as shown in Figure 29. Figure 30 shows that the Bejan number slightly decreases at the lower plate and slightly increases at the upper plate with increasing Prandtl number Pr.


## 5. Conclusion

In this paper, the effects of magnetic field on variable viscosity channel flow with suction/injection together with convective heating/cooling at the walls have been investigated. The nonlinear model problem is tackled numerically using shooting quadrature and fourth-order Runge-Kutta iteration scheme. Based on the results presented above, the following conclusions are deduced.(i)
An increase inε, Ec, G, Pr, and Bi0 increases the velocity profiles, while an increase in Ha and Bi1 decreases the velocity profile along channel centerline region.(ii)
An increase inε, Ec, G, Pr, Re, and Bi0 increases the temperature profiles, while an increase in Ha and Bi1 decreases the temperature profile.(iii)
An increase inε, Bi0, Bi1, G, and BrΩ−1 increases the entropy generation rate. An increase in Re decreases Ns at injection wall, while at suction wall Ns increases. An increase in Ha increases Ns at both walls but decreases Ns at centre of the channel.(iv)
An increase inε, G, BrΩ−1, and Ha decreases the Be with increasing fluid friction and magnetic field irreversibility. An increase in Pr decreases Be at injection wall but increases Be at suction wall. Meanwhile increase in Bi0, Bi1 increases the Be with increasing effect of heat transfer irreversibility.



---

*Source: 630798-2013-05-13.xml*