summaryrefslogtreecommitdiff
path: root/965/CH2/EX2.29/29.sci
blob: 3959570a50b69236293b997afc5f58360eeb4c47 (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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
clc;
clear all;
disp("interface temperature")
La=0.12;//m 
Lb=0.12;//m 
Lc=0.012;//m 

hcf=18;//W/(m^2*C)

kA=1.6;//W/(m*C) 
kB=0.3;//W/(m*C) 
kC=0.14;//W/(m*C)

Rair=0.16;//K/W

thf=1090;// degree C temperature 
tcf=20;// degree C temperature 

Rtotal =(La/kA)+Rair+(Lb/kB)+(Lc/kC)+1/hcf

q=(thf-tcf)/Rtotal;

disp("W",q,"rate of heat flow per m^2 surface area = ")


//q=(thf-t2)/(La/kA)
t2=thf-q*(La/kA);
disp ("degree C",t2,"t2 = ")

//q=(t2-t3)/(Rair)
t3=t2-q*(Rair);
disp ("degree C",t3,"t3 = ")

//q=(t3-t4)/(Lb/B)
t4=t3-q*(Lb/kB);
disp ("degree C",t4,"t4= ")


//q=(t4-t5)/(Lc/kC)
t5=t4-q*(Lc/kC);
disp ("degree C",t5,"t5 = ")