0% found this document useful (0 votes)
54 views5 pages

Lect 3 5

1) One dimensional steady state heat transfer occurs when temperature depends only on one spatial dimension, such as thickness. Under these conditions, heat flux is constant across any plane. 2) The temperature distribution in a slab is linear, following the equation T = C1x + C2. By applying analogous concepts from circuit theory like resistance, more complex multi-slab systems can be analyzed. 3) Various types of boundary conditions like convection, radiation, and contact resistance can be accounted for to model full heat transfer problems using a composite resistance approach.

Uploaded by

api-3737931
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
54 views5 pages

Lect 3 5

1) One dimensional steady state heat transfer occurs when temperature depends only on one spatial dimension, such as thickness. Under these conditions, heat flux is constant across any plane. 2) The temperature distribution in a slab is linear, following the equation T = C1x + C2. By applying analogous concepts from circuit theory like resistance, more complex multi-slab systems can be analyzed. 3) Various types of boundary conditions like convection, radiation, and contact resistance can be accounted for to model full heat transfer problems using a composite resistance approach.

Uploaded by

api-3737931
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

One Dimensional Steady Solutions-I One Dimensional Steady Solutions-II

• As a first class of problem, we now turn to one • This is equivalent to saying


dimensional analysis
dT
−k = cons tan t
• These are useful for large slabs, where transverse dx
dimensions are much larger than thickness
• The same can be very easily obtained from heat
• In such cases, Temperature is only a function of x. equation
• This implies q'' can also be only a function of x. ∂ (ρcT)  ∂  ∂T  ∂  ∂T  ∂  ∂T  
⇒ =  k  +  k  +  k   + q ′′′
• From thermodynamics, at steady state ∂t  ∂x  ∂x  ∂y  ∂y  ∂z  ∂z  
q''in = q''out for any control volume
d  dT   dT  ⇒ T = C1' x + C 2
• This is equivalent to saying that heat ⇒ k =0 ⇒k  = C1
flux is constant across any plane dx  dx   dx 

x • The temperature distribution is linear

One Dimensional Steady Solutions-III One Dimensional Steady Solutions-IV


• Heat flow is analogous to current, as it is same in
• Logically, if we had two slabs series of slabs
• The profile will be a piecewise straight line
• Temperature is analogous to voltage as it drives
• In this case too, the heat flux along any plane the current
will be constant T1 − Ti T − T2
L1 L2 i = q ′′A = k1A = k 2A i
T − Ti T − T2 L1 L2
q ′′ = k1 1 = k2 i
L1 L2 T1 − Ti T − T2
i=q= = i ⇒ q L1 + L2  = T −T
• Thus, knowing T1, T2, k1, k2, L1, and L2, Ti can L1 L2  k 1A k 2 A  1 2
k 1A k 2A
k1 k2 be found and hence the temperature profile
• Thus Resistance is equivalent to L/(kA)
• To simplify analysis of these problems, we
T1 Ti T2 can introduce a principle similar to ohm’s law • q can be computed once, T1, T2, k1, k2, L1, L2, and A are
known. Ti can then be computed once q is known

Convective Resistance Radiation Resistance


• If we also now consider convective boundary on either • If we had a radiation boundary on one side interacting with
sides ambient at T∞ T∞
(TS1 − TS2 ) = h A(T − T )
q = h1A(T∞1 − TS1 ) =
L 2 S2 ∞2 (
q = σεA T∞ 4 − TS 4 )
kA ( )
⇒ q = σε T∞ 2 + TS 2 (T∞ + TS )A(T∞ − TS ) Ts
(T − T ) (T − T ) (T − T∞ 2 )
q = ∞1 S1 = S1 S2 = S2 ⇒ q = h r A(T∞ − TS )
1 L 1
h 1A kA h 2A (
where h r = σε T∞ + TS (T∞ + TS )
2 2
)
• Thus convective resistance is 1/(hA) • Thus we can say radiation resistance can be written as
1/(hrA)
• In a restrictive case of convection and radiation boundary
with the same ambient temperature, the boundary
R eq = 1 +L + 1 resistance shall be 1/((hr+hc)A)
h 1A kA h 2A

1
Contact Resistance Overall Heat Transfer Coefficient
• When we had put slabs together, the interface is assumed
to have one temperature • When we have a series of slabs (contact resistance
neglected)
• In reality, the interface is complex and the heat transfer
q=
(T∞1 − T∞ 2 )
mode is complex too R eq
• It can be idealized to be have a contact resistance with a 1 L L L 1
R eq = + 1 + 2 + 3 +
finite jump in temperature as shown h 1 A k 1A k 2 A k 3 A h 2 A
• This is obtained experimentally and is tabulated in
• In the heat transfer literature
handbooks q
′′ ′′
qR t ,c = ∆Tcontact = AR t ,c = q R t ,c overall heat transfer coefficient,
A
U, is defined as
• Sometimes the handbooks q = UA(T∞1 − T∞ 2 )
list contact conductance
1
⇒ UA =
q = h contact A∆Tcontact R eq

Network Analysis One Dimensional Cylindrical Wall


