1 files changed, 60 insertions, 0 deletions
diff --git a/2090/CH3/EX3.15/Chapter3_Example15.sce b/2090/CH3/EX3.15/Chapter3_Example15.sce
new file mode 100755
index 000000000..6b1117fca
--- /dev/null
+++ b/2090/CH3/EX3.15/Chapter3_Example15.sce
@@ -0,0 +1,60 @@
+clc
+clear
+//Input data
+T=600;//The initial temperature of air in K
+p=1;//The initial pressure of air in atm
+R=8.314;//Universal gas constant in J/molK
+Tr=298;//The temperature of reactants in K
+a=4.503;//Given Constant 
+b=-8.965*10^-3;//Given Constant
+c=37.38*10^-6;//Given Constant
+d=-36.49*10^-9;//Given Constant
+e=12.22*10^-12;//Given Constant
+
+//Calculations
+hfc=-393.52;//Enthalpy of formation for carbondioxide in MJ/kmol fuel
+hfh=-241.82;//Enthalpy of formation for water in MJ/kmol fuel
+hfn=0;//Enthalpy of formation for nitrogen gas
+HP=[1*hfc]+[2*hfh]+[7.52*hfn];//The enthalpy on the products side in MJ/kmol fuel
+hch=[R*[(a*(T-Tr))+((b/2)*(T^2-Tr^2))+((c/3)*(T^3-Tr^3))+((d/4)*(T^4-Tr^4))+((e/5)*(T^5-Tr^5))]]/1000;//The change in enthalpy of the methane in MJ/kmol
+hfc1=-74.87;//The enthalpy of formation for methane in MJ/kmol fuel 
+hfh1=9.247;//The change in enthalpy of the water in MJ/kmol
+hfn1=8.891;//The change in enthalpy of nitrogen in MJ/kmol
+HR=[(hfc1+hch)+(2*hfh1)+(7.52*hfn1)];//The enthalpy on the reactants side in MJ/kmol
+x=HR-HP;//The enthalpy for the remaining gases in the product side in MJ/kmol
+hfn2=x/7.52;//The guess enthalpy for the nitrogen gas in MJ/kmol
+Tc=3700;//The corresponding temperature for the enthalpy of guess nitrogen in K
+T1=2800;//The temperature assumed for the first trail in K
+hco1=140.444;//The change in enthalpy for the assume temp for carbondioxide in MJ/kmol
+hh1=115.294;//The change in enthalpy for the assume temp for water in MJ/kmol
+hn1=85.345;//The change in enthalpy for the assume temp for nitrogen in MJ/kmol
+HP1=hco1+(2*hh1)+(7.52*hn1)+(HR-x);//The total enthalpy on the products side for first trail in MJ/kmol fuel
+T2=2500;//The temperature assumed for the second trail in K
+hco2=121.926;//The change in enthalpy for the assume temp for carbondioxide in MJ/kmol
+hh2=98.964;//The change in enthalpy for the assume temp for water in MJ/kmol
+hn2=74.312;//;//The change in enthalpy for the assume temp for nitrogen in MJ/kmol
+HP2=hco2+(2*hh2)+(7.52*hn2)+(HR-x);//The total enthalpy on the products side for the second trail in MJ/kmol
+T3=2600;//The temperature fo the third trail in K
+hco3=128.085;//The change in enthalpy for the assume temp for carbondioxide in MJ/kmol
+hh3=104.37;//The change in enthalpy for the assume temp for water in MJ/kmol
+hn3=77.973;//The change in enthalpy for the assume temp for nitrogen in MJ/kmol
+HP3=hco3+(2*hh3)+(7.52*hn3)+(HR-x);//The total enthalpy on the products side for the third trail in MJ/kmol
+Ta1=[[(HR-HP2)/(HP3-HP2)]*(T3-T2)]+T2;//The adiabatic temperature for constant pressure process in K
+UR1=HR-(10.52*R*10^-3*T);//The internal energy of reactant in MJ/kmol fuel
+Tc1=3000;//Assume temperature for first trail in K
+hcoa1=146.645;//The change in enthalpy for the assume temp for carbondioxide in MJ/kmol
+hha1=120.813;//The change in enthalpy for the assume temp for carbondioxide in MJ/kmol
+hna1=89.036;//The change in enthalpy for the assume temp for nitrogen in MJ/kmol
+UP1=hcoa1+(2*hha1)+(7.52*hna1)+(HR-x)-(0.08746*Tc1);//The internal energy of products in MJ/kmol fuel
+Tc2=3200;//Assume temperature for the second trail in K
+hcoa2=165.331;//;//The change in enthalpy for the assume temp for carbondioxide in MJ/kmol
+hha2=137.553;//The change in enthalpy for the assume temp for water in MJ/kmol
+hna2=100.161;//The change in enthalpy for the assume temp for nitrogen in MJ/kmol
+UP2=hcoa2+(2*hha2)+(7.52*hna2)+(HR-x)-(0.08746*Tc2);//The internal energy of products in MJ/kmol fuel
+Tu=[[(UR1-UP1)/(UP2-UP1)]*(Tc2-Tc1)]+Tc1;//The adiabatic flame temperature at constant pressure process in K
+
+//Output
+printf('The adiabatic flame temperature at  \n    (a)Constant pressure process is %3.0f K \n    (b)Constant volume process is %3.1f K',Ta1,Tu)
+
+
+