initial commit / add all books

10 files changed, 207 insertions, 0 deletions

diff --git a/1442/CH19/EX19.1/19_1.sce b/1442/CH19/EX19.1/19_1.sce
new file mode 100755
index 000000000..1502b390e
--- /dev/null
+++ b/1442/CH19/EX19.1/19_1.sce
@@ -0,0 +1,22 @@
+clc

+//initialisation of variables

+pN2= 79 //percent

+VN2= 82.3 //m^3

+VCO2= 8 //m^3

+VCO= 0.9 //m^3

+M= 32 //gms

+M1= 28 //gms

+//CALCULATIONS

+P= (pN2/(100-pN2))

+z= VN2/P

+x= VCO2+VCO

+w= VCO2+(VCO/2)+(VCO2/10)

+y= 2*w

+r= y/x

+TO= x+(y/4)

+X= (z/TO)-1

+AF= z*(M+P*M1)/(12*x+y)

+//RESULTS

+printf (' fuel ratio= %.3f  ',r)

+printf (' \n excess air= %.3f  ',X)

+printf (' \n emass air-fuel ratio= %.2f  ',AF)

diff --git a/1442/CH19/EX19.10/19_10.sce b/1442/CH19/EX19.10/19_10.sce
new file mode 100755
index 000000000..5f89c75dd
--- /dev/null
+++ b/1442/CH19/EX19.10/19_10.sce
@@ -0,0 +1,25 @@
+clc

+//initialisation of variables

+T= 25 //C

+R= 8.314 //Jmol K

+k= 1.27

+k1= 1.34

+hf= -393520 //kJ/kmol

+M= 28 //gms

+gf= -394360 //kJ/kmol

+M= 12 //gms

+//CALCULATIONS

+T1= T+(-hf/((R)*((k/(k-1))+(0.2+4.5144)*(k1/(k1-1)))))

+Bin= 0

+dh= (k1*R/(k1-1))*(T1-T)

+dh1= (k1*R/(k1-1))*log((273.15+T1)/(273.15+T))

+H= dh-(273.15+T)*dh1

+h= (k*R/(k-1))*(T1-T)+hf

+h1= (k*R/(k-1))*log((273.15+T1)/(273.15+T))+((hf-gf)/(273.15+T))

+h2= h-(273.15+T)*h1

+Bout= (h2+(0.2+4.5144)*H)/M

+//RESULTS

+printf (' outlet temperature= %.2f C ',T1)

+printf (' \n energy of formation= %.f J ',Bin)

+printf (' \n energy at outlet= %.f kJ/kmol ',H)

+printf (' \n energy of the products= %.f k ',Bout)

diff --git a/1442/CH19/EX19.11/19_11.sce b/1442/CH19/EX19.11/19_11.sce
new file mode 100755
index 000000000..d0dce13d6
--- /dev/null
+++ b/1442/CH19/EX19.11/19_11.sce
@@ -0,0 +1,43 @@
+clc

+//initialisation of variables

+b= 1475.30 //kJ/kg

+b0= 144.44 //kJ/kg

+h2= 3448.6 //kJkg

+h1= 860.5 //kJ/kg

+k= 1.27 

+k1= 1.34

+R= 8.314 //J/mol K

+hf= -393520 //kJ/kmol

+hg= 72596 //kJ/kmol

+Mc= 12 //kg

+n= 1.2 //moles

+n1= 3.76 //moles

+M= 32 //gms

+M1= 28 //gms

+M2= 44 //gms

+n2= 0.2 //moles

+n3= 4.512 //moles

+B1= 25592 //kJ/kmol C

+B2= 394360 //kJ/kmol C

+e= 0.008065

+//CALCULATIONS

+B= b-b0

+Q= h2-h1

+CpCO2= k*R/(k-1)

+CpO2= k1*R/(k1-1)

+Qcoal= (hg+hf)/Mc

