{
 "metadata": {
  "name": "",
  "signature": "sha256:c97b4be5a6b14d1c0295c6b4a641b54921c4d1ce87d0ba32a0413fee26fe8514"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 15.1  page no : 267"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%matplotlib inline\n",
      "\n",
      "from numpy import array,round\n",
      "import math \n",
      "\n",
      "def MCPS(T0,T,A,B,C,D):\n",
      "    t = T/T0;\n",
      "    return (A)+(((B*T0)+(((C*T0*T0)+(D/(t*t*T0*T0)))*(t+1)/2))*((t-1)/math.log(t)))\n",
      "\n",
      "def MCPH(T0,T,A,B,C,D):\n",
      "    t = T/T0;\n",
      "    return (A+((B/2)*T0*(t+1))+((C/3)*T0*T0*((t**2)+t+1))+(D/(t*T0*T0)))\n",
      "\n",
      "# Variables\n",
      "State = ['Supercooled Liquid','Superheated Vapor','Wet Vapor,x = 0.9378','Saturated Liqiud'];\n",
      "T = [318.98, 773.15, 318.98, 318.98];\n",
      "P = [8600, 8600, 10, 10];\n",
      "H = [203.4, 3391.6, 2436, 191.8];\n",
      "S = [0.6580, 6.6858, 7.6846, 0.6493];\n",
      "T0 = 298.15;\n",
      "T1 = 460.;\t\t\t#[K]\n",
      "R = 8.314;\n",
      "T_sigma = T0;\n",
      "#CH4 + 2O2 --> CO2 + 2H2O\n",
      "dH_CO2 = -393509;\n",
      "dH_H2O = -241818;\n",
      "dH_CH4 = -74520;\n",
      "\n",
      "dG_CO2 = -394359;\n",
      "dG_H2O = -228572;\n",
      "dG_CH4 = -50460;\n",
      "\n",
      "# Calculations\n",
      "dH_298 = dH_CO2+(2*dH_H2O)-dH_CH4\n",
      "dG_298 = dG_CO2+(2*dG_H2O)-dG_CH4\n",
      "\n",
      "dS_298 = round((dH_298-dG_298)/T0,3);\n",
      "\n",
      "#Moles Entering\n",
      "ni_O2 = 2*1.25;\n",
      "ni_N2 = round(ni_O2*(79/21),3);\n",
      "ni = ni_O2+ni_N2;\n",
      "\n",
      "#Moles After Combustion\n",
      "n_CO2 = 1;\n",
      "n_H2O = 2;\n",
      "n_O2 = 0.5;\n",
      "n_N2 = ni_N2;\n",
      "n = n_CO2+n_H2O+n_O2+n_N2;\n",
      "m = [n_CO2, n_H2O, n_N2, n_O2];\n",
      "\n",
      "y_CO2 = round(n_CO2/n,4);\n",
      "y_H2O = round(n_H2O/n,4);\n",
      "y_O2 = round(n_O2/n,4);\n",
      "y_N2 = round(n_N2/n,4);\n",
      "\n",
      "y = [y_CO2, y_H2O, y_O2, y_N2];\n",
      "yT = sum(y);\n",
      "\n",
      "#Step(a)\n",
      "dH_a = 0\n",
      "dS_a = round(ni*R*((0.21*math.log(0.21))+(0.79*math.log(0.79))),3)\t\t\t#[J/K]\n",
      "\n",
      "#Step(b)\n",
      "dH_b = dH_298\n",
      "dS_b = dS_298\t\t\t#[J/K]\n",
      "\n",
      "#Step(c)\n",
      "dH_c = 0\n",
      "dS_c = round(-n*R*sum(y*log(y)),3)\t\t\t#[J/K]\n",
      "\n",
      "#Step(d)\n",
      "#For CO2\n",
      "print T0,T1\n",
      "CpH_CO2 = round(R*MCPH(T0,T1,5.457,1.045*(10**-3),0,-1.157*(10**5)),3);\n",
      "#For H2O\n",
      "CpH_H2O = round(R*MCPH(T0,T1,3.470,1.450*(10**-3),0,0.121*(10**5)),3);\n",
      "#For O2\n",
      "CpH_O2 = round(R*MCPH(T0,T1,3.639,0.506*(10**-3),0,-0.227*(10**5)),3);\n",
      "#For N2\n",
      "CpH_N2 = round(R*MCPH(T0,T1,3.280,0.593*(10**-3),0,0.040*(10**5)),3);\n",
      "\n",
      "#For CO2\n",
      "CpS_CO2 = round(R*MCPS(T0,T1,5.457,1.045*(10**-3),0,-1.157*(10**5)),3);\n",
      "#For H2O\n",
      "CpS_H2O = round(R*MCPS(T0,T1,3.470,1.450*(10**-3),0,0.121*(10**5)),3);\n",
      "#For O2\n",
      "CpS_O2 = round(R*MCPS(T0,T1,3.639,0.506*(10**-3),0,-0.227*(10**5)),3);\n",
      "#For N2\n",
      "CpS_N2 = round(R*MCPS(T0,T1,3.280,0.593*(10**-3),0,0.040*(10**5)),3);\n",
      "\n",
      "CpH = array([CpH_CO2, CpH_H2O, CpH_N2, CpH_O2]);\n",
      "CpS = array([CpS_CO2, CpS_H2O, CpS_N2, CpS_O2]);\n",
      "\n",
      "Comp = ['CO2' 'H2O' 'N2' 'O2'];\n",
      "\n",
      "Ans = [CpH,CpS];\n",
      "print '    CpH      CpS',Comp,Ans\n",
      "\n",
      "CpHt = round(sum(m*CpH),3)\t\t\t#[J/K]\n",
      "CpSt = round(sum(m*CpS),3)\t\t\t#[J/K]\n",
      "\n",
      "dH_d = round(CpHt*(T1-T0),0)\t\t\t#[J]\n",
      "dS_d = round((CpSt*math.log(T1/T0)),3)\t\t\t#[J/K]\n",
      "\n",
      "dH = dH_a+dH_b+dH_c+dH_d\t\t\t#[J]\n",
      "dS = dS_a+dS_b+dS_c+dS_d\t\t\t#[J/K]\n",
      "\n",
      "rm = 84.75;\t\t\t#[kg/s]\n",
      "\n",
      "rn_CH4 = round((rm*(H[0]-H[1])*1000)/dH,2)\t\t\t#[mol/s]\n",
      "\n",
      "rW_ideal = round(rn_CH4*((dH/1000)-(T0*dS/1000))/1000,2)*1000\t\t\t#[KW]\n",
      "\n",
      "#(a)  Furnace/Boiler\n",
      "rS_a = round((rn_CH4*dS/1000)+(rm*(S[1]-S[0])),2)\t\t\t#[kJ/s/K]\n",
      "\n",
      "rW_a = round(T_sigma*rS_a/1000,2)*1000\t\t\t#[kW]\n",
      "\n",
      "#(b) Turbine\n",
      "rS_b = round(rm*(S[2]-S[1]),2)\t\t\t#[kW/K]\n",
      "\n",
      "rW_b = round(T_sigma*rS_b/1000,2)*1000\t\t\t#[kW]\n",
      "\n",
      "#(c) Condenser \n",
      "Q_c = H[3]-H[2];\t\t\t#[kJ/kg]\n",
      "rQ_c = round(rm*Q_c/1000,1)*1000\t\t\t#[kJ/s]\n",
      "rS_c = round((rm*(S[3]-S[2]))-(rQ_c/T_sigma),2)\t\t\t#[kW/K]\n",
      "rW_c = round(T_sigma*rS_c/1000,2)*1000\t\t\t#[kW]\n",
      "\n",
      "#(d) Pump\n",
      "rS_d = round(rm*(S[0]-S[3]),2)\t\t\t#[kW/K]\n",
      "rW_d = round(T_sigma*rS_d/1000,2)*1000\t\t\t#[kW]\n",
      "\n",
      "rS = [rS_a, rS_b, rS_c, rS_d];\n",
      "pS = round(rS/sum(rS)*100,1);\n",
      "T = [sum(rS), sum(pS)];\n",
      "Process = ['Furnace/boiler' 'Turbine' 'Condenser' 'Pump'];\n",
      "Ans = [rS,pS];\n",
      "print '   S(kW/K)   %',Process,Ans\n",
      "print (T)\n",
      "rW_ideal = 80000;\n",
      "rW = array([rW_ideal, rW_a ,rW_b, rW_c, rW_d])/1000;\n",
      "pW = round(rW/sum(rW)*100,1);\n",
      "T = [sum(rW) ,sum(pW)];\n",
      "Process = ['Ideal' 'Furnace/boiler' 'Turbine' 'Condenser' 'Pump'];\n",
      "Ans = [rW,pW];\n",
      "print ' W(kW/K)*10**-3   %',Process,Ans\n",
      "print (T)\n",
      "\n",
      "eta = pW[0]\n",
      "\n",
      "print 'Efficiency of the power plant is',eta,'%'\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "298.15 460.0\n",
        "    CpH      CpS ['CO2H2ON2O2'] [array([ 41.649,  34.153,  29.381,  30.473]), array([ 41.377,  34.106,  29.36 ,  30.405])]\n",
        "   S(kW/K)   % ['Furnace/boilerTurbineCondenserPump'] [[578.37, 84.650000000000006, 41.689999999999998, 0.73999999999999999], array([ 82. ,  12. ,   5.9,   0.1])]\n",
        "[705.45000000000005, 100.0]\n",
        " W(kW/K)*10**-3   % ['IdealFurnace/boilerTurbineCondenserPump'] [array([  80.  ,  172.44,   25.24,   12.43,    0.22]), array([ 27.6,  59.4,   8.7,   4.3,   0.1])]\n",
        "[290.33000000000004, 100.09999999999999]\n",
        "Efficiency of the power plant is 27.6 %\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 15.2  page no : 269"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "from numpy import array\n",
      "\n",
      "State = ['Superheated Vapor','Superheated Vapor','Superheated Vapor','Saturated Liqiud','Saturated Vapor','Superheated Vapor'];\n",
      "T = [300 ,300, 207.1, 111.5, 111.5, 295];\n",
      "P = [1, 60, 60, 1, 1, 1];\n",
      "H = [1199.8, 1140, 772, 285.4, 796.9, 1188.9];\n",
      "S = [11.629, 9.359, 7.798, 4.962, 9.523, 11.589];\n",
      "Given = [T,P,H,S];\n",
      "print '   T/K    P/kPa    H/kJ/Kg    S/kJ/kg/K',State,Given\n",
      "T_sigma = 300;\t\t\t#[K]\n",
      "rQ = 5;     \t\t\t#[KJ]\n",
      "rW = 1000;\t    \t\t#[KJ/s]\n",
      "z = round((H[5]-H[1]-rQ)/(H[5]-H[3]),4);\n",
      "\n",
      "#rW_ideal  =  (dH*rm) - (T_sigma(dS*rm))\n",
      "rW_ideal = round(((z*H[3])+((1-z)*H[5])-H[0])-((T_sigma)*((z*S[3])+((1-z)*S[5])-S[0])),1);\n",
      "\n",
      "#(a) Compression/cooling\n",
      "rQ_a = (H[1]-H[0])-rW;\t\t\t#[kJ]\n",
      "rS_a = round((S[1]-S[0])-(rQ_a/T_sigma),4);\t\t\t#[kJ/Kg/K]\n",
      "\n",
      "rW_a = T_sigma*rS_a;\t\t\t#[KJ/Kg]\n",
      "\n",
      "#(b) Exchanger\n",
      "rS_b = round(((S[5]-S[4])*(1-z))+(S[2]-S[1])-(rQ/T_sigma),4);\t\t\t#[kJ/Kg/K]\n",
      "rW_b = T_sigma*rS_b;\t\t\t#[KJ/Kg]\n",
      "\n",
      "#(c) Throttle\n",
      "rS_c = round(((S[3]*z)+(S[4]*(1-z))-S[2]),4);\t\t\t#[KJ/Kg/K]\n",
      "rW_c = T_sigma*rS_c;\t\t\t#[KJ/kg]\n",
      "\n",
      "S = [rS_a, rS_b, rS_c];\n",
      "pS = round((S/sum(S))*100,1);\n",
      "ES = [sum(S) ,sum(pS)];\n",
      "\n",
      "W = round(array([rW_ideal, rW_a ,rW_b, rW_c]),1);\n",
      "pW = round((W/sum(W))*100,1);\n",
      "EW = [sum(W), sum(pW)];\n",
      "Ans = [S,pS];\n",
      "Process = ['Compression/Cooling','Exchanger','Throttle'];\n",
      "\n",
      "print '      Si     %',Process,Ans\n",
      "print 'Sum',ES\n",
      "Ans = [W,pW];\n",
      "\n",
      "Process = ['Ideal','Compression/Cooling','Exchanger','Throttle'];\n",
      "print '      Wi     %'\n",
      "print Process\n",
      "print Ans\n",
      "print 'Sum',EW\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "   T/K    P/kPa    H/kJ/Kg    S/kJ/kg/K ['Superheated Vapor', 'Superheated Vapor', 'Superheated Vapor', 'Saturated Liqiud', 'Saturated Vapor', 'Superheated Vapor'] [[300, 300, 207.1, 111.5, 111.5, 295], [1, 60, 60, 1, 1, 1], [1199.8, 1140, 772, 285.4, 796.9, 1188.9], [11.629, 9.359, 7.798, 4.962, 9.523, 11.589]]\n",
        "      Si     % ['Compression/Cooling', 'Exchanger', 'Throttle'] [[1.2626999999999999, 0.40460000000000002, 1.5033000000000001], array([ 39.8,  12.8,  47.4])]\n",
        "Sum [3.1706000000000003, 100.0]\n",
        "      Wi     %\n",
        "['Ideal', 'Compression/Cooling', 'Exchanger', 'Throttle']\n",
        "[array([  53.8,  378.8,  121.4,  451. ]), array([  5.4,  37.7,  12.1,  44.9])]\n",
        "Sum [1005.0, 100.09999999999999]\n"
       ]
      }
     ],
     "prompt_number": 6
    }
   ],
   "metadata": {}
  }
 ]
}