Scribd
Upload
42 views3 pages

PROBLEM C5.121: Y/2 Y/2 y y 2L 10 MM

This problem involves modeling the transient heat transfer between an acrylic and steel plate. The acrylic and steel plates are initially at different temperatures. Energy balances are written for control volumes within each material to determine the temperature distribution over time. The model is solved numerically to determine the time required for the acrylic surface temperature to reach its softening point and to plot the average temperatures of the acrylic and steel plates and the acrylic surface temperature versus time from 0 to 300 seconds.

Uploaded by

Sadie Hnatow
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
42 views3 pages

PROBLEM C5.121: Y/2 Y/2 y y 2L 10 MM

This problem involves modeling the transient heat transfer between an acrylic and steel plate. The acrylic and steel plates are initially at different temperatures. Energy balances are written for control volumes within each material to determine the temperature distribution over time. The model is solved numerically to determine the time required for the acrylic surface temperature to reach its softening point and to plot the average temperatures of the acrylic and steel plates and the acrylic surface temperature versus time from 0 to 300 seconds.

Uploaded by

Sadie Hnatow
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/ 3

PROBLEMC5.

121
KNOWN: Dimensions and properties of acrylic and steel plates. Initial temperatures.
FIND: Time needed to bring external surface of the acrylic to its softening temperature. Plot of the
average acrylic and steel plate temperatures and acrylic surface temperature for 0 s t s 300 s.
SCHEMATIC:
ASSUMPTIONS: (1) One-dimensional conduction, (2) Constant properties, (3) Negligible contact
resistance.
PROPERTIES: Acrylic (given):
A
= 1990 kg/m
3
, c
A
= 1470 J/kgK and k
A
= 0.21 W/mK. Steel
(given):
B
= 7800 kg/m
3
, c
B
= 500 J/kgK and k
B
= 45 W/mK.
ANALYSIS: We begin by writing energy balances on each of the 20 control volumes using the
implicit method,
Node 1:
( ) ( )
1 1 1
A 1 1 2 1
A A
2
p p p p
T T k T T
y
c
t y

+ + +

A
=
A A
Nodes 2 9:
( ) ( ) ( )
1 1 1 1 1
n n A A 1 1
A A
p p p p p p
n n n n
T T k T T k T T
c y
t y y

+ + + + +
+

A = +
A A A
Node 10:
( ) ( ) ( )
1 1 1 1 1
A 10 10 9 10 11 10
A A
"
p p p p p p
t
T T k T T T T
c y
t y R

+ + + + +

A = +
A A
Node 11:
( ) ( ) ( )
1 1 1 1 1
B 11 11 10 11 12 11
B B
"
p p p p p p
t
T T T T k T T
c y
t y R

+ + + + +

A = +
A A
Nodes 12 19:
( ) ( ) ( )
1 1 1 1 1
n n B B 1 1
B B
p p p p p p
n n n n
T T k T T k T T
c y
t y y

+ + + + +
+

A = +
A A A
Node 20:
( ) ( )
1 1 1
B 20 20 19 20
B B
2
p p p p
T T k T T
y
c
t y

+ + +

A
=
A A
Continued

Material B Material A
Adiabatic
surface
Adiabatic
surface
Interface
Ay/2 Ay/2 Ay
1
2
9 10 11 12 19 20

