adding books

1 files changed, 356 insertions, 0 deletions

diff --git a/Thermodynamics_by_F_P_Durham/chapter8.ipynb b/Thermodynamics_by_F_P_Durham/chapter8.ipynb
new file mode 100755
index 00000000..8fd523c3
--- /dev/null
+++ b/Thermodynamics_by_F_P_Durham/chapter8.ipynb
@@ -0,0 +1,356 @@
+{
+ "metadata": {
+  "name": "",
+  "signature": "sha256:faffd56df5ff320ba3aaf8ef0438e5c3fd03dbfd20ea0e30132ca53319ffe641"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+  {
+   "cells": [
+    {
+     "cell_type": "heading",
+     "level": 1,
+     "metadata": {},
+     "source": [
+      "Chapter 8: Engine Cycles"
+     ]
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 8.1, page no. 154"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialization\n",
+      "ratio = 7.0\n",
+      "Q = 300.0                       #B/lbm\n",
+      "T1 = 60+460.0                   #R\n",
+      "P1 = 14.7                       #lb/in^2\n",
+      "cv = 0.1715                     #B/lvm F\n",
+      "g = 1.4\n",
+      "\n",
+      "#calculation\n",
+      "Tratio = (ratio)**(g-1)\n",
+      "T2 = Tratio*T1\n",
+      "T3 = T2+Q/cv\n",
+      "eta = round(1- 1/Tratio,2)\n",
+      "WbyJ = eta*Q\n",
+      "print WbyJ*778\n",
+      "W = 778*WbyJ\n",
+      "\n",
+      "#result\n",
+      "print \"Final temperature  =  %d R\" %T3\n",
+      "print \"Thermal efficiency  =  %.3f\" %eta\n",
+      "print \"Work done  =  %d ft-lb/lbm\" %W\n",
+      "#difference in answer due to internal rounding off in Python"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "126036.0\n",
+        "Final temperature  =  2881 R\n",
+        "Thermal efficiency  =  0.540\n",
+        "Work done  =  126036 ft-lb/lbm\n"
+       ]
+      }
+     ],
+     "prompt_number": 17
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 8.2, page no. 157"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "import math\n",
+      "\n",
+      "#initialization\n",
+      "cydia = 3.0                     #in\n",
+      "crdia = 5.0                     #in\n",
+      "ratio = 7.0\n",
+      "rpm = 3000.0                    #rpm\n",
+      "hp = 50.0                       #hp\n",
+      "w = 24.2                        #lbm\n",
+      "Q = 18000.0                     #B/lbm\n",
+      "P1 = 14.7                       #lb/in^2\n",
+      "T1 = 60+460.0                   #R\n",
+      "g = 1.4\n",
+      "cv = 0.1715\n",
+      "\n",
+      "#calculation\n",
+      "eta = hp*550*3600/(778*w*Q)\n",
+      "vol = math.pi*(cydia/12)**2 *(crdia/12)*6/4\n",
+      "vdot = vol*rpm/(60*2)\n",
+      "v1 = 53.3*T1/(144*P1)\n",
+      "wdot = vdot/v1\n",
+      "Qdot = w*Q/3600\n",
+      "Qdash = Qdot/wdot\n",
+      "T2 = T1*(ratio)**(g-1)\n",
+      "T3 = T2+Qdash/cv\n",
+      "eta2 = 1- 1/(ratio)**(g-1)\n",
+      "WbyJ = eta2*Qdot\n",
+      "Wdot = WbyJ*778/550.0\n",
+      "\n",
+      "#result\n",
+      "\n",
+      "print \"Part a\"\n",
+      "print \"Thermal efficiency  =  %.3f \" %eta\n",
+      "\n",
+      "print \"part b\"\n",
+      "print \"Temperature at the end of compression  =  %d R\" %T2\n",
+      "print \"Power developed  =  %.1f hp\" %Wdot"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Part a\n",
+        "Thermal efficiency  =  0.292 \n",
+        "part b\n",
+        "Temperature at the end of compression  =  1132 R\n",
+        "Power developed  =  92.6 hp\n"
+       ]
+      }
+     ],
+     "prompt_number": 2
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 8.3, page no. 161"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialization\n",
+      "Pi = 14.0                       #lb/in^2\n",
+      "T1 = 70+460.0                   #F\n",
+      "ratio = 13.0\n",
+      "T3 = 2500+460.0                 #F\n",
+      "cv = 0.171\n",
+      "cp = 0.23\n",
+      "R = 53.3\n",
+      "g = 1.4\n",
+      "\n",
+      "#calculation\n",
+      "T2 = T1*(ratio)**(g-1)\n",
+      "v3ratio = T3/T2\n",
+      "cutoff =  (v3ratio-1)/(ratio-1)\n",
+      "v1ratio = ratio/v3ratio\n",
+      "T4 = T3*(1/v1ratio)**(g-1)\n",
+      "eta = 1.0- cv*(T4-T1)/(T3-T2)/cp\n",
+      "percent = eta*100\n",
+      "\n",
+      "#result\n",
+      "print \"cut off ratio  =  %.4f\" %cutoff\n",
+      "print \"T end expansion  =  %d R\" %T4\n",
+      "print \"Thermal efficiency  =  %.1f\" %percent, \"%\"\n",
+      "#difference in % due to internal rounding off in Python"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "0.563104145624\n",
+        "cut off ratio  =  0.0835\n",
+        "T end expansion  =  1400 R\n",
+        "Thermal efficiency  =  56.3 %\n"
+       ]
+      }
+     ],
+     "prompt_number": 20
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 8.4, page no. 167"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialization\n",
+      "Pratio = 6.0\n",
+      "P = 14.7                        #lb/in^2\n",
+      "Tt1 = 60+460.0                  #R\n",
+      "Tt3 = 1600+460.0                #R\n",
+      "w = 60.0                        #lb/sec\n",
+      "cp = 0.24                       #B/lbm F\n",
+      "g = 1.4\n",
+      "R = 53.3                        #ft-lb/lbm R\n",
+      "\n",
+      "#calculation\n",
+      "Tt2 = Tt1*(Pratio)**((g-1)/g)\n",
+      "Tratio = Tt2/Tt1\n",
+      "Q = cp*(Tt3-Tt2)\n",
+      "eta = 1- 1/Tratio\n",
+      "W = eta*778*Q\n",
+      "Wdot = w*W/550.0\n",
+      "\n",
+      "#result\n",
+      "print \"Thermal efficiency  =  %.3f\" %eta\n",
+      "print \"Horsepower output  =  %d hp\" %Wdot"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thermal efficiency  =  0.401\n",
+        "Horsepower output  =  9731 hp\n"
+       ]
+      }
+     ],
+     "prompt_number": 5
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 8.5, page no. 169"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialization\n",
+      "P  =  14.7        #lb/in^2\n",
+      "T  =  60+460      #R\n",
+      "e1  =  0.8\n",
+      "P2  =  3          #lb/in^2\n",
+      "T2  =  1600+460   #R\n",
+      "Pt4  =  15.6      #lb/in^2\n",
+      "w  =  60          #lbm/sec\n",
+      "e2  =  0.85\n",
+      "\n",
+      "#calculation\n",
+      "# from table 6, initial conditions are\n",
+      "ht1  =  124.3\n",
+      "Prt1  =  1.215\n",
+      "Prt2s  =  6*Prt1\n",
+      "ht2s  =  207.6\n",
+      "ht2 = ht1+(ht2s-ht1)/e1\n",
+      "dht1 = round((ht2s-ht1)/e1, 1)\n",
+      "ht3 = 521.4\n",
+      "Prt3 = 196.2\n",
+      "Pt3 = 6*P-P2\n",
+      "Pratio = Pt3/Pt4\n",
+      "Prt4s = Prt3/Pratio\n",
+      "ht4 = 326.5\n",
+      "dht3 = round(e2*(ht3-ht4), 1)\n",
+      "W = 778*(dht3-dht1)\n",
+      "Q = ht3-ht2\n",
+      "etaf = W/778.0/Q\n",
+      "Wdot = w*W/550.0\n",
+      "\n",
+      "#result\n",
+      "print \"Thermal efficiency  =  %.3f\" %(round(W))\n",
+      "print \"Horsepower output  =  %d hp\" %Wdot\n",
+      "#difference due to internal rounding off in Python."
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thermal efficiency  =  47925.000\n",
+        "Horsepower output  =  5228 hp\n"
+       ]
+      }
+     ],
+     "prompt_number": 11
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Exmaple 8.6, page no. 172"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "\n",
+      "#initialization\n",
+      "g  =  1.4\n",
+      "Tt4  =  2060                      #R\n",
+      "cp  =  0.24\n",
+      "\n",
+      "#calculation\n",
+      "Tt5  =  Tt4/1.67\n",
+      "Tt2  =  868                       #R\n",
+      "Tt3s  =  1234\n",
+      "dTt3  =  (Tt3s-Tt2)/2.0\n",
+      "Tt3  =  Tt2+dTt3\n",
+      "Q  =  cp*(Tt4-Tt3)\n",
+      "eta  =  286*0.401/Q\n",
+      "\n",
+      "#result\n",
+      "\n",
+      "print \"Heat Added is \", round(Q),\"B/lbm\"\n",
+      "print \"eta is \", round(eta, 3)\n",
+      "print \"Improvement is around 6.2 percent in overall efficiency\""
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Heat Added is  242.0 B/lbm\n",
+        "eta is  0.474\n",
+        "Improvement is around 6.2 percent in overall efficiency\n"
+       ]
+      }
+     ],
+     "prompt_number": 1
+    }
+   ],
+   "metadata": {}
+  }
+ ]
+}
\ No newline at end of file