summaryrefslogtreecommitdiff
path: root/965/CH2/EX2.8/8.sci
blob: 77ba284bcd4aab58e1b65b4e8747be38c3420b6d (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
clc;
clear all;
disp("temperature between layers")
La=0.22;//m
Lb=0.15;//m
Lc=0.05;//m
Ld=0.003;//m
kA=4;// kJ/(mh*C)
kB=2.8;// kJ/(mh*C)
kC=0.24;// kJ/(mh*C)
kD=240;// kJ/(mh*C)
t1= 1500;// degree C
t5=90;// degree C
RthA=La/kA;//(m^2)*h*C/kJ
RthB=Lb/kB;// (m^2)*h*C/kJ
RthC=Lc/kC;// (m^2)*h*C/kJ
RthD=Ld/kD;// (m^2)*h*C/kJ
Rtotal= RthA+RthB+RthC+RthD;// (m^2)*h*C/kJ
disp ("(m^2)*h*C/kJ",Rtotal,"total thermal resistance is = ")
q= (t1-t5)/Rtotal;// kJ/(h*m^2)
disp("kJ/(h*m^2)",q,"rate of heat loss per unit area of the wall is =")
// q= (t1-t2)/RthA=  (t3-t4)/RthB
t4=t5+q*RthD;// degree C
disp ("degree C",t4,"the temperature t4 is = ")
t3=t4+q*RthC;// degree C
disp ("degree C",t3,"the temperature t3 is = ")
t2=t3+q*RthB;// degree C
disp ("degree C",t2,"the temperature t2 is = ")