Week 02
Week 02
Week 02
ILO 2: Student can apply Fourier’s law and Newton’s law in fin
theory and thermal resistance context. Further, the student can
analyse 2d heat transfer data in Matlab and formulate an energy
balance for 2d system.
Thermal resistance
Two examples of heat diffusion in 1d
●
Left: initially Gaussian temperature profile diffuses. Amplitude decreases and the
distribution spreads with time. The domain ends are insulated → heat does not escape
from the domain. qleft = qright = 0. Note: for fixed q bc T results.
●
Right: initially constant temperature object is heated from right end.
Temperature diffuses to the left end. Both ends are at fixed temperatures.
Tleft = 293K and Tright = 373K. Note: for fixed T bc q results.
Left and right end insulated i.e. dT/dx = 0 Left and right end fixed temperatures
Early time
Early
Two examples of heat diffusion in 1d
●
Left: initially Gaussian temperature profile diffuses. Amplitude decreases and the
distribution spreads with time. The domain ends are insulated → heat does not escape
from the domain. qleft = qright = 0. Note: for fixed q bc T results.
●
Right: initially constant temperature object is heated from right end.
Temperature diffuses to the left end. Both ends are at fixed temperatures.
Tleft = 293K and Tright = 373K. Note: for fixed T bc q results.
Left and right end insulated i.e. dT/dx = 0 Left and right end fixed temperatures
Fourier
q =k (Tright - Tleft)/L
Early
Derivation of steady state heat rate ([Q]=W=J/s) through a wall.
Convective heat transfer coeff. h (wind&indoor ventilation)
L
Sum: (1) + (3) and note that T2-T1 appears i.e. - (2)
T A−T B +T 2−T 1=Q /(h A A)+Q /(h B A )
Q=UA Δ T , with Δ T =T A −T B
A ΔT
Q=
Σi 1/(k i / Li )+ 1/(h A )+1/(h B)
Energy balance:
Δ T ave (t )
Q IN −Q OUT =c p m
Δt
Steady state:
Q IN =QOUT
Temperature distribution inside the fin in crucial role:
qout
x
qexit
qin T=T(x)
Local energy balance (J/s thinking)
for a single fin
Fourier’s law:
Energy conducts (J/s) into a small volume
dT ( x−dx / 2)
q in =−kA c
dx
A c=wt
Fourier’s law:
Energy conducts (J/s) out of a small volume
dT ( x +dx / 2)
q out =−kA c
dx
Newton’s law:
Energy exits (J/s) from fin to fluid
(strip of area dAs height dx, perimeter P=2(L+d))
q exit =hdA s (T −T ∞ )
dA s=2(w+t )dx
Local energy balance (J/s thinking)
for a single fin
q in −q out −q exit=0
When dx→ 0 we get the heat equation for T=T(x) in the fin but now
the equation has also a heat loss term as heat escapes to the fluid:
2
d T hP
2 = (T −T ∞ )
dx k Ac
2 hP
m = Definition of derivative
k Ac Δx
(
dx
−
dx
)=
2
1 dT (x +dx / 2) dT ( x−dx /2) d T ( x )
dx
2
, when Δ x →0
Incropera: 1d Temperature Distribution and
Heat Loss Along a Fin
2
d T hP
2 = (T −T ∞ )
dx k Ac
Figure 3.18 (Incropera): Conduction and convection in a fin of uniform cross section.
Example: Find temperature distribution and heat rate in a very
long copper rod (diameter D=5mm) with Tb=373K and T∞= 298K
and convection coefficient due to airflow h=100W/m2K
Temperature distribution
hP 4h
Estimate m: m=
√ k Ac
=
kD√≈14.2 m
−1
Heat rate
1) Fourier’s law
qf
ϵf = Usage of fins typically justified if > 2
hA c ,b θ b
Common assumption (not reality but useful)
→ assume that h is unaffected by 1) spatial position, and 2) presence of fins
Example: automobile radiator (fins on the air flow side, hot water in the inside)
A few examples from our own research
Example: liquid cooling an electric circuit by placing
a cooling plate with 3d printed finned microchannels on
top of the circuit
Note:
fins
c p ṁ Δ T =Power
ILO 2: Student can apply Fourier’s law and Newton’s law in fin
theory and thermal resistance context. Further, the student can
analyse 2d heat transfer data in Matlab and formulate an energy
balance for 2d system.
Consider heat conduction in 2d or 3d object (e.g. metal plate).
Divide the object into small elements and carry out energy balance analysis
for 1 of those elements. Assume: no heat losses.
Power exiting
[W/m²]
∆z
∆x
Power entering
[W/m²]
Derivation of heat equation:
Next we apply energy conservation law (“J/s thinking”)
for a small infinitesimal volume assuming conduction only
(e.g.1d metal rod)
x-∆x/2 x x+∆x/2
∆z Power exiting
Power entering
[W/m²] [W/m²]
∆y
∆x
Energy increase of element due to heat fluxes in x-direction during ∆t (J):
∂ T ( x +Δ x / 2 , y , t) ∂ T ( x−Δ x / 2 , y , t)
Δ Q x =[k −k ]Δ y Δ z Δ t
∂x ∂x
Energy increase of element due to heat fluxes in y-direction during ∆t (J):
∂ T ( x , y+Δ y/ 2 , t ) ∂T ( x , y−Δ y / 2 , t )
Δ Q y =[ k −k ]Δ x Δ z Δ t
∂y ∂y
Energy increase of element during ∆t (J):
ρ c p Δ T ( x ,t ) Δ x Δ y Δ z=Δ Q x +Δ Q y +Δ Q z
Then: Divide both sides by ΔxΔyΔzΔt and take the limit when all
Δ-variables → 0 → We get the heat equation.
Heat equation in 2d. Well, it is just thermal
energy conservation law.
General form Terms opened in 2d (assume α = constant)
∂T
=∇⋅α ∇ T ∂T ∂2 T ∂2 T
=α +α
∂t ∂t ∂x
2
∂y
2
What do the partial derivatives represent?
Mathematical interpretation?
General form Terms opened in 2d (assume α = constant)
∂T
=∇⋅α ∇ T ∂T ∂2 T ∂2 T
=α +α
∂t ∂t ∂x
2
∂y
2
∂T
<0
∂t
∂2 T ∂2 T
2
< 0 and 2
<0
In conduction
∂x ∂y
problems heat At a local minimum of a function:
conducts (diffuses)
from hot to cold. ∂2 T ∂2 T
2
> 0 and 2
>0
∂x ∂y
∂T
>0
∂t
On computers, we can solve heat equation by finite
difference methods. We discretize a 2d domain into small
elements.
x=0 x=Lx
T ni , j −1
y=0
Δy
Cell center
Δx
Finite difference discretizations
Δ T ni , j=T ni ,+1j −T ni , j
Numerical solution of temperature distribution
in a heated 2d metal plate
BC 1: BC 2:
- Cool sides and top - Insulated sides
- Hot base - Cool top and hot base
→ T=T(x,y,t) (2d) → T=T(y,t) (1d)
Numerical solution of temperature distribution
with two hot, one cold, and one insulated
boundary
Early time temperature Late time temperature
Zoom to plate
upper left corner
>> T(1:20,1:5)
ans =
269.8795 276.1205 280.3794 283.3261 285.3968
316.1205 309.8795 305.6206 302.6739 300.6032
314.0774 311.9226 309.8685 308.0042 306.3579
313.6455 312.3545 311.0793 309.8467 308.6758
313.4649 312.5351 311.6104 310.7006 309.8141
313.3662 312.6338 311.9038 311.1806 310.4685
313.3043 312.6957 312.0882 311.4845 310.8868
313.2623 312.7377 312.2138 311.6922 311.1742
313.2321 312.7679 312.3041 311.8418 311.3820
313.2095 312.7905 312.3717 311.9540 311.5380
313.1921 312.8079 312.4239 312.0406 311.6586
313.1784 312.8216 312.4651 312.1091 311.7541
313.1673 312.8327 312.4984 312.1645 311.8312
313.1581 312.8419 312.5257 312.2100 311.8948
313.1505 312.8495 312.5486 312.2481 311.9479
313.1440 312.8560 312.5681 312.2805 311.9932
313.1384 312.8616 312.5850 312.3086 312.0325
313.1334 312.8666 312.5999 312.3334 312.0672
313.1289 312.8711 312.6133 312.3557 312.0984
313.1248 312.8752 312.6256 312.3762 312.1270
>> T(1:20,1:5)
Ghost cell row of the top side BC
ans = Note:
The corner 269.8795 276.1205 280.3794 283.3261 285.3968
cell is 316.1205 309.8795 305.6206 302.6739 300.6032 1) the two
redundant 314.0774 311.9226 309.8685 308.0042 306.3579 values are
different →
313.6455 312.3545 311.0793 309.8467 308.6758 not insulated
313.4649 312.5351 311.6104 310.7006 309.8141 Boundary
313.3662 312.6338 311.9038 311.1806 310.4685
313.3043 312.6957 312.0882 311.4845 310.8868 2) the
Ghost cell average
313.2623 312.7377 312.2138 311.6922 311.1742
column of of the two
313.2321 312.7679 312.3041 311.8418 311.3820 values is
the left
313.2095 312.7905 312.3717 311.9540 311.5380 const.
side BC
313.1921 312.8079 312.4239 312.0406 311.6586 → fixed
313.1784 312.8216 312.4651 312.1091 311.7541 Ttop = 293K
313.1673 312.8327 312.4984 312.1645 311.8312
313.1581 312.8419 312.5257 312.2100 311.8948
313.1505 312.8495 312.5486 312.2481 311.9479
313.1440 312.8560 312.5681 312.2805 311.9932
313.1384 312.8616 312.5850 312.3086 312.0325
313.1334 312.8666 312.5999 312.3334 312.0672
313.1289 312.8711 312.6133 312.3557 312.0984
313.1248 312.8752 312.6256 312.3762 312.1270