summaryrefslogtreecommitdiff
path: root/Elementary_heat_Power/Chapter10.ipynb
blob: 3662d932297dc829072d384475ac4570abb32c6d (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
{
 "metadata": {
  "name": "EHP-10"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": "Drafts Fans Blowers and Compressors"
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": "Example 10.1,Page 515"
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nhb=29.0 #in of Hg\nsg=0.491 #specific gravity\nRa=53.3\nTa=460+40.0 #R\nTg=540+460.0 #R\nH=300 #ft\ngam=62.4 #lb/cu ft\n\n#calculations\npb=hb*sg*144\nrhoa=pb/(Ra*Ta) #density\nrhog=pb/(Ra*Tg) #density\ndp=H*(rhoa-rhog)\nD=dp/(gam)\n\n#results\nprint \"Theoretical draft in psf\",round(dp,3)\nprint \" Draft in ft H2O\",round(D,3)",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": "Theoretical draft in psf 11.541\n Draft in ft H2O 0.185\n"
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": "Example 10.2,Page 516"
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nmd=15.0 #lb per lb of coal\nx=0.1\nmss=1.0 #basis\nrea=29.0 #in of Hg\nsg=0.491\nR=53.3\nT=540+460.0 #R\nV=25.0 #fps\ngam=0.038 #lb/ft^3\n\n#calculations\nm=mss-mss*x+md\nms=m\nrhog=rea*0.491*144/(R*T) #density\nA=ms/(gam*V)\n\n#results\nprint \"stack area in sq ft\",round(A,3)",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": "stack area in sq ft 16.737\n"
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": "Example 10.3,Page 526"
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\np=144*29*0.491 #psf\nR=53.3\nT=70+460.0 #R\ngamw=62.4 #lb/ft^3\ngama=0.073 #lb/ft^3\nhw=3.0/12 #ft\nhw2=3.5/12 #ft\nhv=32.2 #ft/s^2\nms=9.0 #lb\ng=32.2 #ft/s^2\n\n#calculations\nrhoa=p/(R*T)\nhs=hw*gamw/gama\nht=hw2*gamw/gama\nhv=ht-hs\nV=sqrt(2*g*hv)\nmsv=ms*V*60.0\nmm=msv*gama\nairhp= ht*mm/33000.0\n\n#results\nprint \"Velocity head in ft of air\",round(hv,3)\nprint \" velocity of air in the duct in fps\",round(V,3)\nprint \" volume in cu ft per min\",round(msv,3)\nprint \" Mass flow rate in lb/min\",round(mm,3)\nprint \" Air horsepower in hp\",round(airhp,3)\nprint \"The answers in the textbook are a bit different due to rounding off error in the textbook Please use a calculator\"",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": "Velocity head in ft of air 35.616\n velocity of air in the duct in fps 47.893\n volume in cu ft per min 25861.99\n Mass flow rate in lb/min 1887.925\n Air hp in hp 14.263\nThe answers in the textbook are a bit different due to rounding off error in the textbook Please use a calculator\n"
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": "Example 10.4,Page 528"
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nA2=9.0 #sq ft\np2=3.0/12 *62.4 #psf\np1=-1.0/12 *62.4 #psf\nms=20000.0 #cfm\nA1=16.0 #sq ft\ngam=0.075 #lb/ft^3\ng=32.2 #ft/s^2\ninp=17.0 #hp\n\n#calculations\nV2=ms/60 *1/A2\nV1=ms/60 *1/A1\nht=(p2-p1)/gam +(V2**2 -V1**2)/(2*g)\nairhp=ht*ms*gam/33000\neta=airhp/inp *100\n\n#results\nprint \"Total head in ft of air\",round(ht,3)\nprint \" Air horsepower in hp\",round(airhp,3)\nprint \" Effifciency in percent\",round(eta,3)",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": "Total head in ft of air 291.894\n Air hp in hp 13.268\n Effifciency in percent 78.047\n"
      }
     ],
     "prompt_number": 12
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": "Example 10.5,Page 532"
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "\n# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nn1=400.0 #rpm\nmv1=10000.0 #lb\nmv2=15000.0 #lb\nh1=2.0 #in of water\nhp1=4.0 #hp\n\n#calculations\nn2=mv2/mv1 *n1 #speed\nh2=h1*(n2/n1)**2 #pressure height\nhp2=hp1 *(n2/n1)**3 #power\n\n#results\nprint  \"The speed in rpm\",round(n2,3)\nprint \" The pressure in inch of water\",round(h2,3)\nprint \" Power in hp\",round(hp2,3)",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": "The speed in rpm 600.0\n The pressure in in of water 4.5\n Power in hp 13.5\n"
      }
     ],
     "prompt_number": 13
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": "Example 10.6,Page 546"
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nm=100000.0 #lb/hr\np1=1.0 #psia\nx=0.8\np2=14.7 #psia\nt2=300.0 #F\n\n#calculations\n#\"from table A3 and A2\"\nh2=1192.8 #Btu/lb\nhf=69.7 #Btu/lb\nhfg=1036.3 #Btu/lb\nh1=hf+x*hfg\nW=h2-h1\npower=m*W\nhp=power/2545.0\n\n#results\nprint \"Power required in hp\",round(hp,3)",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": "11554.42 Power required in hp\n"
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": "Example 10.7,Page 546"
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\np1=14.7 #psia\nt1=60.0 #F\np2=60.0 #psia\nt2=440.0 #F\nm=10.0 #lb/sec\n\n#calculations\n#\"From mollier charts\"\nh2=216.3 #Btu/lb\nh1=124.3 #Btu/lb\nW21=h2-h1\npower=W21*m\nhp=power*3600/2545\ncp=0.237\nW212=cp*(t2-t1)\npower2=W212*m\nhp2=power2*3600.0/2545.0\n\n#results\nprint \"Power required in hp\",round(hp,3)\nprint \" Power required in hp\",round(hp2,3)\nprint \" Work done in Btu/lb\",round(W212,3)",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": "Power required in hp 1301.375\n Power required in hp 1273.933\n Work done in Btu/lb 90.06\n"
      }
     ],
     "prompt_number": 14
    }
   ],
   "metadata": {}
  }
 ]
}