summaryrefslogtreecommitdiff
path: root/1922/CH9
diff options
context:
space:
mode:
Diffstat (limited to '1922/CH9')
-rwxr-xr-x1922/CH9/EX9.1.a/9_1a.sce11
-rwxr-xr-x1922/CH9/EX9.1.b/9_1b.sce23
-rwxr-xr-x1922/CH9/EX9.10.a/9_10a.sce14
-rwxr-xr-x1922/CH9/EX9.10.b/9_10b.sce41
-rwxr-xr-x1922/CH9/EX9.2/9_2.sce16
-rwxr-xr-x1922/CH9/EX9.3/9_3.sce21
-rwxr-xr-x1922/CH9/EX9.4.a/9_4.sce14
-rwxr-xr-x1922/CH9/EX9.4.b/9_4b.sce19
-rwxr-xr-x1922/CH9/EX9.5/9_5.sce10
-rwxr-xr-x1922/CH9/EX9.6.a/9_6a.sce14
-rwxr-xr-x1922/CH9/EX9.6.b/9_6b.sce15
-rwxr-xr-x1922/CH9/EX9.6.c/9_6c.sce15
-rwxr-xr-x1922/CH9/EX9.7/9_7.sce19
-rwxr-xr-x1922/CH9/EX9.9/9_9.sce16
14 files changed, 248 insertions, 0 deletions
diff --git a/1922/CH9/EX9.1.a/9_1a.sce b/1922/CH9/EX9.1.a/9_1a.sce
new file mode 100755
index 000000000..8d022d351
--- /dev/null
+++ b/1922/CH9/EX9.1.a/9_1a.sce
@@ -0,0 +1,11 @@
+clc
+clear
+//Initialization of variables
+g11=178900 //kJ/kmol
+g12=207037 //kJ/kmol
+g21=211852 //kJ/kmol
+g22=228097 //kJ/kmol
+//calculations
+dG=g21-g11
+//results
+printf("Standard free energy change = %d kJ/kmol",dG)
diff --git a/1922/CH9/EX9.1.b/9_1b.sce b/1922/CH9/EX9.1.b/9_1b.sce
new file mode 100755
index 000000000..c82c60efc
--- /dev/null
+++ b/1922/CH9/EX9.1.b/9_1b.sce
@@ -0,0 +1,23 @@
+clc
+clear
+//Initialization of variables
+m1=54.1
+m2=56.1
+m3=2
+cp1=2.122 //kJ/kmol K
+cp2=2.213 //kJ/kmol K
+cp3=14.499 //kJ/kmol K
+hf1=110200 //kJ/kmol
+hf2=-126 //kJ/kmol
+T=700 //K
+Ts=298 //K
+//calculations
+hf=hf1-hf2
+cpn=cp1*m1-cp2*m2+cp3*m3
+h700=hf+ cpn*(T-Ts)
+s298=103.7
+s700 = s298 + cpn*log(T/Ts)
+G700=h700-T*s700
+//results
+printf("Change in gibbs energy = %d kJ/kmol",G700)
+disp("The answer is a bit different due to rounding off error in textbook")
diff --git a/1922/CH9/EX9.10.a/9_10a.sce b/1922/CH9/EX9.10.a/9_10a.sce
new file mode 100755
index 000000000..c528052cf
--- /dev/null
+++ b/1922/CH9/EX9.10.a/9_10a.sce
@@ -0,0 +1,14 @@
+clc
+clear
+//Initialization of variables
+kp=1.09
+feed=[ 1 5 0 0 0 ]
+//calculations
+x=poly(0,"x")
+vec=roots(kp/4^4 /4 *(1-x)*(5-2*x)^2 *(6+2*x)^2 -x^5)
+x=vec(5)
+pro=[1-x 5-2*x x 4*x 0]
+//results
+disp("Equlibrium composition (moles)= ")
+format('v',6);pro
+disp(pro)
diff --git a/1922/CH9/EX9.10.b/9_10b.sce b/1922/CH9/EX9.10.b/9_10b.sce
new file mode 100755
index 000000000..d27c6bd4b
--- /dev/null
+++ b/1922/CH9/EX9.10.b/9_10b.sce
@@ -0,0 +1,41 @@
+clc
+clear
+//Initialization of variables
+kp=1.09
+kp2=0.154
+feed=[ 1 5 0 0 0 ]
+//calculations
+x=poly(0,"x")
+vec=roots(kp/4^4 /4 *(1-x)*(5-2*x)^2 *(6+2*x)^2 -x^5)
+x=vec(5)
+pro=feed- [x 2*x -x -4*x 0]
+y=poly(0,"y")
+vec2=roots(kp2*(0.273-y)*(0.727-y)*(7.454+2*y)^2 - 4*y^2 *(2.908+2*y)^2 *4)
+y=vec2(4)
+pro2=pro-[ y 0 y -2*y -2*y]
+z=poly(0,"z")
+vec3= roots(kp*(0.189-z)*(3.546-2*z)^2 *(7.622+2*z)^2 -(0.643+z)*(3.076+4*z)^4 *4)
+z=vec3(5)
+pro3=pro2 - [z 2*z -z -4*z 0]
+w=poly(0,"w")
+vec4=roots(kp2*(0.229-w)*(0.603-w)*(7.542+2*w) - (2.916+2*w)^2 *(0.168+2*w)^2 *4)
+w=vec4(4)
+w=0.01
+pro4=pro3-[w 0 w -2*w -2*w]
+//results
+disp("feed = ")
+format('v',6);feed
+disp(feed)
+disp("After reactor 1,")
+format('v',6);pro
+disp(pro)
+disp("After reactor 2,")
+format('v',6);pro2
+disp(pro2)
+disp("After reactor 3,")
+format('v',6);pro3
+disp(pro3)
+disp("After reactor 4")
+format('v',6);pro4
+disp(pro4)
+disp("The answers are a bit different due to rounding off error in textbook")
diff --git a/1922/CH9/EX9.2/9_2.sce b/1922/CH9/EX9.2/9_2.sce
new file mode 100755
index 000000000..057d564cf
--- /dev/null
+++ b/1922/CH9/EX9.2/9_2.sce
@@ -0,0 +1,16 @@
+clc
+clear
+//Initialization of variables
+g1=150670 //kJ/kmol
+g2=71500 //kJ/kmol
+R=8.314
+Ts=298 //K
+T=700 //K
+//calculationd
+G=g1-g2
+G2=33875 //kJ/kmol
+K1=exp(-G/R/Ts)
+K2=exp(-G2/R/T)
+//results
+printf("In case 1, equilibrium constant = %.2e",K1)
+printf("\n In case 2, equilibrium constant = %.5f",K2)
diff --git a/1922/CH9/EX9.3/9_3.sce b/1922/CH9/EX9.3/9_3.sce
new file mode 100755
index 000000000..bbc19f1ce
--- /dev/null
+++ b/1922/CH9/EX9.3/9_3.sce
@@ -0,0 +1,21 @@
+clc
+clear
+//Initialization of variables
+R=8.3143
+T1=1273 //K
+T2=2273 //K
+k2=0.0018
+A=123.94
+B=7.554
+C=8.552*10^-3
+D=-13.25e-6
+E=7.002e-9
+F=13.494e-13
+//calculations
+function y= cp(T)
+ y=A/T^2 +B/T +C +D*T +E*T^2 -F*T^3
+endfunction
+lnk=1/R *intg(T1,T2,cp)
+k1=k2/ exp(lnk)
+//results
+printf("Equilibrium constant = %.5f ",k1)
diff --git a/1922/CH9/EX9.4.a/9_4.sce b/1922/CH9/EX9.4.a/9_4.sce
new file mode 100755
index 000000000..129fef3bd
--- /dev/null
+++ b/1922/CH9/EX9.4.a/9_4.sce
@@ -0,0 +1,14 @@
+clc
+clear
+//Initialization of variables
+G=-30050 //kJ/kmol
+R=8.314
+T=573 //K
+//calculations
+lnk=G/(R*T)
+k=exp(lnk)
+x=poly(0,"x")
+vec=roots(4*x^2 - k*(1-x)^2)
+x2=vec(2)
+//results
+printf("Mole fraction of HCN = %.4f",x2)
diff --git a/1922/CH9/EX9.4.b/9_4b.sce b/1922/CH9/EX9.4.b/9_4b.sce
new file mode 100755
index 000000000..8aea25314
--- /dev/null
+++ b/1922/CH9/EX9.4.b/9_4b.sce
@@ -0,0 +1,19 @@
+clc
+clear
+//Initialization of variables
+G=-30050 //kJ/kmol
+R=8.314
+T=573 //K
+phi1=0.980
+phi2=0.915
+phi3=0.555
+//calculations
+lnk=G/(R*T)
+k=exp(lnk)
+kexp= k*phi1*phi2/phi3^2 /4
+x=poly(0,"x")
+vec=roots(x^2 - kexp*(1-x)^2)
+x2=vec(2)
+//results
+printf("Mole fraction of HCN = %.4f",x2)
+
diff --git a/1922/CH9/EX9.5/9_5.sce b/1922/CH9/EX9.5/9_5.sce
new file mode 100755
index 000000000..82ca6ae7a
--- /dev/null
+++ b/1922/CH9/EX9.5/9_5.sce
@@ -0,0 +1,10 @@
+clc
+clear
+//Initialization of variables
+kp=74
+//calculations
+f=poly(0,"f")
+vec=roots(f^2 *(100-6*f) - kp^2 *(1-f)^2 *(9-6*f))
+fn=vec(3)
+//results
+printf("Fractional conversion = %.3f",fn)
diff --git a/1922/CH9/EX9.6.a/9_6a.sce b/1922/CH9/EX9.6.a/9_6a.sce
new file mode 100755
index 000000000..113e6cc48
--- /dev/null
+++ b/1922/CH9/EX9.6.a/9_6a.sce
@@ -0,0 +1,14 @@
+clc
+clear
+//Initialization of variables
+C=3
+phi=3
+R=1
+Sc=0
+function V=fun(C,phi,R,Sc)
+ V=2+C-phi-R-Sc
+endfunction
+//calculations
+V=fun(C,phi,R,Sc)
+//results
+printf("Degrees of freedom = %d ",V)
diff --git a/1922/CH9/EX9.6.b/9_6b.sce b/1922/CH9/EX9.6.b/9_6b.sce
new file mode 100755
index 000000000..568da78e8
--- /dev/null
+++ b/1922/CH9/EX9.6.b/9_6b.sce
@@ -0,0 +1,15 @@
+clc
+clear
+//Initialization of variables
+C=3
+phi=1
+R=1
+Sc=1
+function V=fun(C,phi,R,Sc)
+ V=2+C-phi-R-Sc
+endfunction
+//calculations
+V=fun(C,phi,R,Sc)
+//results
+printf("Degrees of freedom = %d ",V)
+
diff --git a/1922/CH9/EX9.6.c/9_6c.sce b/1922/CH9/EX9.6.c/9_6c.sce
new file mode 100755
index 000000000..5b8d4b5ae
--- /dev/null
+++ b/1922/CH9/EX9.6.c/9_6c.sce
@@ -0,0 +1,15 @@
+clc
+clear
+//Initialization of variables
+C=6
+phi=1
+R=3
+Sc=0
+function V=fun(C,phi,R,Sc)
+ V=2+C-phi-R-Sc
+endfunction
+//calculations
+V=fun(C,phi,R,Sc)
+//results
+printf("Degrees of freedom = %d ",V)
+
diff --git a/1922/CH9/EX9.7/9_7.sce b/1922/CH9/EX9.7/9_7.sce
new file mode 100755
index 000000000..1e5caaa9e
--- /dev/null
+++ b/1922/CH9/EX9.7/9_7.sce
@@ -0,0 +1,19 @@
+clc
+clear
+//Initialization of variables
+a1=0.956
+y=0.014
+x=0.956
+M=18
+z=0.475
+P=8.37 //Mpa
+//calculations
+m=y/(x*M) *10^3
+w=0.0856
+phi1=-0.04
+phi2=0.06
+phi=10^(phi1+ w*phi2)
+f=z*phi*P
+K=m/(f*a1)
+//results
+printf("Equilibrium constant = %.3f",K)
diff --git a/1922/CH9/EX9.9/9_9.sce b/1922/CH9/EX9.9/9_9.sce
new file mode 100755
index 000000000..dae50b85e
--- /dev/null
+++ b/1922/CH9/EX9.9/9_9.sce
@@ -0,0 +1,16 @@
+clc
+clear
+//Initialization of variables
+y=0.18
+z=0.6
+//calculations
+mole=[1-y-z 5-y-2*z y 3*y+4*z z]
+s=sum( mole)
+molef=mole/s
+//results
+disp("Product composition moles = ")
+format('v',6);mole
+disp(mole)
+disp("Mole fraction = ")
+format('v',6);molef
+disp(molef)
generated by cgit (git 2.34.1) at 2025-03-13 10:16:51 +0000