9 files changed, 143 insertions, 0 deletions
diff --git a/2175/CH16/EX16.1/16_1.sce b/2175/CH16/EX16.1/16_1.sce
new file mode 100755
index 000000000..103cc3b9d
--- /dev/null
+++ b/2175/CH16/EX16.1/16_1.sce
@@ -0,0 +1,8 @@
+clc;
+lambda=10^3*0.52;
+x=250;
+t1=40;
+t2=20;
+q=lambda*(t1-t2)/x;
+disp("rate of heat transfer per unit area:");
+disp("W/m^2",q);
diff --git a/2175/CH16/EX16.12/16_12.sce b/2175/CH16/EX16.12/16_12.sce
new file mode 100755
index 000000000..2fb55524f
--- /dev/null
+++ b/2175/CH16/EX16.12/16_12.sce
@@ -0,0 +1,11 @@
+clc;
+delta_p=0.0002;//bar
+d=25;
+rho=7600;//assumed to run program
+c=1.13;
+C=24;
+tou=delta_p*10^5*d/(4*10^3);
+f=tou/(rho*C^2/2);
+alpha=0.125*rho*c*C/(rho*C^2);
+disp("heat transfer coefficient is:");;
+disp("kW/m^2 K",alpha);
diff --git a/2175/CH16/EX16.15/16_15.sce b/2175/CH16/EX16.15/16_15.sce
new file mode 100755
index 000000000..151c739d6
--- /dev/null
+++ b/2175/CH16/EX16.15/16_15.sce
@@ -0,0 +1,6 @@
+clc;
+delta_t=277-17;
+d=0.15;
+alpha=1.32*(delta_t/d)^0.25;
+disp("heat transfer coefficient=");
+disp("W/m^2 K",alpha);
diff --git a/2175/CH16/EX16.16/16_16.sce b/2175/CH16/EX16.16/16_16.sce
new file mode 100755
index 000000000..203bb32db
--- /dev/null
+++ b/2175/CH16/EX16.16/16_16.sce
@@ -0,0 +1,14 @@
+clc;
+Beta=1/303;
+g=9.81;
+l=1;
+delta_t=327-30;
+v=(5.128*10^-5);
+Gr=Beta*g*l^3*delta_t/v^2
+
+alpha=1.31*delta_t^0.33333
+A=1;//m^2
+delta_t=627-27;
+Q=alpha*A*delta_t
+disp("rate of heat loss:");
+disp("kW",Q/1000);
diff --git a/2175/CH16/EX16.18/16_18.sce b/2175/CH16/EX16.18/16_18.sce
new file mode 100755
index 000000000..23b7274c8
--- /dev/null
+++ b/2175/CH16/EX16.18/16_18.sce
@@ -0,0 +1,28 @@
+clc;
+m=3;
+rho=500;
+v=m/rho;
+
+l=4;//m
+r=0.01;
+A=%pi*r^2;
+n=v*l/A;
+disp("number of tumes is:");
+disp(n)
+
+alpha0=260;
+A0=12.7;
+alphai=580;
+Ai=10;
+U=1/[1/alpha0+A0/(alphai*Ai)];
+N=U*%pi*(A0/1000)*l*n/(3*1.5*1000);
+R=3*1.5/(40*1.04);
+
+eta=[1-%e^(-N*(1-R))]/[1-R*%e^(-N*(1-R))]
+disp(eta,"eta is:");
+
+t2=400;
+t1=100;
+tL=eta*(t2-t1)+t1
+disp("exit temperature is :");
+disp(tL);
diff --git a/2175/CH16/EX16.2/16_2.sce b/2175/CH16/EX16.2/16_2.sce
new file mode 100755
index 000000000..2256b536f
--- /dev/null
+++ b/2175/CH16/EX16.2/16_2.sce
@@ -0,0 +1,17 @@
+clc;
+alpha_a=2800;
+lambda=10^3*50;
+x=10;
+alpha_b=11;
+U=1/[1/alpha_a+x/lambda+1/alpha_b];
+
+tA=90;
+tB=15;
+q=(tA-tB)*U;
+disp("rate of heat lost per sq m of surface")
+disp("kW",q)
+
+//part b
+t2=q/alpha_b+tB;
+disp("temperature of outsede surface:");
+disp("C",t2)
diff --git a/2175/CH16/EX16.21/16_21.sce b/2175/CH16/EX16.21/16_21.sce
new file mode 100755
index 000000000..e74fe45c2
--- /dev/null
+++ b/2175/CH16/EX16.21/16_21.sce
@@ -0,0 +1,16 @@
+clc;
+eta=0.4;
+sigma=5.67;
+T1=13.73;
+T2=3.13;
+q=eta*sigma*(T1^4-T2^4);
+disp("heat loss by radiation is:");
+disp("kW",q/1000);
+
+eta2=0.9;
+q1=eta*sigma*T1^4;
+q2=eta2*sigma*T2^4
+q_=q1-q2;
+disp("grey body assumptions overstimates by:");
+pct=(q-q_)/q_
+disp("%",pct*100)
diff --git a/2175/CH16/EX16.25/16_25.sce b/2175/CH16/EX16.25/16_25.sce
new file mode 100755
index 000000000..86d80ec25
--- /dev/null
+++ b/2175/CH16/EX16.25/16_25.sce
@@ -0,0 +1,17 @@
+clc;
+eta=0.8;
+F1_2=5.67*10^-8;
+T1=533;//K
+T2=293;//K
+alpha=eta*F1_2*(T1^2+T2^2)*(T1+T2);
+
+A=%pi*0.6*0.9;
+Q1=alpha*A*(T1-T2);
+
+alpha=8.8;
+A=5;
+Q2=alpha*A*(T1-T2);
+
+Q=Q1+Q2;
+disp("total heat loss is:");
+disp("kW",Q/1000)
diff --git a/2175/CH16/EX16.7/16_7.sce b/2175/CH16/EX16.7/16_7.sce
new file mode 100755
index 000000000..5b8c371d6
--- /dev/null
+++ b/2175/CH16/EX16.7/16_7.sce
@@ -0,0 +1,26 @@
+clc;
+alpha=88.8;
+L=0.05;
+lambda=40;
+Bi=alpha*L/lambda;
+//p1L*[cos(p1L)/sin(p1L)]=1-Bi;
+//from trial and error;
+p1L=0.57;
+
+tou=20*60;
+rho=7600;
+c=0.5*10^3;
+R=0.05;
+F0=lambda*tou/(rho*c*R^2);
+
+tF=20;
+ti=500;
+a=(sin(p1L)-p1L*cos(p1L))*(2*%e^[-(p1L)^2*F0])/(p1L-sin(p1L)*cos(p1L));
+tc=tF+a*(ti-tF);
+disp("temperature of center is:")
+disp("C",tc)
+
+//part b
+tc=tF+[%e^(-3*Bi*F0)]*(ti-tF)
+disp("temperature of center by newtonian cooling is:")
+disp("C",tc)