From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 3516/CH13/EX13.1/Ex13_1.sce | 207 ++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 207 insertions(+)
 create mode 100644 3516/CH13/EX13.1/Ex13_1.sce

(limited to '3516/CH13/EX13.1/Ex13_1.sce')

diff --git a/3516/CH13/EX13.1/Ex13_1.sce b/3516/CH13/EX13.1/Ex13_1.sce
new file mode 100644
index 000000000..18e1019e8
--- /dev/null
+++ b/3516/CH13/EX13.1/Ex13_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
-- 
cgit