diff options
Diffstat (limited to 'Elementary_heat_Power')
| -rwxr-xr-x | Elementary_heat_Power/Chapter1.ipynb | 272 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter10.ipynb | 167 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter11.ipynb | 104 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter12.ipynb | 62 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter13.ipynb | 62 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter2.ipynb | 340 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter3.ipynb | 251 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter5.ipynb | 251 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter6.ipynb | 104 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter7.ipynb | 188 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter8.ipynb | 62 | ||||
| -rwxr-xr-x | Elementary_heat_Power/Chapter9.ipynb | 230 | ||||
| -rwxr-xr-x | Elementary_heat_Power/README.txt | 10 | ||||
| -rwxr-xr-x | Elementary_heat_Power/screenshots/chapter10.png | bin | 0 -> 157093 bytes | |||
| -rwxr-xr-x | Elementary_heat_Power/screenshots/chapter2.png | bin | 0 -> 116158 bytes | |||
| -rwxr-xr-x | Elementary_heat_Power/screenshots/chapter3.png | bin | 0 -> 86505 bytes |
16 files changed, 2103 insertions, 0 deletions
diff --git a/Elementary_heat_Power/Chapter1.ipynb b/Elementary_heat_Power/Chapter1.ipynb new file mode 100755 index 00000000..847dec63 --- /dev/null +++ b/Elementary_heat_Power/Chapter1.ipynb @@ -0,0 +1,272 @@ +{ + "metadata": { + "name": "EHP-1" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Matter and Energy" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.1, Page 17" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng=32.2 #ft/s^2\nm=500.0 #lb\nrate=10.0 #ft/s^2\n\n#calculations\nF1=m/g *rate\nms=m/g\nF2=ms*rate\n\n#results\nprint \"Force in case 1 in lbf\",round(F1,3)\nprint \"Force in case 2 in lbf\",round(F2,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Force in case 1 in lbf 155.28\nForce in case 2 in lbf 155.28\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.2,Page 18" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng=32.2 #ft/s^2\ng2=32.0 #ft/s^2\nrate=10 #ft/s^2\nw1=500 #lbf\n\n#calculations\nfd1=w1*g2/g\nF=fd1/g2 *rate\nms=w1/g\nF2=ms*rate\n#results\nprint \"Net weight of body in case 1 in lbf\",round(F,3)\nprint \"Force in case 2 in lbf\",round(F2,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Net weight of body in case 1 in lbf 155.28\nForce in case 2 in lbf 155.28\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.3, Page 19" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng=32.174 #ft/s^2\nm=500.0 #lbm\nrate=10.0 #ft/s^2\n\n#calculations\nF=1/g *m*rate\n\n#results\nprint \"Force required in lbf\",round(F,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Force required in lbf 155.405\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.4, Page 21" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng1=32.174 #ft/s^2\ngc=g1\ng2=30 #ft/s^2\nm=100 #lbm\n\n#calculations\nw1=g1/gc *m\nw2=g2/gc *m\n\n#results\nprint \"Weight in case 1 in lbf\",round(w1,3)\nprint \" Weight in case 2 in lbf\",round(w2,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Weight in case 1 in lbf 100.0\n Weight in case 2 in lbf 93.243\n" + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.5,Page 22" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nge=32.174 #ft/s^2\ngm=5.47 #ft/s^2\nwe=50 #lbm\n\n#calculations\nwm=we*gm/ge\n\n#results\nprint \"In case a, it will weigh the same, weight in lbm\",round(we,2)\nprint \" In case b, weight in lbf\",round(wm,2)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "In case a, it will weigh the same, weight in lbm 50.0\n In case b, weight in lbf 8.5\n" + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.6,Page 31" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng=32.2 #ft/s^2\np1=100 #psig\np2=29.0 #in of Hg\n\n#calculations\nBP=p2*0.491\nAP=BP+p1\n\n#results\nprint \"Absolute pressure in psia\",round(AP,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Absolute pressure in psia 114.239\n" + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.7,Page 32" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng=32.2 #ft/s^2\nPb=29.5 #in of Hg\nPv=10 #in of Hg\n\n#calculations\nAP=(Pb-Pv)*0.491\n\n#results\nprint \"Absoulte pressure in psia\",round(AP,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Absoulte pressure in psia 9.575\n" + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.8,Page 38" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng=32.2 #ft/s^2\nv1=1 #cu ft\np1=100 #psia\n\n#calculations\nv2=2*v1\nW=144*p1*(v2-v1)\n\n#results\nprint \"Work done in ft-lb\",round(W,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Work done in ft-lb 14400.0\n" + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.9,Page 39" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng=32.2 #ft/s^2\nv1=1 #cu ft\np1= 100.0 #psia\np2=50.0 #psia\nv2=3.0 #cu ft\n\n#calculations\npa=(p1+p2)/2\nW=pa*(v2-v1)*144\n\n#results\nprint \"Work done in ft-lb\",round(W,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Work done in ft-lb 21600.0\n" + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.10,Page 40" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng=32.2 #ft/s^2\np1=100.0 #psia\np2=25.0 #psia\nv2=2.0 #cu ft\n\n#calculations\nW=p1*144*v2 - p2*144*v2\n\n#results\nprint \"Work done in ft-lb\",round(W,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Work done in ft-lb 21600.0\n" + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.11,Page 41" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\ng=32.2 #ft/s^2\nn=100 #rpm\np1=100 #psia\np2=25 #psia\nv2=2 #cu ft\n\n#calculations\nW=p1*144*v2 - p2*144*v2\nHp=W*n/33000\n\n#results\nprint round(Hp,3),\"Horsepower developed in hp\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "65.0 Horsepower developed in hp\n" + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 1.12,Page 46" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "#initialisation of variable\nP=50.0 #hp\nm=30.0 #lb\nE=19000.0 #Btu/lb\n\n#calculations\neta= P*2545/(m*E) *100 #efficiency\n\n#results\nprint \"Efficiency in percent\",round(eta,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Efficiency in percent 22.325\n" + } + ], + "prompt_number": 12 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter10.ipynb b/Elementary_heat_Power/Chapter10.ipynb new file mode 100755 index 00000000..3662d932 --- /dev/null +++ b/Elementary_heat_Power/Chapter10.ipynb @@ -0,0 +1,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": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter11.ipynb b/Elementary_heat_Power/Chapter11.ipynb new file mode 100755 index 00000000..138cf155 --- /dev/null +++ b/Elementary_heat_Power/Chapter11.ipynb @@ -0,0 +1,104 @@ +{ + "metadata": { + "name": "EHP-11" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Feed Water Heaters and Condensers" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 11.1,Page 554" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nm1=1000.0 #lb/hr\nm2=5000.0 #lb/hr\nm3=3000.0 #lb/hr\n#From mollier charts\nh5=196.16 #Btu/lb\nh1=38.04 #Btu/lb\nh2=67.97 #Btu/lb\nh3=117.89 #Btu/lb\nh4=1156.3 #Btu/lb\n\n#calculations\nm4=(m1*h1+m2*h2+m3*h3-(m1+m2+m3)*h5)/(h5-h4) #mass rate\n\n#results\nprint \"Pounds of steam entering the heater in lb/hr\", round(m4,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Pounds of steam entering the heater in lb/hr 1076.801\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 11.2,Page 558" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nP1=100.0 #psia\nT1=400.0 #F\nT2=70.0 #F\n#\"From mollier charts\"\nh1=1227.6 #Btu/lb\nh2=298.4 #Btu/lb\nh3=279.9 #Btu/lb\nh4=38.04 #Btu/lb\n\n#calculations\nm1=(h3-h4)/(h1-h2)\n#results\nprint \"Mass of steam required in lb steam per lb water\",round(m1,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Mass of steam required in lb steam per lb water 0.26\n" + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 11.3,Page 567" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nh0=1260.0 #Btu/lb\nmsr=15.0 #lb\nm4=15.0 #lb per hr per kw\nt2=80.0 #F\nt3=60.0 #F\n\n#calculations\nh1=h0-3413/msr \n#\"from mollier charts\"\nh4=58 #Btu/lb\ndt=t2-t3\nm3=m4*(h1-h4)/dt\n\n#results\nprint \"enthalpy of steam entering the condenser in Btu/lb\",round(h1,3)\nprint \" mass of cooling water in lb per hr per kw\",round(m3,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "enthalpy of steam entering the condenser in Btu/lb 1032.467\n mass of cooling water in lb per hr per kw 730.85\n" + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 11.4,Page 568" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nm4=8*1000000 #lb per hr\ndt=12 #F\n\n#calculations\n#from mollier charts\")\ndh4=950 #Btu/lb\nm3=m4*(dh4)/dt #mass rate\n\n#results\nprint \" mass of cooling water in lb per hr\",round(m3,4)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": " mass of cooling water in lb per hr 633333333.0\n" + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter12.ipynb b/Elementary_heat_Power/Chapter12.ipynb new file mode 100755 index 00000000..c40e0935 --- /dev/null +++ b/Elementary_heat_Power/Chapter12.ipynb @@ -0,0 +1,62 @@ +{ + "metadata": { + "name": "EHP-12" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "The Gas Turbine Power Plant" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 12.1,Page 577" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nT1=80.0 #F\nT2=460.0 #F\nT3=1300.0 #F\nT4=780.0 #F\n\n#calculations\n#\"from mollier charts\"\nh1=129.1 #Btu/lb\nh2 = 221.2 #Btu/lb\nh3= 438.8 #Btu/lb\nh4 = 301.5 #Btu/lb\nwcom=h2-h1\nwcob=h3-h2\nwtur=h3-h4\neta=(wtur-wcom)/wcob *100\n\n#results\nprint \" work done by compressor in btu input as work per lb of air compressed\",round(wcom,3)\nprint \" Heat supplied in the combustor in Btu supplied per lb of air \",round(wcob,3)\nprint \" work done in the turbine in Btu output as work per lb of air\",round(wtur,3)\nprint \" Cycle efficiency in percent\",round(eta,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": " work done by compressor in btu input as work per lb of air compressed 92.1\n Heat supplied in the combustor in Btu supplied per lb of air 217.6\n work done in the turbine in Btu output as work per lb of air 137.3\n Cycle efficiency in percent 20.772\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 12.2,Page 579 " + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nT1=80 #F\nT2=460 #F\nT=700 #F\nT3=1300 #F\nT4=780 #F\n#\"from mollier charts\"\nh1=129.1 #Btu/lb\nh2 = 221.2 #Btu/lb\nh3= 438.8 #Btu/lb\nh4 = 301.5 #Btu/lb\nh=281.1 #Btu/lb\n\n#calculations\nwcom=h2-h1\nwcob=h3-h2\nwtur=h3-h4\noutput=-wcom+wtur\nQ=h3-h\neff=output/Q *100\n\n#results\nprint \" Heat supplied in the combustor in Btu supplied per lb of air \",round(Q,3)\nprint \" Cycle efficiency in percent\",round(eff,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": " Heat supplied in the combustor in Btu supplied per lb of air 157.7\n Cycle efficiency in percent 28.662\n" + } + ], + "prompt_number": 4 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter13.ipynb b/Elementary_heat_Power/Chapter13.ipynb new file mode 100755 index 00000000..9d4ff984 --- /dev/null +++ b/Elementary_heat_Power/Chapter13.ipynb @@ -0,0 +1,62 @@ +{ + "metadata": { + "name": "Chapter 13" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Mechanical Refrigeration" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 13.2,Page 596" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nmr=3.0 #lb\nmj=5.0 #lb\nt2=67.0 #F\nt1=60.0 #lb\nihp=7.25\n#\"From mollier charts\"\nh4=709 #Btu/b\nh3=618 #Btu/lb\n\n#calculations\nenergyin=ihp*2545/60\nenergyout=mr*(h4-h3) + mj*(t2-t1)\n\n#results\nprint \"Energy in Btu/min\",round(energyin,3)\nprint \"Energy out Btu/min\",round(energyout,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Energy in Btu/min 307.521\nEnergy out Btu/min 308.0\n" + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 13.3,Page 602" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nr=3.0 #lb\nhp=10.0 #hp\n\n#calculations\nh3=618 #Btu/lb\nh1=131 #Btu/lb\nQe=mr*(h3-h1)\nwork=hp*2545/60\ncop=Qe/work\n\n#results\nprint \"Coefficient of performance \",round(cop,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Coefficient of performance 3.444\n" + } + ], + "prompt_number": 11 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter2.ipynb b/Elementary_heat_Power/Chapter2.ipynb new file mode 100755 index 00000000..1b4fab13 --- /dev/null +++ b/Elementary_heat_Power/Chapter2.ipynb @@ -0,0 +1,340 @@ +{ + "metadata": { + "name": "EHP-2" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Fuels and Combustion" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.1,Page 76" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nx1=0.135 #moisture \nx2=0.056 #ash content\nveca=array([32.5, 48.4, 5.6, 13.5]) #materials as received\nB1=11788.0 #Btu/lb\n\n#calculations\nvecb=veca/(1-x1) #moisture free\nvecc=veca/(1-x1-x2) #moisture and ash free\nB2=B1/(1-x1)\nB3=B1/(1-x1-x2)\nvecb[3]=0\nvecc[3]=0\nvecc[2]=0\n\n#results\nprint \"In Moisture free case\",vecb\nprint \"In Moisture and Ash free case \",vecc\nprint \"Energy in Moisture free case Btu per lb\",round(B2,3)\nprint \"Energy in Moisture and ash free case Btu per lb\",round(B3,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "In Moisture free case [ 37.57225434 55.95375723 6.47398844 0. ]" + }, + { + "output_type": "stream", + "stream": "stdout", + "text": "\nIn Moisture and Ash free case [ 40.17305315 59.82694685 0. 0. ]\nEnergy in Moisture free case Btu per lb 13627.746\nEnergy in Moisture and ash free case Btu per lb 14571.08\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.2,Page 77" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\ny1=13.5\nx1=0.135 #moisture \nx2=0.056 #ash \nveca=array([66.0, 1.5, 1.1, 5.6, 5.9, 19.9]) # fram table 2.3\nvecb=array([66.0, 1.5, 1.1, 5.6, 5.9, 19.9,13.5]) #with moisture \n\n#calculations\nvecb[4] = vecb[4] - 1.0/9*y1 ;\nvecb[5] = vecb[5] - 8.0/9*y1 ;\nvecc=vecb/(1-x1)\nvecd=vecb/(1-x1-x2) # moisture and ash free\nvecd[3]=0\nvecd[6]=0\nvecc[6]=0\ns1=sum(vecc)\ns2=sum(vecd)\n\n#results\nprint \"With moisture as a separate item\",vecb\nprint \"In Moisture free case, \",vecc\nprint \"In Moisture and Ash free case\",vecd\nprint \"Total Mositure free content percent\",round(s1,3)\nprint \" Total Mositure and ash free content percent\",round(s2,3)\nprint \"slight error due to rounding off error\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "With moisture as a separate item [ 66. 1.5 1.1 5.6 4.4 7.9 13.5]\nIn Moisture free case, [ 76.30057803 1.73410405 1.2716763 6.47398844 5.0867052\n 9.13294798 0. ]\nIn Moisture and Ash free case [ 81.58220025 1.85414091 1.35970334 6.92212608 0. 9.76514215\n 0. ]\nTotal Mositure free content percent 100.0\n Total Mositure and ash free content percent 101.483\n" + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.3,Page 78" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nH=5.9 \nO=19.9 \nH2=4.4 \nO2=7.9 \n\n#calculations\nHa1=H-O/8\nHa2=H2-O2/8\n\n#results\nprint \"Available hydrogen in case 1 percent by weight\",round(Ha1,3)\nprint \" Available hydrogen in case 2 percent by weight\",round(Ha2,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Available hydrogen in case 1 percent by weight 3.413\n Available hydrogen in case 2 percent by weight 3.413\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.4,Page 93" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nH1=0.059\nO1=0.199\nH2=0.044\nO2=0.079\nC=0.66\nS=0.011\n\n#calculations\nQh1= 14600*C+62000*(H1-O1/8)+4050*S #heating value\nQh2=14600*C+62000*(H2-O2/8)+4050*S #heating value\n\n#results\nprint \"Heating value in case 1 Btu/lb \",round(Qh1,3)\nprint \" Heating value in case 2 Btu/lb \",round(Qh2,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Heating value in case 1 Btu/lb 11796.3\n Heating value in case 2 Btu/lb 11796.3\n" + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.5,Page 110" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nH1=0.059\nO1=0.199\nC=0.66\nS=0.011\n\n#calculations\nQh1= 11.52*C+34.56*(H1-O1/8)+4.32*S\n\n#results\nprint \"Theoretical air required in lb of air per lb of coal \",round(Qh1,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Theoretical air required in lb of air per lb of coal 8.83\n" + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.6,Page 110" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nmf=10000 #lb\nmr=700 #lb\nCr=0.20\nCco2=14.1 #% CO2\nCo2=5.1 #% O2\nCco=0.1 #% CO\nCf=0.66 \n\n#calculations\nCn2=100-(Cco2+Co2+Cco) #% N2\nCi=mf*Cf #carbon input\nCa=mr*Cr #carbon in ashpit\nCb=(Ci-Ca)/mf\nCb2=((mf*Cf)-mr*Cr)/(mf)\nMco2=Cco2 *44 #mass of CO2\nMo2=Co2 *32 #mass of O2\nMco=Cco *28 #mass of CO\nMn2=Cn2 *28 #mass of N2\nsumvec=Mo2+Mco2+Mco+Mn2;\nLbc=Cco2*12 + Cco*12\nGc=sumvec/Lbc\nGf=Gc*Cb\n\n#results\nprint \"Carbon in the dry products combustion in lb per lb of fuel\",round(Cb,3)\nprint \" In case 2, Carbon in the dry products combustion in lb per lb of fuel\",round(Cb2,3)\nprint \" Dry gaseous products of combstion per lb of coal in lb \",round(Gf,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Carbon in the dry products combustion in lb per lb of fuel 0.646\n In case 2, Carbon in the dry products combustion in lb per lb of fuel 0.646\n Dry gaseous products of combstion per lb of coal in lb 11.548\n" + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.7,Page 111" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nmf=10000.0 #lb\nmr=700.0 #lb\nCr=0.20\nCco2=14.1\nCo2=5.1\nCco=0.1\nCf=0.66\n\n#calculations\nCn2=100-(Cco2+Co2+Cco)\nCi=mf*Cf\nCa=mr*Cr\nCb=(Ci-Ca)/mf\nCb2=((mf*Cf)-mr*Cr)/(mf)\nMco2=Cco2 *44 #mass of CO2\nMo2=Co2 *32 #mass of O2\nMco=Cco *28 #mass of CO\nMn2=Cn2 *28 #mass of N2\nsumvec=Mn2+Mco2+Mco+Mn2;\nCbb1=Cb*Cco*12/(Cco2*12 + Cco*12)\nCbb2= Cb*(Cco /(Cco + Cco2))\n\n#results\nprint \"In case 1, Carbon burned per lb of fuel in lb per lb of fuel\",round(Cbb1,6)\nprint \" In case 2, Carbon burned per lb of fuel in lb per lb of fuel\",round(Cbb2,6)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "In case 1, Carbon burned per lb of fuel in lb per lb of fuel 0.004549\n In case 2, Carbon burned per lb of fuel in lb per lb of fuel 0.004549\n" + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.8,Page 111" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nH=4.4/100 #hydrogen\nM=13.5/100 #moisture\nH2=0.059 #hydrogen\n\n#calculations\npro=M+9*H\npro2=9*H2\n\n#results\nprint \"In case 1, watervapor present in products in lb\",round(pro,3)\nprint \"In case 2, watervapor present in products in lb\",round(pro2,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "In case 1, watervapor present in products in lb 0.531\nIn case 2, watervapor present in products in lb 0.531\n" + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.9,Page 113" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nGf=11.57;#lb\nf=11.57 #lb per lb of fuel\nH=4.4/100\nM=13.5/100\nmr=700.0 #lb\nmf=10000.0 #lb\nmc=1 #lb\n\n#calculations\npro=M+9*H\nmrf=mr/mf\nAa=Gf+pro+mrf-mc\n\n#results\nprint \"Actual air supplied in lb of air supplied per lb of fuel\",round(Aa,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Actual air supplied in lb of air supplied per lb of fuel 11.171\n" + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.10,Page 114" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nGf=11.57 #lb per lb of fuel\nH=4.4/100\nM=13.5/100\nmr=700.0\nmf=10000.0\nmc=1.0 #lb\n\n#calculations\npro=M+9*H\nmrf=mr/mf\nAa=Gf+pro+mrf-mc\nAt=8.83\nea=(Aa-At)/At *100\n\n#results\nprint \"Excess air in percent\",round(ea,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "26.512 Excess air in percent\n" + } + ], + "prompt_number": 33 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.11a,Page 114" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nGf=11.57 #lb per lb of fuel\ntg=500 #F\nta=70 #F\n#calculations\n\nQ1=0.24*Gf*(tg-ta)\n\n#results\nprint \"Heat loss in Btu per lb of fuel\",round(Q1,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Heat loss in Btu per lb of fuel 1194.024\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.11b,Page 115" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nCo=0.1\nCo2=14.1\nCb=0.646\n\n#calculations\nQ2=Co/(Co+Co2) *Cb*10160\n\n#results\nprint \"Heat loss in Btu per lb of fuel\",round(Q2,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "46.221 Heat loss in Btu per lb of fuel\n" + } + ], + "prompt_number": 35 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.11c ,Page 115" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nmf=10000.0 #lb \nmr=700.0 #lb\nCr=0.2\n\n#calculations\nQ3=mr*Cr/mf *14600\n\n#results\nprint \"Heat loss in Btu per lb of fuel\",round(Q3,2)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "204.4 Heat loss in Btu per lb of fuel\n" + } + ], + "prompt_number": 37 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.11d,Page 115" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nM=0.135\ntg=500.0 #F\nta=70.0 #F\n\n#calculations\nQ4=M*(1089+0.46*tg-ta) #heat loss\n\n#results\nprint \"Heat loss in Btu per lb of fuel\",round(Q4,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "168.615 Heat loss in Btu per lb of fuel\n" + } + ], + "prompt_number": 39 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 2.11e,Page 115" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nPer=0.044 #percentage\ntg=500 #F\nta=70 #F\n\n#calculations\nQ5=9*Per*(1089+0.46*tg-ta) #heat loss\n\n#results\nprint \"Heat loss in Btu per lb of fuel\",round(Q5,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Heat loss in Btu per lb of fuel 494.604\n" + } + ], + "prompt_number": 2 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter3.ipynb b/Elementary_heat_Power/Chapter3.ipynb new file mode 100755 index 00000000..cc592814 --- /dev/null +++ b/Elementary_heat_Power/Chapter3.ipynb @@ -0,0 +1,251 @@ +{ + "metadata": { + "name": "EHP-3" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Internal Combustion Engines" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.1 , Page 131" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nre=6.0\nk=1.4\n\n#calculations\nnt=1-1/re**(k-1)\nntt=nt*100 #efficiency\n\n#results\nprint \"Thermal efficiency in percent\",round(ntt,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Thermal efficiency in percent 51.164\n" + } + ], + "prompt_number": 16 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.2,Page 132" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nd=3.25 #in\nstroke=4.0 #in\nv=6.0 #cu in\n\n#calculations\nDp=d**2 *pi*stroke/4\nr=(Dp+v)/v #compression ratio\n\n#results\nprint \"Compression ratio \",round(r,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Compression ratio 6.531\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.3,Page 133" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nper=20.0\nDp=100.0\n\n#calculations\nr=Dp/per +1\n\n#results\nprint \"Compression ratio \",round(r,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Compression ratio 6.0\n" + } + ], + "prompt_number": 18 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.4,Page 136" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nr=16.0\nrc=4.0\nk=1.4\n\n#calculations\netat=1-1/r**(k-1) *((rc**k -1)/(k*(rc-1)))\neta=etat*100; #efficiency\n\n#results\nprint \"Thermal efficiency in percent\",round(eta,3)\nprint \"The answer is a bit different due to rounding off error in the textbook\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Thermal efficiency in percent 53.154\nThe answer is a bit different due to rounding off error in the textbook\n" + } + ], + "prompt_number": 19 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.5,Page 142" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nF=200 #lb\narea=1.65 #sq. in\nlength=3.5 #in\n\n#calculations\nord=area/length\nmep=ord*F\n\n#results\nprint \"MEP of an engine in psi\",round(mep,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "MEP of an engine in psi 94.286\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.6,Page 143" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nPi=90 #psi\nL=5.0/12.0 #ft\nr=5.0 #in\nx=1.5 #ft\nrpm=1500.0 #rpm\n\n#calculations\nA=pi*x*x\nN=rpm*4/2\nIhp=Pi*L*A*N/33000.0\n\n#results\nprint \"IHP of cylinder \",round(Ihp,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "IHP of cylinder 24.097\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.7,Page 144" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nr=4.0 #ft\nn=300.0 #rpm\nF=60.0 #lb\n\n#calculations\nBhp=2*pi*r*F*n/33000\n\n#results\nprint \"Bhp of the engine \",round(Bhp,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Bhp of the engine 13.709\n" + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.8,Page 147" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nC=1.0/4000.0\nF=125.0 #lb\nn=3500 #rpm\n\n#calculations\nBhp=F*n*C\n\n#results\nprint \"Bhp of the engine \",round(Bhp,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Bhp of the engine 109.375\n" + } + ], + "prompt_number": 23 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.9,Page 151" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nr=1.75 #ft\nF1=72.0 #lb\nF2=24.0 #lb\nn=500.0 #rpm\nm=1.8 #lb\nmi=15.0 #min\nQh=20000.0 #Btu/lb\n\n#calculations\nBhp=2*pi*r*F1*n/33000.0\nFhp=2*pi*r*F2*n/33000.0\nIhp=Bhp+Fhp\nFc=m*60/mi\nBsfc=Fc/Bhp\nIsfc=Fc/Ihp\netam=Bhp/Ihp *100\netabt=Bhp*2545/(Fc*Qh)\netait=Ihp*2545/(Fc*Qh)\n\n#results\nprint \"Thermal efficiency in percent\",round(etam,3)\nprint \" Brake thermal effficiency in percent\",round(etabt*100,3)\nprint \" Indicated thermal effficiency in percent\",round(etait*100,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Thermal efficiency in percent 75.0\n Brake thermal effficiency in percent 21.2\n Indicated thermal effficiency in percent 28.266\n" + } + ], + "prompt_number": 24 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.10,Page 191" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nbore=3.0 #in\nstroke=4.0 #in\nrpm=3000.0 #rpm\nair=110.0 #cu ft per min\n\n#calculations\npdv=bore*bore*pi*stroke*2*bore/4\npde=pdv*rpm /2\nreq=air*1728\neff=req/pde *100; #efficiency\n\n#results\nprint \"Volumetric efficiency in percent\",round(eff,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Volumetric efficiency in percent 74.697\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 3.11, Page 244" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nx1=11.5 # CO2\nx2=4.1 #CO\nx3=0.4 #O2\nx4=2.3 # H2\nx5=0.2 #CH4\nx6=81.5 #N2\nyc=0.842 #lb of C \nyh=0.158 #lb of H2\nbasis=1\nbhp=100\nburn=8.9 #gal/hr\nsg=0.731\nQh=20750.0 #Btu/lbm\nrate=66.0 #gpm\nex=1100.0 #F\nair=70.0 #F\ncp=0.254\nh2=4330.0 #btu/lb\nh4=62000.0 #btu/lb \nh5=23700.0 #btu/lb\n\n#calculations\nc1=x1*44 #mass CO2\nc2=x2*28 #mass CO\nc3=x3*32 #mass O2\nc4=x4*2 #mass H2\nc5=x5*16 #mass CH4\nc6=x6*28 #mass N2\nsumm=c1+c2+c3+c4+c5+c6\ncarbon=x1*12 + x2*12+x5*12\nhydrogen=x4*2+x5*4\nlbdrygas=summ/carbon *yc\nlbfuel=carbon/yc\nlbH=lbfuel*yh\nlbH2=lbH-hydrogen\nlb3=lbH2*9\nlbwater=lb3/lbfuel\nlbair=lbdrygas+lbwater-basis\nbsfc=burn*sg*8.33/bhp\nfuelmin=bsfc*bhp/60\nenergy=2545/bsfc\nper1=energy/Qh #energy supplied\nEc=rate*8.33*10\nEclb=Ec/fuelmin\nper2=Eclb/Qh #energy absorbed\ndryloss=(ex-air)*cp*lbdrygas\nper3=dryloss/Qh #energy lost\nhv2=h2*c2/lbfuel\nhv4=h4*c4/lbfuel\nhv5=h5*c5/lbfuel\nhv=hv2+hv4+hv5\nper4=hv/Qh #energy supplied\neh2=lbwater*(1066+0.5*ex-air)\nper5=eh2/Qh #energy supplied\nrad=1017 \nper6=rad/Qh; #energy supplied\n\n#results\nprint \"Air supplied per lb of fuel in lb air per lb fuel\",round(lbair,3)\nprint \"Percentage of energy supplied utilized in Btu in percent\",round(per1*100,3)\nprint \"Percentage of energy absorbed by coolant in percent\",round(per2*100,3)\nprint \"Energy lost in sensible heat in percent\",round(per3*100,3)\nprint \"Energy supplied in combustiles in exhaust in percent\",round(per4*100,3)\nprint \"Energy supplied in water formed by combustion in percent\",round(per5*100,3)\nprint \"Energy supplied unaccounted for in percent\",round(per6*100,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Air supplied per lb of fuel in lb air per lb fuel 13.189\nPercentage of energy supplied utilized in Btu in percent 22.632\nPercentage of energy absorbed by coolant in percent 29.334\nEnergy lost in sensible heat in percent 16.369\nEnergy supplied in combustiles in exhaust in percent 18.366\nEnergy supplied in water formed by combustion in percent 8.987\nEnergy supplied unaccounted for in percent 4.901\n" + } + ], + "prompt_number": 27 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter5.ipynb b/Elementary_heat_Power/Chapter5.ipynb new file mode 100755 index 00000000..89ee8fb7 --- /dev/null +++ b/Elementary_heat_Power/Chapter5.ipynb @@ -0,0 +1,251 @@ +{ + "metadata": { + "name": "EHP-5" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Steam Generation" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.1,Page 296" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nx=0.98\nvg=26.80\nvf=0.01672\n\n#calculations\nvx=x*vg+(1-x)*vf #specific wet volume\n\n#results\nprint \"Specific volume of wet steam in cu ft per lb\",round(vx,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Specific volume of wet steam in cu ft per lb 26.264\n" + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.2,Page 298" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nhf=167.99 #Btu/lb\nhg=4.5 #Btu/lb\n\n#calculations\nhc=hf+hg\n\n#results\nprint \"Enthalpy of water in Btu/lb\",round(hc,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Enthalpy of water in Btu/lb 172.49\n" + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.3,Page 300" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nx=0.97\nhg=1187.2 #Btu/lb\nhf=298.40 #Btu/lb\nhfg=888.8 #Btu/lb\n\n#calculations\nhx1=x*hg+(1-x)*hf\nhx2=hf+x*hfg\nhx3=hg-(1-x)*hfg\n\n#results\nprint \" In case 1, enthalpy in Btu/lb\",round(hx1,3)\nprint \" In case 2, enthalpy in Btu/lb\",round(hx2,3)\nprint \" In case 3, enthalpy in Btu/lb\",round(hx3,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": " In case 1, enthalpy in Btu/lb 1160.536\n In case 2, enthalpy in Btu/lb 1160.536\n In case 3, enthalpy in Btu/lb 1160.536\n" + } + ], + "prompt_number": 23 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.4,Page 302" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nh1=1172.0 #Btu/lb\nhf1=355.36 #Btu/lb\nhfg1=843.0 #Btu/lb\n\n#calculations\nh2=h1\nx1= (h2-hf1)/hfg1\n\n#results\nprint \"Quality of steam in percent\",round(x1*100,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "96.873 Quality of steam in percent\n" + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.7,Page 306" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nP=200 #psia\nT=260 #F\n\n#calculations\n#\"From mollier chart\"\nhx=1174 #Btu/lb\nx1=2.8 #percent\ny1=100-x1 #percent\n\n#results\nprint \"Quality in percent\",round(y1,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Quality in percent 97.2\n" + } + ], + "prompt_number": 26 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.8,Page 306" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nP=200 #psia\nT=500 #F\n\n#calculations\n#\"From mollier chart\"\nhi=1269 #Btu/lb\nhf=1063 #Btu/lb\ndh=hi-hf\ny1=91\n\n#results\nprint \"Quality in percent\",round(y1,3)\nprint \" Change in enthalpy in Btu/lb\",round(dh,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Quality in percent 91.0\n Change in enthalpy in Btu/lb 206.0\n" + } + ], + "prompt_number": 27 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.9,Page 307" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nP=200 #psia\nTs=260 #F\nTf=220 #F\nm=10000 #lb\nPc=20 #psia\n\n#calculations\n#\"From mollier charts\"\nhf=188.0 #Btu/lb\nh2=1172.0 #Btu/lb\nQ=m*(h2-hf)\n\n#results\nprint \"Heat absorption in Btu/hr\",round(Q,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "9840000.0 Heat absorption in Btu/hr\n" + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.10,Page 307" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nm=150000.0 #lb\nP1=1000.0 #psia\nTs=900.0 #F\nTf=200.0 #F\n\n#calculations\n#\"From mollier charts\"\nh2=1448.2 #Btu/lb\nhf=167.99 #Btu/lb\ncorrec=2.2 #Btu/lb\nhc=hf+correc\nQ=m*(h2-hc)\n\n#results\nprint \"Heat absorption in Btu/hr\",round(Q,3)\nprint\"The answer is a bit different due to rounding off error in textbook\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Heat absorption in Btu/hr 191701500.0\nThe answer is a bit different due to rounding off error in textbook\n" + } + ], + "prompt_number": 28 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.11,Page 308" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nm=150000.0 #lb\nP1=1000.0 #psia\nTs=900.0 #F\nTf=200.0 #F\n\n#calculations\n#\"From mollier charts\"\nh2=1448.2 #Btu/lb\nhf=167.99 #Btu/lb\ncorrec=2.2 #Btu/lb\nhc=hf+correc\nQ=m*(h2-hc)\noutput=Q/1000\n\n#results\nprint \"Output of the steam generating unit kB/hr\",round(output,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Output of the steam generating unit kB/hr 191701.5\n" + } + ], + "prompt_number": 29 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.12,Page 308" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nm=150000 #lb\nP1=1000 #psia\nTs=900 #F\nTf=200 #F\nm2=21000 #lb\nHV=12000 #Btu/lb\n\n#calculations\n#\"From mollier charts\"\nh2=1448.2 #Btu/lb\nhf=167.99 #Btu/lb\ncorrec=2.2 #Btu/lb\nhc=hf+correc\nQ=m*(h2-hc)\noutput=Q\ninpu=m2*HV\neta=output/inpu\n\n#results\nprint \"Efficiency of the steam generating unit in percent\",round(eta*100,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Efficiency of the steam generating unit in percent 76.072\n" + } + ], + "prompt_number": 30 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 5.13,Page 309" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nhv=11780 #Btu/lb\nsteam=55000 #lb/hr\ncoal=6480 #lb\nx1=0.66 # %Carbon\nx2=0.044 #% Hydrogen\nx3=0.079 #% Oxygen\nx4=0.015 #% Nitrogen\nx5=0.11 # % Sulphur\nz1=14.5 # % CO2\nz2=0.2 # CO\nz3=4.4 # O2\nz4=80.9 # N2\nxash=0.076 #% ash\nxmois=0.115 # % moisture\nyc=0.21 # % Carbon in refuse\nrefuse=622 #lb/hr\ncp=0.24 #coeff of heat\ntg=400 #F\nta=70 #F\nQco=10160 #Btu/lb\nQc=14600 #Btu/lb\n\n#calculations\n#\"From steam tables\"\nhf=269.6 #Btu/lbm\nhfg=1.5 #Btu/lbm\nh1=hf+hfg #Btu/lbm\nh2=1196.5\nQb=h2-h1\nh3=1407.7 #Btu/lbm\nQs=h3-h2\nh4=h3-h1\nout=steam*h4/1000\neff=steam*h4/(coal*hv)\n#Energy balance\nCi=coal*x1\nCr=refuse*yc\nCb=(Ci-Cr)/coal\nlbt= z1*44+z2*28+z3*32+z4*28\nlbC=z1*12+z2*12\ndry=lbt/lbC *Cb\nloss1=dry*cp*(tg-ta)\nloss2=z2*12/(lbC) *Cb*Qco\nloss3=Cr*Qc/coal\nloss4=xmois*(1089+0.46*tg-ta)\nloss5=x2*9*(1089+0.46*tg-ta)\nloss6=steam*h4/coal\n\n#results\nprint \"Heat absorbed in the boiler in Btu per lb of steam generated\",round(Qb,3)\nprint \" Heat absorbed in the superheater in Btu/lb of steam\",round(Qs,3)\nprint \" Heat absorbed in steam generating in Btu/lb of steam generated\",round(h4,3)\nprint \" Output of steam generating unit in kB\",round(out,3)\nprint \" Efficiency of steam generating unit in percent\",round(eff*100,3)\nprint \" Carbon burned to CO and CO2 in lb of C per lb of fuel\",round(Cb,3)\nprint \" Dry products of combustion in lb per lb of fuel\",round(dry,3)\nprint \" Loss due to sensible heat in dry gaseous products of combustion in Btu/lb of fuel\",round(loss1,3)\nprint \" Loss due to CO in dry products of combustion in Btu/lb of fuel\",round(loss2,3)\nprint \" Loss due to C in refuse in Btu/lb of fuel\",round(loss3,3)\nprint \" Loss due to evaporating moisture in fuel in Btu/lb of fuel\",round(loss4,3)\nprint \" Loss due to water vapor formed from H in Btu/lb of fuel\",round(loss5,3)\nprint \" Energy absorbed in generating steam in Btu/lb of fuel\",round(loss6,3)\nprint \"The answers are a bit different due to rounding off error in the textbook\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Heat absorbed in the boiler in Btu per lb of steam generated 925.4\n Heat absorbed in the superheater in Btu/lb of steam 211.2\n Heat absorbed in steam generating in Btu/lb of steam generated 1136.6\n Output of steam generating unit in kB 62513.0\n Efficiency of steam generating unit in percent 81.894\n Carbon burned to CO and CO2 in lb of C per lb of fuel 0.64\n Dry products of combustion in lb per lb of fuel 11.062\n Loss due to sensible heat in dry gaseous products of combustion in Btu/lb of fuel 876.078\n Loss due to CO in dry products of combustion in Btu/lb of fuel 88.446\n Loss due to C in refuse in Btu/lb of fuel 294.298\n Loss due to evaporating moisture in fuel in Btu/lb of fuel 138.345\n Loss due to water vapor formed from H in Btu/lb of fuel 476.388\n Energy absorbed in generating steam in Btu/lb of fuel 9647.068\nThe answers are a bit different due to rounding off error in the textbook\n" + } + ], + "prompt_number": 31 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter6.ipynb b/Elementary_heat_Power/Chapter6.ipynb new file mode 100755 index 00000000..b2d3e81f --- /dev/null +++ b/Elementary_heat_Power/Chapter6.ipynb @@ -0,0 +1,104 @@ +{ + "metadata": { + "name": "EHP-6" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Steam Power Plant Cycles" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.1,Page 366" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nP1=200.0 #psia\nT1=600.0 #F\nP2=2.0 #psia\nJ=778.0\n\n#calculations\n#\"from mollier charts,\"\nh1=1322 #Btu/lb\nh2=974 #Btu/lb\nvf2=0.01623 #cu ft per lb\nhf2=94.0 #Btu/lb\nt2=126.0 #F\nWtj=h1-h2\nQout=h2-hf2\nWp=(P1-P2)*vf2\nWpj=Wp/J\nh3=hf2+Wpj\nQin=h1-h3\netat=((h1-h2)-Wpj)/(h1-(hf2+Wpj)) #efficiency\neta=((h1-h2))/(h1-(hf2)) #efficiency\n\n#results\nprint round(etat*100,3),\"Efficiency of rankine cycle in percent\"\nprint round(eta*100,3),\" Efficiency of rankine cycle neglecting boiler feed pump in percent\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "28.339 Efficiency of rankine cycle in percent\n28.339 Efficiency of rankine cycle neglecting boiler feed pump in percent\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.2,Page 366" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nB=70 #F\nP1=140 #psia\nx=0.986\nP2=14.7 #psia\nms=2000.0 #lb/hr\nIhp=80.0\n\n#calculations\n#\"From mollier charts\"\nhc=38 #Btu/lb\nhf=324.82 #Btu/lb\nhfg=868.2 #Btu/lb\nh1=hf+x*hfg\nQin=ms*(h1-hc)\neta=Ihp*2545*100/(Qin) #efficiency\nQw=Ihp*2545 \nQr=Qin-Qw\nper=Qr/Qin *100 #percent\n\n#results\nprint \"Heat input to the boiler in Btu/hr\",round(Qin,3)\nprint \" Cycle efficiency in percent\",round(eta,3)\nprint \" Heat rejected to waste in Btu/hr is\",round(Qr,3),\" or in percent of Qin\",round(per,3)\nprint \"The answer is a bit different due to rounding off error in textbook\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Heat input to the boiler in Btu/hr 2285730.4\n Cycle efficiency in percent 8.907\n Heat rejected to waste in Btu/hr 2082130.4 or in percent of Qin 91.093\nThe answer is a bit different due to rounding off error in textbook\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.3,Page 372" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nB=70.0 #F\nP1=140.0 #psia\nx=0.986\nP2=14.7 #psia\nms=2000 #lb/hr\nIhp=80\n\n#calculations\n#\"From mollier charts\"\nhc=180.0 #Btu/lb\nhf=324.82 #Btu/lb\nhfg=868.2 #Btu/lb\nh1=hf+x*hfg\nQin=ms*(h1-hc)\neta=Ihp*2545*100/(Qin)\nQw=Ihp*2545 #Btu/hr\nQr=Qin-Qw #Btu/hr\nper=Qr/Qin *100 #percent\n\n#results\nprint \"Heat input to the boiler in Btu/hr\",round(Qin,3)\nprint \" Cycle efficiency in percent\",round(eta,3)\nprint \" Heat rejected to waste in Btu/hr is\",round(Qr,3),\" or in percent of Qin\",round(per,3)\nprint \"The answer is a bit different due to rounding off error in textbook\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Heat input to the boiler in Btu/hr 2001730.4\n Cycle efficiency in percent 10.171\n Heat rejected to waste in Btu/hr 1798130.4 in percent of Qin 89.829\nThe answer is a bit different due to rounding off error in textbook\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 6.4,Page 380" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nm=1.24 #lb\nHV=11300 #Btu/lb\n\n#calculations\nHR=m*HV\neff=3413/HR #efficiency\n\n#results\nprint \"Plant heat rate in Btu/kw hr\",round(HR,3)\nprint \" Overall efficiency in percent\",round(eff*100,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Plant heat rate in Btu/kw hr 14012.0\n Overall efficiency in percent 24.358\n" + } + ], + "prompt_number": 10 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter7.ipynb b/Elementary_heat_Power/Chapter7.ipynb new file mode 100755 index 00000000..b3a3267c --- /dev/null +++ b/Elementary_heat_Power/Chapter7.ipynb @@ -0,0 +1,188 @@ +{ + "metadata": { + "name": "EHP-7" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Steam Turbines" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 7.1,Page 390" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nP1=200.0 #psia\nT1=500.0 #psia\nm=1.0 #lb /s\nP4=140.0 #psia\nP11=1.0 #psia\nx=0.808\n\n#calculations\n#\"From mollier charts\"\nh1=1268.9 #Btu/lb\nh4=1234.7 #Btu/lb\nV4=223.8*sqrt(h1-h4)\nv4=3.584 #cu ft/lb\nA4=m*v4/V4\nh11=907.4 #Btu/lb\nV11=223.8*sqrt(h1-h11)\nvf=0.01614 #cu ft/lb\nvg=333.6 #cu ft/lb\nvfg=vg-vf\nv11=x*vg\nA11=m*v11/V11 #area\n\n#results\nprint \"Area of nozzle in sq ft\",round(A4,6)\nprint \" Area of nozzle in sq ft\",round(A11,6)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Area of nozzle in sq ft 0.002738\n Area of nozzle in sq ft 0.063347\n" + } + ], + "prompt_number": 18 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 7.4,Page 410" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nP1=200.0 #psia\nT1=500.0 #F\nP2=1.0 #psia\nalpha=20.0 #degrees\nn=3600.0\ng=32.2 #ft/s^2\n\n#calculations\n#\"From mollier charts\"\nV1=4240.0 #fps\nVb=V1*cos(alpha*pi/180) /2\nR=Vb*60/(n*2*pi)\nwork=1/32.2 *(V1*cos(alpha*pi/180))*Vb\neff=work/(V1**2 /(2*g)) *100 #efficiency\n\n#results\nprint \"Blade velocity in fps\",round(Vb,3)\nprint \" Blade radius in ft\",round(R,3)\nprint \" Work done in per lb of steam\",round(work,3)\nprint \" Blade efficiency in percent\",round(eff,3)\nprint \"The answers are a bit different due to rounding off error in textbook\"\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Blade velocity in fps 1992.148\n Blade radius in ft 5.284\n Work done in per lb of steam 246500.315\n Blade efficiency in percent 88.302\nThe answers are a bit different due to rounding off error in textbook\n" + } + ], + "prompt_number": 19 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 7.5,Page 411" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nP1=200.0 #psia\nT1=500.0 #F\nP2=1.0 #psia\nalpha=20.0 #degrees\nn=3600.0\ng=32.2 #ft/s^2\nVb=1200.0 #fps\n\n#calculations\n#\"From mollier charts\"\nV1=4240.0 #fps\nV1x=3980.0 #fps\nV2x=-1580.0 #fps\nwork=1/32.2 *(V1x - V2x)*Vb\neff=work/(V1**2 /(2*g)) *100 #efficiency\n\n#results\nprint \" Work done in per lb of steam\",round(work,3)\nprint \" Blade efficiency in percent\",round(eff,3)\nprint \"The answers are a bit different due to rounding off error in textbook\"\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": " Work done in per lb of steam 207204.969\n Blade efficiency in percent 74.226\nThe answers are a bit different due to rounding off error in textbook\n" + } + ], + "prompt_number": 20 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 7.6,Page 413" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos, tan\nfrom numpy import *\nP1=200.0 #psia\nT1=500.0 #F\nP2=1.0 #psia\nalpha=20.0 #degrees\nn=3600.0\ng=32.2 #ft/s^2\n\n#calculations\n#\"From mollier charts\")\nV1=2450 #fps\nVb=V1*cos(alpha*pi/180) /2\nR=Vb*60/(n*2*pi)\nwork=1/32.2 *(V1*cos(alpha*pi/180))*Vb\nw3=3*work #workdone\n\n#results\nprint \"Blade velocity in fps\",round(Vb,3)\nprint \"Blade radius in ft\",round(R,3)\nprint \" Work done in per lb of steam\",round(w3,3)\nprint \"The answers are a bit different due to rounding off error in textbook\"\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Blade velocity in fps 1151.123\nBlade radius in ft 3.053\n Work done in per lb of steam 246910.29\nThe answers are a bit different due to rounding off error in textbook\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 7.7, Page 416" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nP1=200.0 #psia\nT1=500.0 #F\nP2=1.0 #psia\nalpha=20.0 #degrees\nn=3600.0\ng=32.2 #ft/s^2\nstage=2.0\n\n#calculations\n#\"From mollier charts\"\nV1=4240.0 #fps\nVb=V1*cos(alpha*pi/180) /(2*stage)\nR=Vb*60/(n*2*pi)\nV1x=3980 #fps\nV2x=-1990 #fps\nwork1=1/g *(V1x-V2x)*Vb #work\nwork2=1/g *(-V2x)*Vb #work\nwt=work1+work2 #work done\n\n#results\nprint \"Blade velocity in fps\",round(Vb,3)\nprint \" Blade radius in ft\",round(R,3)\nprint \" Work done in per lb of steam\",round(wt,3)\nprint \"The answers are a bit different due to rounding off error in textbook\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Blade velocity in fps 996.074\n Blade radius in ft 2.642\n Work done in per lb of steam 246234.486\nThe answers are a bit different due to rounding off error in textbook\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 7.8,Page 422" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan,cos, tan\nfrom numpy import *\nalpha=20.0 #degrees\nn=3600.0\ng=32.2 #ft/s^2\nV1=500 #fps\n\n#calculations\nVb=V1*cos(alpha*pi/180) #blade velocity\nV1x=Vb\nwork=1/32.2 *(V1x)*Vb #work done\n\n#results\nprint \"Blade velocity in fps\",round(Vb,3)\nprint \" Work done in per lb of steam\",round(work,3)\nprint \"The answers are a bit different due to rounding off error in textbook.\"\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Blade velocity in fps 469.846\n Work done in per lb of steam 6855.763\nThe answers are a bit different due to rounding off error in textbook.\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 7.9,Page 436" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\npower=1000.0 #kw\nms=16000.0 #lb/hr\nP=200.0 #psia\nT=540 #F\n\n#calculations\n#\"From mollier charts\"\nh1=1290.0 #Btu/hr\nh2=940.0 #Btu/hr\ndh=h1-h2\nrate=3413.0/dh\nact=ms/power #steam rate\n\n#results\nprint \"Ideal steam rate in per kw hr\",round(rate,3)\nprint \" Actual steam rate in per kw hr\",round(act,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Ideal steam rate in per kw hr 9.751\n Actual steam rate in per kw hr 16.0\n" + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 7.10,Page 439" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nP=200 #psia\nT=540 #F\npower=1000 #kw\nms=16000 #lb/hr\n\n#calculations\n#\"From mollier charts\")\nh1=1290.0 #Btu/hr\nh2=940.0 #Btu/hr\ndh=h1-h2\nhf2=83.0 #Btu/lb\netat=(h1-h2)/(h1-hf2)\nact=power*3413/(ms*(h1-hf2)) #thermal efficiency\netae=act/etat #efficiency of engine\n\n#results\nprint \"Ideal thermal efficiency in percent\",round(etat*100,3)\nprint \"Actual thermal efficiency in percent\",round(act*100,3)\nprint \" Engine efficiency in percent\",round( etae*100,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Ideal thermal efficiency in percent 28.998\nActual thermal efficiency in percent 17.673\n Engine efficiency in percent 60.946\n" + } + ], + "prompt_number": 7 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter8.ipynb b/Elementary_heat_Power/Chapter8.ipynb new file mode 100755 index 00000000..94c8be36 --- /dev/null +++ b/Elementary_heat_Power/Chapter8.ipynb @@ -0,0 +1,62 @@ +{ + "metadata": { + "name": "EHP-8" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Steam Engines" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 8.1,Page 461" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\narea1=2.7\nlength=3.4\nscale=60\narea2=2.75\ndia=12 #ft\nd2=2.5 #ft\nL=15.0/12 #ft\nn=250.0 #rpm\nF=600.0 #lb\nr=3 #ft\n\n#calculations\nAh=dia**2 *pi/4\nAc=(dia**2 -d2**2)*pi/4\nPih=area1/length *scale\nPic=area2/length *scale\nHihp=Pih*L*Ah*n/33000.0\nCihp=Pic*L*Ac*n/33000.0\nTihp=Hihp+Cihp\nBhp=2*pi*r*F*n/33000.0\nFhp=Tihp-Bhp\neff=Bhp/Tihp *100 #efficiency\n\n#results\nprint \"Ihp in ihp\",round(Tihp,3)\nprint \"Bhp in bhp\",round(Bhp,3)\nprint \"Fhp in fhp\",round(Fhp,3)\nprint \"Efficiency in percent\",round(eff,3)\nprint \"The answer is a bit different due to rounding off error in the textbook\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Ihp in ihp 100.749\nBhp in bhp 85.68\nFhp in fhp 15.069\nEfficiency in percent 85.043\nThe answer is a bit different due to rounding off error in the textbook\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 8.2,Page 468" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nIhp=101.1\nBhp=85.7\nmd=3000.0 #lb/hr\nh1=1172.0 #Btu/hr\nh2=180.0 #Btu/hr\nh3=1025.0 #Btu/hr\n\n#calculations\neta1=Ihp*2545/(md*(h1-h2)) *100 #efficiency\neta2=Bhp*2545/(md*(h1-h2)) *100 #efficiency\netat=(h1-h3)/(h1-h2) *100 #efficiency\nengeff=eta1/etat *100 #engine efficiency\nrate1= md/Ihp\nrate2=md/Bhp\nh22=h1-2545/rate1 #enthalpy\n\n#results\nprint \"Actual thermal efficiency based upon Ihp in lb per ihp hr\",round(eta1,3)\nprint \"Actual thermal efficiency based upon Bhp in lb per ihp hr\",round(eta2,3)\nprint \"Ideal thermal efficiency in percent \",round(etat,3)\nprint \"Engine efficiency in percent\",round(engeff,3)\nprint \"Steam rate in lb per ihp hr\",round(rate1,3)\nprint \"Steam rate in lb per bhp hr\",round(rate2,3)\nprint \"Enthalpy of exhaust steam in Btu per lb of steam\",round(h22,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Actual thermal efficiency based upon Ihp in lb per ihp hr 8.646\nActual thermal efficiency based upon Bhp in lb per ihp hr 7.329\nIdeal thermal efficiency in percent 14.819\nEngine efficiency in percent 58.345\nSteam rate in lb per ihp hr 29.674\nSteam rate in lb per bhp hr 35.006\nEnthalpy of exhaust steam in of steam 1086.234\n" + } + ], + "prompt_number": 5 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/Chapter9.ipynb b/Elementary_heat_Power/Chapter9.ipynb new file mode 100755 index 00000000..aa870a60 --- /dev/null +++ b/Elementary_heat_Power/Chapter9.ipynb @@ -0,0 +1,230 @@ +{ + "metadata": { + "name": "EHP-9" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": "Pumps" + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.1,Page 472" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan, cos\nfrom numpy import *\nh=200.0 #ft\ngam=64.0 #lb per cu ft\n\n#calculations\nP=h*gam/144.0 #pressure\n\n#results\nprint \"Pressure in psi\",round(P,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "88.889 Pressure in psi\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.2,Page 472" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nP=20.0 #psi\ngam=62.4 #lb per cu ft\n\n#calculations\nh=P*144/gam #height\n\n#results\nprint \"height in ft\",round(h,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "46.154 height in ft\n" + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.3,Page 472" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nh=3.0/12.0 #ft\ngam=62.4 #lb per cu ft\ngam2=0.075 #lb per cu ft\n\n#calculations\nP=h*gam #pressure\nh2=P/gam2 #height\n\n#results\nprint \"Air height required in ft of air\",round(h2,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Air height required in ft of air 208.0\n" + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.4,Page 477" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\ndif=4.0 #in\ngam=62.4 #lb per cu ft\ndensity=13.6 #g/cc\n\n#calculations\npv=dif*1/12 *density*gam/144 - dif/12 *gam/144 #pressure\nhv=pv*144/gam #height\n\n#results\nprint \"velocity pressure in psi\",round(pv,3)\nprint \" velocity head in ft of water \",round(hv,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "velocity pressure in psi 1.82\n velocity head in ft of water 4.2\n" + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.5,Page 479" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\ndif=4.0 #in\ngam=62.4 #lb per cu ft\ngam2=0.08 #lb per cu ft\n\n#calculations\npv=dif*1/12 *gam/144 #pressure\nhv=pv*144/gam2 #height\n\n#results\nprint \"velocity pressure in psi\",round(pv,3)\nprint \" velocity head in ft of water \",round(hv,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "velocity pressure in psi 0.144\n velocity head in ft of water 260.0\n" + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.6,Page 482" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan,sqrt\nfrom numpy import *\nhw=3.0/12 #ft\ngam1=62.4 #lb/ft^3\ngam2=0.07 #lb/ft^3\ng=32.2 #ft/s^2\n\n#calculations\np=hw*gam1 #pressure\nhg=p/gam2 #height\nV=sqrt(2*g*hg) #velocity\n\n#results\nprint \"velocity of gas in fps\",round(V,3)\n", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "velocity of gas in fps 119.8\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.7,Page 485" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nh=4.0 #in\nden=13.6 #g/cc\nAr=1.0/9\nA1=12 #sq in\ngam=62.4 #lb/ft^3\ng=32.2 #ft/s^2\n\n#calculations\ndh=(h*den-h)/12.0\nVr=1/Ar\nV22=2*g*dh/(1-Ar**2)\nV2=sqrt(V22)\nA2=A1*Ar\nv2=1/gam\nms=A2*V2/(v2*144) #mass flow rate\n\n#results\nprint \"Flow rate of water in lb/sec\",round(ms,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Flow rate of water in lb/sec 9.561\n" + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.8,Page 488" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nmdot=8000.0 #lb/min\nA1=1.0 #sq ft\nA2=3.0/4 #sq ft\nP2=50.0 #psi\nP1=10.0 #psi\ngam=62.4 #lb/ft^3\ny2=-2.0 #ft\ny1=-4.0 #ft\ng=32.2 #ft/s^2\neff=0.7\n\n#calculations\nv=1/gam\ncap=mdot/8.33\nV1=mdot*v/A1 /60.0 #velocity\nV2=mdot*v/A2 /60.0 #velocity\nht= (y2-y1) + (V2**2 -V1**2)/(2*g) + (P2-P1)*144/gam #height\nHhp=mdot*ht/33000.0 #horsepower\nPhp=Hhp/eff #horsepower\n\n#results\nprint \"Capacity in gpm\",round(cap,3)\nprint \" Total dynamic head in ft\",round(ht,3)\nprint \" Hydraulic horsepower in hp\",round(Hhp,3)\nprint \" pump horsepower in hp\",round(Php,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Capacity in gpm 960.384\n Total dynamic head in ft 94.363\n Hydraulic hp in hp 22.876\n pump horsepower in hp 32.68\n" + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.9,Page 491" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nz=12.0 #ft\ngam1=62.4 #lb/ft^3\nsg=0.8\nP2=100.0 #psia\nP1=-10.0 #psia\nmm=10000.0 #lb/min\n\n#calculations\ngam2=sg*gam1\np2g=P2*144/(gam2) +z #pressure\np1g=P1*144*0.491/(gam2) #pressure\nht=p2g-p1g #height\nHhp=mm*ht/33000 #horsepower\n\n#results\nprint \"Total dynamic head in ft of oil\",round(ht,3)\nprint \" Hydraulic horsepower in hp\",round(Hhp,3)", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "Total dynamic head in ft of oil 314.625\n Hydraulic hp in hp 95.341\n" + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": "Example 9.10,Page 510" + }, + { + "cell_type": "code", + "collapsed": false, + "input": "# Initialization of Variable\nfrom math import pi\nfrom math import atan\nfrom numpy import *\nsr=2\n\n#calculations\nhr=sr**2 #ratio\ncapr=sr #ratio\nhpr=sr**3 #ratio\n\n#results\nprint \"head is \",round(hr,3),\"times the original\"\nprint \" capacity is \",round(capr,3),\"times the original\"\nprint \" power is\",round(hpr,3),\"times the original\"", + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": "head is 4.0 times the original\n capacity is 2.0 times the original\n power is 8.0 times the original\n" + } + ], + "prompt_number": 3 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Elementary_heat_Power/README.txt b/Elementary_heat_Power/README.txt new file mode 100755 index 00000000..c837b597 --- /dev/null +++ b/Elementary_heat_Power/README.txt @@ -0,0 +1,10 @@ +Contributed By: Avik Kumar Das +Course: btech +College/Institute/Organization: Indian Institute of Technology Bombay +Department/Designation: Civil Engineering +Book Title: Elementary heat Power +Author: H L Solberg +Publisher: John Wiley and Sons, London +Year of publication: 1986 +Isbn: 0131106724 +Edition: 2
\ No newline at end of file diff --git a/Elementary_heat_Power/screenshots/chapter10.png b/Elementary_heat_Power/screenshots/chapter10.png Binary files differnew file mode 100755 index 00000000..45b6135c --- /dev/null +++ b/Elementary_heat_Power/screenshots/chapter10.png diff --git a/Elementary_heat_Power/screenshots/chapter2.png b/Elementary_heat_Power/screenshots/chapter2.png Binary files differnew file mode 100755 index 00000000..c8f462da --- /dev/null +++ b/Elementary_heat_Power/screenshots/chapter2.png diff --git a/Elementary_heat_Power/screenshots/chapter3.png b/Elementary_heat_Power/screenshots/chapter3.png Binary files differnew file mode 100755 index 00000000..0e9459d3 --- /dev/null +++ b/Elementary_heat_Power/screenshots/chapter3.png |
