diff options
author | Jovina Dsouza | 2014-06-18 12:43:07 +0530 |
---|---|---|
committer | Jovina Dsouza | 2014-06-18 12:43:07 +0530 |
commit | 206d0358703aa05d5d7315900fe1d054c2817ddc (patch) | |
tree | f2403e29f3aded0caf7a2434ea50dd507f6545e2 /Fluid_Mechanics_/Chapter9.ipynb | |
parent | c6f0d6aeb95beaf41e4b679e78bb42c4ffe45a40 (diff) | |
download | Python-Textbook-Companions-206d0358703aa05d5d7315900fe1d054c2817ddc.tar.gz Python-Textbook-Companions-206d0358703aa05d5d7315900fe1d054c2817ddc.tar.bz2 Python-Textbook-Companions-206d0358703aa05d5d7315900fe1d054c2817ddc.zip |
adding book
Diffstat (limited to 'Fluid_Mechanics_/Chapter9.ipynb')
-rw-r--r-- | Fluid_Mechanics_/Chapter9.ipynb | 849 |
1 files changed, 849 insertions, 0 deletions
diff --git a/Fluid_Mechanics_/Chapter9.ipynb b/Fluid_Mechanics_/Chapter9.ipynb new file mode 100644 index 00000000..6c911620 --- /dev/null +++ b/Fluid_Mechanics_/Chapter9.ipynb @@ -0,0 +1,849 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:81c0a87082bc2a398c522de3be6bcb735589adaa60c31ff7051c7c79b83652d0" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Chapter 9 : Turbulent flow in Pipes" + ] + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.1 Page no 308" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Determine Head loss\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "S = 1.26 # specific gravity\n", + "\n", + "mu = 0.826 # kinematic viscosity in Ns/m**2\n", + "\n", + "# for water\n", + "\n", + "rho = 998 # density of water in kg/m**3\n", + "\n", + "mu1 = 1.005*10**-3 # viscosity in Ns/m**2\n", + "\n", + "# for glycerine\n", + "\n", + "rho1 = S*rho # density of glycerine in kg/m**3\n", + "\n", + "Q = 0.1 # discharge in m**3/s\n", + "\n", + "d1 = 0.2 # diameter in m\n", + "\n", + "A = pi*d1**2/4 # area in m**2\n", + "\n", + "g = 9.81 # acceleration due to gravity in m/s**2\n", + "\n", + "l =100 # length of the pipe\n", + "\n", + "# Solution\n", + "\n", + "V = Q/A\n", + "\n", + "R = rho1*V*d1/mu\n", + "\n", + "print \"It is a laminar flow\"\n", + "\n", + "f = 64/R # friction factor\n", + "\n", + "Hf = f*l*V**2/(2*g*d1) # head loss due to friction\n", + "\n", + "print \"(a) Head loss due to flow for glycerine =\",round(Hf,1),\"m \"\n", + "\n", + "R1 = rho*V*d1/mu1\n", + "\n", + "print \"The flow is turbulent\"\n", + "\n", + "e = 0.025\n", + "\n", + "r = e/(d1*100)\n", + "\n", + "f = 0.021\n", + "\n", + "hf = f*l*V**2/(2*g*d1)\n", + "\n", + "print \"(a) Head loss due to flow for water =\",round(hf,2),\"m \"\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "It is a laminar flow\n", + "(a) Head loss due to flow for glycerine = 17.1 m \n", + "The flow is turbulent\n", + "(a) Head loss due to flow for water = 5.42 m \n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.2 Page no 311" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Discharge of water\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "# for water\n", + "\n", + "nu = 1.007*10**-6 # viscosity in m**2/s\n", + "\n", + "e = 0.025 # for cast iron in cm\n", + "\n", + "L = 100 # length of the pipe in m\n", + "\n", + "D = 0.2 # diameter in m\n", + "\n", + "hf = 5.43 # head loss due to friction\n", + "\n", + "r = e/(D*100)\n", + "\n", + "g = 9.81 # acceleration due to gravity in m/s**2\n", + "\n", + "# Solution\n", + "\n", + "A = sqrt(2*g*D*hf/L)\n", + "\n", + "B = D/nu\n", + "\n", + "f = 0.021 # from moodys diagram\n", + "\n", + "V = A/sqrt(f)\n", + "\n", + "print V\n", + "\n", + "R = B*f\n", + "\n", + "A = pi*D**2/4\n", + "\n", + "Q = A*V\n", + "\n", + "print \"Discharge =\",round(Q,2),\"m**3/s\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "3.1853324563\n", + "Discharge = 0.1 m**3/s\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.3 Page no 314" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Size of the case iron\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "Q =0.1 # discharge in m**3/s\n", + "\n", + "hf = 5.43 # friction loss head in m\n", + "\n", + "L = 100 # length of pipe\n", + "\n", + "nu = 1.00*10**-6 # viscosity in m**2/s\n", + "\n", + "e = 0.025 # for cast iron in cm\n", + "\n", + "g = 9.81 # acceleration due to gravity in m/s**2\n", + "\n", + "# Solution\n", + "\n", + "A = 8*L*Q**2/(hf*g*pi**2)\n", + "\n", + "B = 4*Q/(pi*nu)\n", + "\n", + "# for D = 0.172 ; f=0.01\n", + "D = 0.172\n", + "\n", + "r = e/D\n", + "\n", + "Re = B/D\n", + "\n", + "f = 0.022 # for Re and r\n", + "\n", + "# for D1=0.199 ; f=0.021\n", + "\n", + "D1 = 0.199\n", + "\n", + "r1 = e/D1\n", + "\n", + "R = B/D1\n", + "\n", + "f = 0.021 # for R and r\n", + "\n", + "print \"Hence the convergence is attained, D=\",round(D1,1),\"m\"\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Hence the convergence is attained, D= 0.2 m\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.4 Page no 318" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Head loss \n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "L = 500 # length of the pipe in ft\n", + "\n", + "D= 9*2.54/100 # diameter in cm\n", + "\n", + "C = 100 # constant\n", + "\n", + "S = 0.004\n", + "\n", + "# Solution\n", + "\n", + "Hf = S*L\n", + "\n", + "print \"Head loss =\",round(Hf,0),\"ft\"\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Head loss = 2.0 ft\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.5 Page no 319" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Head loss of water\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "Q = 0.1 # water flow rate in m**3/s\n", + "\n", + "d = 30 # diameter in m\n", + "\n", + "l = 500 # length in m\n", + "\n", + "e = 0.025 # for cast iron\n", + "\n", + "g = 9.81 # acceleration due to gravity in m/s**2\n", + "\n", + "# Solution\n", + "\n", + "r = log(d/e,10)\n", + "\n", + "K = (pi/4)*sqrt(2*g)*(2*r+1.14)*(0.3)**(2.5)\n", + " \n", + "S = (Q/K)**2\n", + "\n", + "hf = S*l\n", + "\n", + "print \"Head loss of water =\",round(hf,1),\"m\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Head loss of water = 3.2 m\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.6 Page no 319" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Head loss by conveyance method\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "Q = 0.1 # water flow rate in m**3/s\n", + "\n", + "d = 20 # diameter in cm\n", + "\n", + "l = 500 # length in m\n", + "\n", + "e = 0.025 # for cast iron\n", + "\n", + "g = 9.81 # acceleration due to gravity in m/s**2\n", + "\n", + "S = 5.43 \n", + "\n", + "# Solution\n", + "\n", + "r = log(d/e,10)\n", + "\n", + "K = (pi/4)*sqrt(2*g)*(2*r+1.14)*(0.2)**2.5\n", + "\n", + "Q=K*sqrt(S/100)\n", + "\n", + "print \"Head loss of water =\",round(Q,2),\"m**3/s\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Head loss of water = 0.1 m**3/s\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.7 Page no 320" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Solve using the conveyence method\n", + "\n", + "from math import *\n", + "\n", + "# given\n", + "\n", + "eps = 0.025*10**-2 # for cast iron epsilon = 0.0025 cm\n", + "\n", + "# we get the value of K = 0.432 m**2/s\n", + "# we need to do trial and error to find the value of D\n", + "\n", + "# we use the value of D = 0.2 m\n", + "\n", + "D = 0.2 # value in m\n", + "\n", + "g = 9.81\n", + "\n", + "# Solution\n", + "\n", + "K = (pi/4)*sqrt(2*g)*(2*log10(D/(eps))+1.14)*D**(2.5)\n", + "\n", + "print \"K = \",round(K,3),\" from trial and error\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "K = 0.432 from trial and error\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.8 Page no 326" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Determine head loss\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "d = 0.1 # diameter of the pipe\n", + "\n", + "Q= 0.075 # discharge in m**3/s\n", + "\n", + "L = 30 # length in m\n", + "\n", + "A = pi*d**2/4\n", + "\n", + "g = 9.81 # acceleration due to gravity in m/s**2\n", + "\n", + "# for water\n", + "\n", + "nu = 1.007*10**-6 # viscosity in m**2/s\n", + "\n", + "e = 0.025\n", + "\n", + "r = e/(10*d)\n", + "\n", + "# Solution\n", + "\n", + "V = Q/A\n", + "\n", + "Re = V*d/nu\n", + "\n", + "f = 0.025 # firction factor from moodys diagram\n", + "\n", + "hf = f*L*V**2/(2*g*d) \n", + "\n", + "K= 0.5 # contraction constant\n", + "\n", + "hc = K*V**2/(2*g) \n", + "\n", + "K1 =10 # loss of the globe valve\n", + "\n", + "hg = K1*V**2/(2*g)\n", + "\n", + "Th = hf+hc+hg\n", + "\n", + "print \"Total head loss =\",round(Th,1),\"m\"\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Total head loss = 83.7 m\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.9 Page no 328" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# discharge through the pipe\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "# for water\n", + "\n", + "nu = 1.007*10**-6 # viscosity in m**2/s\n", + "\n", + "d1 = 0.3 # diameter of pipe 1 in m\n", + "\n", + "d2 = 0.15 # diameter of pipe 2 in m\n", + "\n", + "d3 = 0.08 # diameter of pipe 3 in m\n", + "\n", + "g = 9.81 # acclelration due to gravity in m/s**2\n", + "\n", + "e = 0.025 # for cast iron\n", + "\n", + "f1 = 0.019 # foe e/d1\n", + "\n", + "f2 = 0.022 # foe e/d2\n", + "\n", + "# Solution\n", + "\n", + "V3 = sqrt(2*g*100/((8.4*(f1)+268.85*(f2)+4.85)))\n", + "\n", + "V1 = (d3/d1)**2*V3\n", + "\n", + "V2 = (d3/d2)**2*V3\n", + "\n", + "# reynolds number for pipe BC\n", + "\n", + "R1 = V1*d1/nu\n", + "\n", + "R2 = V2*d2/nu\n", + "\n", + "Q = V3*pi*d3**2/4\n", + "\n", + "print \"Discharge through the pipe =\",round(Q,3),\"m**3/s\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge through the pipe = 0.067 m**3/s\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.10 Page no 332" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Replace the flow system\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "D = 0.2 # diameter of pipe 1\n", + "\n", + "D1 = 0.15 # diameter of pipe 2\n", + "\n", + "Q = 0.1 # discharge in m**3/s\n", + "\n", + "nu = 1.007*10**-6 # viscosity in m**2/s\n", + "\n", + "e = 0.025 # e for cast iron\n", + "\n", + "r = e/(100*D) \n", + "\n", + "# Solution\n", + "\n", + "V = Q/(pi*(0.2)**2/4)\n", + "\n", + "R = V*D/nu\n", + "\n", + "f = 0.021 # from moodys law\n", + "\n", + "r2 = e/D1\n", + "\n", + "V1 = Q/(pi*D1**2/4)\n", + "\n", + "R1 = V*D1/nu\n", + "\n", + "f2 = 0.023 # from moodys law\n", + "\n", + "L2 = 28562*D1**5/f2\n", + "\n", + "print \"Replacement of the flow system =\",round(L2,2),\"m\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Replacement of the flow system = 94.3 m\n" + ] + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.11 Page no 335" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Discharge through each branch\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "e = 0.025 # in cm\n", + "\n", + "nu = 1.007*10**-6 # viscosity in m**2/s\n", + "\n", + "Q1 = 0.5 # discharge in m**3/s\n", + "\n", + "D1 = 50\n", + "\n", + "L1 = 500 # length in m\n", + "\n", + "g = 9.81\n", + "\n", + "# Pipe 1\n", + "\n", + "r1 = e/D # r1 for pipe 1\n", + "\n", + "V1 = Q1/(pi*(0.5)**2/4)\n", + "\n", + "R = V*(0.5)/nu\n", + "\n", + "f1 = 0.018 # for the reynolds no\n", + "\n", + "hf1 = f*L1*V1**2/(2*g*D1)\n", + "\n", + "# pipe 2\n", + "\n", + "hf2 = hf1\n", + "\n", + "L2 =200 # length in m\n", + "\n", + "D2 = 0.3 # diameter in m\n", + "\n", + "r2 = e/D2\n", + "\n", + "f2 = 0.02 \n", + "\n", + "V2 = 0.419/sqrt(f2) \n", + "\n", + "R2 = V2*D2/nu\n", + "\n", + "Q2 = V2*(pi*D2**2/4)\n", + "\n", + "#pipe 3\n", + "\n", + "hf3=hf1\n", + "\n", + "L3 = 300 # length of pipe 3 in m\n", + "\n", + "D3 =0.15 # diameter of pipe 3 in m\n", + "\n", + "r3 = e/D3 \n", + "\n", + "f = 0.022 # from moody's law\n", + "\n", + "V3 = 0.242/sqrt(f2)\n", + "\n", + "R3 = V3*D3/nu\n", + "\n", + "Q3 = V3*(pi*D3**2/4)\n", + "\n", + "Td = Q1+Q2+Q3\n", + "\n", + "q1 = Q1*(2.0/Td)\n", + "\n", + "q2 = Q2*(2.0/Td)\n", + "\n", + "q3 = Q3*(2.0/Td)\n", + "\n", + "print \"Discharge through branch 1 =\",round(q1,2),\"m**3/s\"\n", + "\n", + "print \"Discharge through branch 2 =\",round(q2,3),\"m**3/s\"\n", + "\n", + "print \"Discharge through branch 3 =\",round(q3,3),\"m**3/s\"\n", + "\n", + "# Actual head loss\n", + "\n", + "d = 0.5\n", + "\n", + "v1 = q1/(pi*(d)**2/4)\n", + "\n", + "R4 = v1*d/nu\n", + "\n", + "r4 = 0.0005 # ratio of e/D\n", + "\n", + "f = 0.018\n", + "\n", + "Hf1 = f*L1*v1**2/(2*g*d)\n", + "\n", + "print \"It is found that hf1=hf2=hf3 =\",round(Hf1,1),\"The distribution od discharge is correct\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge through branch 1 = 1.35 m**3/s\n", + "Discharge through branch 2 = 0.566 m**3/s\n", + "Discharge through branch 3 = 0.082 m**3/s\n", + "It is found that hf1=hf2=hf3 = 43.5 The distribution od discharge is correct\n" + ] + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 3, + "metadata": {}, + "source": [ + "Example 9.14 Page no 349" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Find minimum depth below the ridge\n", + "\n", + "from math import *\n", + "\n", + "# Given\n", + "\n", + "e = 0.00015 # from moody's chart\n", + "\n", + "D = 2 # depth in ft\n", + "\n", + "r = e/D\n", + "\n", + "z1 = 100 # elevation in ft\n", + "\n", + "mu = 1.084*10**-5 # viscosity in Ns/ft**2\n", + "\n", + "p1 = 34 # pressure head in ft\n", + "\n", + "p2 = 10 # pressure head in ft\n", + "\n", + "g = 32.2 # acclelration due to gravity in ft/s**2\n", + "\n", + "L = 1000 # length in ft\n", + "\n", + "# Solution\n", + "\n", + "f = 0.011 # assume\n", + "\n", + "V = sqrt(100/(10000/(2*2*g)))/sqrt(f)\n", + "\n", + "R = V*D/mu\n", + "\n", + "V1 = 10.15\n", + "\n", + "f1 = 0.0125\n", + "\n", + "Q = V1*pi*D**2/4\n", + "\n", + "x = p1-p2-(V1**2/(2*g))-(f1*L*V1**2/(2*g*D))\n", + "\n", + "Dp = 30 - x\n", + "\n", + "print \"Minimum depth =\",round(Dp,2),\"ft\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Minimum depth = 17.6 ft\n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [], + "language": "python", + "metadata": {}, + "outputs": [] + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file |