Next Article in Journal
Gibbs Paradox and the Concepts of Information, Symmetry, Similarity and Their Relationship
Previous Article in Journal
Time in Quantum Measurement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermodynamic Analysis of Gravity-driven Liquid Film along an Inclined Heated Plate with Hydromagnetic and Viscous Dissipation Effects

1
Professional Training Institute, Saïd Otba, 30 000 Ouargla-Algeria
2
Faculty of Sciences and Engineering Sciences, University of Ouargla, 30 000 Ouargla- Algeria
*
Author to whom correspondence should be addressed.
Submission received: 21 December 2005 / Accepted: 13 September 2006 / Published: 23 October 2006

Abstract

:
The purpose of this work is to investigate the entropy generation in a laminar, gravity-driven conducting liquid film with fully developed velocity flowing along an incline heated plate in the presence of a transverse magnetic field. The upper surface of the liquid film is considered free and adiabatic. The effect of heat generation by viscous dissipation is included in the analysis. The influence of the applied magnetic field and the viscous dissipation on velocity, temperature and entropy generation is examined.

Introduction

Entropy generation is closely associated with thermodynamic irreversibility, which is encountered in all heat transfer processes. Different sources are responsible for generation of entropy such as heat transfer and viscous dissipation [1982; 1996]. The analysis of entropy generation rate in a circular duct with imposed heat flux at the wall and its extension to determine the optimum Reynolds number as function of the Prandtl number and the duty parameter were presented by Bejan [1996; 1979]. Sahin [1998] introduced the second law analysis to a viscous fluid in circular duct with isothermal boundary conditions. In another paper, Sahin [1999] presented the effect of variable viscosity on entropy generation rate for heated circular duct. A comparative study of entropy generation rate inside duct of different shapes and the determination of the optimum duct shape subjected to isothermal boundary condition were done by Sahin [1998]. Narusawa [1998] gave an analytical and numerical analysis of the second law for flow and heat transfer inside a rectangular duct. In a more recent paper, Mahmud and Fraser [2003; 2002] applied the second law analysis to fundamental convective heat transfer problems and to non-Newtonian fluid flow through channel made of two parallel plates. The study of entropy generation in a falling liquid film along an inclined heated plate was carried out by Saouli and Aïboud-Saouli [2004]. As far as the effect of a magnetic field on the entropy generation is concerned, Mahmud et al. [2003] studied the case of mixed convection in a channel.
The purpose of this article is thermodynamics analysis of a fully developed liquid film flowing along an inclined heated plate in the presence of a transverse magnetic field. The effect of heat generation by viscous dissipation is included in the analysis. Expressions for dimensionless velocity and temperature, entropy generation number are obtained.

Problem formulation and analytical solution

