1 files changed, 65 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