y
2L = 10 mm
PROBLEM5.121 (Cont.)
where
"
A B
/ 2 / 2
t
y y
R
k k
A A
= +
Also, note that the average temperature of each material may be written as
0
1
( )
L
y
T T y dy
L
=
=
}
or, in finite difference form,
A 1 2 3 4 5 6 7 8 9 10
A
1
2
y
T T T y T y T y T y T y T y T y T y T y
L
A (
= + A + A + A + A + A + A + A + A + A
(

and
B 11 12 13 14 15 16 17 18 19 20
B
1
2
y
T T y T y T y T y T y T y T y T y T y T
L
A (
= A + A + A + A + A + A + A + A + A +
(

The preceding equations were solved using IHT code that is listed in the Comments. The spatially-
averaged temperatures of the two plates, as well as the external temperature of the acrylic, T
1
, are
shown in the plot below. <
From the simulation, we also find that the surface of the acrylic reaches T
soft
= 90C at t = 87 s. <
COMMENTS: (1) Ultimately, the temperatures of the two plates will reach the same steady-state
value. The steady-state temperature may be found by recognizing the energy gained by the acrylic is
lost by the steel, L
A
c
A
(T
ss
T
i,A
) = L
B
c
B
(T
i,B
T
ss
) yielding T
ss
= 180C, as evident in the plot of the
average temperatures. (2) If the surface of the acrylic in contact with the metal is assumed to be of
constant temperature and equal to 300C, the external surface of the acrylic reaches the softening
temperature at t = 74 s. (3) The IHT code is shown on the next page.
Continued
Steel average temperature
Acrylic average
temperature
Acrylic surface temperature
Average Material and Acrylic Surface Temperatures
0 75 150 225 300
Time (s)
0
50
100
150
200
250
300
350
T
e
m
p
e
r
a
t
u
r
e
(
C
)
PROBLEM5.121 (Cont.)
//Geometry and Discretization
LA = 0.005 //m
LB = 0.005 //m
delyA = LA/9.5 //m
delyB = LB/9.5 //m
//Properties
kA = 0.21 //W/mK
kB = 45 //W/mK
cA = 1470 //J/kgK
cB = 500 //J/kgK
rhoA = 1990 //kg/m^3
rhoB = 7800 //kg/m^3
//Initial Temperatures
TiA = 20 //C
TiB = 300 //C
//Interface Resistance
Rt = delyA/kA/2 + delyB/kB/2
//Node 1
rhoA*cA*(delyA/2)*der(T1,t) = (kA/delyA)*(T2 - T1)
//Nodes 2 through 9
rhoA*cA*delyA*der(T2,t) = (kA/delyA)*(T1 - T2) + (kA/delyA)*(T3 - T2)
rhoA*cA*delyA*der(T3,t) = (kA/delyA)*(T2 - T3) + (kA/delyA)*(T4 - T3)
rhoA*cA*delyA*der(T4,t) = (kA/delyA)*(T3 - T4) + (kA/delyA)*(T5 - T4)
rhoA*cA*delyA*der(T5,t) = (kA/delyA)*(T4 - T5) + (kA/delyA)*(T6 - T5)
rhoA*cA*delyA*der(T6,t) = (kA/delyA)*(T5 - T6) + (kA/delyA)*(T7 - T6)
rhoA*cA*delyA*der(T7,t) = (kA/delyA)*(T6 - T7) + (kA/delyA)*(T8 - T7)
rhoA*cA*delyA*der(T8,t) = (kA/delyA)*(T7 - T8) + (kA/delyA)*(T9 - T8)
rhoA*cA*delyA*der(T9,t) = (kA/delyA)*(T8 - T9) + (kA/delyA)*(T10 - T9)
//Node 10
rhoA*cA*delyA*der(T10,t) = (kA/delyA)*(T9 - T10) + (1/Rt)*(T11 - T10)
//Node 11
rhoB*cB*delyB*der(T11,t) = (1/Rt)*(T10 - T11) + (kB/delyB)*(T12 - T11)
//Nodes 12 through 19
rhoB*cB*delyB*der(T12,t) = (kB/delyB)*(T11 - T12) + (kB/delyB)*(T13 - T12)
rhoB*cB*delyB*der(T13,t) = (kB/delyB)*(T12 - T13) + (kB/delyB)*(T14 - T13)
rhoB*cB*delyB*der(T14,t) = (kB/delyB)*(T13 - T14) + (kB/delyB)*(T15 - T14)
rhoB*cB*delyB*der(T15,t) = (kB/delyB)*(T14 - T15) + (kB/delyB)*(T16 - T15)
rhoB*cB*delyB*der(T16,t) = (kB/delyB)*(T15 - T16) + (kB/delyB)*(T17 - T16)
rhoB*cB*delyB*der(T17,t) = (kB/delyB)*(T16 - T17) + (kB/delyB)*(T18 - T17)
rhoB*cB*delyB*der(T18,t) = (kB/delyB)*(T17 - T18) + (kB/delyB)*(T19 - T18)
rhoB*cB*delyB*der(T19,t) = (kB/delyB)*(T18 - T19) + (kB/delyB)*(T20 - T19)
//Node 20
rhoB*cB*(delyB/2)*der(T20,t) = (kB/delyB)*(T19 - T20)
//Average Temperatures
TavgA*LA = T1*delyA/2 + T2*delyA + T3*delyA + T4*delyA + T5*delyA + T6*delyA + T7*delyA +
T8*delyA + T9*delyA + T10*delyA
TavgB*LB = T11*delyB + T12*delyB + T13*delyB + T14*delyB + T15*delyB + T16*delyB + T17*delyB
+ T18*delyB + T19*delyB + T20*delyB/2
//Steady State Temperature
LA*cA*rhoA*(TiA - Tss) = LB*cB*rhoB*(Tss - TiB)

You might also like