summaryrefslogtreecommitdiff
path: root/Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter8_1.ipynb
diff options
context:
space:
mode:
authorTrupti Kini2016-03-08 23:30:20 +0600
committerTrupti Kini2016-03-08 23:30:20 +0600
commitdf4be1d705b512d28c07ff063ba4045423e019fb (patch)
tree3bf559595e10d22a6c494eaa61884c72059600a0 /Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter8_1.ipynb
parent50b3a7d796c572462e5a5e6a79fee66414f20a1d (diff)
downloadPython-Textbook-Companions-df4be1d705b512d28c07ff063ba4045423e019fb.tar.gz
Python-Textbook-Companions-df4be1d705b512d28c07ff063ba4045423e019fb.tar.bz2
Python-Textbook-Companions-df4be1d705b512d28c07ff063ba4045423e019fb.zip
Added(A)/Deleted(D) following books
A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter10_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter11_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter12_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter13_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter14_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter15_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter16_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter17_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter18_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter19_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter1_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter2_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter4_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter5_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter7_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter8_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter9_1.ipynb A Basic_Engineering_Thermodynamics_by_Rayner_Joel/screenshots/chap11_1.png A Basic_Engineering_Thermodynamics_by_Rayner_Joel/screenshots/chap14_1.png A Basic_Engineering_Thermodynamics_by_Rayner_Joel/screenshots/chap7_1.png A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch1.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch11.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch12.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch15.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch16.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch3.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch5.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch6.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch7.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch8.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/Ch9.ipynb A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/screenshots/euclideanAlgo11.png A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/screenshots/gcd11.png A Discrete_Mathematics_by_S._Lipschutz,_M._Lipson_And_V._H._Patil/screenshots/primeFactor11.png A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter10_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter11_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter1_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter2_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter3_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter4_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter5_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter6_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter8_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/Chapter9_1.ipynb A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/screenshots/chap1_1.png A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/screenshots/chap2_1.png A Introduction_to_Electrical_Engineering_by_Er._J.P._Navani_&_Er._Sonal_Sapra/screenshots/chap3_1.png A sample_notebooks/UmangAgarwal/Sample_Notebook_Umang_1.ipynb
Diffstat (limited to 'Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter8_1.ipynb')
-rw-r--r--Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter8_1.ipynb1584
1 files changed, 1584 insertions, 0 deletions
diff --git a/Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter8_1.ipynb b/Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter8_1.ipynb
new file mode 100644
index 00000000..9280416c
--- /dev/null
+++ b/Basic_Engineering_Thermodynamics_by_Rayner_Joel/Chapter8_1.ipynb
@@ -0,0 +1,1584 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 8 - Combustion"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 1: pg 198"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.1\n",
+ " The stoichiometric mass of air required to burn 1 kg the fuel should be (kg/kg fuel) = 14.8\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 198\n",
+ "print('Example 8.1');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the stoichiometric mass of air required to burn 1 kg the fuel \n",
+ "\n",
+ "# Given values\n",
+ "C = .72;# mass fraction of C; [kg/kg]\n",
+ "H2 = .20;# mass fraction of H2;, [kg/kg]\n",
+ "O2 = .08;# mass fraction of O2, [kg/kg]\n",
+ "aO2=.232;# composition of oxygen in air\n",
+ "\n",
+ "# solution\n",
+ "# for 1kg of fuel\n",
+ "mO2 = 8./3*C+8*H2-O2;# mass of O2, [kg]\n",
+ "\n",
+ "# hence stoichiometric mass of O2 required is\n",
+ "msO2 = mO2/aO2;# [kg]\n",
+ "\n",
+ "#results\n",
+ "print ' The stoichiometric mass of air required to burn 1 kg the fuel should be (kg/kg fuel) = ',round(msO2,1)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 2: pg 198"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.2\n",
+ " The stoichiometric mass of air required to burn 1 kg the fuel should be (kg/kg fuel) = 15.18\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 198\n",
+ "print('Example 8.2');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the stoichiometric mass of air required to burn 1 kg of heptane\n",
+ "\n",
+ "# Given values\n",
+ "C = .84;# mass fraction of C; [kg/kg]\n",
+ "H2 = .16;# mass fraction of H2;, [kg/kg]\n",
+ "x1 = 11.5;# mass fraction of O2, [kg/kg]\n",
+ "x2 = 34.5;# composition of oxygen in air\n",
+ "\n",
+ "# solution\n",
+ "# for 1kg of fuel\n",
+ "mO2 = x1*C + x2*H2 ;# mass of O2, [kg]\n",
+ "\n",
+ "mO22 = ((11*32)+100)/100\n",
+ "mO2x = (mO22-1)/.232\n",
+ "#results\n",
+ "print ' The stoichiometric mass of air required to burn 1 kg the fuel should be (kg/kg fuel) = ',round(mO2,2)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 3: pg 199"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.3\n",
+ " The stoichiometric mass of air is (kg/kg fuel) = 13.52\n",
+ " CO2 produced = 3.01 kg/kg fuel, percentage composition = 20.7 ,\n",
+ " H2O produced = 1.08 kg/kg fuel, percentage composition = 7.43 ,\n",
+ " SO2 produced = 0.02 kg/kg fuel, percentage composition = 0.14 ,\n",
+ " N2 produced = 10.43 kg/kg fuel, percentage composition = 71.74\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 199\n",
+ "print('Example 8.3');\n",
+ "\n",
+ "# aim : To determine \n",
+ "# the stoichiometric mass of air\n",
+ "# the products of combustion both by mass and as percentage \n",
+ "\n",
+ "# Given values\n",
+ "C = .82;# mass composition C\n",
+ "H2 = .12;# mass composition of H2\n",
+ "O2 = .02;# mass composition of O2\n",
+ "S = .01;# mass composition of S\n",
+ "N2 = .03;# mass composition of N2\n",
+ "\n",
+ "# solution\n",
+ "# for 1kg fuel\n",
+ "mo2 = 8./3*C+8*H2-O2+S*1;# total mass of O2 required, [kg]\n",
+ "sa = mo2/.232;# stoichimetric air, [kg]\n",
+ "print ' The stoichiometric mass of air is (kg/kg fuel) = ',round(sa,2)\n",
+ "\n",
+ "# for one kg fuel\n",
+ "mCO2 = C*11/3;# mass of CO2 produced, [kg]\n",
+ "mH2O = H2*9;# mass of H2O produced, [kg]\n",
+ "mSO2 = S*2;# mass of SO2 produce, [kg]\n",
+ "mN2 = C*8.84+H2*26.5-O2*.768/.232+S*3.3+N2;# mass of N2 produced, [kg]\n",
+ "\n",
+ "mt = mCO2+mH2O+mSO2+mN2;# total mass of product, [kg]\n",
+ "\n",
+ "x1 = mCO2/mt*100;# %age mass composition of CO2 produced\n",
+ "x2 = mH2O/mt*100;# %age mass composition of H2O produced\n",
+ "x3 = mSO2/mt*100;# %age mass composition of SO2 produced\n",
+ "x4 = mN2/mt*100;# %age mass composition of N2 produced\n",
+ "\n",
+ "print ' CO2 produced = ',round(mCO2,2),' kg/kg fuel, percentage composition = ',round(x1,1),',\\n H2O produced = ',mH2O,' kg/kg fuel, percentage composition = ',round(x2,2),',\\n SO2 produced = ',mSO2,' kg/kg fuel, percentage composition = ',round(x3,2),',\\n N2 produced = ',round(mN2,2),' kg/kg fuel, percentage composition = ',round(x4,2)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 4: pg 202"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.4\n",
+ " The stoichiometric volume of air required for complete combustion is (m^3/m^3 fuel) = 0.952\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 202\n",
+ "print('Example 8.4');\n",
+ "\n",
+ "# aim : To determine \n",
+ "# the stoichiometric volume of air required for complete combution of 1 m^3 of the gas\n",
+ "\n",
+ "# Given values\n",
+ "H2 = .14;# volume fraction of H2\n",
+ "CH4 = .02;# volume fraction of CH4\n",
+ "CO = .22;# volume fraction of CO\n",
+ "CO2 = .05;# volume fraction of CO2\n",
+ "O2 = .02;# volume fraction of O2\n",
+ "N2 = .55;# volume fraction of N2\n",
+ "\n",
+ "# solution\n",
+ "# for 1 m^3 of fuel\n",
+ "Va = .5*H2+2*CH4+.5*CO-O2;# [m^3]\n",
+ "\n",
+ "# stoichiometric air required is\n",
+ "Vsa = Va/.21;# [m^3]\n",
+ "#results\n",
+ "print ' The stoichiometric volume of air required for complete combustion is (m^3/m^3 fuel) = ',round(Vsa,3)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 5: pg 203"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.5\n",
+ " The volume of air required is(m^3/m^3 fuel) = 3.451\n",
+ "Result in the book is misprinted\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 203\n",
+ "print('Example 8.5');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the volume of the air required \n",
+ "\n",
+ "# Given values\n",
+ "H2 = .45;# volume fraction of H2\n",
+ "CO = .40;# volume fraction of CO\n",
+ "CH4 = .15;# volume fraction of CH4\n",
+ "\n",
+ "# solution\n",
+ "V = 2.38*(H2+CO)+9.52*CH4;# stoichimetric volume of air, [m^3]\n",
+ "#results\n",
+ "print ' The volume of air required is(m^3/m^3 fuel) = ',V\n",
+ "\n",
+ "print 'Result in the book is misprinted'\n",
+ "\n",
+ "# End\n",
+ "\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 6: pg 203"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 6,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.6\n",
+ " Stoichiometric volume of air required is (m^3/m^3 fuel) = 3.15\n",
+ " N2 in products = 2.516 m^3/m^3 fuel, percentage composition = 66.8 ,\n",
+ " CO2 in products = 0.561 m^3/m^3 fuel, percentage composition = 14.9 ,\n",
+ " H2O in products = 0.689 m^3/m^3 fuel, percentage composition = 18.3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 203\n",
+ "print('Example 8.6');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the stoichiometric volume of air for the complete combustion\n",
+ "# the products of combustion\n",
+ "\n",
+ "# given values\n",
+ "CH4 = .142;# volumetric composition of CH4\n",
+ "CO2 = .059;# volumetric composition of CO2\n",
+ "CO = .360;# volumetric composition of CO\n",
+ "H2 = .405;# volumetric composition of H2\n",
+ "O2 = .005;# volumetric composition of O2\n",
+ "N2 = .029;# volumetric composition of N2\n",
+ "\n",
+ "aO2 = .21;# O2 composition into air by volume\n",
+ "\n",
+ "# solution\n",
+ "svO2 = CH4*2+CO*.5+H2*.5-O2;# stroichiometric volume of O2 required, [m^3/m^3 fuel]\n",
+ "svair = svO2/aO2;# stroichiometric volume of air required, [m^3/m^3 fuel]\n",
+ "print ' Stoichiometric volume of air required is (m^3/m^3 fuel) = ',svair\n",
+ "\n",
+ "# for one m^3 fuel\n",
+ "vN2 = CH4*7.52+CO*1.88+H2*1.88-O2*.79/.21+N2;# volume of N2 produced, [m^3]\n",
+ "vCO2 = CH4*1+CO2+CO*1;# volume of CO2 produced, [m^3]\n",
+ "vH2O = CH4*2+H2*1;# volume of H2O produced, [m^3]\n",
+ "\n",
+ "vt = vN2+vCO2+vH2O;# total volume of product, [m^3]\n",
+ "\n",
+ "x1 = vN2/vt*100;# %age composition of N2 in product,\n",
+ "x2 = vCO2/vt*100;# %age composition of CO2 in product\n",
+ "x3 = vH2O/vt*100;# %age composition of H2O in product\n",
+ "\n",
+ "print ' N2 in products = ',round(vN2,3),' m^3/m^3 fuel, percentage composition = ',round(x1,1),',\\n CO2 in products = ',vCO2,' m^3/m^3 fuel, percentage composition =',round(x2,1),',\\n H2O in products = ',vH2O,'m^3/m^3 fuel, percentage composition =',round(x3,1)\n",
+ "\n",
+ "# End \n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7: pg 206"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.7\n",
+ " Mass percentage of CO2 is = 27.8 \n",
+ "\n",
+ " Mass percentage of N2 is = 62.0 \n",
+ "\n",
+ " Mass percentage of O2 is = 10.1\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 206\n",
+ "print('Example 8.7');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the percentage analysis of the gas by mass\n",
+ "\n",
+ "# Given values\n",
+ "CO2 = 20.;# percentage volumetric composition of CO2\n",
+ "N2 = 70.;# percentage volumetric composition of N2\n",
+ "O2 = 10.;# percentage volumetric composition of O2\n",
+ "\n",
+ "mCO2 = 44.;# moleculer mas of CO2\n",
+ "mN2 = 28.;# moleculer mass of N2\n",
+ "mO2 = 32.;# moleculer mass of O2\n",
+ "\n",
+ "# solution\n",
+ "mgas = CO2*mCO2+N2*mN2+O2*mO2;# moleculer mass of gas \n",
+ "m1 = CO2*mCO2/mgas*100;# percentage composition of CO2 by mass \n",
+ "m2 = N2*mN2/mgas*100;# percentage composition of N2 by mass \n",
+ "m3 = O2*mO2/mgas*100;# percentage composition of O2 by mass \n",
+ "#results\n",
+ "print ' Mass percentage of CO2 is = ',round(m1,1),' \\n\\n Mass percentage of N2 is = ',round(m2,1),' \\n\\n Mass percentage of O2 is = ',round(m3,1)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 8: pg 206"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 14,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.8\n",
+ " The percentage composition of CO by volume is = 6.9 \n",
+ ",\n",
+ "The percentage composition of N2 by volume is = 4.6 \n",
+ "\n",
+ "The percentage composition of CH4 by volume is = 6.0 \n",
+ "\n",
+ "The percentage composition of H2 by volume is = 80.5 \n",
+ "\n",
+ "The percentage composition of O2by volume is= 2.0\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 206\n",
+ "print('Example 8.8');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the percentage composition of the gas by volume\n",
+ "\n",
+ "# given values\n",
+ "CO = 30.;# %age mass composition of CO\n",
+ "N2 = 20.;# %age mass composition of N2\n",
+ "CH4 = 15.;# %age mass composition of CH4\n",
+ "H2 = 25.;# %age mass composition of H2\n",
+ "O2 = 10.;# %age mass composition of O2\n",
+ "\n",
+ "mCO = 28.;# molculer mass of CO\n",
+ "mN2 = 28.;# molculer mass of N2\n",
+ "mCH4 = 16.;# molculer mass of CH4\n",
+ "mH2 = 2.;# molculer mass of H2\n",
+ "mO2 = 32.;# molculer mass of O2\n",
+ "\n",
+ "# solution\n",
+ "vg = CO/mCO+N2/mN2+CH4/mCH4+H2/mH2+O2/mO2;\n",
+ "v1 = CO/mCO/vg*100;# %age volume composition of CO\n",
+ "v2 = N2/mN2/vg*100;# %age volume composition of N2\n",
+ "v3 = CH4/mCH4/vg*100;# %age volume composition of CH4\n",
+ "v4 = H2/mH2/vg*100;# %age volume composition of H2\n",
+ "v5 = O2/mO2/vg*100;# %age volume composition of O2\n",
+ "#results\n",
+ "print ' The percentage composition of CO by volume is = ',round(v1,1),' \\n,\\nThe percentage composition of N2 by volume is =',round(v2,1), '\\n\\nThe percentage composition of CH4 by volume is = ',round(v3,1),'\\n\\nThe percentage composition of H2 by volume is = ',round(v4,1),' \\n\\nThe percentage composition of O2by volume is= ',round(v5,1)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {
+ "collapsed": true
+ },
+ "source": [
+ "## Example 9: pg 209"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 15,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.9\n",
+ " The mass of air supplied per/kg of fuel burnt is (kg) = 21.1\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 209\n",
+ "print('Example 8.9');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the mass of air supplied per kilogram of fuel burnt\n",
+ "\n",
+ "# given values\n",
+ "CO2 = 8.85;# volume composition of CO2\n",
+ "CO = 1.2;# volume composition of CO\n",
+ "O2 = 6.8;# volume composition of O2\n",
+ "N2 = 83.15;# volume composition of N2 \n",
+ "\n",
+ "# composition of gases in the fuel oil\n",
+ "C = .84;# mass composition of carbon \n",
+ "H = .14;# mass composition of hydrogen\n",
+ "o2 = .02;# mass composition of oxygen\n",
+ "\n",
+ "mC = 12.;# moleculer mass of Carbon\n",
+ "mCO2 = 44.;# molculer mass of CO2\n",
+ "mCO = 28.;# molculer mass of CO\n",
+ "mN2 = 28.;# molculer mass of N2\n",
+ "mO2 = 32.;# molculer mass of O2\n",
+ "aO2 = .23;# mass composition of O2 in air\n",
+ "\n",
+ "# solution\n",
+ "ma = (8./3*C+8*H-o2)/aO2;# theoretical mass of air/kg fuel, [kg]\n",
+ "\n",
+ "mgas = CO2*mCO2+CO*mCO+N2*mN2+O2*mO2;# total mass of gas/kg fuel, [kg]\n",
+ "x1 = CO2*mCO2/mgas;# composition of CO2 by mass \n",
+ "x2 = CO*mCO/mgas;# composition of CO by mass\n",
+ "x3 = O2*mO2/mgas;# composition of O2 by mass \n",
+ "x4 = N2*mN2/mgas;# composition of N2 by mass \n",
+ "\n",
+ "m1 = x1*mC/mCO2+x2*mC/mCO;# mass of C/kg of dry flue gas, [kg]\n",
+ "m2 = C;# mass of C/kg fuel, [kg]\n",
+ "mf = m2/m1;# mass of dry flue gas/kg fuel, [kg]\n",
+ "mo2 = mf*x3;# mass of excess O2/kg fuel, [kg]\n",
+ "mair = mo2/aO2;# mass of excess air/kg fuel, [kg]\n",
+ "m = ma+mair;# mass of excess air supplied/kg fuel, [kg]\n",
+ "#results\n",
+ "print ' The mass of air supplied per/kg of fuel burnt is (kg) = ',round(m,1)\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 10: pg 210"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 16,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.10\n",
+ "The percentage composition of CO2 by volume is = 10.98 \n",
+ ",\n",
+ "The percentage composition of H2O by volume is = 10.72 \n",
+ ",\n",
+ "The percentage composition of O2 by volume is = 3.27 \n",
+ ",\n",
+ "The percentage composition of N2 by volume is = 75.03\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 210\n",
+ "print('Example 8.10');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# volumetric composition of the products of combustion\n",
+ "\n",
+ "# given values\n",
+ "C = .86;# mass composition of carbon\n",
+ "H = .14;# mass composition of hydrogen\n",
+ "Ea = .20;# excess air for combustion\n",
+ "O2 = .23;# mass composition of O2 in air \n",
+ "\n",
+ "MCO2 = 44.;# moleculer mass of CO2\n",
+ "MH2O = 18.;# moleculer mass of H2O\n",
+ "MO2 = 32.;# moleculer mass of O2\n",
+ "MN2 = 28.;# moleculer mass of N2,\n",
+ "\n",
+ "\n",
+ "# solution\n",
+ "sO2 = (8./3*C+8*H);# stoichiometric O2 required, [kg/kg petrol]\n",
+ "sair = sO2/O2;# stoichiometric air required, [kg/kg petrol]\n",
+ "# for one kg petrol\n",
+ "mCO2 = 11./3*C;# mass of CO2,[kg]\n",
+ "mH2O = 9*H;# mass of H2O, [kg]\n",
+ "mO2 = Ea*sO2;# mass of O2, [kg]\n",
+ "mN2 = 14.84*(1+Ea)*(1-O2);# mass of N2, [kg]\n",
+ "\n",
+ "mt = mCO2+mH2O+mO2+mN2;# total mass, [kg]\n",
+ "# percentage mass composition\n",
+ "x1 = mCO2/mt*100;# mass composition of CO2\n",
+ "x2 = mH2O/mt*100;# mass composition of H2O\n",
+ "x3 = mO2/mt*100;# mass composition of O2\n",
+ "x4 = mN2/mt*100;# mass composition of N2\n",
+ "\n",
+ "vt = x1/MCO2+x2/MH2O+x3/MO2+x4/MN2;# total volume of petrol\n",
+ "v1 = x1/MCO2/vt*100;# %age composition of CO2 by volume\n",
+ "v2 = x2/MH2O/vt*100;# %age composition of H2O by volume\n",
+ "v3 = x3/MO2/vt*100;# %age composition of O2 by volume\n",
+ "v4 = x4/MN2/vt*100;# %age composition of N2 by volume\n",
+ " #results\n",
+ "print 'The percentage composition of CO2 by volume is =',round(v1,2),'\\n,\\nThe percentage composition of H2O by volume is = ',round(v2,2),' \\n,\\nThe percentage composition of O2 by volume is = ',round(v3,2),'\\n,\\nThe percentage composition of N2 by volume is = ',round(v4,2)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 11: pg 211"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 17,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.11\n",
+ " The energy carried away by the dry flue gas/kg is (kg) = 5000.0\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 211\n",
+ "print('Example 8.11');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the energy carried away by the dry flue gas/kg of fuel burned\n",
+ "\n",
+ "# given values\n",
+ "C = .78;# mass composition of carbon\n",
+ "H2 = .06;# mass composition of hydrogen\n",
+ "O2 = .09;# mass composition of oxygen\n",
+ "Ash = .07;# mass composition of ash\n",
+ "Ea = .50;# excess air for combustion\n",
+ "aO2 = .23;# mass composition of O2 in air \n",
+ "Tb = 273.+20;# boiler house temperature, [K]\n",
+ "Tf = 273.+320;# flue gas temperature, [K]\n",
+ "c = 1.006;# specific heat capacity of dry flue gas, [kJ/kg K]\n",
+ "\n",
+ "# solution\n",
+ "# for one kg of fuel\n",
+ "sO2 = (8./3*C+8*H2);# stoichiometric O2 required, [kg/kg fuel]\n",
+ "sO2a = sO2-O2;# stoichiometric O2 required from air, [kg/kg fuel]\n",
+ "sair = sO2a/aO2;# stoichiometric air required, [kg/kg fuel]\n",
+ "ma = sair*(1+Ea);# actual air supplied/kg of fuel, [kg]\n",
+ "# total mass of flue gas/kg fuel is\n",
+ "mf = ma+1;# [kg]\n",
+ "mH2 = 9*H2;#H2 produced, [kg] \n",
+ "# hence, mass of dry flue gas/kg coall is\n",
+ "m = mf-mH2;# [kg]\n",
+ "Q = m*c*(Tf-Tb);# energy carried away by flue gas, [kJ]\n",
+ "#results\n",
+ "print ' The energy carried away by the dry flue gas/kg is (kg) = ',round(Q)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 12: pg 212"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 18,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.12\n",
+ " (a) The stoichiometric volume of air for the complete combustion (m^3/m^3 gas) = 9.702\n",
+ " (b) The percentage volumetric composition of CO2 in produced is = 9.6 \n",
+ ",\n",
+ " The percentage volumetric composition of H2O in produced is = 18.8 \n",
+ ",\n",
+ " The percentage volumetric composition of N2 in produced is = 71.6\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 212\n",
+ "print('Example 8.12');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# (a) the stoichiometric volume of air for the complete combustion of 1 m^3\n",
+ "# (b) the percentage volumetric analysis of the products of combustion\n",
+ "\n",
+ "# given values\n",
+ "N2 = .018;# volumetric composition of N2\n",
+ "CH4 = .94;# volumetric composition of CH4\n",
+ "C2H6 = .035;# volumetric composition of C2H6\n",
+ "C3H8 = .007;# volumetric composition of C3H8\n",
+ "aO2 = .21;# O2 composition in air\n",
+ "\n",
+ "# solution\n",
+ "# (a)\n",
+ "# for CH4\n",
+ "# CH4 +2 O2= CO2 + 2 H2O\n",
+ "sva1 = 2./aO2;# stoichiometric volume of air, [m^3/m^3 CH4]\n",
+ "svn1 = sva1*(1-aO2);# stoichiometric volume of nitrogen in the air, [m^3/m^3 CH4]\n",
+ "\n",
+ "# for C2H6\n",
+ "# 2 C2H6 +7 O2= 4 CO2 + 6 H2O\n",
+ "sva2 = 7./2/aO2;# stoichiometric volume of air, [m^3/m^3 C2H6]\n",
+ "svn2 = sva2*(1-aO2);# stoichiometric volume of nitrogen in the air, [m^3/m^3 C2H6]\n",
+ "\n",
+ "# for C3H8\n",
+ "# C3H8 +5 O2=3 CO2 + 4 H2O\n",
+ "sva3 = 5/aO2;# stoichiometric volume of air, [m^3/m^3 C3H8]\n",
+ "svn3 = sva3*(1-aO2);# stoichiometric volume of nitrogen in the air, [m^3/m^3 C3H8]\n",
+ "\n",
+ "Sva = CH4*sva1+C2H6*sva2+C3H8*sva3;# stoichiometric volume of air required, [m^3/m^3 gas]\n",
+ "print ' (a) The stoichiometric volume of air for the complete combustion (m^3/m^3 gas) = ',round(Sva,3)\n",
+ "\n",
+ "# (b)\n",
+ "# for one m^3 of natural gas\n",
+ "vCO2 = CH4*1+C2H6*2+C3H8*3;# volume of CO2 produced, [m^3]\n",
+ "vH2O = CH4*2+C2H6*3+C3H8*4;# volume of H2O produced, [m^3]\n",
+ "vN2 = CH4*svn1+C2H6*svn2+C3H8*svn3+N2;# volume of N2 produced, [m^3]\n",
+ "\n",
+ "vg = vCO2+vH2O+vN2;# total volume of gas, [m^3]\n",
+ "x1 = vCO2/vg*100;# volume percentage of CO2 produced\n",
+ "x2 = vH2O/vg*100;# volume percentage of H2O produced\n",
+ "x3 = vN2/vg*100;# volume percentage of N2 produced\n",
+ "\n",
+ "print ' (b) The percentage volumetric composition of CO2 in produced is = ',round(x1,1),' \\n,\\n The percentage volumetric composition of H2O in produced is = ',round(x2,1),' \\n,\\n The percentage volumetric composition of N2 in produced is = ',round(x3,1)\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 13: pg 214"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 19,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.13\n",
+ " (a) Volume of air taken by fan is (m^3/s) = 2.51\n",
+ " (b) Percentage mass composition of CO2 is (percent) = 18.77\n",
+ " Percentage mass composition of O2 is (percent) = 5.24\n",
+ " Percentage mass composition of N2 is (percent) = 75.99\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 214\n",
+ "print('Example 8.13');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# (a) the volume of air taken by the fan\n",
+ "# (b) the percentage composition of dry flue gas\n",
+ "\n",
+ "# gien values\n",
+ "C = .82;# mass composition of carbon\n",
+ "H = .08;# mass composition of hydrogen\n",
+ "O = .03;# mass composition of oxygen\n",
+ "A = .07;# mass composition of ash\n",
+ "mc = .19;# coal uses, [kg/s] \n",
+ "ea = .3;# percentage excess air of oxygen in the air required for combustion\n",
+ "Oa = .23;# percentage of oxygen by mass in the air\n",
+ "\n",
+ "# solution\n",
+ "# (a)\n",
+ "P = 100.;# air pressure, [kN/m^2]\n",
+ "T = 18.+273;# air temperature, [K]\n",
+ "R = .287;# [kJ/kg K]\n",
+ "# basis one kg coal\n",
+ "sO2 = 8./3*C+8*H;# stoichiometric O2 required, [kg]\n",
+ "aO2 = sO2-.03;# actual O2 required, [kg]\n",
+ "tO2 = aO2/Oa;# theoretical O2 required, [kg]\n",
+ "Aa = tO2*(1+ea);# actual air supplied, [kg]\n",
+ "m = Aa*mc;# Air supplied, [kg/s]\n",
+ "\n",
+ "# now using P*V=m*R*T\n",
+ "V = m*R*T/P;# volume of air taken ,[m^3/s]\n",
+ "print ' (a) Volume of air taken by fan is (m^3/s) = ',round(V,2)\n",
+ "\n",
+ "# (b)\n",
+ "mCO2 = 11./3*C;# mass of CO2 produced, [kg]\n",
+ "mO2 = aO2*.3;# mass of O2 produces, [kg]\n",
+ "mN2 = Aa*.77;# mass of N2 produced, [[kg]\n",
+ "mt = mCO2+mO2+mN2;# total mass, [kg]\n",
+ "\n",
+ "print ' (b) Percentage mass composition of CO2 is (percent) = ',round(mCO2/mt*100,2)\n",
+ "print ' Percentage mass composition of O2 is (percent) = ',round(mO2/mt*100,2)\n",
+ "print ' Percentage mass composition of N2 is (percent) = ',round(mN2/mt*100,2)\n",
+ "\n",
+ "\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 14: pg 215"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 21,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.14\n",
+ " (a) Mass of fuel used per cycle is (g) = 0.444\n",
+ " (b) The mass of air supplied per cycle is (kg) = 0.0136\n",
+ " (c) The volume of air taken in per cycle is (m^3) = 0.0114\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 215\n",
+ "print('Example 8.14');\n",
+ "\n",
+ "# aim : To determine \n",
+ "# (a) the mass of fuel used per cycle\n",
+ "# (b) the actual mass of air taken in per cycle\n",
+ "# (c) the volume of air taken in per cycle\n",
+ "\n",
+ "# given values\n",
+ "W = 15.;# work done, [kJ/s]\n",
+ "N = 5.;# speed, [rev/s]\n",
+ "C = .84;# mass composition of carbon\n",
+ "H = .16;# mass composition of hydrogen\n",
+ "ea = 1.;# percentage excess air supplied \n",
+ "CV = 45000.;# calorificvalue of fuel, [kJ/kg]\n",
+ "n_the = .3;# thermal efficiency\n",
+ "P = 100.;# pressuer, [kN/m^2]\n",
+ "T = 273.+15;# temperature, [K]\n",
+ "R = .29;# gas constant, [kJ/kg K]\n",
+ "\n",
+ "# solution\n",
+ "# (a)\n",
+ "E = W*2/N/n_the;# energy supplied, [kJ/cycle]\n",
+ "mf = E/CV;# mass of fuell used, [kg]\n",
+ "print ' (a) Mass of fuel used per cycle is (g) = ',round(mf*10**3,3)\n",
+ "\n",
+ "# (b)\n",
+ "# basis 1 kg fuel\n",
+ "mO2 = C*8./3+8*H;# mass of O2 requirea, [kg]\n",
+ "smO2 = mO2/.23;# stoichiometric mass of air, [kg]\n",
+ "ma = smO2*(1+ea);# actual mass of air supplied, [kg]\n",
+ "m = ma*mf;# mass of air supplied, [kg/cycle]\n",
+ "print ' (b) The mass of air supplied per cycle is (kg) = ',round(m,4)\n",
+ "\n",
+ "# (c)\n",
+ "V = m*R*T/P;# volume of air, [m^3]\n",
+ "print ' (c) The volume of air taken in per cycle is (m^3) = ',round(V,4)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 15: pg 216"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 22,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.15\n",
+ " (a) The mass of coal used per hour is (kg) = 95.7\n",
+ " (b) The mass of air supplied per hour is (kg) = 1590.0\n",
+ " (c) The mass percentage composition of CO2 = 20.83 ,\n",
+ " The mass percentage composition of H2O = 3.08 ,\n",
+ " The mass percentage composition of O2 = 3.89 ,\n",
+ " The mass percentage composition of N2 = 72.2\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 216\n",
+ "print('Example 8.15');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# (a) the mass of coal used per hour\n",
+ "# (b) the mass of air used per hour\n",
+ "# (c) the percentage analysis of the flue gases by mass\n",
+ "\n",
+ "# given values\n",
+ "m = 900.;# mass of steam boiler generate/h, [kg]\n",
+ "x = .96;# steam dryness fraction\n",
+ "P = 1400.;# steam pressure, [kN/m^2]\n",
+ "Tf = 52.;# feed water temperature, [C]\n",
+ "BE = .71;# boiler efficiency\n",
+ "CV = 33000.;# calorific value of coal, [kJkg[\n",
+ "ea = .22;# excess air supply\n",
+ "aO2 = .23;# oxygen composition in air\n",
+ "c = 4.187;# specific heat capacity of water, [kJ/kg K]\n",
+ "\n",
+ "# coal composition\n",
+ "C = .83;# mass composition of carbon\n",
+ "H2 = .05;# mass composition of hydrogen\n",
+ "O2 = .03;# mass composition of oxygen\n",
+ "ash = .09;# mass composition of ash\n",
+ "\n",
+ "# solution\n",
+ "# from steam table at pressure P\n",
+ "hf = 830.1;# specific enthalpy, [kJ/kg]\n",
+ "hfg = 1957.1;# specific enthalpy, [kJ/kg]\n",
+ "hg = 2728.8;# specific enthalpy, [kJ/kg]\n",
+ "\n",
+ "# (a)\n",
+ "h = hf+x*hfg;# specific enthalpy of steam generated by boiler, [kJ/kg]\n",
+ "hfw = c*Tf;# specific enthalpy of feed water, [kJ/kg]\n",
+ "Q = m*(h-hfw);# energy to steam/h, [kJ]\n",
+ "Qf = Q/BE;# energy required from fuel/h, [kJ]\n",
+ "mc = Qf/CV;# mass of coal/h,[kg]\n",
+ "print ' (a) The mass of coal used per hour is (kg) = ',round(mc,1)\n",
+ "\n",
+ "# (b)\n",
+ "# for one kg coal\n",
+ "mO2 = 8./3*C+8*H2-O2;# actual mass of O2 required, [kg]\n",
+ "mta = mO2/aO2;# theoretical mass of air, [kg]\n",
+ "ma = mta*(1+ea);# mass of air supplied, [kg]\n",
+ "mas = ma*116;# mass of air supplied/h, [kg]\n",
+ "print ' (b) The mass of air supplied per hour is (kg) = ',round(mas)\n",
+ "\n",
+ " \n",
+ "# (c)\n",
+ "# for one kg coal\n",
+ "mCO2 = 11./3*C;# mass of CO2 produced, [kg]\n",
+ "mH2O = 9*H2;# mass of H2O produced, [kg]\n",
+ "mO2 = mO2*ea;# mass of excess O2 in flue gas, [kg]\n",
+ "mN2 = ma*(1-aO2);# mass of N2 in flue gas, [kg]\n",
+ "\n",
+ "mt = mCO2+mH2O+mO2+mN2;# total mass of gas\n",
+ "x1 = mCO2/mt*100;# mass percentage composition of CO2\n",
+ "x2 = mH2O/mt*100;# mass percentage composition of H2O\n",
+ "x3 = mO2/mt*100;# mass percentage composition of O2\n",
+ "x4 = mN2/mt*100;# mass percentage composition of N2\n",
+ "\n",
+ "print ' (c) The mass percentage composition of CO2 = ',round(x1,2),' ,\\n The mass percentage composition of H2O = ',round(x2,2),' ,\\n The mass percentage composition of O2 = ',round(x3,2),' ,\\n The mass percentage composition of N2 =',round(x4,2)\n",
+ "# mass of coal taken in part (b) is wrong so answer is not matching\n",
+ "\n",
+ "# End\n",
+ "\n",
+ "\n",
+ "\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 16: pg 223"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 23,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.16\n",
+ " (a) The volume of the gas is (m^3) = 22.4\n",
+ " (b)(1) The average moleculer mass of air is (g/mol) = 28.84\n",
+ " (2) The value of R is (kJ/kg K) = 0.288\n",
+ " (3) The mass of one cubic metre of air at STP is (kg/m^3) = 1.287\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 223\n",
+ "print('Example 8.16');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# (a) volume of gas\n",
+ "# (b) (1) the average molecular mass of air\n",
+ "# (2) the value of R\n",
+ "# (3) the mass of 1 m^3 of air at STP\n",
+ "\n",
+ "# given values\n",
+ "n = 1.;# moles of gas, [kmol]\n",
+ "P = 101.32;# standard pressure, [kN/m^2]\n",
+ "T = 273.;# gas tempearture, [K]\n",
+ "\n",
+ "O2 = 21.;# percentage volume composition of oxygen in air\n",
+ "N2 = 79.;# percentage volume composition of nitrogen in air\n",
+ "R = 8.3143;# molar gas constant, [kJ/kg K]\n",
+ "mO2 = 32.;# moleculer mass of O2\n",
+ "mN2 = 28.;# moleculer mass of N2\n",
+ "\n",
+ "# solution\n",
+ "# (a)\n",
+ "V = n*R*T/P;# volume of gas, [m^3]\n",
+ "print ' (a) The volume of the gas is (m^3) = ',round(V,1)\n",
+ "\n",
+ "# (b)\n",
+ "#(1)\n",
+ "Mav = (O2*mO2+N2*mN2)/(O2+N2);# average moleculer mass of air\n",
+ "print ' (b)(1) The average moleculer mass of air is (g/mol) = ',Mav\n",
+ "\n",
+ "# (2)\n",
+ "Rav = R/Mav;# characteristic gas constant, [kJ/kg k]\n",
+ "print ' (2) The value of R is (kJ/kg K) = ',round(Rav,3)\n",
+ "\n",
+ "# (3)\n",
+ "rho = Mav/V;# density of air, [kg/m^3]\n",
+ "print ' (3) The mass of one cubic metre of air at STP is (kg/m^3) = ',round(rho,3)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 17: pg 223"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 24,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.17\n",
+ " (a)The partial pressure of O2 is (kN/m^2) = 104.0 ,\n",
+ " The partial pressure of N2 is (kN/m^2) = 104.0 \n",
+ " The partial pressure of CO2 is (kN/m^2) = 208.0\n",
+ " (b) The volume of the container is (m^3) = 6.655\n",
+ " (c) The new total pressure in the vessel is (kN/m^2) = 626.0\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 223\n",
+ "print('Example 8.17');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# (a) the partial pressure of each gas in the vessel\n",
+ "# (b) the volume of the vessel\n",
+ "# (c) the total pressure in the gas when temperature is raised to228 C\n",
+ "\n",
+ "# given values\n",
+ "MO2 = 8.;# mass of O2, [kg]\n",
+ "MN2 = 7.;# mass of N2, [kg]\n",
+ "MCO2 = 22.;# mass of CO2, [kg]\n",
+ "\n",
+ "P = 416.;# total pressure in the vessel, [kN/m^2]\n",
+ "T = 273.+60;# vessel temperature, [K]\n",
+ "R = 8.3143;# gas constant, [kJ/kmol K]\n",
+ "\n",
+ "mO2 = 32.;# molculer mass of O2 \n",
+ "mN2 = 28.;# molculer mass of N2\n",
+ "mCO2 = 44.;# molculer mass of CO2\n",
+ "\n",
+ "# solution\n",
+ "# (a)\n",
+ "n1 = MO2/mO2;# moles of O2, [kmol]\n",
+ "n2 = MN2/mN2;# moles of N2, [kmol]\n",
+ "n3 = MCO2/mCO2;# moles of CO2, [kmol]\n",
+ "\n",
+ "n = n1+n2+n3;# total moles in the vessel, [kmol]\n",
+ "# since,Partial pressure is proportinal, so\n",
+ "P1 = n1*P/n;# partial pressure of O2, [kN/m^2]\n",
+ "P2 = n2*P/n;# partial pressure of N2, [kN/m^2]\n",
+ "P3 = n3*P/n;# partial pressure of CO2, [kN/m^2]\n",
+ "\n",
+ "print ' (a)The partial pressure of O2 is (kN/m^2) = ',P1,',\\n The partial pressure of N2 is (kN/m^2) = ',P2,'\\n The partial pressure of CO2 is (kN/m^2) = ',P3\n",
+ "\n",
+ "# (b)\n",
+ "# assuming ideal gas \n",
+ "V = n*R*T/P;# volume of the container, [m^3]\n",
+ "print ' (b) The volume of the container is (m^3) = ',round(V,3)\n",
+ "\n",
+ "# (c)\n",
+ "T2 = 273.+228;# raised vessel temperature, [K]\n",
+ "# so volume of vessel will constant , P/T=constant\n",
+ "P2 = P*T2/T;# new pressure in the vessel , [kn/m62]\n",
+ "print ' (c) The new total pressure in the vessel is (kN/m^2) = ',round(P2)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 18: pg 225"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 25,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.18\n",
+ " The actual mass of air supplied is (kg/kg coal) = 13.64\n",
+ " The velocity of flue gas is (m/s) = 3.99\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 225\n",
+ "print('Example 8.18');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the actual mass of air supplied/kg coal\n",
+ "# the velocity of flue gas\n",
+ "\n",
+ "# given values\n",
+ "mc = 635;# mass of coal burn/h, [kg]\n",
+ "ea = .25;# excess air required\n",
+ "C = .84;# mass composition of carbon\n",
+ "H2 = .04;# mass composition of hydrogen\n",
+ "O2 = .05;# mass composition of oxygen\n",
+ "ash = 1-(C+H2+O2);# mass composition of ash\n",
+ "\n",
+ "P1 = 101.3;# pressure, [kJn/m^2]\n",
+ "T1 = 273;# temperature, [K]\n",
+ "V1 = 22.4;# volume, [m^3]\n",
+ "\n",
+ "T2 = 273.+344;# gas temperature, [K]\n",
+ "P2 = 100.;# gas pressure, [kN/m^2]\n",
+ "A = 1.1;# cross section area, [m^2]\n",
+ "aO2 = .23;# composition of O2 in air\n",
+ "\n",
+ "mCO2 = 44.;# moleculer mass of carbon\n",
+ "mH2O = 18.;# molecular mass of hydrogen\n",
+ "mO2 = 32.;# moleculer mas of oxygen\n",
+ "mN2 = 28.;# moleculer mass of nitrogen\n",
+ "\n",
+ "# solution\n",
+ "mtO2 = 8./3*C+8*H2-O2;# theoretical O2 required/kg coal, [kg]\n",
+ "msa= mtO2/aO2;# stoichiometric mass of air supplied/kg coal, [kg]\n",
+ "mas = msa*(1+ea);# actual mass of air supplied/kg coal, [kg]\n",
+ "\n",
+ "m1 = 11./3*C;# mass of CO2/kg coal produced, [kg]\n",
+ "m2 = 9*H2;# mass of H2/kg coal produced, [kg]\n",
+ "m3 = mtO2*ea;# mass of O2/kg coal produced, [kg]\n",
+ "m4 = mas*(1-aO2);# mass of N2/kg coal produced, [kg]\n",
+ "\n",
+ "mt = m1+m2+m3+m4;# total mass, [kg]\n",
+ "x1 = m1/mt*100;# %age mass composition of CO2 produced\n",
+ "x2 = m2/mt*100;# %age mass composition of H2O produced\n",
+ "x3 = m3/mt*100;# %age mass composition of O2 produced\n",
+ "x4 = m4/mt*100;# %age mass composition of N2 produced\n",
+ "\n",
+ "vt = x1/mCO2+x2/mH2O+x3/mO2+x4/mN2;# total volume\n",
+ "v1 = x1/mCO2/vt*100;# %age volume composition of CO2\n",
+ "v2 = x2/mH2O/vt*100;# %age volume composition of H2O\n",
+ "v3 = x3/mO2/vt*100;# %age volume composition of O2\n",
+ "v4 = x4/mN2/vt*100;# %age volume composition of N2\n",
+ "\n",
+ "Mav = (v1*mCO2+v2*mH2O+v3*mO2+v4*mN2)/(v1+v2+v3+v4);# average moleculer mass, [kg/kmol]\n",
+ "# since no of moles is constant so PV/T=constant\n",
+ "V2 = P1*V1*T2/(P2*T1);#volume, [m^3]\n",
+ "\n",
+ "mp = mt*mc/3600.;# mass of product of combustion/s, [kg]\n",
+ "\n",
+ "V = V2*mp/Mav;# volume of flowing gas /s,[m^3]\n",
+ "\n",
+ "v = V/A;# velocity of flue gas, [m/s]\n",
+ "print ' The actual mass of air supplied is (kg/kg coal) = ',round(mas,2)\n",
+ "print ' The velocity of flue gas is (m/s) = ',round(v,2)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 19: pg 227"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 26,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.19\n",
+ " (a) The temperature of the gas after compression is (C) = 354.1\n",
+ " (b) The density of air-gas mixture is (kg/m^3) = 1.133\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 227\n",
+ "print('Example 8.19');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# (a) the temperature of the gas after compression\n",
+ "# (b) the density of the air-gas mixture\n",
+ "\n",
+ "# given values\n",
+ "CO = 26.;# %age volume composition of CO \n",
+ "H2 = 16.;# %age volume composition of H2\n",
+ "CH4 = 7.;# %age volume composition of CH4 \n",
+ "N2 = 51.;# %age volume composition of N2\n",
+ "\n",
+ "P1 = 103.;# gas pressure, [kN/m^2]\n",
+ "T1 = 273.+21;# gas temperature, [K]\n",
+ "rv = 7.;# volume ratio\n",
+ "\n",
+ "aO2 = 21.;# %age volume composition of O2 in the air\n",
+ "c = 21.;# specific heat capacity of diatomic gas, [kJ/kg K]\n",
+ "cCH4 = 36.;# specific heat capacity of CH4, [kJ/kg K]\n",
+ "R = 8.3143;# gas constant, [kJ/kg K]\n",
+ "\n",
+ "mCO = 28.;# moleculer mass of carbon\n",
+ "mH2 = 2.;# molecular mass of hydrogen\n",
+ "mCH4 = 16.;# moleculer mas of methane\n",
+ "mN2 = 28.;# moleculer mass of nitrogen\n",
+ "mO2 = 32.;# moleculer mass of oxygen\n",
+ "\n",
+ "# solution\n",
+ "# (a)\n",
+ "Cav = (CO*c+H2*c+CH4*cCH4+N2*c+100*2*c)/(100.+200);# heat capacity, [kJ/kg K]\n",
+ "\n",
+ "Gama = (Cav+R)/Cav;# heat capacity ratio\n",
+ "# rv = V1/V2\n",
+ "# process is polytropic, so\n",
+ "T2 = T1*(rv)**(Gama-1);# final tempearture, [K]\n",
+ "print ' (a) The temperature of the gas after compression is (C) = ',round(T2-273.15,1)\n",
+ "\n",
+ "# (b)\n",
+ "\n",
+ "Mav = (CO*mCO+H2*mH2+CH4*mCH4+N2*mN2+42*mO2+158*mN2)/(100.+200)\n",
+ "\n",
+ "# for 1 kmol of gas\n",
+ "V = R*T1/P1;# volume of one kmol of gas, [m^3]\n",
+ "# hence\n",
+ "rho = Mav/V;# density of gas, [kg/m^3]\n",
+ "\n",
+ "print ' (b) The density of air-gas mixture is (kg/m^3) = ',round(rho,3)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 20: pg 228"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 27,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.20\n",
+ "The required stoichiometric equation is = \n",
+ "2 H2+ O2+3.76 N2 = 2 H2O+3.76 N2\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 228\n",
+ "print('Example 8.20');\n",
+ "\n",
+ "# aim : to determine \n",
+ "# stoichiometric equation for combustion of hydrogen\n",
+ "\n",
+ "# solution\n",
+ "# equation with algebric coefficient is\n",
+ "# H2+aO2+79/21*aN2=bH2O+79/21*aN2\n",
+ "# by equating coefficients\n",
+ "b = 1;\n",
+ "a = b/2.;\n",
+ "# so equation becomes\n",
+ "# 2 H2+ O2+3.76 N2=2 H2O+3.76 N2\n",
+ "#results\n",
+ "print('The required stoichiometric equation is = ');\n",
+ "print('2 H2+ O2+3.76 N2 = 2 H2O+3.76 N2');\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 22: pg 229"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 28,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.22\n",
+ " Percentage gravimetric composition of CO2 = 12.65 \n",
+ " ,\n",
+ " Percentage gravimetric composition of H2O = 11.7 \n",
+ "\n",
+ " Percentage gravimetric composition of O2 = 2.77 \n",
+ "\n",
+ " Percentage gravimetric composition of N2 = 72.89\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 229\n",
+ "print('Example 8.22');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the percentage gravimetric analysis of the total products of combustion\n",
+ "\n",
+ "# given values\n",
+ "CO = 12.;# %age volume composition of CO\n",
+ "H2 = 41.;# %age volume composition of H2\n",
+ "CH4 = 27.;# %age volume composition of CH4\n",
+ "O2 = 2.;# %age volume composition of O2\n",
+ "CO2 = 3.;# %age volume composition of CO2\n",
+ "N2 = 15.;# %age volume composition of N2\n",
+ "\n",
+ "mCO2 = 44.;# moleculer mass of CO2,[kg/kmol]\n",
+ "mH2O = 18.;# moleculer mass of H2O, [kg/kmol]\n",
+ "mO2 = 32.;# moleculer mass of O2, [kg/kmol]\n",
+ "mN2 = 28.;# moleculer mass of N2, [kg/kmol]\n",
+ " \n",
+ "ea = 15.;# %age excess air required\n",
+ "aO2 = 21.;# %age air composition in the air\n",
+ "\n",
+ "# solution\n",
+ "# combustion equation by no. of moles\n",
+ "# 12CO + 41H2 + 27CH4 + 2O2 + 3CO2 + 15N2 + aO2+79/21*aN2 = bCO2 + dH2O + eO2 + 15N2 +79/21*aN2\n",
+ "# equating C coefficient\n",
+ "b = 12.+27+3;# [mol]\n",
+ "# equatimg H2 coefficient\n",
+ "d = 41.+2*27;# [mol]\n",
+ "# O2 required is 15 % extra,so\n",
+ "# e/(e-a)=.15 so e=.13a\n",
+ "# equating O2 coefficient\n",
+ "# 2+3+a=b+d/2 +e\n",
+ "\n",
+ "a = (b+d/2.-5)/(1-.13);\n",
+ "e = .13*a;# [mol]\n",
+ "\n",
+ "# gravimetric analysis of product\n",
+ "v1 = b*mCO2;# gravimetric volume of CO2 \n",
+ "v2 = d*mH2O ;# gravimetric volume of H2O \n",
+ "v3 = e*mO2;# gravimetric volume of O2\n",
+ "v4 = 15*mN2 +79./21*a*mN2;# gravimetric volume of N2 \n",
+ "\n",
+ "vt = v1+v2+v3+v4;# total\n",
+ "x1 = v1/vt*100;# percentage gravimetric of CO2\n",
+ "x2 = v2/vt*100;# percentage gravimetric of H2O\n",
+ "x3 = v3/vt*100;# percentage gravimetric of O2\n",
+ "x4 = v4/vt*100;# percentage gravimetric of N2\n",
+ "#results\n",
+ "print ' Percentage gravimetric composition of CO2 = ',round(x1,2),' \\n ,\\n Percentage gravimetric composition of H2O = ',round(x2,2),' \\n\\n Percentage gravimetric composition of O2 = ',round(x3,2),' \\n\\n Percentage gravimetric composition of N2 = ',round(x4,2)\n",
+ "\n",
+ "# End \n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 23: pg 231"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 29,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.23\n",
+ " (a) The mass of actual air supplied per kg of fuel is (kg) = 31.2\n",
+ " (b) The volumetric efficiency of the engine is (percent) = 89.9\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 231\n",
+ "print('Example 8.23');\n",
+ "import math\n",
+ "# aim : To determine\n",
+ "# (a) the actual quantity of air supplied/kg of fuel\n",
+ "# (b) the volumetric efficiency of the engine\n",
+ "\n",
+ "# given values\n",
+ "d = 300.*10**-3;# bore,[m]\n",
+ "L = 460.*10**-3;# stroke,[m]\n",
+ "N = 200.;# engine speed, [rev/min]\n",
+ "\n",
+ "C = 87.;# %age mass composition of Carbon in the fuel\n",
+ "H2 = 13.;# %age mass composition of H2 in the fuel\n",
+ "\n",
+ "mc = 6.75;# fuel consumption, [kg/h]\n",
+ "\n",
+ "CO2 = 7.;# %age composition of CO2 by volume\n",
+ "O2 = 10.5;# %age composition of O2 by volume\n",
+ "N2 = 7.;# %age composition of N2 by volume\n",
+ "\n",
+ "mC = 12.;# moleculer mass of CO2,[kg/kmol]\n",
+ "mH2 = 2.;# moleculer mass of H2, [kg/kmol]\n",
+ "mO2 = 32.;# moleculer mass of O2, [kg/kmol]\n",
+ "mN2 = 28.;# moleculer mass of N2, [kg/kmol]\n",
+ "\n",
+ "T = 273.+17;# atmospheric temperature, [K]\n",
+ "P = 100;# atmospheric pressure, [kn/m**2]\n",
+ "R =.287;# gas constant, [kJ/kg k]\n",
+ "\n",
+ "# solution\n",
+ "# (a)\n",
+ "# combustion equation by no. of moles\n",
+ "# 87/12 C + 13/2 H2 + a O2+79/21*a N2 = b CO2 + d H2O + eO2 + f N2\n",
+ "# equating coefficient\n",
+ "b = 87./12;# [mol]\n",
+ "a = 22.7;# [mol]\n",
+ "e = 10.875;# [mol]\n",
+ "f = 11.8*b;# [mol]\n",
+ "# so fuel side combustion equation is\n",
+ "# 87/12 C + 13/2 H2 +22.7 O2 +85.5 N2\n",
+ "mair = ( 22.7*mO2 +85.5*mN2)/100;# mass of air/kg fuel, [kg]\n",
+ "print ' (a) The mass of actual air supplied per kg of fuel is (kg) = ',round(mair,2)\n",
+ "\n",
+ "# (b)\n",
+ "m = mair*mc/60;# mass of air/min, [kg]\n",
+ "V = m*R*T/P;# volumetric flow of air/min, [m**3]\n",
+ "SV = math.pi/4*d**2*L*N/2;# swept volume/min, [m**3]\n",
+ "\n",
+ "VE = V/SV;# volumetric efficiency\n",
+ "print ' (b) The volumetric efficiency of the engine is (percent) = ',round(VE*100,1)\n",
+ "\n",
+ "# End\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 24: pg 232"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 30,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Example 8.24\n",
+ " The mass of air supplied per kg of fuel is (kg) = 21.07\n"
+ ]
+ }
+ ],
+ "source": [
+ "#pg 232\n",
+ "print('Example 8.24');\n",
+ "\n",
+ "# aim : To determine\n",
+ "# the mass of air supplied/kg of fuel burnt\n",
+ "\n",
+ "# given values\n",
+ "# gas composition in the fuel\n",
+ "C = 84.;# %age mass composition of Carbon in the fuel\n",
+ "H2 = 14.;# %age mass composition of H2 in the fuel\n",
+ "O2f = 2.;# %age mass composition of O2 in the fuel\n",
+ "\n",
+ "# exhaust gas composition\n",
+ "CO2 = 8.85;# %age composition of CO2 by volume\n",
+ "CO = 1.2# %age composition of CO by volume\n",
+ "O2 = 6.8;# %age composition of O2 by volume\n",
+ "N2 = 83.15;# %age composition of N2 by volume\n",
+ "\n",
+ "mC = 12.;# moleculer mass of CO2,[kg/kmol]\n",
+ "mH2 = 2.;# moleculer mass of H2, [kg/kmol]\n",
+ "mO2 = 32.;# moleculer mass of O2, [kg/kmol]\n",
+ "mN2 = 28.;# moleculer mass of N2, [kg/kmol]\n",
+ "\n",
+ "# solution\n",
+ "# combustion equation by no. of moles\n",
+ "# 84/12 C + 14/2 H2 +2/32 O2 + a O2+79.3/20.7*a N2 = b CO2 + d CO2+ eO2 + f N2 +g H2\n",
+ "# equating coefficient and given condition\n",
+ "b = 6.16;# [mol]\n",
+ "a = 15.14;# [mol]\n",
+ "d = .836;# [mol]\n",
+ "f = 69.3*d;# [mol]\n",
+ "# so fuel side combustion equation is\n",
+ "# 84/12 C + 14/2 H2 +2/32 O2 + 15.14 O2 +85.5 N2\n",
+ "mair = ( a*mO2 +f*mN2)/100;# mass of air/kg fuel, [kg]\n",
+ "#results\n",
+ "print ' The mass of air supplied per kg of fuel is (kg) = ',round(mair,2)\n",
+ "\n",
+ "# End\n"
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.11"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}