+mcoal= Q/(-Qcoal)

+ncoal= mcoal/Mc

+r= (n*M+n1*M1)/Mc

+r1= (M2+n2*M+n3*M1)/Mc

+mair= r*mcoal

+mgas= r1*mcoal

+Bfuel= (B1-B2)*e

+Bnet= Bfuel+B

+p= B*100/(-Bfuel)

+//RESULTS

+printf (' change in energy= %.2f kJ/kg ',B)

+printf (' \n amount of air= %.3f kg/kg ',mair)

+printf (' \n amount of gas= %.3f kg/kg ',mgas)

+printf (' \n net change in energy= %.2f kg/kg steam ',Bnet)

+printf (' \n percent energy in original fuel= %.2f percent ',p)

diff --git a/1442/CH19/EX19.2/19_2.sce b/1442/CH19/EX19.2/19_2.sce
new file mode 100755
index 000000000..3dc62b6ec
--- /dev/null
+++ b/1442/CH19/EX19.2/19_2.sce
@@ -0,0 +1,28 @@
+clc

+//initialisation of variables

+m1= 24 //kg

+M1= 32 //kg

+m2= 28 //kg

+M2= 28 //kg

+e= 0.5

+T3= 1800 //C

+T0= 25 //C

+T1= 25 //C

+T2= 100 //C

+R= 8.314 //Jmol K

+cp= 4.57 //J/mol K

+cp1= 3.5 //J/mol K

+cp2= 3.5 //J/mol K

+hCO2= -393522 //J

+hCO= -110529 //J

+//CALCULATIONS

+n1= m1/M1

+n2= m2/M2

+N= n1-0.5*e

+N1= n2-e

+N2= e

+N3= N+N1+N2

+y1= N/N3

+Q= ((N*cp+N1*cp1+N2*cp2)*R*(T3-T0)-(n1*cp*(T1-T0)+n2*cp2*(T2-T1))+N*(hCO2-hCO))/60

+//RESULTS

+printf (' Heat interaction= %.f kW ',Q)

diff --git a/1442/CH19/EX19.3/19_3.sce b/1442/CH19/EX19.3/19_3.sce
new file mode 100755
index 000000000..7e2472763
--- /dev/null
+++ b/1442/CH19/EX19.3/19_3.sce
@@ -0,0 +1,13 @@
+clc

+//initialisation of variables

+hCO2= -393520 //kJ/kg mol

+hH2O= -285840 //kJ/kg mol

+hC7H16= -187820 //kJ/kg mol

+M= 100

+hH2O1= -241830 //kJkg mol

+//CALCULATIONS

+HHV= -(7*hCO2+8*hH2O-hC7H16)/M

+LLV= -(7*hCO2+8*hH2O1-hC7H16)/M

+//RESULTS

+printf (' Higher heating vlue= %.f kJ/kg mol ',HHV)

+printf (' \n Lower heating vlue= %.f kJ/kg mol ',LLV)

diff --git a/1442/CH19/EX19.4/19_4.sce b/1442/CH19/EX19.4/19_4.sce
new file mode 100755
index 000000000..fdfbdf29b
--- /dev/null
+++ b/1442/CH19/EX19.4/19_4.sce
@@ -0,0 +1,14 @@
+clc

+//initialisation of variables

+T0= 25 //C

+T1= 220 //C

+hCO2= -393520 //kJ/kg

+hH2O= -241830 //kJ/kg

+hC3H8= -103850 //kJ/kg= 1.4

+R= 8.314 //Jmol K

+k= 1.4

+k1= 1.29

+//CALCULATIONS

+T= T0+((15*(R*(k/(k-1)))*4.762*(T1-T0)-(3*hCO2+4*hH2O-hC3H8))/(R*((3+4)*(k1/(k1-1))+(10+56.43)*(k/(k-1)))))

+//RESULTS

+printf (' adiabatic flame temperature= %.f C ',T)