The problem as shown in Fig. 1 concerns a fully developed Newtonian, laminar, gravity-driven liquid film of thickness δ flowing along an inclined heated plate in the presence of a transverse uniform magnetic field B . The magnetic Reynolds number R e is assumed to be small, so that the induced magnetic field is neglected and the Hall effect of magnetohydrodynamics is ignored.
Figure 1. Schematic diagram of the problem.
Figure 1. Schematic diagram of the problem.
Entropy 08 00188 g001
Neglecting the inertia terms in the momentum equation compared to the body force and the magnetic term, the momentum equation is then:
μ 2 u ( y ) y 2 σ B 2 u ( y ) + ρ g s i n θ = 0
where μ is the dynamic viscosity, ρ the fluid density, σ the electric conductivity and g the gravity acceleration.
The boundary conditions are:
No-slip condition                                            u ( 0 ) = 0
Free surface                                                   u ( δ ) y = 0
The velocity profile is obtained by integrating Eq. (1) and using the boundary conditions given by Eq. (2). It may be written:
u ( y ) = ρ g s i n θ σ B 2 ( 1 c o s h ( B σ μ ( δ y ) ) c o s h ( B δ σ μ ) )
Introducing the following dimensionless variables for the velocity and the transverse distance U ( Y ) = u ( y ) u m , Y = y δ , the dimensionless velocity becomes:
U ( Y ) = c o s h ( H a ) c o s h ( H a ( 1 Y ) ) c o s h ( H a ) 1
where H a is the Hartman number defined as:
H a = B δ σ μ
and
u m = ρ g s i n θ σ B 2 ( c o s h ( H a ) c o s h ( H a ) 1 )
The liquid mass flow rate Q is:
Q = 0 δ ρ u ( y ) d y
Substituting Eq. (3) into Eq. (7), the liquid mass flow rate is:
Q = ρ 2 g s i n θ δ 3 μ [ 1 H a 2 1 H a 3 t a n h ( H a ) ]
In absence of the magnetic field, the liquid mass flow rate is [10]:
Q 0 = ρ 2 g s i n θ δ 3 3 μ
Therefore the liquid mass flow rate may be written:
Q Q 0 = 3 [ 1 H a 2 1 H a 3 t a n h ( H a ) ]
The energy equation for the present problem is:
u ( y ) T ( x , y ) x = λ ρ C P 2 T ( x , y ) y 2 + σ B 2 ρ C P u 2 ( y ) + μ ρ C P ( u ( y ) y ) 2
The boundary conditions are:
Inlet temperature                                     T ( 0 , y ) = T 0
Constant heat flux at the wall             λ T ( x , 0 ) y = q
Adiabatic upper surface                         T ( x , δ ) y = 0
Using the following dimensionless variables:
X = λ x ρ u m C P δ 2 , Y = y δ U ( Y ) = u ( y ) u m Θ ( X , Y ) = T ( x , y ) T 0 Δ T
where Δ T is a reference temperature difference defined as:
Δ T = q δ λ
The energy equation can be written in the following dimensionless form:
U ( Y ) Θ ( X , Y ) X = 2 Θ ( X , Y ) Y 2 + B r H a 2 U 2 ( Y ) + B r ( U ( Y ) Y ) 2
subjected to the following boundary conditions:
Θ ( 0 , Y ) = 0
Θ ( X , 0 ) Y = 1
Θ ( X , 1 ) Y = 0
B r = μ u m 2 C P 2 / λ Δ T is the Brinkman number.
To get a solution of Eq. (15), a separation of variables solution is assumed in the following form [10]:
Θ ( X , Y ) = Θ 1 ( X ) Θ 2 ( Y ) + Θ 1 ( X ) + Θ 2 ( Y )
The first term in the right-hand side of Eq. (17) is significant for decaying initial transition and entrance effect, the second term is significant for axial temperature rise due to accumulated wall heat flux and the third term is significant for transverse temperature variation to wall heat flux into fluid. Neglecting entrance effect and assuming that the system already passed the decaying initial transition. Then the first term at the right-hand side of Eq. (17) will disappear [8, 9]. Combination of Eq. (15) and Eq. (17) leaves two separated ordinary equations [10]. The solution of these two ordinary equations is:
Θ ( X , Y ) = α X + α ( c o s h ( H a ) 1 ) [ Y 2 2 c o s h ( H a ) c o s h ( H a ( 1 Y ) ) H a 2 ] B r H a 2 ( c o s h ( H a ) 1 ) 2 [ Y 2 2 c o s h 2 ( H a ) 2 H a 2 c o s h ( H a ) c o s h ( H a ( 1 Y ) ) + 1 4 H a 2 c o s h ( 2 H a ( 1 Y ) ) ] + C 1 Y + C
where α , C 2 and C are constants of integration.
Using the boundary conditions (16b) and (16c), it is found that:
α = A 3 A 4 A 1 A 2 C 1 = A 1 A 4 A 2 A 3 A 1 A 2
In the above expression A 1 , A 2 , A 3 and A 4 can be defined by:
A 1 = s i n h ( H a ) H a ( c o s h ( H a ) 1 ) , A 2 = c o s h ( H a ) c o s h ( H a ) 1 A 3 = B r H a 2 ( c o s h ( H a ) 1 ) 2 ( 2 H a c o s h ( H a ) s i n h ( H a ) 1 2 H a s i n h ( 2 H a ) ) 1 A 4 = B r H a 2 c o s h 2 ( H a ) ( c o s h ( H a ) 1 ) 2
To evaluate the constant of integration C , the bulk mean temperature given in Eq. (21) is examined:
Θ b ( X ) = 1 A A Θ ( X , Y ) d A
where the area element d A and the area A are:
d A = δ d Y A = 0 1 δ d Y = δ
Using Eq. (21), the bulk mean temperature is:
Θ b ( X ) = 0 1 Θ ( X , Y ) d Y
The boundary conditions defined by Eq. (16a) leads the following condition on the bulk mean temperature:
Θ b ( 0 ) = 0
Substituting Eq. (18) in Eq. (23) and using Eq. (24), the constant of integration is:
C = α ( c o s h H a 1 ) [ s i n h ( H a ) H a 3 ] B r H a 2 ( c o s h ( H a ) 1 ) 2 [ 2 H a 3 c o s h ( H a ) s i n h ( H a ) 2 16 H a 3 s i n h ( 2 H a ) ] α 6 ( c o s h ( H a ) 1 ) + B r H a 2 c o s h 2 ( H a ) 6 ( c o s h ( H a ) 1 ) 2 C 2 2
According to Woods [1975], the entropy generation rate is:
S G = λ T 0 2 [ ( T ( x , y ) x ) 2 + ( T ( x , y ) y ) 2 ] + μ T 0 ( u ( y ) y ) 2 + σ B 2 T 0 u 2 ( y )
The dimensionless entropy generation number is defined by the following relationship:
N S = λ T 0 2 q 2 S G
using the dimensionless velocity and temperature, Eq.(27) can be rewritten as:
N S = 1 P e 2 ( Θ ( X , Y ) X ) 2 + ( Θ ( X , Y ) Y ) 2 + B r Ω ( U ( Y ) Y ) 2 + B r H a 2 Ω U 2 ( Y )
N S = N C + N Y + N F + N B
where P e = ρ u m C P δ / λ and Ω = Δ T / T 0 are respectively the Peclet number and the dimensionless temperature difference. N C and N Y , are respectively the entropy generation numbers due to the conductive heat in the axial and the transverse directions. N F is the entropy generation number due to the fluid friction and N B is the entropy generation due to the hydromagnetic effect.

