diff options
author | Thomas Stephen Lee | 2015-09-04 22:04:10 +0530 |
---|---|---|
committer | Thomas Stephen Lee | 2015-09-04 22:04:10 +0530 |
commit | 41f1f72e9502f5c3de6ca16b303803dfcf1df594 (patch) | |
tree | f4bf726a3e3ce5d7d9ee3781cbacfe3116115a2c /Fluid_Mechanics/Chapter_25.ipynb | |
parent | 9c9779ba21b9bedde88e1e8216f9e3b4f8650b0e (diff) | |
download | Python-Textbook-Companions-41f1f72e9502f5c3de6ca16b303803dfcf1df594.tar.gz Python-Textbook-Companions-41f1f72e9502f5c3de6ca16b303803dfcf1df594.tar.bz2 Python-Textbook-Companions-41f1f72e9502f5c3de6ca16b303803dfcf1df594.zip |
add/remove/update books
Diffstat (limited to 'Fluid_Mechanics/Chapter_25.ipynb')
-rwxr-xr-x | Fluid_Mechanics/Chapter_25.ipynb | 328 |
1 files changed, 0 insertions, 328 deletions
diff --git a/Fluid_Mechanics/Chapter_25.ipynb b/Fluid_Mechanics/Chapter_25.ipynb deleted file mode 100755 index 9b08e6f3..00000000 --- a/Fluid_Mechanics/Chapter_25.ipynb +++ /dev/null @@ -1,328 +0,0 @@ -{ - "metadata": { - "name": "", - "signature": "sha256:2bebadfad4a18ede042a136ada8984e181f57290a23ad4ed8003a1f83445a447" - }, - "nbformat": 3, - "nbformat_minor": 0, - "worksheets": [ - { - "cells": [ - { - "cell_type": "heading", - "level": 1, - "metadata": {}, - "source": [ - "Chapter 25: Machine\u2013Network Interactions" - ] - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 25.4, Page 893" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from __future__ import division\n", - "import math\n", - "\n", - " #Initializing the variables\n", - "Pa_P1 = -200; # From previous Question\n", - "Q = 1.4311 ; # From previous questions.\n", - "\n", - " #Calculations\n", - "DpSys = Pa_P1 + 98.9*Q**2;\n", - "print \"System Operating point (m^3/s):\",round(DpSys,2)" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "System Operating point (m^3/s): 2.55\n" - ] - } - ], - "prompt_number": 1 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 25.7, Page 906" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from __future__ import division\n", - "import math\n", - "import sympy\n", - "from sympy import solve,symbols\n", - " \n", - "\n", - " #Initializing the variables\n", - "Vo = 25.3; #Outlet velocity\n", - "D = 10 ; # Mean hydraulic diameter\n", - "f = 0.008; # friction factor\n", - "X = 1000; # Length of road\n", - "P = 12600; # Absorbing power\n", - "Va = 300; # Tunnel air flow\n", - "K1 = 0.96;\n", - "K2 = 0.9;\n", - "T = 590; #Thrust\n", - "rho = 1.2; # Air density \n", - "\n", - " #Calculations\n", - "alpha = (1/D)**2;\n", - "A = math.pi*D**2/4; # Area of tunnel\n", - "Vt = Va/A;\n", - "W = Vo/Vt; #Omega\n", - "E = (1-alpha*W);\n", - "C = (1-alpha*W)*(1-E)**2 + E**2 - 1;\n", - " # Manipulating equation 25.20;\n", - "LHS = f*X*(E+1)**2/D + C + 1 ;\n", - "\n", - "n1 = symbols('n1')\n", - "result=solve(K1*(2*((alpha*W**2 + (1-alpha)*E**2-1)+(n1-1)*(alpha*W*(W-1)-C/2)))-LHS)\n", - "\n", - "n=result[0]\n", - "\n", - "\n", - " # Alternative approach using equation 25.22\n", - "n2 = (rho*((4*f*X*Vt**2)/(2*D) + 1.5*Vt**2/2))*A/(K1*K2*T); \n", - "Pt = round(n2)*P;\n", - "\n", - "print \"Number of fans required :\",round(n2)\n", - "print \"Total power consumed (KW) :\",Pt/1000" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "Number of fans required : 6.0\n", - "Total power consumed (KW) : 75.6\n" - ] - } - ], - "prompt_number": 1 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 25.8, Page 907" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from __future__ import division\n", - "import math\n", - "import sympy\n", - "from sympy import solve,symbols\n", - " #Initializing the variables\n", - "f = 0.008;\n", - "T = 290;\n", - "L = 750;\n", - "Dt = 9; # Diameter Tunnel\n", - "Df = 0.63; # Diameter fan\n", - "K1 = 0.98;\n", - "K2 = 0.92;\n", - "Vo = 27.9;\n", - "n = 10;\n", - "A=math.pi*Dt**2/4\n", - "rho=1.2\n", - "X=750\n", - " #Calculations\n", - "alpha = (Df/Dt)**2;\n", - " # equation 25.20 becomes when E = 1 nad C = 0\n", - "W=symbols('W')\n", - "omega = solve(2*K1* (alpha*W**2 +(n-1)*alpha*W*(W-1)) - 4*f*L/Dt -1)\n", - " \n", - "\n", - "for i in range(1,len(omega)): # since omega is always positive and real\n", - " if omega[i]>0:\n", - " w = round(omega[i],1);\n", - "Vt = Vo/w;\n", - "\n", - "# by equation 25.22\n", - "VT=(n*(K1*K2*T)/(A*(rho*((4*f*X)/(2*Dt) + 1.5/2))))**0.5\n", - "print \"Tunnel Velocity (m/s) :\",round(VT,2)" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "Tunnel Velocity (m/s) : 4.05\n" - ] - } - ], - "prompt_number": 2 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 25.9, Page 914" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from __future__ import division\n", - "import math\n", - "\n", - "\n", - " #Initializing the variables\n", - "Ws = 0.45;\n", - "Ks = 3.2;\n", - "H = 152;\n", - "h = 0;\n", - "Hatm = 10.3;\n", - "Pv = 350; #vapour pressure\n", - "g = 9.81;\n", - "rho = 1000;\n", - " \n", - " #Calculations\n", - "Ht1 = 152*(Ws/Ks)**(4/3); # the value of Ht1 is 11.12 and in book it is taken as 11.2 so there will be a difference in final answer\n", - "Hvap = round(Pv/(rho*g),3);\n", - "Z = Hatm -h -Hvap -Ht1;\n", - "print \"Elevation of pump (m):\",round(Z,3)\n", - " " - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "Elevation of pump (m): -0.851\n" - ] - } - ], - "prompt_number": 4 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 25.11, Page 927" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from __future__ import division\n", - "import math\n", - "import sympy\n", - "from sympy import symbols,solve\n", - "import numpy as np\n", - " \n", - "\n", - " #Initializing the variables\n", - "Co = 0;\n", - "Qc = 0.0024;\n", - "V = 5400;\n", - "c = 10;\n", - " #Calculations\n", - "#####--------------------PART(A)-----------------#######\n", - "n1=symbols('n1')\n", - "def partA(n1):\n", - " Ci = 10;\n", - " # t = infinity so e^(-nt) = 0\n", - " Q=10000*Qc/(c-Co)\n", - " n1 = Q*3600/V; \n", - " return n1\n", - "ans=partA(n1)\n", - "\n", - "print \"Part(A) : number of air changes per hour if the garage is in continuous use and the maximum permissible concentration of carbon monoxide is 0.1 per cent. :\",ans,\"\\n\"\n", - "\n", - "#####--------------------PART(B)-----------------#######\n", - "n=symbols('n')\n", - "def partB(n):\n", - " Ci = 0; \n", - " n=[1.5,1.2,0.9,1.0] \n", - " t=1 # time in hours\n", - " error=[]\n", - " mini=100\n", - " ans=0\n", - " for i in range(4): \n", - " Q = V/3600; \n", - " A = 10000*Qc/Q; # as Co=0 \n", - " error.append(abs((A*(1-math.e**(-n[i]*t))/c)-n[i]));\n", - " if(error[i]<mini):\n", - " mini=error[i]\n", - " ans=n[i]\n", - " return ans \n", - "ans=partB(n)\n", - "print \"Part(B) : number of air changes per hour if this maximum level is reached after 1 hour and the garage is out of use :\",ans,'\\n'\n", - " \n", - "#####--------------------PART(C)-----------------#######\n", - "c=symbols('c')\n", - "def partC(c):\n", - " Ci = 0;\n", - " n = 1; \n", - " t = 0.333333 # 20 minutes in hours\n", - " Q = V*n/3600;\n", - " y = (Co + 10000*Qc/Q)*(1-math.e**(-n*t)) + Ci*math.e**(-n*t) ; \n", - " return y\n", - "ans=partC(c)\n", - "print \"Part(C) :the concentration after 20 minutes (Parts per 10000) :\",round(ans,3),'\\n'\n", - "#####--------------------PART(D)-----------------#######\n", - "t=symbols('t')\n", - "def partD(t):\n", - " Ci = 10;\n", - " n = 1; \n", - " c = 0.1;\n", - " t=np.log(100) \n", - " return round(t,2)\n", - "ans=partD(t) \n", - "print \"Part(D) : time necessary to run the ventilation system at the rate calculated in (b) to reduce the concentration to 0.001 per cent (in hours) :\",ans,\"hours\"" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "Part(A) : number of air changes per hour if the garage is in continuous use and the maximum permissible concentration of carbon monoxide is 0.1 per cent. : 1.6 \n", - "\n", - "Part(B) : number of air changes per hour if this maximum level is reached after 1 hour and the garage is out of use : 1.0 \n", - "\n", - "Part(C) :the concentration after 20 minutes (Parts per 10000) : 4.535 \n", - "\n", - "Part(D) : time necessary to run the ventilation system at the rate calculated in (b) to reduce the concentration to 0.001 per cent (in hours) : 4.61 hours\n" - ] - } - ], - "prompt_number": 3 - } - ], - "metadata": {} - } - ] -}
\ No newline at end of file |