diff options
Diffstat (limited to 'Basic_Engineering_Thermodynamics/ch12.ipynb')
| -rwxr-xr-x | Basic_Engineering_Thermodynamics/ch12.ipynb | 763 |
1 files changed, 763 insertions, 0 deletions
diff --git a/Basic_Engineering_Thermodynamics/ch12.ipynb b/Basic_Engineering_Thermodynamics/ch12.ipynb new file mode 100755 index 00000000..896ca7c0 --- /dev/null +++ b/Basic_Engineering_Thermodynamics/ch12.ipynb @@ -0,0 +1,763 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:53a735279bd0b4b47d350239e84293c3d52753ef579aa2d9f892dff4596cafc6" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 12 : Fuels and Combustion" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.1 Page No : 412" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "\n", + "xCO2 = 0.03 \t\t\t#mole fraction of CO2\n", + "xCO = 0.19 \t\t\t#mole fraction of CO\n", + "xH2 = 0.41 \t\t\t#mole fraction of H2\n", + "xCH4 = 0.25 \t\t\t#mole fraction of CH4\n", + "xC2H4 = 0.09 \t\t\t#mole fraction of C2H4\n", + "xN2 = 0.03 \t\t\t#mole fraction of N2\n", + "\n", + "\t\t\t\n", + "# Calculations and Results\n", + "#Part(a)\n", + "print \"Parta\";\n", + "M = xCO2*44 + xCO*28 + xH2*2 + xC2H4*28 +xCH4*16 + xN2*28 \t\t\t#kg/mol \t\t\t#average molar mass\n", + "yCO2 = xCO2*(44/M)\n", + "yCO = xCO*(28/M)\n", + "yH2 = xH2*(2/M)\n", + "yCH4 = xCH4*(16/M)\n", + "yC2H4 = xC2H4*(28/M)\n", + "yN2 = xN2*(28/M)\n", + "print \"yCO2 = %.1f%%\"%(yCO2*100)\n", + "print \"yCO = %.1f%%\"%(yCO*100)\n", + "print \"yH2 = %.1f%%\"%(yH2*100)\n", + "print \"yCH4 = %.0f%%\"%(yCH4*100)\n", + "print \"yC2H4 = %.1f%%\"%(yC2H4*100)\n", + "print \"yN2 = %.1f%%\"%(yN2*100)\n", + "\n", + "#Part(b)\n", + "print \"Partb\";\n", + "nC = xCO2 + xCO + xCH4 + xC2H4*2 \t\t\t#number of moles of C\n", + "nH2 = xH2 + xCH4*2 + xC2H4*2 \t\t\t#number of moles of H2\n", + "nO2 = xCO2 + 0.5*xCO \t\t\t#number of moles of O2\n", + "nN2 = xN2 \t\t\t#number of moles of N2\n", + "\n", + "mC = nC*12/M\n", + "mH2 = nH2*2/M\n", + "mO2 = nO2*32/M\n", + "mN2 = nN2*28/M\n", + "print \"mC = %.1f%%\"%(mC*100)\n", + "print \"mH2 = %.1f%%\"%(mH2*100)\n", + "print \"mO2 = %.0f%%\"%(mO2*100)\n", + "print \"mN2 = %.1f%%\"%(mN2*100)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Parta\n", + "yCO2 = 8.9%\n", + "yCO = 35.9%\n", + "yH2 = 5.5%\n", + "yCH4 = 27%\n", + "yC2H4 = 17.0%\n", + "yN2 = 5.7%\n", + "Partb\n", + "mC = 52.6%\n", + "mH2 = 14.7%\n", + "mO2 = 27%\n", + "mN2 = 5.7%\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.2 Page No : 416" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "#Molar masses of O2,H2,N2,C and S respectively\n", + "MO2 = 32. \t\t\t#g/mol\n", + "MH2 = 2. \t\t\t#g/mol\n", + "MN2 = 28. \t\t\t#g/mol\n", + "MC = 12. \t\t\t#g/mol\n", + "MS = 32. \t\t\t#g/mol\n", + "#Part(a)\n", + "print \"Parta\"\n", + "print \"Stoichiometric airCarbon) = %.2f kg/kg carbon\"%((MO2 + 3.76*MN2)/MC)\n", + "#Part(b)\n", + "print \"Partb\"\n", + "print \"Stoichiometric airHydrogen) = %.1f kg/kg hydrogen\"%(0.5*(MO2 + 3.76*MN2)/MH2)\n", + "#Part(c)\n", + "print \"Partc\"\n", + "print \"Stoichiometric airSulphur) = %.1f kg/kg sulphur\"%((MO2 + 3.76*MN2)/MS)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Parta\n", + "Stoichiometric airCarbon) = 11.44 kg/kg carbon\n", + "Partb\n", + "Stoichiometric airHydrogen) = 34.3 kg/kg hydrogen\n", + "Partc\n", + "Stoichiometric airSulphur) = 4.3 kg/kg sulphur\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.3 Page No : 417" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "xC3H8 = 0.2 \t\t\t#mole fraction of propane\n", + "xC4H10 = 0.8 \t\t\t#mole fraction of bumath.tane\n", + "#Molar masses of O2,H2,N2 & C respectively\n", + "MO2 = 32. \t\t\t#g/mol\n", + "MH2 = 2. \t\t\t#g/mol\n", + "MN2 = 28. \t\t\t#g/mol\n", + "MC = 12. \t\t\t#g/mol\n", + "\n", + "\n", + "# Calculations and Results\n", + "#C balance\n", + "b = xC3H8*3+xC4H10*4\n", + "#H2 balance\n", + "d = xC3H8*4 + xC4H10*5\n", + "#O2 balance\n", + "a = b + d/2\n", + "#N2 balance\n", + "c = 3.76*a\n", + "\n", + "Stoichiometric_air = a*(MO2 + 3.76*MN2)/(xC3H8*(MC+MO2)+xC4H10*(MC*4+MH2*5))\n", + "print \"Stoichiometric air = %.2f kg/kg cooking gas\"%(Stoichiometric_air)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Stoichiometric air = 15.42 kg/kg cooking gas\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.4 Page No : 418" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "xC3H8 = 0.2 \t\t\t#mole fraction of propane\n", + "xC4H10 = 0.8 \t\t\t#mole fraction of bumath.tane\n", + "\n", + "xO2 = 0.21 \t\t\t#mole fraction of O2\n", + "xN2 = 0.79 \t\t\t#mole fraction of N2\n", + "#Molar masses of O2,H2,N2 & C respectively\n", + "MO2 = 32. \t\t\t#g/mol\n", + "MH2 = 2. \t\t\t#g/mol\n", + "MN2 = 28. \t\t\t#g/mol\n", + "MC = 12. \t\t\t#g/mol\n", + "Stoichiometric_air = 15.42 \t\t\t#kg/kg\n", + "\n", + "\n", + "# Calculations and Results\n", + "#C balance\n", + "b = xC3H8*3+xC4H10*4\n", + "#H2 balance\n", + "e = xC3H8*4 + xC4H10*5\n", + "#O2 balance\n", + "d = 40*xO2-b -0.5*e\n", + "#N2 balance\n", + "c = 40*xN2\n", + "\n", + "actual_air = 40*(xO2*MO2 + xN2*MN2)/(xC3H8*(MC*3+MH2*4)+xC4H10*(MC*4+MH2*5)) \t\t\t#kg/kg gas\n", + "\n", + "excess_air = (actual_air - Stoichiometric_air)/Stoichiometric_air*100\n", + "print \"Excess air = %.1f%%\"%(excess_air)\n", + "\n", + "theoritical_air = 100+excess_air\n", + "print \"Theoritical air = %.1f%%\"%(theoritical_air)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Excess air = 35.5%\n", + "Theoritical air = 135.5%\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.5 Page No : 420" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "nCO2 = 3.8 \t\t\t#Number of moles of CO2\n", + "nN2 = 31.6 \t\t\t#Number of moles of N2\n", + "nO2 = 2.2 \t\t\t#Number of moles of O2\n", + "\n", + "\t\t\t\n", + "# Calculations and Results\n", + "n = nCO2 + nN2 + nO2\n", + "xCO2 = nCO2/n \t\t\t#mole fraction of CO2\n", + "xN2 = nN2/n \t\t\t#mole fraction of N2\n", + "xO2 = nO2/n \t\t\t#mole fraction of O2\n", + "\n", + "print \"xCO2 = %.1f %% \"%(xCO2*100)\n", + "print \"xN2 = %.1f %% \"%(xN2*100)\n", + "print \"xO2 = %.1f %% \"%(xO2*100)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "xCO2 = 10.1 % \n", + "xN2 = 84.0 % \n", + "xO2 = 5.9 % \n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.6 Page No : 421" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "nCO2 = 3.8 \t\t\t#Number of moles of CO2\n", + "nN2 = 31.6 \t\t\t#Number of moles of N2\n", + "nO2 = 2.2 \t\t\t#Number of moles of O2\n", + "nH2O = 4.8 \t\t\t#Number of moles of H2O\n", + "\n", + "\t\t\t\n", + "# Calculations and Results\n", + "\t\t\t\n", + "# Calculations and Results\n", + "n = nCO2 + nN2 + nO2 + nH2O\n", + "xCO2 = nCO2/n \t\t\t#mole fraction of CO2\n", + "xN2 = nN2/n \t\t\t#mole fraction of N2\n", + "xO2 = nO2/n \t\t\t#mole fraction of O2\n", + "xH2O = nH2O/n \t\t\t#mole fraction of H2O\n", + "\n", + "print \"Volume fraction of CO2 = %.1f%%\"%(xCO2*100)\n", + "print \"Volume fraction of N2 = %.1f%%\"%(xN2*100)\n", + "print \"Volume fraction of O2 = %.1f%%\"%(xO2*100)\n", + "print \"Volume fraction of H2O = %.1f%%\"%(xH2O*100)\n", + "\n", + "M = xCO2*44 + xN2*28 + xO2*32 + xH2O*18 \t\t\t#Mass of combustion product\n", + "\n", + "yCO2 = xCO2*(44/M)\n", + "yH2O = xH2O*(18/M)\n", + "yN2 = xN2*(28/M)\n", + "yO2 = xO2*(32/M)\n", + "\n", + "print \"Mass fraction of CO2 = %.1f%%\"%(yCO2*100)\n", + "print \"Mass fraction of N2 = %.1f%%\"%(yN2*100)\n", + "print \"Mass fraction of O2 = %.1f%%\"%(yO2*100)\n", + "print \"Mass fraction of H2O = %.1f%%\"%(yH2O*100)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Volume fraction of CO2 = 9.0%\n", + "Volume fraction of N2 = 74.5%\n", + "Volume fraction of O2 = 5.2%\n", + "Volume fraction of H2O = 11.3%\n", + "Mass fraction of CO2 = 13.8%\n", + "Mass fraction of N2 = 73.2%\n", + "Mass fraction of O2 = 5.8%\n", + "Mass fraction of H2O = 7.1%\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.7 Page No : 422" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Part(a)\n", + "print \"Parta\"\n", + "mC = 0.65 \t\t\t#kg \t\t\t#mass of C per kg coal\n", + "mA = 0.15 \t\t\t#kg \t\t\t#mass of Ash per kg coal\n", + "CR = 0.05 \t\t\t#kg \t\t\t#mass of carbon in solid refuse per kg coal\n", + "mR = 0.2 \t\t\t#kg \t\t\t#mass of refuse per kg coal\n", + "m = mC- CR \t\t\t#kg \t\t\t#mass of carbon burnt per kg coal\n", + "\n", + "#By C balance\n", + "x = (14 + 1)*(12/0.6) \t\t\t#kg \t\t\t#mass of burnt coal\n", + "#By H2 balance\n", + "b = x*(0.06/2)\n", + "#By O2 Balance\n", + "a = (14 + 0.5 + 3.5 + 4.5)-(x*0.1/32) \n", + "actual_air = a*(32+3.76*28)/x \t\t\t#kg/kg coal\n", + "print \"Actual air = %.3f kg/kg coal\"%(actual_air)\n", + "\n", + "#Part(b)\n", + "print \"Partb\"\n", + "Stoichiometric_air = (0.6*11.45+0.06*34.3)-(0.1/0.232) \t\t\t#kg\n", + "excess_air = (actual_air - Stoichiometric_air)/Stoichiometric_air*100 \n", + "print \"Excess air = %.1f%%\"%(excess_air)\n", + "\n", + "#Part(c)\n", + "print \"Partc\";\n", + "print \"Percentage theoritical air = %.1f%%\"%(100+excess_air)\n", + "\n", + "#Part(d)\n", + "print \"Partd\"\n", + "m = 14*44 + 1*28 +3.5*32 +81.5*28 +9*18 \t\t\t#kg \t\t\t#mass of combustion product\n", + "print \"Mass fraction of CO2 = %.2f%%\"%(14*44/m*100)\n", + "print \"Mass fraction of CO = %.3f%%\"%(1*28/m*100)\n", + "print \"Mass fraction of O2 = %.2f%%\"%(3.5*32/m*100)\n", + "print \"Mass fraction of N2 = %.2f%%\"%(81.5*28/m*100)\n", + "print \"Mass fraction of H2O = %.2f%%\"%(9*18/m*100)\n", + "\n", + "#Part(e)\n", + "print \"Parte\"\n", + "xH2O = 9/(14+1+3.5+81.5+9) \t\t\t#molfraction of H2O\n", + "pH2O = xH2O*1e5 \t\t\t#Pa \t\t\t#partial pressure\n", + "#From steam table\n", + "tdp = 42.5 \t\t\t#\u00b0C\n", + "print \"Dew point temperature = %.1f \u00b0C\"%(tdp)\n", + "\n", + "\n", + "# note : part b is calculated wrong in book. so answers might be different." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Parta\n", + "Actual air = 9.867 kg/kg coal\n", + "Partb\n", + "Excess air = 16.1%\n", + "Partc\n", + "Percentage theoritical air = 116.1%\n", + "Partd\n", + "Mass fraction of CO2 = 19.25%\n", + "Mass fraction of CO = 0.875%\n", + "Mass fraction of O2 = 3.50%\n", + "Mass fraction of N2 = 71.31%\n", + "Mass fraction of H2O = 5.06%\n", + "Parte\n", + "Dew point temperature = 42.5 \u00b0C\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.8 Page No : 425" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "xCO2 = 9.7 \t\t\t#mole percent CO2\n", + "xCO = 1.1 \t\t\t#mole percent CO\n", + "xO2 = 4 \t\t\t#mole percent O2\n", + "xN2 = 85.2 \t\t\t#mole percent N2\n", + "\n", + "\n", + "# Calculations and Results\n", + "#by C balance\n", + "b = 2\n", + "#by H2 balance\n", + "d = 2\n", + "#by O2 balance\n", + "a = b+d*.5\n", + "#by N2 balance\n", + "c = 3.76*a\n", + "\n", + "Stoichiometric_air = a*(32+3.76*28)/28 \t\t\t#kg/kg ethylene\n", + "\n", + "#by C balance\n", + "x = (xCO2+xCO)/2 \t\t\t#kmol of ehtylene be burnt\n", + "#by H2 balance\n", + "q = 2*x\n", + "#by O2 balance\n", + "p = xCO2 + xCO/2 + xO2 + q/2\n", + "\n", + "actual_air = p*(32+3.76*28)/(x*28) \t\t\t#kg/kg ethylene\n", + "excess_air = (actual_air - Stoichiometric_air)/Stoichiometric_air*100\n", + "print \"Actual air = %.1f kg/kg ethylene \"%(actual_air)\n", + "print \"Excess air = %.0f%%\"%(excess_air)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Actual air = 17.8 kg/kg ethylene \n", + "Excess air = 21%\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.9 Page No : 433" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# variables\n", + "excess_air = 10. \t\t\t#%\n", + "tR = 30. \t\t\t#\u00b0C \t\t\t#temperature of reacmath.tant\n", + "tP = 120. \t\t\t#\u00b0C \t\t\t#temperature of product\n", + "delta_H = -802.3*1000 \t\t\t#kJ/kmol CH4\n", + "Cp_CH4 = 36. \t\t\t#kJ/lmolK\n", + "Cp_O2 = 29.5 \t\t\t#kJ/lmolK\n", + "Cp_N2 = 29.2 \t\t\t#kJ/lmolK\n", + "Cp_CO2 = 37. \t\t\t#kJ/lmolK\n", + "Cp_H2O = 33.7 \t\t\t#kJ/lmolK\n", + "tA = 30. \t\t\t#\u00b0C\n", + "tX = 25. \t\t\t#\u00b0C\n", + "tY = tX \t\t\t#\u00b0C\n", + "tB = 120. \t\t\t#\u00b0C\n", + "\t\t\t\n", + "# Calculations and Results\n", + "Q_AB = (1*Cp_CO2 + 0.2*Cp_O2 + 8.272*Cp_N2 + 2*Cp_H2O)*(tB-tX)+ delta_H + (1*Cp_CH4 + 2.2*Cp_O2 +8.272*Cp_N2)*(tY-tA)\n", + "print \"The heat transfer from the combustor = %.1f kJ/kg CH4 \"%(Q_AB/16)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The heat transfer from the combustor = -48161.7 kJ/kg CH4 \n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.10 Page No : 435" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "deltaH_gasgas = -2651.4e3 \t\t\t#kJ/kmol bumath.tane\n", + "hfg_butane = 370 \t\t\t#kJ/kg \t\t\t#enthalpies of vaporisation of butance\n", + "hfg_water = 2442 \t\t\t#kJ/kg \t\t\t#enthalpies of vaporisation of water\n", + "M_butane = 58 \t\t\t#g/mol\n", + "M_water = 18 \t\t\t#g/mol\n", + "\n", + "\t\t\t\n", + "# Calculations and Results\n", + "deltaH_liqliq = deltaH_gasgas + M_butane*hfg_butane - 5*M_water*hfg_water\n", + "print \"Enthalpy of reaction = %.1f kJ/kg\"%(deltaH_liqliq/M_butane)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Enthalpy of reaction = -49133.1 kJ/kg\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.11 Page No : 437" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "#Formation Enthalpies\n", + "hf_CO2 = -393510. \t\t\t#kJ/kmol\n", + "hf_H2O = -285838. \t\t\t#kJ/kmol\n", + "hf_C3H8 = -104680. \t\t\t#kJ/kmol\n", + "hf_O2 = 0. \t\t\t#kJ/kmol\n", + "\n", + "\t\t\t\n", + "# Calculations and Results\n", + "Hp = 3*hf_CO2 + 4*hf_H2O \n", + "Hr = hf_C3H8 + 5*hf_O2\n", + "deltaH_std = Hp-Hr\n", + "print \"Standard change in enthalpy for the reaction = %.0f kJ/kg\"%(deltaH_std/44)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Standard change in enthalpy for the reaction = -50436 kJ/kg\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.12 Page No : 438" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "Cp_C2H4 = 28*1.548 \t\t\t#kJ/lmolK\n", + "Cp_O2 = 32*0.922 \t\t\t#kJ/lmolK\n", + "Cp_N2 = 28*1.042 \t\t\t#kJ/lmolK\n", + "Cp_CO2 = 44*0.842 \t\t\t#kJ/lmolK\n", + "Cp_H2O = 18*1.86 \t\t\t#kJ/lmolK\n", + "\n", + "\t\t\t\n", + "# Calculations and Results\n", + "deltaH_BX = (2*Cp_CO2 + 2*Cp_H2O + 0.3*Cp_O2 + 12.408*Cp_N2)*(120-25)\n", + "deltaH_YA = (Cp_C2H4 + 3.3*Cp_O2 + 12.408*Cp_N2)*(25-50)\n", + "hf_CO2 = -393510. \t\t\t#kJ/kmol\n", + "hf_H2O = -241820. \t\t\t#kJ/kmol\n", + "hf_C2H4 = 52283. \t\t\t#kJ/kmol\n", + "deltaH_XY = 2*hf_CO2 + 2*hf_H2O - hf_C2H4 \t\t\t#kJ/kmol\n", + "deltaH_BA = deltaH_BX + deltaH_YA + deltaH_XY \t\t\t#kJ/kmol\n", + "\n", + "print \"The heat transferred from the combustor per kg ethylene = %.0f kJ/kg ethane\"%(deltaH_BA/28)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The heat transferred from the combustor per kg ethylene = -45960 kJ/kg ethane\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.13 Page No : 441" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Part(a)\n", + "print \"Parta\"\n", + "#Picking up various Cp values from Table 12.1\n", + "tB = (16*50010)/(44*0.842+2*18*1.86+3*32*0.922+22.56*28*1.042)+25\n", + "print \"Adiabatic combustion temperature using average Cp values)= %.1f K\"%(tB+273)\n", + "\n", + "#Part(b)\n", + "print \"Partb\"\n", + "tb1 = 1000. \t\t\t#K \t\t\t#first guess temperature\n", + "tb2 = 1200. \t\t\t#K second guess temperature\n", + "tb = (tb1 - tb2)/(637617-836847)*(800160-836847) + tb2\n", + "print \"Adiabatic combustion temperature using ideal gas enthalpy = %.0f K\"%(tb)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Parta\n", + "Adiabatic combustion temperature using average Cp values)= 1238.6 K\n", + "Partb\n", + "Adiabatic combustion temperature using ideal gas enthalpy = 1163 K\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.14 Page No : 443" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "delta_H_std = -45101. \t\t\t#kJ/kg\n", + "hfg = 2442. \t\t\t#kJ/kg \t\t\t#enthalpy of vaporisation\n", + "\n", + "\t\t\t\n", + "# Calculations and Results\n", + "LCV = -1*delta_H_std \t\t\t# kJ/kg hexane\n", + "print \"LCV of gaseous hexane = %.1f kJ/kg hexane\"%(LCV)\n", + "\n", + "m = 7.*18./86 \t\t\t#mass of H2O per kg hexane\n", + "HCV = LCV+m*hfg \t\t\t#kJ/kg hexane\n", + "print \"HCV of gaseous hexane = %d kJ/kg hexane\"%(HCV)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "LCV of gaseous hexane = 45101.0 kJ/kg hexane\n", + "HCV of gaseous hexane = 48678 kJ/kg hexane\n" + ] + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file |