From f270f72badd9c61d48f290c3396004802841b9df Mon Sep 17 00:00:00 2001 From: kinitrupti Date: Fri, 12 May 2017 18:53:46 +0530 Subject: Removed duplicates --- .../Chapter_10.ipynb | 457 +++++++++++++++ .../Chapter_11.ipynb | 307 ++++++++++ .../Chapter_12.ipynb | 109 ++++ .../Chapter_13.ipynb | 623 +++++++++++++++++++++ .../Chapter_14.ipynb | 258 +++++++++ .../Chapter_15.ipynb | 171 ++++++ .../Chapter_16.ipynb | 178 ++++++ .../Chapter_17.ipynb | 289 ++++++++++ .../Chapter_18.ipynb | 131 +++++ .../Chapter_19.ipynb | 118 ++++ .../Chapter_2.ipynb | 141 +++++ .../Chapter_21.ipynb | 261 +++++++++ .../Chapter_22.ipynb | 256 +++++++++ .../Chapter_23.ipynb | 302 ++++++++++ .../Chapter_24.ipynb | 76 +++ .../Chapter_25.ipynb | 294 ++++++++++ .../Chapter_3.ipynb | 240 ++++++++ .../Chapter_5.ipynb | 225 ++++++++ .../Chapter_6.ipynb | 154 +++++ .../Chapter_8.ipynb | 244 ++++++++ .../README.txt | 10 + .../screenshots/chap23.png | Bin 0 -> 51627 bytes .../screenshots/chap24.png | Bin 0 -> 48124 bytes .../screenshots/chap25.png | Bin 0 -> 41706 bytes 24 files changed, 4844 insertions(+) create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_10.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_11.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_12.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_13.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_14.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_15.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_16.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_17.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_18.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_19.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_2.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_21.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_22.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_23.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_24.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_25.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_3.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_5.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_6.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_8.ipynb create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/README.txt create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap23.png create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap24.png create mode 100755 Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap25.png (limited to 'Mechanical_Engineering_Thermodynamics_by_Mooney_D_A') diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_10.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_10.ipynb new file mode 100755 index 00000000..aa1e31ac --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_10.ipynb @@ -0,0 +1,457 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 10 - Tabulated Properties : Steam Tables" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 155" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the work and heat of vaporization of water\n", + "#initialization of varaibles\n", + "P=100. #psia\n", + "hfg=888.8 #B/lb\n", + "#calculations\n", + "print '%s' %(\"From steam tables,\")\n", + "vg=4.432 #cu ft/lb\n", + "vf=0.001774 #cu ft/lb\n", + "W=P*(vg-vf)*144\n", + "ufg=807.1 #B/lb\n", + "W=hfg-ufg\n", + "sfg=1.1286\n", + "Q=788*sfg #Heat\n", + "#results\n", + "print '%s %.1f %s' %(\"Work done = \",W,\"B/lb\")\n", + "print '%s %d %s' %(\"\\n Heat of vaporization of water =\",Q,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From steam tables,\n", + "Work done = 81.7 B/lb\n", + "\n", + " Heat of vaporization of water = 889 B/lb\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 156" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#caclulate the pressure, temperature and enthalpy\n", + "#initialization of varaibles\n", + "s=1.6315 #B/lb R\n", + "#calculations\n", + "print '%s' %(\"From table 1 \")\n", + "h=1180.6 #B/lb\n", + "t=302.92 #F\n", + "p=70 #psia\n", + "#results\n", + "print '%s %d %s' %(\"Pressure = \",p,\"psia\")\n", + "print '%s %.2f %s' %(\"\\n Temperature =\",t,\"F\")\n", + "print '%s %.1f %s' %(\"\\n Enthalpy =\",h,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 1 \n", + "Pressure = 70 psia\n", + "\n", + " Temperature = 302.92 F\n", + "\n", + " Enthalpy = 1180.6 B/lb\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 156" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the Internal energy of the gas\n", + "#initialization of varaibles\n", + "T=250. #F\n", + "print '%s' %(\"From table 1,\")\n", + "p=29.825 #psia\n", + "hg=1164 #B/lb\n", + "vg=13.821 #cu ft/lb\n", + "#calculations\n", + "ug=hg-(p)*144*vg/778.\n", + "#results\n", + "print '%s %.1f %s' %(\"Internal energy of the gas = \",ug,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 1,\n", + "Internal energy of the gas = 1087.7 B/lb\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4 - Pg 158" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the specific volume, Enthalpy and Entropy of the gas\n", + "#initialization of varaibles\n", + "x=0.4\n", + "P=100. #psia\n", + "#calculations\n", + "y=1-x\n", + "print '%s' %(\"From table 2,\")\n", + "vf=0.01774\n", + "vg=4.432\n", + "vx=x*vf+y*vg\n", + "hf=298.4\n", + "hfg=888.8\n", + "hx=hf+y*hfg\n", + "sg=1.6026\n", + "sfg=1.1286\n", + "sx=sg-x*sfg\n", + "#results\n", + "print '%s %.3f %s' %(\"Specific volume = \",vx,\"cu ft/lb\")\n", + "print '%s %.1f %s' %(\"\\n Enthalpy =\",hx,\"B/lb\")\n", + "print '%s %.4f %s' %(\"\\n Entropy = \",sx,\"B/lb R\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 2,\n", + "Specific volume = 2.666 cu ft/lb\n", + "\n", + " Enthalpy = 831.7 B/lb\n", + "\n", + " Entropy = 1.1512 B/lb R\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5 - Pg 158" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the enthalpy of the gas\n", + "#initialization of varaibles\n", + "x=0.97\n", + "P=100. #psia\n", + "#calculations\n", + "print '%s' %(\"From table 2,\")\n", + "hf=298.4\n", + "hfg=888.8\n", + "hx=hf+x*hfg\n", + "hg=1187.2\n", + "hx2=hg-(1-x)*hfg\n", + "#results\n", + "print '%s %d %s' %(\"\\n Enthalpy =\",hx,\"B/lb\")\n", + "print '%s %.1f %s' %(\"\\n Accurate Enthalpy =\",hx2,\"B/lb\")" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 2,\n", + "\n", + " Enthalpy = 1160 B/lb\n", + "\n", + " Accurate Enthalpy = 1160.5 B/lb\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6 - Pg 159" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the specific volume and enthalpy of the process\n", + "#initialization of varaibles\n", + "s=1.7050 #B/lb R\n", + "#calculations\n", + "print '%s' %(\"From table 2,\")\n", + "sx=1.7050\n", + "sg=1.7549\n", + "sfg=1.4415\n", + "dx=(sg-sx)/sfg\n", + "hg=1150.8\n", + "hfg=969.7\n", + "hx=hg-dx*hfg\n", + "vg=26.29\n", + "vfg=26.27\n", + "vx=vg-dx*vfg\n", + "#results\n", + "print '%s %.2f %s' %(\"Specific volume = \",vx,\"cu ft/lb\")\n", + "print '%s %.1f %s' %(\"\\n Enthalpy =\",hx,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 2,\n", + "Specific volume = 25.38 cu ft/lb\n", + "\n", + " Enthalpy = 1117.2 B/lb\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7 - Pg 159" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the enthalpy required\n", + "#initialization of varaibles\n", + "P=150 #psia\n", + "T=400 #F\n", + "#calculations\n", + "print '%s' %(\"From table 3,\")\n", + "h=1219.4 #B/lb\n", + "#results\n", + "print '%s %.1f %s' %(\"Enthalpy =\",h,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 3,\n", + "Enthalpy = 1219.4 B/lb\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8 - Pg 159" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the temperature of steam\n", + "#initialization of varaibles\n", + "en=1303.7 #B/lb\n", + "P=300 #psia\n", + "#calculations\n", + "print '%s' %(\"Given hg is less than h, steam is superheated. T=580 F\")\n", + "T=580 #F\n", + "#results\n", + "print '%s %d %s' %(\"Temperature =\",T,\"F\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Given hg is less than h, steam is superheated. T=580 F\n", + "Temperature = 580 F\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9 - Pg 160" + ] + }, + { + "cell_type": "raw", + "metadata": {}, + "source": [ + "This is a theoretical example." + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10 - Pg 161" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the Enthalpy and Volume of the gas\n", + "#initialization of varaibles\n", + "import math\n", + "T=100. #F\n", + "P=1000. #psia\n", + "#calculations\n", + "print '%s' %(\"From table 4\")\n", + "dvf=-5.1*math.pow(10,-5)\n", + "dhf=2.7\n", + "vf=0.01613\n", + "hf=67.97\n", + "v=vf+dvf\n", + "h=hf+dhf\n", + "#results\n", + "print '%s %.2f %s' %(\"Enthalpy =\",h,\"B/lb\")\n", + "print '%s %.5f %s' %(\"\\n Volume = \",v,\"cu ft/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 4\n", + "Enthalpy = 70.67 B/lb\n", + "\n", + " Volume = 0.01608 cu ft/lb\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 11 - Pg 164" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the Quality of steam\n", + "#initialization of varaibles\n", + "h1=1183.2 #B/lb\n", + "hg=1198.4 #B/lb\n", + "hfg=843.\n", + "#calculations\n", + "x=1- (hg-h1)/hfg\n", + "#results\n", + "print '%s %.3f' %(\"Quality =\",x)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Quality = 0.982\n" + ] + } + ], + "prompt_number": 11 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_11.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_11.ipynb new file mode 100755 index 00000000..fc43bf01 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_11.ipynb @@ -0,0 +1,307 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 11 - Properties of Gases" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 172" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the volume of the appratus\n", + "#initialization of varaibles\n", + "P1=15. #psia\n", + "T1=80+460. #R\n", + "dm=3. #lb\n", + "T2=75+460. #R\n", + "P2=25. #psia\n", + "#calculations\n", + "mratio=P1*T2/(P2*T1)\n", + "m2=dm/(1-mratio)\n", + "V2=m2*55.16*T2/(P2*144.)\n", + "#results\n", + "print '%s %.1f %s' %(\"Volume of the apparatus =\",V2,\"cu ft\")\n", + "print '%s' %(\"The answer is a bit different due to rounding off error in textbook\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Volume of the apparatus = 60.6 cu ft\n", + "The answer is a bit different due to rounding off error in textbook\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 176" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the specific heats at constant pressure and volume\n", + "#initialization of varaibles\n", + "R=48.3 #ft lb/lb R\n", + "k=1.4\n", + "#calculations\n", + "dc=R/778.\n", + "cp=k*dc/(k-1)\n", + "cv=cp/k\n", + "#results\n", + "print '%s %.3f %s' %(\"Specific heat at constant volume = \",cv,\"B/lb R\")\n", + "print '%s %.3f %s' %(\"\\n Specific heat at constant pressure = \",cp,\"B/lb R\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Specific heat at constant volume = 0.155 B/lb R\n", + "\n", + " Specific heat at constant pressure = 0.217 B/lb R\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 177" + ] + }, + { + "cell_type": "raw", + "metadata": {}, + "source": [ + "This is a theoretical example." + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + " Example 4 - Pg 180" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the enthalpy generated\n", + "#initialization of varaibles\n", + "import math\n", + "P1=100. #psia\n", + "P2=10. #psia\n", + "T1=140. +460 #R\n", + "g=1.4\n", + "cp=0.248\n", + "#calculations\n", + "dh=g*55.16*T1*(math.pow((P2/P1),((g-1)/g)) -1)/(g-1)*0.01286 #0.01286 is the conversion factor\n", + "T2=T1*math.pow((P2/P1),((g-1)/g))\n", + "dh2=cp*(T2-T1)\n", + "#results\n", + "print '%s %d %s' %(\"In method 1, Enthalpy = \",dh,\"Btu/lb\")\n", + "print '%s %.1f %s' %(\"\\n In method 2, Enthalpy =\",dh2,\"ft lb/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "In method 1, Enthalpy = -718 Btu/lb\n", + "\n", + " In method 2, Enthalpy = -71.7 ft lb/lb\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5a - Pg 186" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the changes in Enthalpy and Specific Volume\n", + "#initialization of varaibles\n", + "import math\n", + "P1=100. #psia\n", + "T1=2000.+460 #R\n", + "P2=15. #psia\n", + "g=1.4\n", + "cp=0.24\n", + "#calculations\n", + "v1=53.34*T1/(P1*144.)\n", + "v2=53.34*T1*math.pow((P1/P2),(1/g)) /(P1*144)\n", + "T2=T1*P2*v2/(P1*v1)\n", + "dh=cp*(T2-T1)\n", + "dv=v2-v1\n", + "#results\n", + "print '%s %d %s' %(\"Change in enthalpy =\",dh,\"B/lb\")\n", + "print '%s %.1f %s' %(\"\\n Specific volume change = \",dv,\"cu ft/lb\")\n", + "print '%s' %(\"The answer is a bit different due to rounding off error in textbook\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Change in enthalpy = -247 B/lb\n", + "\n", + " Specific volume change = 26.2 cu ft/lb\n", + "The answer is a bit different due to rounding off error in textbook\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5b - Pg 186" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the changes in Enthalpy and Specific Volume\n", + "#initialization of varaibles\n", + "import math\n", + "P1=100. #psia\n", + "T1=2000.+460 #R\n", + "P2=15. #psia\n", + "g=1.4\n", + "cp=0.276\n", + "cv=0.207\n", + "T2=1520. #R\n", + "#calculations\n", + "k=cp/cv\n", + "v1=53.34*T1/(P1*144)\n", + "v2=v1*math.pow((P1/P2),(1/k))\n", + "dh=cp*(T2-T1)\n", + "dv=v2-v1\n", + "#results\n", + "print '%s %d %s' %(\"Enthalpy change =\",dh,\"B/lb\")\n", + "print '%s %.1f %s' %(\"\\n Volume change =\",dv,\"cu t/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Enthalpy change = -259 B/lb\n", + "\n", + " Volume change = 28.7 cu t/lb\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5c - Pg 186" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the changes in Enthalpy and Specific Volume\n", + "#initialization of varaibles\n", + "import math\n", + "P1=100. #psia\n", + "T1=2000.+460 #R\n", + "P2=15. #psia\n", + "g=1.4\n", + "cp=0.276\n", + "cv=0.207\n", + "T2=1520. #R\n", + "#calculations\n", + "h1=634.4\n", + "pr1=407.3\n", + "pr2=pr1*P2/P1\n", + "print '%s' %(\"From table 1,\")\n", + "T2=1535. #R\n", + "h2=378.44\n", + "dh=h2-h1\n", + "v1=53.34*T1/(P1*144)\n", + "v2=53.34*T2/(P2*144.)\n", + "dv=v2-v1\n", + "#results\n", + "print '%s %.2f %s' %(\"Enthalpy change = \",dh,\"B/lb\")\n", + "print '%s %.1f %s' %(\"\\n Volume change =\",dv,\"cu ft/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 1,\n", + "Enthalpy change = -255.96 B/lb\n", + "\n", + " Volume change = 28.8 cu ft/lb\n" + ] + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_12.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_12.ipynb new file mode 100755 index 00000000..82670034 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_12.ipynb @@ -0,0 +1,109 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 12 - Properties of Gaseous mixtures" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 201" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the partial pressure and specific humidity\n", + "#initialization of varaibles\n", + "P=15. #psia\n", + "T2=70.+460 #R\n", + "T1=55.+460 #R\n", + "#calculations\n", + "pw=0.2141\n", + "pA=P-pw\n", + "mratio=pA*29/(pw*18.)\n", + "mAbym=mratio/(1+mratio)\n", + "mwbym=1/(1+mratio)\n", + "pg=0.3631 #psia\n", + "phi=pw/pg\n", + "gamma=1/mratio\n", + "#results\n", + "print '%s %.2f %s' %(\"Partial pressure of water vapor = \",pA,\"psia\")\n", + "print '%s %.4f %s' %(\"\\n Specific humidity = \",gamma,\"lb vapor/lb air\")\n", + "print '%s %.2f' %(\"\\n Relative humidity = \",phi)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Partial pressure of water vapor = 14.79 psia\n", + "\n", + " Specific humidity = 0.0090 lb vapor/lb air\n", + "\n", + " Relative humidity = 0.59\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 201" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the mass of watervapor per pound of dry air\n", + "#initialization of varaibles\n", + "rh=0.75\n", + "pg=0.5069\n", + "inc=10 #in\n", + "pA=29.50 #psia\n", + "#calculations\n", + "pw=rh*pg\n", + "p=(29.50+ inc/13.6)*0.491\n", + "pA=p-pw\n", + "mratio=pw*18/(pA*29.)\n", + "#results\n", + "print '%s %.4f %s' %(\"Pounds of water vapor enter the surface per pound of dry air =\",mratio,\"lb vapor/lb air\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Pounds of water vapor enter the surface per pound of dry air = 0.0163 lb vapor/lb air\n" + ] + } + ], + "prompt_number": 3 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_13.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_13.ipynb new file mode 100755 index 00000000..bcedae7c --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_13.ipynb @@ -0,0 +1,623 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 13 - Process Calculations for Stationary Systems" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1a - Pg 203" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the final temperature and heat transferred from the steam\n", + "#initialization of varaibles\n", + "P1=100. #psia\n", + "T1=500+460. #R\n", + "v=10. #cu ft\n", + "P2=50. #psia\n", + "cv=0.172\n", + "R=53.34\n", + "m=2.81 #lb\n", + "#calculations\n", + "T2=T1*P2/P1\n", + "Q1=P1*144*v*cv*(T2-T1)/(R*T1)\n", + "u1=165.26 #Btu/lb\n", + "u2=81.77 #Btu/lb\n", + "du=u2-u1\n", + "Q2=m*du\n", + "#results\n", + "print '%s' %(\"Case 1,\")\n", + "print '%s %d %s' %(\"\\n Final temperature of the steam = \",T2,\"R\")\n", + "print '%s %.1f %s' %(\"\\n Heat transferred =\",Q1,\"Btu\")\n", + "print '%s %.1f %s' %(\"\\n Heat transferred in case 2 =\",Q2,\"Btu\")\n", + "print '%s' %(\"\\n The answer may be a bit different due to rounding off error in the textbook\")" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Case 1,\n", + "\n", + " Final temperature of the steam = 480 R\n", + "\n", + " Heat transferred = -232.2 Btu\n", + "\n", + " Heat transferred in case 2 = -234.6 Btu\n", + "\n", + " The answer may be a bit different due to rounding off error in the textbook\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1b - Pg 204" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the final temperature and heat transferred\n", + "#initialization of varaibles\n", + "P1=100. #psia\n", + "T1=500.+460 #R\n", + "V=10. #cu ft\n", + "P2=50. #psia\n", + "cv=0.172\n", + "R=53.34\n", + "v=5.589 #cu ft/lb\n", + "#calculations\n", + "m=V/v\n", + "x2=(v-0.017)/8.498\n", + "print '%s' %(\"From table 2,\")\n", + "T2=281.01 #F\n", + "h1=1279.1\n", + "u1=h1-144*P1*v/778.\n", + "uf=249.93\n", + "ufg=845.4\n", + "u2=uf+x2*ufg\n", + "Q=m*(u2-u1)\n", + "#results\n", + "print '%s %.2f %s' %(\"Final temperature =\",T2,\"F\")\n", + "print '%s %d %s' %(\"\\n Heat transferred =\",Q,\"Btu\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 2,\n", + "Final temperature = 281.01 F\n", + "\n", + " Heat transferred = -664 Btu\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2a - Pg 205" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the final temperature, Enthalpy and heat transferred in the process\n", + "#initialization of varaibles\n", + "T1=350+460 #R\n", + "v1=6 #cu ft/lb\n", + "m=1 #lb\n", + "R=53.34\n", + "v2=2*v1\n", + "cp=0.24\n", + "#calculations\n", + "P=R*T1/(v1*144)\n", + "W=P*144*(v2-v1)\n", + "T2=T1*v2/v1\n", + "Q=cp*(T2-T1)\n", + "h1=194.25\n", + "h2=401.09\n", + "dh=h2-h1\n", + "#results\n", + "print '%s %d %s' %(\"Final temperature =\",T2-460,\"F\")\n", + "print '%s %.2f %s' %(\"\\n Enthalpy =\",dh,\"B/lb\")\n", + "print '%s %d %s' %(\"\\n Heat =\",Q,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Final temperature = 1160 F\n", + "\n", + " Enthalpy = 206.84 B/lb\n", + "\n", + " Heat = 194 B/lb\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2b - Pg 206" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the final temperature, Enthalpy and heat transferred in the process\n", + "#initialization of varaibles\n", + "T1=350+460. #R\n", + "v1=6 #cu ft/lb\n", + "m=1 #lb\n", + "R=53.34\n", + "v2=2*v1\n", + "cp=0.24\n", + "#calculations\n", + "print '%s' %(\"From steam tables,\")\n", + "vg=3.342 #cu ft/lb\n", + "P1=77.5 #psia\n", + "P2=P1\n", + "h1=1204.8 #B/lb\n", + "v2=2*v1\n", + "T2=1106 #F\n", + "h2=1586.7 #B/lb\n", + "Q=h2-h1\n", + "W=P1*144*(v2-v1)\n", + "#results\n", + "print '%s %d %s' %(\"Final temperature =\",T2,\"F\")\n", + "print '%s %d %s' %(\"\\n Work =\",W,\"ft lb/lb\")\n", + "print '%s %.1f %s' %(\"\\n Heat =\",Q,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From steam tables,\n", + "Final temperature = 1106 F\n", + "\n", + " Work = 66960 ft lb/lb\n", + "\n", + " Heat = 381.9 B/lb\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3a - Pg 206" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the final pressure, work done and change in internal energy of the process\n", + "#initialization of varaibles\n", + "import math\n", + "T1=400+460. #R\n", + "P1=50. #psia\n", + "ratio=1/10.\n", + "R=53.34\n", + "#calculations\n", + "P2=P1/ratio\n", + "W=R*T1*math.log(ratio)\n", + "du=0\n", + "#results\n", + "print '%s %d %s' %(\"Final pressure =\",P2,\"psia\")\n", + "print '%s %.1f %s' %(\"\\n Work done =\",W,\"B/lb\")\n", + "print '%s %d %s' %(\"\\n Change in Internal energy \",du,\"Btu/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Final pressure = 500 psia\n", + "\n", + " Work done = -105625.1 B/lb\n", + "\n", + " Change in Internal energy 0 Btu/lb\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3b - Pg 207" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the final pressure, work done and change in internal energy of the process\n", + "#initialization of varaibles\n", + "T1=400+460. #R\n", + "P1=50. #psia\n", + "ratio=1/10.\n", + "R=53.34\n", + "v1=10.065 #cu ft/lb\n", + "vfg=1.8447 #cu ft/lb\n", + "vg=1.8633 #cu ft/lb\n", + "#calculations\n", + "v2=v1*ratio\n", + "dx=(v2-vg)/vfg\n", + "P2=247.3 #psia\n", + "print '%s' %(\"From steam tables,\")\n", + "u2=773. #B/lb\n", + "u1=1141.6 #B/lb\n", + "du=u2-u1\n", + "s1=1.7349 #B/lb R\n", + "s2=1.082 #B/lb R\n", + "W=T1*(s2-s1) - du\n", + "#results\n", + "print '%s %.1f %s' %(\"Final pressure = \",P2,\"psia\")\n", + "print '%s %d %s' %(\"\\n Work done =\",W,\"B/lb\")\n", + "print '%s %d %s' %(\"\\n Change in Internal energy = \",du,\"B/lb \")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From steam tables,\n", + "Final pressure = 247.3 psia\n", + "\n", + " Work done = -192 B/lb\n", + "\n", + " Change in Internal energy = -368 B/lb \n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4a - Pg 208" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate final specific volume and work per pound of fluid in the problem\n", + "#initialization of varaibles\n", + "import math\n", + "P1=150. #psia\n", + "T1=400.+460 #R\n", + "P2=15. #psia\n", + "g=1.4\n", + "R=53.34\n", + "#calculations\n", + "Tratio=math.pow((P2/P1),((g-1)/g))\n", + "W=53.34*T1*(Tratio-1)/(1-g)\n", + "T2=T1*Tratio\n", + "v2=R*T2/(P2*144)\n", + "u1=147.50\n", + "Pr1=7.149\n", + "Pr2=Pr1*P2/P1\n", + "print '%s' %(\"From tables,\")\n", + "Pr=0.7149\n", + "T2=447. #R\n", + "u2=76.13 #B/lb\n", + "W=-(u2-u1)\n", + "v2=R*T2/(P2*144)\n", + "#results\n", + "print '%s %.1f %s' %(\"Final specific volume =\",v2,\"cu ft/lb\")\n", + "print '%s %.1f %s' %(\"\\n Work per pound of fluid =\",W,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From tables,\n", + "Final specific volume = 11.0 cu ft/lb\n", + "\n", + " Work per pound of fluid = 71.4 B/lb\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4b - Pg 209" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate final specific volume and work per pound of fluid in the problem\n", + "#initialization of varaibles\n", + "print '%s' %(\"From Steam tables,\")\n", + "h1=1219.4\n", + "P1=150 #psia\n", + "v1=0.59733 #cu ft/lb\n", + "s1=1.5995 #B/lb R\n", + "#calculations\n", + "u1=h1-P1*v1\n", + "sg=1.7549\n", + "sfg=1.4415\n", + "s2=s1\n", + "dx=(sg-s2)/sfg\n", + "u2=981.3\n", + "W=u1-u2\n", + "v2=23.48\n", + "#results\n", + "print '%s %.2f %s' %(\"Final specific volume =\",v2,\"cu ft/lb\")\n", + "print '%s %.1f %s' %(\"\\n Work per pound of fluid = \",W,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From Steam tables,\n", + "Final specific volume = 23.48 cu ft/lb\n", + "\n", + " Work per pound of fluid = 148.5 B/lb\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5a - Pg 210" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the final specific volume, temperature and the net heat transferred in the process\n", + "#initialization of varaibles\n", + "import math\n", + "P1=150. #psia\n", + "T1=400.+460 #R\n", + "P2=15. #psia\n", + "n=1.15\n", + "cv=0.172\n", + "R=53.34\n", + "#calculations\n", + "v2=R*T1*math.pow((P1/P2),(1/n)) /(P1*144.)\n", + "v1=R*T1/(P1*144.)\n", + "T2=T1*P2*v2/(P1*v1)\n", + "Q=(cv - 0.458)*(T2-T1)\n", + "#results\n", + "print '%s %.1f %s' %(\"Final specific volume = \",v2,\"cu ft/lb\")\n", + "print '%s %d %s' %(\"\\n Final temperature =\",T2,\"R\")\n", + "print '%s %.1f %s' %(\"\\n Heat transferred =\",Q,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Final specific volume = 15.7 cu ft/lb\n", + "\n", + " Final temperature = 636 R\n", + "\n", + " Heat transferred = 63.8 B/lb\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5b - Pg 211" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the final specific volume, temperature and the net heat transferred in the process\n", + "#initialization of varaibles\n", + "import math\n", + "print '%s' %(\"From table 3,\")\n", + "v1=3.223 #cu ft/lb\n", + "P1=150. #psia\n", + "T1=400.+460 #R\n", + "P2=15. #psia\n", + "n=1.15\n", + "#calculations\n", + "v2=v1*math.pow((P1/P2),(1/n))\n", + "T2=213 #F\n", + "W=144*(P2*v2-P1*v1)*0.00129/(1-n)\n", + "u1=1129.8 #B/lb\n", + "v2=23.9\n", + "vg=26.29\n", + "vfg=26.27\n", + "dx=(vg-v2)/vfg\n", + "u2=996.1\n", + "Q=(u2-u1)+W\n", + "#results\n", + "print '%s %.1f %s' %(\"Final specific volume =\",v2,\"cu ft/lb\")\n", + "print '%s %d %s' %(\"\\n Final temperature =\",T2,\" F\")\n", + "print '%s %.1f %s' %(\"\\n Heat transferred =\",Q,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 3,\n", + "Final specific volume = 23.9 cu ft/lb\n", + "\n", + " Final temperature = 213 F\n", + "\n", + " Heat transferred = 21.6 B/lb\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6a - Pg 212" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the final specific volume, temperature and the net work done in the process\n", + "#initialization of varaibles\n", + "v2=15.7 #cu ft/lb\n", + "T2=640 #R\n", + "cv=0.172\n", + "T1=400+460. #R\n", + "#calculations\n", + "du=cv*(T2-T1)\n", + "W=-du\n", + "#results\n", + "print '%s %.1f %s' %(\"Final specific volume =\",v2,\"cu ft/lb\")\n", + "print '%s %d %s' %(\"\\n Final temperature =\",T2,\"R \")\n", + "print '%s %.1f %s' %(\"\\n Work done =\",W,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Final specific volume = 15.7 cu ft/lb\n", + "\n", + " Final temperature = 640 R \n", + "\n", + " Work done = 37.8 B/lb\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6b - Pg 213" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the final specific volume, temperature and the net work done in the process\n", + "#initialization of varaibles\n", + "print '%s' %(\"From steam tables,\")\n", + "T2=213 #F\n", + "v2=23.9 #cu ft/lb\n", + "W=133.7 #B/lb\n", + "#results\n", + "print '%s %.1f %s' %(\"Final specific volume =\",v2,\"cu ft/lb\")\n", + "print '%s %d %s' %(\"\\n Final temperature =\",T2,\"F\")\n", + "print '%s %.1f %s' %(\"\\n Work done = \",W,\" B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From steam tables,\n", + "Final specific volume = 23.9 cu ft/lb\n", + "\n", + " Final temperature = 213 F\n", + "\n", + " Work done = 133.7 B/lb\n" + ] + } + ], + "prompt_number": 12 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_14.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_14.ipynb new file mode 100755 index 00000000..67743161 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_14.ipynb @@ -0,0 +1,258 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 14 - Vapor Cycles : The Rankine Cycle" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 219" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the heat supplied, work done, efficiency and steam rate in the process\n", + "#initialization of varaibles\n", + "P1=200. #psia\n", + "T1=750+460. #R\n", + "P2=1. #psia\n", + "#calculations\n", + "print '%s' %(\"From steam tables,\")\n", + "h1=1399.2\n", + "h2=976\n", + "h3=69.70\n", + "v3=0.01614\n", + "dh3=v3*(P1-P2)*144/778.\n", + "h4=dh3+h3\n", + "Q1=h1-h4\n", + "Wt=h1-h2\n", + "Wp=h4-h3\n", + "eta=(Wt-Wp)/Q1\n", + "w=2545./Wt\n", + "#results\n", + "print '%s %.1f %s' %(\"Heat supplied = \",Q1,\"B/lb\")\n", + "print '%s %d %s' %(\"\\n Turbine work =\",Wt,\"B/lb\")\n", + "print '%s %.3f %s' %(\"\\n Pump work =\",Wp,\"B/lb\")\n", + "print '%s %.3f' %(\"\\n Efficiency =\",eta)\n", + "print '%s %.2f %s' %(\"\\n Steam rate =\",w,\"lb steam per hr\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From steam tables,\n", + "Heat supplied = 1328.9 B/lb\n", + "\n", + " Turbine work = 423 B/lb\n", + "\n", + " Pump work = 0.594 B/lb\n", + "\n", + " Efficiency = 0.318\n", + "\n", + " Steam rate = 6.01 lb steam per hr\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 228" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the engine efficiency\n", + "#initialization of varaibles\n", + "h1=1399.2 #B/lb\n", + "h2s=976 #B/lb\n", + "wt=8. #lb /hp hr\n", + "#calculations\n", + "Wt=2545/wt\n", + "etaT=Wt/(h1-h2s)\n", + "h2=h1-Wt\n", + "#results\n", + "print '%s %.3f' %(\"Engine efficiency =\",etaT)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Engine efficiency = 0.752\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3a - Pg 234" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the backwork, efficiency and steam rate for both the cycles\n", + "#initialization of varaibles\n", + "P1=200. #psia\n", + "P2=1. #psia\n", + "e=0.7\n", + "#calculations\n", + "h1=1198.4\n", + "h2s=863.5\n", + "h3r=69.7\n", + "h4r=70.3\n", + "h3c=300.7\n", + "h4c=355.4\n", + "print '%s' %(\"For Rankine cycle, \")\n", + "Wtr=h1-h2s\n", + "Q1r=h1-h4r\n", + "Wpr=h4r-h3r\n", + "Wnetr=Wtr-Wpr\n", + "eta1=(Wtr-Wpr)/Q1r\n", + "wr=2545./Wtr\n", + "print '%s %.1f %s' %(\"Back work =\",Wnetr,\"B/lb\")\n", + "print '%s %.3f' %(\"\\n Efficiency =\",eta1)\n", + "print '%s %.1f %s' %(\"\\n Steam rate =\",wr,\"lb/hp hr\")\n", + "print '%s' %(\"\\n For carnot cycle,\")\n", + "Wtc=h1-h2s\n", + "Q1c=h1-h4c\n", + "Wpc=h4c-h3c\n", + "Wnetc=Wtc-Wpc\n", + "eta2=(Wtc-Wpc)/Q1c\n", + "wc=9.1\n", + "print '%s %.1f %s' %(\"\\n Back work =\",Wnetc,\"B/lb\")\n", + "print '%s %.3f' %(\"\\n Efficiency = \",eta2)\n", + "print '%s %.2f %s' %(\"\\n Steam rate =\",wc,\"lb/hp hr\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For Rankine cycle, \n", + "Back work = 334.3 B/lb\n", + "\n", + " Efficiency = 0.296\n", + "\n", + " Steam rate = 7.6 lb/hp hr\n", + "\n", + " For carnot cycle,\n", + "\n", + " Back work = 280.2 B/lb\n", + "\n", + " Efficiency = 0.332\n", + "\n", + " Steam rate = 9.10 lb/hp hr\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3b - Pg 234" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the backwork, efficiency and steam rate for both the cycles\n", + "#initialization of varaibles\n", + "P1=200. #psia\n", + "P2=1. #psia\n", + "e=0.7\n", + "#calculations\n", + "h1=1198.4\n", + "h2s=863.5\n", + "h3r=69.7\n", + "h4r=70.3\n", + "h3c=300.7\n", + "h4c=355.4\n", + "print '%s' %(\"For Rankine cycle with actual machines, \")\n", + "Wtr=e*(h1-h2s)\n", + "Q1r=(h1-h4r)\n", + "Wpr=(h4r-h3r)/e\n", + "Wnetr=Wtr-Wpr\n", + "eta1=(Wtr-Wpr)/Q1r\n", + "wr=2545./Wtr\n", + "print '%s %.1f %s' %(\"Back work =\",Wnetr,\"B/lb\")\n", + "print '%s %.3f '%(\"\\n Efficiency = \",eta1)\n", + "print '%s %.1f %s' %(\"\\n Steam rate =\",wr,\"lb/hp hr\")\n", + "print '%s' %(\"For carnot cycle,\")\n", + "Wtc=e*(h1-h2s)\n", + "Q1c=h1-h4c\n", + "Wpc=(h4c-h3c)/e\n", + "Wnetc=Wtc-Wpc\n", + "eta2=(Wtc-Wpc)/Q1c\n", + "wc=16.2\n", + "print '%s %.1f %s' %(\"Back work =\",Wnetc,\"B/lb\")\n", + "print '%s %.3f' %(\"\\n Efficiency = \",eta2)\n", + "print '%s %.1f %s' %(\"\\n Steam rate =\",wc,\"lb/hp hr\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For Rankine cycle with actual machines, \n", + "Back work = 233.6 B/lb\n", + "\n", + " Efficiency = 0.207 \n", + "\n", + " Steam rate = 10.9 lb/hp hr\n", + "For carnot cycle,\n", + "Back work = 156.3 B/lb\n", + "\n", + " Efficiency = 0.185\n", + "\n", + " Steam rate = 16.2 lb/hp hr\n" + ] + } + ], + "prompt_number": 4 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_15.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_15.ipynb new file mode 100755 index 00000000..546b25c5 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_15.ipynb @@ -0,0 +1,171 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 15 - Vapor Cycles : More Efficient Cycles" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 243" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the improvement in efficiency\n", + "#initialization of varaibles\n", + "e=0.85\n", + "print '%s' %(\"From Mollier chart and table 3,\")\n", + "h1=1474.5 #B/lb\n", + "s1=1.5603 #B/lb R\n", + "h2s=1277.5 #B/lb\n", + "#calculations\n", + "h2=h1-e*(h1-h2s)\n", + "h3=1522.4 #B/lb\n", + "h4s=948 #B/lb\n", + "h4=h3-e*(h3-h4s)\n", + "h5=47.6 #B/lb\n", + "h6=53.5 #B/lb\n", + "h7s=840 #B/lb\n", + "h7=h1-e*(h1-h7s)\n", + "h8=1493.2 #B/lb\n", + "h9s=866 #B/lb\n", + "h9=h8-e*(h8-h9s)\n", + "h11=51.5 #B/lb\n", + "eta1=0.401\n", + "eta2=0.375\n", + "eta3=0.366\n", + "IE1=(eta1-eta2)/eta2\n", + "IE2=(eta1-eta3)/eta3\n", + "#results\n", + "print '%s %.1f %s' %(\"Improvement in efficiency =\",IE1*100,\"percent\")\n", + "print '%s %.1f %s' %(\"\\nImprovement in efficiency in case 2=\",IE2*100,\"percent\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From Mollier chart and table 3,\n", + "Improvement in efficiency = 6.9 percent\n", + "\n", + "Improvement in efficiency in case 2= 9.6 percent\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2a - Pg 257" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the fraction of energy supplied \n", + "#initialization of varaibles\n", + "print '%s' %(\"From mollier chart and table 3,\")\n", + "h1=1371. #B/lb\n", + "h2s=1149. #B/lb\n", + "h3=118. #B/lb\n", + "e=0.9\n", + "print '%s' %(\"Neglecting pump work,\")\n", + "Q1=h1-h3\n", + "W=156. #B/lb\n", + "eta1=W/Q1\n", + "Q=h1-W-h3\n", + "UE=W+e*Q\n", + "fraction = UE/Q1\n", + "#results\n", + "print '%s %.2f' %(\"Fraction supplied = \",fraction)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From mollier chart and table 3,\n", + "Neglecting pump work,\n", + "Fraction supplied = 0.91\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2b - Pg 258" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the fraction of energy supplied \n", + "#initialization of varaibles\n", + "print '%s' %(\"From mollier chart and table 3,\")\n", + "h1=1371. #B/lb\n", + "h2s=1149. #B/lb\n", + "h3=118. #B/lb\n", + "e=0.23\n", + "e2=0.9\n", + "print '%s' %(\"Neglecting pump work,\")\n", + "Q1=h1-h3\n", + "W=156. #B/lb\n", + "eta1=W/Q1\n", + "Q=h1-W-h3\n", + "We=W/e\n", + "UE=We+Q\n", + "UE1=W+e2*Q\n", + "Q2=Q+We\n", + "fraction = UE1/UE\n", + "#results\n", + "print '%s %.2f' %(\"Fraction supplied = \",fraction)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From mollier chart and table 3,\n", + "Neglecting pump work,\n", + "Fraction supplied = 0.64\n" + ] + } + ], + "prompt_number": 3 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_16.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_16.ipynb new file mode 100755 index 00000000..c5f31a4f --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_16.ipynb @@ -0,0 +1,178 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 16 - Gas Cycles" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1a - Pg 279" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the efficiency and work done per pound of fluid\n", + "#initialization of varaibles\n", + "import math\n", + "Pb=75. #psia\n", + "Pc=15. #psia\n", + "k=1.4\n", + "Td=550. #R\n", + "Tb=1700. #R\n", + "cp=0.24\n", + "#calculations\n", + "print '%s' %(\"Gas law solution\")\n", + "Pratio=Pb/Pc\n", + "Ta=Td*math.pow((Pratio),((k-1)/k))\n", + "Tc=Tb/math.pow((Pratio),((k-1)/k))\n", + "Q1=cp*(Tb-Ta)\n", + "Q2=cp*(Tc-Td)\n", + "Wnet=Q1-Q2\n", + "eta=Wnet/Q1\n", + "eta2=1-Td/Ta\n", + "#results\n", + "print '%s %.3f' %(\"Efficiency in 1= \",eta)\n", + "print '%s %.2f' %(\"\\n Efficiency in 2 = \",eta2)\n", + "print '%s %d %s' %(\"\\n Work per pound of fluid =\",Wnet,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Gas law solution\n", + "Efficiency in 1= 0.369\n", + "\n", + " Efficiency in 2 = 0.37\n", + "\n", + " Work per pound of fluid = 73 B/lb\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1b - Pg 280" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the efficiency and work done per pound of fluid\n", + "#initialization of variables\n", + "import math\n", + "Pb=75. #psia\n", + "Pc=15. #psia\n", + "k=1.4\n", + "Td=550. #R\n", + "Tb=1700. #R\n", + "cp=0.24\n", + "#calculations\n", + "Prd=1.4779\n", + "hd=131.46 #B/lb\n", + "Prb=90.95\n", + "hb=422.59 #B/lb\n", + "Pratio=Pb/Pc\n", + "Pra=Pratio*(Prd)\n", + "Ta=868 #R\n", + "ha=208.41\n", + "Prc=Prb/Pratio\n", + "Tc=1113 #R\n", + "hc=269.27\n", + "Q1=hb-ha\n", + "Q2=hc-hd\n", + "Wnet=Q1-Q2\n", + "eta=Wnet/Q1\n", + "#results\n", + "print '%s %.3f' %(\"Efficiency = \",eta)\n", + "print '%s %.2f %s' %(\"\\n Work per pound of fluid =\",Wnet,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Efficiency = 0.357\n", + "\n", + " Work per pound of fluid = 76.37 B/lb\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 282" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the work done and efficiency of the process\n", + "#initialization of varaibles\n", + "e=0.75\n", + "Ta=870. #R\n", + "Tc=1075. #R\n", + "cp=0.24\n", + "Td=550. #R\n", + "Tb=1700. #R\n", + "#calculations\n", + "Tadash=e*(Tc-Ta) +Ta\n", + "Tcdash=Tc+Ta-Tadash\n", + "Q1=cp*(Tb-Tadash)\n", + "Q2=cp*(Tcdash-Td)\n", + "Wnet=Q1-Q2\n", + "eta=Wnet/Q1\n", + "#results\n", + "print '%s %d %s' %(\"Net work done =\",Wnet,\"B/lb\")\n", + "print '%s %.2f' %(\"\\n efficiency = \",eta)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Net work done = 73 B/lb\n", + "\n", + " efficiency = 0.45\n" + ] + } + ], + "prompt_number": 3 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_17.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_17.ipynb new file mode 100755 index 00000000..3b6e3958 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_17.ipynb @@ -0,0 +1,289 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 17 - Fluid Flow : Nozzles and Orifices" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1a - Pg 299" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the exit area\n", + "#initialization of varaibles\n", + "import math\n", + "w=1 #lb/sec\n", + "v2=36.4\n", + "h1=1279.1 #B/lb\n", + "h2=1091.7 #B/lb\n", + "V1=100 #fps\n", + "#calculations\n", + "a2=w*v2/(math.sqrt(2*32.2*778*(h1-h2) + V1*V1)) #sq ft\n", + "a2=1.705 #sq in\n", + "#results\n", + "print '%s %.3f %s' %(\"Exit area =\",a2,\"sq. in\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Exit area = 1.705 sq. in\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1b - Pg 300" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the throat area\n", + "#initialization of varaibles\n", + "import math\n", + "w=1 #lb/sec\n", + "k=1.3\n", + "P=100 #psia\n", + "#calculations\n", + "Pratio=math.pow((2/(k+1)),(k/(k-1)))\n", + "Pt=Pratio*P\n", + "print '%s' %(\"From table 3,\")\n", + "ht=1221.5 #B/lb\n", + "vt=8.841 #cu ft/lb\n", + "at=w*vt/1700.\n", + "#results\n", + "print '%s %.4f %s' %(\"Throat area =\",at,\"sq ft\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 3,\n", + "Throat area = 0.0052 sq ft\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 309" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the Throat and exit areas for both frictionless and frictioned nozzles\n", + "#initialization of varaibles\n", + "import math\n", + "k=1.3\n", + "P=250 #psia\n", + "h0=1263.4 #B/lb\n", + "w=10000. #lb\n", + "cv=0.949\n", + "vts=3.415 #cu ft/lb\n", + "#calculations\n", + "Pratio=math.pow((2/(k+1)),(k/(k-1)))\n", + "Pt=Pratio*P\n", + "hts=1208.2 #B/lb\n", + "h2s=891 #B/lb\n", + "Vts=math.sqrt(2*32.2*778*(h0-hts))\n", + "w=w/3600. #lb/sec\n", + "at=w*vts/(Vts)\n", + "V2=cv*math.sqrt(2*32.2*778*(h0-h2s))\n", + "etan=cv*cv\n", + "h2=928. #B/lb\n", + "print '%s' %(\"From table 3,\")\n", + "v2=276. #cu ft/lb\n", + "a2=w*v2/V2\n", + "a2s=0.17 #ft^2\n", + "Cw=0.98\n", + "at2=at/Cw\n", + "#results\n", + "print '%s %.5f %s' %(\"\\n Throat area =\",at,\"ft^2\")\n", + "print '%s %.3f %s' %(\"\\n Exit area =\",a2,\"ft^2\")\n", + "print '%s %.3f %s' %(\"\\n For frictionless nozzle =\",a2s,\"ft^2\")\n", + "print '%s %.5f %s' %(\"\\n Changed throat area =\",at2,\"ft^2 and exit area is unchanged\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From table 3,\n", + "\n", + " Throat area = 0.00570 ft^2\n", + "\n", + " Exit area = 0.187 ft^2\n", + "\n", + " For frictionless nozzle = 0.170 ft^2\n", + "\n", + " Changed throat area = 0.00582 ft^2 and exit area is unchanged\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3a - Pg 310" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the Exit velocity, throat and exit areas\n", + "#initialization of varaibles\n", + "import math\n", + "w=1 #lb/sec\n", + "Pratio=0.53\n", + "k=1.4\n", + "T0=800. #R\n", + "cp=0.24\n", + "P0=150. #psia\n", + "P2=15. #psia\n", + "#calculations\n", + "Pt=Pratio*P0\n", + "Tratio=math.pow((Pratio),((k-1)/k))\n", + "Tts=T0*Tratio\n", + "Vts=math.sqrt(2*32.2*778*cp*(T0-Tts))\n", + "vts=53.34*Tts/(Pt*144.)\n", + "at=w*vts/(Vts)\n", + "T2s=T0*math.pow((Pt/P0),((k-1)/k))\n", + "T2=460. #R\n", + "V2=math.sqrt(2*32.2*cp*778*(T0-T2))\n", + "v2=53.34*T2/(144.*P2)\n", + "a2=w*v2/V2\n", + "#results\n", + "print '%s %d %s' %(\"Exit velocity =\",Vts,\"fps\")\n", + "print '%s %.5f %s' %(\"\\n Throat area =\",at,\"ft^2\")\n", + "print '%s %.5f %s' %(\"\\n Exit area =\",a2,\"ft^2\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Exit velocity = 1263 fps\n", + "\n", + " Throat area = 0.00246 ft^2\n", + "\n", + " Exit area = 0.00562 ft^2\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3b - Pg 310" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the Exit velocity, throat and exit areas\n", + "#initialization of varaibles\n", + "import math\n", + "h0=191.81 #B/lb\n", + "Pr0=5.526\n", + "w=1 #lb/sec\n", + "Pratio=0.53\n", + "k=1.4\n", + "T0=800. #R\n", + "cp=0.24\n", + "P0=150. #psia\n", + "P2=15. #psia\n", + "Pt=79.5 #psia\n", + "#calculations\n", + "Prt=Pratio*Pr0\n", + "print '%s' %(\"From keenan and kaye steam tables,\")\n", + "Pr=2.929\n", + "Tts=668 #R\n", + "hts=159.9 #B/lb\n", + "Vts=math.sqrt(2*32.2*778*(h0-hts))\n", + "vts=53.34*Tts/(Pt*144.)\n", + "at=w*vts/(Vts)\n", + "Pr2=P2*Pr0/P0\n", + "T2s=415 #R\n", + "h2s=99.13 #B/lb\n", + "h2=110.25 #B/lb\n", + "T2=462 #R\n", + "V2=math.sqrt(2*32.2*778*(h0-h2))\n", + "v2=53.34*T2/(144.*P2)\n", + "a2=w*v2/V2\n", + "#results\n", + "print '%s %d %s' %(\"Exit velocity =\",Vts,\"fps\")\n", + "print '%s %.5f %s' %(\"\\n Throat area =\",at,\"ft^2\")\n", + "print '%s %.5f %s' %(\"\\n Exit area =\",a2,\"ft^2\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From keenan and kaye steam tables,\n", + "Exit velocity = 1264 fps\n", + "\n", + " Throat area = 0.00246 ft^2\n", + "\n", + " Exit area = 0.00564 ft^2\n" + ] + } + ], + "prompt_number": 5 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_18.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_18.ipynb new file mode 100755 index 00000000..7131f2a7 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_18.ipynb @@ -0,0 +1,131 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 18 - Turbines" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 336" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate work done, axial thrust and nozzle bucket efficiency\n", + "#initialization of varaibles\n", + "import math\n", + "drop=50. #B/lb\n", + "cv=0.95\n", + "Vb=700. #fps\n", + "alpha=20.*math.pi/180. #radians\n", + "beta=30.*math.pi/180. #radians\n", + "Cb=0.95\n", + "#calculations\n", + "V1=cv*math.sqrt(2*32.2*778*drop)\n", + "y1=V1*math.cos(alpha)\n", + "z1=V1*math.sin(alpha)\n", + "y1R=y1-Vb\n", + "V1R=math.sqrt(y1R*y1R + z1*z1)\n", + "V2R=Cb*V1R\n", + "y2R=-V2R*math.cos(beta)\n", + "z2=V2R*math.sin(beta)\n", + "Wx=(y1R-y2R)*Vb/32.2\n", + "Fa=(z1-z2)/32.2\n", + "Vc=1582.77\n", + "etanb=Wx/(Vc*Vc /(2*32.2))\n", + "#results\n", + "print '%s %d %s' %(\"Work per pound of fluid =\",Wx,\"ft lbf/lbm\")\n", + "print '%s %.1f %s' %(\"\\n Axial thrust =\",Fa,\"lbf/lbm/sec\")\n", + "print '%s %.2f' %(\"\\n Nozzle bucket efficiency = \",etanb)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Work per pound of fluid = 31221 ft lbf/lbm\n", + "\n", + " Axial thrust = 3.0 lbf/lbm/sec\n", + "\n", + " Nozzle bucket efficiency = 0.80\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 349" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the internal efficiency and the reheat factor\n", + "#initialization of varaibles\n", + "ha=1187.2 #B/lb\n", + "sa=1.6026 #B/lb R\n", + "h3s=895 #B/lb\n", + "h1s=1090 #B/lb\n", + "p1=28 #psia\n", + "h2s=993 #B/lb\n", + "p2=6.2 #psia\n", + "n=0.65\n", + "#calculations\n", + "print '%s' %(\"From Table 3,\")\n", + "h1=ha-n*(ha-h1s)\n", + "s1=1.65 #B/lb R\n", + "h2dash=1024 #B/lb\n", + "h2=h1-n*(h1-h2dash)\n", + "s2=1.706 #B/lb R\n", + "h3dash=953. #B/lb\n", + "h3=h2-n*(h2-h3dash)\n", + "etaT=(ha-h3)/(ha-h3s)\n", + "reheat=etaT/n\n", + "#results\n", + "print '%s %.3f' %(\"Internal efficiency = \",etaT)\n", + "print '%s %.2f' %(\"\\n Reheat factor = \",reheat)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From Table 3,\n", + "Internal efficiency = 0.675\n", + "\n", + " Reheat factor = 1.04\n" + ] + } + ], + "prompt_number": 2 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_19.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_19.ipynb new file mode 100755 index 00000000..7bea3b13 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_19.ipynb @@ -0,0 +1,118 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 19 - Reciprocating Expanders and Compressors" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 370" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the Horsepower output and mean effective pressure\n", + "#initialization of varaibles\n", + "print '%s' %(\"From tables,\")\n", + "h1=1185.3 #B/lb\n", + "v1=4.896 #cu ft/lb\n", + "v2=23.66 #cu ft/lb\n", + "h2=1054.3 #B/lb\n", + "Pd1=1 #cu ft\n", + "Pd2=0.98 #cu ft\n", + "N=300 #rpm\n", + "#calculations\n", + "Wx=h1-h2\n", + "Pd=Pd1+Pd2\n", + "Cl=0.05\n", + "mf=Pd*(1-Cl*(v2/v1 - 1))/v2\n", + "P=Wx*mf*N/(2545./60.)\n", + "mep=P*33000./(N*Pd)\n", + "#results\n", + "print '%s %.3f %s' %(\"Horsepower output =\",P,\"hp\")\n", + "print '%s %d %s' %(\"\\n Mean effective pressure =\",mep,\"psf\")\n", + "#The answers in the book are a bit different due to round off error.\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From tables,\n", + "Horsepower output = 62.679 hp\n", + "\n", + " Mean effective pressure = 3482 psf\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 370" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the mass flow rate when the clearences are 3% and 6%\n", + "#initialization of varaibles\n", + "R=53.34\n", + "T1=540. #R\n", + "P1=15. #psia\n", + "T2=720. #R\n", + "P2=60. #psia\n", + "PD=150. #cu ft/min\n", + "p1=0.03\n", + "p2=0.06\n", + "#calculations\n", + "v1=R*T1/(P1*144.)\n", + "vratio=T1*P2/(P1*T2)\n", + "Nmf=PD*(1-p1*(vratio-1))/v1\n", + "Nmf2=PD*(1-p2*(vratio-1))/v1\n", + "#results\n", + "print '%s %.1f %s' %(\"For clearance of 3 percent, Mass per min =\",Nmf,\"lb/min\")\n", + "print '%s %.1f %s' %(\"\\n For clearance of 6 percent, Mass per min =\",Nmf2,\"lb/min\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For clearance of 3 percent, Mass per min = 10.6 lb/min\n", + "\n", + " For clearance of 6 percent, Mass per min = 9.9 lb/min\n" + ] + } + ], + "prompt_number": 2 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_2.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_2.ipynb new file mode 100755 index 00000000..97405b03 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_2.ipynb @@ -0,0 +1,141 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 2 - Work" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 16" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the net work done\n", + "#Initialization of variables\n", + "from numpy import math\n", + "g=1.4 #Gamma\n", + "P=100. #psia\n", + "V1=3. #cu ft\n", + "Pf=20. #psia\n", + "#calculations\n", + "V2=V1*math.pow((P/Pf),(1/g))\n", + "W=(Pf*V2-P*V1)*144./(1-g)\n", + "#results\n", + "print '%s %d %s' %(\"Net work done =\",W,\"ft lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Net work done = 39810 ft lb\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 17" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the net work done\n", + "#Initialization of variables\n", + "Wb=-33000. #ft-lb\n", + "V2=3. #cu ft\n", + "V1=1. #cu ft\n", + "P=69.4 #psia\n", + "#calculations\n", + "Wa=P*(V2-V1)*144.\n", + "W=Wa+Wb\n", + "#results\n", + "print '%s %d %s' %(\"Net work done =\",W,\"ft-lb\")\n", + "print '%s' %(\"The answer is a bit different due to rounding off error in textbook\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Net work done = -13012 ft-lb\n", + "The answer is a bit different due to rounding off error in textbook\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 20" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the net work done\n", + "import math\n", + "#Initialization of variables\n", + "b=11. #in\n", + "s=15. #in\n", + "l=2.4 #in\n", + "k=80. #psi per in\n", + "#calculations\n", + "a=math.pi*b*b /4\n", + "L=s/12.\n", + "Pm=1.6/l *k\n", + "W=Pm*a*L\n", + "#results\n", + "print '%s %d %s' %(\"Net work done =\",W,\"ft lb\")\n", + "print '%s' %(\"The answer is a bit different due to rounding off error in textbook\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Net work done = 6335 ft lb\n", + "The answer is a bit different due to rounding off error in textbook\n" + ] + } + ], + "prompt_number": 4 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_21.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_21.ipynb new file mode 100755 index 00000000..540caef4 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_21.ipynb @@ -0,0 +1,261 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 21 - Gas Compression" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 387" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the work done and the net heat transferred\n", + "#initialization of varaibles\n", + "import math\n", + "R=53.34\n", + "T1=540 #R\n", + "n=1.4\n", + "g=n\n", + "n2=1.3\n", + "P2=90. #psia\n", + "P1=15. #psia\n", + "cv=0.171\n", + "cp=0.24\n", + "#calculations\n", + "pv=R*T1\n", + "Wk=n*R*T1*(math.pow((P2/P1),((g-1)/g)) -1) /(n-1)\n", + "Wn=n2*R*T1*(math.pow((P2/P1),((n2-1)/n2)) -1) /(n2-1)\n", + "Wt=R*T1*math.log(P2/P1)\n", + "Q=cv*(n-n2)*778*T1*(math.pow((P2/P1),((n2-1)/n2)) -1) /(1-n2)*0.001305\n", + "#0.001305 is the conversion factor\n", + "#results\n", + "print '%s %d %s' %(\"\\n Work in case 1 =\",Wk,\"ft lb/lb\")\n", + "print '%s %d %s' %(\"\\n Work in case 2 =\",Wn,\"ft lb/lb\")\n", + "print '%s %d %s' %(\"\\n Work in case 3 =\",Wt,\"ft lb/lb\")\n", + "print '%s %.1f %s' %(\"\\n Heat transferred =\",Q,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "\n", + " Work in case 1 = 67394 ft lb/lb\n", + "\n", + " Work in case 2 = 63914 ft lb/lb\n", + "\n", + " Work in case 3 = 51609 ft lb/lb\n", + "\n", + " Heat transferred = -16.0 B/lb\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 389" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the adiabatic and isothermal efficiencies\n", + "#initialization of varaibles\n", + "import math\n", + "R=53.34\n", + "T1=540 #R\n", + "n=1.4 #gamma\n", + "g=n\n", + "n2=1.3 #gamma\n", + "P2=90. #psia\n", + "P1=15. #psia\n", + "cv=0.171\n", + "#calculations\n", + "pv=R*T1\n", + "Wk=n*R*T1*(math.pow((P2/P1),((g-1)/g)) -1) /(n-1)\n", + "Wn=n2*R*T1*(math.pow((P2/P1),((n2-1)/n2)) -1) /(n2-1)\n", + "Wt=R*T1*math.log(P2/P1)\n", + "eta1=Wt/Wn\n", + "eta2=Wk/Wn\n", + "#results\n", + "print '%s %.2f' %(\"Adiabatic efficiency = \",eta2)\n", + "print '%s %.2f' %(\"\\n Isothermal efficiency = \",eta1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Adiabatic efficiency = 1.05\n", + "\n", + " Isothermal efficiency = 0.81\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 389" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the heat transferred\n", + "#initialization of varaibles\n", + "import math\n", + "R=53.34\n", + "T1=540 #R\n", + "n=1.4 #gamma\n", + "g=n\n", + "n2=1.3 #gamma\n", + "P2=90 #psia\n", + "P1=15 #psia\n", + "cv=0.171\n", + "eta=0.95\n", + "cp=0.24\n", + "#calculations\n", + "pv=R*T1\n", + "Wk=n*R*T1*(math.pow((P2/P1),((g-1)/g)) -1) /(n-1)\n", + "Wn=n2*R*T1*(math.pow((P2/P1),((n2-1)/n2)) -1) /(n2-1)\n", + "Wt=R*T1*math.log(P2/P1)\n", + "Wx=-Wk/eta\n", + "dh=cp*T1*(1.52 - 1)\n", + "Q=dh+Wx/778.\n", + "#results\n", + "print '%s %.1f %s' %(\"Heat transferred =\",Q,\"B/lb\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Heat transferred = -23.8 B/lb\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4 - Pg 395" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the volumetric efficiency\n", + "#initialization of varaibles\n", + "import math\n", + "n=1.3\n", + "P1=15. #psia\n", + "P2=75. #psia\n", + "eta=0.5\n", + "eta2=0\n", + "#calculations\n", + "Pr=math.pow((P2/P1),(1/n))\n", + "Cl=(1-eta)/(Pr-1)\n", + "Cl2=(1-eta2)/(Pr-1)\n", + "#results\n", + "print '%s %.3f' %(\"For volumetric efficiency to be 0.5, Clearance = \",Cl)\n", + "print '%s %.3f' %(\"\\n For volumetric efficiency to be 0, Clearance = \",Cl2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For volumetric efficiency to be 0.5, Clearance = 0.204\n", + "\n", + " For volumetric efficiency to be 0, Clearance = 0.408\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5 - Pg 398" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the single stage and two stage efficiencies\n", + "#initialization of varaibles\n", + "import math\n", + "P1=5 #psia\n", + "P2=83.5 #psia\n", + "n=1.25\n", + "per=0.03\n", + "#calculations\n", + "nv1=1- per*(math.pow((P2/P1),(1/n)) -1)\n", + "nv2=1-per*(math.pow((math.sqrt(P2/P1)),(1/n)) -1)\n", + "#results\n", + "print '%s %.3f' %(\"For single stage machine = \",nv1)\n", + "print '%s %.3f' %(\"\\n For Two stage machine = \",nv2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For single stage machine = 0.745\n", + "\n", + " For Two stage machine = 0.937\n" + ] + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_22.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_22.ipynb new file mode 100755 index 00000000..4a1a7d0c --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_22.ipynb @@ -0,0 +1,256 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 22 - Combustion Processes : First law analysis" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5 - Pg 419" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the heat transferred from the system\n", + "#initialization of varaibles\n", + "mO=1.33\n", + "CO=0.155\n", + "mC=3.67\n", + "CC=0.165\n", + "t2=1000. #F\n", + "tb=68. #F\n", + "t1=300. #F\n", + "mC2=1.\n", + "CC2=0.17\n", + "mO2=4.\n", + "CO2=0.155\n", + "H=-14087 #B/lb\n", + "#calculations\n", + "dE2=mO*CO*(t2-tb) + mC*CC*(t2-tb)\n", + "dE1=mO2*CO2*(tb-t1) + mC2*CC2*(tb-t1)\n", + "Q=dE2+dE1+H\n", + "#results\n", + "print '%s %d %s' %(\"Heat transfer from the system =\",Q,\"Btu\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Heat transfer from the system = -13513 Btu\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6 - Pg 422" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the cp value of methane\n", + "#initialization of varaibles\n", + "H1=17889 #Cal/g\n", + "H2=-94052 #Cal/g\n", + "H3=2* -68317 #Cal/g\n", + "#calculations\n", + "x=H1+H2+H3\n", + "#results\n", + "print '%s %d %s' %(\"Constant pressure heating value of methane =\",x,\" cal/gm formula wt.\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Constant pressure heating value of methane = -212797 cal/gm formula wt.\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7 - Pg 423" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the cv value of methane\n", + "#initialization of varaibles\n", + "HV=4344 #B/lb\n", + "xC=56 #lb\n", + "R=1.986\n", + "T=530 #R\n", + "MC=56 #g/mol\n", + "#calculations\n", + "HR=xC*HV\n", + "Eb=-HR -R*T*(2-3)\n", + "HV=-Eb/MC\n", + "#results\n", + "print '%s %d %s' %(\"COnstant volume heating value =\",HV,\"B/lb \")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "COnstant volume heating value = 4325 B/lb \n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8 - Pg 426" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the products temperature\n", + "#initialization of varaibles\n", + "dH2=14087 #B/lb\n", + "xc=3.67 #lb\n", + "xN=8.78 #lb\n", + "tb=100 #F\n", + "#calculations\n", + "dt2=dH2/(xc*0.196 + xN*0.248)\n", + "t2=dt2+tb\n", + "#results\n", + "print '%s %d %s' %(\"products temperature =\",t2,\"F\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "products temperature = 4963 F\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9 - Pg 430" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the efficiency\n", + "#initialization of varaibles\n", + "Heat=14087 #Btu/lb\n", + "x1=0.9 #lb\n", + "x2=0.05 #lb\n", + "x3=0.05 #lb\n", + "Heat2=3952 #Btu/lb\n", + "#calculations\n", + "h1=x1*Heat\n", + "h2=x2*Heat2\n", + "e=(h1+h2)/Heat\n", + "#results\n", + "print '%s %.2f' %(\"Efficiency = \",e)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Efficiency = 0.91\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10 - Pg 431" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the efficiency of the process\n", + "#initialization of varaibles\n", + "print '%s' %(\"From data and steam tables,\")\n", + "Q=10240000. #B/hr\n", + "w=700. #lb/hr\n", + "h=19500. #B/lb\n", + "#calculations\n", + "HV=w*h\n", + "e=Q/HV\n", + "#results\n", + "print '%s %.2f' %(\"Efficiency = \",e)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From data and steam tables,\n", + "Efficiency = 0.75\n" + ] + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_23.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_23.ipynb new file mode 100755 index 00000000..1ecba1a2 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_23.ipynb @@ -0,0 +1,302 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 23 - Internal Combustion" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 437" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the efficiency, mean effective pressure and work per machine cycle\n", + "#initialization of varaibles\n", + "print '%s' %(\"from chart\")\n", + "T6=2600 #R\n", + "mratio=0.05\n", + "V6d=82 #cu ft\n", + "E6d=465 #Btu\n", + "H6d=655 #Btu\n", + "T6d=2480 #R\n", + "Hs=58 #Btu\n", + "LHV=19256\n", + "#calculations\n", + "H1=mratio*H6d + (1-mratio)*Hs\n", + "dV=22-3.67\n", + "PD=0.12\n", + "Work=446*PD/dV\n", + "pm=Work*778/(144.*PD)\n", + "eta=446./((1-mratio)*(LHV*0.0665))\n", + "#results\n", + "print '%s %.3f' %(\"Efficiency = \",eta)\n", + "print '%s %d %s' %(\"\\n Mean effective pressure =\",pm,\"psi\")\n", + "print '%s %.2f %s' %(\"\\n Work per machine cycle =\",Work,\"Btu\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "from chart\n", + "Efficiency = 0.367\n", + "\n", + " Mean effective pressure = 131 psi\n", + "\n", + " Work per machine cycle = 2.92 Btu\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 441" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the efficiency of the machine\n", + "#initialization of varaibles\n", + "f=0.03\n", + "T6=1500. #R\n", + "print '%s' %(\"from air tables,\")\n", + "hi=131.46 #B/lb\n", + "h6=381 #B/lb\n", + "vratio=1/15.\n", + "v1r=120.7\n", + "P1=15. #psi\n", + "T1=580. #R\n", + "x=0.5\n", + "Tb=520. #R\n", + "H=18500 #B/lb\n", + "mh=0.0345\n", + "m3=1.065\n", + "#calculations\n", + "h1=f*h6+(1-f)*hi\n", + "v2r=v1r*vratio\n", + "T2=1615 #R\n", + "u2=289.05 #B/lb\n", + "P2=T2*1/vratio *P1/T1\n", + "theo=0.069 #lb/lb of air\n", + "m=theo*x\n", + "h3B=0.242*Tb\n", + "m3=1+0.03+0.0345\n", + "h3=(638+284)/1.065 +h3B\n", + "T3=3520 #R\n", + "P3=626 #psi\n", + "v3=53.34*T3/(P3*144)\n", + "v3p=v3*m3\n", + "v1=53.35*T1/(144*P1)\n", + "v2=14.7/P1\n", + "m1=1.03\n", + "h3=992.\n", + "h4=531.\n", + "T3=3520. #R\n", + "T4=2030. #R\n", + "W12=m1*(98.9-289.05)\n", + "W23=P2*(v3p-v2)*144/778.\n", + "W34=m3*(h3-h4-53.4*(T3-T4)/778.)\n", + "W=W12+W23+W34\n", + "eta=W/(mh*H)\n", + "#results\n", + "print '%s %.3f' %(\"Efficiency = \",eta)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "from air tables,\n", + "Efficiency = 0.516\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 448" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the efficiency, Air rate and back work ratio of the process\n", + "#initialization of varaibles\n", + "print '%s' %(\"Using air tables,\")\n", + "h1=124.27\n", + "pr1=1.2147\n", + "p2byp1=6.\n", + "p1=15.\n", + "p4=15.\n", + "eta=0.8\n", + "#calculations\n", + "pr2=p2byp1*pr1\n", + "h2s=197.5\n", + "h2=h1+(h2s-h1)/eta\n", + "h2B=124.3\n", + "dhB=-18500 #B/lb\n", + "dh2=h2B-h2\n", + "T3=1910 #R\n", + "h3=479.85\n", + "pr3=144.53\n", + "h3B=h2B\n", + "dh3=h3-h3B\n", + "wratio=(-dh3-dh2)/(dh3+dhB)\n", + "pr4=28.91\n", + "h4s=306.9\n", + "h4=h3-eta*(h3-h4s)\n", + "Wt=(1+wratio)*(h3-h4)\n", + "Wc=(h2-h2B)\n", + "Wnet=Wt-Wc\n", + "E=Wnet/(wratio*-dhB)\n", + "rate=2545./Wnet\n", + "BWratio=Wc/Wnet\n", + "#results\n", + "print '%s %.3f' %(\"Efficiency = \",E)\n", + "print '%s %.1f %s' %(\"\\n Air rate =\",rate,\"lb air/hp hr\")\n", + "print '%s %.2f' %(\"\\n Back work ratio = \",BWratio)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Using air tables,\n", + "Efficiency = 0.182\n", + "\n", + " Air rate = 52.1 lb air/hp hr\n", + "\n", + " Back work ratio = 1.87\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4 - Pg 451" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the specific impulse, thrust and efficiency of the process\n", + "#initialization of varaibles\n", + "import math\n", + "V1=587. #fps\n", + "etaD=0.9\n", + "etaC=0.8\n", + "h1=114.69\n", + "P1=10. #psia\n", + "P6=P1\n", + "P3=626 #psia\n", + "dhB=-19100. #B/lb\n", + "T1=480. #R\n", + "#calculations\n", + "h2s=etaD*V1*V1 /(778.*2*32.2) +h1\n", + "print '%s' %(\"From tables,\")\n", + "Pr2s=1.104\n", + "Pr1=0.9182\n", + "P2=P1*Pr2s/Pr1\n", + "h2=h1+(h2s-h1)/etaD\n", + "T2=509. #R\n", + "Pr2=1.127\n", + "Pr3s=Pr2*P3/P2\n", + "Pr3s=6.76\n", + "h3s=203.2\n", + "h3=(h3s-h2)/etaC +h2\n", + "T3=930. #R\n", + "P3=6*P2\n", + "T4=2160. #R\n", + "h4=549.35\n", + "Pr4=238.\n", + "h4B=126.66\n", + "dh4=422.7\n", + "h3B=h4B\n", + "dh3=h3-h3B\n", + "cp=0.5\n", + "Ta=480. #R\n", + "Tb=530. #R\n", + "dhf=cp*(Tb-Ta)\n", + "wratio=(-dh4+dh3)/(dh4+dhf+dhB)\n", + "h5s=425.3\n", + "Pr5s=93.1\n", + "P5=27.6\n", + "T5=1801 #R\n", + "Pr5=114.28\n", + "Pr6s=Pr5*P6/P5\n", + "h5=450.\n", + "h6=351.\n", + "V6=math.sqrt(2*32.2*778*(h5-h6))\n", + "SI=((1+wratio)*V6 -V1)/(32.2)\n", + "v1=53.34*T1/(P1*144.)\n", + "wa=V1/v1\n", + "thrust = wa*SI\n", + "SC=wa*0.0174*3600/1840.\n", + "eff=2545/(SC*-dhB)\n", + "#results\n", + "print '%s %.1f %s' %(\"Specific impulse =\",SI,\"lb/lb per sec of air\")\n", + "print '%s %d %s' %(\"\\n Thrust =\",thrust,\"lb\")\n", + "print '%s %.3f' %(\"\\n Efficiency = \",eff)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From tables,\n", + "Specific impulse = 52.1 lb/lb per sec of air\n", + "\n", + " Thrust = 1721 lb\n", + "\n", + " Efficiency = 0.119\n" + ] + } + ], + "prompt_number": 4 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_24.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_24.ipynb new file mode 100755 index 00000000..bc3bc67c --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_24.ipynb @@ -0,0 +1,76 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 24 - Refrigeration" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 461" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the cop in both cases and also the piston displacement\n", + "#initialization of varaibles\n", + "print '%s' %(\"From tables,\")\n", + "h1=611.8 #B/lb\n", + "h2=704.4 #B/lb\n", + "h3=127.4 #B/lb\n", + "h4=h3\n", + "T2=460. #R\n", + "T1=76+460. #R\n", + "W=10000. #B/hr\n", + "e=0.7\n", + "#calculations\n", + "Qe=h1-h4\n", + "Wc=h2-h1\n", + "CP=Qe/Wc\n", + "CP2=T2/(T1-T2)\n", + "w=W/(Qe*60.)\n", + "v1=9.116 #cu ft/lb\n", + "PD=w*v1/(e)\n", + "#results\n", + "print '%s %.2f' %(\"Coefficient of performance in case 1 = \",CP)\n", + "print '%s %.2f' %(\"\\n Coefficient of performance in case 2 = \",CP2)\n", + "print '%s %.2f %s' %(\"\\n Piston displacement =\",PD,\"cu ft/min\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From tables,\n", + "Coefficient of performance in case 1 = 5.23\n", + "\n", + " Coefficient of performance in case 2 = 6.05\n", + "\n", + " Piston displacement = 4.48 cu ft/min\n" + ] + } + ], + "prompt_number": 1 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_25.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_25.ipynb new file mode 100755 index 00000000..8b36d41a --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_25.ipynb @@ -0,0 +1,294 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 25 - Air - Water Vapor Mixtures" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 482" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the specific humidity and dew temperature\n", + "#initialization of varaibles\n", + "Pg=0.4298 #steam tables psi\n", + "phi=0.5\n", + "P=14.7 #psi\n", + "#calculations\n", + "pw=phi*Pg\n", + "Pa=P-pw\n", + "gamma=0.622*pw/Pa\n", + "T=55 #F from dew point tables \n", + "#results\n", + "print '%s %.5f %s' %(\"Specific humidity =\",gamma,\"lb water/lb dry air\")\n", + "print '%s %d %s' %(\"\\n Dew temperature =\",T,\"F\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Specific humidity = 0.00923 lb water/lb dry air\n", + "\n", + " Dew temperature = 55 F\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 486" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the enthalpy, relative and specific humidities, Dew point temperature of the liquid\n", + "#initialization of varaibles\n", + "print '%s' %(\"From psychrometric chart ,\")\n", + "hgdp=1061.8 \n", + "cpw=0.44\n", + "tdb=72 #F\n", + "cp=0.24\n", + "g=0.0071\n", + "#calculations\n", + "rh=0.42\n", + "sp=g\n", + "tdp=58 #F\n", + "hw=hgdp+cpw*tdb\n", + "h=cp*tdb+g*hw\n", + "#results\n", + "print '%s %.2f %s' %(\"Enthalpy =\",h,\" B/lb dry air\")\n", + "print '%s %.2f' %(\"\\n relative humidity = \",rh)\n", + "print '%s %.2f' %(\"\\n specific humidity = \",sp)\n", + "print '%s %d %s' %(\"\\n Dew point temperature =\",tdp,\"F\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From psychrometric chart ,\n", + "Enthalpy = 25.04 B/lb dry air\n", + "\n", + " relative humidity = 0.42\n", + "\n", + " specific humidity = 0.01\n", + "\n", + " Dew point temperature = 58 F\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 488" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the water and heat to be supplied\n", + "#initialization of varaibles\n", + "print '%s' %(\"From the psychrometric chart,\")\n", + "ha=12.9 #B/lb\n", + "g1=0.0032 #lb water/ lb dry air\n", + "g2=0.0078 #lb water/ lb dry air\n", + "hl=13 #B/lb\n", + "hd=25.33 #B/lb\n", + "p=14.7 #psia\n", + "phi=0.6\n", + "cp=0.24\n", + "t2=70 #F\n", + "#calculations\n", + "wl=g2-g1\n", + "Q=hd-ha-wl*hl\n", + "pg=0.1217 #psia\n", + "pa=p-pg\n", + "G1=0.622*pg*phi/pa\n", + "G2=0.00788\n", + "wl2=G2-G1\n", + "t1=40 #F\n", + "hw1=1061.8 + 0.44*t1\n", + "hw2=1092.6 #B/lb\n", + "Q2=cp*(t2-t1) + G2*hw2 -G1*hw1 - wl2*hl\n", + "#results\n", + "print '%s' %(\"Method 1\")\n", + "print '%s %.4f %s' %(\"\\n Water to be supplied =\",wl,\"lb/lb of dry air\")\n", + "print '%s %.1f %s' %(\"\\n heat supplied =\",Q,\"B/lb of dry air\")\n", + "print '%s' %(\"\\n Method 2\")\n", + "print '%s %.5f %s' %(\"\\n Water to be supplied =\",wl2,\"lb/lb of dry air\")\n", + "print '%s %.1f %s' %(\"\\n heat supplied =\",Q2,\"B/lb of dry air\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From the psychrometric chart,\n", + "Method 1\n", + "\n", + " Water to be supplied = 0.0046 lb/lb of dry air\n", + "\n", + " heat supplied = 12.4 B/lb of dry air\n", + "\n", + " Method 2\n", + "\n", + " Water to be supplied = 0.00476 lb/lb of dry air\n", + "\n", + " heat supplied = 12.4 B/lb of dry air\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4 - Pg 489" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the air supplied, temperature and humidity of air\n", + "#initialization of varaibles\n", + "print '%s' %(\"From psychrometric charts,\")\n", + "e=0.7\n", + "phi=0.5\n", + "g1=0.0131 #lb water/lb dry air\n", + "h1=32.36 #B/lb of dry air\n", + "g3=0.0073\n", + "h3=24.26 #B/lb\n", + "pg=0.3390 #psia\n", + "T3=528 #R\n", + "V3=1000\n", + "Rw=85.8\n", + "#calculations\n", + "pw3=phi*pg\n", + "ww3=pw3*144*V3/(Rw*T3)\n", + "wa3=ww3/g3\n", + "wa1=phi*wa3\n", + "wa2=phi*wa3\n", + "ww1=g1*wa1\n", + "ww2=ww3-ww1\n", + "g2=ww2/wa2\n", + "h2=(wa3*h3-wa1*h1)/wa2\n", + "tdb=61 #F\n", + "#results\n", + "print '%s %.3f %s' %(\"Air supplied =\",ww2,\"lb/min\")\n", + "print '%s %d %s' %(\"\\n temperature =\",tdb,\"F\")\n", + "print '%s %.5f %s' %(\"\\n Humidity =\",g2,\"lb water/lb dry air\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From psychrometric charts,\n", + "Air supplied = 0.055 lb/min\n", + "\n", + " temperature = 61 F\n", + "\n", + " Humidity = 0.00150 lb water/lb dry air\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5 - Pg 493" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the cooling temperature, heat transfer and fraction of heat removed\n", + "#initialization of varaibles\n", + "print '%s' %(\"From psychrometric charts,\")\n", + "g1=0.0131 #lb water/lb dry air\n", + "g2=0.0093 #lb water/lb dry air\n", + "h1=32.36 #B/lb dry air\n", + "h2=27.03\n", + "hd2=23.4 #B/lb dry air\n", + "hf=23.4 #B/lb dry air\n", + "hw1=1094.5\n", + "#calculations\n", + "tdp=55.3 #F\n", + "wratio=g1-g2\n", + "Qc=hd2-h1+wratio*hf\n", + "Qh=h2-hd2\n", + "Heat=wratio*(hw1-hf)\n", + "frac=-Heat/Qc\n", + "#results\n", + "print '%s %.1f %s' %(\"Cooling temperature = \",tdp,\"F\")\n", + "print '%s %.2f %s' %(\"\\n heat transfer = \",Heat,\"B/lb dry air\")\n", + "print '%s %.2f' %(\"\\n Fraction of heat removed =\",frac)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "From psychrometric charts,\n", + "Cooling temperature = 55.3 F\n", + "\n", + " heat transfer = 4.07 B/lb dry air\n", + "\n", + " Fraction of heat removed = 0.46\n" + ] + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_3.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_3.ipynb new file mode 100755 index 00000000..f5793d78 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_3.ipynb @@ -0,0 +1,240 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 3 - Temperature and Heat" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 33" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the mass of water required per pound of iron \n", + "import math\n", + "#Initialization of variables\n", + "T1=500 #F\n", + "T2=100 #F\n", + "Tf=75 #F\n", + "cpi=0.120 #B/lb F\n", + "cpw=1.0 #B/lb F\n", + "#calculations\n", + "Qw=1*cpw*(T2-Tf)\n", + "Qi=-1*cpi*(T2-T1)\n", + "mw=Qi/Qw\n", + "#results\n", + "print '%s %.2f %s' %(\"Mass of water = \",mw,\"lb water/lb iron\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Mass of water = 1.92 lb water/lb iron\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 34" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the net heat transferred\n", + "import math\n", + "import scipy\n", + "from scipy import integrate\n", + "#Initialization of variables\n", + "m=5 #lb\n", + "T1=1540+460 #R\n", + "T2=540+460 #R\n", + "#calculations\n", + "def q(T):\n", + "\tcp=m*(0.248+0.448*math.pow(10,-8) *T*T)\n", + "\treturn cp;\n", + "\n", + "Q=scipy.integrate.quad(q,T1,T2)\n", + "#results\n", + "print '%s %d %s' %(\"Heat transferred =\",Q[0],\"Btu\")\n", + "print '%s' %(\"The answer is a bit different due to rounding off error in textbook\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Heat transferred = -1292 Btu\n", + "The answer is a bit different due to rounding off error in textbook\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 36" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the heat required for the process\n", + "#Initialization of variables\n", + "Tm=235 #F\n", + "Tb=832 #F\n", + "T=70 #F\n", + "cps=0.18 #B/lb F\n", + "cpl=0.235 #B/lb F\n", + "Lf=15.8 #B/lb\n", + "Lv=120 #B/lb\n", + "m=10 #lb\n", + "#calculations\n", + "Qa=m*cps*(Tm-T)\n", + "Qb=m*Lf\n", + "Qc=m*cpl*(Tb-Tm)\n", + "Qd=m*Lv\n", + "Q=Qa+Qb+Qc+Qd\n", + "#results\n", + "print '%s %.1f %s' %(\"Heat required =\",Q,\"Btu\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Heat required = 3057.9 Btu\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4 - Pg 36" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the mass of ice required for the process\n", + "#Initialization of variables\n", + "T1=22 #F\n", + "T2=32 #F\n", + "T3=40 #F\n", + "T4=70 #F\n", + "cps=0.501 #B/lb F\n", + "cpw=1 #B/lb F\n", + "Lf=143.3 #B/lb\n", + "m=40 #lb\n", + "#calculations\n", + "Qa=cps*(T2-T1)\n", + "Qb=Lf\n", + "Qc=cpw*(T3-T2)\n", + "Qd=m*cpw*(T3-T4)\n", + "mi=-Qd/(Qa+Qb+Qc)\n", + "#results\n", + "print '%s %.2f %s' %(\"Mass of ice required =\",mi,\"lb ice\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Mass of ice required = 7.68 lb ice\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5 - Pg 37" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the extra mass of ice required\n", + "#Initialization of variables\n", + "T1=22 #F\n", + "T2=32 #F\n", + "T3=40 #F\n", + "T4=70 #F\n", + "cps=0.501 #B/lb F\n", + "cpw=1 #B/lb F\n", + "Lf=143.3 #B/lb\n", + "m=40 #lb\n", + "cp=0.092\n", + "mc=10\n", + "#calculations\n", + "Qa=cps*(T2-T1)\n", + "Qb=Lf\n", + "Qc=cpw*(T3-T2)\n", + "Qe=mc*cp*(T3-T4)\n", + "mi=-Qe/(Qa+Qb+Qc)\n", + "#results\n", + "print '%s %.3f %s' %(\"Extra Mass of ice required =\",mi,\"lb ice\")" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Extra Mass of ice required = 0.177 lb ice\n" + ] + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_5.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_5.ipynb new file mode 100755 index 00000000..ef287258 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_5.ipynb @@ -0,0 +1,225 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 5 - First Law of Thermodynamics" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 59" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the change in internal energy\n", + "#initialization of varaibles\n", + "V1=10 #cu ft\n", + "P1=15 #psia\n", + "V2=5 #cu ft\n", + "H=34.7 #Btu\n", + "#calculations\n", + "W=P1*(V2-V1)*144\n", + "dE=-H-W/778.\n", + "#results\n", + "print '%s %.1f %s' %(\"Internal energy change =\",dE,\"Btu\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Internal energy change = -20.8 Btu\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 59" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the work done and the change in Internal energy\n", + "#initialization of varaibles\n", + "dT=35 #F\n", + "H=34 #Btu\n", + "cv=1.2 #B/lb F\n", + "m= 2 #lb\n", + "#calculations\n", + "U=cv*dT*m\n", + "W=H-U\n", + "#results\n", + "print '%s %d %s' %(\"Work done =\",W,\"Btu\")\n", + "print '%s %.1f %s' %(\"\\n Internal energy change =\",U,\"Btu\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Work done = -50 Btu\n", + "\n", + " Internal energy change = 84.0 Btu\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 60" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the change in internal energy in both cases\n", + "#initialization of varaibles\n", + "p=500 #psia\n", + "V2=0.9278 #cu ft\n", + "V1=0.0197 #cu ft\n", + "h2=1204.4 #B/lb\n", + "h1=449.4 #B/lb\n", + "#calculations\n", + "W=p*(V2-V1)*144.\n", + "du=h2-h1-144*p*(V2-V1)/778.\n", + "du2=h2-h1-W/778.\n", + "#results\n", + "print '%s %.1f %s' %(\"Change in internal energy = \",du,\"Btu\")\n", + "print '%s %.1f %s' %(\"\\n Method 2, Internal energy change = \",du2,\"Btu\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Change in internal energy = 671.0 Btu\n", + "\n", + " Method 2, Internal energy change = 671.0 Btu\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4a - Pg 61" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the heat liberated\n", + "#initialization of varaibles\n", + "import math\n", + "P1=75. #psia\n", + "P2=15. #psia\n", + "V1=6. #cu ft\n", + "g=1.2 #gamma\n", + "m=3. #lb\n", + "#calculations\n", + "V2=V1*math.pow((P1/P2),(1/g))\n", + "U=0.48*(P2*V2-P1*V1)\n", + "W=(P2*V2-P1*V1)*144./((1-g)*778.)\n", + "Q=U+W\n", + "#results\n", + "print '%s %.3f %s' %(\"Heat = \",Q,\"Btu\")\n", + "#The answer given in textbook is wrong. please check using a calculator\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Heat = 47.162 Btu\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4b - Pg 62" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the net work done\n", + "#initialization of varaibles\n", + "import math\n", + "P1=75. #psia\n", + "P2=15. #psia\n", + "V1=6. #cu ft\n", + "g=1.2 #gamma\n", + "m=3 #lb\n", + "#calculations\n", + "Q=30 #Btu\n", + "V2=V1*math.pow((P1/P2),(1/g))\n", + "U=0.48*(P2*V2-P1*V1)\n", + "W=Q-U\n", + "#results\n", + "print '%s %.1f %s' %(\"Work done =\",W,\" Btu\")\n", + "#The answer given in textbook is wrong. please check using a calculator\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Work done = 80.8 Btu\n" + ] + } + ], + "prompt_number": 4 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_6.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_6.ipynb new file mode 100755 index 00000000..49800538 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_6.ipynb @@ -0,0 +1,154 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 6 - Flow processes : First law analysis" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 76" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the power output of the machine\n", + "#initialization of varaibles\n", + "import math\n", + "u1=1111.9 #Btu/lb\n", + "P1=170. #psia\n", + "v1=2.675 #cu ft/lb\n", + "v2=100.9 #cu ft/lb\n", + "z1=10. #ft\n", + "V1=6000./60. #ft/sec\n", + "Q=-1000.\n", + "u2=914.6 #B/lb\n", + "P2=3. #psia\n", + "V2=300. #ft/sec\n", + "rate=2500. #lb/hr\n", + "#calculations\n", + "Wx=rate*(u1-u2 + (P1*v1-P2*v2)*144/778 +(V1*V1 -V2*V2)/(2*32.2*778.) +z1/778. +Q/rate)\n", + "f=3.92*math.pow(10,-4)\n", + "power = Wx*f\n", + "#results\n", + "print '%s %d %s' %(\"Power output of turbine = \",Wx,\"B/hr\")\n", + "print '%s %.1f %s' %(\"\\n Power output in hp =\",power,\"hp\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Power output of turbine = 558647 B/hr\n", + "\n", + " Power output in hp = 219.0 hp\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 78" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the flow rate\n", + "#initialization of varaibles\n", + "w1=500. #lb/min\n", + "h1=132.9 #lb/min\n", + "h2=1150. #B/lb\n", + "h3=180. #B/lb\n", + "#calculations\n", + "w2=(w1*h1-w1*h3)/(h3-h2)\n", + "#results\n", + "print '%s %.1f %s' %(\"Flow rate =\",w2,\"lb/min\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Flow rate = 24.3 lb/min\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 79" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the average velocity and rate of flow\n", + "#initialization of varaibles\n", + "import math\n", + "v2=5.434 #cu ft/lb\n", + "v1=4.937 #cu ft/lb\n", + "h1=1227.6\n", + "h2=1223.9\n", + "A1=math.pi/144\n", + "#calculations\n", + "Vratio=v2/v1\n", + "V1=math.sqrt(64.4*(h1-h2)*778./(Vratio*Vratio -1))\n", + "V2=V1*Vratio\n", + "w=A1*V1/v1\n", + "#results\n", + "print '%s %d %s' %(\"Average velocity at 1 =\",V1,\"fps\")\n", + "print '%s %d %s' %(\"\\n Average velocity at 2 =\",V2,\"fps\")\n", + "print '%s %.2f %s' %(\"\\n Rate of flow = \",w,\"lb/sec\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Average velocity at 1 = 936 fps\n", + "\n", + " Average velocity at 2 = 1030 fps\n", + "\n", + " Rate of flow = 4.14 lb/sec\n" + ] + } + ], + "prompt_number": 3 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_8.ipynb b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_8.ipynb new file mode 100755 index 00000000..fe970912 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/Chapter_8.ipynb @@ -0,0 +1,244 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 8 - Basic Applications of the Second Law" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1 - Pg 130" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the Max Efficiency possible\n", + "#initialization of varaibles\n", + "T1=85+460. #R\n", + "T2=50+460. #R\n", + "#calculations\n", + "eta=(T1-T2)/T1\n", + "#results\n", + "print '%s %.1f %s' %(\"Max. efficiency =\",eta*100,\"percent\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Max. efficiency = 6.4 percent\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2 - Pg 130" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the Max Efficiency possible\n", + "#initialization of varaibles\n", + "T1=1050+460. #R\n", + "T2=90+460. #R\n", + "#calculations\n", + "eta=(T1-T2)/T1\n", + "#results\n", + "print '%s %d %s' %(\"Max. possible efficiency = \",eta*100,\"percent\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Max. possible efficiency = 63 percent\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3 - Pg 130" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the change in entropy\n", + "#initialization of varaibles\n", + "import math\n", + "T1=50+460. \t#R\n", + "T2=150+460. #R\n", + "m=1.\n", + "cp=0.240\n", + "#calculations\n", + "ds=m*cp*(math.log(T2) - math.log(T1))\n", + "#results\n", + "print '%s %.4f %s' %(\"Change in entropy =\",ds,\"B/ F abs\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Change in entropy = 0.0430 B/ F abs\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4 - Pg 131" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate the change in entropy\n", + "#initialization of varaibles\n", + "import math\n", + "T1=50+460. #R\n", + "T2=150+460. #R\n", + "m=1\n", + "cp=0.240\n", + "#calculations\n", + "ds=m*cp*(math.log(T2) - math.log(T1))\n", + "#results\n", + "print '%s %.4f %s' %(\"Change in entropy = \",ds,\"B/ F abs\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Change in entropy = 0.0430 B/ F abs\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5 - Pg 131" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the total change in entropy\n", + "#initialization of varaibles\n", + "import math\n", + "Q=826. #B/lb\n", + "T=860. #R\n", + "T1=2000+460. #R\n", + "T2=1000+460. #R\n", + "#calculations\n", + "ds=Q/T\n", + "dsgas=Q*(math.log(T2)-math.log(T1))/(T1-T2)\n", + "dst=ds+dsgas\n", + "#results\n", + "print '%s %.3f %s' %(\"Total entropy change =\",dst,\"B/R\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Total entropy change = 0.530 B/R\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6 - Pg 132" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#Calculate the loss in energy and increase in unavailable energy\n", + "#initialization of varaibles\n", + "T0=540. #R\n", + "Q=826. #B/lb\n", + "ds=0.534\n", + "ds2=0.431\n", + "#calculations\n", + "tds=T0*ds\n", + "tds2=T0*ds2\n", + "H=Q-tds2\n", + "Loss=tds/H\n", + "#results\n", + "print '%s %.1f %s' %(\"Increase in unavailable energy = \",tds,\"Btu\")\n", + "print '%s %.1f %s' %(\"Loss =\",Loss*100,\"percent\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Increase in unavailable energy = 288.4 Btu\n", + "Loss = 48.6 percent\n" + ] + } + ], + "prompt_number": 4 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/README.txt b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/README.txt new file mode 100755 index 00000000..2046585c --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/README.txt @@ -0,0 +1,10 @@ +Contributed By: Chaitanya Potti +Course: btech +College/Institute/Organization: IITB +Department/Designation: Chemical engineering +Book Title: Mechanical Engineering Thermodynamics +Author: Mooney D A +Publisher: Prentice Hall, London +Year of publication: 1953 +Isbn: 1466511796 +Edition: 1 \ No newline at end of file diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap23.png b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap23.png new file mode 100755 index 00000000..0eed9ce1 Binary files /dev/null and b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap23.png differ diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap24.png b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap24.png new file mode 100755 index 00000000..4ff3b58c Binary files /dev/null and b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap24.png differ diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap25.png b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap25.png new file mode 100755 index 00000000..850e306c Binary files /dev/null and b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap25.png differ -- cgit