{
 "metadata": {
  "name": "",
  "signature": "sha256:96d9d4935cd6f3c8e019e8d452a9980bb4f36e96fb606a03b9b38b11670ead8c"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 8: Balances on Nonreactive Processes"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.1-1, page no. 362"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "basis=100. #mol/s\n",
      "ninAcv=66.9 #mol/s\n",
      "noutAcv=3.35 #mol/s\n",
      "noutAcl=63.55 #mol/s\n",
      "ninN2=33.1 #mol/s\n",
      "noutN2=ninN2\n",
      "HinAcv=35.7 #Kj/mol\n",
      "HoutAcv=32.0 #Kj/mol\n",
      "HinN2=1.16 #Kj/mol\n",
      "HoutN2= -0.10 #Kj/mol\n",
      "#Calculations and printing:\n",
      "deltaH=noutAcv*HoutAcv+ noutN2*HoutN2-ninN2*HinN2-ninAcv*HinAcv\n",
      "Q=deltaH\n",
      "print '%s %.2f' %(\"Heat transferred from condenser to achieve required cooling and condensation at the rate of (kW) \",Q)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Heat transferred from condenser to achieve required cooling and condensation at the rate of (kW)  -2322.84\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 2,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.3-1, page no. 367"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "import scipy\n",
      "from scipy import integrate\n",
      "mass=200.0 #kg\n",
      "Ti=20.0 #C\n",
      "Tf=150.0 #C\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \\n \")\n",
      "def Cv(T):\n",
      "    Cv=0.855+ T*9.42*math.pow(10,-4)\n",
      "    return Cv\n",
      "\n",
      "Ucap, err= scipy.integrate.quad(Cv,Ti,Tf)  #integrate.quad is an inbult function used for definite integration\n",
      "Q=mass*Ucap\n",
      "print '%s %.3f' %(\"Heat Required (KJ) = \",Q)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \n",
        " \n",
        "Heat Required (KJ) =  24311.820\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 3,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.3-2, page no. 369"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "import scipy\n",
      "from scipy import integrate\n",
      "T1=20.0 #C\n",
      "T2=100.0 #C\n",
      "T3=90.0 #C\n",
      "T4=30.0 #C\n",
      "P=3.0 #bar\n",
      "V=5.0 #L\n",
      "R=0.08314 #L.bar/mol.K\n",
      "T=363.0 #K\n",
      "ndot=100.0 #mol/min\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \\n \")\n",
      "print(\"part 1\")\n",
      "def fun1(T):\n",
      "    fun1=0.02900+ T*0.2199*math.pow(10,-5) + math.pow(T,2) * 0.5723 *math.pow(10,-8) - math.pow(T,3) * 2.871 * math.pow(10,-12)\n",
      "    return fun1\n",
      "\n",
      "deltaH, err=scipy.integrate.quad(fun1,T1,T2) #scipy.integrate.quad is an inbult function used for definite integration\n",
      "Qdot=ndot*deltaH\n",
      "print '%s %.3f' %(\"Heat Transferred (KJ/min) = \",Qdot)\n",
      "print(\"part2\")\n",
      "def fun2(T):\n",
      "    fun2=0.02900+ T*0.2199*math.pow(10,-5) + math.pow(T,2) * 0.5723 *math.pow(10,-8) - math.pow(T,3) * 2.871 * math.pow(10,-12)-8.14* math.pow(10,-3)\n",
      "    return fun2\n",
      "\n",
      "deltaU, err2=scipy.integrate.quad(fun2,T3,T4) #scipy.integrate.quad is an inbult function used for definite integration\n",
      "n=P*V/(R*T)\n",
      "Q=n*deltaU\n",
      "print '%s %.3f' %(\"Heat transferred (KJ) = \",Q)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \n",
        " \n",
        "part 1\n",
        "Heat Transferred (KJ/min) =  233.238\n",
        "part2\n",
        "Heat transferred (KJ) =  -0.627\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 4,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.3-3, page no. 371"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "import scipy\n",
      "from scipy import integrate\n",
      "T1=430.0 #C\n",
      "T2=100.0 #C\n",
      "ndot=15.0 #Kmol/min\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \\n \")\n",
      "def fun(T):\n",
      "    fun=0.02894 + T* 0.4147 *math.pow(10,-5) + math.pow(T,2) * 0.3191 * math.pow(10,-8) - math.pow(T,3) * 1.965 * math.pow(10,-12)\n",
      "    return fun\n",
      "\n",
      "deltaH, err=scipy.integrate.quad(fun,T1,T2)  #scipy.integrate.quad is an inbult function used for definite integration\n",
      "Qdot=ndot*deltaH*math.pow(10,3) /60.\n",
      "print '%s %.3f' %(\" \\n Rate of heat removal (KW) = \",Qdot)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \n",
        " \n",
        " \n",
        " Rate of heat removal (KW) =  -2494.904\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 5,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.3-4, page no. 373"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "import scipy\n",
      "from scipy import integrate\n",
      "x=0.6\n",
      "T1=0.0\n",
      "T2=400.0\n",
      "ndot=150.0 #mol/h\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \\n \")\n",
      "def fun1(T):\n",
      "    fun1=0.04937 + T*13.92*math.pow(10,-5) - math.pow(T,2) *5.816*math.pow(10,-8) + math.pow(T,3) *7.280 * math.pow(10,-12)\n",
      "    return fun1\n",
      "\n",
      "def fun2(T):\n",
      "    fun2=0.06803 +T*22.59*math.pow(10,-5) - math.pow(T,2) *13.11*math.pow(10,-8) + math.pow(T,3) *31.71 * math.pow(10,-12)\n",
      "    return fun2\n",
      "\n",
      "def fun(T):\n",
      "    fun=x*fun1(T)+ (1-x)*fun2(T)\n",
      "    return fun\n",
      "\n",
      "deltaH , err=scipy.integrate.quad(fun,T1,T2) #scipy.integrate.quad is an inbult function used for definite integration\n",
      "print '%s %.3f' %(\" \\n Heat capacity of Mixture (KJ/mol) = \",deltaH)\n",
      "Qdot=ndot*deltaH\n",
      "print '%s %.3f' %(\" \\n Heat Required (KJ/h) = \",Qdot)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \n",
        " \n",
        " \n",
        " Heat capacity of Mixture (KJ/mol) =  34.890\n",
        " \n",
        " Heat Required (KJ/h) =  5233.495\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 6,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.3-5, page no. 374"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "import scipy\n",
      "from scipy import integrate\n",
      "x=0.1 #CH4\n",
      "T1=20.0\n",
      "T2=300.0\n",
      "Vdot=2000.0 #L/min\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \\n \")\n",
      "ndot=Vdot/22.4\n",
      "def fun(T):\n",
      "    fun=0.03431 + T*5.469*math.pow(10,-5) + math.pow(T,2) *0.3661*math.pow(10,-8) + math.pow(T,3) *11*math.pow(10,-12)\n",
      "    return fun\n",
      "\n",
      "H1 ,err=scipy.integrate.quad(fun,T1,T2)  #scipy.integrate.quad is an inbult function used for definite integration\n",
      "print '%s %.3f' %(\"H1 (Kj/mol) = \",H1)\n",
      "print(\"From Tables H2= -0.15 Kj/mol   ,   H3=8.17 Kj/mol\")\n",
      "H2= -0.15\n",
      "H3=8.17\n",
      "Qdot=(ndot*x*H1+ ndot*(1-x)*(H3-H2))/60.\n",
      "print '%s %.3f' %(\"Heat Input (KW) = \",Qdot)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \n",
        " \n",
        "H1 (Kj/mol) =  12.112\n",
        "From Tables H2= -0.15 Kj/mol   ,   H3=8.17 Kj/mol\n",
        "Heat Input (KW) =  12.945\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 7,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.4-1, page no. 379"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "\n",
      "mdot=1500.0 #g/min\n",
      "M=32.0 #g/mol\n",
      "deltaHv=35.3 #Kj/mol\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \\n \")\n",
      "Qdot=mdot*deltaHv/(M*60.)\n",
      "print '%s %.3f' %(\"Rate of Heat transfer (KW) = \",Qdot)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook \n",
        " \n",
        "Rate of Heat transfer (KW) =  27.578\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 9,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.4-2, page no. 379"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "import scipy\n",
      "from scipy import integrate\n",
      "deltaHv=28.85  #Kj/mol at 69 C\n",
      "T1=25.0 #C\n",
      "T2=69.0 #C\n",
      "Cp=0.2163 #Kj/mol C\n",
      "V=1.0 #L\n",
      "P=7.0 #bar\n",
      "D=0.659 #KG/L\n",
      "M=86.17 #Kg\n",
      "ndot=100.0 #mol/h\n",
      "T3=300.0 #C\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\")\n",
      "print(\"We have deltaHv at 69C hence We follow path ADG\")\n",
      "deltaHA=Cp*(T2-T1) + V*(1.013-P)*M/(D*math.pow(10,4))\n",
      "print '%s %.3f' %(\" \\n deltaHA (Kj/mol) =\",deltaHA)\n",
      "deltaHD=deltaHv\n",
      "print '%s %.3f' %(\" \\n deltaHD (Kj/Kg) =\",deltaHD)\n",
      "def fun1(T):\n",
      "    fun1=0.13744 + T*40.85*math.pow(10,-5) - math.pow(T,2) *23.92*math.pow(10,-8) + math.pow(T,3) *57.66*math.pow(10,-12)\n",
      "    return fun1\n",
      "\n",
      "deltaHG, err=scipy.integrate.quad(fun1,T2,T3)  #scipy.integrate.quad is an inbult function used for definite integration\n",
      "print '%s %.3f' %(\" \\n deltaHG (KJ/mol) = \",deltaHG)\n",
      "Qdot=ndot*(deltaHA+deltaHD+deltaHG)/3600\n",
      "print '%s %.3f' %(\" \\n rate of Heat supply (KJ/mol) = \",Qdot)\n",
      "print(\" \\n In this problem we neglected V*deltaP as it is negligible \")\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\n",
        "We have deltaHv at 69C hence We follow path ADG\n",
        " \n",
        " deltaHA (Kj/mol) = 9.439\n",
        " \n",
        " deltaHD (Kj/Kg) = 28.850\n",
        " \n",
        " deltaHG (KJ/mol) =  47.149\n",
        " \n",
        " rate of Heat supply (KJ/mol) =  2.373\n",
        " \n",
        " In this problem we neglected V*deltaP as it is negligible \n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 10,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.4-3, page no. 382"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "\n",
      "T1=337.9 #K\n",
      "T2=473.0 #K\n",
      "Tc=513.2 #K\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\")\n",
      "print(\"using Trouton rule,\")\n",
      "deltaHvT1=0.109*T1\n",
      "print(\"In this case, trouton rule gives a better estimate\")\n",
      "print(\"using Watson correction\")\n",
      "deltaHvT2=36.8*math.pow(((Tc-T2)/(Tc-T1)),0.38)\n",
      "print '%s %.3f' %(\"Estimated value using Trouton rule (Kj/mol) = \",deltaHvT1)\n",
      "print '%s %.3f' %(\" \\n Estimated value using watson correction (Kj/mol) = \",deltaHvT2)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\n",
        "using Trouton rule,\n",
        "In this case, trouton rule gives a better estimate\n",
        "using Watson correction\n",
        "Estimated value using Trouton rule (Kj/mol) =  36.831\n",
        " \n",
        " Estimated value using watson correction (Kj/mol) =  21.029\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 11,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.4-4, page no. 383"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "\n",
      "basis=1.0 #mol feed\n",
      "x=0.684 #mole fraction Of B\n",
      "y=0.4\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\")\n",
      "A=([[1, 1],[x ,y]])\n",
      "b=([[basis],[basis/2]])\n",
      "C=numpy.dot(linalg.inv(A),b)\n",
      "#Here We solved two linear equations simultaneously\n",
      "nV=C[0,0]\n",
      "nL=C[1,0]\n",
      "H1=5.332\n",
      "H2=6.340\n",
      "H3=37.52\n",
      "H4=42.93\n",
      "Q=nV*x*H1 + nV*(1-x)*H2 + nL*y*H3 + nL*(1-y)*H4\n",
      "print '%s %.3f' %(\" \\n Heat transferred (KJ) = \",Q)\n",
      "print(\"The answer for this problem in Text is wrong\")\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\n",
        " \n",
        " Heat transferred (KJ) =  28.401\n",
        "The answer for this problem in Text is wrong\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 1,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.4-6, page no. 390"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "basis=1. #lbm dry air\n",
      "Tin=80. #F Inlet temperature\n",
      "Tout=51. #F Outlet temperature\n",
      "\n",
      "\n",
      "\n",
      "#Calculations and printing:\n",
      "print(\"For 80F and 80 percent RH, ha=0.018 lbm H20/lbm DA H1=38.8 Btu/lbm DA\")\n",
      "ha=0.018 #lbm H20/lbm DA\n",
      "H1=38.8 #Btu/lbm DA\n",
      "m1=basis*ha\n",
      "print(\"\\n For 51F and saturated, ha=0.0079 lbm H20/lbm DA H2=20.9 Btu/lbm DA\")\n",
      "ha2=0.0079 #lbm H20/lbm DA\n",
      "H2=20.9 #Btu/lbm DA\n",
      "m2=basis*ha2\n",
      "print(\"\\n Balance on H2O\")\n",
      "m3=m1-m2\n",
      "Fraction=m3/m1\n",
      "print '%s %.3f' %(\"\\n Fraction H2O condensed = \",Fraction)\n",
      "H3=basis*(Tout-32.)\n",
      "Q=basis*H2+m3*H3-basis*H1\n",
      "print(\"\\n From psychrometric chart, Vh=13.0 ft^3/lbm DA\")\n",
      "Vbasis=basis*13.0\n",
      "Qdot=Q*1000./Vbasis\n",
      "print '%s %.2f' %(\"Rate at which heat must be removed (Btu/min) = \",Qdot)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "For 80F and 80 percent RH, ha=0.018 lbm H20/lbm DA H1=38.8 Btu/lbm DA\n",
        "\n",
        " For 51F and saturated, ha=0.0079 lbm H20/lbm DA H2=20.9 Btu/lbm DA\n",
        "\n",
        " Balance on H2O\n",
        "\n",
        " Fraction H2O condensed =  0.561\n",
        "\n",
        " From psychrometric chart, Vh=13.0 ft^3/lbm DA\n",
        "Rate at which heat must be removed (Btu/min) =  -1362.16\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 13,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 13
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.5-1, page no. 397"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "import scipy\n",
      "from scipy import integrate\n",
      "M1=35.6 #g/mol\n",
      "M2=18.0 #g/mol\n",
      "x=0.2\n",
      "mdot=1000.0 #kg/h\n",
      "T1=25.0 #C\n",
      "T2=100.0 #C\n",
      "T3=40.0 #C\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\")\n",
      "nHCl=mdot*x*math.pow(10,3) /M1\n",
      "nH2O=mdot*(1-x)*math.pow(10,3) /M2\n",
      "def fun(T):\n",
      "    fun=29.13*math.pow(10,-3) - T*0.1341*math.pow(10,-5) +math.pow(T,2) *0.9715*math.pow(10,-8) - math.pow(T,3) *4.335*math.pow(10,-12)\n",
      "    return fun\n",
      "\n",
      "H1, err=scipy.integrate.quad(fun,T1,T2)  #scipy.integrate.quad is an inbult function used for definite integration\n",
      "r=nH2O/nHCl\n",
      "print(\"From table B.11, deltaHa= -67.4 Kj/mol HCl\")\n",
      "deltaHa= -67.4\n",
      "y=nHCl/(nHCl+nH2O)\n",
      "Cp=0.73*mdot*4.184/nHCl\n",
      "deltaHb=Cp*(T3-T1)\n",
      "H2=deltaHa+deltaHb\n",
      "Qdot=nHCl*(H2-H1)\n",
      "print '%s %.3E' %(\"Heat (kJ/h) = \",Qdot)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\n",
        "From table B.11, deltaHa= -67.4 Kj/mol HCl\n",
        "Heat (kJ/h) =  -3.451E+05\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 14,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 14
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.5-2, page no. 400"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "\n",
      "x=0.05\n",
      "y=0.4\n",
      "mdot=1000.0 #lbm/h\n",
      "Hf=10.0 #Btu/lbm\n",
      "Hl= -17 #Btu/lbm\n",
      "Hv=1138.0 #Btu/lbm\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\")\n",
      "print(\"\\n sulphuric acid balance\")\n",
      "m2=x*mdot/y\n",
      "print(\"\\n Total Mass balance\")\n",
      "m1=mdot-m2\n",
      "Qdot=m1*Hv+m2*Hl-mdot*Hf\n",
      "print '%s %.3f' %(\" Rate of Heat transfer (Btu/h) = \",Qdot)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\n",
        "\n",
        " sulphuric acid balance\n",
        "\n",
        " Total Mass balance\n",
        " Rate of Heat transfer (Btu/h) =  983625.000\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 15,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 15
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.5-5, page no. 405"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "\n",
      "#Initialization of variables\n",
      "import math\n",
      "import numpy\n",
      "from numpy import linalg\n",
      "\n",
      "basis=100.0 #lbm/h\n",
      "Hv=728.0 #Btu/lbm\n",
      "Hl=45.0 #Btu/lbm\n",
      "HF=100.0 #Btu/lbm\n",
      "T=120.0 #F\n",
      "xF=0.30\n",
      "\n",
      "#Calculations and printing :\n",
      "\n",
      "print(\" All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\")\n",
      "print(\"From figure 8.5-2, \")\n",
      "xL=0.185\n",
      "xV=0.89\n",
      "mL=basis*((xV-xF)/(xV-xL))\n",
      "mV=basis-mL\n",
      "Qdot=mV*Hv + mL*Hl - basis*HF\n",
      "print '%s %.3f' %(\"Rate of heat transfer (Btu/h) =\",Qdot)\n",
      "raw_input('press enter key to exit')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " All the values in the textbook are Approximated hence the values in this code differ from those of Textbook\n",
        "From figure 8.5-2, \n",
        "Rate of heat transfer (Btu/h) = 5641.135\n"
       ]
      },
      {
       "name": "stdout",
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "press enter key to exit\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 16,
       "text": [
        "''"
       ]
      }
     ],
     "prompt_number": 16
    }
   ],
   "metadata": {}
  }
 ]
}