{
 "metadata": {
  "name": "",
  "signature": "sha256:b48de029dcc8d0a67a0785422dda0e79a8ad980134c0bbc3f290af08b9a74588"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 16 - Combustion"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1 - Pg 299"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the Molecule\n",
      "#Initialization of variables\n",
      "per=85.\n",
      "#calculations\n",
      "a=per/12.\n",
      "b=100-per\n",
      "ad=1.13*a\n",
      "bd=1.13*b+1\n",
      "#results\n",
      "print '%s %d %s %d' %(\"Molecule is C\",ad,\"H\",bd)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Molecule is C 8 H 17\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 2 - Pg 299"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the Oxygen and Nitrogen content\n",
      "#Initialization of variables\n",
      "per=0.071\n",
      "#calculations\n",
      "O2=8.74\n",
      "N2=per/2 + 3.76*O2\n",
      "#results\n",
      "print '%s %.2f %s %.2f' %(\"Oxygen =\",O2,\"and Nitrogen =\",N2)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Oxygen = 8.74 and Nitrogen = 32.90\n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 4 - Pg 302"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the air fuel ratio\n",
      "#Initialization of variables\n",
      "N2=78.1\n",
      "M=29.\n",
      "co2=8.7\n",
      "co=8.9\n",
      "x4=0.3\n",
      "x5=3.7\n",
      "x6=14.7\n",
      "#calculations\n",
      "O2=N2/3.76\n",
      "Z=(co2+co+x4)/8\n",
      "AF=(O2+N2)*M/(Z*113)\n",
      "#results\n",
      "print '%s %.1f %s' %(\"Air fuel ratio =\",AF,\"lbm air/lbm fuel\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Air fuel ratio = 11.3 lbm air/lbm fuel\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 5 - Pg 303"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the air fuel ratio\n",
      "#Initialization of variables\n",
      "N2=78.1\n",
      "M=29\n",
      "ba=2.12\n",
      "x4=0.3\n",
      "x5=3.7\n",
      "x6=14.7\n",
      "#calculations\n",
      "O2=N2/3.76\n",
      "O2=N2/3.76\n",
      "Z=(x4*4+x5*2+x6*2)/17\n",
      "AF=(O2+N2)*M/(Z*113)\n",
      "#results\n",
      "print '%s %.1f %s' %(\"Air fuel ratio =\",AF,\"lbm air/lbm fuel\")\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Air fuel ratio = 11.4 lbm air/lbm fuel\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 6 - Pg 303"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the air fuel ratio\n",
      "#Initialization of variables\n",
      "N2=78.1\n",
      "M=29\n",
      "ba=2.12\n",
      "x4=0.3\n",
      "x5=3.7\n",
      "x6=14.7\n",
      "#calculations\n",
      "O2=N2/3.76\n",
      "c=14.7\n",
      "b= x4*4 + x5*2 + x6*2\n",
      "a=b/ba\n",
      "AF=(O2+N2)*M/(a*12. + b)\n",
      "#results\n",
      "print '%s %.1f %s' %(\"Air fuel ratio =\",AF,\"lbm air/lbm fuel\")"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Air fuel ratio = 11.3 lbm air/lbm fuel\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 7 - Pg 304"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the air fuel ratio\n",
      "#Initialization of variables\n",
      "N2=78.1\n",
      "M=29\n",
      "ba=2.12\n",
      "co2=8.7\n",
      "co=8.9\n",
      "x4=0.3\n",
      "x5=3.7\n",
      "x6=14.7\n",
      "#calculations\n",
      "O2=N2/3.76\n",
      "c=14.7\n",
      "Z=2.238\n",
      "X=(Z*17-x4*4-x5*2)/2\n",
      "a=co2+co/2+x4+x6/2\n",
      "b=3.764*a\n",
      "AF=(O2+N2)*M/(Z*113)\n",
      "#results\n",
      "print '%s %.1f %s' %(\"Air fuel ratio =\",AF,\"lbm air/lbm fuel\")\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Air fuel ratio = 11.3 lbm air/lbm fuel\n"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8 - Pg 305"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the air fuel ratio\n",
      "#Initialization of variables\n",
      "x1=8.7\n",
      "x2=8.9\n",
      "x3=0.3\n",
      "N=78.1\n",
      "z=113\n",
      "M=29\n",
      "#calculations\n",
      "co2=(x1+x2+x3)*100/(N+x1+x2+x3)\n",
      "a=2.325\n",
      "AF=103*M/(a*z)\n",
      "#results\n",
      "print '%s %.2f' %(\"Air fuel ratio = \",AF)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Air fuel ratio =  11.37\n"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 9 - Pg 308"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the Higher heating value\n",
      "#Initialization of variables\n",
      "dH=-2369859 #Btu\n",
      "r=1.986\n",
      "dn=5.5\n",
      "T=536.7 #R\n",
      "#calculations\n",
      "dQ=dH+dn*r*T\n",
      "#results\n",
      "print '%s %d %s' %(\"Higher heating value =\",dQ,\" Btu\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Higher heating value = -2363996  Btu\n"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 10 - Pg 308"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the Lower heating value\n",
      "#Initialization of variables\n",
      "y=13\n",
      "x=12\n",
      "M2=18\n",
      "M=170\n",
      "p=0.4593\n",
      "vfg=694.9\n",
      "J=778.2\n",
      "m=9*18\n",
      "u1=-2363996 #Btu\n",
      "#calculations\n",
      "z=y*M2/M\n",
      "hfg=1050.4 #Btu/lbm\n",
      "ufg= hfg- p*vfg*144/J\n",
      "dU=ufg*m \n",
      "Lhv=u1+dU\n",
      "#results\n",
      "print '%s %d %s' %(\"Lower heating value =\",Lhv,\"Btu/lbm\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Lower heating value = -2203398 Btu/lbm\n"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 11 - Pg 309"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the Heat of the reaction\n",
      "#Initialization of variables\n",
      "n1=8\n",
      "n2=9\n",
      "n3=1\n",
      "n4=12.5\n",
      "U11=3852\n",
      "U12=115\n",
      "U21=3009\n",
      "U22=101\n",
      "U31=24773\n",
      "U32=640\n",
      "U41=2539\n",
      "U42=83\n",
      "H=-2203389\n",
      "#calculations\n",
      "dU1=n1*(U11-U12)+n2*(U21-U22)\n",
      "dU2=n3*(U31-U32)+n4*(U41-U42)\n",
      "Q=H+dU1-dU2\n",
      "#results\n",
      "print '%s %d %s' %(\"Heat of reaction =\",Q,\"Btu\")\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Heat of reaction = -2202154 Btu\n"
       ]
      }
     ],
     "prompt_number": 12
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12 - Pg 310"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the final temperature\n",
      "#Initialization of variables\n",
      "n1=8\n",
      "n2=9\n",
      "n3=47\n",
      "h1=118\n",
      "h2=104\n",
      "h3=82.5\n",
      "Q=2203279 #Btu\n",
      "#calculations\n",
      "U11=n1*h1+n2*h2+n3*h3\n",
      "U12=U11+Q\n",
      "T2=5271 #R\n",
      "#results\n",
      "print '%s %d %s' %(\"Upon interpolating, T2 =\",T2,\" R\")\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Upon interpolating, T2 = 5271  R\n"
       ]
      }
     ],
     "prompt_number": 13
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 13 - Pg 313"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the percentage of dissociation\n",
      "#Initialization of variables\n",
      "import numpy\n",
      "from numpy import roots\n",
      "kp=5. \n",
      "#calculations\n",
      "p1=[24.,0,3,-2]\n",
      "vec=numpy.roots(p1)\n",
      "x=numpy.real(vec[2]*100.)\n",
      "p2=[249.,0,3,-2]\n",
      "vec2=numpy.roots(p2)\n",
      "y=numpy.real(vec2[2]*100.)\n",
      "#results\n",
      "print '%s %.1f %s' %(\"percentage of dissociation =\",x,\" percent\")\n",
      "print '%s %.d %s' %(\"\\n If pressure =10 . degree of dissociation =\",y,\"percent\")\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "percentage of dissociation = 34.3  percent\n",
        "\n",
        " If pressure =10 . degree of dissociation = 18 percent\n"
       ]
      }
     ],
     "prompt_number": 14
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 14 - Pg 314"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#calculate the Extent of reaction\n",
      "#Initialization of variables\n",
      "import numpy\n",
      "from numpy import roots\n",
      "p=[24.,48,7,-4]\n",
      "vec=numpy.roots(p)\n",
      "x=numpy.real(vec[2] *100)\n",
      "#results\n",
      "print '%s %d %s' %(\"Extent of reaction=\",100-x,\"percent\")\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Extent of reaction= 78 percent\n"
       ]
      }
     ],
     "prompt_number": 15
    }
   ],
   "metadata": {}
  }
 ]
}