diff options
Diffstat (limited to 'Power_Plant_Engineering_by_P._K._Nag/Ch8.ipynb')
| -rw-r--r-- | Power_Plant_Engineering_by_P._K._Nag/Ch8.ipynb | 282 |
1 files changed, 282 insertions, 0 deletions
diff --git a/Power_Plant_Engineering_by_P._K._Nag/Ch8.ipynb b/Power_Plant_Engineering_by_P._K._Nag/Ch8.ipynb new file mode 100644 index 00000000..37aa3970 --- /dev/null +++ b/Power_Plant_Engineering_by_P._K._Nag/Ch8.ipynb @@ -0,0 +1,282 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 8 : Condenser Feedwater and Circulating water systems" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Ex: 8.1 Pg: 591 " + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + " (a) The rate of flow of cooling water = 5817.7 kg/s \n", + " (b) The rate of air leakage into the condenser shell = 5.615 kg/s \n", + " (c) The length of tubes = 20.74 m \n", + " (d) The number of tubes = 7849 \n" + ] + } + ], + "source": [ + "from math import pi,log\n", + "#Input data\n", + "ws=250##The amount of steam received by the surface condenser in t/h\n", + "tsat=40##The saturated temperature in degree centigrade\n", + "m=12##The amount of moisture present in the steam in percentage\n", + "tc1=32##The inlet temperature of cooling water in degree centigrade\n", + "tc2=38##The outlet temperature of cooling water in degree centigrade\n", + "p=0.078##The pressure inside the condenser in bar\n", + "V=1.8##velocity of circulating water in m/s\n", + "do=0.0254##The outer diameter of the condenser tubes in m\n", + "T=0.00125##The thickness of the condenser tubes in m\n", + "pi=3.141##Mathematical constant of pi\n", + "U=2.6##The overall heat transfer coefficient in kW/m**2K\n", + "Cpc=4.187##The specific heat of water in kJ/kgK\n", + "R=0.287##Real gas constant in kJ/kgK\n", + "P=1000##The density of water in kg/m**3\n", + "\n", + "#Calculations\n", + "x2=0.88##The quality of the steam\n", + "hfg=2407##The enthalpy of evaporation at 40 degreeC in kJ/kg\n", + "h=x2*hfg##The change in enthalpy in kJ/kg\n", + "di=do-(2*T)##The inner diameter of the condenser tubes in m\n", + "wc=(((ws*1000)/3600)*h)/(Cpc*(tc2-tc1))##Mass flow rate of water in kg/s\n", + "psat=0.07375##The saturated pressure at 40 degree centigrade in bar\n", + "pair=(p-psat)*100##The pressure of air in kPa\n", + "vf=0.001008##Specific volume of saturated liquid in m**3/kg\n", + "vfg=19.544##Specific volume of vapour in m**3/kg\n", + "v2=vf+(x2*vfg)##Specific volume at 40 degree centigrade in m**3/kg\n", + "wair=(pair*((ws*1000)/3600)*v2)/(R*(tsat+273))##Mass flow rate of air in kg/s\n", + "t1m=(8-2)/log(8/2)##The mean temperature in degree centigrade\n", + "Ao=(((ws*1000)/3600)*h)/(U*t1m)##The area of the tubes in m**2\n", + "n=(wc*(4/pi)*(1/di**2)*(1/(P*V)))##The number of tubes\n", + "l=Ao/(pi*do*n)##Yhe length of tubes in m\n", + "\n", + "#Output\n", + "print \" (a) The rate of flow of cooling water = %3.1f kg/s \\n (b) The rate of air leakage into the condenser shell = %3.3f kg/s \\n (c) The length of tubes = %3.2f m \\n (d) The number of tubes = %3.0f \"%(wc,wair,l,n)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Ex: 8.2 Pg: 592" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + " (a) The rate of saving of condensate and the rate of saving in the heat supply in the boiler due to seperate air extraction pump = 1.105 kW \n", + " (b) The percentage reduction in air ejector load due to this seperate air extraction method = 80.8 percent \n", + " (c) The rate of cooling water flow = 533 kg/s \n" + ] + } + ], + "source": [ + "#Input data\n", + "ws=20##The amount of dry saturated steam received by a surface condencer in t/h\n", + "tsh=40##The temperature of dry saturated steam in degree centigrade \n", + "wa=(0.35/1000)##The mass flow rate of air per 1000 kg of steam in kg\n", + "tc=38##The temperature at which condensate leaves the temperature in degree centigrade\n", + "tm=10##The temperature at which makeup water is supplied in degree centigrade\n", + "tc1=32##The inlet temperature of cooling water in degree centigrade\n", + "tc2=38##The outlet temperature of cooling water in degree centigrade\n", + "tas=27##The temperature of air along with steam in degree centigrade\n", + "psat=0.07384##The saturated pressure at 40 degree C in bar\n", + "vg=19.52##The specific volume at 40 degree C in m**3/kg\n", + "R=0.287##Real gas constant in kJ/kgK\n", + "Cpc=4.187##The specific heat of water in kJ/kgK\n", + "Cp=1.005##The specific heat of air in kJ/kgK\n", + "\n", + "#Calculations\n", + "pair=((wa*R*(tsh+273)*1000)/vg)*10**-5##The pressure of air in bar\n", + "psat1=0.06624##The saturated pressure at 38 degree C in bar\n", + "vg1=21.63##The specific volume at 38 degree C in m**3/kg\n", + "pair1=psat-psat1##The pressure of air in bar\n", + "wa1=(ws*1000)*wa##Mass of air removed per hour in kg/h\n", + "V1=((wa1*R*(273+tc2)*1000))/(pair1*10**5)##Volume of air remove per hour\n", + "ws1=V1/vg1##The mass of steam accompanying air in kg/h\n", + "psat2=0.03564##The saturated pressure at 27 degree C in bar\n", + "vg2=38.81##The specific volume at 27 degree C in m**3/kg\n", + "pair2=psat-psat2##The pressure of air in bar\n", + "V2=(wa1*R*1000*(tas+273))/(pair2*10**5)##Volume of air removed per hr in m**3/hr\n", + "ws2=V2/vg2##The mass of steam accompanying air in kg/h\n", + "ws3=ws1-ws2##Saving mass of steam by using seperate extraction in kg/hr\n", + "Q3=(ws3*Cpc*(tc-tm))/3600##Saving in heat supply in the boiler in kW\n", + "V=((V1-V2)/V1)*100##Percentage reduction in air ejector load in %\n", + "hc=159.3##Enthalpy at 38 degree C in kJ/kg\n", + "hs1=2574.3##Enthalpy at 40 degree C in kJ/kg\n", + "hs2=2550.3##Enthalpy at 27 degree C in kJ/kg\n", + "Q2=(((ws*1000)*(hs1-hc))-((wa1*(Cp*(tsh-tas))))-(ws3*hs2))/3600##The amount of heat in kW\n", + "wc=Q2/(Cpc*(tc2-tc1))##The mass flow rate of water in kg/s\n", + "\n", + "#Output\n", + "print \" (a) The rate of saving of condensate and the rate of saving in the heat supply in the boiler due to seperate air extraction pump = %3.3f kW \\n (b) The percentage reduction in air ejector load due to this seperate air extraction method = %3.1f percent \\n (c) The rate of cooling water flow = %3.0f kg/s \"%(Q3,V,wc)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Ex: 8.3 Pg: 593" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + " (a) The temperature of water leaving the tower = 21.8 degree centigrade \n", + " (b) The fraction of water evaporated = 0.0125 kg/kg dry air \n", + " (c) The approach of the cooling tower = 6.6 degree centigrade \n", + " The Range of the cooling tower = 8.2 degree centigrade \n" + ] + } + ], + "source": [ + "#Input data\n", + "tw3=30##The inlet temperature of water in degree centigrade\n", + "wc=1.15##Mass flow rate of cooling water in kg per kg air\n", + "tdb1=20##The dry bulb temperature of air in degree centigrade\n", + "R1=60##Relative humidity of air while entering in percentage\n", + "tdb2=28##The dry bulb temperature while leaving in degree centigrade\n", + "R2=90##Relative humidity of air while leaving in percentage\n", + "tm=20##The temperature of makeup water in degree centigrade\n", + "Cpc=4.187##The specific heat of water in kJ/kgK\n", + "G=1##Mass flow rate of dry air in kg/s\n", + "\n", + "#Calculations\n", + "twb1=15.2## from psychrometric chart The wet bulb temperature while entering in degree centigrade\n", + "twb2=26.7## from psychrometric chart The wet bulb temperature while leaving in degree centigrade\n", + "h1=43##The enthalpy from chart for dry air in kJ/kg dry air\n", + "h2=83.5##The enthalpy from chart in kJ/kg dry air\n", + "W1=0.0088##Humidity in kg water vapour/kg dry air\n", + "W2=0.0213##Humidity in kg water vapour/kg dry air\n", + "hw3=125.8##Enthalpy of water entering the tower in kJ/kg\n", + "hw=84##Enthalpy of makeup water in kJ/kg\n", + "hwc=((G/wc)*((h2-h1)-(W2-W1)*hw))##The change in enthalpy of water in kJ/kg\n", + "tw4=tw3-(hwc/Cpc)##The exit temperature of water in degree centigrade\n", + "ta=tw4-twb1##The approach temperature in degree centigrade\n", + "tr=tw3-tw4##The range temperature in degree centigrade\n", + "x=G*(W2-W1)##Fraction of water evaporated in kg/kg dry air\n", + "\n", + "#Output\n", + "print \" (a) The temperature of water leaving the tower = %3.1f degree centigrade \\n (b) The fraction of water evaporated = %3.4f kg/kg dry air \\n (c) The approach of the cooling tower = %3.1f degree centigrade \\n The Range of the cooling tower = %3.1f degree centigrade \"%(tw4,x,ta,tr)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Ex: 8.4 Pg: 593" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + " (a) The final temperature of the water = 16.60 degree centigrade \n", + " (b)The amount of makeup water required per second = 0.2605 kg/s \n" + ] + } + ], + "source": [ + "#Input data\n", + "tw3=45##The temperature of warm water in degree centigrade\n", + "wc1=6##The cooling water inflow in kg/s\n", + "V=10##volume flow of ID fan in m**3/s\n", + "Ws=4.90##Heat absorbed by air in kW\n", + "ti=20##The temperature of air entering the tower in degree centigrade\n", + "R=60##The relative humidity in percentage\n", + "to=26##The temperature of air leaving the tower in degree centigrade\n", + "p=1.013##The constant pressure throughout the tower in bar\n", + "r=0.287##Real gas constant in kJ/kgK\n", + "Cpc=4.187##The specific heat of water in kJ/kgK\n", + "Cp=1.005##The specific heat of air in kJ/kgK\n", + "\n", + "#Calculations\n", + "ps=0.0234##The pressure at 20 degreec in bar\n", + "ps1=(R/100)*ps##The pressure of water vapour in bar\n", + "pa1=p-ps1##The pressure of air in bar\n", + "G1=(pa1)/(r*10**-3*(ti+273))##The mass flow rate of dry air in kg/s\n", + "w1=(ps1*10**5*V)/(0.4619*10**3*(ti+273))##Mass flow rate of vopour in kg/s\n", + "W1=w1/G1##Moisture flow in kg vap/kg dry air\n", + "ps2=0.0336##The pressure at 26 degree C at exit in bar\n", + "pw2=0.0336##The pressure of water vapour at 26 degree C at exit in bar\n", + "W2=(0.622)*(pw2/(1-pw2))##oisture flow in kg vap/kg dry air\n", + "G2=G1##The mass flow rate of dry air in kg/s\n", + "w2=W2*G2##Moisture flow at exit in kg/s\n", + "wm=w2-w1##Makeup water required in kg/s\n", + "wc2=wc1-wm##The cooling water outflow in kg/s\n", + "hw3=Cpc*tw3##The enthalpy of warm water in kJ/kg\n", + "hg=2538.1##The enthalpy of gas at 20 degree C in kJ/kg\n", + "tsat=12##The saturation temperature in degree centigrade\n", + "pw1=0.01404##The pressure at 12 degree C in bar\n", + "hw1=hg+(1.88*(ti-tsat))##The enthalpy of warm water in kJ/kg\n", + "hw2=2548.4##The enthalpy of evaporation at 26 degree C in kJ/kg\n", + "hw4=((G1*((Cp*(to-ti))+W2*hw2-W1*hw1)-Ws)-(wc1*hw3))/-wc2##The enthapy of warm water at point 4 in kJ/kg\n", + "E=hw4/Cpc##Exit water temperature in degree centigrade\n", + "\n", + "#Output\n", + "print \" (a) The final temperature of the water = %3.2f degree centigrade \\n (b)The amount of makeup water required per second = %3.4f kg/s \"%(E,wm)" + ] + } + ], + "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 +} |