diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /1328/CH13/EX13.2/13_2.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '1328/CH13/EX13.2/13_2.sce')
| -rw-r--r-- | 1328/CH13/EX13.2/13_2.sce | 59 |
1 files changed, 59 insertions, 0 deletions
diff --git a/1328/CH13/EX13.2/13_2.sce b/1328/CH13/EX13.2/13_2.sce new file mode 100644 index 000000000..387a23df0 --- /dev/null +++ b/1328/CH13/EX13.2/13_2.sce @@ -0,0 +1,59 @@ +printf("\t example 13.2 \n");
+printf("\t approximate values are mentioned in the book \n");
+printf("\t bubble point at 95F and 14.7psia \n");
+x(1)=0.077; // mole fraction of C4
+x(2)=0.613; // mole fraction of C5
+x(3)=0.310; // mole fraction of C6
+K(1)=3.13; // fig 7
+K(2)=0.92; // fig 7
+K(3)=0.3; // fig 7
+a(1)=3.4; // a= alpha
+a(2)=1;
+a(3)=0.326;
+i=1;
+while(i<4)
+ Z(i)=(a(i)*x(i));
+ i=i+1;
+end
+Zt=Z(1)+Z(2)+Z(3);
+printf("\t Zt is : %.3f \n",Zt);
+i=1;
+while(i<4)
+ y(i)=(a(i)*x(i)/(Zt));
+ printf(" \n x(i) K(i) a(i) Z(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(a(i))+" "+string(Z(i))+" "+string(y(i))+" \n");
+ i=i+1;
+end
+yt=y(1)+y(2)+y(3);
+printf("\t yt is : %.3f \n",yt);
+K2=(y(2)/x(2));
+printf("\t K2 is : %.3f \n",K2);
+printf("\t bubble point is 102 \n"); // from fig 7 , comparing K2 value
+printf("\t dew point at 130F and 14.7psia \n");
+y(1)=0.077; // mole fraction of C4
+y(2)=0.613; // mole fraction of C5
+y(3)=0.310; // mole fraction of C6
+K(1)=5; // fig 7
+K(2)=1.65; // fig 7
+K(3)=0.62; // fig 7
+a(1)=3.03; // a= alpha
+a(2)=1;
+a(3)=0.376;
+i=1;
+while(i<4)
+ Z(i)=(y(i)/a(i));
+ i=i+1;
+end
+Zt=Z(1)+Z(2)+Z(3);
+printf("\t Zt is : %.3f \n",Zt);
+i=1;
+while(i<4)
+ x(i)=(Z(i)/Zt);
+ printf(" \n y(i) K(i) a(i) Z(i) x(i) \n "+string(y(i))+" "+string(K(i))+" "+string(a(i))+" "+string(Z(i))+" "+string(x(i))+" \n");
+ i=i+1;
+end
+xt=x(1)+x(2)+x(3);
+printf("\t xt is : %.0f \n",xt);
+K2=(y(2)/x(2));
+printf("\t K2 is : %.2f \n",K2);
+printf("\t dew point is 122F \n"); // from fig 7, comparing K2 value
+// end
|