summaryrefslogtreecommitdiff
path: root/Hydraulics_by_J._Lal/Chapter_2.ipynb
diff options
context:
space:
mode:
Diffstat (limited to 'Hydraulics_by_J._Lal/Chapter_2.ipynb')
-rwxr-xr-xHydraulics_by_J._Lal/Chapter_2.ipynb399
1 files changed, 399 insertions, 0 deletions
diff --git a/Hydraulics_by_J._Lal/Chapter_2.ipynb b/Hydraulics_by_J._Lal/Chapter_2.ipynb
new file mode 100755
index 00000000..8353a544
--- /dev/null
+++ b/Hydraulics_by_J._Lal/Chapter_2.ipynb
@@ -0,0 +1,399 @@
+{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:e77c0ecf71b6b19d4b74c1b5653eedac5bf7ec259a0d03cbae3bef3229b3cdac"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter 2 : Hydrodynamics"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 2.1 pageno : 41"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "#initialisation of variables\n",
+ "\n",
+ "import math \n",
+ "\n",
+ "Q = 0.8 \t\t\t#ft**3/sec\n",
+ "w = 62.4 \t\t\t#lb/sec\n",
+ "d1 = 3. \t\t\t#in\n",
+ "d2 = 1.5 \t\t\t#in\n",
+ "#CALCULATIONS\n",
+ "Q1 = Q*w*60/10.\n",
+ "a1 = math.pi*(d1/12.)**2/4.\n",
+ "a2 = math.pi*(d2/12.)**2/4.\n",
+ "v1 = Q/a1\n",
+ "v2 = Q/a2\n",
+ "\n",
+ "#RESULTS\n",
+ "print 'Q = %.f gpm'%Q1\n",
+ "print 'v1 = %.1f ft/sec '%(v1)\n",
+ "print 'v2 = %.1f ft/sec '%(v2)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Q = 300 gpm\n",
+ "v1 = 16.3 ft/sec \n",
+ "v2 = 65.2 ft/sec \n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 2.2 page no : 45"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "#initialisation of variables\n",
+ "import math \n",
+ "d1 = 12. \t\t\t#in\n",
+ "d2 = 9. \t\t\t#in\n",
+ "z1 = 10. \t\t\t#ft\n",
+ "z2 = 10. \t\t\t#ft\n",
+ "p1 = 15. \t\t\t#lb/in**2\n",
+ "w = 62.4 \t\t\t#lb/ft**3\n",
+ "Q = 2. \t\t\t#cuses\n",
+ "g = 32.2 \t\t\t#ft/sec**2\n",
+ "#CALCULATIONS\n",
+ "v1 = Q/(math.pi*(d1/12.)**2/4.)\n",
+ "v2 = Q/(math.pi*(d2/12.)**2/4.)\n",
+ "p2 = w*(z1-z2+(p1*144/w)+(v1**2/(2*g))-(v2**2/(2*g)))/144\n",
+ "#RESULTS\n",
+ "print 'p2 = %.2f lb/in**2 '%(p2)\n",
+ "\n",
+ "# note : answer is slightly differet because of rounding error."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "p2 = 14.91 lb/in**2 \n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 2.3 page no : 46"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "\n",
+ "#initialisation of variables\n",
+ "import math \n",
+ "d0 = 4. \t\t\t#ft\n",
+ "d2 = 2. \t\t\t#ft\n",
+ "z0 = 0 \t\t\t#ft\n",
+ "z1 = 5. \t\t\t#ft\n",
+ "z2 = 13. \t\t\t#ft\n",
+ "h = 9.5 \t\t\t#in\n",
+ "w = 62.4 \t\t\t#lb/ft**3\n",
+ "w1 = 30. \t\t\t#lb/ft**3\n",
+ "g = 32.2 \t\t\t#ft/sec**2\n",
+ "r = 0.1\n",
+ "#CALCULATIONS\n",
+ "p2 = -h*34/w1\n",
+ "v2 = math.sqrt(2*g*(z1-p2-z2)/(1+r))\n",
+ "Q = math.pi*(d2/12)**2*v2*w*60/(10*4)\n",
+ "#RESULTS\n",
+ "print 'Discharge = %.f gpm '%(Q)\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Discharge = 104 gpm \n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 2.4 pageno :47"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "\n",
+ "#initialisation of variables\n",
+ "\n",
+ "import math \n",
+ "d1 = 2 \t\t\t#ft\n",
+ "d2 = 3 \t\t\t#ft\n",
+ "v1 = 20 \t\t\t#ft/sec\n",
+ "z1 = 20 \t\t\t#ft\n",
+ "z2 = 0 \t\t\t#ft\n",
+ "h = 5 \t\t\t#ft\n",
+ "w = 62.4 \t\t\t#lb/ft**3\n",
+ "g = 32.2 \t\t\t#ft/sec**2\n",
+ "#CALCULATIONS\n",
+ "Hl = v1**2*0.15/(2*g)\n",
+ "a1 = math.pi*d1**2/4\n",
+ "a2 = math.pi*d2**2/4\n",
+ "v2 = a1*v1/a2\n",
+ "p1 = ((h-z1+(v2**2)/(2*g))-(0.85*v1**2/(2*g)))\n",
+ "#RESULTS\n",
+ "print 'water pressure at top = %.2f ft of water '%(p1)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "water pressure at top = -19.05 ft of water \n"
+ ]
+ }
+ ],
+ "prompt_number": 7
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 2.5 page no : 50"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "\n",
+ "#initialisation of variables\n",
+ "\n",
+ "import math \n",
+ "d1 = 15. \t\t\t#in\n",
+ "d2 = 6.\t\t\t#in\n",
+ "h = 10. \t\t\t#in of mercury\n",
+ "C = 0.98\n",
+ "sm = 13.6\n",
+ "w = 12.\n",
+ "g = 32.2 \t\t\t#ft/sec**2\n",
+ "#CALCULATIONS\n",
+ "a1 = math.pi*(d1/12)**2/4\n",
+ "a2 = math.pi*(d2/12)**2/4\n",
+ "h1 = h*(sm-1)/w\n",
+ "Q = round(C*(a1*a2/(math.sqrt(a1**2-a2**2)))*math.sqrt(2*g)*math.sqrt(h1)*6.24*60*60,-2)\n",
+ "\n",
+ "#RESULTS\n",
+ "print 'Discharge = %.f gph '%(Q)\n",
+ "\n",
+ "# answer is different because of rounding error. Please calculate manually."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Discharge = 113900 gph \n"
+ ]
+ }
+ ],
+ "prompt_number": 12
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 2.6 page no : 52"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "\n",
+ "#initialisation of variables\n",
+ "import math \n",
+ "d1 = 8. \t\t\t#in\n",
+ "d2 = 4. \t\t\t#in\n",
+ "h = 10. \t\t\t#in of mercury\n",
+ "Cd = 0.98\n",
+ "g = 32.2 \t\t\t#ft/sec**2\n",
+ "sm = 13.56\n",
+ "#CALCULATIONS\n",
+ "a1 = math.pi*(d1/12)**2/4\n",
+ "a2 = math.pi*(d2/12)**2/4\n",
+ "h1 = h*(sm-1)/12\n",
+ "Q = a1*a2*math.sqrt(2*g)*math.sqrt(h1)/math.sqrt(a1**2-a2**2)\n",
+ "Qactual = Cd*Q\n",
+ "#RESULTS\n",
+ "print 'Actual discharge = %.2f cusecs '%(Qactual)\n",
+ "\n",
+ "# Note : Answer is slightly different because of rounding error."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Actual discharge = 2.29 cusecs \n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 2.7 page no : 54"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "\n",
+ "#initialisation of variables\n",
+ "\n",
+ "import math \n",
+ "h = 6.8 \t\t\t#in of mercury\n",
+ "sm = 13.6\n",
+ "ssw = 1.026\n",
+ "g = 32.2 \t\t\t#ft/sec**2\n",
+ "#CALCULATIONS\n",
+ "V = math.sqrt(2*g*h*(sm-ssw)/12)*3600/5280\n",
+ "#RESULTS\n",
+ "print 'speed of submarine = %.1f miles per hour '%(V)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "speed of submarine = 14.6 miles per hour \n"
+ ]
+ }
+ ],
+ "prompt_number": 17
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 2.8 page no : 58"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "\n",
+ "#initialisation of variables\n",
+ "\n",
+ "import math \n",
+ "g = 32.2 \t\t\t#ft/sec**2\n",
+ "d1 = 2. \t\t\t#in\n",
+ "d2 = 12. \t\t\t#in\n",
+ "r = 1.4\n",
+ "n = 0.905\n",
+ "Q = 2995. \t\t\t#lb/ft**2\n",
+ "w = 0.083 \t\t\t#lb/ft**3\n",
+ "#CALCULATIONS\n",
+ "V1 = round(1/w,2)\n",
+ "n1 = round(n**((r-1)/r),2)\n",
+ "n2 = n**(2/r)\n",
+ "Q = math.pi*(d1/12)**2*math.sqrt(2*g*Q*(1-n1)*r/((r-1)*n2*(1-(d1/d2)**2)))\n",
+ "\n",
+ "#RESULTS\n",
+ "print 'Volume of air passing through the Venturimeter = %.1f cuses '%(Q)\n",
+ "\n",
+ "# note : answer is different because of rounding error."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Volume of air passing through the Venturimeter = 13.5 cuses \n"
+ ]
+ }
+ ],
+ "prompt_number": 7
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [],
+ "language": "python",
+ "metadata": {},
+ "outputs": []
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+} \ No newline at end of file