diff options
Diffstat (limited to 'Elementary_Heat_Power_by_H_L_Solberg')
12 files changed, 3570 insertions, 0 deletions
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/1-Matter_and_energy.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/1-Matter_and_energy.ipynb new file mode 100644 index 0000000..7a4ff72 --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/1-Matter_and_energy.ipynb @@ -0,0 +1,388 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 1: Matter and energy" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.10: Example_10.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g=32.2 //ft/s^2\n", +"p1=100 //psia\n", +"p2=25 //psia\n", +"v2=2 //cu ft\n", +"//calculations\n", +"W=p1*144*v2 - p2*144*v2\n", +"//results\n", +"printf('Work done = %d ft-lb',W)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.11: Example_11.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g=32.2 //ft/s^2\n", +"n=100 //rpm\n", +"p1=100 //psia\n", +"p2=25 //psia\n", +"v2=2 //cu ft\n", +"//calculations\n", +"W=p1*144*v2 - p2*144*v2\n", +"Hp=W*n/33000\n", +"//results\n", +"printf('Horsepower developed = %.1f hp',Hp)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.12: Example_12.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P=50 //hp\n", +"m=30 //lb\n", +"E=19000 //Btu/lb\n", +"//calculations\n", +"eta= P*2545/(m*E) *100\n", +"//results\n", +"printf('Efficiency = %.1f percent',eta)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g=32.2 //ft/s^2\n", +"m=500 //lb\n", +"rate=10 //ft/s^2\n", +"//calculations\n", +"F1=m/g *rate\n", +"ms=m/g\n", +"F2=ms*rate\n", +"//results\n", +"printf('Force in case 1 = %.1f lbf',F1)\n", +"printf('\n Force in case 2 = %.1f lbf',F2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g=32.2 //ft/s^2\n", +"g2=32.0 //ft/s^2\n", +"rate=10 //ft/s^2\n", +"w1=500 //lbf\n", +"//calculations\n", +"fd1=w1*g2/g\n", +"F=fd1/g2 *rate\n", +"ms=w1/g\n", +"F2=ms*rate\n", +"//results\n", +"printf('Net weight of body in case 1 = %.1f lbf',F)\n", +"printf('\n Force in case 2 = %.1f lbf',F2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.3: Example_3.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g=32.174 //ft/s^2\n", +"m=500 //lbm\n", +"rate=10 //ft/s^2\n", +"//calculations\n", +"F=1/g *m*rate\n", +"//results\n", +"printf('Force required = %.1f lbf',F)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.4: Example_4.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g1=32.174 //ft/s^2\n", +"gc=g1\n", +"g2=30 //ft/s^2\n", +"m=100 //lbm\n", +"//calculations\n", +"w1=g1/gc *m\n", +"w2=g2/gc *m\n", +"//results\n", +"printf('Weight in case 1 = %d lbf',w1)\n", +"printf('\n Weight in case 2 = %.1f lbf',w2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.5: Example_5.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"ge=32.174 //ft/s^2\n", +"gm=5.47 //ft/s^2\n", +"we=50 //lbm\n", +"//calculations\n", +"wm=we*gm/ge\n", +"//results\n", +"printf('In case a, it will weigh the same, weight = %d lbm',we)\n", +"printf('\n In case b, weight = %.1f lbf',wm)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.6: Example_6.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g=32.2 //ft/s^2\n", +"p1=100 //psig\n", +"p2=29.0 //in of Hg\n", +"//calculations\n", +"BP=p2*0.491\n", +"AP=BP+p1\n", +"//results\n", +"printf('Absolute pressure = %.2f psia',AP)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.7: Example_7.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g=32.2 //ft/s^2\n", +"Pb=29.5 //in of Hg\n", +"Pv=10 //in of Hg\n", +"//calculations\n", +"AP=(Pb-Pv)*0.491\n", +"//results\n", +"printf('Absoulte pressure = %.2f psia',AP)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.8: Example_8.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g=32.2 //ft/s^2\n", +"v1=1 //cu ft\n", +"p1=100 //psia\n", +"//calculations\n", +"v2=2*v1\n", +"W=144*p1*(v2-v1)\n", +"//results\n", +"printf('Work done = %d ft-lb',W)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.9: Example_9.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"g=32.2 //ft/s^2\n", +"v1=1 //cu ft\n", +"p1= 100 //psia\n", +"p2=50 //psia\n", +"v2=3 //cu ft\n", +"//calculations\n", +"pa=(p1+p2)/2\n", +"W=pa*(v2-v1)*144\n", +"//results\n", +"printf('Work done = %d ft-lb',W)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/10-Drafts_fans_blowers_and_compressors.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/10-Drafts_fans_blowers_and_compressors.ipynb new file mode 100644 index 0000000..14df2c4 --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/10-Drafts_fans_blowers_and_compressors.ipynb @@ -0,0 +1,303 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 10: Drafts fans blowers and compressors" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"hb=29 //in of Hg\n", +"sg=0.491\n", +"Ra=53.3\n", +"Ta=460+40 //R\n", +"Tg=540+460 //R\n", +"H=300 //ft\n", +"gam=62.4 //lb/cu ft\n", +"//calculations\n", +"pb=hb*sg*144\n", +"rhoa=pb/(Ra*Ta)\n", +"rhog=pb/(Ra*Tg)\n", +"dp=H*(rhoa-rhog)\n", +"D=dp/(gam)\n", +"//results\n", +"printf('Theoretical draft = %.1f psf',dp)\n", +"printf('\n Draft = %.2f ft H2O',D)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"md=15 //lb per lb of coal\n", +"x=0.1\n", +"mss=1 //basis\n", +"rea=29 //in of Hg\n", +"sg=0.491\n", +"R=53.3\n", +"T=540+460 //R\n", +"V=25 //fps\n", +"gam=0.038 //lb/ft^3\n", +"//calculations\n", +"m=mss-mss*x+md\n", +"ms=m\n", +"rhog=rea*0.491*144/(R*T)\n", +"A=ms/(gam*V)\n", +"//results\n", +"printf('stack area = %.1f sq ft',A)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10.3: Example_3.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"p=144*29*0.491 //psf\n", +"R=53.3\n", +"T=70+460 //R\n", +"gamw=62.4 //lb/ft^3\n", +"gama=0.073 //lb/ft^3\n", +"hw=3/12 //ft\n", +"hw2=3.5/12 //ft\n", +"hv=32.2 //ft/s^2\n", +"ms=9 //lb\n", +"g=32.2 //ft/s^2\n", +"//calculations\n", +"rhoa=p/(R*T)\n", +"hs=hw*gamw/gama\n", +"ht=hw2*gamw/gama\n", +"hv=ht-hs\n", +"V=sqrt(2*g*hv)\n", +"msv=ms*V*60\n", +"mm=msv*gama\n", +"airhp= ht*mm/33000\n", +"//results\n", +"printf('Velocity head = %d ft of air',hv)\n", +"printf('\n velocity of air in the duct = %.1f fps',V)\n", +"printf('\n volume = %d cu ft per min',msv)\n", +"printf('\n Mass flow rate = %d lb/min',mm)\n", +"printf('\n Air hp = %.1f hp',airhp)\n", +"disp('The answers in the textbook are a bit different due to rounding off error in the textbook. Please use a calculator')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10.4: Example_4.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"A2=9 //sq ft\n", +"p2=3/12 *62.4 //psf\n", +"p1=-1/12 *62.4 //psf\n", +"ms=20000 //cfm\n", +"A1=16 //sq ft\n", +"gam=0.075 //lb/ft^3\n", +"g=32.2 //ft/s^2\n", +"inp=17 //hp\n", +"//calculations\n", +"V2=ms/60 *1/A2\n", +"V1=ms/60 *1/A1\n", +"ht=(p2-p1)/gam +(V2^2 -V1^2)/(2*g)\n", +"airhp=ht*ms*gam/33000\n", +"eta=airhp/inp *100\n", +"//results\n", +"printf('Total head = %.1f ft of air',ht)\n", +"printf('\n Air hp = %.1f hp',airhp)\n", +"printf('\n Effifciency = %.1f percent',eta)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10.5: Example_5.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"n1=400 //rpm\n", +"mv1=10000 //lb\n", +"mv2=15000 //lb\n", +"h1=2 //in of water\n", +"hp1=4 //hp\n", +"//calculations\n", +"n2=mv2/mv1 *n1\n", +"h2=h1*(n2/n1)^2\n", +"hp2=hp1 *(n2/n1)^3\n", +"//results\n", +"printf('The speed = %d rpm',n2)\n", +"printf('\n The pressure = %.1f in of water',h2)\n", +"printf('\n Power = %.1f hp',hp2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10.6: Example_6.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"m=100000 //lb/hr\n", +"p1=1 //psia\n", +"x=0.8\n", +"p2=14.7 //psia\n", +"t2=300 //F\n", +"//calculations\n", +"disp('from table A3 and A2')\n", +"h2=1192.8 //Btu/lb\n", +"hf=69.7 //Btu/lb\n", +"hfg=1036.3 //Btu/lb\n", +"h1=hf+x*hfg\n", +"W=h2-h1\n", +"power=m*W\n", +"hp=power/2545\n", +"//results\n", +"printf('Power required = %d hp',hp)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10.7: Example_7.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"p1=14.7 //psia\n", +"t1=60 //F\n", +"p2=60 //psia\n", +"t2=440 //F\n", +"m=10 //lb/sec\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h2=216.3 //Btu/lb\n", +"h1=124.3 //Btu/lb\n", +"W21=h2-h1\n", +"power=W21*m\n", +"hp=power*3600/2545\n", +"cp=0.237\n", +"W212=cp*(t2-t1)\n", +"power2=W212*m\n", +"hp2=power2*3600/2545\n", +"//results\n", +"printf('Power required = %d hp',hp)\n", +"printf('\n Power required = %d hp',hp2)\n", +"printf('\n Work done = %.1f Btu/lb',W212)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/11-Feed_water_heaters_and_condensers.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/11-Feed_water_heaters_and_condensers.ipynb new file mode 100644 index 0000000..d19e3c7 --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/11-Feed_water_heaters_and_condensers.ipynb @@ -0,0 +1,162 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 11: Feed water heaters and condensers" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"m1=1000 //lb/hr\n", +"m2=5000 //lb/hr\n", +"m3=3000 //lb/hr\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h5=196.16 //Btu/lb\n", +"h1=38.04 //Btu/lb\n", +"h2=67.97 //Btu/lb\n", +"h3=117.89 //Btu/lb\n", +"h4=1156.3 //Btu/lb\n", +"m4=(m1*h1+m2*h2+m3*h3-(m1+m2+m3)*h5)/(h5-h4)\n", +"//results\n", +"printf('Pounds of steam entering the heater = %d lb/hr',m4)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P1=100 //psia\n", +"T1=400 //F\n", +"T2=70 //F\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h1=1227.6 //Btu/lb\n", +"h2=298.4 //Btu/lb\n", +"h3=279.9 //Btu/lb\n", +"h4=38.04 //Btu/lb\n", +"m1=(h3-h4)/(h1-h2)\n", +"//results\n", +"printf('Mass of steam required = %.2f lb steam per lb water',m1)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.3: Example_3.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"h0=1260 //Btu/lb\n", +"msr=15 //lb\n", +"m4=15 //lb per hr per kw\n", +"t2=80 //F\n", +"t3=60 //F\n", +"//calculations\n", +"h1=h0-3413/msr\n", +"disp('from mollier charts,')\n", +"h4=58 //Btu/lb\n", +"dt=t2-t3\n", +"m3=m4*(h1-h4)/dt\n", +"//results\n", +"printf('enthalpy of steam entering the condenser = %d Btu/lb',h1)\n", +"printf('\n mass of cooling water = %d lb per hr per kw',m3)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.4: Example_4.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"m4=8*1000000 //lb per hr\n", +"dt=12 //F\n", +"//calculations\n", +"disp('from mollier charts,')\n", +"dh4=950 //Btu/lb\n", +"m3=m4*(dh4)/dt\n", +"//results\n", +"printf('\n mass of cooling water = %.3e lb per hr',m3)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/12-The_Gas_turbine_power_plant.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/12-The_Gas_turbine_power_plant.ipynb new file mode 100644 index 0000000..54491f3 --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/12-The_Gas_turbine_power_plant.ipynb @@ -0,0 +1,113 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 12: The Gas turbine power plant" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"T1=80 //F\n", +"T2=460 //F\n", +"T3=1300 //F\n", +"T4=780 //F\n", +"//calculations\n", +"disp('from mollier charts,')\n", +"h1=129.1 //Btu/lb\n", +"h2 = 221.2 //Btu/lb\n", +"h3= 438.8 //Btu/lb\n", +"h4 = 301.5 //Btu/lb\n", +"wcom=h2-h1\n", +"wcob=h3-h2\n", +"wtur=h3-h4\n", +"eta=(wtur-wcom)/wcob *100\n", +"//results\n", +"printf('\n work done by compressor = %.1f btu input as work per lb of air compressed',wcom)\n", +"printf('\n Heat supplied in the combustor = %.1f Btu supplied per lb of air ',wcob)\n", +"printf('\n work done in the turbine = %.1f Btu output as work per lb of air',wtur)\n", +"printf('\n Cycle efficiency = %.1f percent',eta)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"T1=80 //F\n", +"T2=460 //F\n", +"T=700 //F\n", +"T3=1300 //F\n", +"T4=780 //F\n", +"//calculations\n", +"disp('from mollier charts,')\n", +"h1=129.1 //Btu/lb\n", +"h2 = 221.2 //Btu/lb\n", +"h3= 438.8 //Btu/lb\n", +"h4 = 301.5 //Btu/lb\n", +"wcom=h2-h1\n", +"wcob=h3-h2\n", +"wtur=h3-h4\n", +"output=-wcom+wtur\n", +"h=281.1 //Btu/lb\n", +"Q=h3-h\n", +"eff=output/Q *100\n", +"//results\n", +"printf('\n Heat supplied in the combustor = %.1f Btu supplied per lb of air ',Q)\n", +"printf('\n Cycle efficiency = %.1f percent',eff)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/13-Mechanical_Refrigeration.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/13-Mechanical_Refrigeration.ipynb new file mode 100644 index 0000000..17f0ba3 --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/13-Mechanical_Refrigeration.ipynb @@ -0,0 +1,134 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 13: Mechanical Refrigeration" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"\n", +"disp('From mollier diagram from ammonia, values are found')\n", +"disp('part a')\n", +"h1=65 //Btu/lb\n", +"printf('enthalpy in case a = %d Btu/lb',h1)\n", +"h2=99 //Btu/lb\n", +"v2=0.93 //ft^3/lb\n", +"printf('\n In case 2, enthalpy and specific volume are %d Btu/lb and %.2f ft^3/lb respectively',h2,v2)\n", +"h3=583 //Btu/lb\n", +"v3=8.8 //ft^3/lb\n", +"s3=1.275\n", +"printf('\n In case 3, enthalpy, specific volume and entropy are %d Btu/lb, %.2f ft^3/lb and %.3f respectively',h3,v3,s3)\n", +"h4=720 //Btu/lb\n", +"v4=10.4 //ft^3/lb\n", +"s4=1.50\n", +"printf('\n In case 4, enthalpy, specific volume and entropy are %d Btu/lb, %.2f ft^3/lb and %.3f respectively',h4,v4,s4)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"mr=3 //lb\n", +"mj=5 //lb\n", +"t2=67 //F\n", +"t1=60 //lb\n", +"ihp=7.25\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h4=709 //Btu/b\n", +"h3=618 //Btu/lb\n", +"energyin=ihp*2545/60\n", +"energyout=mr*(h4-h3) + mj*(t2-t1)\n", +"//results\n", +"printf('Energy in = %.1f Btu/min',energyin)\n", +"printf('\n Energy out = %.1f Btu/min',energyout)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.3: Example_3.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"mr=3 //lb\n", +"hp=10 //hp\n", +"//calculations\n", +"h3=618 //Btu/lb\n", +"h1=131 //Btu/lb\n", +"Qe=mr*(h3-h1)\n", +"work=hp*2545/60\n", +"cop=Qe/work\n", +"//results\n", +"printf('Coefficient of performance = %.2f',cop)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/2-Fuels_and_Combustion.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/2-Fuels_and_Combustion.ipynb new file mode 100644 index 0000000..59b6918 --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/2-Fuels_and_Combustion.ipynb @@ -0,0 +1,547 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 2: Fuels and Combustion" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.10: Example_10.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"Gf=11.57 //lb per lb of fuel\n", +"H=4.4/100\n", +"M=13.5/100\n", +"mr=700\n", +"mf=10000\n", +"mc=1 //lb\n", +"//calculations\n", +"pro=M+9*H\n", +"mrf=mr/mf\n", +"Aa=Gf+pro+mrf-mc\n", +"At=8.83\n", +"ea=(Aa-At)/At *100\n", +"//results\n", +"printf('Excess air = %.1f percent',ea)\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.11_a: Example_11.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear\n", +"//Initialization of variables\n", +"Gf=11.57 //lb per lb of fuel\n", +"tg=500 //F\n", +"ta=70 //F\n", +"//calculations\n", +"Q1=0.24*Gf*(tg-ta)\n", +"//results\n", +"printf('Heat loss = %d Btu per lb of fuel',Q1)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.11_b: Example_12.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"Co=0.1\n", +"Co2=14.1\n", +"Cb=0.646\n", +"//calculations\n", +"Q2=Co/(Co+Co2) *Cb*10160\n", +"//results\n", +"printf('Heat loss = %d Btu per lb of fuel',Q2)\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.11c: Example_13.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"mf=10000 //lb \n", +"mr=700 //lb\n", +"Cr=0.2\n", +"//calculations\n", +"Q3=mr*Cr/mf *14600\n", +"//results\n", +"printf('Heat loss = %d Btu per lb of fuel',Q3)\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.11d: Example_14.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"M=0.135\n", +"tg=500 //F\n", +"ta=70 //F\n", +"//calculations\n", +"Q4=M*(1089+0.46*tg-ta)\n", +"//results\n", +"printf('Heat loss = %.1f Btu per lb of fuel',Q4)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.11e: Example_15.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"Per=0.044 //percentage\n", +"tg=500 //F\n", +"ta=70 //F\n", +"//calculations\n", +"Q5=9*Per*(1089+0.46*tg-ta)\n", +"//results\n", +"printf('Heat loss = %.1f Btu per lb of fuel',Q5)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"x1=0.135\n", +"x2=0.056\n", +"veca=[32.5 48.4 5.6 13.5]\n", +"B1=11788\n", +"//calculations\n", +"vecb=veca/(1-x1)\n", +"vecc=veca/(1-x1-x2)\n", +"B2=B1/(1-x1)\n", +"B3=B1/(1-x1-x2)\n", +"vecb(4)=0\n", +"vecc(3)=0\n", +"vecc(4)=0\n", +"//results\n", +"printf('In Moisture free case, ')\n", +"format('v',6);vecb\n", +"disp(vecb)\n", +"printf('In Moisture and Ash free case, ')\n", +"format('v',6);vecc\n", +"disp(vecc)\n", +"printf('Energy in Moisture free case = %d Btu per lb',B2)\n", +"printf('\n Energy in Moisture and ash free case = %d Btu per lb',B3)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"y1=13.5\n", +"x1=0.135\n", +"x2=0.056\n", +"veca=[66 1.5 1.1 5.6 5.9 19.9]\n", +"//calculations\n", +"vecb=[veca y1]\n", +"vecb(5) = vecb(5) - 1/9*y1\n", +"vecb(6) = vecb(6) - 8/9*y1\n", +"vecc=vecb/(1-x1)\n", +"vecd=vecb/(1-x1-x2)\n", +"vecd(4)=0\n", +"vecd(7)=0\n", +"vecc(7)=0\n", +"s1=sum(vecc)\n", +"s2=sum(vecd)\n", +"//results\n", +"printf('With moisture as a separate item, ')\n", +"format ('v',6);vecb\n", +"disp(vecb)\n", +"printf('In Moisture free case, ')\n", +"format('v',4);vecc\n", +"disp(vecc)\n", +"printf('In Moisture and Ash free case, ')\n", +"format('v',5);vecd\n", +"disp(vecd)\n", +"printf('Total Mositure free content = %.1f percent',s1)\n", +"printf('\n Total Mositure and ash free content = %.1f percent',s2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.3: Example_3.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"H=5.9\n", +"O=19.9\n", +"H2=4.4\n", +"O2=7.9\n", +"//calculations\n", +"Ha1=H-O/8\n", +"Ha2=H2-O2/8\n", +"//results\n", +"printf('Available hydrogen in case 1 = %.1f percent by weight',Ha1)\n", +"printf('\n Available hydrogen in case 1 = %.1f percent by weight',Ha2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.4: Example_4.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"H1=0.059\n", +"O1=0.199\n", +"H2=0.044\n", +"O2=0.079\n", +"C=0.66\n", +"S=0.011\n", +"//calculations\n", +"Qh1= 14600*C+62000*(H1-O1/8)+4050*S\n", +"Qh2=14600*C+62000*(H2-O2/8)+4050*S\n", +"//results\n", +"printf('Heating value in case 1 = %d Btu/lb ',Qh1)\n", +"printf('\n Heating value in case 2 = %d Btu/lb ',Qh2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.5: Example_5.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"H1=0.059\n", +"O1=0.199\n", +"C=0.66\n", +"S=0.011\n", +"//calculations\n", +"Qh1= 11.52*C+34.56*(H1-O1/8)+4.32*S\n", +"//results\n", +"printf('Theoretical air required = %.2f lb of air per lb of coal ',Qh1)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.6: Example_6.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"mf=10000 //lb\n", +"mr=700 //lb\n", +"Cr=0.20\n", +"Cco2=14.1\n", +"Co2=5.1\n", +"Cco=0.1\n", +"Cf=0.66\n", +"//calculations\n", +"Cn2=100-(Cco2+Co2+Cco)\n", +"Ci=mf*Cf\n", +"Ca=mr*Cr\n", +"Cb=(Ci-Ca)/mf\n", +"Cb2=((mf*Cf)-mr*Cr)/(mf)\n", +"veca=[Cco2 Co2 Cco Cn2]\n", +"vecb=veca\n", +"vecb(1)=vecb(1) *44\n", +"vecb(2)=vecb(2) *32\n", +"vecb(3)=vecb(3) *28\n", +"vecb(4)=vecb(4) *28\n", +"sumvec=sum(vecb)\n", +"Lbc=Cco2*12 + Cco*12\n", +"Gc=sumvec/Lbc\n", +"Gf=Gc*Cb\n", +"//results\n", +"printf('Carbon in the dry products combustion = %.3f lb per lb of fuel',Cb)\n", +"printf('\n In case 2, Carbon in the dry products combustion = %.3f lb per lb of fuel',Cb2)\n", +"printf('\n Dry gaseous products of combstion per lb of coal = %.2f lb ',Gf)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.7: Example_7.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"mf=10000 //lb\n", +"mr=700 //lb\n", +"Cr=0.20\n", +"Cco2=14.1\n", +"Co2=5.1\n", +"Cco=0.1\n", +"Cf=0.66\n", +"//calculations\n", +"Cn2=100-(Cco2+Co2+Cco)\n", +"Ci=mf*Cf\n", +"Ca=mr*Cr\n", +"Cb=(Ci-Ca)/mf\n", +"Cb2=((mf*Cf)-mr*Cr)/(mf)\n", +"veca=[Cco2 Co2 Cco Cn2]\n", +"vecb=veca\n", +"vecb(1)=vecb(1) *44\n", +"vecb(2)=vecb(2) *32\n", +"vecb(3)=vecb(3) *28\n", +"vecb(4)=vecb(4) *28\n", +"Cbb1=Cb*Cco*12/(Cco2*12 + Cco*12)\n", +"Cbb2= Cb*(veca(3) /(veca(3) + veca(1)))\n", +"//results\n", +"printf('In case 1, Carbon burned per lb of fuel = %.5f lb per lb of fuel',Cbb1)\n", +"printf('\n In case 2, Carbon burned per lb of fuel = %.5f lb per lb of fuel',Cbb2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.8: Example_8.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"H=4.4/100\n", +"M=13.5/100\n", +"H2=0.059\n", +"//calculations\n", +"pro=M+9*H\n", +"pro2=9*H2\n", +"//results\n", +"printf('In case 1, watervapor present in products = %.3f lb',pro)\n", +"printf('\n In case 2, watervapor present in products = %.3f lb',pro2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.9: Example_9.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"Gf=11.57 //lb per lb of fuel\n", +"H=4.4/100\n", +"M=13.5/100\n", +"mr=700\n", +"mf=10000\n", +"mc=1 //lb\n", +"//calculations\n", +"pro=M+9*H\n", +"mrf=mr/mf\n", +"Aa=Gf+pro+mrf-mc\n", +"//results\n", +"printf('Actual air supplied = %.2f lb of air supplied per lb of fuel',Aa)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/3-Internal_Combustion_Engines.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/3-Internal_Combustion_Engines.ipynb new file mode 100644 index 0000000..f982fbe --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/3-Internal_Combustion_Engines.ipynb @@ -0,0 +1,420 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 3: Internal Combustion Engines" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.10: Example_10.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"bore=3 //in\n", +"str=4 //in\n", +"rpm=3000 //rpm\n", +"air=110 //cu ft per min\n", +"//calculations\n", +"pdv=bore*bore*%pi*str*2*bore/4\n", +"pde=pdv*rpm /2\n", +"req=air*1728\n", +"eff=req/pde *100\n", +"//results\n", +"printf('Volumetric efficiency = %.1f percent',eff)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.11: Example_11.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"x1=11.5\n", +"x2=4.1\n", +"x3=0.4\n", +"x4=2.3\n", +"x5=0.2\n", +"x6=81.5\n", +"yc=0.842\n", +"yh=0.158\n", +"basis=1\n", +"bhp=100\n", +"burn=8.9 //gal/hr\n", +"sg=0.731\n", +"Qh=20750 //Btu/lbm\n", +"rate=66 //gpm\n", +"ex=1100 //F\n", +"air=70 //F\n", +"cp=0.254\n", +"h2=4330\n", +"h4=62000\n", +"h5=23700\n", +"//calculations\n", +"c1=x1*44\n", +"c2=x2*28\n", +"c3=x3*32\n", +"c4=x4*2\n", +"c5=x5*16\n", +"c6=x6*28\n", +"summ=c1+c2+c3+c4+c5+c6\n", +"carbon=x1*12 + x2*12+x5*12\n", +"hydrogen=x4*2+x5*4\n", +"lbdrygas=summ/carbon *yc\n", +"lbfuel=carbon/yc\n", +"lbH=lbfuel*yh\n", +"lbH2=lbH-hydrogen\n", +"lb3=lbH2*9\n", +"lbwater=lb3/lbfuel\n", +"lbair=lbdrygas+lbwater-basis\n", +"bsfc=burn*sg*8.33/bhp\n", +"fuelmin=bsfc*bhp/60\n", +"energy=2545/bsfc\n", +"per1=energy/Qh\n", +"Ec=rate*8.33*10\n", +"Eclb=Ec/fuelmin\n", +"per2=Eclb/Qh\n", +"dryloss=(ex-air)*cp*lbdrygas\n", +"per3=dryloss/Qh\n", +"hv2=h2*c2/lbfuel\n", +"hv4=h4*c4/lbfuel\n", +"hv5=h5*c5/lbfuel\n", +"hv=hv2+hv4+hv5\n", +"per4=hv/Qh\n", +"eh2=lbwater*(1066+0.5*ex-air)\n", +"per5=eh2/Qh\n", +"rad=1017 \n", +"per6=rad/Qh\n", +"//results\n", +"printf('Air supplied per lb of fuel = %.1f lb air per lb fuel',lbair)\n", +"printf('\n Percentage of energy supplied utilized in Btu = %.2f percent',per1*100)\n", +"printf('\n Percentage of energy absorbed by coolant = %.2f percent',per2*100)\n", +"printf('\n Energy lost in sensible heat = %.2f percent',per3*100)\n", +"printf('\n Energy supplied in combustiles in exhaust = %.2f percent',per4*100)\n", +"printf('\n Energy supplied in water formed by combustion = %.2f percent',per5*100)\n", +"printf('\n Energy supplied unaccounted for = %.2f percent',per6*100)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"re=6\n", +"k=1.4\n", +"//calculations\n", +"nt=1-1/re^(k-1)\n", +"ntt=nt*100\n", +"//results\n", +"printf('Thermal efficiency = %.1f percent',ntt)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"d=3.25 //in\n", +"str=4 //in\n", +"v=6 //cu in\n", +"//calculations\n", +"Dp=d^2 *%pi*str/4\n", +"r=(Dp+v)/v\n", +"//results\n", +"printf('Compression ratio = %.2f',r)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.3: Example_3.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"per=20\n", +"Dp=100\n", +"//calculations\n", +"r=Dp/per +1\n", +"//results\n", +"printf('Compression ratio = %d ',r)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.4: Example_4.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"r=16\n", +"rc=4\n", +"k=1.4\n", +"//calculations\n", +"etat=1-1/r^(k-1) *((rc^k -1)/(k*(rc-1)))\n", +"eta=etat*100\n", +"//results\n", +"printf('Thermal efficiency = %.1f percent',eta)\n", +"disp('The answer is a bit different due to rounding off error in the textbook')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.5: Example_5.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"F=200 //lb\n", +"area=1.65 //sq. in\n", +"len=3.5 //in\n", +"//calculations\n", +"ord=area/len\n", +"mep=ord*F\n", +"//results\n", +"printf('MEP of an engine = %.1f psi',mep)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.6: Example_6.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"Pi=90 //psi\n", +"L=5/12 //ft\n", +"r=5 //in\n", +"x=1.5 //ft\n", +"rpm=1500 //rpm\n", +"//calculations\n", +"A=%pi*x*x\n", +"N=rpm*4/2\n", +"Ihp=Pi*L*A*N/33000\n", +"//results\n", +"printf('IHP of cylinder = %d',Ihp)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.7: Example_7.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"r=4 //ft\n", +"n=300 //rpm\n", +"F=60 //lb\n", +"//calculations\n", +"Bhp=2*%pi*r*F*n/33000\n", +"//results\n", +"printf('Bhp of the engine = %.1f',Bhp)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.8: Example_8.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"C=1/4000\n", +"F=125 //lb\n", +"n=3500 //rpm\n", +"//calculations\n", +"Bhp=F*n*C\n", +"//results\n", +"printf('Bhp developed by the engine = %.1f',Bhp)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.9: Example_9.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"r=1.75 //ft\n", +"F1=72 //lb\n", +"F2=24 //lb\n", +"n=500 //rpm\n", +"m=1.8 //lb\n", +"mi=15 //min\n", +"Qh=20000 //Btu/lb\n", +"//calculations\n", +"Bhp=2*%pi*r*F1*n/33000\n", +"Fhp=2*%pi*r*F2*n/33000\n", +"Ihp=Bhp+Fhp\n", +"Fc=m*60/mi\n", +"Bsfc=Fc/Bhp\n", +"Isfc=Fc/Ihp\n", +"etam=Bhp/Ihp *100\n", +"etabt=Bhp*2545/(Fc*Qh)\n", +"etait=Ihp*2545/(Fc*Qh)\n", +"//results\n", +"printf('Thermal efficiency = %d percent',etam)\n", +"printf('\n Brake thermal effficiency = %.1f percent',etabt*100)\n", +"printf('\n Indicated thermal effficiency = %.1f percent',etait*100)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/5-Steam_Generation.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/5-Steam_Generation.ipynb new file mode 100644 index 0000000..0c02a5b --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/5-Steam_Generation.ipynb @@ -0,0 +1,496 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 5: Steam Generation" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.10: Example_10.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"m=150000 //lb\n", +"P1=1000 //psia\n", +"Ts=900 //F\n", +"Tf=200 //F\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h2=1448.2 //Btu/lb\n", +"hf=167.99 //Btu/lb\n", +"correc=2.2 //Btu/lb\n", +"hc=hf+correc\n", +"Q=m*(h2-hc)\n", +"//results\n", +"printf('Heat absorption = %d Btu/hr',Q)\n", +"disp('The answer is a bit different due to rounding off error in textbook')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.11: Example_11.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"m=150000 //lb\n", +"P1=1000 //psia\n", +"Ts=900 //F\n", +"Tf=200 //F\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h2=1448.2 //Btu/lb\n", +"hf=167.99 //Btu/lb\n", +"correc=2.2 //Btu/lb\n", +"hc=hf+correc\n", +"Q=m*(h2-hc)\n", +"output=Q/1000\n", +"//results\n", +"printf('Output of the steam generating unit = %d kB/hr',output)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.12: Example_12.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"m=150000 //lb\n", +"P1=1000 //psia\n", +"Ts=900 //F\n", +"Tf=200 //F\n", +"m2=21000 //lb\n", +"HV=12000 //Btu/lb\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h2=1448.2 //Btu/lb\n", +"hf=167.99 //Btu/lb\n", +"correc=2.2 //Btu/lb\n", +"hc=hf+correc\n", +"Q=m*(h2-hc)\n", +"output=Q\n", +"inpu=m2*HV\n", +"eta=output/inpu\n", +"//results\n", +"printf('Efficiency of the steam generating unit = %.1f percent',eta*100)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.13: Example_13.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"hv=11780 //Btu/lb\n", +"steam=55000 //lb/hr\n", +"coal=6480 //lb\n", +"x1=0.66\n", +"x2=0.044\n", +"x3=0.079\n", +"x4=0.015\n", +"x5=0.11\n", +"z1=14.5\n", +"z2=0.2\n", +"z3=4.4\n", +"z4=80.9\n", +"xash=0.076\n", +"xmois=0.115\n", +"yc=0.21\n", +"refuse=622 //lb/hr\n", +"cp=0.24\n", +"tg=400 //F\n", +"ta=70 //F\n", +"Qco=10160 //Btu/lb\n", +"Qc=14600 //Btu/lb\n", +"//calculations\n", +"disp('From steam tables,')\n", +"hf=269.6 //Btu/lbm\n", +"hfg=1.5 //Btu/lbm\n", +"h1=hf+hfg\n", +"h2=1196.5\n", +"Qb=h2-h1\n", +"h3=1407.7 //Btu/lbm\n", +"Qs=h3-h2\n", +"h4=h3-h1\n", +"out=steam*h4/1000\n", +"eff=steam*h4/(coal*hv)\n", +"//Energy balance\n", +"Ci=coal*x1\n", +"Cr=refuse*yc\n", +"Cb=(Ci-Cr)/coal\n", +"lbt= z1*44+z2*28+z3*32+z4*28\n", +"lbC=z1*12+z2*12\n", +"dry=lbt/lbC *Cb\n", +"loss1=dry*cp*(tg-ta)\n", +"loss2=z2*12/(lbC) *Cb*Qco\n", +"loss3=Cr*Qc/coal\n", +"loss4=xmois*(1089+0.46*tg-ta)\n", +"loss5=x2*9*(1089+0.46*tg-ta)\n", +"loss6=steam*h4/coal\n", +"//results\n", +"printf('Heat absorbed in the boiler = %.2f Btu per lb of steam generated',Qb)\n", +"printf('\n Heat absorbed in the superheater = %.2f Btu/lb of steam',Qs)\n", +"printf('\n Heat absorbed in steam generating = %.2f Btu/lb of steam generated',h4)\n", +"printf('\n Output of steam generating unit = %d kB',out)\n", +"printf('\n Efficiency of steam generating unit = %.1f percent',eff*100)\n", +"printf('\n Carbon burned to CO and CO2 = %.2f lb of C per lb of fuel',Cb)\n", +"printf('\n Dry products of combustion = %.2f lb per lb of fuel',dry)\n", +"printf('\n Loss due to sensible heat in dry gaseous products of combustion = %d Btu/lb of fuel',loss1)\n", +"printf('\n Loss due to CO in dry products of combustion = %.1f Btu/lb of fuel',loss2)\n", +"printf('\n Loss due to C in refuse = %.1f Btu/lb of fuel',loss3)\n", +"printf('\n Loss due to evaporating moisture in fuel = %.1f Btu/lb of fuel',loss4)\n", +"printf('\n Loss due to water vapor formed from H = %.1f Btu/lb of fuel',loss5)\n", +"printf('\n Energy absorbed in generating steam = %d Btu/lb of fuel',loss6)\n", +"disp('The answers are a bit different due to rounding off error in the textbook')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"x=0.98\n", +"vg=26.80\n", +"vf=0.01672\n", +"//calculations\n", +"vx=x*vg+(1-x)*vf\n", +"//results\n", +"printf('Specific volume of wet steam = %.6f cu ft per lb',vx)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"hf=167.99 //Btu/lb\n", +"hg=4.5 //Btu/lb\n", +"//calculations\n", +"hc=hf+hg\n", +"//results\n", +"printf('Enthalpy of water = %.1f Btu/lb',hc)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.3: Example_3.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"x=0.97\n", +"hg=1187.2 //Btu/lb\n", +"hf=298.40 //Btu/lb\n", +"hfg=888.8 //Btu/lb\n", +"//calculations\n", +"hx1=x*hg+(1-x)*hf\n", +"hx2=hf+x*hfg\n", +"hx3=hg-(1-x)*hfg\n", +"//results\n", +"printf('\n In case 1, enthalpy = %.1f Btu/lb',hx1)\n", +"printf('\n In case 2, enthalpy = %.1f Btu/lb',hx2)\n", +"printf('\n In case 3, enthalpy = %.1f Btu/lb',hx3)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.4: Example_4.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"h1=1172 //Btu/lb\n", +"hf1=355.36 //Btu/lb\n", +"hfg1=843 //Btu/lb\n", +"//calculations\n", +"h2=h1\n", +"x1= (h2-hf1)/hfg1\n", +"//results\n", +"printf('Quality of steam = %.1f percent',x1*100)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.5: Example_5.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P=200 //psia\n", +"x=0.95\n", +"m=1//lb\n", +"//calculations\n", +"disp('From mollier chart,')\n", +"hx=1156 //Btu/lb\n", +"sx=1.495 //Btu/lb F\n", +"//results\n", +"printf('Enthalpy = %d Btu/lb',hx)\n", +"printf('\n entropy = %.3f Btu/lb F',sx)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.6: Example_6.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P=200 //psia\n", +"T=600 //F\n", +"m=1 //lb\n", +"//calculations\n", +"disp('From mollier chart,')\n", +"hx=1322 //Btu/lb\n", +"sx=1.676 //Btu/lb F\n", +"//results\n", +"printf('Enthalpy = %d Btu/lb',hx)\n", +"printf('\n entropy = %.3f Btu/lb F',sx)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.7: Example_7.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P=200 //psia\n", +"T=260 //F\n", +"//calculations\n", +"disp('From mollier chart,')\n", +"hx=1174 //Btu/lb\n", +"x1=2.8\n", +"y1=100-x1\n", +"//results\n", +"printf('Quality = %.1f percent',y1)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.8: Example_8.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P=200 //psia\n", +"T=500 //F\n", +"//calculations\n", +"disp('From mollier chart,')\n", +"hi=1269 //Btu/lb\n", +"hf=1063 //Btu/lb\n", +"dh=hi-hf\n", +"y1=91\n", +"//results\n", +"printf('Quality = %.1f percent',y1)\n", +"printf('\n Change in enthalpy = %d Btu/lb',dh)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.9: Example_9.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P=200 //psia\n", +"Ts=260 //F\n", +"Tf=220 //F\n", +"m=10000 //lb\n", +"Pc=20 //psia\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"hf=188 //Btu/lb\n", +"h2=1172 //Btu/lb\n", +"Q=m*(h2-hf)\n", +"//results\n", +"printf('Heat absorption = %d Btu/hr',Q)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/6-Steam_power_plant_cycles.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/6-Steam_power_plant_cycles.ipynb new file mode 100644 index 0000000..8c343df --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/6-Steam_power_plant_cycles.ipynb @@ -0,0 +1,189 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 6: Steam power plant cycles" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P1=200 //psia\n", +"T1=600 //F\n", +"P2=2 //psia\n", +"J=778\n", +"//calculations\n", +"disp('from mollier charts,')\n", +"h1=1322 //Btu/lb\n", +"h2=974 //Btu/lb\n", +"vf2=0.01623 //cu ft per lb\n", +"hf2=94 //Btu/lb\n", +"t2=126 //F\n", +"Wtj=h1-h2\n", +"Qout=h2-hf2\n", +"Wp=(P1-P2)*vf2\n", +"Wpj=Wp/J\n", +"h3=hf2+Wpj\n", +"Qin=h1-h3\n", +"etat=((h1-h2)-Wpj)/(h1-(hf2+Wpj))\n", +"eta=((h1-h2))/(h1-(hf2))\n", +"//results\n", +"printf('Efficiency of rankine cycle = %.1f percent',etat*100)\n", +"printf('\n Efficiency of rankine cycle neglecting boiler feed pump = %.1f percent',eta*100)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"B=70 //F\n", +"P1=140 //psia\n", +"x=0.986\n", +"P2=14.7 //psia\n", +"ms=2000 //lb/hr\n", +"Ihp=80\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"hc=38 //Btu/lb\n", +"hf=324.82 //Btu/lb\n", +"hfg=868.2 //Btu/lb\n", +"h1=hf+x*hfg\n", +"Qin=ms*(h1-hc)\n", +"eta=Ihp*2545*100/(Qin)\n", +"Qw=Ihp*2545\n", +"Qr=Qin-Qw\n", +"per=Qr/Qin *100\n", +"//results\n", +"printf('Heat input to the boiler = %d Btu/hr',Qin)\n", +"printf('\n Cycle efficiency = %.1f percent',eta)\n", +"printf('\n Heat rejected to waste = %d Btu/hr or %.1f percent of Qin',Qr,per)\n", +"disp('The answer is a bit different due to rounding off error in textbook')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.3: Example_3.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"B=70 //F\n", +"P1=140 //psia\n", +"x=0.986\n", +"P2=14.7 //psia\n", +"ms=2000 //lb/hr\n", +"Ihp=80\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"hc=180 //Btu/lb\n", +"hf=324.82 //Btu/lb\n", +"hfg=868.2 //Btu/lb\n", +"h1=hf+x*hfg\n", +"Qin=ms*(h1-hc)\n", +"eta=Ihp*2545*100/(Qin)\n", +"Qw=Ihp*2545\n", +"Qr=Qin-Qw\n", +"per=Qr/Qin *100\n", +"//results\n", +"printf('Heat input to the boiler = %d Btu/hr',Qin)\n", +"printf('\n Cycle efficiency = %.2f percent',eta)\n", +"printf('\n Heat rejected to waste = %d Btu/hr or %.2f percent of Qin',Qr,per)\n", +"disp('The answer is a bit different due to rounding off error in textbook')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.4: Example_4.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"m=1.24 //lb\n", +"HV=11300 //Btu/lb\n", +"//calculations\n", +"HR=m*HV\n", +"eff=3413/HR\n", +"//results\n", +"printf('Plant heat rate = %d Btu/kw hr',HR)\n", +"printf('\n Overall efficiency = %.1f percent',eff*100)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/7-Steam_turbines.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/7-Steam_turbines.ipynb new file mode 100644 index 0000000..aad6b7a --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/7-Steam_turbines.ipynb @@ -0,0 +1,343 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 7: Steam turbines" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.10: Example_10.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P=200 //psia\n", +"T=540 //F\n", +"pow=1000 //kw\n", +"ms=16000 //lb/hr\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h1=1290 //Btu/hr\n", +"h2=940 //Btu/hr\n", +"dh=h1-h2\n", +"hf2=83 //Btu/lb\n", +"etat=(h1-h2)/(h1-hf2)\n", +"act=pow*3413/(ms*(h1-hf2))\n", +"etae=act/etat\n", +"//results\n", +"printf('Ideal thermal efficiency = %.1f percent',etat*100)\n", +"printf('\n Actual thermal efficiency = %.1f percent',act*100)\n", +"printf('\n Engine efficiency = %.1f percent',etae*100)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P1=200 //psia\n", +"T1=500 //psia\n", +"m=1 //lb /s\n", +"P4=140 //psia\n", +"P11=1 //psia\n", +"x=0.808\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h1=1268.9 //Btu/lb\n", +"h4=1234.7 //Btu/lb\n", +"V4=223.8*sqrt(h1-h4)\n", +"v4=3.584 //cu ft/lb\n", +"A4=m*v4/V4\n", +"h11=907.4 //Btu/lb\n", +"V11=223.8*sqrt(h1-h11)\n", +"vf=0.01614 //cu ft/lb\n", +"vg=333.6 //cu ft/lb\n", +"vfg=vg-vf\n", +"v11=x*vg\n", +"A11=m*v11/V11\n", +"//results\n", +"printf('Area of nozzle = %.5f sq ft',A4)\n", +"printf('\n Area of nozzle = %.4f sq ft',A11)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.4: Example_4.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P1=200 //psia\n", +"T1=500 //F\n", +"P2=1 //psia\n", +"alpha=20 //degrees\n", +"n=3600\n", +"g=32.2 //ft/s^2\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"V1=4240 //fps\n", +"Vb=V1*cosd(alpha) /2\n", +"R=Vb*60/(n*2*%pi)\n", +"work=1/32.2 *(V1*cosd(alpha))*Vb\n", +"eff=work/(V1^2 /(2*g)) *100\n", +"//results\n", +"printf('Blade velocity = %d fps',Vb)\n", +"printf('\n Blade radius = %.1f ft',R)\n", +"printf('\n Work done = %d ft-lb per lb of steam',work)\n", +"printf('\n Blade efficiency = %.1f percent',eff)\n", +"disp('The answers are a bit different due to rounding off error in textbook.')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.5: Example_5.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P1=200 //psia\n", +"T1=500 //F\n", +"P2=1 //psia\n", +"alpha=20 //degrees\n", +"n=3600\n", +"g=32.2 //ft/s^2\n", +"Vb=1200 //fps\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"V1=4240 //fps\n", +"V1x=3980 //fps\n", +"V2x=-1580 //fps\n", +"work=1/32.2 *(V1x - V2x)*Vb\n", +"eff=work/(V1^2 /(2*g)) *100\n", +"//results\n", +"printf('\n Work done = %d ft-lb per lb of steam',work)\n", +"printf('\n Blade efficiency = %.1f percent',eff)\n", +"disp('The answers are a bit different due to rounding off error in textbook.')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.6: Example_6.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P1=200 //psia\n", +"T1=500 //F\n", +"P2=1 //psia\n", +"alpha=20 //degrees\n", +"n=3600\n", +"g=32.2 //ft/s^2\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"V1=2450 //fps\n", +"Vb=V1*cosd(alpha) /2\n", +"R=Vb*60/(n*2*%pi)\n", +"work=1/32.2 *(V1*cosd(alpha))*Vb\n", +"w3=3*work\n", +"//results\n", +"printf('Blade velocity = %d fps',Vb)\n", +"printf('\n Blade radius = %.2f ft',R)\n", +"printf('\n Work done = %d ft-lb per lb of steam',w3)\n", +"disp('The answers are a bit different due to rounding off error in textbook.')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.7: Example_7.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P1=200 //psia\n", +"T1=500 //F\n", +"P2=1 //psia\n", +"alpha=20 //degrees\n", +"n=3600\n", +"g=32.2 //ft/s^2\n", +"stage=2\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"V1=4240 //fps\n", +"Vb=V1*cosd(alpha) /(2*stage)\n", +"R=Vb*60/(n*2*%pi)\n", +"V1x=3980 //fps\n", +"V2x=-1990 //fps\n", +"work1=1/g *(V1x-V2x)*Vb\n", +"work2=1/g *(-V2x)*Vb\n", +"wt=work1+work2\n", +"//results\n", +"printf('Blade velocity = %d fps',Vb)\n", +"printf('\n Blade radius = %.2f ft',R)\n", +"printf('\n Total Work done = %d ft-lb per lb of steam',wt)\n", +"disp('The answers are a bit different due to rounding off error in textbook.')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.8: Example_8.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"alpha=20 //degrees\n", +"n=3600\n", +"g=32.2 //ft/s^2\n", +"V1=500 //fps\n", +"//calculations\n", +"Vb=V1*cosd(alpha)\n", +"V1x=Vb\n", +"work=1/32.2 *(V1x)*Vb\n", +"//results\n", +"printf('Blade velocity = %d fps',Vb)\n", +"printf('\n Work done = %d ft-lb per lb of steam',work)\n", +"disp('The answers are a bit different due to rounding off error in textbook.')\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.9: Example_9.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"pow=1000 //kw\n", +"ms=16000 //lb/hr\n", +"P=200 //psia\n", +"T=540 //F\n", +"//calculations\n", +"disp('From mollier charts,')\n", +"h1=1290 //Btu/hr\n", +"h2=940 //Btu/hr\n", +"dh=h1-h2\n", +"rate=3413/dh\n", +"act=ms/pow\n", +"//results\n", +"printf('Ideal steam rate = %.2f lb per kw hr',rate)\n", +"printf('\n Actual steam rate = %d lb per kw hr',act)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/8-Steam_engines.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/8-Steam_engines.ipynb new file mode 100644 index 0000000..b87ae6b --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/8-Steam_engines.ipynb @@ -0,0 +1,122 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 8: Steam engines" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"area1=2.7\n", +"len=3.4\n", +"scale=60\n", +"area2=2.75\n", +"dia=12 //ft\n", +"d2=2.5 //ft\n", +"L=15/12 //ft\n", +"n=250 //rpm\n", +"F=600 //lb\n", +"r=3 //ft\n", +"//calculations\n", +"Ah=dia^2 *%pi/4\n", +"Ac=(dia^2 -d2^2)*%pi/4\n", +"Pih=area1/len *scale\n", +"Pic=area2/len *scale\n", +"Hihp=Pih*L*Ah*n/33000\n", +"Cihp=Pic*L*Ac*n/33000\n", +"Tihp=Hihp+Cihp\n", +"Bhp=2*%pi*r*F*n/33000\n", +"Fhp=Tihp-Bhp\n", +"eff=Bhp/Tihp *100\n", +"//results\n", +"printf('Ihp = %.1f ihp',Tihp)\n", +"printf('\n Bhp = %.1f bhp',Bhp)\n", +"printf('\n Fhp = %.1f fhp',Fhp)\n", +"printf('\n Efficiency = %.1f percent',eff)\n", +"disp('The answer is a bit different due to rounding off error in the textbook.')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"Ihp=101.1\n", +"Bhp=85.7\n", +"md=3000 //lb/hr\n", +"h1=1172 //Btu/hr\n", +"h2=180 //Btu/hr\n", +"h3=1025 //Btu/hr\n", +"//calculations\n", +"eta1=Ihp*2545/(md*(h1-h2)) *100\n", +"eta2=Bhp*2545/(md*(h1-h2)) *100\n", +"etat=(h1-h3)/(h1-h2) *100\n", +"engeff=eta1/etat *100\n", +"rate1= md/Ihp\n", +"rate2=md/Bhp\n", +"h22=h1-2545/rate1\n", +"//results\n", +"printf('Actual thermal efficiency based upon Ihp = %.2f lb per ihp hr',eta1)\n", +"printf('\n Actual thermal efficiency based upon Bhp = %.2f lb per ihp hr',eta2)\n", +"printf('\n Ideal thermal efficiency = %.2f percent ',etat)\n", +"printf('\n Engine efficiency = %.1f percent',engeff)\n", +"printf('\n Steam rate = %.2f lb per ihp hr',rate1)\n", +"printf('\n Steam rate = %.2f lb per bhp hr',rate2)\n", +"printf('\n Enthalpy of exhaust steam = %d Btu/lb of steam',h22)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Elementary_Heat_Power_by_H_L_Solberg/9-Pumps.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/9-Pumps.ipynb new file mode 100644 index 0000000..b61182a --- /dev/null +++ b/Elementary_Heat_Power_by_H_L_Solberg/9-Pumps.ipynb @@ -0,0 +1,353 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 9: Pumps" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.10: Example_10.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"sr=2\n", +"//calculations\n", +"hr=sr^2\n", +"capr=sr\n", +"hpr=sr^3\n", +"//results\n", +"printf('head is %d times the original',hr)\n", +"printf('\n capacity is %d times the original',capr)\n", +"printf('\n power is %d times the original',hpr)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.1: Example_1.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"h=200 //ft\n", +"gam=64 //lb per cu ft\n", +"//calculations\n", +"P=h*gam/144\n", +"//results\n", +"printf('Pressure = %.1f psi',P)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.2: Example_2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"P=20 //psi\n", +"gam=62.4 //lb per cu ft\n", +"//calculations\n", +"h=P*144/gam\n", +"//results\n", +"printf('height = %.1f ft',h)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.3: Example_3.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"h=3/12 //ft\n", +"gam=63.4 //lb per cu ft\n", +"gam2=0.075 //lb per cu ft\n", +"//calculations\n", +"P=h*gam\n", +"h2=P/gam2\n", +"//results\n", +"printf('Air height required = %d ft of air',h2)\n", +"disp('The answer is a bit different due to roundoff error in textbook.')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.4: Example_4.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"dif=4 //in\n", +"gam=62.4 //lb per cu ft\n", +"density=13.6 //g/cc\n", +"//calculations\n", +"pv=dif*1/12 *density*gam/144 - dif/12 *gam/144\n", +"hv=pv*144/gam\n", +"//results\n", +"printf('velocity pressure = %.2f psi',pv)\n", +"printf('\n velocity head = %.1f ft of water ',hv)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.5: Example_5.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"dif=4 //in\n", +"gam=62.4 //lb per cu ft\n", +"gam2=0.08 //lb per cu ft\n", +"//calculations\n", +"pv=dif*1/12 *gam/144\n", +"hv=pv*144/gam2\n", +"//results\n", +"printf('velocity pressure = %.3f psi',pv)\n", +"printf('\n velocity head = %.1f ft of air ',hv)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.6: Example_6.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"hw=3/12 //ft\n", +"gam1=62.4 //lb/ft^3\n", +"gam2=0.07 //lb/ft^3\n", +"g=32.2 //ft/s^2\n", +"//calculations\n", +"p=hw*gam1\n", +"hg=p/gam2\n", +"V=sqrt(2*g*hg)\n", +"//results\n", +"printf('velocity of gas = %.1f fps',V)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.7: Example_7.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"h=4 //in\n", +"den=13.6 //g/cc\n", +"Ar=1/9\n", +"A1=12 //sq in\n", +"gam=62.4 //lb/ft^3\n", +"g=32.2 //ft/s^2\n", +"//calculations\n", +"dh=(h*den-h)/12\n", +"Vr=1/Ar\n", +"V22=2*g*dh/(1-Ar^2)\n", +"V2=sqrt(V22)\n", +"A2=A1*Ar\n", +"v2=1/gam\n", +"ms=A2*V2/(v2*144)\n", +"//results\n", +"printf('Flow rate of water = %.1f lb/sec',ms)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.8: Example_8.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"mdot=8000 //lb/min\n", +"A1=1 //sq ft\n", +"A2=3/4 //sq ft\n", +"P2=50 //psi\n", +"P1=10 //psi\n", +"gam=62.4 //lb/ft^3\n", +"y2=-2 //ft\n", +"y1=-4 //ft\n", +"g=32.2 //ft/s^2\n", +"eff=0.7\n", +"//calculations\n", +"v=1/gam\n", +"cap=mdot/8.33\n", +"V1=mdot*v/A1 /60\n", +"V2=mdot*v/A2 /60\n", +"ht= (y2-y1) + (V2^2 -V1^2)/(2*g) + (P2-P1)*144/gam\n", +"Hhp=mdot*ht/33000\n", +"Php=Hhp/eff\n", +"//results\n", +"printf('Capacity = %d gpm',cap)\n", +"printf('\n Total dynamic head = %.1f ft',ht)\n", +"printf('\n Hydraulic hp = %.1f hp',Hhp)\n", +"printf('\n pump hp = %.1f hp',Php)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.9: Example_9.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//Initialization of variables\n", +"z=12 //ft\n", +"gam1=62.4 //lb/ft^3\n", +"sg=0.8\n", +"P2=100 //psia\n", +"P1=-10 //psia\n", +"mm=10000 //lb/min\n", +"//calculations\n", +"gam2=sg*gam1\n", +"p2g=P2*144/(gam2) +z\n", +"p1g=P1*144*0.491/(gam2)\n", +"ht=p2g-p1g\n", +"Hhp=mm*ht/33000\n", +"//results\n", +"printf('Total dynamic head = %.1f ft of oil',ht)\n", +"printf('\n Hydraulic hp = %.1f hp',Hhp)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} |