path: root/Engineering_Thermodynamics_by_P._K._Nag/ch12.ipynb

{
 "metadata": {
  "name": "",
  "signature": "sha256:abd75e2d0144cfdf6dd59935ddb4f9b5824ec80059ba3e73a72bf780292d0a06"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 12 : Vapour Power Cycles"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.1 Page No : 475"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "# Part (a)\n",
      "P1 = 1e05;\n",
      "P2 = 10e05;\n",
      "vf = 0.001043;\n",
      "\n",
      "# Calculation and Results\n",
      "Wrev = vf*(P1-P2);\n",
      "print \"The work required in saturated liquid form is\",Wrev/1000,\"kJ/kg\"\n",
      "\n",
      "# Part (b)\n",
      "h1 = 2675.5; \n",
      "s1 = 7.3594;\n",
      "s2 = s1;\n",
      "h2 = 3195.5;\n",
      "Wrev1 = h1-h2;\n",
      "print \"The work required in saturated vapour form is\",Wrev1,\"kJ/kg\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The work required in saturated liquid form is -0.9387 kJ/kg\n",
        "The work required in saturated vapour form is -520.0 kJ/kg\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.2 Page No : 476"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 3159.3;\n",
      "s1 = 6.9917;\n",
      "h3 = 173.88; \n",
      "s3 = 0.5926; \n",
      "sfp2 = s3; \n",
      "hfp2 = h3;\n",
      "hfgp2 = 2403.1; \n",
      "sgp2 = 8.2287;\n",
      "vfp2 = 0.001008; \n",
      "sfgp2 = 7.6361;\n",
      "\n",
      "# Calculation and Results\n",
      "x2s = (s1-sfp2)/(sfgp2);\n",
      "h2s = hfp2+(x2s*hfgp2);\n",
      "\n",
      "# Part (a)\n",
      "P1 = 20e02; \n",
      "P2 = 0.08e02;\n",
      "h4s = vfp2*(P1-P2)+h3 ; \n",
      "Wp = h4s-h3;\n",
      "Wt = h1-h2s;\n",
      "Wnet = Wt-Wp;\n",
      "Q1 = h1-h4s;\n",
      "n_cycle = Wnet/Q1;\n",
      "print \"Net work per kg of steam is %.2f kJ/Kg\"%Wnet\n",
      "print \"Cycle efficiency is %.3f\"%n_cycle\n",
      "\n",
      "# Part (b)\n",
      "n_p = 0.8; \n",
      "n_t = 0.8;\n",
      "Wp_ = Wp/n_p;\n",
      "Wt_ = Wt*n_t;\n",
      "Wnet_ = Wt_-Wp_;\n",
      "P = 100*((Wnet-Wnet_)/Wnet) ;\n",
      "n_cycle_ = Wnet_/Q1;\n",
      "P_ = 100*((n_cycle-n_cycle_)/n_cycle);\n",
      "print \"Percentage reduction in net work per kg of steam is %.1f %%\"%P\n",
      "print \"Percentage reduction in cycle efficiency is %.1f %%\"%P_\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Net work per kg of steam is 969.60 kJ/Kg\n",
        "Cycle efficiency is 0.325\n",
        "Percentage reduction in net work per kg of steam is 20.1 %\n",
        "Percentage reduction in cycle efficiency is 20.1 %\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.3 Page No : 477"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "P1 = 0.08; \t\t\t# in bar\n",
      "sf = 0.5926; \n",
      "x2s = 0.85; \n",
      "sg = 8.2287; \n",
      "\n",
      "# Calculation\n",
      "s2s = sf+(x2s*(sg-sf));\n",
      "s1 = s2s;\n",
      "P2 = 16.832; \t\t\t# by steam table opposite to s1 in bar\n",
      "h1 = 3165.54;\n",
      "h2s = 173.88 + (0.85*2403.1);\n",
      "h3 = 173.88;\n",
      "vfp2 = 0.001;\n",
      "h4s = h3 + (vfp2*(P2-P1)*100);\n",
      "Q1 = h1-h4s;\n",
      "Wt = h1-h2s;\n",
      "Wp = h4s-h3;\n",
      "n_cycle = 100*((Wt-Wp)/Q1);\n",
      "Tm = (h1-h4s)/(s2s-sf);\n",
      "\n",
      "# Results\n",
      "print \"The greatest allowable steam pressure at the turbine inlet is\",P2,\"bar\"\n",
      "print \"Rankine cycle efficiency is %.2f %%\"%n_cycle\n",
      "print \"Mean temperature of heat addition is %.2f K\"%Tm\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The greatest allowable steam pressure at the turbine inlet is 16.832 bar\n",
        "Rankine cycle efficiency is 31.68 %\n",
        "Mean temperature of heat addition is 460.66 K\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.4 Page No : 478"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 3465.;\n",
      "h2s = 3065.; \n",
      "h3 = 3565.;\n",
      "h4s = 2300.; \n",
      "x4s = 0.88; \n",
      "h5 = 191.83;\n",
      "v = 0.001;\n",
      "P = 150.; \t\t\t# in bar\n",
      "\n",
      "# Calculation\n",
      "Wp = v*P*100;\n",
      "h6s = 206.83;\n",
      "Q1 = (h1-h6s)+(h3-h2s);\n",
      "Wt = (h1-h2s)+(h3-h4s);\n",
      "Wnet = Wt-Wp;\n",
      "n_cycle = 100*Wnet/Q1;\n",
      "sr = 3600/Wnet;\n",
      "\n",
      "# Results\n",
      "print \"Quality at turbine exhaust is\",0.88\n",
      "print \"Cycle efficiency is %.1f %%\"%n_cycle\n",
      "print \"steam rate is %.2f Kg/kW h\"%sr\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Quality at turbine exhaust is 0.88\n",
        "Cycle efficiency is 43.9 %\n",
        "steam rate is 2.18 Kg/kW h\n"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.5 Page No : 479"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 3230.9; \n",
      "s1 = 6.9212; \n",
      "s2 = s1; \n",
      "s3 = s1;\n",
      "h2 = 2796; \n",
      "sf = 0.6493; \n",
      "sfg = 7.5009;\n",
      "\n",
      "# Calculation\n",
      "x3 = (s3-sf)/sfg;\n",
      "h3 = 191.83 + x3*2392.8;\n",
      "h4 = 191.83; h5 = h4;\n",
      "h6 = 640.23; h7 = h6;\n",
      "m = (h6-h5)/(h2-h5);\n",
      "Wt = (h1-h2)+(1-m)*(h2-h3);\n",
      "Q1 = h1-h6;\n",
      "n_cycle = 100*Wt/Q1;\n",
      "sr = 3600/Wt;\n",
      "s7 = 1.8607; s4 = 0.6493;\n",
      "Tm = (h1-h7)/(s1-s7);\n",
      "Tm1 = (h1-h4)/(s1-s4); \t\t\t# With out regeneration\n",
      "dT = Tm-Tm1;\n",
      "Wt_ = h1-h3;\n",
      "sr_ = 3600/Wt_;\n",
      "dsr = sr-sr_;\n",
      "n_cycle_ = 100*(h1-h3)/(h1-h4);\n",
      "dn = n_cycle-n_cycle_;\n",
      "\n",
      "# Results\n",
      "print \"Efficiency of the cycle is %.2f %%\"%n_cycle\n",
      "print \"Steam rate of the cycle is %.2f kg/kW h\"%sr\n",
      "print \"Increase in temperature due to regeneration is %.1f degree centigrade\"%dT\n",
      "print \"Increase in steam rate due to regeneration is %.1f kg/kW h\"%dsr\n",
      "print \"Increase in Efficiency of the cycle due to regeneration is %.2f %%\"%dn\n",
      "\n",
      "# note: rounding error is there because of decimal points."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Efficiency of the cycle is 36.07 %\n",
        "Steam rate of the cycle is 3.85 kg/kW h\n",
        "Increase in temperature due to regeneration is 27.4 degree centigrade\n",
        "Increase in steam rate due to regeneration is 0.4 kg/kW h\n",
        "Increase in Efficiency of the cycle due to regeneration is 1.90 %\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.6 Page No : 481"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 3023.5; \n",
      "s1 = 6.7664; \n",
      "s2 = s1; \n",
      "s3 = s1; \n",
      "s4 = s1;\n",
      "t_sat_20 = 212.;\n",
      "t_sat_1 = 46.;\n",
      "\n",
      "# Calculation\n",
      "dt = t_sat_20-t_sat_1;\n",
      "n =3; \t\t\t# number of heaters\n",
      "t = dt/n;\n",
      "t1 = t_sat_20-t;\n",
      "t2 = t1-t;\n",
      "\n",
      "# 0.1 bar\n",
      "hf = 191.83; \n",
      "hfg = 2392.8; \n",
      "sf = 0.6493;\n",
      "sg = 8.1502;\n",
      "# At 100 degree\n",
      "hf100 = 419.04; \n",
      "hfg100 = 2257.0; \n",
      "sf100 = 1.3069; \n",
      "sg100 = 7.3549;\n",
      "# At 150 degree\n",
      "hf150 = 632.20; \n",
      "hfg150 = 2114.3; \n",
      "sf150 = 1.8418; \n",
      "sg150 = 6.8379;\n",
      "x2 = (s1-sf150)/4.9961;\n",
      "h2 = hf150+(x2*hfg150);\n",
      "x3 = (s1-sf100)/6.0480;\n",
      "h3 = hf100+(x3*hfg100);\n",
      "x4 = (s1-sf)/7.5010;\n",
      "h4 = hf+(x4*hfg);\n",
      "h5 = hf; h6 = h5;\n",
      "h7 = hf100; h8 = h7;\n",
      "h9 = 632.2; h10 = h9;\n",
      "m1 = (h9-h7)/(h2-h7);\n",
      "m2 = ((1-m1)*(h7-h6))/(h3-h6);\n",
      "Wt = 1*(h1-h2)+(1-m1)*(h2-h3)+(1-m1-m2)*(h3-h4);\n",
      "Q1 = h1-h9;\n",
      "Wp = 0 ; \t\t\t# Pump work is neglected\n",
      "n_cycle = 100*(Wt-Wp)/Q1;\n",
      "sr = 3600/(Wt-Wp);\n",
      "\n",
      "# Results \n",
      "print \"Net work per kg of stem is %.2f kJ/Kg\"%Wt\n",
      "print \"Cycle efficiency is %.2f %%\"%n_cycle\n",
      "print \"Stream rate is %.2f Kg/kW h\"%sr\n",
      "\n",
      "# rounding off error would be there."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Net work per kg of stem is 798.64 kJ/Kg\n",
        "Cycle efficiency is 33.40 %\n",
        "Stream rate is 4.51 Kg/kW h\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.7 Page No : 483"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "# Variables\n",
      "Ti = 2000.;\n",
      "Te = 450.;\n",
      "T0 = 300.;\n",
      "Q1_dot = 100e03; \t\t\t# in kW\n",
      "cpg = 1.1;\n",
      "\n",
      "# Calculation\n",
      "wg = Q1_dot/(cpg*(Ti-Te));\n",
      "af1 = wg*cpg*T0*((Ti/T0)-1-math.log(Ti/T0));\n",
      "af2 = wg*cpg*T0*((Te/T0)-1-math.log(Te/T0));\n",
      "afi = af1-af2;\n",
      "h1 = 2801; h3 = 169; h4 = 172.8; h2 = 1890.2;\n",
      "s1 = 6.068; s2 = s1; s3 = 0.576; s4 = s3;\n",
      "Wt = h1-h2;\n",
      "Wp = h4-h3;\n",
      "Q1 = h1-h4;\n",
      "Q2 = h2-h3;\n",
      "Wnet = Wt-Wp;\n",
      "ws = Q1_dot/2628;\n",
      "afu = 38*(h1-h4-T0*(s1-s3));\n",
      "I_dot = afi-afu;\n",
      "Wnet_dot = ws*Wnet;\n",
      "afc = ws*(h2-h3-T0*(s2-s3));\n",
      "n2 = 100*Wnet_dot/af1;\n",
      "\n",
      "# Results\n",
      "print \"The second law efficiency is %.1f %%\"%n2\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The second law efficiency is 47.3 %\n"
       ]
      }
     ],
     "prompt_number": 17
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.8 Page No : 485"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "# Variables\n",
      "h1 = 2758.;      # kJ/kg\n",
      "h2 = 1817.; \n",
      "h3 = 192.\n",
      "h4 = 200.;\n",
      "cpg = 1.1\n",
      "\n",
      "# Calculation and Results\n",
      "Wt = h1-h2; \n",
      "Wp = h4-h3;\n",
      "Q1 = h1-h4; \n",
      "Wnet = Wt-Wp;\n",
      "n1 = Wnet/Q1;\n",
      "WR = Wnet/Wt;\n",
      "Q1_ = 100.;\n",
      "PO = n1*Q1_;\n",
      "cp = 1000.;\n",
      "wg = (Q1_/(833.-450));\n",
      "EIR = 59.3\n",
      "n2 = PO/EIR ;\n",
      "print (\"Part (a)\")\n",
      "print \"n1 is %.1f %%\"%(n1*100)\n",
      "print \"n2 is %.1f %%\"%(n2*100)\n",
      "print \"Work ratio is %.3f\"%WR\n",
      "\n",
      "# Part (b)\n",
      "h1b = 3398.; \n",
      "h2b = 2130.\n",
      "h3b = 192.; \n",
      "h4b = 200.;\n",
      "Wtb = 1268.; \n",
      "Wpb = 8.; \n",
      "Q1b = 3198.;\n",
      "Wt = 1268.\n",
      "Wp = 8.\n",
      "n1b = (Wt-Wp)/Q1b;\n",
      "WRb = (Wt-Wp)/Wtb;\n",
      "EIRb = 59.3;\n",
      "Wnetb = Q1b*n1b;\n",
      "n2b = 36.5/EIRb;\n",
      "print (\"\\nPart (b)\")\n",
      "print \"n1 is %.1f %%\"%(n1b*100)\n",
      "print \"n2 is %.1f %%\"%(n2b*100)\n",
      "print \"Work ratio is %.3f\"%WRb\n",
      "\n",
      "# Part (c)\n",
      "h1c = 3398.; \n",
      "h2c = 2761.; \n",
      "h3c = 3482.; \n",
      "h4c = 2522.; \n",
      "h5c = 192.; \n",
      "h6c = 200.;\n",
      "Wt1 = 637.; \n",
      "Wt2 = 960.; \n",
      "Wtc = Wt1+Wt2; Wpc = 8.;\n",
      "Wnetc = 1589.; \n",
      "Q1c = 3198+721.;\n",
      "n1c = Wnetc/Q1c;\n",
      "WRc = Wnetc/Wtc;\n",
      "POc = Q1_*n1c;\n",
      "EIRc = 59.3;\n",
      "n2c = POc/EIRc;\n",
      "print (\"\\nPart (c)\")\n",
      "print \"n1 is %.1f %%\"%(n1c*100)\n",
      "print \"n2 is %.1f %%\"%(n2c*100)\n",
      "print \"Work ratio is %.3f\"%WRc\n",
      "\n",
      "# Part (d)\n",
      "T3 = 318.8; T1 = 568.;\n",
      "n1d = 1-(T3/T1);\n",
      "Q1d = 2758-1316.;\n",
      "Wnet = Q1d*n1d;\n",
      "Wpd = 8.; \n",
      "Wtd = 641.;\n",
      "WRd = (Wtd-Wpd)/Wtd;\n",
      "POd = Q1_*0.439;\n",
      "EIRd = (Q1_/(833-593.))*cpg*((833.-300)-300*(math.log(833./300)));\n",
      "n2d = POd/94.6;\n",
      "print (\"\\nPart (d)\")\n",
      "print \"n1 is %.1f %%\"%(n1d*100)\n",
      "print \"n2 is %.1f %%\"%(n2d*100)\n",
      "print \"Work ratio is %.3f\"%WRd\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Part (a)\n",
        "n1 is 36.5 %\n",
        "n2 is 61.5 %\n",
        "Work ratio is 0.991\n",
        "\n",
        "Part (b)\n",
        "n1 is 39.4 %\n",
        "n2 is 61.6 %\n",
        "Work ratio is 0.994\n",
        "\n",
        "Part (c)\n",
        "n1 is 40.5 %\n",
        "n2 is 68.4 %\n",
        "Work ratio is 0.995\n",
        "\n",
        "Part (d)\n",
        "n1 is 43.9 %\n",
        "n2 is 46.4 %\n",
        "Work ratio is 0.988\n"
       ]
      }
     ],
     "prompt_number": 25
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.9 Page No : 488"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "hfg = 2202.6;\n",
      "Qh = 5.83;\n",
      "ws = Qh/hfg;\n",
      "eg = 0.9; \t\t\t# efficiency of generator\n",
      "P = 1000;\n",
      "Wnet = 1000/0.9;\n",
      "nbrake = 0.8;\n",
      "\n",
      "# Calculation\n",
      "h1_2s = Wnet/(ws*nbrake); \t\t\t# h1-h2s\n",
      "n_internal = 0.85;\n",
      "h12 = n_internal*h1_2s;\n",
      "hg = 2706.3; \n",
      "h2 = hg;\n",
      "h1 = h12+h2;\n",
      "h2s = h1-h1_2s;\n",
      "hf = 503.71;\n",
      "x2s = (h2s-hf)/hfg;\n",
      "sf = 1.5276; \n",
      "sfg = 5.6020;\n",
      "s2s = sf+(x2s*sfg);\n",
      "s1 = s2s;\n",
      "P1 = 22.5; \t\t\t# in bar from Moiller chart\n",
      "t1 = 360;\n",
      "\n",
      "# Results\n",
      "print \"Temperature of the steam is\",t1,\"degree\"\n",
      "print \"Pressure of the steam is \",P1,\"bar\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Temperature of the steam is 360 degree\n",
        "Pressure of the steam is  22.5 bar\n"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.10 Page No : 489"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 3037.3; h2 = 561+(0.96*2163.8);\n",
      "s2 = 1.6718+(0.96*5.3201);\n",
      "s3s = s2;\n",
      "\n",
      "# Calculation\n",
      "x3s = (s3s-0.6493)/7.5009;\n",
      "h3s = 191.83+(x3s*2392.8);\n",
      "h23 = 0.8*(h2-h3s); \t\t\t# h2-h3\n",
      "h3 = h2-h23;\n",
      "h5 = 561.47; h4 = 191.83;\n",
      "Qh = 3500.; \t\t\t# in kJ/s\n",
      "w = Qh/(h2-h5);\n",
      "Wt = 1500.;\n",
      "ws = (Wt+w*(h2-h3))/(h1-h3); \n",
      "ws_ = 3600.*ws ; \t\t\t# in kg/h\n",
      "h6 = ((ws-w)*h4+w*h5)/ws;\n",
      "h7 = h6;\n",
      "n_boiler = 0.85;\n",
      "CV = 44000.; \t\t\t# in kJ/kg\n",
      "wf = (1.1*ws_*(h1-h7))/(n_boiler*CV);\n",
      "\n",
      "# Results\n",
      "print \"Fuel buring rate is %.2f Kg/h\"%wf\n",
      "print \"Wf = %.2f tonnes/day\"%(wf*24/1000)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Fuel buring rate is 756.64 Kg/h\n",
        "Wf = 18.16 tonnes/day\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.11 Page No : 491"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 3285.; \n",
      "h2s = 3010.\n",
      "h3 = 3280.; \n",
      "h4s = 3030.;\n",
      "\n",
      "# Calculation\n",
      "h4 = h3-0.83*(h3-h4s);\n",
      "h5s = 2225;\n",
      "h5 = h4-0.83*(h4-h5s);\n",
      "h6 = 162.7; h7 = h6;\n",
      "h8 = 762.81;\n",
      "h2 = h1-0.785*(h1-h2s);\n",
      "m = (h8-h7)/(h4-h7);\n",
      "n_cycle = ((h1-h2)+(h3-h4)+(1-m)*(h4-h5))/((h1-h8)+(h3-h2))\n",
      "\n",
      "# Results\n",
      "print \"Steam flow at turbine inlet is %.3f Kg/s\"%m\n",
      "print \"cycle efficiency is %.2f %%\"%(n_cycle*100)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Steam flow at turbine inlet is 0.206 Kg/s\n",
        "cycle efficiency is 35.92 %\n"
       ]
      }
     ],
     "prompt_number": 21
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 12.12 Page No : 493"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# From table and graph \n",
      "\n",
      "# Variables\n",
      "h1 = 2792.2;\n",
      "h4 = 122.96;\n",
      "hb = 254.88;\n",
      "hc = 29.98;\n",
      "ha = 355.98;\n",
      "hd = hc;\n",
      "h2 = 1949.27;\n",
      "\n",
      "# Calculation and Results\n",
      "m = (h1-h4)/(hb-hc); \t\t\t# Amount of mercury circulating\n",
      "Q1t = m*(ha-hd);\n",
      "W1t = m*(ha-hb) + (h1-h2);\n",
      "Nov = W1t/Q1t ;\n",
      "print \"Overall efficiency of the cycle %.2f %%\"%(Nov*100)\n",
      "\n",
      "S = 50000.; \t\t\t# Stem flow rate through turbine in kg/h\n",
      "wm = S*m;\n",
      "print \"Flow through the mercury turbine is %.2e Kg/h\"%wm\n",
      "\n",
      "Wt = W1t*S/3600;\n",
      "print \"Useful work done in binary vapour cycle is %.2e kW\"%Wt\n",
      "\n",
      "nm = 0.85; \t\t\t# Internal efficiency of mercury turbine\n",
      "ns = 0.87; \t\t\t# Internal efficiency of steam turbine\n",
      "WTm = nm*(ha-hb);\n",
      "hb_ = ha-WTm; \t\t\t# hb'\n",
      "m_ = (h1-h4)/(hb_-hc); \t\t\t# m'\n",
      "h1_ = 3037.3; \t\t\t# h'\n",
      "Q1t = m_*(ha-hd)+(h1_-h1);\n",
      "x2_ = (6.9160-0.4226)/(8.47-0.4226);\n",
      "h2_ = 121+(0.806*2432.9);\n",
      "WTst = ns*(h1_-h2_);\n",
      "WTt = m_*(ha-hb_)+WTst;\n",
      "Nov = WTt/Q1t;\n",
      "\n",
      "print \"Overall efficiency is %.1f %%\"%(Nov*100)\n",
      "\n",
      "# note : rounding off error."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Overall efficiency of the cycle 52.80 %\n",
        "Flow through the mercury turbine is 5.93e+05 Kg/h\n",
        "Useful work done in binary vapour cycle is 2.84e+04 kW\n",
        "Overall efficiency is 46.2 %\n"
       ]
      }
     ],
     "prompt_number": 4
    }
   ],
   "metadata": {}
  }
 ]
}