diff options
Diffstat (limited to 'Hydraulics_by_J._Lal/Chapter_2.ipynb')
-rwxr-xr-x | Hydraulics_by_J._Lal/Chapter_2.ipynb | 399 |
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 |