Results and discussion

The velocity profiles U ( Y ) are represented in Fig. 2 for various values of the Hartman number H a . As it can be seen, the action of the applied magnetic field B is to flatten the velocity profile near the free surface of the liquid film. The influence of the Hartman number on the liquid mass flow rate is illustrated in Fig. 3. As the Hartman number increases, the liquid mass fluid rate decreases. This means that the mean velocity of the liquid decreases. The application of the magnetic field induces a resistive force acting in the opposite direction of the flow, thus causing its deceleration.
Figure 2. Velocity Profiles as function of the transverse distance for different Hartman numbers.
Figure 2. Velocity Profiles as function of the transverse distance for different Hartman numbers.
Entropy 08 00188 g002
Fig. 4 presents a comparison between the temperature profiles Θ ( X , Y ) across the liquid film without the hydomagnetic effect, with hydromagnetic effect and with hydromagnetic and viscous dissipation effects. The temperature decreases along the transverse direction. For given transverse and axial directions, the temperature increases with the application of the magnetic field and the presence of viscous dissipation because of the heat generated by the hydromagnetic and viscous dissipation effects.
For a given axial distance X and Brinkman number B r , the effect of the Hartman number Ha on the temperature profile Θ ( X , Y ) is illustrated in Fig. 5. An increase of the Hartman number yields higher temperature profiles because of the heat dissipation due the action of the magnetic field.
Figure 3. Influence of the Hartman number on the liquid mass flow rate.
Figure 3. Influence of the Hartman number on the liquid mass flow rate.
Entropy 08 00188 g003
The action of the Brinkman number B r on the temperature profiles for given axial position X and Hartman number H a is illustrated in Fig. 6. As the Brinkman number increases, the temperature increases consequently because of the heat generated by viscous dissipation.
Figure 4. Comparison between the temperature profiles across the liquid film without the hydomagnetic effect, with hydromagnetic effect and with hydromagnetic and viscous dissipation effects.
Figure 4. Comparison between the temperature profiles across the liquid film without the hydomagnetic effect, with hydromagnetic effect and with hydromagnetic and viscous dissipation effects.
Entropy 08 00188 g004
Figure 5. Temperature profiles as function of the transverse distance at different Hartman numbers.
Figure 5. Temperature profiles as function of the transverse distance at different Hartman numbers.
Entropy 08 00188 g005
The effect of the Hartman number H a on the transverse distribution of the entropy generation number is plotted in Fig. 7. As the Hartman number increases the entropy generation number increases in the transverse direction and a minimum in the entropy generation number appears near the heated plate. At the upper surface where both velocity and temperature are maximum (or minimum) which cause zero velocity and temperature gradients leaving no contribution to the entropy generation number ( second and third term of Eq. (28)), the entropy generation number is most sensitive to the Hartman number which is proportional to the magnetic field. The presence of the magnetic field creates additional entropy (fourth term of Eq. (28)).
Figure 6. Temperature profiles as function of the transverse distance at different Brinkman numbers.
Figure 6. Temperature profiles as function of the transverse distance at different Brinkman numbers.
Entropy 08 00188 g006
Fig. 8 illustrates the effect of the Brinkman number B r , for fixed Hartman number and dimensionless group, on the transverse distribution of the entropy generation number, in which there is a minimum near the heated plate. For a given transverse position, the entropy generation number is higher for higher Brinkman number. The augmentation of the Brinkman number increases the contribution of the entropy generation number due to the conductive heat in the transverse direction because of the viscous dissipation. In all cases the heated plate acts as a strong source of irreversibility.
Figure 7. Entropy generation number as function of thetransverse distance at different Hartman number.
Figure 7. Entropy generation number as function of thetransverse distance at different Hartman number.
Entropy 08 00188 g007
Figure 8. Entropy generation number as function of thetransverse distance at different Brinkman numbers.
Figure 8. Entropy generation number as function of thetransverse distance at different Brinkman numbers.
Entropy 08 00188 g008
The effect of the dimensionless group B r Ω 1 , on the transverse distribution of the entropy generation number is depicted in Fig. 9. The dimensionless group determines the relative importance of viscous effect. For small dimensionless group, the entropy generation number decreases along the transverse distance. For higher dimensionless group, the entropy generation number decreases first, then increases with increasing transverse distance. For a given transverse position, the entropy generation number is higher for higher dimensionless group. This is due to the fact that for high dimensionless group, the entropy generation numbers due to the fluid friction and to the magnetic field increase (third and fourth term of Eq. (28)).
Figure 9. Entropy generation number as function of the transversedistance at different dimensionless group.
Figure 9. Entropy generation number as function of the transversedistance at different dimensionless group.
Entropy 08 00188 g009

