Lecture 3
Lecture 3
Lecture 3
Heat Transfer
Semester: January 2022
Dimensional Analysis
The equations of convective heat
Convection ..
transfer
1. Continuity equation u v w
+ + =0
x y z
2. Momentum equations
u u u u p 2u 2u 2u
+ u + v + w = − + 2 + 2 + 2 + Fx
t x y z x x y z
v v v v p 2v 2v 2v
+ u + v + w = − + 2 + 2 + 2 + Fy
t x y z y x y z
w w w w p 2w 2w 2w
+u +v + w = − + 2 + 2 + 2 + Fz
t x y z z x y z
3. Energy equation
T T T T 2T 2T 2T
C p +u +v +w =k 2 + 2 + 2 + q +
t x y z x y z
Three-dimensional, incompressible flow, Newtonian fluid with constant
properties in Cartesian coordinate system. [Pressure-velocity formulation]
Navier-Stokes and Energy
Convection ..
equations
➢ Three dependent variables
(
1. Velocity vector V = uiˆ + vjˆ + wkˆ )
2. Pressure (P)
3. Temperature (T)
➢ Four independent variables (Cartesian coordinate)
1. x
2. y (space)
3. z
4. t (time)
➢ Four known parameters (fluid properties)
1. ρ (density) Dimensional
2. (dynamic viscosity) Governing
3. k (thermal conductivity) parameters
4. Cp (specific heat at constant pressure)
All are dimensional quantities and therefore, the total number of inputs are 4.
Non-dimensional or Dimensionless
Convection ..
Parameters
➢ Some of the important parameters for convection problem
are listed below
1. Reynolds number (Re)
2. Prandtl number (Pr)
3. Grashof number (Gr) Governing parameters
4. Rayleigh number (Ra)
5. Peclet number (Pe)
6. Richardson number (Ri)
7. Nusselt number (Nu)
8. Biot number (Bi) Performance parameters
9. Stanton number (St)
Gr = =
viscous force 2
where,
= kinematic viscosity of fluid, m2/s
g = gravitational acceleration, m/s2
Franz Grashof
β = coefficient of volume expansion, 1/K
ΔT = temperature difference, K
Lref = characteristic or reference length of the geometry, m
Gr is a measure of the relative magnitudes of the buoyancy force
and the opposing viscous force acting on the fluid. The flow
regime in natural convection is governed by Grashof number. It
provides the main criterion in determining whether the fluid flow is
laminar or turbulent in natural convection.
Convection .. Rayleigh Number
( uU ) ( uV )
+ =0
( LX ) ( LY )
u U V
+ =0
L X Y
U V
+ =0
X Y
This is the non-dimensional continuity equation.
Dimensional Analysis:
Convection ..
Forced Convection
➢ x-momentum equation
u u p 2u 2u
u + v = − + 2 + 2
x y x x y
Using the definition, x = LX , y = LY ,u = uU ,v = uV , p = u2 P
1 2 2
Energy U +V = 2+ 2
X Y ReL Pr X Y
➢ Two independent variables (X, Y), four dependent variables (U,
V, P, Θ) and two known parameters or inputs (ReL and Pr).
Convection .. Boussinesq approximation
− T − T T − T
1 − β = 1/T∞ is the
T T volume expansion
1 − (T − T ) coefficient at
constant pressure.
T T 2T 2T
Energy C p u +v =k 2 + 2
x y x y
Dimensional Analysis: Free
Convection ..
or Mixed convection
➢ Change of nondimensional governing equations due to change
of reference scales:
Free Convection Case Reference Velocity, uref Reference Pressure, Pref
1 α / Lref ρα uref/ Lref
2 α / Lref ρν uref/ Lref
3 α / Lref uref g TLref
4 ν / Lref ρα uref/ Lref
5 ν / Lref ρν uref/ Lref
6 ν / Lref uref g TLref
7 g TLref ρα uref/ Lref
8 g TLref ρν uref/ Lref
9 g TLref uref g TLref
Mixed convection um ρu2ref
Dimensional Analysis:
Convection ..
Free convection (All Cases)
➢ Generalized non-dimensional governing equations become
U V
Continuity + =0
X Y
U U P 2U 2U
X-momentum ND U +V =− + ND 2 + 2
X Y X X Y
V V P 2V 2V
Y-momentum ND U +V =− + ND 2 + 2 + FND
X Y Y X Y
2 2
Energy U +V = k ND 2 + 2
X Y X Y
where
uref
2
uref Lref
ND = , ND = ,FND = g (T − Tref ) ,k ND =
Pref Pref Lref Pref uref Lref
Dimensional Analysis:
Convection ..
Free convection (case 1)
➢ Two independent variables (x, y), four dependent variables (u,
v, p, T) and four know parameters or inputs (ρ, , k, Cp).
➢ Let us choose the following reference or characteristic scales:
1. Reference length, Lref = L
2. Reference velocity, uref = α / Lref = α / L
3. Reference pressure, Pref = ραuref /L = ρα2/L2
4. Reference temperature, Tref = T∞
➢ Hence the non-dimensional parameters are
u 2
uref g TL3
ND = ref
= 1, ND = = Pr,k ND = =1 Gr =
Pref Pref Lref uref Lref 2
Lref Pr =
FND = g (T − Tref ) = Ra Pr ,
Pref
Ra = Gr Pr
Dimensional Analysis:
Convection ..
Free convection (case 1)
➢ Considering two-dimensional steady state problem with
buoyancy forces acting in y-direction, no internal heat generation
and negligible viscous dissipation, nondimensional governing
equations are U V
2 2
Energy U +V = + 2
X Y X 2
Y
➢ Two independent variables (X, Y), four dependent variables (U,
V, P, Θ) and two known parameters or inputs (Ra and Pr).
Dimensional Analysis:
Convection ..
Mixed convection
➢ Let us choose the following reference or characteristic scales:
1. Reference length, Lref = L
2. Reference velocity, uref = um
3. Reference pressure, Pref = ρu2ref
4. Reference temperature, Tref = T∞
➢ Hence the non-dimensional governing equations become
U V
Continuity + =0 g TL3
X Y Gr =
2
U U P 1 2U 2U
X-momentum U +V =− + 2+ 2 Pr =
X Y X ReL X Y
V V P 1 2V 2V um L
Y-momentum U +V =− + 2 + 2 + Ri Re =
X Y Y ReL X Y
1 2 2 Gr
+V = Ri = 2
Energy U 2+ 2 Re
X Y ReL Pr X Y
Dimensional Analysis: Forced
Convection ..
Convection
Question: What do you mean by dimensional analysis for forced
convective heat transfer?
Dimensional analysis for forced convective heat transfer is a
method which helps us to qualitatively relate the heat transfer
coefficient with fewer dimensionless variables that are required to
describe the forced convection problem.
For finding the heat transfer coefficient during forced convection,
the following factors are considered,
h = f(ρ, L, u∞, μ, Cp, k)
or, f(h, ρ, L, u∞, μ, Cp, k) = 0
No. of physical quantities :7 [h, ρ, L, u∞, μ, Cp, k]
No. of fundamental dimensions : 4 [MLTθ]
No. of nondimensional π-terms : 7 – 4 = 3
Buckingham π-theorem:
f(π1, π2, π3) = 0
Select three variables containing all four dimensions: ρ, h, C
Dimensional Analysis: Forced
Convection ..
Convection
So 1 = L −1u = Re
Dimensional Analysis: Forced
Convection ..
Convection
So 2 = Lk −1h = Nu
Dimensional Analysis: Forced
Convection ..
Convection
( L) ( MLT ) ( ML T ) ( LT ) ( )
j −1 k −1 l
−3 −1 −1
2 −2 −1
= LT
i 0 0 0 0
L T M
So 3 = k −1C p = Pr
f ( Re, Nu, Pr ) = 0 Nu = f ( Re, Pr )
Convection .. Problem #1: Dimensional Analysis
Given, 12 m/s
Bottom plate temperature, T1 = 15oC
Top plate temperature, T2 = 40oC T2 = 40C
Top plate speed, V = 12 m/s L = 0.7 mm Oil
Distance between plates, L = 0.7 mm
T1 = 15C
d 2u
x-momentum: 2
=0
dy
Integrating it twice gives u ( y ) = C1 y + C2
T T T T
2 2 u v v u
2 2
2
C p u +v = k 2 + 2 + 2 + + +
x y x y x y x y
2
T u
2
V y
2
d 2T
0 = k 2 + k 2 = − u( y) = V
y y dy L L
Dividing both sides by k and integrating twice give
y
2
T ( y) = − V + C3 y + C4
2k L
Applying the boundary conditions T(0) = T1 and T(L) = T2 gives
the temperature distribution to be
T2 − T1 V 2 y y 2
T ( y) = y + T1 + − 2
L 2k L L
Convection .. Problem #1: Dimensional Analysis
T2 − T1 V2 y y2
Tmax = T (0.00038) = y + T1 + − 2
L 2k L L
(40 − 15)C (0.605 N s/m 2 )(12 m/s) 2 0.00038 m (0.00038 m)2
= (0.00038 m) + 15C + − 2
0.0007 m 2(0.145 W/m C) 0.0007 m (0.0007 m)
= 103.12C
c) Heat flux at the plates is determined from the definition of heat
flux,
dT T2 − T1 V2 T2 − T1 V 2
q0 = −k = −k −k (1 − 0 ) = −k −
dy y =0
L 2kL L 2L
(40 − 15)C (0.605 N s/m 2 )(12 m/s) 2 1 W
= −(0.145 W/m.C) − = − 6.74 10 4
W / m2
dT T2 − T1 V2 T2 − T1 V 2
qL = − k = −k −k (1 − 2 ) = −k +
dy y=L
L 2kL L 2L
(40 − 15)C (0.605 N s/m 2 )(12 m/s) 2 1 W
= −(0.145 W/m.C) + = 5.705 10 W / m
4 2