initial commit / add all books

8 files changed, 125 insertions, 0 deletions

diff --git a/1709/CH11/EX11.1/11_1.sce b/1709/CH11/EX11.1/11_1.sce
new file mode 100755
index 000000000..6a708fe94
--- /dev/null
+++ b/1709/CH11/EX11.1/11_1.sce
@@ -0,0 +1,28 @@
+clc

+//Initialization of variables

+x=1.5

+P=14.696 //psia

+m=28.96

+//calculations

+mf=114 // lbm/mol fuel

+ma=x*12.5*(1+3.76)*m

+AF=ma/mf

+n1=8

+n2=9

+n3=(x-1)*12.5 

+n4= x*3.76*12.5

+np=n1+n2+n3+n4

+x1=n1/np

+x2=n2/np

+x3=n3/np

+x4=n4/np

+ph=x2*P

+Td=113.5 //F

+//results

+printf("Air fuel ratio = %.1f lbm air/lbm fuel",AF)

+printf("\n Mole fraction of CO2 = %.2f percent",x1)

+printf("\n Mole fraction of H2O = %.2f percent",x2)

+printf("\n Mole fraction of O2 = %.2f percent",x3)

+printf("\n Mole fraction of N2 = %.2f percent",x4)

+disp("From tables of saturation pressure")

+printf("Dew point  = %.1f F",Td)

diff --git a/1709/CH11/EX11.2/11_2.sce b/1709/CH11/EX11.2/11_2.sce
new file mode 100755
index 000000000..e2185acfd
--- /dev/null
+++ b/1709/CH11/EX11.2/11_2.sce
@@ -0,0 +1,15 @@
+clc

+//Initialization of variables

+x1=9

+x2=1.2

+x3=1.5

+x4=88.3

+//calculations

+a=x1+x2

+b=2*a

+xO=(2*x1 + x2+ 2*x3 + b)/2

+xN=x4/3.76

+ratio=xO/a

+percent=ratio/2 *100

+//results

+printf("Percent theoretical air = %.1f percent",percent)

diff --git a/1709/CH11/EX11.3/11_3.sce b/1709/CH11/EX11.3/11_3.sce
new file mode 100755
index 000000000..249fcc810
--- /dev/null
+++ b/1709/CH11/EX11.3/11_3.sce
@@ -0,0 +1,11 @@
+clc

+//Initialization of variables

+T=440 //F

+//calculations

+disp("From steam tables,")

+h1=-169290

+h2=7597.6

+h3=4030.2

+ht=h1+h2-h3

+//results

+printf("Molal enthalpy of CO2 = %d Btu/lbm mole",ht)

diff --git a/1709/CH11/EX11.4/11_4.sce b/1709/CH11/EX11.4/11_4.sce
new file mode 100755
index 000000000..7bd462a14
--- /dev/null
+++ b/1709/CH11/EX11.4/11_4.sce
@@ -0,0 +1,11 @@
+clc

+//Initialization of variables

+T=77 //F

+//calculations

+Hr=-36420 //B

+hc=-169290 //B/lb mol

+hh=-122970 //B/lb mol

+Hp=2*hc+3*hh

+Q=Hp-Hr

+//results

+printf("Heat transfer = %d B/mol fuel",Q)

diff --git a/1709/CH11/EX11.5/11_5.sce b/1709/CH11/EX11.5/11_5.sce
new file mode 100755
index 000000000..ca8eb84bb
--- /dev/null
+++ b/1709/CH11/EX11.5/11_5.sce
@@ -0,0 +1,20 @@
+clc

+//Initialization of variables

+T2=440 //F

+T1=77 //F

+Mch4=16

+Mw=18

+//calculations

+h77=3725.1

+ht=6337.9

+ht2=7597.6

+h772=4030.2

+hwt=1260.3

+h77w=45.02

+hr77=-383040 //B/lbm mol

+dHR=1*Mch4*0.532*(T1-T2) + 2*(h77-ht)

+dHp=1*(ht2-h772) + 2*Mw*(hwt - h77w)

+hrp=dHp+hr77+dHR

+//results

+printf("Enthalpy of combustion of gaseous methane = %d B/lbm mol fuel",hrp)

+//The calculation in textbook is wrong Please check it using a calculator.

diff --git a/1709/CH11/EX11.6/11_6.sce b/1709/CH11/EX11.6/11_6.sce
new file mode 100755
index 000000000..366cdd795
--- /dev/null
+++ b/1709/CH11/EX11.6/11_6.sce
@@ -0,0 +1,7 @@
+clc

+//Initialization of variables

+Hr=-107530 //B/mol fuel

+disp("By iteration of temperatures, T=2700 R")

+T=2700 //R

+//results

+printf("Adiabatic flame temperature = %d R",T)

diff --git a/1709/CH11/EX11.7/11_7.sce b/1709/CH11/EX11.7/11_7.sce
new file mode 100755
index 000000000..83ef3c77d
--- /dev/null
+++ b/1709/CH11/EX11.7/11_7.sce
@@ -0,0 +1,17 @@
+clc

+//Initialization of variables

+Kp=0.668

+y=Kp^2

+//calculations

+x=poly(0,"x")

+vec=roots(x^3 + y*x^3 + 2*y*x^2 -y*x -2*y)

+eps=vec(1)

+x1=(1-eps)/(1+ eps/2)

+x2=eps/(1+eps/2)

+x3=eps/2/(1+ eps/2)

+//results

+printf("degree of reaction = %.3f ",eps)

+printf("\n Equilibrium concentration of CO2 = %.3f ",x1)

+printf("\n Equilibrium concentration of CO = %.3f ",x2)

+printf("\n Equilibrium concentration of O2 = %.3f ",x3)

+//the answers are a bit different due to approximation in textbook

diff --git a/1709/CH11/EX11.8/11_8.sce b/1709/CH11/EX11.8/11_8.sce
new file mode 100755
index 000000000..39fe1e0e3
--- /dev/null
+++ b/1709/CH11/EX11.8/11_8.sce
@@ -0,0 +1,16 @@
+clc

+//Initialization of variables

+Kp=15.63

+y=Kp

+//calculations

+x=poly(0,"x")

+vec=roots(x^2 + y*x^2 - y)

+eps=vec(1)

+x1=(1-eps)/(1+eps)

+x2=eps/(1+eps)

+x3=eps/(1+eps)

+//results

+printf(" Equilibrium concentration of Cs = %.4f ",x1)

+printf("\n Equilibrium concentration of Cs+ = %.4f ",x2)

+printf("\n Equilibrium concentration of e- = %.4f ",x3)

+//the answers are a bit different due to approximation in textbook