summaryrefslogtreecommitdiff
path: root/1328/CH13/EX13.1/13_1.sce
diff options
context:
space:
mode:
Diffstat (limited to '1328/CH13/EX13.1/13_1.sce')
-rw-r--r--1328/CH13/EX13.1/13_1.sce207
1 files changed, 207 insertions, 0 deletions
diff --git a/1328/CH13/EX13.1/13_1.sce b/1328/CH13/EX13.1/13_1.sce
new file mode 100644
index 000000000..18e1019e8
--- /dev/null
+++ b/1328/CH13/EX13.1/13_1.sce
@@ -0,0 +1,207 @@
+printf("\t example 13.1 \n");
+// at atmospheric pressure,Pt=760 mm Hg
+printf("\t approximate values are mentioned in the book \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
+printf("\t for T 100 F \n");
+Pp(1)=3170; // vapour pressure of C4, from fig 13.3
+Pp(2)=790; // vapour pressure of C5,from fig 13.3
+Pp(3)=250; // vapour pressure of C6,from fig 13.3
+i=1;
+while(i<4)
+ p(i)=(Pp(i)*x(i));
+ printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
+i=i+1;
+end
+pt=p(1)+p(2)+p(3);
+printf("\t total pressure is : %.1f mm Hg \n",pt);
+printf("\t pressure is too high \n");
+printf("\t for T 96 F \n");
+Pp(1)=2990; // vapour pressure of C4, from fig 13.3
+Pp(2)=725; // vapour pressure of C5,from fig 13.3
+Pp(3)=229; // vapour pressure of C6,from fig 13.3
+i=1;
+while(i<4)
+ p(i)=(Pp(i)*x(i));
+ printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
+i=i+1;
+end
+pt=p(1)+p(2)+p(3);
+printf("\t total pressure is : %.1f mm Hg \n",pt);
+printf("\t pressure is too low \n");
+printf("\t for T 97 F \n");
+Pp(1)=3040; // vapour pressure of C4, from fig 13.3
+Pp(2)=740; // vapour pressure of C5,from fig 13.3
+Pp(3)=234; // vapour pressure of C6,from fig 13.3
+i=1;
+while(i<4)
+ p(i)=(Pp(i)*x(i));
+ printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
+i=i+1;
+end
+pt=p(1)+p(2)+p(3);
+printf("\t total pressure is : %.1f mm Hg \n",pt);
+i=1;
+while(i<4)
+ y(i)=(Pp(i)*x(i)/pt);
+ printf("\n x(i) y(i) \n "+string(x(i))+" "+string(y(i))+" \n");
+ i=i+1;
+end
+printf("\t solution for b \n");
+// Similarly at what temperature will the mixture start to boil if the system is under a pressure of 35 psia
+printf("\t for T 150 F \n");
+Pp(1)=6100; // vapour pressure of C4, from fig 13.3
+Pp(2)=1880; // vapour pressure of C5,from fig 13.3
+Pp(3)=680; // vapour pressure of C6,from fig 13.3
+i=1;
+while(i<4)
+ p(i)=(Pp(i)*x(i));
+ printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
+i=i+1;
+end
+pt=p(1)+p(2)+p(3);
+printf("\t total pressure is : %.0f mm Hg \n",pt);
+printf("\t pressure is too high \n");
+printf("\t for T 149F \n");
+Pp(1)=6050; // vapour pressure of C4, from fig 13.3
+Pp(2)=1850; // vapour pressure of C5,from fig 13.3
+Pp(3)=670; // vapour pressure of C6,from fig 13.3
+i=1;
+while(i<4)
+ p(i)=(Pp(i)*x(i));
+ printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
+i=i+1;
+end
+pt=p(1)+p(2)+p(3);
+printf("\t total pressure is : %.0f mm Hg \n",pt);
+i=1;
+while(i<4)
+ y(i)=(Pp(i)*x(i)/pt);
+ printf("\n x(i) y(i) \n "+string(x(i))+" "+string(y(i))+" \n");
+ i=i+1;
+end
+printf("\t solution for c \n");
+printf("\t for T 95F \n");
+K(1)=3.13; // fig 7
+K(2)=0.92; // fig 7
+K(3)=0.30; // fig 7
+i=1;
+while(i<4)
+ y(i)=(K(i)*x(i));
+ printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
+ i=i+1;
+end
+yt=y(1)+y(2)+y(3);
+printf("\t yt is : %.3f \n",yt);
+printf("\t yt is too low \n");
+printf("\t for T 100F \n");
+K(1)=3.35; // fig 7
+K(2)=1; // fig 7
+K(3)=0.335; // fig 7
+i=1;
+while(i<4)
+ y(i)=(K(i)*x(i));
+ printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
+ i=i+1;
+end
+yt=y(1)+y(2)+y(3);
+printf("\t yt is : %.3f \n",yt);
+printf("\t yt is too low \n");
+printf("\t for T 102F \n");
+K(1)=3.45; // fig 7
+K(2)=1.02; // fig 7
+K(3)=0.35; // fig 7
+i=1;
+while(i<4)
+ y(i)=(K(i)*x(i));
+ printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
+ i=i+1;
+end
+yt=y(1)+y(2)+y(3);
+printf("\t yt is : %.3f \n",yt);
+printf("\t solution for d \n");
+// The use of K values gives y, directly and permits use of the total mol fraction of yt = 1.00 as the criterion for equilibrium
+printf("\t for T 150F \n");
+K(1)=2.8; // fig 7
+K(2)=1.01; // fig 7
+K(3)=0.4; // fig 7
+i=1;
+while(i<4)
+ y(i)=(K(i)*x(i));
+ printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
+ i=i+1;
+end
+yt=y(1)+y(2)+y(3);
+printf("\t yt is : %.3f \n",yt);
+printf("\t yt is too low \n");
+printf("\t for T 153F \n");
+K(1)=2.90; // fig 7
+K(2)=1.06; // fig 7
+K(3)=0.415; // fig 7
+i=1;
+while(i<4)
+ y(i)=(K(i)*x(i));
+ printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
+ i=i+1;
+end
+yt=y(1)+y(2)+y(3);
+printf("\t yt is : %.3f \n",yt);
+printf("\t solution for e at pt=760mm Hg \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
+printf("\t for T 130F \n");
+K(1)=5; // fig 7
+K(2)=1.65; // fig 7
+K(3)=0.62; // fig 7
+i=1;
+while(i<4)
+ x(i)=(y(i)/K(i));
+ printf("\n y(i) K(i) x(i) \n "+string(y(i))+" "+string(K(i))+" "+string(x(i))+" \n");
+ i=i+1;
+end
+xt=x(1)+x(2)+x(3);
+printf("\t xt is : %.3f \n",xt);
+printf("\t xt is too low \n");
+printf("\t for T 120F \n");
+K(1)=4.4; // fig 7
+K(2)=1.4; // fig 7
+K(3)=0.51; // fig 7
+i=1;
+while(i<4)
+ x(i)=(y(i)/K(i));
+ printf("\n y(i) K(i) x(i) \n "+string(y(i))+" "+string(K(i))+" "+string(x(i))+" \n");
+ i=i+1;
+end
+xt=x(1)+x(2)+x(3);
+printf("\t xt is : %.3f \n",xt);
+printf("\t xt is high \n");
+printf("\t for T 123F \n");
+K(1)=4.6; // fig 7
+K(2)=1.49; // fig 7
+K(3)=0.545; // fig 7
+i=1;
+while(i<4)
+ x(i)=(y(i)/K(i));
+ printf("\n y(i) K(i) x(i) \n "+string(y(i))+" "+string(K(i))+" "+string(x(i))+" \n");
+ i=i+1;
+end
+xt=x(1)+x(2)+x(3);
+printf("\t xt is : %.3f \n",xt);
+printf("\t dew point at 760mm is 123F \n");
+printf("\t dew point at 35psia \n");
+printf("\t for T 174F \n");
+K(1)=3.7; // fig 7
+K(2)=1.38; // fig 7
+K(3)=0.58; // fig 7
+i=1;
+while(i<4)
+ x(i)=(y(i)/K(i));
+ printf("\n y(i) K(i) x(i) \n "+string(y(i))+" "+string(K(i))+" "+string(x(i))+" \n");
+ i=i+1;
+end
+xt=x(1)+x(2)+x(3);
+printf("\t xt is : %.3f \n",xt);
+printf("\t dew point is 174F \n");
+// end
generated by cgit (git 2.34.1) at 2024-12-19 11:13:03 +0000