{
 "metadata": {
  "name": ""
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h1>Chapter 5: CONTROL VOLUME ANALYSIS USING ENERGY</h1>"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.01, page: 98</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "# Initialization  of  Variable\n",
      "P1 = 7 #bar\n",
      "T1 = 200 #deg C\n",
      "m1dot = 40 #kg/s\n",
      "P2 = 7 #bar\n",
      "T2 = 40 #degC\n",
      "A2 = 25 #cm2\n",
      "P3 = 7 #bar\n",
      "vdot = 0.06 #m3/s\n",
      "\n",
      "#calculations:\n",
      "#from Table T-3, specific vol\n",
      "v3 = 1.108E-3 #m3/kg\n",
      "m3dot = vdot/v3\n",
      "\n",
      "m2dot = m3dot - m1dot\n",
      "\n",
      "#from table T-2,\n",
      "v2 = 1.0078E-3 #m3/kg\n",
      "V2 = m2dot*v2*10000/A2\n",
      "\n",
      "#Results\n",
      "print  \"mass  flow  rate at inlet 2 is\", round(m2dot,2),\"kg/s and at exit is\",round(m3dot,2),\"kg/s\"\n",
      "print  \"velocity at inlet 2 is\", round(V2,1),\"m/s\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "mass  flow  rate at inlet 2 is 14.15 kg/s and at exit is 54.15 kg/s\n",
        "velocity at inlet 2 is 5.7 m/s\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.02, page: 99</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "%pylab inline\n",
      "\n",
      "# Initialization  of  Variable\n",
      "midot = 30 #lb/s\n",
      "A = 3 #ft2\n",
      "p = 62.4 #lb/ft2\n",
      "\n",
      "#calculations:\n",
      "#creating empty lists for plotting\n",
      "L = []\n",
      "t = []\n",
      "\n",
      "for h in range(1200):\n",
      "    t.append((h-1)/10)\n",
      "    k=(h-1)/10\n",
      "    L.append(3.33*(1 - math.e**(-9*k/(p*A))))\n",
      "\n",
      "#Results\n",
      "# plots\n",
      "fig  = plt.figure()\n",
      "ax = fig.add_subplot(1, 1, 1)\n",
      "ax.plot(t,L)\n",
      "xlabel('Time, s')\n",
      "ylabel('Height, ft')\n",
      "title('Liquid Height vs Time')\n",
      "show()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Populating the interactive namespace from numpy and matplotlib\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAEZCAYAAABvpam5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlcVPX+P/DXKJigpOICCrgkyOLC4kL3poERmbiRWy4p\nV6lL5hJ5LW6/W4m3XMisNM2sTDO7knlNKUckl3FHUtG6YrEEOiySiCSo7Of3x/k6OuzbzGcOvJ6P\nx3nMzJkz57zG6LznfD7nc45KkiQJRERE9dBKdAAiIlIeFg8iIqo3Fg8iIqo3Fg8iIqo3Fg8iIqo3\nFg8iIqo3Fg8S4vjx43BxcWmSdc2bNw/vvPNOte+3atUKv//+e5Nsq6Kvv/4ao0aNqtOyW7duxYgR\nIwySo6kEBATgq6++Eh2DFIDFgwyqd+/eOHToUKX5I0aMwK+//tok29i4cSPeeOONBn3W19cXmzdv\n1pun0Wjg4OBQp8/PnDkTBw4caNC265KlqfXv3x9WVlawsrKCmZkZLCwsdK9XrVoFtVqNWbNmGTQD\nNQ9mogNQ86ZSqaBSqUTHqJYp5TNGjkuXLumejxw5ErNmzcLcuXMNvl1qfnjkQUJU/HUfHx8PLy8v\nPPzww5g2bRqmTZuGN998E0DVzT0PNkX97W9/0y0LAKtXr0aPHj1gb2+PL774otFZ//zzTwQHB+vW\n+eabb6K8vLzKbDExMXB2dkbHjh0xf/58+Pj4VDqaePXVV2FtbY1HHnkE0dHRAIB//etfOH78OBYs\nWAArKyssWrSoUo7Ro0djw4YNevPc3d2xZ88eAMArr7wCGxsbdOjQAYMGDdIrFNWpeIGJB49+tm7d\nisceewyLFy9Gp06d4OjoiFOnTmHLli3o2bMnbGxssG3bNt1ni4qKsGTJEvTq1Qu2traYN28eCgsL\na81AysTiQcIVFxcjMDAQQUFBuHnzJqZMmYLdu3fX+Zf4g0cP0dHRWLNmDQ4ePIjExEQcPHiw1s/X\ndoWev/3tb2jTpg1SUlIQHx+PmJgYfP7555WWy8nJwZQpUxAREYHc3Fw4Ozvj9OnTet/jzJkzcHFx\nwY0bN/Daa68hODgYALB8+XKMGDECGzZsQH5+PtatW1dp/TNmzMCOHTt0rxMSEnD16lWMGTMGBw4c\nwPHjx5GUlIQ///wT3377LTp37lzrd6/4b1zxSCwuLg7u7u7Izc3F9OnTMXXqVJw/fx4pKSnYvn07\nFixYgDt37gAA/vnPfyI5ORkXL15EcnIyMjIy8O9//7vWDKRMLB4kXGxsLEpLS/Hyyy+jdevWmDRp\nEoYOHdqgde3cuRNz586Fm5sbLC0tsWzZshqXlyQJixYtQqdOnXTTuHHjdDvQ7Oxs7N+/Hx988AEs\nLCzQtWtXhIaGIjIystK61Go1BgwYgMDAQLRq1QqLFi2Cra2t3jK9evVCcHAwVCoVZs+ejaysLPzx\nxx96eaoTGBiICxcuQKvVApA76ydNmgRzc3OYm5sjPz8fly9fRnl5OZydnSttuyH69OmDoKAgqFQq\nTJ06FZmZmXjrrbdgbm4Of39/tGnTBsnJyZAkCZ999hnef/99dOzYEe3bt8frr79e5b8TNQ8sHiRc\nZmYm7Ozs9Ob16tWrQevKysrSaw7r2bNnjcurVCp89NFHuHnzpm764YcfdDvxK1euoKSkBN27d9cV\nlxdffBHXr1+v8nvY29vrzav4+sEduqWlJQCgoKBAL091rKysMGbMGN3RR2RkJGbOnAkAeOKJJ7Bg\nwQLMnz8fNjY2CAkJQX5+fo3fvS5sbGx0zy0sLAAAXbt21ZtXUFCA69ev486dOxg8eLDu32n06NHI\nyclpdAYyTSweJFz37t2RkZGhN+/KlSu65+3atdM1jQDAtWvXalzX1atXda8ffF5XD/76d3BwwEMP\nPYQbN27oisuff/6JX375pdLnevTogfT0dL31PPi6NnVppps+fTp27NiB06dPo7CwECNHjtS9t3Dh\nQpw9exYJCQlITEzE6tWr67ztxurSpQssLCyQkJCg+3fKy8vDrVu3jJaBjIvFgwyuuLgYhYWFuqms\nrEzv/b/85S8wMzPDunXrUFJSgt27d+Onn37Sve/u7o5Lly7h4sWLKCwsRHh4uN7nJUnS7fCnTp2K\nrVu34vLly7hz506tzVb3Pl+d7t2746mnnsLixYuRn5+P8vJypKSk4NixY5WWDQgIwC+//IK9e/ei\ntLQUGzZsqLHQVWRjY4OUlJQalwkICMCVK1ewdOlSTJs2TTf/7NmzOHPmDEpKSmBpaYm2bduidevW\ntW6zqe7I0KpVK7zwwgsIDQ3VHZVlZGQgJiamSdZPpofFgwwuICAAlpaWumnZsmV6HbNt2rTB7t27\nsXXrVnTu3Bk7d+7ExIkTdTu2fv364a233sKTTz4JZ2dnjBgxQu9X+oPrevrppxEaGoonnngC/fr1\ng5+fX62/6Kt6/8F527ZtQ3FxMdzc3GBtbY0pU6boisKD2+7SpQu+/fZbvPbaa+jSpQsuX76MIUOG\n4KGHHqq0bFXbefnll7Fr1y5YW1sjNDS0yqxt2rTBxIkTcejQIcyYMUM3/9atW/j73/8Oa2tr9O7d\nG126dMGrr75a4/eu7rs/+F5NeSuKiIiAo6MjHn30UXTo0AH+/v5ITEysNQMpk0rUzaAKCwvh4+OD\noqIiFBcXY8KECVi5cqXeMhqNBhMmTMAjjzwCAJg0aVKDB4ORssyZMwf29vZ4++23RUdpsPLycjg4\nOOA///kPfHx8RMchalLCBgm2bdsWR44cgaWlJUpLSzF8+HCcOHECw4cP11vOx8cHUVFRglKSKEq9\nwWVMTAyGDRsGCwsLXZ/Do48+KjgVUdMT2mx172yT4uJilJWVwdrautIySt2JUOOY0sjv+jh9+jQc\nHR3RtWtX7Nu3D3v27NE1WxE1J8KarQD5sN7LywspKSmYN28e3n33Xb33jx49iokTJ8Le3h52dnZ4\n77334ObmJigtERHdI/TIo1WrVrhw4QLS09Nx7NgxaDQavfe9vLyg1Wpx8eJFLFy4EIGBgWKCEhGR\nHqFHHg96++23YWFhgSVLllS7TJ8+fXDu3LlKzVuOjo61nuJIRET39e3bF8nJyQ3+vLAjj5ycHOTl\n5QEA7t69ix9//BGenp56y2RnZ+v6POLi4iBJUpX9IikpKbpz/ZU2LV26VHgG5hefg/mVOSk5f2N/\ncAs72yorKwtBQUEoLy9HeXk5Zs2aBT8/P2zatAkAEBISgl27dmHjxo0wMzODpaUlr5NDRGQihBWP\ngQMH4vz585Xmh4SE6J7Pnz8f8+fPN2YsIiKqA44wF8zX11d0hEZhfrGYXyyl528Mk+kwbwyVSoVm\n8DWIiIymsftNHnkQEVG9sXgQEVG9sXgQEVG9CTvbioioOZEkoKwMKCrSn4qLa55XUgKUlt5/NMTz\nsrL7j/emxmLxIKJmS5KAwkLg1i0gPx8oKADu3NGf7t6t3+vCwuqLgkoFPPSQ/tSmTeV5D843N5cn\nM7P7j1U9b9u26vdrm2dmBrRuXXlyd2/cvy3PtiIik1ReLu/0c3OBmzfvP968CeTlycUgP/9+Yaju\nuZkZYGUFPPww0L490K4dYGEBWFren+r62sJCnqorCnW4eaPJaOx+k8WDiAyuvFze4f/xR9VTbq5+\nccjNlXf+7dsDnTrJk7X1/ccOHeRiYGV1f6rutbm56G9vmlg8wOJBJEppKXDtGpCRAaSny4/3pgeL\nQ06O/Iu/W7fKU9euQOfO+sWhUye5QJixYd1gWDzA4kFkKHl5QGrq/SktTb9IXL8OdOkC2NnJk729\n/NijB2Bre79AdOkiN+uQ6WDxAIsHUUOVl8vF4LffgKQk4Pff9YtFaSnQpw/wyCPyY69egIPD/WJh\na8tmIaVi8QCLB1Ft7tyRC8Svv+o/JiXJfQMuLoCT0/0icW/q3Fk+g4iaHxYPsHgQ3SNJwJUrwM8/\nAxcvyo8//wxcvQo4OspFwtn5/qOzs1w8qOVh8QCLB7VMkgQkJwNxcfJ0/rxcKNq1k8/hHzTo/qOz\nM5uXSB+LB1g8qGX44w/gzJn7xeKnn+RTUYcNkycvL7lQdO0qOikpAYsHWDyoedJqgWPH7k9ZWYC3\n9/1iMXSo3GFN1BAsHmDxoObh2jUgJgY4fBg4elS+lMbjj9+fBg1S1ghmMm0sHmDxIGUqKgJOnAAO\nHJCLxpUrgJ8f8OSTgI+P3KnNM53IUBRbPAoLC+Hj44OioiIUFxdjwoQJWLlyZaXlFi1ahP3798PS\n0hJbt26Fp6dnpWVYPEgpsrOB778H9u6Vjy769weeegoYNUpuiuKIajKWxu43hf2ptm3bFkeOHIGl\npSVKS0sxfPhwnDhxAsOHD9cto1arkZycjKSkJJw5cwbz5s1DbGysqMhEDZKSAuzZI0+//AI8/TQw\nYwbw5ZfypTiIlEjo7xxLS0sAQHFxMcrKymBd4f+kqKgoBAUFAQC8vb2Rl5eH7Oxs2NjYGD0rUX1c\nuQLs2CFP2dnAhAnA//t/wBNP8DId1DwIvZNgeXk5PDw8YGNjg5EjR8LNzU3v/YyMDDg4OOhe29vb\nIz093dgxieokJwfYuBEYMQIYPFi+DtS6dfI1oDZtAkaPZuGg5kNo8WjVqhUuXLiA9PR0HDt2DBqN\nptIyFdvkVOxBJBNSWir3YUyYII/gPn4cCAsDMjOBTz6RO755hhQ1RybRPdehQweMGTMGZ8+eha+v\nr26+nZ0dtFqt7nV6ejrs7OyqXEd4eLjuua+vr956iJpaWhqweTOwZYt8ocAXXgC+/lq+/wSRKdJo\nNFX+QG8oYWdb5eTkwMzMDB07dsTdu3cxatQoLF26FH5+frpl1Go11q9fD7VajdjYWISGhlbZYc6z\nrcgYysuB6Gi5Keqnn4DnngOefx4YOFB0MqL6U+zZVllZWQgKCkJ5eTnKy8sxa9Ys+Pn5YdOmTQCA\nkJAQBAQEQK1Ww9HREe3atcOWLVtExaUW7PZtYNs2YO1a+RakoaHAd9/Jz4laKg4SJKpGVpZcMDZv\nBoYPl4vG449z4B41D43dbwrtMCcyRVotsGCBPIDvzh0gNlY+0vDxYeEguofFg+j//P673PHt4SFf\n1vzyZbl/o29f0cmITA+LB7V4GRly0Rg2TL5KbWIiEBEBcCwqUfVYPKjFyssDXn9dvlqttbV8S9a3\n35ZvvUpENWPxoBanqAhYswbo1w+4fl2+XWtEBNCpk+hkRMphEoMEiYxFrZbPmnJ2BjQaoMIVcYio\njlg8qEVISZGLxm+/yaffjh4tOhGRsrHZipq1oiJg6VL59q2PPSZfEp2Fg6jxeORBzdbp00BwsNy3\nceECYG8vOhFR88HiQc1OQQHwxhvAN9/I4zQmT+bgPqKmxmYralaOHpVPvb15E/jf/4ApU1g4iAyB\nRx7ULBQXy30bX34JfPYZMGaM6EREzRuLByneb78BM2cC3bvLfRvduolORNT8sdmKFO3eFW+Dg4Go\nKBYOImPhkQcp0t27wPz5wJkzwLFjgKur6ERELQuPPEhxkpOBv/wFKCyUiwcLB5HxsXiQokRFAX/9\nK+8ZTiQam61IESQJWLUK+PhjuYA8+qjoREQtG4sHmbyiIuDvfwcuXZLv6mdnJzoREQlrttJqtRg5\nciT69++PAQMGYN26dZWW0Wg06NChAzw9PeHp6Yl33nlHQFIS6fp1wM8PuH1b7hhn4SAyDcKOPMzN\nzfHBBx/Aw8MDBQUFGDx4MPz9/eFaoffTx8cHUVFRglKSSElJwNNPA9OmyTdpasUeOiKTIex/R1tb\nW3h4eAAA2rdvD1dXV2RmZlZaTpIkY0cjE3D+PODjA/zzn8Dy5SwcRKbGJP6XTEtLQ3x8PLy9vfXm\nq1QqnDp1Cu7u7ggICEBCQoKghGRMhw/LRxwbNshnVRGR6RHeYV5QUIDJkydj7dq1aF/hvEsvLy9o\ntVpYWlpi//79CAwMRGJiYpXrCQ8P1z339fWFr6+vAVOToezaBbz0EvDtt/KRBxE1DY1GA41G02Tr\nU0kC24VKSkowduxYjB49GqGhobUu36dPH5w7dw7W1tZ681UqFZu3moGvvgLCwuRbxf5fiyYRGUhj\n95vCmq0kSUJwcDDc3NyqLRzZ2dm6LxcXFwdJkioVDmoevvxS7t84eJCFg0gJhDVbnTx5Etu3b8eg\nQYPg6ekJAFixYgWuXr0KAAgJCcGuXbuwceNGmJmZwdLSEpGRkaLikgF98QXw1lvAoUOAi4voNERU\nF0KbrZoKm62U6/PPgWXL5MLRr5/oNEQtR2P3m8I7zKnl2r5dLhyHDwNOTqLTEFF98MiDhIiKki85\ncvgw4OYmOg1Ry8MjD1KcI0eA558H9u1j4SBSKhYPMqq4OGDqVHkcx9ChotMQUUOZxAhzahmSk4EJ\nE+Rbx3IMJ5GysXiQUdy4AQQEAOHhwPjxotMQUWOxw5wMrrAQ8PeX7wAYESE6DREBjd9vsniQQZWX\nA889B5SWApGRvDoukang2VZk0sLDgbQ0eRAgCwdR88HiQQbz3//K16z66SfAwkJ0GiJqSmy2IoP4\n3/+AkSOB6Ghg8GDRaYioIsVeVZear5s3gcBA4P33WTiImiseeVCTKisDxoyRr4774Yei0xBRdXjk\nQSYlPBwoKgJWrxadhIgMiR3m1GQOHpTvzXH+PGBuLjoNERkSiwc1iWvXgNmz5VvJ2tiITkNEhsZm\nK2q0sjJ5IODzzwN+fqLTEJExsHhQo61aBZSUyLeSJaKWgc1W1CgnTwIffQScOweY8a+JqMUQduSh\n1WoxcuRI9O/fHwMGDMC6deuqXG7RokVwcnKCu7s74uPjjZySalJQIPdzbNoE2NmJTkNExiTst6K5\nuTk++OADeHh4oKCgAIMHD4a/vz9cXV11y6jVaiQnJyMpKQlnzpzBvHnzEBsbKyoyVfCPfwA+PvI9\nOoioZRFWPGxtbWFrawsAaN++PVxdXZGZmalXPKKiohAUFAQA8Pb2Rl5eHrKzs2HD03mEU6uBAweA\nn38WnYSIRDCJDvO0tDTEx8fD29tbb35GRgYcHBx0r+3t7ZGenm7seFRBTg7wwgvA1q3Aww+LTkNE\nIgjv4iwoKMDkyZOxdu1atG/fvtL7FYfPq1SqKtcTHh6ue+7r6wtf3ufUYF56CXj2Wd5KlkhJNBoN\nNBpNk61P6LWtSkpKMHbsWIwePRqhoaGV3n/xxRfh6+uLadOmAQBcXFxw9OjRSs1WvLaV8Xz3HRAW\nBly8yMusEymZYq9tJUkSgoOD4ebmVmXhAIDx48dj27ZtAIDY2Fh07NiR/R0C5eUBCxcCn33GwkHU\n0gk78jhx4gQef/xxDBo0SNcUtWLFCly9ehUAEBISAgBYsGABoqOj0a5dO2zZsgVeXl6V1sUjD+N4\n8UVAkuRTc4lI2XgPc7B4GMPx48D06fJNnjp2FJ2GiBpLsc1WpByFhfLZVR99xMJBRLJai8fatWvr\nNI+ar5UrATc34JlnRCchIlNRa7OVp6dnpcuCeHh44MKFCwYNVh9stjKclBTA2xuIjwceGHJDRArX\n2P1mteM8duzYgf/85z9ITU3FuHHjdPPz8/PRuXPnBm+QlCU0FFiyhIWDiPRVWzy8vb3RvXt35OTk\nYMmSJboKZWVlBXd3d6MFJHF++AFITAR27RKdhIhMTbXFY8qUKTh37hwsLCzg4+NjzExkAgoL5aOO\nDRuAhx4SnYaITE21xaOsrAzLly9HYmIi3n//fb22MZVKhcWLFxslIImxejUwaBAwapToJERkiqot\nHpGRkdizZw/KysqQn59vzEwkmFYLfPghcP686CREZKpqPdtKrVYjICDAWHkahGdbNa2gILmD/J13\nRCchIkPhCHOweDSl8+eBMWPkjnIrK9FpiMhQOMKcmowkyaflLl3KwkFENWPxIB21GsjKAp5/XnQS\nIjJ19S4eP/30EzIzMw2RhQQqLQVefVU+y8pM+C3CiMjU1bt4fPTRRxgzZgyeffZZQ+QhQTZvBmxt\n5f4OIqLaNLjD/NatW3jYRG5gzQ7zxrl7F3B0BPbuBYYMEZ2GiIzB4B3mfn5+Vc4zlcJBjbdxo3zx\nQxYOIqqralu37969izt37uD69evIzc3Vzb916xYyMjKMEo4Mr6AAePdd4McfRSchIiWptnhs2rQJ\na9euRWZmJgYPHqybb2VlhQULFhglHBne2rXAE08AAweKTkJESlJrn8e6deuwaNEiY+VpEPZ5NMzN\nm0C/fsDJk/IjEbUcRhlhfurUKaSlpaG0tFQ3b/bs2Q3e6D1z587Fvn370K1bN/zyyy+V3tdoNJgw\nYQIeeeQRAMCkSZPwxhtvVFqOxaNh3nwTyMyUz7QiopbFYDeDuue5557D77//Dg8PD7Ru3Vo3vymK\nx5w5c7Bw4cIa1+Xj44OoqKhGb4v05eQAH38MnDsnOgkRKVGtxePcuXNISEiASqVq8o2PGDECaWlp\nNS7DIwrD+OADYMoUoHdv0UmISIlqPVV3wIAByMrKMkaWSlQqFU6dOgV3d3cEBAQgISFBSI7mJi8P\n2LQJCAsTnYSIlKraI4979y0vKCiAm5sbhg0bhof+75ZyKpXKKE1JXl5e0Gq1sLS0xP79+xEYGIjE\nxMQqlw0PD9c99/X1ha+vr8HzKdWGDUBAANCnj+gkRGQsGo0GGo2mydZXbYd5bRtpqp1zWloaxo0b\nV2WHeUV9+vTBuXPnYG1trTefHeZ1d/s28MgjgEYDuLqKTkNEohisw9wUfrlnZ2ejW7duUKlUiIuL\ngyRJlQoH1c9nnwHDh7NwEFHj1NphblXFjR06dOiAoUOHYs2aNbrTaBti+vTpOHr0KHJycuDg4IBl\ny5ahpKQEABASEoJdu3Zh48aNMDMzg6WlJSIjIxu8LQKKioD33pOvYUVE1Bi1jvN444034ODggOnT\npwOQ722ekpICT09PfPLJJ03ahtZQbLaqm08/Bb77Dti/X3QSIhLN4IMEBw0ahJ9//llvnoeHBy5c\nuAB3d3dcvHixwRtvKiwetSsrA5ydgS1bgBEjRKchItEMflVdS0tLfPPNNygvL0d5eTl27tyJtm3b\n6jZOyhAVBXTpIvd3EBE1Vq3F4+uvv8ZXX32Fbt26oVu3bti2bRu2b9+Ou3fvYv369cbISE3g/feB\nxYsB1nsiagoNvhmUKWGzVc3i4oCpU4HkZN5ilohkBjtVNyIiAmFhYVi4cGGVG123bl2DN0rG9f77\nwMsvs3AQUdOpdnfi5uYGAHr38rhXqdjXoRxXrsg3evr0U9FJiKg5qXOz1e3bt9GuXTtD52kQNltV\nb8kS+fG998TmICLTYvCzrU6dOgU3Nze4uLgAAC5cuICXXnqpwRsk47l1Sz41t4qWRyKiRqm1eISG\nhiI6OhpdunQBII/xOHr0qMGDUeNt2QI8+STQq5foJETU3NSpC7Vnz576H2LPq8krL5evnrtli+gk\nRNQc1VoFevbsiZMnTwIAiouLsW7dOrjyqnom7/BhwMIC+OtfRSchouao1marjRs3YsOGDcjIyICd\nnR3i4+OxYcMGY2SjRtiwAXjpJQ4KJCLD4CDBZkirBTw85NN027cXnYaITJHBBgk+ODiw4kY4SNC0\nbdoEzJzJwkFEhlNt8Rg8eLCuaCxduhT//ve/dQWEgwRNV3Ex8PnnwJEjopMQUXNWp2YrT09PxMfH\nGyNPg7DZ6r4dO+TiceiQ6CREZMoMPkiQlOXjj+WOciIiQ2LxaEYSEoCUFGD8eNFJiKi5q7bPo337\n9rq+jbt37+rdy1ylUuHWrVuGT0f1snkzEBQEmJuLTkJEzV21Rx4FBQXIz89Hfn4+SktLdc/z8/Ob\npHDMnTsXNjY2GDhwYLXLLFq0CE5OTnB3dzfpPhdTUFwMfPUVMHeu6CRE1BIIa7aaM2cOoqOjq31f\nrVYjOTkZSUlJ+PTTTzFv3jwjplOe778HXF0BJyfRSYioJRBWPEaMGIFOnTpV+35UVBSCgoIAAN7e\n3sjLy0N2drax4inOF18AwcGiUxBRS2GyHeYZGRlwcHDQvba3t0d6errARKYrIwM4fRqYPFl0EiJq\nKUz68rgVz0GuaXBieHi47rmvry98fX0NlMr0bN0q36Pc0lJ0EiIyVRqNBhqNpsnWZ7LFw87ODlqt\nVvc6PT0ddnZ21S7/YPFoScrL5SaryEjRSYjIlFX8Ub1s2bJGrc9km63Gjx+Pbdu2AQBiY2PRsWNH\n2NjYCE5leo4eBdq1A4YMEZ2EiFoSYUce06dPx9GjR5GTkwMHBwcsW7YMJSUlAICQkBAEBARArVbD\n0dER7dq1wxbe1ahKW7bIHeW83BgRGRMvya5gt28DdnZAYiLQrZvoNESkJLy2VQu2Zw/w2GMsHERk\nfCweCrZ9O/Dcc6JTEFFLxGYrhcrOBlxc5DEePEWXiOqLzVYt1I4dwIQJLBxEJAaLh0KxyYqIRGLx\nUKDLl4GsLGDkSNFJiKilYvFQoO3bgRkzgNatRSchopaKHeYKI0lAnz7A3r2Au7voNESkVOwwb2HO\nnJE7yQcNEp2EiFoyFg+F2blTvoIuL0dCRCKx2UpBysuBXr2AAwcANzfRaYhIydhs1YLExgIdO7Jw\nEJF4LB4K8s03cpMVEZFobLZSiPJywN4eOHIEcHYWnYaIlI7NVi3EyZPy1XNZOIjIFLB4KASbrIjI\nlLDZSgHKyuSbPh0/Djg5iU5DRM0Bm61agOPHgR49WDiIyHSweCjAd98BkyaJTkFEdJ/Q4hEdHQ0X\nFxc4OTkhIiKi0vsajQYdOnSAp6cnPD098c477whIKZYkybebDQwUnYSI6D4zURsuKyvDggULcPDg\nQdjZ2WHo0KEYP348XF1d9Zbz8fFBVFSUoJTinT8PPPQQBwYSkWkRduQRFxcHR0dH9O7dG+bm5pg2\nbRr27t1babnm3BFeF3v2AM88w2tZEZFpEVY8MjIy4ODgoHttb2+PjIwMvWVUKhVOnToFd3d3BAQE\nICEhwdgxhfvuOzZZEZHpEdZsparDT2kvLy9otVpYWlpi//79CAwMRGJiYpXLhoeH6577+vrC19e3\niZKKk5QJHWrOAAANT0lEQVQE5OYC3t6ikxCR0mk0Gmg0miZbn7BxHrGxsQgPD0d0dDQAYOXKlWjV\nqhXCwsKq/UyfPn1w7tw5WFtb681vruM83n0XSE0FNm4UnYSImhvFjvMYMmQIkpKSkJaWhuLiYnzz\nzTcYP3683jLZ2dm6LxcXFwdJkioVjuaMZ1kRkakS1mxlZmaG9evXY9SoUSgrK0NwcDBcXV2xadMm\nAEBISAh27dqFjRs3wszMDJaWloiMjBQV1+iysoBffwVGjhSdhIioMl6exER98ok8svzrr0UnIaLm\nSLHNVlSzqChgwgTRKYiIqsYjDxN0+zZgawukpwMdOohOQ0TNEY88mqHDh4GhQ1k4iMh0sXiYoB9+\nAMaMEZ2CiKh6bLYyMZIEODgAhw7xroFEZDhstmpmLl4E2rYF+vUTnYSIqHosHiZm3z5g7FheCJGI\nTBuLh4lhfwcRKQH7PEzI9euAoyPwxx/yPTyIiAyFfR7NyP79gJ8fCwcRmT4WDxNyr7+DiMjUsdnK\nRJSWAl27AgkJQPfuotMQUXPHZqtm4swZoHdvFg4iUgYWDxMREwOMGiU6BRFR3bB4mIgDB4CnnhKd\ngoiobtjnYQJyc+Umq+vXeaYVERkH+zyagUOHgBEjWDiISDlYPExATAybrIhIWYQWj+joaLi4uMDJ\nyQkRERFVLrNo0SI4OTnB3d0d8fHxRk5oeJIk93ews5yIlERY8SgrK8OCBQsQHR2NhIQE7NixA5cv\nX9ZbRq1WIzk5GUlJSfj0008xb948QWkN57ff5Edefp2IlERY8YiLi4OjoyN69+4Nc3NzTJs2DXv3\n7tVbJioqCkFBQQAAb29v5OXlITs7W0Rcg7l31MGr6BKRkggrHhkZGXBwcNC9tre3R0ZGRq3LpKen\nGy2jMbC/g4iUSFjxUNXxp3bFU8nq+jklKCoCjh+XL4ZIRKQkZqI2bGdnB61Wq3ut1Wphb29f4zLp\n6emws7Orcn3h4eG6576+vvD19W3SvIZw8iTg5gZYW4tOQkTNnUajgUajabL1CRskWFpaCmdnZxw6\ndAg9evTAsGHDsGPHDri6uuqWUavVWL9+PdRqNWJjYxEaGorY2NhK61LqIME//gCuXgWGDBGdhIha\nmsbuN4UdeZiZmWH9+vUYNWoUysrKEBwcDFdXV2zatAkAEBISgoCAAKjVajg6OqJdu3bYsmWLqLgG\n0a2bPBERKQ0vT0JE1ALx8iRERGR0LB5ERFRvLB5ERFRvLB5ERFRvLB5ERFRvLB5ERFRvLB5ERFRv\nLB5ERFRvLB5ERFRvLB5ERFRvLB5ERFRvLB5ERFRvLB5ERFRvLB5ERFRvLB5ERFRvLB5ERFRvLB5E\nRFRvLB5ERFRvLB5ERFRvZiI2mpubi2effRZXrlxB7969sXPnTnTs2LHScr1798bDDz+M1q1bw9zc\nHHFxcQLSEhFRRUKOPFatWgV/f38kJibCz88Pq1atqnI5lUoFjUaD+Pj4Zls4NBqN6AiNwvxiMb9Y\nSs/fGEKKR1RUFIKCggAAQUFB2LNnT7XLSpJkrFhCKP2Pj/nFYn6xlJ6/MYQUj+zsbNjY2AAAbGxs\nkJ2dXeVyKpUKTz75JIYMGYLPPvvMmBGJiKgGBuvz8Pf3x7Vr1yrNX758ud5rlUoFlUpV5TpOnjyJ\n7t274/r16/D394eLiwtGjBhhkLxERFQPkgDOzs5SVlaWJEmSlJmZKTk7O9f6mfDwcOm9996r8r2+\nfftKADhx4sSJUx2nvn37Nmo/LuRsq/Hjx+PLL79EWFgYvvzySwQGBlZa5s6dOygrK4OVlRVu376N\nmJgYLF26tMr1JScnGzoyERE9QCVJxu+Rzs3NxdSpU3H16lW9U3UzMzPxwgsvYN++ffj9998xceJE\nAEBpaSlmzpyJ119/3dhRiYioCkKKBxERKZtiR5i/+uqrcHV1hbu7OyZOnIg///xT997KlSvh5OQE\nFxcXxMTECExZs+joaLi4uMDJyQkRERGi49RKq9Vi5MiR6N+/PwYMGIB169YBkI8k/f390a9fPzz1\n1FPIy8sTnLR6ZWVl8PT0xLhx4wAoK3teXh4mT54MV1dXuLm54cyZM4rKv3LlSvTv3x8DBw7EjBkz\nUFRUZNL5586dCxsbGwwcOFA3r6a8prbfqSp/k+43G9VjIlBMTIxUVlYmSZIkhYWFSWFhYZIkSdKl\nS5ckd3d3qbi4WEpNTZX69u2rW86UlJaWSn379pVSU1Ol4uJiyd3dXUpISBAdq0ZZWVlSfHy8JEmS\nlJ+fL/Xr109KSEiQXn31VSkiIkKSJElatWqV7r+FKVqzZo00Y8YMady4cZIkSYrKPnv2bGnz5s2S\nJElSSUmJlJeXp5j8qampUp8+faTCwkJJkiRp6tSp0tatW006/7Fjx6Tz589LAwYM0M2rLq8p7neq\nyt+U+03FFo8H7d69W5o5c6YkSZK0YsUKadWqVbr3Ro0aJZ0+fVpUtGqdOnVKGjVqlO71ypUrpZUr\nVwpMVH8TJkyQfvzxR8nZ2Vm6du2aJElyganL2XMiaLVayc/PTzp8+LA0duxYSZIkxWTPy8uT+vTp\nU2m+UvLfuHFD6tevn5SbmyuVlJRIY8eOlWJiYkw+f2pqqt7Ot7q8prrfqZj/QY3dbyq22epBX3zx\nBQICAgAAmZmZsLe3171nb2+PjIwMUdGqlZGRAQcHB91rU81ZnbS0NMTHx8Pb27vOgz5Fe+WVV7B6\n9Wq0anX/z14p2VNTU9G1a1fMmTMHXl5eeOGFF3D79m3F5Le2tsY//vEP9OzZEz169EDHjh3h7++v\nmPz3VJdXKfudBzV2v2nSxcPf3x8DBw6sNH3//fe6ZZYvX442bdpgxowZ1a6nukGIIpliproqKCjA\npEmTsHbtWlhZWem9V9OgT5F++OEHdOvWDZ6entVe8sZUswPyGYfnz5/HSy+9hPPnz6Ndu3aVrgln\nyvlTUlLw4YcfIi0tDZmZmSgoKMD27dv1ljHl/FWpLa8pf5em2G8KGedRVz/++GON72/duhVqtRqH\nDh3SzbOzs4NWq9W9Tk9Ph52dncEyNlTFnFqtVq/ym6qSkhJMmjQJs2bN0o3PsbGxwbVr12Bra4us\nrCx069ZNcMrKTp06haioKKjVahQWFuLWrVuYNWuWIrID8i9Be3t7DB06FAAwefJkrFy5Era2torI\nf/bsWfz1r39F586dAQATJ07E6dOnFZP/nur+XpSy3wGabr9p0kceNYmOjsbq1auxd+9etG3bVjd/\n/PjxiIyMRHFxMVJTU5GUlIRhw4YJTFq1IUOGICkpCWlpaSguLsY333yD8ePHi45VI0mSEBwcDDc3\nN4SGhurm3xv0CaDaQZ+irVixAlqtFqmpqYiMjMQTTzyBr776ShHZAcDW1hYODg5ITEwEABw8eBD9\n+/fHuHHjFJHfxcUFsbGxuHv3LiRJwsGDB+Hm5qaY/PdU9/eilP1Ok+43m6hfxugcHR2lnj17Sh4e\nHpKHh4c0b9483XvLly+X+vbtKzk7O0vR0dECU9ZMrVZL/fr1k/r27SutWLFCdJxaHT9+XFKpVJK7\nu7vu333//v3SjRs3JD8/P8nJyUny9/eXbt68KTpqjTQaje5sKyVlv3DhgjRkyBBp0KBB0jPPPCPl\n5eUpKn9ERITk5uYmDRgwQJo9e7ZUXFxs0vmnTZsmde/eXTI3N5fs7e2lL774osa8prbfqZh/8+bN\nTbrf5CBBIiKqN8U2WxERkTgsHkREVG8sHkREVG8sHkREVG8sHkREVG8sHkREVG8sHkRVuHHjBjw9\nPeHp6Ynu3bvD3t4enp6esLKywoIFC0THIxKO4zyIarFs2TJYWVlh8eLFoqMQmQweeRDVwb3fWBqN\nRncjqfDwcAQFBeHxxx9H7969sXv3bixZsgSDBg3C6NGjUVpaCgA4d+4cfH19MWTIEDz99NO4du1a\njds6evSo7qjHy8sLBQUFhv1yRA3A4kHUCKmpqThy5AiioqLw3HPPwd/fHz///DMsLCywb98+lJSU\nYOHChfjvf/+Ls2fPYs6cOfjXv/5V4zrXrFmDjz/+GPHx8Thx4gQsLCyM9G2I6s6kr6pLZMpUKhVG\njx6N1q1bY8CAASgvL8eoUaMAAAMHDkRaWhoSExNx6dIlPPnkkwDk2+D26NGjxvU+9thjeOWVVzBz\n5kxMnDjRZK/OSi0biwdRI7Rp0wYA0KpVK5ibm+vmt2rVCqWlpZAkCf3798epU6fqvM6wsDCMHTsW\n+/btw2OPPYYDBw7A2dm5ybMTNQabrYgaqC7nmjg7O+P69euIjY0FIN8PJSEhAQCwfv16bNiwodJn\nUlJS0L9/f7z22msYOnQofvvtt6YNTtQEWDyI6uDeXdUevHtcxTvJVbzzmkqlgrm5OXbt2oWwsDB4\neHjA09MTp0+fBgD8+uuv6NKlS6VtrV27FgMHDoS7uzvatGmD0aNHG+prETUYT9UlEmTcuHH47rvv\nYGbG1mNSHhYPIiKqNzZbERFRvbF4EBFRvbF4EBFRvbF4EBFRvbF4EBFRvbF4EBFRvbF4EBFRvf1/\n2tzBm6RkOo8AAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0x9808c88>"
       ]
      }
     ],
     "prompt_number": 12
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.03, page: 105</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "# Initialization  of  Variable\n",
      "P1 = 40 #bar\n",
      "T1 = 400 #degC\n",
      "V1 = 10 #m/s\n",
      "P2 = 15 #bar\n",
      "V2 = 665 #m/s\n",
      "mdot = 2 #kg/s\n",
      "\n",
      "#calculations:\n",
      "#specific potential energy at 1 is\n",
      "h1 = 3213.6 #kJ/kg\n",
      "#at 2\n",
      "h2 = h1 + (V1**2 - V2**2)/2\n",
      "\n",
      "#From Table T-4\n",
      "v2 = 0.1627 #m3/kg\n",
      "\n",
      "#exit area:\n",
      "A2 = mdot*v2/V2\n",
      "\n",
      "#Results\n",
      "print  \"Exit Area is\", round(A2,6),\"m2\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Exit Area is 0.000489 m2\n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.04, page: 107</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "# Initialization  of  Variable\n",
      "mdot = 4600 #kg/h\n",
      "Wdot = 1000 #kW\n",
      "P1 = 60 #bar\n",
      "T1 = 400 #degC\n",
      "V1 = 40 #m/s\n",
      "P2 = 0.1 #bar\n",
      "x2 = 0.9\n",
      "V2 = 50 #m/s\n",
      "\n",
      "#calculations:\n",
      "#from table T-4, specific enthalpy\n",
      "h1 = 3177.2 #kJ/kg\n",
      "#From table T-3,\n",
      "hf2 = 191.83 #kJ/kg\n",
      "hg2 = 2584.63 #kJ/kg\n",
      "\n",
      "h2 = hf2 + x2*(hg2-hf2)\n",
      "\n",
      "#rate of heat transfer\n",
      "Qdot = Wdot + mdot*(h2 - h1 + (V2**2 - V1**2)/2000)/3600\n",
      "\n",
      "#Results\n",
      "print  \"rate of heat transfer is\", round(Qdot,1),\"kW\"\n",
      "#answer wrong in book\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "rate of heat transfer is -62.3 kW\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.05, page: 109</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "# Initialization  of  Variable\n",
      "P1 = 1 #bar\n",
      "T1 = 290 #K\n",
      "V1 = 6 #m/s\n",
      "A1 = 0.1 #m2\n",
      "P2 = 7 #bar\n",
      "T2 = 450 #K\n",
      "V2 = 2 #m/s\n",
      "Qdot = -180 #kJ/min\n",
      "R = 8314/28.97 #J/kg-K\n",
      "\n",
      "#calculations:\n",
      "#specific volume at 1\n",
      "v1 = R*T1/P1\n",
      "#mass flow rate\n",
      "mdot = A1*V1*1E5/v1\n",
      "\n",
      "#from table T-9\n",
      "h1 = 290.16 #kJ/kg\n",
      "h2 = 451.8 #kJ/kg\n",
      "#power required by the compressor\n",
      "Wdot = Qdot/60 + mdot*((h1-h2) + (V1**2 - V2**2)/2000)\n",
      "\n",
      "#Results\n",
      "print  \"power required by the compressor is\", round(Wdot,1),\"kW\" \n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "power required by the compressor is -119.5 kW\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.06, page: 110</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "# Initialization  of  Variable\n",
      "T1 = 20 #degC\n",
      "P1 = 1 #atm\n",
      "Vdot = 0.1 #lt/s\n",
      "D1 = 2.5 #cm\n",
      "T2 = 23 #degC\n",
      "V2 = 50 #m/s\n",
      "z1 = 0 #m\n",
      "z2 = 5 #m\n",
      "x2 = 0.1\n",
      "c = 4.18 #kJ/kg-K\n",
      "g = 9.81 #m/s2\n",
      "\n",
      "#cqlculations:\n",
      "#from Table T-2, Specific volume at 1\n",
      "v = 1.0018E-3 #m3/kg\n",
      "#mass flow rate\n",
      "mdot = Vdot*1E-3/v\n",
      "#inlet Area\n",
      "A1 = math.pi*(D1**2)/4\n",
      "#inlet velocity\n",
      "V1 = Vdot/A1\n",
      "\n",
      "# power\n",
      "Wdot = (mdot/(1 - x2))*(c*(T1-T2) + (V1**2 - V2**2)/2000 + g*(z1 - z2)/1000)\n",
      "\n",
      "#Results\n",
      "print  \"power input is\", round(Wdot,2),\"kW\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "power input is -1.53 kW\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.07, page: 111</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "# Initialization  of  Variable\n",
      "P1 = 0.1 #bar\n",
      "x1 = 0.95\n",
      "P2 = 0.1 #bar\n",
      "T2 = 45 #deg C\n",
      "T3 = 20 #degC\n",
      "T4 = 35 #degC\n",
      "\n",
      "#calculations:\n",
      "#assumption\n",
      "Wdot = 0\n",
      "#from table T-3,\n",
      "hf1 = 191.83 #kJ/kg\n",
      "hf2 = 188.45 #kJ/kg\n",
      "hg1 = 2584.7 #kJ/kg\n",
      "h4_3 = 62.7 #kJ/kg\n",
      "#specific enthalpy\n",
      "h1 = hf1 + x1*(hg1-hf1)\n",
      "h2 = hf2\n",
      "\n",
      "#ratio of mdots\n",
      "m3dot_m1dot = (h1 - h2)/h4_3\n",
      "\n",
      "#rate of energy transfer from the condensing steam to the cooling water\n",
      "Qdotm = h2 - h1\n",
      "\n",
      "#Results\n",
      "print  \"a) ratio of the mass flow rate of the cooling water to the mass flow rate of the condensing stream is\", round(m3dot_m1dot,1)\n",
      "print  \"b) rate of energy transfer from the condensing steam to the cooling water is\",round(Qdotm,1),\"kJ/kg\" \n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "a) ratio of the mass flow rate of the cooling water to the mass flow rate of the condensing stream is 36.3\n",
        "b) rate of energy transfer from the condensing steam to the cooling water is -2276.6 kJ/kg\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.08, page: 113</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "# Initialization  of  Variable\n",
      "T1 = 20 #degC\n",
      "P1 = 1 #bar\n",
      "V1 = 1.3 #m/s\n",
      "T2 = 32 #degC\n",
      "Pc = -80 #W\n",
      "Pf = -18 #W\n",
      "Cp = 1.005 #kJ/kg-K\n",
      "R = 8314/28.97 #J/kg-K\n",
      "\n",
      "#calculations:\n",
      "#assumptions\n",
      "Qdot = 0\n",
      "Wdot = Pc + Pf\n",
      "#inlet Area\n",
      "A1 = (1/V1)*(-1*Wdot/(Cp*(T2 - T1)*1000))*R*(T1+273)/101325\n",
      "\n",
      "D1 = ((4*A1/math.pi)**0.5)*100\n",
      "\n",
      "\n",
      "\n",
      "#Results\n",
      "print  \"smallest fan inlet diameter is\", round(D1,0),\"cm\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "smallest fan inlet diameter is 8.0 cm\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.09, page: 115</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "# Initialization  of  Variable\n",
      "P1 = 300 #lbf/in2\n",
      "P2 = 14.7 #lbf/in2\n",
      "T2 = 250 #degF\n",
      "\n",
      "#calculations:\n",
      "#From  TAble T-3E\n",
      "hf1 = 394.1 #Btu/lb\n",
      "hg1 = 1203.9 #Btu/lb\n",
      "#from Table T-4E\n",
      "h2 = 1168.8 #Btu/lb\n",
      "#Quality factor\n",
      "x1 = (h2-hf1)/(hg1-hf1)\n",
      "\n",
      "#Results\n",
      "print  \"quality in the line is\", round(x1*100,1),\"%\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "quality in the line is 95.7 %\n"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h2>Example 5.10, page: 117</h2>"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "# Initialization  of  Variable\n",
      "AV1 = 2E5 #ft3/min\n",
      "T1 = 400 #degF\n",
      "P1 = 1 #atm\n",
      "T2 = 260 #degF\n",
      "P2 = 1 #atm\n",
      "P3 = 40 #lbf/in2\n",
      "T3 = 102 #degF\n",
      "m3dot = 275 #lb/min\n",
      "P5 = 1 #lbf/in2\n",
      "x5 = 0.93\n",
      "cost = 0.08 #$/W\n",
      "t = 8000 #hrs\n",
      "R = 1545/28.97 #ft.lbf/lb-degR\n",
      "\n",
      "#calculations:\n",
      "#in Rankine\n",
      "T1r = (((T1-32)*5/9)+273)*1.8\n",
      "#mass flow rate at 1\n",
      "m1dot = AV1*(P1*14.7)*144/(R*T1r)\n",
      "#from Table T-9E\n",
      "h1 = 206.46 #Btu/lb\n",
      "h2 = 172.39 #Btu/lb\n",
      "#From Table T-2E\n",
      "h3 = 70 #Btu/lb\n",
      "#From Table T-3E\n",
      "hf5 = 69.74 #Btu/lb\n",
      "hg5 = 1105.74 #Btu/lb\n",
      "#enthapy at 5\n",
      "h5 = hf5 + x5*(hg5-hf5)\n",
      "#Power\n",
      "Wdot = m1dot*(h1-h2) + m3dot*(h3 - h5)\n",
      "\n",
      "#for temp T4\n",
      "h4 = h3 + (m1dot/m3dot)*(h1-h2)\n",
      "#from Table T-4E\n",
      "P4 = 40 #lbf/in2\n",
      "T4 = 354 #degF\n",
      "\n",
      "#total cost\n",
      "AnnVal = (Wdot*60/3413)*t*cost\n",
      "\n",
      "\n",
      "#From  TAble T-3E\n",
      "hf1 = 394.1 #Btu/lb\n",
      "hg1 = 1203.9 #Btu/lb\n",
      "#from Table T-4E\n",
      "h2 = 1168.8 #Btu/lb\n",
      "#Quality factor\n",
      "x1 = (h2-hf1)/(hg1-hf1)\n",
      "\n",
      "#Results\n",
      "print \"a) power developed by the turbine is\", round(Wdot,0),\"Btu/min\"\n",
      "print \"b) Turbine Inlet temperature is\", T4,\"degF\"\n",
      "print \"c) Annual Value of power developed is\", round(AnnVal,0),\"$/year\"\n",
      "#Answer wrong in book\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "a) power developed by the turbine is 49822.0 Btu/min\n",
        "b) Turbine Inlet temperature is 354 degF\n",
        "c) Annual Value of power developed is 560548.0 $/year\n"
       ]
      }
     ],
     "prompt_number": 9
    }
   ],
   "metadata": {}
  }
 ]
}