summaryrefslogtreecommitdiff
path: root/617/CH4/EX4.3/Example4_3.sci
blob: bc85a2de5c53a45689c98c315863594237a9f526 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
clc();
clear;

// To find the tempearure and heat low in case of sudden heat change

t = 10;              // time elapsed in hr
Ti= 70;              // tempearature of wall initially in degF
Ts = 1500;           // temperature of surface when suddenly changed in degF
a = 0.03;            // thermal diffusivity in ft^2/hr
k = 0.5;             // thermal conductivity in Btu/hr-ft-degF
A = 10;              // area of wall in sq ft
x = 7/12;            // distance from surface where tempearture is to be found in ft
f = x/(2*sqrt(a*t)); 
// From gaussian error function table erf can be found
errorf = 0.55;       // Referred from table

T = Ts+(Ti-Ts)*errorf;
printf("Temperaure at a distance of 7/12ft from surface is %.1f degF \n",T);
q = -k*A*(Ti-Ts)*exp(-x^2/(4*a*t))/sqrt(t*%pi*a);   // heat flow rate at a distance
qtot = -k*A*(Ti-Ts)*2*sqrt(t/(%pi*a));              // total heat flowing after 10 hrs in Btu
printf(" Heat flowing at a distance of 7/12 ft from surface is %d Btu/hr\n",q); 
printf(" Total heat flow after 10hrs is %f Btu",%pi);