blob: ab297d0dff5bf2eb2a9aef129767b7ab2718030a (
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
|
clc();
clear;
// To find the temperature at the free end is made of copper iron and glass
D = 3/48; // diameter in ft
L = 9/12; // Length of steam vessel in ft
T1 = 210; // Vessel temperature in degF
T2 = 80; // Air temperature in degF
th0 = T1-T2; // Temperature difference in degF
h = 1.44; // Assumed heat coefficient in Btu/hr-ft^2-degF
C = %pi*D; // Circumference of vessel in ft
A = %pi*D*D/4; // Area of vessel in ft^2
// For copper
k1 = 219; // Heat conductivity of copper in Btu/hr-ft-degF
m1 = sqrt(h*C/(k1*A)); // in /ft
th1 = th0*2/(exp(m1*L)+exp(-m1*L));
Tl1 = round(th1+T2); // The temperaure at the free end in degF
printf("Temperature at free end of the copper rod is %d degF \n",Tl1);
// For iron
k2 = 36; // heat conductivity of copper in Btu/hr-ft-degF
m2 = sqrt(h*C/(k2*A)); // in /ft
th2 = th0*2/(exp(m2*L)+exp(-m2*L));
Tl2 = th2+T2; // The temperaure at the free end in degF
printf(" Temperature at free end of the iron rod is %.2f degF \n",Tl2);
// For glass
k3 = 0.64; // Heat conductivity of copper in Btu/hr-ft-degF
m3 = sqrt(h*C/(k3*A)); // in /ft
th3 = th0*2/(exp(m3*L)+exp(-m3*L));
Tl3 = th3+T2; // The temperaure at the free end in degF
printf(" Temperature at free end of the glass rod is %.2f degF \n",Tl3);
|