adding book

1 files changed, 219 insertions, 0 deletions

diff --git a/Introduction_To_Chemical_Engineering/ch9.ipynb b/Introduction_To_Chemical_Engineering/ch9.ipynb
new file mode 100644
index 00000000..ebdb2cdc
--- /dev/null
+++ b/Introduction_To_Chemical_Engineering/ch9.ipynb
@@ -0,0 +1,219 @@
+{
+ "metadata": {
+  "name": ""
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+  {
+   "cells": [
+    {
+     "cell_type": "heading",
+     "level": 1,
+     "metadata": {},
+     "source": [
+      "Chapter 9 : Computers and their application"
+     ]
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "example 9.1 page number 384"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "#to find the pressure drop in the coil\n",
+      "\n",
+      "import math \n",
+      "# Variables\n",
+      "D = 38.*10**-3;    #in m\n",
+      "U = 1.   #in m/s\n",
+      "density = 998.   #in kg/cubic m\n",
+      "viscosity = 8.*10**-4   #in Pa-s\n",
+      "DC = 1.   #in m\n",
+      "N = 10.\n",
+      "e = 4.*10**-6;   #in m\n",
+      "\n",
+      "# Calculations and Results\n",
+      "Re = (density*U*D)/viscosity;\n",
+      "print  \"Reynolds number = %f\"%(Re)\n",
+      "\n",
+      "f = (4*math.log10((e/D)/3.7+(6.81/Re)**0.9))**-2;\n",
+      "print \"friction factor = %f\"%(f);\n",
+      "\n",
+      "L = 3.14*DC*N;\n",
+      "\n",
+      "delta_Pstr = (2*f*U*density*L)/D;\n",
+      "print \"pressure drop through straight pipe = %f Pa\"%(delta_Pstr)\n",
+      "\n",
+      "S = 1+3.54*(D/DC);\n",
+      "print \"correction factor = %f\"%(S)\n",
+      "\n",
+      "delta_P = S*delta_Pstr\n",
+      "print \"pressure drop of coil = %f Pa\"%(delta_P)\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Reynolds number = 47405.000000\n",
+        "friction factor = 0.005330\n",
+        "pressure drop through straight pipe = 8791.184173 Pa\n",
+        "correction factor = 1.134520\n",
+        "pressure drop of coil = 9973.774268 Pa\n"
+       ]
+      }
+     ],
+     "prompt_number": 1
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "example 9.2 page number 384\n"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "#to find the shell side pressure drop in heat exchanger\n",
+      "\n",
+      "import math \n",
+      "# Variables\n",
+      "U = 0.5   #in m/s\n",
+      "N = 19.;\n",
+      "DT = 0.026   #in m\n",
+      "L = 2.7   #in m\n",
+      "DS = 0.2   #in m\n",
+      "e = 0.0002   #in m\n",
+      "density = 836.  #in kg/cu m\n",
+      "viscosity = 0.00032   #in Pa s\n",
+      "Pr = 6.5;\n",
+      "Prw = 7.6;\n",
+      "\n",
+      "# Calculations and Results\n",
+      "HYDIA = (DS**2-N*DT**2)/(DS+N*DT);\n",
+      "Re = HYDIA*U*density/viscosity;\n",
+      "print  \"Reynolds number = %f\"%(Re)\n",
+      "\n",
+      "f = (4*math.log10((e/HYDIA)/3.7+(6.81/Re)**0.9))**-2;\n",
+      "print \"friction factor = %f\"%(f);\n",
+      "\n",
+      "L = 3.14*DT*N;\n",
+      "\n",
+      "delta_Pstr = (2*f*U*density*L)/HYDIA;\n",
+      "print \"pressure drop through straight pipe = %f Pa\"%(delta_Pstr)\n",
+      "\n",
+      "S = (Prw/Pr)**0.33;\n",
+      "print \"correction factor = %f\"%(S)\n",
+      "\n",
+      "delta_P = S*delta_Pstr\n",
+      "print \"pressure drop of coil = %f Pa\"%(delta_P)\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Reynolds number = 51113.148415\n",
+        "friction factor = 0.008158\n",
+        "pressure drop through straight pipe = 270.362537 Pa\n",
+        "correction factor = 1.052948\n",
+        "pressure drop of coil = 284.677794 Pa\n"
+       ]
+      }
+     ],
+     "prompt_number": 2
+    },
+    {
+     "cell_type": "heading",
+     "level": 3,
+     "metadata": {},
+     "source": [
+      "example 9.3 page number 385\n"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "# Variables\n",
+      "MH = 10.   #in kg/s\n",
+      "MC = 12.5  #in kg/s\n",
+      "CPH = 4.2   #in kJ/kg\n",
+      "CPC = 4.2   #in kJ/kg\n",
+      "THI = 353.   #in K\n",
+      "THO = 333.   #in K\n",
+      "TCI = 300.   #in K\n",
+      "U = 1.8    #in kW/sq m K\n",
+      "\n",
+      "# Calculations and Results\n",
+      "Q = MH*CPH*(THI-THO);\n",
+      "print \"heat load = %f J\"%(Q)\n",
+      "\n",
+      "TCO = Q/(MC*CPC)+TCI;\n",
+      "print \"cold fluid outlet temperature = %f K\"%(TCO)\n",
+      "\n",
+      "#for co current flow\n",
+      "\n",
+      "DT1 = THI-TCO;\n",
+      "DT2 = THO-TCO;\n",
+      "\n",
+      "LMTD = (DT1-DT2)/math.log(DT1/DT2);\n",
+      "\n",
+      "A = Q/(U*LMTD);\n",
+      "print \"for co current flow area = %f sq m\"%(A);\n",
+      "\n",
+      "#for counter current flow\n",
+      "\n",
+      "DT1 = THI-TCO;\n",
+      "DT2 = THO-TCI;\n",
+      "\n",
+      "LMTD = (DT1-DT2)/math.log(DT1/DT2);\n",
+      "\n",
+      "A = Q/(U*LMTD);\n",
+      "print \"for counter current flow area = %f sq m\"%(A);\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "heat load = 840.000000 J\n",
+        "cold fluid outlet temperature = 316.000000 K\n",
+        "for co current flow area = 18.146440 sq m\n",
+        "for counter current flow area = 13.347874 sq m\n"
+       ]
+      }
+     ],
+     "prompt_number": 3
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [],
+     "language": "python",
+     "metadata": {},
+     "outputs": []
+    }
+   ],
+   "metadata": {}
+  }
+ ]
+}
\ No newline at end of file