{
 "metadata": {
  "name": "",
  "signature": "sha256:b83edffd06279953e9d9c22648c2ce548f84e9427ec05c98afa676ab5960d29e"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 3 : Macroscopic properties of pure substances"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 3.1 pg : 58"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\t\t\t\n",
      "# Variables\n",
      "V = 1.  \t\t\t#ft**3\n",
      "m = 30. \t\t\t#lbm\n",
      "\t\t\t\n",
      "# Calculations\n",
      "v = V/m\n",
      "vf1 = 0.01665\n",
      "vfg1 = 32.38 \t\t\t#ft**3/lbm\n",
      "x1 = 0.000515\n",
      "uf1 = 169.92\n",
      "ufg1 = 904.8\n",
      "u1 = uf1+x1*ufg1\n",
      "vfg = 0.0216\n",
      "vfg2 = 0.4240\n",
      "v2 = v\n",
      "x2 = 0.0277\n",
      "uf2 = 538.4\n",
      "ufg2 = 571.\n",
      "u2 = uf2+x2*ufg2\n",
      "Q = m*(u2-u1)\n",
      "\t\t\t\n",
      "# Results\n",
      "print \"Heat transfer  =  %d Btu\"%(round(Q,-2))\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Heat transfer  =  11500 Btu\n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 3.2 pg : 59"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\t\t\t\n",
      "# Variables\n",
      "V2 = 2.5 \t\t\t#ft**3\n",
      "V1 = 0.5 \t\t\t#ft**3\n",
      "P = 100. \t\t\t#psia\n",
      "x1 = 0.5\n",
      "\t\t\t\n",
      "# Calculations\n",
      "W = -P*144*(V2-V1)\n",
      "vf1 = 0.01774\n",
      "vfg1 = 4.414\n",
      "v1 = vf1+x1*vfg1\n",
      "m = V1/v1\n",
      "v2 = V2/m\n",
      "uf1 = 298.08\n",
      "ufg1 = 807.1\n",
      "u1 = uf1+x1*ufg1\n",
      "h2 = 1747.9\n",
      "u2 = h2-P*144*v2/778\n",
      "Q = m*(u2-u1)\n",
      "\t\t\t\n",
      "# Results\n",
      "print \"Amount of heat  =  %d Btu\"%(Q)\n",
      "\t\t\t#The answer for u2 is given wrong in the textbook. Please use a calculator to find it\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Amount of heat  =  188 Btu\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 3.3 pg : 60"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\t\t\t\n",
      "# Variables\n",
      "V1 = 1.735*10**-4 \t\t\t#ft**3\n",
      "v1 = 0.016080    \t\t\t#ft**3/lbm\n",
      "h1 = 70.61 \t\t    \t    #B/lbm\n",
      "P1 = 100. \t\t        \t#psia\n",
      "V2 = 1. \t\t    \t    #ft**3\n",
      "\t\t\t\n",
      "# Calculations\n",
      "u1 = h1-P1*v1*144/778.\n",
      "m = V1/v1\n",
      "v2 = V2/m\n",
      "vf2 = 0.01613\n",
      "vfg2 = 350.3\n",
      "x2 = (v2-vf2)/vfg2\n",
      "hf2 = 67.97\n",
      "hfg2 = 1037.2\n",
      "h2 = hf2+x2*hfg2\n",
      "P2 = 0.9492\n",
      "u2 = h2- P2*144*v2/778.\n",
      "Q = m*(u2-u1)\n",
      "\t\t\t\n",
      "# Results\n",
      "print \"Enthalpy change  =  %.2f Btu\"%(Q)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Enthalpy change  =  2.76 Btu\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 3.4 pg : 64"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\t\t\t\n",
      "# Variables\n",
      "P = 20. \t\t\t#psia\n",
      "V = 1.   \t\t\t#ft**3\n",
      "T = 560. \t\t\t#R\n",
      "cv = 0.1715\n",
      "Q = 10. \t\t\t#Btu\n",
      "\t\t\t\n",
      "# Calculations\n",
      "m = P*144*V/(53.35*T)\n",
      "T2 = Q/(m*cv) +T\n",
      "P2 = m*53.35*T2/V\n",
      "\t\t\t\n",
      "# Results\n",
      "print \"Fina pressure  =  %d lbf/ft**2\"%(P2)\n",
      "\n",
      "# note : rounding off error."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Fina pressure  =  5990 lbf/ft**2\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 3.5 pg : 66"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from scipy.integrate import quad\n",
      "\t\t\t\n",
      "# Variables\n",
      "T1 = 560 \t\t\t#R\n",
      "T2 = 3460  \t\t\t#R\n",
      "m = 28.02 \t\t\t#lb\n",
      "cv = 0.248\n",
      "\t\t\t\n",
      "# Calculations\n",
      "def fun(T):\n",
      "    return 9.47 - 3.29*10**3 /T +1.07*10**6 /T**2\n",
      "\n",
      "Q1 = quad(fun,T1,T2)[0]\n",
      "Q2 = m*cv*(T2-T1)\n",
      "Error = (Q1-Q2)/Q1\n",
      "\t\t\t\n",
      "# Results\n",
      "print \"Percentage error  =  %.1f percent\"%(Error*100)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Percentage error  =  12.7 percent\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 3.6 pg : 66"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\t\t\t\n",
      "# Variables\n",
      "import math \n",
      "rate = 20. \t\t\t#gal/min\n",
      "P1 = 20. \t\t\t#psia\n",
      "P2 = 1000. \t\t\t#psia\n",
      "T = 100.+460 \t\t\t#R\n",
      "\t\t\t\n",
      "# Calculations\n",
      "vf = 0.01613\n",
      "dv = -5.2*10**-5 \t\t\t#ft**3/lbm\n",
      "K = -dv/(vf*P2*144)\n",
      "wt = K*vf*(P2**2 - P1**2)*144*144*10**4 /2\n",
      "m = rate*8.33\n",
      "Wt = wt*m\n",
      "Wthp = Wt/33000\n",
      "\t\t\t\n",
      "# Results\n",
      "print \"Pump power required  =  %d hp\"%(Wthp)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Pump power required  =  188 hp\n"
       ]
      }
     ],
     "prompt_number": 8
    }
   ],
   "metadata": {}
  }
 ]
}