Conclusion

This paper presents the application of the second law of thermodynamics to a gravity-driven liquid film along an inclined heated plate in the presence of a transverse magnetic field and viscous dissipation effects. The velocity and temperature profiles are obtained and used to evaluate the entropy generation number. The effects of the Hartman number, Brinkman number and the dimensionless group on velocity, temperature and entropy generation number are discussed.
From the results the following conclusions could be drawn:
Higher Hartman number causes flattened velocity profiles because the magnetic field slows down the movement of the fluid along the plate.
Temperature profiles shift to higher temperatures with increasing Hartman and Brinkman numbers because of the heat generated by magnetic and viscous dissipation.
The entropy generation number increases with Hartman number, Brinkman number and dimensionless group. As the Hartman number, Brinkman number and dimensionless group increases, the entropy generation number due respectively to the magnetic field, the conductive heat in the transverse direction and the fluid friction increases.

Nomenclature

A
area, (m2)
B
magnetic induction, (Wb.m-2)
B r
Brinkman number, μ u m 2 C P 2 / λ Δ T
C P
specific heat, (J.kg-1.K-1)
H a
Hartman number, B δ σ / μ
N B
entropy generation number, magnetic induction
N C
entropy generation, axial conduction
N F
entropy generation, fluid friction
N S
entropy generation number, total
N Y
entropy generation number, transverse conduction
P e
Peclet number, ρ u m C P δ / λ
q
wall heat flux, (W.m-2)
Q
liquid mass flow rate, (kg.m-1.s-1)
Q 0
liquid mass flow rate in absence of the magnetic field, (kg.m-1.s-1)
R e
Reynolds magnetic number, η σ u m δ
S G
entropy generation rate, (W.m-3.K-1)
T
temperature, (K)
u
axial velocity, (m.s-1)
U
dimensionless axial velocity
x
axial distance, (m)
X
dimensionless axial distance
y
transverse distance, (m)
Y
dimensionless transverse distance