• Although one • The heat equation in cylindrical coordinate for a
dimensional analysis is constant property system can be written as
very restrictive,
∂ (T )  1 ∂  ∂T  1 ∂ 2 T ∂ 2 T 
sometimes it is extended ρc = k r + 2 2
+ 2  + q ′′′
to a network ∂t  r ∂r  ∂r  r ∂φ ∂z 
• This methodology is fast • In steady state with no heat generation and if the
becoming obsolete due to system can be assumed such that T varies only along r
the progress in numerical we can write d  dT  dT
heat transfer  rk =0 ⇒ rk =C
dr  dr  dr
or ⇒ 2πrLq ′r′ = C = q r
*
⇒ rq ′r′ = C
This implies that the total heat crossing
L any radial plane in the cylinder is constant

Electrical Analogy-I Electrical Analogy-II


dT qr
⇒ 2πrL(− k ) = qr ⇒ −dT = dr
dr 2πrLk

Integrating between Ri and Ro, we have


To Ro R
qr q
∫ − dT = ∫ 2πrLk dr ⇒ Ti − To = 2πLk ln R
r o

Ti Ri i

Ro
ln
Ri
Thermal resistance ⇒ R cond =
2πLk

2
Computation of Temperature One dimensional Spherical wall
• To compute temperature at any radius r within a • The heat equation in spherical coordinate for a
cylindrical wall with To, Ti, Ro and Ri, we proceed constant property system can be written as
as follows ∂ (T ) 1 ∂  2 ∂T  1 ∂  ∂T  1 ∂  ∂T 
ρc = 2  kr + k +  k sin θ  + q ′′′
qr 2πrL( − k ) dT
dr ⇒ T − T = − r dT ln R o ∂t r ∂r  ∂r  r 2 sin 2 θ ∂φ  ∂φ  r 2 sin θ ∂θ  ∂θ 
Ti − To = ⇒ Ti − To = i o
2πLk 2πLk dr Ri
R Ro
ln o ln 1 ∂  2 ∂T  ∂T
Ri Ri ⇒  kr =0 ⇒ r2k =C ⇒ r 2 q ′r′ = C
• Now we integrate between Ri and R as follows: r 2 ∂r  ∂r  ∂r

Ti − To R dr T R
This implies that the total heat crossing
⇒ ln
R o R∫ r
= − ∫ dT
T −T Ri ⇒ 4πr 2 q ′r′ = C* = q r any radial plane in the sphere is constant
ln i Ti ⇒ i =
Ri Ti − To R
ln o
Ri

Electrical Analogy Conduction with Heat Generation


 dT  q r dr • We had analyzed one-dimensional cases in slab,
4πr 2 q ′r′ = q r ⇒ 4πr 2  − k  = qr ⇒ −dT = cylinder and sphere cases without heat generation
 dr  4πk r 2
Integrating between Ri and Ro, we have • In nuclear and other applications conduction in
To Ro Ro heat generating body is of important significance
qr dr T q r −1
⇒ − ∫ dT = ∫ ⇒ −T To = • The heat generation can be uniform (electrical
Ti 4πk Ri r2 i 4πk r R i
heating) or non-uniform (Radiation shields)
1 1 
 −  • The current crossing a plane is no longer a
r ro
⇒ Ti − To = q r  i  Thermal resistance constant
4πk 1 1 
 r − r 
• The heat generation affects the temperature
=  i 
T − Ti distribution
Similarly, temperature at ⇒
any R can be written as To − Ti  1 1 
 r − r 
 o i 

1-D Slab with Heat Generation-I 1-D Slab with Heat Generation-II
x=-L x=L dT −q′′′
T=Ts T=Ts If k can be taken as constant, then ⇒ = x + C1*
∂ (ρcT)  ∂  ∂T  ∂  ∂T  ∂  ∂T   dx k
=  k  +  k  +  k   + q′′′ −q ′′′ 2
∂t  ∂x  ∂x  ∂y  ∂y  ∂z  ∂z   Integrating once more ⇒T= x + C1* x + C 2
2k
d  dT  Boundary conditions imply
⇒ k  + q′′′ = 0 For -L≤x≤L −q ′′′ 2
dx  dx  At x = -L Ts = L + C1* (− L ) + C 2
2k
Boundary Conditions; At x=-L, T=Ts; At x=L, T=Ts −q′′′ 2
x At x = L ⇒ Ts = L + C1*L + C 2
d  dT  2k
⇒ k  = −q ′′′ We shall first restrict to constant
dx  dx  volumetric heat generation rate Subtracting the above equations leads to C1*=0
Integration of the above equation leads to q′′′ 2
Adding and simplifying leads to C 2 = Ts + L
2k
 dT  q′′′ 2
⇒ k  = −q ′′′x + C1
 dx 
⇒T=
2k
(
L − x 2 + Ts ) Parabolic distribution x

