adding book

1 files changed, 457 insertions, 0 deletions

diff --git a/Mechanics_Of_Fluids/ch4.ipynb b/Mechanics_Of_Fluids/ch4.ipynb
new file mode 100755
index 00000000..8c70a7d4
--- /dev/null
+++ b/Mechanics_Of_Fluids/ch4.ipynb
@@ -0,0 +1,457 @@
+{
+ "metadata": {
+  "name": ""
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+  {
+   "cells": [
+    {
+     "cell_type": "heading",
+     "level": 1,
+     "metadata": {},
+     "source": [
+      "Chapter 4 : Bernoulli's Equation and Measurement of Flow of Incompressible Fluids"
+     ]
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.1.1 page no : 74"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialisation of variables\n",
+      "\n",
+      "H= 33. \t\t\t#ft lbf/lbf\n",
+      "Q= 100. \t\t\t#ft**3/min\n",
+      "w= 62.4 \t\t\t#lbf/ft**3\n",
+      "s= 0.8\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "P= s*w*Q*H/33000.\n",
+      "\t\t\t\n",
+      "#RESULTS\n",
+      "print  ' power required= %.2f h.p'%(P)\n",
+      "    "
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " power required= 4.99 h.p\n"
+       ]
+      }
+     ],
+     "prompt_number": 1
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.2.2 page no : 78"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialisation of variables\n",
+      "\n",
+      "import math \n",
+      "g= 32.2 \t\t\t#ft/sec**2\n",
+      "h= 1. \t\t\t    #in\n",
+      "ww= 62.4 \t\t\t#lbf/ft**3\n",
+      "w= 0.0764 \t\t\t#lbf/ft**3\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "u= math.sqrt(2*g*h*(1./12)*(ww/w))\n",
+      "\t\t\t\n",
+      "#RESULTS\n",
+      "print  ' speed of air through the tunnel= %.1f ft/sec'%(u)\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " speed of air through the tunnel= 66.2 ft/sec\n"
+       ]
+      }
+     ],
+     "prompt_number": 2
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.3.1. page no : 82"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialisation of variables\n",
+      "import math \n",
+      "za= 0. \t\t\t#ft\n",
+      "zb= 12. \t\t\t#ft\n",
+      "w= 62.3 \t\t\t#lbf/ft**2\n",
+      "pa= 750. \t\t\t#lbf/in**2\n",
+      "p= 700. \t\t\t#lbf/in**2\n",
+      "ua= 3. \t\t\t#ft/sec\n",
+      "g= 32.2 \t\t\t#ft/sec**2\n",
+      "d= 2. \t\t\t#in**2\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "ub= 4*ua\n",
+      "Hl= (za-zb)+((pa-p)*144/w)+(ua**2-ub**2)/(2*g)\n",
+      "P= (w*ua*(math.pi/4)*Hl*d**2)/(144*550.)\n",
+      "\t\t\t\n",
+      "#RESULTS\n",
+      "print  ' horse-power expended in over coming losses= %.2f h.p'%(P)\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " horse-power expended in over coming losses= 0.75 h.p\n"
+       ]
+      }
+     ],
+     "prompt_number": 3
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.4.1 pageno : 85"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialisation of variables\n",
+      "import math \n",
+      "d= 1. \t\t\t#in\n",
+      "d1= 3. \t\t\t#in\n",
+      "h= 9. \t\t\t#in\n",
+      "p= 3. \t\t\t#percent\n",
+      "g= 32.2 \t\t\t#ft/sec**2\n",
+      "s= 13.6 \t\t\t#gm/cm**3\n",
+      "a= 0.97\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "Ka= 1./(1-(d/d1)**2)\n",
+      "C= Ka*math.pi*(d/2)**2*math.sqrt(2*g*(s-1))/144\n",
+      "C1= a*C\n",
+      "Q= C1*h/12.\n",
+      "\t\t\t\n",
+      "#RESULTs\n",
+      "print  ' flow rate = %.2f ft**3/sec '%(Q)\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " flow rate = 0.13 ft**3/sec \n"
+       ]
+      }
+     ],
+     "prompt_number": 5
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.4.2 page no : 86"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialisation of variables\n",
+      "import math \n",
+      "Q= 1.4 \t\t\t#ft**3/sec\n",
+      "d= 6. \t\t\t#in\n",
+      "d1= 3. \t\t\t#in\n",
+      "h= 9. \t\t\t#in\n",
+      "s= 13.6/0.78\n",
+      "Cd= 0.96\n",
+      "g= 32.2 \t\t\t#ft/sec**2\n",
+      "w= 62.3 \t\t\t#lb/ft**3\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "h1= (Q*4.*12**2/(Cd*math.pi*d1**2))**2*(1-(d1/(2*d))**2)/(2*g*(s-1))\n",
+      "dpbyw= (h/12.)+((s)-1)*h1\n",
+      "dp= dpbyw*h1*w/144.\n",
+      "\t\t\t\n",
+      "#RESULTS\n",
+      "print  ' pressure difference = %.2f lbf/in**2 '%(dp)\n",
+      "\n",
+      "# Note : answer may be vary because of rounding error."
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " pressure difference = 4.60 lbf/in**2 \n"
+       ]
+      }
+     ],
+     "prompt_number": 1
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.5.1 page no : 88"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialisation of variables\n",
+      "import math \n",
+      "C= 0.6\n",
+      "s= 0.0767 \t\t\t#lbf/ft**3\n",
+      "g= 32.2 \t\t\t#ft/sec**2\n",
+      "w= 62.4 \t\t\t#lbf/ft**3\n",
+      "Hw= 0.7 \t\t\t#in\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "Ha= Hw*w/(s*12)\n",
+      "Q= C*math.pi*math.sqrt(2*g*Ha)/144.\n",
+      "\t\t\t\n",
+      "#RESULTS\n",
+      "print  ' volumetric flow rate = %.3f ft**3/sec '%(Q)\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " volumetric flow rate = 0.724 ft**3/sec \n"
+       ]
+      }
+     ],
+     "prompt_number": 7
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.5.3 pageno : 89"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "\n",
+      "#initialisation of variables\n",
+      "import math \n",
+      "g= 32.2 \t\t#ft/sec**2\n",
+      "h= 5. \t\t\t#ft\n",
+      "Q= 0.6 \t\t\t#ft**3/sec\n",
+      "Cd= 0.6\n",
+      "d= 2. \t\t\t#in\n",
+      "Q1= 0.315 \t\t#ft**3/sec\n",
+      "h1= 8. \t\t\t#ft\n",
+      "h2= 2. \t\t\t#ft\n",
+      "A= 9. \t\t\t#ft**2\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "H= Q1**2./((Cd*math.pi*(d/24.)**2.)**2.*2*g)\n",
+      "T= A*2*(math.sqrt(h1)-math.sqrt(h2))/(Cd*math.sqrt(2.*g)*60.*(d/24.)**2)\n",
+      "dhbyt= (Q-Cd*math.pi*(d/24.)**2*math.sqrt(2*g*h))*60*12/(math.pi*A)\n",
+      "\t\t\t\n",
+      "#RESULTS\n",
+      "print  'depth of the water = %.f ft '%(H)\n",
+      "print  'time taken = %.1f min '%(T)\n",
+      "print  'rate of rise in water = %.1f in/min '%(dhbyt)\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "depth of the water = 9 ft \n",
+        "time taken = 12.7 min \n",
+        "rate of rise in water = 9.3 in/min \n"
+       ]
+      }
+     ],
+     "prompt_number": 9
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.6.1 page no : 91"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "\n",
+      "#initialisation of variables\n",
+      "import math \n",
+      "d= 8. \t\t\t#in\n",
+      "d1= 1.5 \t\t#in\n",
+      "Cd= 0.65\n",
+      "w= 62.3 \t\t#lbf.ft**3\n",
+      "W= 25. \t\t\t#tonf\n",
+      "u= 5. \t\t\t#miles/hour\n",
+      "u1= 20.5        #miles/hour\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "ds= W*2240*d1**4*Cd**2*math.log(u1/u)/(w*d**4*math.pi*(d/24)**2)\n",
+      "T= W*2240*d1**4*Cd**2*((5/(u*7.33))-(20/(u1*29.35)))/(w*d**4*math.pi*(d/24)**2)\n",
+      "\n",
+      "\n",
+      "#RESULTS\n",
+      "print  ' Distance that piston moves= %.2f ft '%(ds)\n",
+      "print  ' time taken = %.4f sec '%(T)\n",
+      "print \" Answers may vary because of rounding error.\""
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " Distance that piston moves= 1.90 ft \n",
+        " time taken = 0.1388 sec \n",
+        " Answers may vary because of rounding error.\n"
+       ]
+      }
+     ],
+     "prompt_number": 6
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.7.1 page no : 94"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialisation of variables\n",
+      "\n",
+      "c= 0.002378 \t\t\t#slug/ft**3\n",
+      "u= 420. \t\t\t#mile/hour.\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "P= 0.5*c*u*616**2/420.\n",
+      "\t\t\t\n",
+      "#RESULTS\n",
+      "print  ' Dynamic pressure= %.f lbf/ft**2 '%(P)\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " Dynamic pressure= 451 lbf/ft**2 \n"
+       ]
+      }
+     ],
+     "prompt_number": 11
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "Example 4.8.2 page no : 97"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialisation of variables\n",
+      "import math \n",
+      "g= 32.2 \t\t\t#ft/sec**2\n",
+      "A= 13. \t\t\t#in**2\n",
+      "l= 10. \t\t\t#in**1.5\n",
+      "\t\t\t\n",
+      "#CALCULATIONS\n",
+      "Q= 2*math.pi*1.05*math.sqrt(2*g*12.)*A*l/1728.\n",
+      "\t\t\t\n",
+      "#RESULTS\n",
+      "print  ' Rate of flow= %.1f ft**3/sec '%(Q)\n",
+      "print 'Answer may vary because of rounding error'"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " Rate of flow= 13.8 ft**3/sec \n",
+        "Answer may vary because of rounding error\n"
+       ]
+      }
+     ],
+     "prompt_number": 13
+    }
+   ],
+   "metadata": {}
+  }
+ ]
+}
\ No newline at end of file