diff --git a/1442/CH19/EX19.6/19_6.sce b/1442/CH19/EX19.6/19_6.sce
new file mode 100755
index 000000000..9769b578c
--- /dev/null
+++ b/1442/CH19/EX19.6/19_6.sce
@@ -0,0 +1,14 @@
+clc

+//initialisation of variables

+T= 25 //C

+hfT= -241820 //kJ/kmol

+R= 8.314 //J/mol K

+k= 1.4

+cpH2O= 4.45

+cpO2= 3.5

+T1= 1000 //C

+//CALCULATIONS

+S= (cpH2O-k*cpO2)

+hfT1= hfT+S*(T1-T)

+//RESULTS

+printf (' enthalpy formation= %.f kJ/kmol ',hfT1)

diff --git a/1442/CH19/EX19.7/19_7.sce b/1442/CH19/EX19.7/19_7.sce
new file mode 100755
index 000000000..74e370684
--- /dev/null
+++ b/1442/CH19/EX19.7/19_7.sce
@@ -0,0 +1,15 @@
+clc

+//initialisation of variables

+R= 8.314 //J/mol K

+T= 25 //C

+gf= 16590 //kJ/kmol

+T1= 500 //C

+Cp= 4.157 //J/mol K

+hf= -46190 //kJ/kmol

+//CALCULATIONS

+K= %e^(gf/(R*(273.15+T)))

+r= (1-((273.15+T)/(273.15+T1)))*((hf/(R*(273.15+T)))+(R/Cp))-2*log((273.15+T1)/(273.15+T))+0.6

+KT1= K*%e^r

+//RESULTS

+printf (' equilibrium constant= %.1f bar^-1 ',K)

+printf (' \n equilibrium constant= %.5f bar^-1 ',KT1)

diff --git a/1442/CH19/EX19.8/19_8.sce b/1442/CH19/EX19.8/19_8.sce
new file mode 100755
index 000000000..5751cfd3e
--- /dev/null
+++ b/1442/CH19/EX19.8/19_8.sce
@@ -0,0 +1,24 @@
+clc

+//initialisation of variables

+uCO2= -394374 //J/mol

+uCO= -137150 //J/mol

+uO2= 0

+R= 8.314 //J/mol K

+T= 25 //C

+cpCO2= 4.57 //J/mol K

+cpCO= 3.5 //J/mol K

+cpO2= 3.5 //J/mol K

+T1= 1500 //C

+hf= -393522 //kJ/kmol

+gf= -110529 //kJ/kmol

+T2= 2500 //C

+//CALCULATIONS

+r= -(uCO2-uCO-0.5*uO2)/(R*(273.15+T))

+s= (cpCO2-cpCO-0.5*cpO2)

+r1= (1-((273.15+T)/(273.15+T1)))*((hf-gf)/(R*(273.15+T))-s)+s*log((273.15+T1)/(273.15+T))

+KT1= %e^(r+r1)

+r2= (1-((273.15+T)/(273.15+T2)))*((hf-gf)/(R*(273.15+T))-s)+s*log((273.15+T2)/(273.15+T))

+KT2= %e^(r+r2)

+//RESULTS

+printf (' equilibrium constant at T1= %.3f C ',KT1)

+printf (' \n equilibrium constant at T2= %.3f C ',KT2)

diff --git a/1442/CH19/EX19.9/19_9.sce b/1442/CH19/EX19.9/19_9.sce
new file mode 100755
index 000000000..80fd6fe80
--- /dev/null
+++ b/1442/CH19/EX19.9/19_9.sce
@@ -0,0 +1,9 @@
+clc

+//initialisation of variables

+Wc= 12 //kg

+hf= -393520 //kJ/kmol

+gf= -394360 //kJ/kmol

+//CALCULATIONS

+Wmax= -gf/Wc

+//RESULTS

+printf (' maximum work= %.f kJ/kg of carbon ',Wmax)