summaryrefslogtreecommitdiff
path: root/Applied_Thermodynamics/Chapter13.ipynb
diff options
context:
space:
mode:
Diffstat (limited to 'Applied_Thermodynamics/Chapter13.ipynb')
-rwxr-xr-xApplied_Thermodynamics/Chapter13.ipynb916
1 files changed, 916 insertions, 0 deletions
diff --git a/Applied_Thermodynamics/Chapter13.ipynb b/Applied_Thermodynamics/Chapter13.ipynb
new file mode 100755
index 00000000..1d963cc2
--- /dev/null
+++ b/Applied_Thermodynamics/Chapter13.ipynb
@@ -0,0 +1,916 @@
+{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:3c23e07bbff73d50716bb3f0344f5f19803ca8de755edfb991f464d2a6411a44"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter 13: Nozzles"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1, page no. 584"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "from __future__ import division\n",
+ " \n",
+ "#Variable Declaration: \n",
+ "p1 = 10 #Pressure of dry steam(in bar):\n",
+ "C1 = 100 #Velocity of steam entering(in m/s):\n",
+ "C2 = 300 #Velocity of steam leaving the nozzle(in m/s):\n",
+ "p2 = 5 #Pressure of steam at exit(in bar):\n",
+ "m = 16 #Mass flow rate(in kg/s):\n",
+ "q = 10 #Heat loss to surroundings(in kJ/kg):\n",
+ "#From steam tables:\n",
+ "h1 = 2778.1 #kJ/kg\n",
+ "hf = 640.23 #kJ/kg\n",
+ "hfg = 2108.5 #kJ/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "dh = (q*10**3+(C1**2-C2**2)/2)/1000\t#Heat drop in the nozzle(in kJ/kg):\n",
+ "dQ = -dh*m #Total heat drop(in kJ/s):\n",
+ "h2 = h1+dh #Enthalpy at state 2(in kJ/kg):\n",
+ "x2 = (h2-hf)/hfg #Dryness fraction at state 2:\n",
+ "\n",
+ "#Results: \n",
+ "print \"Heat drop in the nozzle: \",round(-dh,2),\"kJ/kg\" \n",
+ "print \"Total heat drop: \",round(dQ,2),\"kJ/s\"\n",
+ "print \"Dryness fraction at exit: \",round(x2,4)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Heat drop in the nozzle: 30.0 kJ/kg\n",
+ "Total heat drop: 480.0 kJ/s\n",
+ "Dryness fraction at exit: 0.9997\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 2, page no. 585"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "from __future__ import division\n",
+ " \n",
+ "from math import sqrt\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "p1 = 10 #Steam entering at pressure(in bar):\n",
+ "p2 = 6 #Pressure at which steam leaves(in bar):\n",
+ "A2 = 20 #Cross-section area of exit of nozzle(in cm**2):\n",
+ "#From steam tables:\n",
+ "h1 = 3478.5 #kJ/kg \n",
+ "s1 = 7.7622 #kJ/kg.K\n",
+ "T2 = 418.45 #C(by interpolation)\n",
+ "h2 = 3309.51 #kJ/kg\n",
+ "v2 = 0.5281 #m**3/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "s2 = s1\n",
+ "C2 = sqrt(2*(h1-h2)*10**3) #Velocity at exit(in m/s):\n",
+ "m = A2*10**(-4)*C2/v2 #Mass flow rate(in kg/s):\n",
+ "\n",
+ "#Result: \n",
+ "print \"Mass flow rate: \",round(m,3),\"kg/s\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Mass flow rate: 2.202 kg/s\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3, page no. 587"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "from __future__ import division\n",
+ "\n",
+ "from math import sqrt \n",
+ "\n",
+ "#Variable Declaration: \n",
+ "p1 = 12 #Pressure of steam entering(in bar):\n",
+ "p2 = 6 #Pressure at exit(in bar):\n",
+ "m1 = 5 #Mass flow rate(in kg/s):\n",
+ "m2 = m1\n",
+ "m3 = m1\n",
+ "C3a = 500 #Exit velocity(in m/s):\n",
+ "#From steam tables:\n",
+ "h1 = 3045.8 #kJ/kg \n",
+ "h2 = 2900.05 #kJ/kg\n",
+ "s2 = 7.0317 #kJ/kg.K\n",
+ "v2 = 0.3466 #m**3/kg\n",
+ "h3 = 2882.55 #kJ/kg\n",
+ "v3 = 0.3647 #m**3/kg\n",
+ "n = 1.3 #For superheated steam:\n",
+ "\n",
+ "#Calculations:\n",
+ "s1 = s2\n",
+ "s3 = s2\n",
+ "p2 = p1*(2/(n+1))**(n/(n-1)) #Pressue at state 2(in bar):\n",
+ "C2 = sqrt(2*(h1-h2)*10**3) #Velocity at throat(in m/s):\n",
+ "A2 = m2*v2/C2 #Cross-sectional area at throat(in m**2):\n",
+ "C3 = sqrt(2*(h1-h3)*10**3) #Ideal velocity at exit(in m/s):\n",
+ "A3 = m3*v3/C3a #Cross-sectional area at exit(in m**2): \n",
+ "r = C3a/C3 #Coefficient of velocity:\n",
+ "\n",
+ "#Results: \n",
+ "print \"Cross-sectional area at throat: \",round(A2*10**3,3),\" x 10^-3 m^2\"\n",
+ "print \"Cross-sectional area at exit: \",round(A3*10**3,3),\" x 10^-3 m^2\"\n",
+ "print \"Coefficient of velocity: \",round(r,3)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Cross-sectional area at throat: 3.21 x 10^-3 m^2\n",
+ "Cross-sectional area at exit: 3.647 x 10^-3 m^2\n",
+ "Coefficient of velocity: 0.875\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4, page no. 588"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "from __future__ import division\n",
+ "\n",
+ "from math import sqrt\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "p1 = 16 #Pressure of steam entering(in bar):\n",
+ "p3 = 5 #Pressure at exit(in bar):\n",
+ "m1 = 1 #Mass flow rate(in kg/s):\n",
+ "m2 = m1\n",
+ "m3 = m1\n",
+ "#From steam tables:\n",
+ "#For case 1:\n",
+ "h1 = 3034.8 #kJ/kg\n",
+ "s1 = 6.8844 #kJ/kg.K\n",
+ "v1 = 0.15862 #m**3/kg\n",
+ "n = 1.3\n",
+ "h2 = 2891.39 #kJ/kg\n",
+ "h3 = 2777 #kJ/kg\n",
+ "v2 = 0.2559 #m**3/kg\n",
+ "v3 = 0.3882 #m**3/kg\n",
+ "#For case 2:\n",
+ "h2a = 2905.73 #kJ/kg\n",
+ "v2a = 0.2598 #m**3/kg\n",
+ "v3a = 0.40023 #m**3/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "p2 = p1*(2/(n+1))**(n/(n-1)) #Pressure at the throat of nozzle(in bar):\n",
+ "q12 = h1-h2 #Heat drop up to throat section(in kJ/kg):\n",
+ "C2 = sqrt(2*(h1-h2)*10**3) #Velocity at throat(in m/s):\n",
+ "q23 = h2-h3 #Heat drop from exit(in kJ/kg):\n",
+ "C3 = sqrt(2*(h2-h3)*10**3+C2**2) #Velocity at exit(in m/s):\n",
+ "A2 = m2*v2/C2 #Throat area(in m**2):\n",
+ "A3 = m3*v3/C3 #Exit area(in m**2):\n",
+ "q12a = 0.9*q12 #Considering expansion to have 10% friction loss:\n",
+ "C2a = sqrt(2*q12a*10**3) #Actual velocity at throat(in m/s):\n",
+ "A2a = m2*v2a/C2a #Actual throat area(in m**2):\n",
+ "q23a = 0.9*q23 #Actual drop at the exit of the nozzle(in kJ/kg):\n",
+ "h3a = h2a-q23a #Actual enthalpy at state 3(in kJ/kg):\n",
+ "C3a = sqrt(2*q23a*10**3+C2a**2) #Actual velocity at exit(in m/s):\n",
+ "A3a = m3*v3a/C3a #Actual area at exit(in m**2):\n",
+ "\n",
+ "#Results:\n",
+ "print \"Throat area, without friction consideration: \",round(A2*10**4,2),\"cm**2\"\n",
+ "print \"Exit area, without friction consideration: \",round(A3*10**4,2),\"cm**2\"\n",
+ "print \"Throat area, with friction consideration: \",round(A2a*10**4,2),\"cm**2\"\n",
+ "print \"Exit area, with friction consideration: \",round(A3a*10**4,3),\"cm**2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Throat area, without friction consideration: 4.78 cm**2\n",
+ "Exit area, without friction consideration: 5.41 cm**2\n",
+ "Throat area, with friction consideration: 5.11 cm**2\n",
+ "Exit area, with friction consideration: 5.875 cm**2\n"
+ ]
+ }
+ ],
+ "prompt_number": 9
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5, page no. 590"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "from __future__ import division\n",
+ " \n",
+ "from math import sqrt,pi\n",
+ "#Variable Declaration: \n",
+ "P = 1 #Power of turbine(in MW):\n",
+ "p1 = 20 #Pressure of steam entering(in bar):\n",
+ "m = 8 #Steam consumption rate(in kg/kW.h):\n",
+ "p3 = 0.2 #Pressure at which steam leaves(in bar):\n",
+ "d = 0.01 #Throat diameter(in m):\n",
+ "#From Mollier diagram:\n",
+ "q12 = 142 #kJ/kg \n",
+ "v2 = 0.20 #m**3/kg\n",
+ "q13 = 807 #kJ/kg\n",
+ "v3 = 7.2 #m**3/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "C2 = sqrt(2*q12*10**3) #Velocity at throat(in m/s):\n",
+ "m2 = pi*d**2/4*C2/v2 #Mass flow rate:\n",
+ "m3 = m2\n",
+ "n = 10**3*m/(3600*m2) #Number of nozzles:\n",
+ "q13a = 0.90*q13 #Useful heat drop:\n",
+ "C3 = sqrt(2*10**3*q13a)#Velocity at exit(in m/s):\n",
+ "A3 = m3*v3/C3 #Area at exit(in m**2):\n",
+ "\n",
+ "#Results: \n",
+ "print \"Number of nozzles required: \",round(n) \n",
+ "print \"Area at exit: \",round(A3*10**4,2),\"cm**2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Number of nozzles required: 11.0\n",
+ "Area at exit: 12.5 cm**2\n"
+ ]
+ }
+ ],
+ "prompt_number": 8
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6, page no. 591"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "\n",
+ "from math import pi, atan, sqrt\n",
+ "#Variable Declaration: \n",
+ "p1 = 0.7 #Pressure at which steam is supplied(in MPa):\n",
+ "l = 0.06 #Length of diverging nozzle(in m):\n",
+ "d = 0.005 #Throat diameter(in mm):\n",
+ "p3 = 0.1 #Pressure at which steam leaves the nozzle(in MPa):\n",
+ "#From Mollier diagram:\n",
+ "q12 = 138 #kJ/kg\n",
+ "v2 = 0.58 #m**3/kg\n",
+ "T = 203 #\u00b0C\n",
+ "q23 = 247 #kJ/kg\n",
+ "q23a = 209.95 #kJ/kg\n",
+ "v3a = 1.7 #m**3/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "C2 = sqrt(2*q12*10**3) #Velocity at throat(in m/s):\n",
+ "m1 = pi*d**2/4*C2/v2 #Mass flow rate(in kg/s):\n",
+ "m2 = m1\n",
+ "m3 = m1\n",
+ "q = q12+q23a #Total heat drop(in kJ/kg):\n",
+ "C3 = sqrt(2*10**3*q) #Velocity at exit(in m/s):\n",
+ "A3 = m3*v3a/C3 #Area at exit(in m**2):\n",
+ "d1 = (sqrt(A3*4/pi))*10**3 #Diameter at exit(in mm):\n",
+ "a = atan((d1-d*10**3)/(2*60))*180/pi\n",
+ "\n",
+ "#Results: \n",
+ "print \"With no losses, temperature at throat: \",round(T,2),\"\u00b0C\"\n",
+ "print \"Velocity at throat: \",round(C2,2),\"m/s\"\n",
+ "print \"With losses, cone angle: \",round(2*a,2),\"\u00b0\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "With no losses, temperature at throat: 203.0 \u00b0C\n",
+ "Velocity at throat: 525.36 m/s\n",
+ "With losses, cone angle: 1.71 \u00b0\n"
+ ]
+ }
+ ],
+ "prompt_number": 7
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7, page no. 593"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "\n",
+ "from math import sin,pi,sqrt\n",
+ "\n",
+ "#Variable Declaration:\n",
+ "P = 5000 #Power of the turbine(in hp):\n",
+ "m = P*6/3600 #Steam required(in kg of steam/hp-hr):\n",
+ "n = 0.90 #Efficiency of nozzle:\n",
+ "a = 12 #Nozzle angle:\n",
+ "p = 5 #Pitch(in cm):\n",
+ "t = 0.3 #Thickness(in cm):\n",
+ "#From steam tables:\n",
+ "h1 = 2794 #kJ/kg\n",
+ "s1 = 6.4218 #kJ/kg.K\n",
+ "x2 = 0.9478\n",
+ "h2 = 2662.2 #kJ/kg\n",
+ "x2a = 0.9542\n",
+ "v2a = 0.2294 #m**3/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "s2 = s1\n",
+ "h12 = h1-h2 #Change in enthalpy(in kJ/kg):\n",
+ "h12a = n*h12 #Actual change(in kJ/kg):\n",
+ "C2 = sqrt(2*h12a*10**3) #Velocity at inlet(in m/s):\n",
+ "A2 = m*v2a/C2*10**4 #Area at exit of nozzle(in cm**2):\n",
+ "l = 60*pi/3 #Approximate length of the nozzle(in cm):\n",
+ "n = int(l/p)+1 #Number of nozzles:\n",
+ "l1 = n*p #Correct length of nozzle arc:\n",
+ "h = A2/((p*sin(a*pi/180)-t)*n)#Radial height of nozzle(in cm):\n",
+ "\n",
+ "#Results:\n",
+ "print \"Length of nozzle: \",round(l1,2),\"cm\"\n",
+ "print \"Radial height of nozzle: \",round(h,2),\"cm\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Length of nozzle: 65.0 cm\n",
+ "Radial height of nozzle: 4.08 cm\n"
+ ]
+ }
+ ],
+ "prompt_number": 10
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8, page no. 594"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "\n",
+ "from math import sqrt\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "p1 = 13 #Pressure at which steam enters(in bar):\n",
+ "p2 = 6 #Pressure at which steam leaves(in bar):\n",
+ "T1 = 150+273 #Temperature of steam entering(in K):\n",
+ "r = 1.4 #Adibatic insex of compression:\n",
+ "\n",
+ "#Calculations:\n",
+ "T2 = T1*(p2/p1)**((r-1)/r) #Final temperature of steam(in K):\n",
+ "C2 = sqrt(2*1.005*(T1-T2)) #Exit velocity(in m/s):\n",
+ "\n",
+ "#Results: \n",
+ "print \"Exit velocity: \",round(C2,2),\"m/s\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Exit velocity: 12.98 m/s\n"
+ ]
+ }
+ ],
+ "prompt_number": 12
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 9, page no. 595"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "\n",
+ "from math import sqrt\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "F = 350 #Force on the plate(in N):\n",
+ "p1 = 8 #Initial pressure(in bar):\n",
+ "p3 = 1 #Final pressure(in bar):\n",
+ "A2 = 5*10**(-4) #Throat cross-sectional area(in m**2):\n",
+ "#From steam tables:\n",
+ "h1 = 2769.1 #kJ/kg\n",
+ "s1 = 6.6628 #kJ/kg.K\n",
+ "x2 = 0.9717\n",
+ "h2 = 2685.17 #kJ/kg\n",
+ "v2 = 0.3932 #m**3/kg\n",
+ "x3 = 0.8238\n",
+ "h3 = 2277.6 #kJ/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "s2 = s1\n",
+ "s3 = s1\n",
+ "h12 = h1-h2 #Enthalpy change(in kJ/kg):\n",
+ "C2 = sqrt(2*h12*10**3) #Velocity at throat(in m/s):\n",
+ "m = A2*C2/v2 #Discharge at throat(in kg/s):\n",
+ "C3a = F/m #Actual exit velocity(in m/s):\n",
+ "h23 = h2-h3 #Theoretical enthalpy drop(in kJ/kg):\n",
+ "n = C3a**2/(2*h23*10**3) #Nozzle efficiency:\n",
+ "\n",
+ "#Results:\n",
+ "print \"Discharge at throat: \",round(m,3),\"kg/s\"\n",
+ "print \"Nozzle efficiency: \",round(n*100,2),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Discharge at throat: 0.521 kg/s\n",
+ "Nozzle efficiency: 55.37 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 10, page no. 597"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ " \n",
+ "from math import sqrt,pi\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "m = 5/60 #Mass flow rate(in kg/s):\n",
+ "p3 = 1 #Pressure at which steam is discharged(in bar):\n",
+ "p1 = 10 #Initial pressure(in bar):\n",
+ "T1 = 200+273 #Initial temperature(in K)\n",
+ "n = 1.3 #Adiabatic index of compression:\n",
+ "\n",
+ "#From steam tables:\n",
+ "h1 = 2827.9 #kJ/kg\n",
+ "s1 = 6.6940 #kJ/kg.K\n",
+ "v1 = 0.2060 #m**3/kg\n",
+ "h2a = 2711.23 #kJ/kg\n",
+ "s2a = 6.6749 #kJ/kg.K\n",
+ "h3 = 2420.08 #kJ/kg\n",
+ "v3 = 1.5025 #m**3/kg\n",
+ "psat = 3.44 #bar (at T = 138.18 \u00b0C)\n",
+ "Tsat = 155.12 #C (at p = 5.45 bar)\n",
+ "\n",
+ "#Calculations:\n",
+ "s3 = s2a\n",
+ "p2 = p1*(2/(n+1))**(n/(n-1)) #Pressure at throat(in bar):\n",
+ "C3 = sqrt(2*(h1-h3)*10**3) #Velocity at exit(in m/s):\n",
+ "A3 = m*v3/C3 #Exit area(in m**2):\n",
+ "d = sqrt(A3*4/pi) #Diameter of nozzle at exit(in m):\n",
+ "T2 = T1*(p2/p1)**((n-1)/n) #Temperature at throat(in K):\n",
+ "d1 = p2/psat #Degree of supersaturation:\n",
+ "u = Tsat-(T2-273) #Amount of undercooling(in \u00b0C):\n",
+ "\n",
+ "#Results: \n",
+ "print \"Degree of supersaturation: \",round(d1,2) \n",
+ "print \"Amount of undercooling: \",round(u,2),\"\u00b0C\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Degree of supersaturation: 1.59\n",
+ "Amount of undercooling: 16.82 \u00b0C\n"
+ ]
+ }
+ ],
+ "prompt_number": 16
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11, page no. 599"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "p1 = 4 #Initial pressure(in bar):\n",
+ "T1 = 180+273 #Initial temperature(in K):\n",
+ "p2 = 1.5 #Final pressure(in bar):\n",
+ "n = 1.3 #Index of compression:\n",
+ "nn = 0.95 #Efficiency due to heat loss:\n",
+ "C = 2.174 #Specific heat(in kJ/kg.K):\n",
+ "#From steam tables:\n",
+ "v1 = 0.5088 #m**3/kg \n",
+ "Tsat = 111.37+273 #K (at p = 1.5 bar)\n",
+ "\n",
+ "#Calculations:\n",
+ "h1 = p1*v1*10**2+2614 #Enthalpy at state 1(in kJ/kg):\n",
+ "v2 = v1*(p1/p2)**(1/n) #Specific volume at state 2(in m**3/kg):\n",
+ "h2 = p2*v2*10**2+2614 #Enthalpy at state 2(in kJ/kg):\n",
+ "dh = nn*(h1-h2) #Actual heat drop(in kJ/kg):\n",
+ "T2 = T1*(p2/p1)**((n-1)/n) #Temperature at state 2(in K):\n",
+ "dT = (1-nn)*(h1-h2)/C #Temperature rise due to supersaturation:\n",
+ "T2a = T2+dT #Actual temperature at state 2(in K):\n",
+ "u = Tsat-T2a #Amount of undercooling(in \u00b0C):\n",
+ "\n",
+ "#Results:\n",
+ "print \"Actual heat drop: \",round(dh,2),\"kJ/kg\"\n",
+ "print \"Amount of undercooling: \",round(u,2),\"\u00b0C\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Actual heat drop: 39.16 kJ/kg\n",
+ "Amount of undercooling: 22.18 \u00b0C\n"
+ ]
+ }
+ ],
+ "prompt_number": 19
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 12, page no. 600"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ " \n",
+ "from math import sqrt\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "p1 = 14 #Initial pressure(in bar):\n",
+ "T1 = 400+273 #Initial temperature(in K):\n",
+ "N = 16 #Number of nozzles:\n",
+ "p2 = 10 #Final pressure(in bar):\n",
+ "m = 5 #Discharge(in kg/s):\n",
+ "nn = 0.90 #Nozzle efficiency:\n",
+ "C1 = 100 #Inlet velocity(in m/s):\n",
+ "n = 1.3 #Insex of compression:\n",
+ "\n",
+ "#From steam tables:\n",
+ "h1 = 3257.5 #kJ/kg \n",
+ "s1 = 7.3026 #kJ/kg.K\n",
+ "T2 = 350.46 #\u00b0C\n",
+ "h2 = 3158.7 #kJ/kg\n",
+ "v2 = 0.2827 #m**3/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "h12 = (h1-h2)*nn #Actual enthalpy change(inn kJ/kg):\n",
+ "C2 = sqrt(2*h12*10**3) #Velocity at exit(in m/s):\n",
+ "A2 = m*v2/(C2*N)*10**4 #Cross-sectional area at exit(in cm**2):\n",
+ "C2a = sqrt(2*h12*10**3+C1**2)#Modified velocity at nozzle exit(in m/s):\n",
+ "ma = 16*2.13*433.41*10**(-4)/0.2827#ma = A2*C2a*N/v2*10**(-4) \n",
+ "p = (ma-m)/m*100 #% increase in discharge:\n",
+ "\n",
+ "#Results: \n",
+ "print \"Cross-sectional area at exit of nozzle: \",round(A2,2),\"cm**2\"\n",
+ "print \"Percentage increase in discharge: \",round(p,2),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Cross-sectional area at exit of nozzle: 2.09 cm**2\n",
+ "Percentage increase in discharge: 4.5 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 21
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 13, page no. 602"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "p1 = 20 #Initial pressure(in bar):\n",
+ "p3 = 5 #Final pressure(in bar):\n",
+ "n = 1.3\n",
+ "#From steam tables:\n",
+ "T1 = 212.42+273 #K\n",
+ "Tsat = 186.43+273 #K (at 11.6 bar)\n",
+ "psat = 5.452 #bar (at 155.14 \u00b0C)\n",
+ "h1 = 2799.5 #kJ/kg\n",
+ "v1 = 0.009963 #m**3/kg\n",
+ "s1 = 6.3409 #kJ/kg.K\n",
+ "h2aa = 2693.98 #kJ/kg\n",
+ "s2a = 6.5484 #kJ/kg.K\n",
+ "h3a = 2632.76 #kJ/kg\n",
+ "h3 = 2544.21 #kJ/kg\n",
+ "p2 = p1*0.58 #Pressure at throat(in bar):\n",
+ "\n",
+ "#Calculations:\n",
+ "s2aa = s1\n",
+ "s3a = s2a\n",
+ "s3 = s1\n",
+ "T2 = T1*(p2/p1)**((n-1)/n)\t#Temperature at state 2(in K):\n",
+ "d = p2/psat #Degree of supersaturation:\n",
+ "d1 = Tsat-T2 #Degree of undercooling:\n",
+ "h12 = (n/(n-1))*p1*10**2*v1*(1-(T2/T1)) #Isentropic enthalpy drop:\n",
+ "h2 = h1-h12 #Enthalpy at state 2(in kJ/kg):\n",
+ "h12aa = h1-h2aa #Heat drop with no saturation(in kJ/kg):\n",
+ "L = h12aa-h12 #Loss of available heat drop(in kJ/kg):\n",
+ "s12a = L/Tsat #Increase in entropy(in kJ/kg.K):\n",
+ "L1 = h3a-h3 #Loss due to undercooling(in kJ/kg):\n",
+ "p = L1/(h1-h3)*100 #Percentage loss:\n",
+ "\n",
+ "#Results: \n",
+ "print \"Degree of supersaturation: \",round(d,2)\n",
+ "print \"Degree of undercooling: \",round(d1,2),\"\u00b0C\"\n",
+ "print \"Entropy change: \",round(s12a,4),\"kJ/kg.K\"\n",
+ "print \"Loss due to undercooling: \",round(L1,2),\"kJ/kg\"\n",
+ "print \"Percentage loss: \",round(p,2)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Degree of supersaturation: 2.13\n",
+ "Degree of undercooling: 31.35 \u00b0C\n",
+ "Entropy change: 0.2075 kJ/kg.K\n",
+ "Loss due to undercooling: 88.55 kJ/kg\n",
+ "Percentage loss: 34.69\n"
+ ]
+ }
+ ],
+ "prompt_number": 24
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 14, page no. 604"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ " \n",
+ "from math import sqrt,pi\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "m1 = 150/60 #Mass flow rate(in kg/s):\n",
+ "H = 5 #Height of water level from the axis of injector(in m):\n",
+ "p4 = 20 #Pressuer at which steam is injected(in bar):\n",
+ "Z4 = 0.8 #Water level in boiler from the injector(in m):\n",
+ "x1 = 0.95 #Dryness fraction at state 1:\n",
+ "C4 = 20 #Velocity in delivery pipe(in m/s):\n",
+ "p3 = 1.013 #Atmospheric pressure(in bar):\n",
+ "d = 10**3 #Density(in kg/m**3):\n",
+ "g = 9.81 #Acceleration due to gravity(in m/s**2):\n",
+ "Cps = 3.18 #Specific heat of steam(in kJ/kg.K):\n",
+ "Cpw = 4.18 #Specific heat of water(in kJ/kg.K):\n",
+ "#From steam tables:\n",
+ "T1 = 212.42 #\u00b0C\n",
+ "Tw = 25 #\u00b0C\n",
+ "p2 = 0.7*p4\n",
+ "h1 = 2704.95 #kJ/kg\n",
+ "hfg1 = 1890.7 #kJ/kg\n",
+ "s1 = 6.1462 #kJ/kg.K\n",
+ "x2 = 0.923\n",
+ "h2 = 2639.10 #kJ/kg\n",
+ "v2 = 0.13 #m**3/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "s2 = s1\n",
+ "C2 = sqrt(2*(h1-h2)*10**3) #Velocity of steam at throat(in m/s):\n",
+ "C3 = sqrt(2*(g*Z4+p4*10**5/d+C4**2/2-p3*10**5/d))#Velocity at state 3(in m/s):\n",
+ "m = (C2-C3)/(sqrt(2*g*H)+C3) #Mass of water pumped per kg of steam(in kg):\n",
+ "m3 = m1+m1/m #Mass of mixture passing through state 3(in kg/s):\n",
+ "A3 = m3/(d*C3)*10**4 #Area of throat of mixing nozzle(in cm**2):\n",
+ "d3 = sqrt(A3*4/pi) #Diameter of throat of the mixing nozzle(in cm):\n",
+ "ms = m1/m #Mass of steam required for given flow rate(in kg/s):\n",
+ "A2 = ms*v2/C2*10**4 #Area at state 2(in cm**2):\n",
+ "d2 = sqrt(A2*4/pi) #Diameter of throat of steam nozzle(in cm):\n",
+ "T3 = (x1*hfg1+Cps*T1+m*Cpw*Tw)/(m*Cpw+Cps)#Temperature of water coming out of the injector(in C):\n",
+ "\n",
+ "#Results: \n",
+ "print \"Mass of water pumped per kg of steam: \",round(m,2),\"kg\"\n",
+ "print \"Diameter of throat of the mixing nozzle: \",round(d3,3),\"cm\"\n",
+ "print \"Diameter of throat of steam nozzle: \",round(d2,2),\"cm\"\n",
+ "print \"Temperature of water coming out of the injector: \",round(T3,2),\"\u00b0C\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Mass of water pumped per kg of steam: 3.98 kg\n",
+ "Diameter of throat of the mixing nozzle: 0.783 cm\n",
+ "Diameter of throat of steam nozzle: 1.69 cm\n",
+ "Temperature of water coming out of the injector: 145.63 \u00b0C\n"
+ ]
+ }
+ ],
+ "prompt_number": 27
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 15, page no. 607"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "\n",
+ "from math import sqrt\n",
+ "\n",
+ "#Variable Declaration: \n",
+ "p4 = 20 #Pressure at which steam is generated(in bar):\n",
+ "p1 = 1.5 #Pressure at inlet(in bar):\n",
+ "x1 = 0.9 #Dryness fraction:\n",
+ "M = 5000 #Mass of water taken from feed water tank(in kg/hr):\n",
+ "d = 10**3 #Density(in kg/m**3):\n",
+ "#From steam tables:\n",
+ "h1 = 2470.96 #kJ/kg \n",
+ "s1 = 6.6443 #kJ/kg.K\n",
+ "s2 = s1\n",
+ "x2 = 0.88\n",
+ "h2 = 2396.72 #kJ/kg\n",
+ "v2 = 1.7302 #m**3/kg\n",
+ "\n",
+ "#Calculations:\n",
+ "C2 = sqrt(2*(h1-h2)*10**3) #Steam velocity(in m/s):\n",
+ "C3 = sqrt(1.2*p4*2*10**5/d) #Velocity at 3(in m/s):\n",
+ "m = C2/C3-1 #Mass entrained per kg of steam:\n",
+ "ms = M/(3600*m) #Mass of steam supplied per second(in kg/s):\n",
+ "A2 = ms*v2/C2*10**4 #Area of steam nozzle(in cm**2):\n",
+ "D = M/3600+ms #Total discharge from injector(in kg/s):\n",
+ "A = D/(C3*d)*10**4 #Area of discharge orifice(in cm**2):\n",
+ "\n",
+ "#Results: \n",
+ "print \"Mass of water pumped per kg of steam: \",round(m,2),\"kg water/kg of steam\"\n",
+ "print \"Area of steam nozzle: \",round(A2,2),\"cm**2\"\n",
+ "print \"Area of discharge orifice: \",round(A,3),\"cm**2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Mass of water pumped per kg of steam: 4.56 kg water/kg of steam\n",
+ "Area of steam nozzle: 13.67 cm**2\n",
+ "Area of discharge orifice: 0.244 cm**2\n"
+ ]
+ }
+ ],
+ "prompt_number": 29
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+} \ No newline at end of file