initial commit / add all books

1 files changed, 50 insertions, 0 deletions

diff --git a/611/CH12/EX12.5/Chap12_Ex5_R1.sce b/611/CH12/EX12.5/Chap12_Ex5_R1.sce
new file mode 100755
index 000000000..bbcd878e4
--- /dev/null
+++ b/611/CH12/EX12.5/Chap12_Ex5_R1.sce
@@ -0,0 +1,50 @@
+// Y.V.C.Rao ,1997.Chemical Engineering Thermodynamics.Universities Press,Hyderabad,India.

+

+//Chapter-12,Example 5,Page 435

+//Title: VLE data using the van Laar model

+//================================================================================================================

+clear 

+clc

+

+//INPUT

+//For convenience chloroform is taken as 1 and methanol is taken as 2

+P=760;//pressure in Torr at which chloroform and methanol form an azeotrope

+T=53.5;//temperature in degree celsius at which chloroform and methanol form an azeotrope

+x1=0.65;//mole fraction of chloroform in the liquid phase (no unit) (corresponding to azeotropic composition)

+antoine_const_chloroform=[6.95465;1170.966;226.232];//Antoine's constants for acetone from Table A.7

+antoine_const_methanol=[8.08097;1582.271;239.726];//Antoine's constants for acetone from Table A.7

+

+//CALCULATION

+//The form of the Antoine's equation used is logP=A-(B/(t+C)), where P is in Torr and t is in degree celsius

+x2=1-x1;//calculation of the mole fraction of methanol in the liquid phase (no unit) (corresponding to azeotropic composition)

+//calculation of saturation pressure of chloroform at T in Torr

+P1_s=10^(antoine_const_chloroform(1,:)-(antoine_const_chloroform(2,:)/(T+antoine_const_chloroform(3,:))));

+//calculation of saturation pressure of methanol at T in Torr

+P2_s=10^(antoine_const_methanol(1,:)-(antoine_const_methanol(2,:)/(T+antoine_const_methanol(3,:))));

+//At the azeotropic conditions, the activity coefficients are determined using Eq.(12.15 and 12.16)

+gaamma1=P/P1_s;//calculation of activity coefficient using Eq.(12.15) (no unit)

+gaamma2=P/P2_s;//calculation of activity coefficient using Eq.(12.16) (no unit)

+A=log(gaamma1)*(1+((x2*log(gaamma2))/(x1*log(gaamma1))))^2;//calculation of the van Laar parameter (no unit) using Eq.(11.84)

+B=log(gaamma2)*(1+((x1*log(gaamma1))/(x2*log(gaamma2))))^2;//calculation of the van Laar parameter (no unit) using Eq.(11.85)

+x1=0.1:0.1:0.9;//taking the values of mole fraction of chloroform in the liquid phase to compute the VLE data (no unit)

+l=length(x1);//iteration parameter

+j=1;//iteration parameter

+while j<l|j==l

+    ln_gaamma1(j)=A/(1+((A*x1(:,j))/(B*(1-x1(:,j)))))^2;//calculation of ln(activity coefficient) (no unit)

+    ln_gaamma2(j)=B/(1+((B*(1-x1(:,j)))/(A*x1(:,j))))^2;//calculation of ln(activity coefficient) (no unit)

+    gaamma1(j)=exp(ln_gaamma1(j));//calculation of the activity coefficient (no unit)

+    gaamma2(j)=exp(ln_gaamma2(j));//calculation of the activity coefficient (no unit)

+    P(j)=(gaamma1(j)*x1(:,j)*P1_s)+(gaamma2(j)*(1-x1(:,j))*P2_s);//calculation of pressure in Torr

+    y1(j)=(gaamma1(j)*x1(:,j)*P1_s)/P(j);//calculation of mole fraction of chloroform in vapour phase (no unit)

+    j=j+1;

+end

+

+//OUTPUT

+mprintf('VLE data: \n\n');

+i=1;

+mprintf('x1 \tgamma1 \t\t gamma2  \t P (Torr) \t y1 \n\n');

+for i=1:l

+    mprintf('%0.1f \t %f \t %f \t %f \t %f \n',x1(i),gaamma1(i),gaamma2(i),P(i),y1(i));

+end

+

+//===============================================END OF PROGRAM===================================================