9 files changed, 267 insertions, 0 deletions
diff --git a/1328/CH18/EX18.1/18_1.sce b/1328/CH18/EX18.1/18_1.sce
new file mode 100644
index 000000000..6d31d4ba8
--- /dev/null
+++ b/1328/CH18/EX18.1/18_1.sce
@@ -0,0 +1,65 @@
+printf("\t example 18.1 \n");
+// specific gravity of benzene is 0.88
+// specific heat of benzene is 0.48 Btu/(lb)*(F)
+U=50;
+A=400;
+T1=400;
+t1=100;
+t2=300;
+c=0.48;
+w=40000;
+C=0.60;
+W=10000;
+printf("\t values are approximately mentioned in the book \n");
+printf("\t for a \n");
+M=(7500*8.33*0.88);
+printf("\t weight of benzene is : %.1e lb \n",M);
+Q1=(w*c);
+printf("\t Q1 is : %.2e Btu/(hr)*(F) \n",Q1);
+Q2=(W*C);
+printf("\t Q2 is : %.0e Btu/(hr)*(F) \n",Q2);
+Ks=((%e)^(U*A*((1/Q1)-(1/Q2)))); // eq 18.16
+printf("\t Ks is : %.3f \n",Ks);
+Z=log((T1-t1)/(T1-t2));
+printf("\t Z is : %.3f \n",Z);
+theta=((M*(Z)*(Ks*6000-(19200)))/((Ks-1)*40000*6000));
+printf("\t theta is : %.1f hr \n",theta);
+printf("\t for b \n");
+R=(Q1/Q2);
+printf("\t R is : %.1f \n",R);
+KT=((%e)^(U*(A/Q1)*(1+R^2)^(1/2)));
+printf("\t KT is : %.0f \n",KT);
+S=((2*(KT-1))/((KT*(R+1+(1+R^2)^(1/2)))-(R+1-(1+R^2)^(1/2)))); // eq 18.24
+printf("\t S is : %.3f \n",S);
+theta1=((M*Z)/(0.266*40000)); // eq 18.25
+printf("\t theta1 is : %.2f hr \n",theta1);
+printf("\t for c \n");
+U1=100;
+A1=200;
+K8=((%e)^(U*(A/(2*Q1))*(1+R^2)^(1/2))); // eq 18.32
+S1=((2*(K8-1)*(1+((1-0.266)*(1-(3.2*0.266)))^(1/2)))/(((K8-1)*(3.2+1))+((K8+1)*(1+3.2^2)^(1/2)))); // eq 18.31
+printf("\t K8 is : %.2f \n",K8);
+printf("\t S1 is : %.3f \n",S1);
+theta2=((M*Z)/(0.282*40000)); // eq 18.25
+printf("\t theta2 is : %.2f hr \n",theta2);
+printf("\t for d \n");
+K9=((%e)^(U*(A/(Q1))*(R-1)));
+S2=((K9-1)/((K9*R)-1)); // eq 18.36
+printf("\t K9 is : %.2f \n",K9);
+printf("\t S2 is : %.2f \n",S2);
+t=100;
+t1=t+(S2*(T1-t)); // 18.37
+printf("\t t1 is : %.0f F \n",t1);
+t2=t1+(S2*(T1-t1));
+printf("\t t2 is : %.0f F \n",t2);
+t3=t2+(S2*(T1-t2));
+printf("\t t3 is : %.0f F \n",t3);
+t4=t3+(S2*(T1-t3));
+printf("\t t4 is : %.0f F \n",t4);
+x=0.23;
+printf("\t fractional circulation is : %.2f \n",x);
+N=3+x;
+printf("\t total fractional circulation : %.2f \n",N);
+theta3=(N*(M/w));
+printf("\t theta3 is : %.2f \n",theta3);
+// end
diff --git a/1328/CH18/EX18.2/18_2.sce b/1328/CH18/EX18.2/18_2.sce
new file mode 100644
index 000000000..1ed2ebc5d
--- /dev/null
+++ b/1328/CH18/EX18.2/18_2.sce
@@ -0,0 +1,21 @@
+printf("\t example 18.2 \n");
+tav=500; // F
+Ts=1000; 
+t0=100;
+c=0.12; // Btu/(lb)*(F)
+k=24; // Btu/(hr)*(ft^2)*(F/ft)
+row=488; // lb/ft^3
+alpha=0.41; // alpha=(k/(c*row)), ft^2/hr
+x=0.333; // ft
+theta=4;
+printf("\t values are approximately mentioned in the book \n");
+X=(x/(2*(alpha*theta)^(1/2)));
+printf("\t X is : %.2f \n",X);
+Y=0.142; // Y=f1(X) from fig 18.7 
+t=Ts+(t0-Ts)*(Y); // eq 18.43
+printf("\t t si : %.0f F \n",t);
+q=((k*(Ts-t0))/(3.14*alpha*theta)^(1/2)); // q=(Q/A),from eq 18.47
+printf("\t q is : %.0f Btu/(hr)*(ft^2) \n",q);
+q1=(2*k*(Ts-t0)*(theta/(3.14*alpha))^(1/2)); // q=(Q1/A). eq 18.49
+printf("\t The total heat which flowed through a square foot of wall in the 4 hr is : %.1e Btu/ft^2 \n",q1);
+// end
diff --git a/1328/CH18/EX18.3/18_3.sce b/1328/CH18/EX18.3/18_3.sce
new file mode 100644
index 000000000..7f1efed0d
--- /dev/null
+++ b/1328/CH18/EX18.3/18_3.sce
@@ -0,0 +1,12 @@
+printf("\t example 18.3 \n");
+Ts=1000; 
+t0=100;
+alpha=0.41; // alpha=(k/(c*row)), ft^2/hr
+theta=15/60;
+l=1; // ft
+X=(4*alpha*theta)/(l^2);
+printf("\t X is : %.2f \n",X);
+Y=0.155; // Y=f3*(X)from fig 18.9 when L=infinity
+t=Ts+(t0-Ts)*(Y); // eq 18.52
+printf("\t t si : %.1e F \n",t);
+// end
diff --git a/1328/CH18/EX18.4/18_4.sce b/1328/CH18/EX18.4/18_4.sce
new file mode 100644
index 000000000..ca2c87da4
--- /dev/null
+++ b/1328/CH18/EX18.4/18_4.sce
@@ -0,0 +1,41 @@
+printf("\t example 18.4 \n");
+T1=1100; // F
+T2=70; // F
+t1=T1+460; // R
+t2=T2+460; // R
+k=27; // from appendix
+c=0.14; // from appendix
+row=490; // from appendix
+alpha=0.394;
+theta=4;
+l=10/12; // ft
+x=0.173*10^(-8); // stefan constant
+e=0.7; // emmisivity
+printf("\t values are approximately mentioned in the book \n");
+printf("\t for a \n");
+// Assume the temperature is 500°F after 4 hr. The coefficient from plate to air is the· sum of the radiation and convection coefficients
+hri=(e*x*(t1^4-t2^4))/(T1-T2);
+printf("\t radiation coefficient is : %.1f Btu/(hr)*(ft^2)*(F) \n",hri); // eq 4.32
+hci=(0.3*(T1-T2)^(1/4)); // eq 10.10
+printf("\t convection coefficient is : %.1f Btu/(hr)*(ft^2)*(F) \n",hci);
+hti=hri+hci;
+printf("\t total intial coefficient is : %.1f Btu/(hr)*(ft^2)*(F) \n",hti);
+// For the 4-hr coefficient at 500°F
+hr=2.2; // Btu/(hr)*(ft^2)*(F)
+hc=1.35; // Btu/(hr)*(ft^2)*(F)
+ht=hr+hc;
+printf("\t total intial coefficient is : %.1f Btu/(hr)*(ft^2)*(F) \n",ht);
+h=(hti+ht)/2;
+printf("\t mean coefficient is : %.1f Btu/(hr)*(ft^2)*(F) \n",h);
+X=(4*alpha*theta)/(l^2);
+Y=(h*l)/(2*k);
+printf("\t X is : %.1f \n",X);
+printf("\t Y is : %.3f \n",Y);
+Z=0.42; // Z=f3(X,Y), from fig 18.10
+t=T2+((T1-T2)*Z); // eq 18.53
+printf("\t t is : %.0f F \n",t);
+printf("\t for b \n");
+Z1=0.43; // Z=f4(X,Y), from fig 18.11
+t1=T2+((T1-T2)*Z1); // eq 18.53
+printf("\t temperature of center plane is : %.0f F \n",t1);
+// end
diff --git a/1328/CH18/EX18.5/18_5.sce b/1328/CH18/EX18.5/18_5.sce
new file mode 100644
index 000000000..0aaaedb89
--- /dev/null
+++ b/1328/CH18/EX18.5/18_5.sce
@@ -0,0 +1,18 @@
+printf("\t example 18.5 \n");
+Ts=400;
+t0=200;
+k=25; // from appendix
+c=0.12; // from appendix
+row=490; // from appendix
+alpha=0.45; // alpha=(k/(c*row))
+theta=15/60;
+l=8/12; // ft
+h=50;
+X=(4*alpha*theta)/(l^2);
+Z=(2*k)/(h*l);
+printf("\t X is : %.2f \n",X);
+printf("\t Z is : %.1f \n",Z);
+Y=0.31; // Y=(Ts-t)/(Ts-t0), from fig 18.13
+t=Ts+(t0-Ts)*(Y); // eq 18.43
+printf("\t t is : %.0f F \n",t);
+//end
diff --git a/1328/CH18/EX18.6/18_6.sce b/1328/CH18/EX18.6/18_6.sce
new file mode 100644
index 000000000..382778d6a
--- /dev/null
+++ b/1328/CH18/EX18.6/18_6.sce
@@ -0,0 +1,52 @@
+printf("\t example 18.6 \n");
+Ts=300; 
+t0=70;
+c=0.25; // Btu/(lb)*(F)
+k=0.3; // Btu/(hr)*(ft^2)*(F/ft)
+row=103; // lb/ft^3
+alpha=0.01164; // alpha=(k/(c*row)), ft^2/hr
+theta=1;
+lx=9/12;
+ly=4.5/12;
+lz=2.5/12;
+h=4.1;
+printf("\t values are approximately mentioned in the book \n")
+X1=(4*alpha*theta)/(lx^2);
+Z1=(2*k)/(h*lx);
+printf("\t X1 is : %.4f \n",X1);
+printf("\t Z1 is : %.3f \n",Z1);
+X2=(4*alpha*theta)/(ly^2);
+Z2=(2*k)/(h*ly);
+printf("\t X2 is : %.4f \n",X2);
+printf("\t Z2 is : %.3f \n",Z2);
+X3=(4*alpha*theta)/(lz^2);
+Z3=(2*k)/(h*lz);
+printf("\t X3 is : %.3f \n",X3);
+printf("\t Z3 is : %.3f \n",Z3);
+printf("\t at centre (2*x/l) is zero \n");
+Yx=0.98; // fig 18.12
+Yy=0.75; // fig 18.12
+Yz=0.43; // fig 18.12
+printf("\t at surface (2*x/l) is one \n");
+Yx1=0.325; // fig 18.12
+Yy1=0.29; // fig 18.12
+Yz1=0.245; // fig 18.12
+printf("\t center of brick \n");
+t1=Ts-(Yx*Yy*Yz*(Ts-t0));
+printf("\t t1 is : %.1f F \n",t1);
+printf("\t corner of brick \n");
+t2=Ts-(Yx1*Yy1*Yz1*(Ts-t0));
+printf("\t t2 is : %.1f F \n",t2);
+printf("\t center of 9 by 4.5in face \n");
+t3=Ts-(Yx*Yy*Yz1*(Ts-t0));
+printf("\t t3 is : %.1f F \n",t3);
+printf("\t center of 9 by 2.5in face \n");
+t4=Ts-(Yx*Yy1*Yz*(Ts-t0));
+printf("\t t4 is : %.0f F \n",t4);
+printf("\t center of 4.5 by 2.5in face \n");
+t5=Ts-(Yx1*Yy*Yz*(Ts-t0));
+printf("\t t5 is : %.1f F \n",t5);
+printf("\t middle of long edge \n");
+t6=Ts-(Yx*Yy1*Yz1*(Ts-t0));
+printf("\t t6 is : %.0f F \n",t6);
+//end
diff --git a/1328/CH18/EX18.7/18_7.sce b/1328/CH18/EX18.7/18_7.sce
new file mode 100644
index 000000000..cda06e521
--- /dev/null
+++ b/1328/CH18/EX18.7/18_7.sce
@@ -0,0 +1,11 @@
+printf("\t example 18.7 \n");
+t=20; // min
+alpha=0.40; // ft^2/hr
+delx=0.167; // ft
+// From the conditions of Eq. (18.61) take time increments such that alpha(deltheta/delx^2)=1/2
+printf("\t approximate values are mentioned in the book \n");
+deltheta=(delx^2/(2*alpha));
+printf("\t deltheta is : %.3f hr \n",deltheta);
+N=(t/(deltheta*60));
+printf("\t number of steps required : %.1f \n",N);
+// end
diff --git a/1328/CH18/EX18.8/18_8.sce b/1328/CH18/EX18.8/18_8.sce
new file mode 100644
index 000000000..03350bed6
--- /dev/null
+++ b/1328/CH18/EX18.8/18_8.sce
@@ -0,0 +1,26 @@
+printf("\t example 18.8 \n");
+k=0.3;
+row=103;
+c=0.25;
+alpha=0.01164;
+f=1/24;
+t1=120;
+t2=60;
+printf("\t approximate values are mentioned in the book \n");
+printf("\t temperature lag 6in below the surface \n");
+x=6/12;
+theta=(x/2)*(1/(3.14*f*alpha))^(1/2); // eq 18.65
+printf("\t theta is : %.2f hr \n",theta);
+printf("\t amplitude \n");
+deltom=(t1-t2)/2;
+printf("\t deltom is : %.0f F \n",deltom);
+delt=(deltom)*(%e)^(-x*(3.14*f/alpha)^(1/2)); // eq 18.67
+printf("\t delt is : %.1f F \n",delt); // calculation mistake in book
+printf("\t temperature deviation after 2 hr \n");
+theta1=2; // hr
+deltx=(deltom)*((%e)^(-x*(3.14*f/alpha)^(1/2)))*cos((2*3.14*f*theta1)-(x*(3.14*f/alpha)^(1/2))); // eq 18.69
+printf("\t deltx is : %.1f F \n",deltx);
+printf("\t heat flow during the half period \n");
+q=(k*deltom*(2/(3.14*f*alpha))^(1/2)); // eq 18.70
+printf("\t heat flow is : %.0f Btu/(hr)*(ft^2) \n",q);
+// end
diff --git a/1328/CH18/EX18.9/18_9.sce b/1328/CH18/EX18.9/18_9.sce
new file mode 100644
index 000000000..9c2aa3e96
--- /dev/null
+++ b/1328/CH18/EX18.9/18_9.sce
@@ -0,0 +1,21 @@
+printf("\t example 18.9 \n");
+G=60; // lb/(hr)*(ft^2)
+De=1/12; // ft
+theta=6; // hr
+cs=41.3; // Btu/(ft^3)*(F)
+c=0.0191; // Btu/(ft^3)*(F)
+f=0.45; // void fraction
+T=90;
+T1=200;
+t0=50;
+h=(0.79*(G/De)^0.7); // eq 18.90
+printf("\t h is : %.1f \n",h);
+X=(h*theta/(cs*(1-f)));
+Y=(T-t0)/(T1-t0);
+printf("\t X is : %.0f \n",X);
+printf("\t Y is : %.3f \n",Y);
+row=0.0807; // lb/(ft^3) air
+Z=24.5; // Z=(h*x*row/(c*G)), by comparing X an Y in fig 18.21
+x=24.5*(c*G/(h*row));
+printf("\t x is : %.1f ft \n",x);
+// end