Greek symbols

α
scalar constant
δ
thickness of the liquid film (m)
Δ T
reference temperature difference, Δ T = q δ λ
η
magnetic permeability, (H.m-1)
μ
dynamic viscosity, (kg.m-1.s-1)
λ
thermal conductivity, (W.m-1.K-1)
Θ
dimensionless temperature, ( T ( x , y ) T 0 ) / Δ T
Ω
dimensionless temperature difference, Δ T / T 0
ρ
density of the fluid, (kg.m-3)
σ
electric conductivity, (Ω-1.m-1)

Subscripts 

b
bulk value
m
maximum value
0
inlet value, reference value

References

  1. Bejan, A. Second-law analysis in heat transfer and thermal design. Adv. Heat Transfer. 1982, 15, 1–58. [Google Scholar]
  2. Bejan, A. Entropy generation minimization. CRC Press: New York; Boca Raton, 1996. [Google Scholar]
  3. Bejan, A. A study of entropy generation in fundamental convective heat transfer. J. Heat Transfer 1979, 101, 718–725. [Google Scholar] [CrossRef]
  4. Sahin, A.Z. Second law analysis of laminar viscous flow through a duct subjected to constant wall temperature. J. Heat Transfer 1998, 120, 76–83. [Google Scholar] [CrossRef]
  5. Sahin, A.Z. Effect of variable viscosity on the entropy generation and pumping power in a laminar fluid flow through a duct subjected to constant heat flux. Heat Mass Transfer 1999, 35, 499–506. [Google Scholar]
  6. Sahin, A.Z. A second law comparison for optimum shape of duct subjected to constant wall temperature and laminar flow. Heat Mass Transfer 1998, 33, 425–430. [Google Scholar]
  7. Narusawa, U. The second-law analysis of mixed convection in rectangular ducts. Heat Mass Transfer 1998, 37, 197–203. [Google Scholar] [CrossRef]
  8. Mahmud, S.; Fraser, R.A. The second law analysis in fundamental convective heat transfer problems. Int. J. Therm. Sci. 2003, 42, 177–186. [Google Scholar] [CrossRef]
  9. Mahmud, S.; Fraser, R.A. Thermodynamic analysis of flow and heat transfer inside channel with two parallel plates. Exergy 2002, 2, 140–146. [Google Scholar] [CrossRef]
  10. Saouli, S.; Aïboud-Saouli, S. Second law analysis of laminar falling liquid film along an inclined heated plate. Int. Communications in Heat and Mass Transfer 2004, 31, 879–886. [Google Scholar] [CrossRef]
  11. Mahmud, S.; Tasnim, S.H.; Mamun, S.A.A. Thermodynamic analysis of mixed convection in a channel with transverse hydromagnetic effect. Int. J. Therm. Sci. 2003, 42, 731–740. [Google Scholar] [CrossRef]
  12. Woods, L.C. Thermodynamics of fluid systems. Oxford University Press: Oxford, 1975. [Google Scholar]

Share and Cite

MDPI and ACS Style

Aïboud-Saouli, S.; Saouli, S.; Settou, N.; Meza, N. Thermodynamic Analysis of Gravity-driven Liquid Film along an Inclined Heated Plate with Hydromagnetic and Viscous Dissipation Effects. Entropy 2006, 8, 188-199. https://0-doi-org.brum.beds.ac.uk/10.3390/e8040188

AMA Style

Aïboud-Saouli S, Saouli S, Settou N, Meza N. Thermodynamic Analysis of Gravity-driven Liquid Film along an Inclined Heated Plate with Hydromagnetic and Viscous Dissipation Effects. Entropy. 2006; 8(4):188-199. https://0-doi-org.brum.beds.ac.uk/10.3390/e8040188

Chicago/Turabian Style

Aïboud-Saouli, Soraya, Salah Saouli, Noureddine Settou, and Nouredine Meza. 2006. "Thermodynamic Analysis of Gravity-driven Liquid Film along an Inclined Heated Plate with Hydromagnetic and Viscous Dissipation Effects" Entropy 8, no. 4: 188-199. https://0-doi-org.brum.beds.ac.uk/10.3390/e8040188

Article Metrics

Back to TopTop