{
 "metadata": {
  "name": "",
  "signature": "sha256:0acf769d31306f5bc291e2574308a74de8651b4a4723b52ee59419306e6b7112"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 1 : Introduction"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.1 Page No : 20"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "d_r = 13640.; \t\t\t# Density of mercury in kg/m3\n",
      "g = 9.79; \t\t\t    # Acceleration due to gravity in m/s2\n",
      "z = 562e-03; \t\t\t# Difference in height in m\n",
      "z0 = 761e-03; \t\t\t# Reading of barometer in m\n",
      "\n",
      "# Calculation\n",
      "P = (d_r*g*(z+z0))*(0.987/1e05); \t\t\t# Gas Pressure in bar\n",
      "\n",
      "# Results\n",
      "print \"Gas Pressure is %.3f bar\"%P\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Gas Pressure is 1.744 bar\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.2 Page No : 21"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "d_r = 13.6e03; \t\t\t# Density of mercury in kg/m3\n",
      "g = 9.81; \t\t\t# Acceleration due to gravity in m/s2\n",
      "z = 710e-03; \t\t\t# Stean flow pressure in m\n",
      "z0 = 772e-03; \t\t\t# Reading of barometer in m\n",
      "P = 1.4e06; \t\t\t# Gauge pressure of applied steam in Pa\n",
      "\n",
      "# Calculation\n",
      "P0 = d_r*g*z0; \t\t\t# Atmospheric pressure in Pa\n",
      "Pi = P+P0 ; \t\t\t# Inlet steam pressure in Pa\n",
      "Pc = d_r*g*(z0-z); \t\t\t# Condenser pressure in Pa\n",
      "\n",
      "# Results\n",
      "print \"Inlet steam pressure is %.3e Pa\"%Pi\n",
      "print \"Condenser pressure is %.2e Pa\"%Pc\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Inlet steam pressure is 1.503e+06 Pa\n",
        "Condenser pressure is 8.27e+03 Pa\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.3 Page No : 21"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "z = 0.760; \t\t\t# Barometer reading in m\n",
      "# Part (a)\n",
      "h1 = 40e-02; \t\t\t# Mercury height in vaccume in m\n",
      "d_r = 13.6e03; \t\t\t# Density of mercury in kg/m3\n",
      "g = 9.80; \t\t\t# Acceleration due to gravity in m/s2\n",
      "\n",
      "# Calculation and Results\n",
      "Patm = z*d_r*g; \t\t\t# Atmospheric pressure in Pas\n",
      "Pv = h1*d_r*g; \t\t\t# Pressue in vaccume in Pa\n",
      "Pabst = Patm-Pv; \t\t\t# Absolute pressure in Pa\n",
      "Pabs = 101.325 - Pv/1000\n",
      "\n",
      "print \"Pvaccum = %.2f kPa\"%(Pv/1000)\n",
      "print \"Pabsolute = %.f\"%(Pabst/1000),\"kPa\"\n",
      "print \"40 cmHg vacuum : %.3f kPa\"%Pabs\n",
      "\n",
      "# Part (b)\n",
      "h2 = 90e-02; \t\t\t# Mercury height in gauge in m\n",
      "Pg = h2*d_r*g; \t\t\t# Gauge Pressure in Pa\n",
      "Pabs1 = Patm + Pg ; \t\t\t# Absolute pressure in Pa\n",
      "\n",
      "print \"\\nPgauge = %.f kPa\"%(Pg/1000)\n",
      "print \"90cmHg gauge is %.3f\"%(Pabs1/1000),\"kPa\"\n",
      "\n",
      "# Part(c)\n",
      "d_w = 1e03 ; \t\t\t# Density of water in kg/m3\n",
      "h3 = 1.2 ; \t\t\t# Gauge Pressure water height in m\n",
      "Pga = d_w*h3*g; \t\t\t# Gauge Pressure in Pa\n",
      "Pabs3 = Patm + Pga ; \t\t\t# Absolute pressure in Pa\n",
      "print \"\\n1.2 m H2O gauge is %.3f\"%(Pabs3/1000),\"kPa\"\n",
      "\n",
      "# rounding off error"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Pvaccum = 53.31 kPa\n",
        "Pabsolute = 48 kPa\n",
        "40 cmHg vacuum : 48.013 kPa\n",
        "\n",
        "Pgauge = 120 kPa\n",
        "90cmHg gauge is 221.245 kPa\n",
        "\n",
        "1.2 m H2O gauge is 113.053 kPa\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.4 Page No : 22"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from scipy.integrate import quad\n",
      "\n",
      "# Variables\n",
      "Pr = 1.033e05; \t\t\t# Required Pressure in bar\n",
      "\n",
      "# Calculation\n",
      "def pressure(p):\n",
      "    return p**(-0.714);\n",
      "g = 9.81; \t\t\t# Acceleration due to gravity in m/s2\n",
      "H = ((2.5e05**0.714)/g)* quad(pressure,0,Pr)[0]; \t\t\t# Depth of atmosphere required in m\n",
      "\n",
      "# Results\n",
      "print \"The depth of atmosphere required is %.3f Km\"%(H/1000)\n",
      "  \n",
      "# note : there will be rounding off error because of quad function"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The depth of atmosphere required is 69.203 Km\n"
       ]
      }
     ],
     "prompt_number": 23
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.5 Page No : 22"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "m = 68. ; \t\t\t# Astronaut mass in Kg\n",
      "g = 9.806; \t\t\t# Acceleration due to gravity in m/s2\n",
      "\n",
      "# Calculation\n",
      "a = 10*g ; \t\t\t# Lift off acceleration in m/s2\n",
      "F = m*a; \t\t\t# Net vertical force in N\n",
      "\n",
      "# Results\n",
      "print \"Net vertical force experienced by astronaut is %.0f N\"%F\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Net vertical force experienced by astronaut is 6668 N\n"
       ]
      }
     ],
     "prompt_number": 6
    }
   ],
   "metadata": {}
  }
 ]
}