diff options
| author | kinitrupti | 2017-05-12 18:53:46 +0530 |
|---|---|---|
| committer | kinitrupti | 2017-05-12 18:53:46 +0530 |
| commit | 6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d (patch) | |
| tree | 22789c9dbe468dae6697dcd12d8e97de4bcf94a2 /Mechanical_Engineering_Thermodynamics_by_Mooney_D_A | |
| parent | d36fc3b8f88cc3108ffff6151e376b619b9abb01 (diff) | |
| download | Python-Textbook-Companions-6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d.tar.gz Python-Textbook-Companions-6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d.tar.bz2 Python-Textbook-Companions-6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d.zip | |
Removed duplicates
Diffstat (limited to 'Mechanical_Engineering_Thermodynamics_by_Mooney_D_A')
24 files changed, 4844 insertions, 0 deletions
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 Binary files differnew file mode 100755 index 00000000..0eed9ce1 --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap23.png diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap24.png b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap24.png Binary files differnew file mode 100755 index 00000000..4ff3b58c --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap24.png diff --git a/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap25.png b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap25.png Binary files differnew file mode 100755 index 00000000..850e306c --- /dev/null +++ b/Mechanical_Engineering_Thermodynamics_by_Mooney_D_A/screenshots/chap25.png |
