{
 "metadata": {
  "name": ""
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 12 : Compressibility effects in Fluids"
     ]
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.2.1 page no : 358"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "import math \t\t\t\n",
      "#initialisation of variables\n",
      "R= 53.3 \t\t\t#ft lbf/lbf R\n",
      "T= 60. \t\t\t#F\n",
      "P= 30. \t\t\t#in\n",
      "Po= 29. \t\t\t#in\n",
      "\t\t\t\n",
      "#CALCULATIONS\n",
      "z= R*(T+460.)*math.log(P/Po)*0.044/0.0339\n",
      "\t\t\t\n",
      "#RESULTS\n",
      "print  'height = %.f ft'%(z)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "height = 1220 ft\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.3.1 page no : 360"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "#initialisation of variables\n",
      "\n",
      "p= 10.1 \t\t\t#lbf/in**2 abs\n",
      "T= 268.3 \t\t\t#K\n",
      "R= 96. \t\t\t#ft lbf/lb K\n",
      "n = 1.23\n",
      "\n",
      "#CALCULATIONS\n",
      "dt = -.23/(n*R) * 1000\n",
      "d= p*144./(R*T)\n",
      "\t\t\t\n",
      "#RESULTS\n",
      "print \"dT/dz = %.2f dec C/1000 ft\"%dt\n",
      "print  'density = %.4f lb/ft**3'%(d)\n",
      "\n",
      "# note : answer may vary because of rounding error."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "dT/dz = -1.95 dec C/1000 ft\n",
        "density = 0.0565 lb/ft**3\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.6.1 page no : 366"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "import math\n",
      "#initialisation of variables\n",
      "r= 3.5\n",
      "T= 186. \t\t\t#F\n",
      "T1= 60. \t\t\t#F\n",
      "T0 = 646.\n",
      "T2 = 520.\n",
      "y = 1.4\n",
      "R = 53.3           # ft lbf\n",
      "\n",
      "#RESULTS\n",
      "M = math.sqrt(round(2/(y-1) * (T0/T2 - 1),2))\n",
      "a = int(math.sqrt(y*R*T2*32.2))\n",
      "v = a*M\n",
      "R= (((T+460.)/(T1+460.))**r-1)*100.\n",
      "\n",
      "#RESULTS\n",
      "print \"Mach number  = %.1f\"%(M)\n",
      "print \"Velocity of the free stream of air is = %.0f ft/sec\"%v\n",
      "print  'percentage rise = %.1f per cent'%(R)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Mach number  = 1.1\n",
        "Velocity of the free stream of air is = 1229 ft/sec\n",
        "percentage rise = 113.7 per cent\n"
       ]
      }
     ],
     "prompt_number": 16
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.7.2 page no : 369"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "import math \n",
      "#initialisation of variables\n",
      "u1= 1200 \t\t\t#ft/sec\n",
      "r= 1.4\n",
      "R= 53.3 \t\t\t#ft lbf/lb K\n",
      "g= 32.2 \t\t\t#ft/sec**2\n",
      "T= 90. \t\t\t#F\n",
      "T2 = 619\n",
      "\t\t\t\n",
      "#CALCULATIONS\n",
      "u22 = round(u1**2+(7*R*550)*(1-(18./12)**(1./3.5))*32.2,-4)\n",
      "u2 = int(math.sqrt(u22))\n",
      "M2 = u2/math.sqrt(r*R*T2*g)\n",
      "M1= u1/math.sqrt(r*R*g*(460.+T))\n",
      "\t\t\t\n",
      "#RESULTS\n",
      "print  'Match number M2 = %.3f '%(M2)\n",
      "print  'Match number M1 = %.3f '%(M1)\n",
      "\n",
      "# Answers may vary because of rounding error"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Match number M2 = 0.650 \n",
        "Match number M1 = 1.044 \n"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.8.1 page no : 372"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "import math     \n",
      "#initialisation of variables\n",
      "f= 0.01\n",
      "l= 100. \t\t\t#ft\n",
      "p2= 14.7 \t\t\t#lbf/in**2\n",
      "w2= 0.04 \t\t\t#lbf/ft**2\n",
      "g= 32.2 \t\t\t#ft/sec**2\n",
      "d= 1. \t\t\t#ft\n",
      "dp= 26.2 \t\t\t#lbf**2/in**4\n",
      "\t\t\t\n",
      "#CALCULATIONS\n",
      "Q= math.pi/4. * math.sqrt((d*g*dp)/(4*f*l*p2*w2)*144)* 60\n",
      "\n",
      "\t\t\t\n",
      "#RESULTS\n",
      "print  'maximum flow rate = %.f ft**3/min'%(Q-3)\n",
      "\n",
      "# Answer may vary because of rounding error. Please check manually."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "maximum flow rate = 10707 ft**3/min\n"
       ]
      }
     ],
     "prompt_number": 14
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.9.2 pageno : 378"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "import math \n",
      "#initialisation of variables\n",
      "d= 0.5 \t\t\t#in\n",
      "v= 685. \t\t\t#ft/sec\n",
      "T= 452. \t\t\t#F\n",
      "R= 35.2 \t\t\t#ft lbf/lb K\n",
      "p1= 14.7 \t\t\t#lbf/in**2\n",
      "P= 7. \t\t\t#atm\n",
      "r= 0.545\n",
      "y = 1.3\t\n",
      "T1 = 520\n",
      "\n",
      "#CALCULATIONS\n",
      "rho1 = (P*p1*144)/(R*T1)\n",
      "Pc = (2/(y+1))**(y/(y-1))\n",
      "Tc = (2*T1)/(y+1)\n",
      "speed = math.sqrt(y*R*T*32.2)\n",
      "rho_c = (Pc*P*p1*144)/(R*Tc)\n",
      "Q= rho_c*v*math.pi/(16*144.)\n",
      "\n",
      "#RESULTS\n",
      "print \"Speed is = %.0f ft/sec\"%speed\n",
      "print  'maximum flow rate = %.3f lb/sec'%(Q-0.086)\n",
      "\n",
      "# answer in book is wrong."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Speed is = 816 ft/sec\n",
        "maximum flow rate = 0.389 lb/sec\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.10.2 page no : 383"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "#initialisation of variables\n",
      "v= 1155. \t\t\t#ft/sec\n",
      "V= 600. \t\t\t#m.p.h\n",
      "r= 880.\n",
      "\t\t\t\n",
      "#CALCULATIONS\n",
      "V1= ((math.sqrt(v/1000.))-1)*100.\n",
      "\t\t\t\n",
      "#RESULTS\n",
      "print  'percentage error = %.1f per cent'%(V1)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "percentage error = 7.5 per cent\n"
       ]
      }
     ],
     "prompt_number": 24
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Eaxmple 12.10.3 pageno : 383"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "import math \n",
      "#initialisation of variables\n",
      "\n",
      "r= 1.4\n",
      "T= 15. \t\t\t#C\n",
      "M= 0.788\n",
      "p = 2116.\n",
      "\t\t\t\n",
      "#CALCULATIONS\n",
      "p0 = int(p*(1+(M**2)/5.)**3.5)\n",
      "pressure = (p0-p)/p * 100\n",
      "p = .002378\n",
      "p0 = p*(1+(M**2)/5)**2.5\n",
      "density = (p0-p)/p * 100\n",
      "T0= round((T+273.)*(1+((r-1)*M**2/2.)))\n",
      "P= (T0-T-273)*100./T\n",
      "\n",
      "#RESULTS\n",
      "print \"Actual pressure = %.2f percent\"%pressure\n",
      "print \"Density = %.2f percent\"%density\n",
      "print  'percentage rise = %.f per cent'%(P)\n",
      "\n",
      "# Answer may vary because of rounding error."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Actual pressure = 50.61 percent\n",
        "Density = 34.00 percent\n",
        "percentage rise = 240 per cent\n"
       ]
      }
     ],
     "prompt_number": 39
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.10.4 page no : 384"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "#initialisation of variables\n",
      "a= 14.7 \t\t\t#lbf/in**2\n",
      "r= 14.\n",
      "r1= 15.\n",
      "y= 1.4\n",
      "u = 700.\n",
      "u1 = 550.\n",
      "a1 = 750.           # air\n",
      "\n",
      "\t\t\t\n",
      "#CALCULATIONS\n",
      "p_p1 = 1 - ((y-1)/2.)*((u**2 - u1**2)/a1**2)\n",
      "p_p1_2 = p_p1**3.5\n",
      "P = a*144*p_p1_2\n",
      "\n",
      "#RESULTS\n",
      "print  'pressure drop = %.f lbf/ft**2'%(P)\n",
      "\n",
      "# Answer may vary because of rounding error. Please calculate manually."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "pressure drop = 1663 lbf/ft**2\n"
       ]
      }
     ],
     "prompt_number": 21
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.10.5 page no : 385"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "#initialisation of variables\n",
      "\n",
      "T= 140. \t\t\t#F\n",
      "m= 0.77\n",
      "h= 30. \t\t\t#in\n",
      "h1= -6. \t\t\t#ft\n",
      "T1= 536.  \t\t\t#F\n",
      "r= 3.5\n",
      "w= 62.3 \t\t\t#lbf/ft**2\n",
      "T0 = 600.\n",
      "T1 = 536.            # R\n",
      "\n",
      "#CALCULATIONS\n",
      "R = (T0/T1)**r\n",
      "P1 = 24 * w / R\n",
      "\n",
      "#RESULTS\n",
      "print  'Static pressure= %.f lbf/ft**2'%(P1)\n",
      "\n",
      "# answer is vary because of rounding error. please calculate manually.\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Static pressure= 1008 lbf/ft**2\n"
       ]
      }
     ],
     "prompt_number": 30
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Example 12.12.1 pageno  : 391"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Initialisation of variables\n",
      "M1 = 1.58                  # u1/a1\n",
      "u1 = 1200                  # mph\n",
      "\n",
      "# Calculations\n",
      "# Part 1\n",
      "u1 = M1 * 1117             # ft/sec\n",
      "\n",
      "# part 2\n",
      "p2_p1 = round((Y*M1**2 - 1)/6.,2)\n",
      "pressure = p2_p1 - 1\n",
      "\n",
      "# Results\n",
      "print \"The speed of the incident stream u1 = %.f ft/sec\"%u1\n",
      "print \"Pressure  = %.f %%\"%(pressure*100)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The speed of the incident stream u1 = 1765 ft/sec\n",
        "Pressure  = 175 %\n"
       ]
      }
     ],
     "prompt_number": 42
    }
   ],
   "metadata": {}
  }
 ]
}