3
1-D Slab with Heat Generation-III 1-D Slab with Heat Generation-IV
• Now if we add a clad to it we can just modify the • The temperature drop in clad is equal to
results in the following manner t q ′′′Lt q ′′′Lt
∆Tclad = q ′′′AL = ⇒ Tci = Tco +
kcA kc kc
• Now the surface of the clad temperature
will be the boundary condition for the • The temperature anywhere inside the clad can be
composite geometry computed using linear interpolation as temperature
• Since no heat is generated in the clad, profile is linear in the clad
the heat rate crossing every plane will be • If we neglect the contact resistance,
constant and we can use the previous then Ts of fuel is same as that of Tci.
concepts derived • The fuel temperature can be now estimated
• The total heat rate crossing any clad plane is equal from the equation derived previously
q ′′′ 2
to heat rate generated in one half of the fuel
⇒T=
2k f
(
L − x 2 + Ts)
q = q ′′′AL

1-D Slab with Heat Generation-V 1-D Slab with Heat Generation-VI
• If we have clad, contact resistance and fluid
• The centerline temperature To can be computed as cooling the fuel, then we can simply write
q ′′′Lt q′′′L2
⇒ To = Tco + +
kc 2k f
 1 t R ′′ L 
• We can manipulate the above to get some interesting ⇒ To − T∞ = q ′′′LA + + contact + 

 hA k c A A 2Ak f 
results
 t L 
⇒ To − Tco = q ′′′LA + 
 • The Temperature profiles can be found in fuel
 k c A 2Ak f  once Ts is found, by using
∆Toverall q Rclad Rfuel
q ′′′ 2
• This method is valid only to get To and valid only
⇒T=
2k f
(
L − x 2 + Ts)
if kf is constant

1-D cylinder with heat addition-I 1-D cylinder with heat addition-II
• The analysis of cylinder is similar r=R First Boundary condition implies C1 = 0
T=Ts
∂ (ρcT) 1 ∂  ∂T  1 ∂  ∂T  ∂  ∂T  dT r dT r
k  + k ⇒k = −q ′′′ ⇒ = −q ′′′
=  rk +  + q ′′′
∂t r ∂r  ∂r  r 2 ∂φ  ∂φ  ∂z  ∂z  dr 2 dr 2k
If k can be taken as constant, and −q′′′ 2
1 d  dT  ⇒T= r + C2
⇒  rk  + q ′′′ = 0 For 0≤r≤R we integrate once more 4k
r dr  dr 
Second Boundary condition −q ′′′ 2
Boundary Conditions; At r=0, dT/dr=0; At r=R, T=Ts implies Ts = R + C2
r 4k
d  dT  Restricting to constant q′′′ 2
⇒ r=0 Or C 2 = Ts + R
 rk  = −q′′′r volumetric heat generation rate 4k
dr  dr  dT/dr=0
q′′′ 2 2
Integration of the above equation leads to ⇒T=
4k
(
R − r + Ts )
2 2
dT r q ′′′R
⇒ rk = −q ′′′ + C1 ⇒ To = + Ts Parabolic distribution r
dr 2 4k

4
1-D cylinder with heat addition-III 1-D cylinder with heat addition-IV
• Now if we add a clad to it we can just modify the • If we have clad, contact resistance and fluid
results in the following manner cooling the fuel, then we can simply write
• The total heat rate crossing any clad radial  
ln  co 
R
 R ci  R ′contact
′ 
+ 
plane is 1 1 
⇒ To − T∞ = q ′′′πR ci 2 L + +
q = q′′′πR 2 L  h 2πR co L 2πLk c 2πR ci L 4πk f L 
 
 
• The temperature drop in clad is equal to
ln  o 
R
2  Ri  Note that L, being
∆Tclad = q ′′′πR L Clad temperature Profile
2πLk c arbitrary cancels out
R
ln o 
R ln
Ri  Tci − T R ci
⇒ Tci = Tco + q′′′R 2  =
Tci − Tco R
2k c ln co
R ci

1-D sphere with heat addition-I 1-D sphere with heat addition-II
• The analysis of sphere is similar First Boundary condition implies C1 = 0
dT r dT r
∂ (T ) 1 ∂  2 ∂T  1 ∂  ∂T  1 ∂  ∂T  ⇒k = −q ′′′ ⇒ = −q ′′′
ρc = 2  kr + k +  k sin θ  + q ′′′ dr 3 dr 3k
∂t r ∂r  ∂r  r 2 sin 2 θ ∂φ  ∂φ  r 2 sin θ ∂θ  ∂θ 
If k can be taken as constant, and −q′′′ 2
1 d  2 dT  ⇒T= r + C2
⇒ r k  + q ′′′ = 0 For 0≤r≤R we integrate once more 6k
r 2 dr  dr 
Second Boundary condition −q ′′′ 2
Boundary Conditions; At r=0, dT/dr=0; At r=R, T=Ts Ts = R + C2
implies 6k
R
q ′′′ 2
d dT  Restricting to constant Or C 2 = Ts + R q"'
⇒  r 2k  = −q′′′r
2
6k
dr  dr  volumetric heat generation rate
q ′′′ 2 2
Integration of the above equation leads to ⇒T=
6k
(
R − r + Ts )
3 2
dT r q ′′′R
⇒ r 2k = −q′′′ + C1 ⇒ To = + Ts Parabolic distribution
dr 3 6k

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy