summaryrefslogtreecommitdiff
path: root/Introduction_To_Chemical_Engineering_Thermodynamics_by_G._Halder/Ch7.ipynb
diff options
context:
space:
mode:
Diffstat (limited to 'Introduction_To_Chemical_Engineering_Thermodynamics_by_G._Halder/Ch7.ipynb')
-rwxr-xr-xIntroduction_To_Chemical_Engineering_Thermodynamics_by_G._Halder/Ch7.ipynb597
1 files changed, 0 insertions, 597 deletions
diff --git a/Introduction_To_Chemical_Engineering_Thermodynamics_by_G._Halder/Ch7.ipynb b/Introduction_To_Chemical_Engineering_Thermodynamics_by_G._Halder/Ch7.ipynb
deleted file mode 100755
index 6e32feac..00000000
--- a/Introduction_To_Chemical_Engineering_Thermodynamics_by_G._Halder/Ch7.ipynb
+++ /dev/null
@@ -1,597 +0,0 @@
-{
- "cells": [
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Chapter 7 - Thermodynamics to flow processes"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.1 Page: 256"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 2,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.1 - Page: 256\n",
- "\n",
- "\n",
- "Flow rate is 0.1237 m/s\n",
- "\n",
- "Velocity of water at the outlet is 3.938 m/s\n"
- ]
- }
- ],
- "source": [
- "from __future__ import division\n",
- "from math import pi\n",
- "print \"Example: 7.1 - Page: 256\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "# *****Data******#\n",
- "d1 = 0.15## [inlet dia, m]\n",
- "d2 = 0.20## [outlet dia, m]\n",
- "U1 = 7## [inlet velocity, m/s]\n",
- "#****************#\n",
- "\n",
- "# From Fig. 7.2 (Pg 256)\n",
- "# At the inlet:\n",
- "A1 = (pi/4)*d1**2## [square m]\n",
- "# At the outlet:\n",
- "A2 = (pi/4)*d2**2## [square m]\n",
- "Q = A1*U1## [cubic m/s]\n",
- "print \"Flow rate is %.4f m/s\\n\"%(Q)#\n",
- "# Using Continuity Eqn.\n",
- "# density1*U1*A1 = Density2*U2*A2\n",
- "# For water: Density1 = Density2. Therefore:\n",
- "U2 = A1*U1/A2#\n",
- "print \"Velocity of water at the outlet is %.3f m/s\"%(U2)#"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.2 Page: 257"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 3,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.2 - Page: 257\n",
- "\n",
- "\n",
- "Discharge through the 10 cm pipe is 0.0942 cubic m/sec\n",
- "\n",
- "Average velocity in the 15 cm pipe is 6.38 m/s\n"
- ]
- }
- ],
- "source": [
- "print \"Example: 7.2 - Page: 257\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "#*****Data******#\n",
- "d1 = 0.2## [m]\n",
- "d2 = 0.15## [m]\n",
- "d3 = 0.1## [m]\n",
- "U1 = 3## [m/s]\n",
- "U2 = 2.5## [m/s]\n",
- "#**************#\n",
- "\n",
- "# From Fig. 7.3 (Pg: 257)\n",
- "# For pipe I:\n",
- "A1 = (pi/4)*d1**2## [square m]\n",
- "Q1 = A1*U1## [cubic m/s]\n",
- "# For pipe II:\n",
- "A2 = (pi/4)*d2**2## [square m]\n",
- "Q2 = A2*U2## [cubic m/s]\n",
- "# For pipe III:\n",
- "A3 = (pi/4)*d3**2## [square m]\n",
- "# From continuity Eqn.:\n",
- "Q3 = Q1 - Q2## [cubic m/s]\n",
- "U3 = Q3/A3## [m/s]\n",
- "print \"Discharge through the 10 cm pipe is %.4f cubic m/sec\\n\"%(Q1)#\n",
- "print \"Average velocity in the 15 cm pipe is %.2f m/s\"%(U3)#"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.3 Page: 262"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 4,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.3 - Page: 262\n",
- "\n",
- "\n",
- "Pressure at section 2 is 5.20 bar\n"
- ]
- }
- ],
- "source": [
- "print \"Example: 7.3 - Page: 262\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "#*****Data******#\n",
- "d1 = 0.3## [m]\n",
- "d2 = 0715##[m]\n",
- "Q = 40/1000## [cubic m/s]\n",
- "Z1 = 8## [m]\n",
- "Z2 = 6## [m]\n",
- "P1 = 5*10**5## [Pa]\n",
- "density = 1000## [kg/cubic m]\n",
- "g = 9.81## [m/square s]\n",
- "#*************#\n",
- "\n",
- "# From Fig. 7.3 (Pg: 262)\n",
- "A1 = (pi/4)*d1**2## [square m]\n",
- "A2 = (pi/4)*d2**2## [square m]\n",
- "U1 = Q/A1## [m/s]\n",
- "U2 = Q/A2## [m/s]\n",
- "# Applying Bernoulli's equations at sections 1 & 2:\n",
- "P2 = ((U1**2/(2*g) + Z1 + P1/(density*g)) - (U2**2/(2*g) + Z2))*(density*g)## [Pa]\n",
- "print \"Pressure at section 2 is %.2f bar\"%(P2/10**5)#"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.4 Page: 268"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 5,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.4 - Page: 268\n",
- "\n",
- "\n",
- "Power Requirement of the compressor is 3.79 kW\n"
- ]
- }
- ],
- "source": [
- "print \"Example: 7.4 - Page: 268\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "#*****Data******#\n",
- "P1 = 100## [kPa]\n",
- "T1 = 320## [K]\n",
- "P2 = 600## [kPa]\n",
- "T2 = 430## [K]\n",
- "m_dot = 0.03## [kg/s]\n",
- "Qout = 15## [kJ/kg]\n",
- "#*************#\n",
- "\n",
- "# The energy balance around the compressor:\n",
- "# dE_System/dt = Ein - Eout\n",
- "# Since it is a steady state process:\n",
- "# dE_Sysytem/dt = 0\n",
- "# Ein = Eout\n",
- "# Win + m_dot*H1 = Qout + m_dot*H2\n",
- "# Since, Qout = Qout/m\n",
- "# Win = m_dot*(Qout + (H2 - H1))\n",
- "# From enthalpy chart of air:\n",
- "H1 = 320.20## [Enthalpy of air at 320 K, kJ/kg]\n",
- "H2 = 431.43## [Enthalpy of air at 430 K, kJ/kg]\n",
- "Win = m_dot*(Qout + (H2 - H1))## [kW]\n",
- "print \"Power Requirement of the compressor is %.2f kW\"%(Win)#"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.5 Page: 269"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 1,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.5 - Page: 269\n",
- "\n",
- "\n",
- "Work done by reversible adiabatic compression when gama = 1.4 is -96.53 J/g\n",
- "\n",
- "Work done by isothermal compression is -76.03 J/g\n",
- "\n",
- "Work done in single stage compression is -19.49 kW\n"
- ]
- }
- ],
- "source": [
- "from math import log\n",
- "print \"Example: 7.5 - Page: 269\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "#*****Data******#\n",
- "P1 = 100## [kPa]\n",
- "T1 = 250## [K]\n",
- "Q = 0.1## [cubic m/s]\n",
- "P2 = 500## [kPa]\n",
- "M = 44## [g/mol]\n",
- "R = 8.314## [J/mol K]\n",
- "#****************#\n",
- "\n",
- "# Solution (a)\n",
- "# Work done by reversible adiabatic compression, gama = 1.4#\n",
- "gama = 1.4#\n",
- "T2 = T1*(P2/P1)**((gama - 1)/gama)## [K]\n",
- "Wad = (gama*R/(gama - 1))*(T1 - T2)## [J/mol]\n",
- "Wad = Wad/M## [J/g]\n",
- "print \"Work done by reversible adiabatic compression when gama = 1.4 is %.2f J/g\\n\"%(Wad)#\n",
- "\n",
- "# Solution (b)\n",
- "# Work done by isothermal compression:\n",
- "Wiso = - (R*T1)*log(P2/P1)## [J/mol]\n",
- "Wiso = Wiso/M## [J/g]\n",
- "print \"Work done by isothermal compression is %.2f J/g\\n\"%(Wiso)#\n",
- "\n",
- "# Solution (c)\n",
- "# Work done in single stage compression, gama = 1.3:\n",
- "gama = 1.3#\n",
- "V1 = Q## [cubic m]\n",
- "Wsingle_stage = (gama*P1*V1/(gama - 1))*(1-(P2/P1)**((gama - 1)/gama))## [kW]\n",
- "print \"Work done in single stage compression is %.2f kW\"%(Wsingle_stage)#"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.7 Page: 274"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 7,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.7 - Page: 274\n",
- "\n",
- "\n",
- "The mass flow rate of water is 448.70 kg/min\n",
- "\n",
- "The rate of heat transfer is 18777.00 kJ/min\n"
- ]
- }
- ],
- "source": [
- "print \"Example: 7.7 - Page: 274\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "#*****Data******#\n",
- "T_steam1 = 50## [OC]\n",
- "T_steam2 = 30## [OC]\n",
- "msteam_dot = 10## [kg/min]\n",
- "T_water1 = 15## [OC]\n",
- "T_water2 = 25## [OC]\n",
- "#***************#\n",
- "\n",
- "# Solution (a)\n",
- "# From the Stem Table:\n",
- "H1 = 2645.9## [kJ/kg, At 50 OC]\n",
- "H2 = 768.2## [kJ/kg, At 30 OC]\n",
- "H3 = 62.982## [kJ/kg, At 15 OC]\n",
- "H4 = 104.83## [kJ/kg, At 25 OC]\n",
- "# The mass & Energy balance of the above flow gives:\n",
- "mwater_dot = msteam_dot*(H1 - H2)/(H4 - H3)## [kg/min]\n",
- "print \"The mass flow rate of water is %.2f kg/min\\n\"%(mwater_dot)#\n",
- "\n",
- "# Solution (b)\n",
- "Qdot = mwater_dot*(H4 - H3)## [kJ/min]\n",
- "print \"The rate of heat transfer is %.2f kJ/min\"%(Qdot)#"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.8 Page: 279"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 9,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.8 - Page: 279\n",
- "\n",
- "\n",
- "Outlrt velocity is 603.3 m/s\n",
- "\n"
- ]
- }
- ],
- "source": [
- "from math import sqrt\n",
- "print \"Example: 7.8 - Page: 279\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "#*****Data******#\n",
- "P1 = 500## [kPa]\n",
- "T1 = 623## [K]\n",
- "mdot = 12## [kg/s]\n",
- "P2 = 500## [kPa]\n",
- "T2 = 523## [K]\n",
- "Qdot = -120## [kW]\n",
- "H1 = 3168## [kJ/kg]\n",
- "H2 = 2976## [kJ/kg]\n",
- "#************#\n",
- "\n",
- "Q = Qdot/mdot## [kJ/kg]\n",
- "# By energy balance:\n",
- "# (deltaU**2/2) + g*deltaZ + deltaH = Q - Ws\n",
- "# Considering negligible change in P.E., deltaZ = 0 & Ws = 0.\n",
- "# (U2**2 - U1**2)/2 + deltaH = Q\n",
- "deltaH = H2 - H1## [kJ/kg]\n",
- "# Neglecting inlet velocity.\n",
- "U2 = sqrt(2*(Q - deltaH)*1000)## [m/s]\n",
- "print \"Outlrt velocity is %.1f m/s\\n\"%(U2)#"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.9 Page: 279"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 10,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.9 - Page: 279\n",
- "\n",
- "\n",
- "Critical Ratio is 0.53\n",
- "\n",
- "The discharge velocity is 1174.46 m/s\n"
- ]
- }
- ],
- "source": [
- "print \"Example: 7.9 - Page: 279\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "#*****Data******#\n",
- "Pin = 1000## [kPa]\n",
- "Tin = 600## [K]\n",
- "Uin = 50## [m/s]\n",
- "gama = 1.4#\n",
- "M = 17## [g/mol]\n",
- "R = 8314## [kJ/mol K]\n",
- "MachNumber = 2#\n",
- "#***************#\n",
- "\n",
- "# Solution (i)\n",
- "# Using Eqn. (7.36):\n",
- "Critical_Ratio = (2/(gama + 1))**(gama/(gama - 1))#\n",
- "print \"Critical Ratio is %.2f\\n\"%(Critical_Ratio)#\n",
- "\n",
- "# Solution (ii)\n",
- "PV_in = R*Tin/M## [square m]\n",
- "Uthroat = sqrt(Uin**2 + (2*gama*PV_in/(gama - 1))*(1-(Critical_Ratio)**((gama - 1)/gama)))## [m/s]\n",
- "Uout = MachNumber*Uthroat## [m/s]\n",
- "print \"The discharge velocity is %.2f m/s\"%(Uout)#"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.10 Page: 280"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 11,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.10 - Page: 280\n",
- "\n",
- "\n",
- "(a) Final Velocity is 901.11 m/s\n",
- "\n",
- "(b) Final Velocity is 1553.06 m/s\n",
- "\n"
- ]
- }
- ],
- "source": [
- "print \"Example: 7.10 - Page: 280\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "#*****Data******#\n",
- "P1 = 800## [kPa]\n",
- "T1 = 773## [K]\n",
- "H1 = 3480## [kJ/kg]\n",
- "P2 = 100## [kPa]\n",
- "T2 = 573## [K]\n",
- "H2 = 3074## [kJ/kg]\n",
- "#***************#\n",
- "\n",
- "# Solution (a)\n",
- "# Velocity of the fluid exiting the nozzle:\n",
- "# U2 = sqrt(U1**2 + 2*(H1 - H2))\n",
- "# Neglecting initial velocity:\n",
- "U2 = sqrt(2*(H1 - H2)*1000)## [m/s]\n",
- "print \"(a) Final Velocity is %.2f m/s\\n\"%(U2)#\n",
- "\n",
- "# Solution (b)\n",
- "U1 = 40## [m/s]\n",
- "U2 = sqrt((U1**2 + 2*(H1 - H2))*1000)## [m/s]\n",
- "print \"(b) Final Velocity is %.2f m/s\\n\"%(U2)#"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "## Example: 7.11 Page: 281"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 12,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Example: 7.11 - Page: 281\n",
- "\n",
- "\n",
- "Mass flow rate of the steam is 17.50 kg/s\n",
- "\n",
- "The temperature of the steam leaving the outlet is 393.38 K\n",
- "\n",
- "Area of diffuser outlet is 0.28 square m\n",
- "\n"
- ]
- }
- ],
- "source": [
- "print \"Example: 7.11 - Page: 281\\n\\n\"\n",
- "\n",
- "# Solution\n",
- "\n",
- "#*****Data******#\n",
- "P1 = 100## [kPa]\n",
- "T1 = 200## [OC]\n",
- "U1 = 190## [m/s]\n",
- "A1 = 2000/10**4## [square m]\n",
- "U2 = 70## [m/s]\n",
- "P2 = 200## [kPa]\n",
- "Qdot = 100## [kW]\n",
- "V1 = 2.172## [cubic m/kg]\n",
- "H1 = 2875.3## [kJ/kg]\n",
- "#***************#\n",
- "\n",
- "# Solution (a)\n",
- "mdot = U1*A1/V1## [kg/s]\n",
- "print \"Mass flow rate of the steam is %.2f kg/s\\n\"%(mdot)#\n",
- "\n",
- "# Solution (b)\n",
- "# Amount of heat transferred to the surrounding per unit steam:\n",
- "Q = Qdot/mdot## [kJ/kg]\n",
- "# The Enthalpy at the diffuser outlet can be estimated as:\n",
- "H2 = Q + H1 + (U1**2 - U2**2)/2## [kJ/kg]\n",
- "# From the steam table:\n",
- "T2 = 393.38## [K]\n",
- "V2 = 1.123## [cubic m/kg]\n",
- "print \"The temperature of the steam leaving the outlet is %.2f K\\n\"%(T2)#\n",
- "\n",
- "# Solution (c)\n",
- "A2 = V2*mdot/U2## [square m]\n",
- "print \"Area of diffuser outlet is %.2f square m\\n\"%(A2)#"
- ]
- }
- ],
- "metadata": {
- "kernelspec": {
- "display_name": "Python 2",
- "language": "python",
- "name": "python2"
- },
- "language_info": {
- "codemirror_mode": {
- "name": "ipython",
- "version": 2
- },
- "file_extension": ".py",
- "mimetype": "text/x-python",
- "name": "python",
- "nbconvert_exporter": "python",
- "pygments_lexer": "ipython2",
- "version": "2.7.9"
- }
- },
- "nbformat": 4,
- "nbformat_minor": 0
-}