3 files changed, 67 insertions, 0 deletions
diff --git a/1394/CH12/EX12.2.1/Ex12_2_1.sce b/1394/CH12/EX12.2.1/Ex12_2_1.sce
new file mode 100755
index 000000000..54bab483b
--- /dev/null
+++ b/1394/CH12/EX12.2.1/Ex12_2_1.sce
@@ -0,0 +1,15 @@
+
+clc
+//initialization of variables
+l = 1.22 // length of tower
+Gflow = 0.026 // mol/sec
+GbyL = 0.07
+dia = 0.088 // m
+pl = 1.1/100// pl = 1-yl
+p0 = 0.04/100 // p0 = 1-y0
+//Calculations
+A = %pi*(dia^2)/4 // cross sectional of tower in m^2
+G = Gflow/A // Gas flux in mol/m^2-sec
+Kya = (G/l)*(1/(1-GbyL))*(log(pl/p0))// Mass transfer per volume in mol/m^3-sec
+//Results
+printf("The mass transfer per volume is %.1f mol/m^3-sec",Kya)
diff --git a/1394/CH12/EX12.2.2/Ex12_2_2.sce b/1394/CH12/EX12.2.2/Ex12_2_2.sce
new file mode 100755
index 000000000..10b317dbc
--- /dev/null
+++ b/1394/CH12/EX12.2.2/Ex12_2_2.sce
@@ -0,0 +1,27 @@
+
+clc
+//Initialization of variables
+x1=0.99
+x2=0.99
+y1=0.95
+y2=0.95
+alpha=1.5
+m=0.42
+l=2
+HTU=0.34
+//calculations
+y1s= (y1-0.58)/m
+xrd= (x2-y2)/(x1-y1s)
+Rd=xrd/(1-xrd)
+Rds=alpha*Rd
+xl= ((Rds+1)*y1 - x1)/(Rds)
+function z= ystar(y)
+    z=0.58+0.42*y
+endfunction
+NTU = log((ystar(xl) -y1)/(ystar(x1)-x1)) /(1- m*(Rds+1)/Rds)
+NTU2=l/HTU
+xd2=(ystar(y1)-y1)/%e^(NTU2*(1-m))
+xd=(0.58-xd2)/(1-m)
+//results
+printf("In case 1, NTU = %.2f",NTU)
+printf("\n In case 2, xd = %.3f",xd)
diff --git a/1394/CH12/EX12.4.1/Ex12_4_1.sce b/1394/CH12/EX12.4.1/Ex12_4_1.sce
new file mode 100755
index 000000000..fccc568ce
--- /dev/null
+++ b/1394/CH12/EX12.4.1/Ex12_4_1.sce
@@ -0,0 +1,25 @@
+
+clc
+//initialization of variables
+F=3500 //mol/hr
+xf=0.4
+x1=0.98
+y1=0.97
+y2=0.625
+x1=0.97
+x2=0.4
+ratio=1.5
+HTU=0.2
+//calculations
+A=[1 1;x1 1-x1]
+B=[F;xf*F]
+C=A\B
+DA=C(1)
+BA=C(2)
+Rds=(y1-y2)/(x1-x2)
+Rd=Rds/(1-Rds)
+Rdreq=ratio*Rd
+NTU=13.9
+l=HTU*NTU
+//results
+printf("length of the tower = %.1f m",l)