{
"cells": [
 {
		   "cell_type": "markdown",
	   "metadata": {},
	   "source": [
       "# Chapter 8: Vapour Power Cycles"
	   ]
	},
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.10: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_10.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 10')\n",
"disp('at inlet to HP turbine,h2=3287.1 KJ/kg,s2=6.6327 KJ/kg K')\n",
"h2=3287.1;\n",
"s2=6.6327;\n",
"disp('By interpolation state 3 i.e. for isentropic expansion betweeen 2-3 lies at 328.98oc at 30 bar.h3=3049.48 KJ/kg')\n",
"h3=3049.48;\n",
"disp('actual enthapy at 3_a,h3_a=h2-0.80*(h2-h3)in KJ/kg')\n",
"h3_a=h2-0.80*(h2-h3)\n",
"disp('enthalpy at inlet to LP turbine,h4=3230.9 KJ/kg,s4=6.9212 KJ K')\n",
"h4=3230.9;\n",
"s4=6.9212;\n",
"disp('for ideal expansion from 4-6,s4=s6.Let dryness fraction at state 6 be x6.')\n",
"s6=s4;\n",
"disp('s6=6.9212=sf at 0.075 bar+x6* sfg at 0.075 bar in KJ/kg K')\n",
"disp('from steam tables,at 0.075 bar,sf=0.5764 KJ/kg K,sfg=7.6750 KJ/kg K')\n",
"sf=0.5764;\n",
"sfg=7.6750;\n",
"disp('so x6=(s6-sf)/sfg')\n",
"x6=(s6-sf)/sfg\n",
"x6=0.827;//approx.\n",
"disp('h6=hf at 0.075 bar+x6*hfg at 0.075 bar in KJ/kg K')\n",
"disp('from steam tables,at 0.075 bar,hf=168.79 KJ/kg,hfg=2406.0 KJ/kg')\n",
"hf=168.79;\n",
"hfg=2406.0;\n",
"h6=hf+x6*hfg\n",
"disp('for actual expansion process in LP turbine.')\n",
"disp('h6_a=h4-0.85*(h4-h6) in KJ/kg')\n",
"h6_a=h4-0.85*(h4-h6)\n",
"disp('Ideally,enthalpy at bleed point can be obtained by locating state 5 using s5=s4.The pressure at bleed point shall be saturation pressure corresponding to the 140oc i.e from steam tables.Let dryness fraction at state 5 be x5.')\n",
"p5=3.61;\n",
"s5=s4;\n",
"disp('s5_a=6.9212=sf at 140oc+x5*sfg at 140oc')\n",
"disp('from steam tables,at 140oc,sf=1.7391 KJ/kg K,sfg=5.1908 KJ/kg K')\n",
"sf=1.7391;\n",
"sfg=5.1908;\n",
"disp('so x5=(s5-sf)/sfg')\n",
"x5=(s5-sf)/sfg\n",
"x5=0.99;//approx.\n",
"disp('h5=hf at 140oc+x5*hfg at 140oc in KJ/kg')\n",
"disp('from steam tables,at 140oc,hf=589.13 KJ/kg,hfg=2144.7 KJ/kg')\n",
"hf=589.13;\n",
"hfg=2144.7;\n",
"h5=hf+x5*hfg\n",
"disp('actual enthalpy,h5_a=h4-0.85*(h4-h5)in KJ/kg')\n",
"h5_a=h4-0.85*(h4-h5)\n",
"disp('enthalpy at exit of open feed water heater,h9=hf at 30 bar=1008.42 KJ/kg')\n",
"h9=1008.42;\n",
"disp('specific volume at inlet of CEP,v7=0.001008 m^3/kg')\n",
"v7=0.001008;\n",
"disp('enthalpy at inlet of CEP,h7=168.79 KJ/kg')\n",
"h7=168.79;\n",
"disp('for pumping process 7-8,h8=h7+v7*(3.61-0.075)*10^2 in KJ/kg')\n",
"h8=h7+v7*(3.61-0.075)*10^2\n",
"disp('Applying energy balance at open feed water heater.Let mass of bled steam be m kg per kg of steam generated.')\n",
"disp('m*h5+(1-m)*h8=h9')\n",
"disp('so m=(h9-h8)/(h5-h8) in kg /kg of steam generated')\n",
"m=(h9-h8)/(h5-h8)\n",
"disp('For process on feed pump,9-1,v9=vf at 140oc=0.00108 m^3/kg')\n",
"v9=0.00108;\n",
"disp('h1=h9+v9*(70-3.61)*10^2 in KJ/kg')\n",
"h1=h9+v9*(70-3.61)*10^2 \n",
"disp('Net work per kg of steam generated,W_net=(h2-h3_a)+(h4-h5_a)+(1-m)*(h5_a-h6_a)-{(1-m)*(h8-h7)+(h1-h9)}in KJ/kg steam generated')\n",
"W_net=(h2-h3_a)+(h4-h5_a)+(1-m)*(h5_a-h6_a)-{(1-m)*(h8-h7)+(h1-h9)}\n",
"disp('heat added per kg of steam generated,q_add=(h2-h1)+(h4-h3_a)in KJ/kg of steam generated')\n",
"q_add=(h2-h1)+(h4-h3_a)\n",
"disp('Thermal efficiency,n=W_net/q_add')\n",
"n=W_net/q_add\n",
"disp('in percentage')\n",
"n=n*100\n",
"disp('so thermal efficiency=39.03%%')\n",
"disp('NOTE=>In this question there is some calculation mistake while calculating W_net and q_add in book, which is corrected above and the answers may vary.')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.11: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_11.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 11')\n",
"disp('Enthalpy of steam entering ST1,h2=3308.6 KJ/kg,s2=6.3443 KJ/kg K')\n",
"h2=3308.6;\n",
"s2=6.3443;\n",
"disp('for isentropic expansion 2-3-4-5,s2=s3=s4=s5')\n",
"s3=s2;\n",
"s4=s3;\n",
"s5=s4;\n",
"disp('Let dryness fraction of states 3,4 and 5 be x3,x4 and x5')\n",
"disp('s3=6.3443=sf at 10 bar+x3*sfg at 10 bar')\n",
"disp('so x3=(s3-sf)/sfg')\n",
"disp('from steam tables,at 10 bar,sf=2.1387 KJ/kg K,sfg=4.4478 KJ/kg K')\n",
"sf=2.1387;\n",
"sfg=4.4478;\n",
"x3=(s3-sf)/sfg\n",
"x3=0.945;//approx.\n",
"disp('h3=hf+x3*hfg in KJ/kg')\n",
"disp('from steam tables,hf=762.81 KJ/kg,hfg=2015.3 KJ/kg')\n",
"hf=762.81;\n",
"hfg=2015.3;\n",
"h3=hf+x3*hfg\n",
"disp('s4=6.3443=sf at 1.5 bar+x4*sfg at 1.5 bar')\n",
"disp('so x4=(s4-sf)/sfg')\n",
"disp('from steam tables,at 1.5 bar,sf=1.4336 KJ/kg K,sfg=5.7897 KJ/kg K')\n",
"sf=1.4336;\n",
"sfg=5.7897;\n",
"x4=(s4-sf)/sfg\n",
"x4=0.848;//approx.\n",
"disp('so h4=hf+x4*hfg in KJ/kg')\n",
"disp('from steam tables,at 1.5 bar,hf=467.11 KJ/kg,hfg=2226.5 KJ/kg')\n",
"hf=467.11;\n",
"hfg=2226.5;\n",
"h4=hf+x4*hfg\n",
"disp('s5=6.3443=sf at 0.05 bar+x5*sfg at 0.05 bar')\n",
"disp('so x5=(s5-sf)/sfg')\n",
"disp('from steam tables,at 0.05 bar,sf=0.4764 KJ/kg K,sfg=7.9187 KJ/kg K')\n",
"sf=0.4764;\n",
"sfg=7.9187;\n",
"x5=(s5-sf)/sfg\n",
"x5=0.739;//approx.\n",
"disp('h5=hf+x5*hfg in KJ/kg')\n",
"disp('from steam tables,at 0.05 bar,hf=137.82 KJ/kg,hfg=2423.7 KJ/kg')\n",
"hf=137.82;\n",
"hfg=2423.7;\n",
"h5=hf+x5*hfg \n",
"disp('h6=hf at 0.05 bar=137.82 KJ/kg')\n",
"h6=137.82;\n",
"disp('v6=vf at 0.05 bar=0.001005 m^3/kg')\n",
"v6=0.001005;\n",
"disp('h7=h6+v6*(1.5-0.05)*10^2 in KJ/kg')\n",
"h7=h6+v6*(1.5-0.05)*10^2\n",
"disp('h8=hf at 1.5 bar=467.11 KJ/kg')\n",
"h8=467.11; \n",
"disp('v8=0.001053 m^3/kg=vf at 1.5 bar')\n",
"v8=0.001053;\n",
"disp('h9=h8+v8*(150-1.5)*10^2 in KJ/kg')\n",
"h9=h8+v8*(150-1.5)*10^2\n",
"disp('h10=hf at 150 bar=1610.5 KJ/kg')\n",
"h10=1610.5; \n",
"disp('v10=0.001658 m^3/kg=vf at 150 bar')\n",
"v10=0.001658;\n",
"disp('h12=h10+v10*(150-10)*10^2 in KJ/kg')\n",
"h12=h10+v10*(150-10)*10^2\n",
"disp('Let mass of steam bled out at 10 bar,1.5 bar be m1 and m2 per kg of steam generated.')\n",
"disp('Heat balance on closed  feed water heater yields,')\n",
"disp('m1*h3+(1-m)*h9=m1*h10+(1-m1)*4.18*150')\n",
"disp('so m1=(4.18*150-h9)/(h3-h9-h10+4.18*150)in kg/kg of steam generated.')\n",
"m1=(4.18*150-h9)/(h3-h9-h10+4.18*150)\n",
"disp('heat balance on open feed water can be given as under,')\n",
"disp('m2*h4+(1-m1-m2)*h7=(1-m1)*h8')\n",
"disp('so m2=((1-m1)*(h8-h7))/(h4-h7)in kg/kg of steam')\n",
"m2=((1-m1)*(h8-h7))/(h4-h7)\n",
"disp('for mass flow rate of 300 kg/s=>m1=36 kg/s,m2=39 kg/s')\n",
"disp('For mixing after closed feed water heater,')\n",
"disp('h1=(4.18*150)*(1-m1)+m1*h12 in KJ/kg')\n",
"h1=(4.18*150)*(1-m1)+m1*h12\n",
"disp('Net work output per kg of steam generated=W_ST1+W_ST2+W_ST3-{W_CEP+W_FP+W_FP2}')\n",
"disp('W_net=(h2-h3)+(1-m1)*(h3-h4)+(1-m1-m2)*(h4-h5)-{(1-m1-m2)*(h7-h6)+(1-m1)*(h9-h8)+(m1*(h12-h10))}in KJ/kg of steam generated.')\n",
"W_net=(h2-h3)+(1-m1)*(h3-h4)+(1-m1-m2)*(h4-h5)-{(1-m1-m2)*(h7-h6)+(1-m1)*(h9-h8)+(m1*(h12-h10))}\n",
"disp('heat added per kg of steam generated,q_add=(h2-h1) in KJ/kg')\n",
"q_add=(h2-h1)\n",
"disp('cycle thermal efficiency,n=W_net/q_add')\n",
"n=W_net/q_add\n",
"disp('in percentage')\n",
"n=n*100\n",
"disp('Net power developed in KW=1219*300 in KW')\n",
"1219*300\n",
"disp('cycle thermal efficiency=47.6%')\n",
"disp('Net power developed=365700 KW')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.12: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_12.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 12')\n",
"P=100*10^3;//net power output in KW\n",
"disp('At inlet to HPT,h2=3373.7 KJ/kg,s2=6.5966 KJ/kg K')\n",
"h2=3373.7;\n",
"s2=6.5966;\n",
"disp('For isentropic expansion between 2-3-4-5,s2=s3=s4=s5')\n",
"s3=s2;\n",
"s4=s3;\n",
"s5=s4;\n",
"disp('state 3 lies in superheated region as s3>sg at 20 bar.By interpolation from superheated steam table,T3=261.6oc.Enthalpy at 3,h3=2930.57 KJ/kg')\n",
"T3=261.6;\n",
"h3=2930.57;\n",
"disp('since s4<sg at 4 bar so state 4 and 5 lies in wet region.')\n",
"disp('Let dryness fraction at state 4 ans 5 be x4 and x5.')\n",
"disp('s4=6.5966=sf at 4 bar+x4*sfg at 4 bar')\n",
"disp('from steam tables,at 4 bar,sf=1.7766 KJ/kg K,sfg=5.1193 KJ/kg K')\n",
"sf=1.7766;\n",
"sfg=5.1193;\n",
"disp('x4=(s4-sf)/sfg')\n",
"x4=(s4-sf)/sfg\n",
"x4=0.941;//approx.\n",
"disp('h4=hf at 4 bar+x4*hfg at 4 bar in KJ/kg')\n",
"disp('from steam tables,at 4 bar,hf=604.74 KJ/kg,hfg=2133.8 KJ/kg')\n",
"hf=604.74;\n",
"hfg=2133.8;\n",
"h4=hf+x4*hfg\n",
"disp('for state 5,')\n",
"disp('s5=6.5966=sf at 0.075 bar+x5*sfg at 0.075 bar')\n",
"disp('from steam tables,at 0.075 bar,sf=0.5764 KJ/kg K,sfg=7.6750 KJ/kg K')\n",
"sf=0.5764;\n",
"sfg=7.6750;\n",
"disp('x5=(s5-sf)/sfg')\n",
"x5=(s5-sf)/sfg\n",
"x5=0.784;//approx.\n",
"disp('h5=hf at 0.075 bar+x5*hfg at 0.075 bar in KJ/kg')\n",
"disp('from steam tables,at 0.075 bar,hf=168.76 KJ/kg,hfg=2406.0 KJ/kg')\n",
"hf=168.76;\n",
"hfg=2406.0;\n",
"h5=hf+x5*hfg\n",
"disp('Let mass of steam bled at 20 bar be m1 and m2 per kg of steam generated.')\n",
"disp('h10=hf at 20 bar=908.76 KJ/kg,h8=hf at 4 bar=604.74 KJ/kg')\n",
"h10=908.76;\n",
"h8=604.74;\n",
"disp('At trap h10=h11=908.79 KJ/kg')\n",
"h11=h10;\n",
"disp('At condensate extraction pump,(CEP),h7-h6=v6*(4-0.075)*10^2 in KJ/kg')\n",
"disp('here v6=vf at 0.075 bar=0.001008 m^3/kg,h6=hf at 0.075 bar=168.79 KJ/kg')\n",
"v6=0.001008;\n",
"h6=168.79;\n",
"disp('so h7=h6+v6*(4-0.075)*10^2 in KJ/kg')\n",
"h7=h6+v6*(4-0.075)*10^2\n",
"disp('At feed pump,(FP),h9-h8=v8*(20-4)*10^2 in KJ/kg')\n",
"disp('here v8=vf at 4 bar=0.001084 m^3/kg,h8=hf at 4 bar=604.74 KJ/kg')\n",
"v8=0.001084;\n",
"h8=604.74;\n",
"disp('so h9=h8+v8*(20-4)*10^2 in KJ/kg')\n",
"h9=h8+v8*(20-4)*10^2\n",
"disp('Let us apply heat balance at closed feed water heater,')\n",
"disp('m1*h3+h9=m1*h10+4.18*200')\n",
"disp('so m1=(4.18*200-h9)/(h3-h10) in kg')\n",
"m1=(4.18*200-h9)/(h3-h10)\n",
"m1=0.114;//approx.\n",
"disp('Applying heat balance at open feed water,')\n",
"disp('m1*h11+m2*h4+(1-m1-m2)*h7=h8')\n",
"disp('so m2=(h8-m1*h11-h7+m1*h7)/(h4-h7) in kg')\n",
"m2=(h8-m1*h11-h7+m1*h7)/(h4-h7)\n",
"m2=0.144;//approx.\n",
"disp('Net work per kg steam generated,')\n",
"disp('w_net=(h2-h3)+(1-m1)*(h3-h4)+(1-m1-m2)*(h4-h5)-{(1-m1-m2)*(h7-h6)+(h9-h8)} in KJ/kg')\n",
"w_net=(h2-h3)+(1-m1)*(h3-h4)+(1-m1-m2)*(h4-h5)-{(1-m1-m2)*(h7-h6)+(h9-h8)} \n",
"disp('Heat added per kg steam generated,q_add=(h2-h1) in KJ/kg')\n",
"h1=4.18*200;\n",
"q_add=(h2-h1)\n",
"disp('Thermal efficiency=w_net/q_add')\n",
"w_net/q_add\n",
"disp('in percentage')\n",
"w_net*100/q_add\n",
"disp('steam generation rate=P/w_net in kg/s')\n",
"P/w_net\n",
"disp('so thermal efficiency=44.78%')\n",
"disp('steam generation rate=87.99 kg/s')\n",
"disp('a> For the reheating introduced at 20 bar up to 400oc.The modified cycle representation is shown on T-S diagram by 1-2-3-3_a-4_a-5_a-6-7-8-9-10-11')\n",
"disp('At state 2,h2=3373.7 KJ/kg,s2=6.5966 KJ/kg K')\n",
"h2=3373.7;\n",
"s2=6.5966;\n",
"disp('At state 3,h3=2930.57 KJ/kg')\n",
"h3=2930.57;\n",
"disp('At state 3_a,h3_a=3247.6 KJ/kg,s3_a=7.1271 KJ/kg K')\n",
"h3_a=3247.6;\n",
"s3_a=7.1271;\n",
"disp('At state 4_a and 5_a,s3_a=s4_a=s5_a=7.1271 KJ/kg K')\n",
"s4_a=s3_a;\n",
"s5_a=s4_a;\n",
"disp('From steam tables by interpolation state 4_a is seen to be at 190.96oc at 4 bar,h4_a=2841.02 KJ/kg')\n",
"h4_a=2841.02;\n",
"disp('Let dryness fraction at state 5_a be x5,')\n",
"disp('s5_a=7.1271=sf at 0.075 bar+x5_a*sfg at 0.075 bar')\n",
"disp('from steam tables,at 0.075 bar,sf=0.5764 KJ/kg K,sfg=7.6750 KJ/kg K')\n",
"disp('so x5_a=(s5_a-sf)/sfg')\n",
"x5_a=(s5_a-sf)/sfg\n",
"x5_a=0.853;//approx.\n",
"disp('h5_a=hf at 0.075 bar+x5_a*hfg at 0.075 bar in KJ/kg')\n",
"disp('from steam tables,at 0.075 bar,hf=168.76 KJ/kg,hfg=2406.0 KJ/kg')\n",
"h5_a=hf+x5_a*hfg\n",
"disp('Let mass of bled steam at 20 bar and 4 bar be m1_a,m2_a per kg of steam generated.Applying heat balance at closed feed water heater.')\n",
"disp('m1_a*h3+h9=m1*h10+4.18*200')\n",
"disp('so m1_a=(4.18*200-h9)/(h3-h10) in kg')\n",
"m1_a=(4.18*200-h9)/(h3-h10)\n",
"m1_a=0.114;//approx.\n",
"disp('Applying heat balance at open feed water heater,')\n",
"disp('m1_a*h11+m2_a*h4_a+(1-m1_a-m2_a)*h7=h8')\n",
"disp('so m2_a=(h8-m1_a*h11-h7+m1_a*h7)/(h4_a-h7) in kg')\n",
"m2_a=(h8-m1_a*h11-h7+m1_a*h7)/(h4_a-h7)\n",
"m2_a=0.131;//approx.\n",
"disp('Net work per kg steam generated')\n",
"disp('w_net=(h2-h3)+(1-m1_a)*(h3_a-h4_a)+(1-m1_a-m2_a)*(h4_a-h5_a)-{(1-m1_a-m2_a)*(h7-h6)+(h9-h8)}in KJ/kg')\n",
"w_net=(h2-h3)+(1-m1_a)*(h3_a-h4_a)+(1-m1_a-m2_a)*(h4_a-h5_a)-{(1-m1_a-m2_a)*(h7-h6)+(h9-h8)}\n",
"disp('Heat added per kg steam generated,q_add=(h2-h1)+(1-m1_a)*(h3_a-h3)in KJ/kg')\n",
"q_add=(h2-h1)+(1-m1_a)*(h3_a-h3)\n",
"disp('Thermal efficiency,n=w_net/q_add')\n",
"n=w_net/q_add\n",
"disp('in percentage')\n",
"n=n*100\n",
"disp('% increase in thermal efficiency due to reheating=(0.4503-0.4478)*100/0.4478')\n",
"(0.4503-0.4478)*100/0.4478\n",
"disp('so thermal efficiency of reheat cycle=45.03%')\n",
"disp('% increase in efficiency due to reheating=0.56%')\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
""
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.13: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_13.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 13')\n",
"disp('For mercury cycle,')\n",
"disp('insentropic heat drop=349-234.5 in KJ/kg Hg')\n",
"349-234.5\n",
"disp('actual heat drop=0.85*114.5 in KJ/kg Hg')\n",
"0.85*114.5\n",
"disp('Heat rejected in condenser=(349-97.325-35)in KJ/kg')\n",
"(349-97.325-35)\n",
"disp('heat added in boiler=349-35 in KJ/kg')\n",
"349-35\n",
"disp('For steam cycle,')\n",
"disp('Enthalpy of steam generated=h at 40 bar,0.98 dry=2767.13 KJ/kg')\n",
"h=2767.13;\n",
"disp('Enthalpy of inlet to steam turbine,h2=h at 40 bar,450oc=3330.3 KJ/kg')\n",
"h2=3330.3;\n",
"disp('Entropy of steam at inlet to steam turbine,s2=6.9363 KJ/kg K')\n",
"s2=6.9363;\n",
"disp('Therefore,heat added in condenser of mercury cycle=h at 40 bar,0.98 dry-h_feed at 40 bar in KJ/kg')\n",
"h-4.18*150\n",
"disp('Therefore,mercury required per kg of steam=2140.13/heat rejected in condenser in kg per kg of steam')\n",
"2140.13/216.675\n",
"disp('for isentropic expansion,s2=s3=s4=s5=6.9363 KJ/kg K')\n",
"s3=s2;\n",
"s4=s3;\n",
"s5=s4;\n",
"disp('state 3 lies in superheated region,by interpolation the state can be given by,temperature 227.07oc at 8 bar,h3=2899.23 KJ/kg')\n",
"h3=2899.23;\n",
"disp('state 4 lies in wet region,say with dryness fraction x4')\n",
"disp('s4=6.9363=sf at 1 bar+x4*sfg at 1 bar')\n",
"disp('so x4=(s4-sf)/sfg')\n",
"disp('from steam tables,at 1 bar,sf=1.3026 KJ/kg K,sfg=6.0568 KJ/kg K')\n",
"sf=1.3026;\n",
"sfg=6.0568;\n",
"x4=(s4-sf)/sfg\n",
"x4=0.93;//approx.\n",
"disp('h4=hf at 1 bar+x4*hfg at 1 bar in KJ/kg')\n",
"disp('from steam tables,at 1 bar,hf=417.46 KJ/kg,hfg=2258.0 KJ/kg')\n",
"hf=417.46;\n",
"hfg=2258.0;\n",
"h4=hf+x4*hfg\n",
"disp('Let state 5 lie in wet region with dryness fraction x5,')\n",
"disp('s5=6.9363=sf at 0.075 bar+x5*sfg at 0.075 bar')\n",
"disp('so x5=(s5-sf)/sfg')\n",
"disp('from steam tables,at 0.075 bar,sf=0.5764 KJ/kg K,sfg=7.6750 KJ/kg K')\n",
"sf=0.5764;\n",
"sfg=7.6750;\n",
"x5=(s5-sf)/sfg\n",
"x5=0.828;//approx.\n",
"disp('h5=hf+x5*hfg in KJ/kg')\n",
"disp('from steam tables,at 0.075 bar,hf=168.79 KJ/kg,hfg=2406.0 KJ/kg')\n",
"hf=168.79;\n",
"hfg=2406.0;\n",
"h5=hf+x5*hfg\n",
"disp('Let mass of steam bled at 8 bar and 1 bar be m1 and m2 per kg of steam generated.')\n",
"disp('h7=h6+v6*(1-0.075)*10^2 in KJ/kg')\n",
"disp('from steam tables,at 0.075 bar,h6=hf=168.79 KJ/kg,v6=vf=0.001008 m^3/kg')\n",
"h6=168.79;\n",
"v6=0.001008;\n",
"h7=h6+v6*(1-0.075)*10^2\n",
"disp('h9=hf at 1 bar=417.46 KJ/kg,h13=hf at 8 bar=721.11 KJ/kg')\n",
"h9=417.46;\n",
"h13=721.11;\n",
"disp('Applying heat balance on CFWH1,T1=150oc and also T15=150oc')\n",
"T1=150;\n",
"T15=150;\n",
"disp('m1*h3+(1-m1)*h12=m1*h13+(4.18*150)*(1-m1)')\n",
"disp('(m1-2899.23)+(1-m1)*h12=(m1*721.11)+627*(1-m1)')\n",
"disp('Applying heat balance on CFEH2,T11=90oc')\n",
"T11=90;\n",
"disp('m2*h4+(1-m1-m2)*h7=m2*h9+(1-m1-m2)*4.18*90')\n",
"disp('(m2*2517.4)+(1-m1-m2)*168.88=(m2*417.46)+376.2*(1-m1-m2)')\n",
"disp('Heat balance at mixing between CFWH1 and CFWH2,')\n",
"disp('(1-m1-m2)*4.18*90+m2*h10=(1-m1)*h12')\n",
"disp('376.2*(1-m1-m2)+m2*h10=(1-m1)*h12')\n",
"disp('For pumping process,9-10,h10=h9+v9*(8-1)*10^2 in KJ/kg')\n",
"disp('from steam tables,h9=hf at 1 bar=417.46 KJ/kg,v9=vf at 1 bar=0.001043 m^3/kg')\n",
"h9=417.46;\n",
"v9=0.001043;\n",
"h10=h9+v9*(8-1)*10^2 \n",
"disp('solving above equations,we get')\n",
"disp('m1=0.102 kg per kg steam generated')\n",
"disp('m2=0.073 kg per kg steam generated')\n",
"m1=0.102;\n",
"m2=0.073;\n",
"disp('pump work in process 13-14,h14-h13=v13*(40-8)*10^2')\n",
"disp('so h14-h13 in KJ/kg')\n",
"v13=0.001252;\n",
"v13*(40-8)*10^2\n",
"disp('Total heat supplied(q_add)=(9.88*314)+(3330.3-2767.13) in KJ/kg of steam')\n",
"q_add=(9.88*314)+(3330.3-2767.13)\n",
"disp('net work per kg of steam,w_net=w_mercury+w_steam')\n",
"disp('w_net=(9.88*97.325)+{(h2-h3)+(1-m1)*(h3-h4)+(1-m1-m2)*(h4-h5)-(1-m1-m2)*(h4-h6)-m2*(h10-h9)-m1*(h14-h13)} in KJ/kg')\n",
"w_net=(9.88*97.325)+{(h2-h3)+(1-m1)*(h3-h4)+(1-m1-m2)*(h4-h5)-(1-m1-m2)*(h7-h6)-m2*(h10-h9)-m1*4.006}\n",
"disp('thermal efficiency of binary vapour cycle=w_net/q_add')\n",
"w_net/q_add\n",
"disp('in percentage')\n",
"w_net*100/q_add\n",
"disp('so thermal efficiency=55.36%')\n",
"disp('NOTE=>In this question there is some mistake in formula used for w_net which is corrected above.')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.14: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_14.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 14')\n",
"n=0.8;//efficiency of both HP and LP turbine\n",
"P=2500;//output in hp\n",
"disp('This is a mixed pressure turbine so the output of turbine shall be sum of the contributions by HP and LP steam streams.')\n",
"disp('For HP:at inlet of HP steam=>h1=3023.5 KJ/kg,s1=6.7664 KJ/kg K')\n",
"h1=3023.5;\n",
"s1=6.7664;\n",
"disp('ideally, s2=s1=6.7664 KJ/kg K')\n",
"s2=s1;\n",
"disp('s2=sf at 0.075 bar +x3* sfg at 0.075 bar')\n",
"disp('from steam tables,at 0.075 bar,sf=0.5764 KJ/kg K,sfg=7.6750 KJ/kg K')\n",
"sf=0.5764;\n",
"sfg=7.6750;\n",
"disp('so x3=(s2-sf)/sfg')\n",
"x3=(s2-sf)/sfg\n",
"x3=0.806;//approx.\n",
"disp('h_3HP=hf at 0.075 bar+x3*hfg at 0.075 bar in KJ/kg')\n",
"disp('from steam tables,at 0.075 bar,hf=168.79 KJ/kg,hfg=2406.0 KJ/kg')\n",
"hf=168.79;\n",
"hfg=2406.0; \n",
"h_3HP=hf+x3*hfg\n",
"disp('actual enthalpy drop in HP(h_HP)=(h1-h_3HP)*n in KJ/kg')\n",
"h_HP=(h1-h_3HP)*n\n",
"disp('for LP:at inlet of LP steam')\n",
"disp('h2=2706.7 KJ/kg,s2=7.1271 KJ/kg K')\n",
"h2=2706.7;\n",
"s2=7.1271;\n",
"disp('Enthalpy at exit,h_3LP=2222.34 KJ/kg')\n",
"h_3LP=2222.34;\n",
"disp('actual enthalpy drop in LP(h_LP)=(h2-h_3LP)*n in KJ/kg')\n",
"h_LP=(h2-h_3LP)*n\n",
"disp('HP steam consumption at full load=P*0.7457/h_HP in kg/s')\n",
"P*0.7457/h_HP\n",
"disp('HP steam consumption at no load=0.10*(P*0.7457/h_HP)in kg/s')\n",
"0.10*(P*0.7457/h_HP)\n",
"disp('LP steam consumption at full load=P*0.7457/h_LP in kg/s')\n",
"P*0.7457/h_LP\n",
"disp('LP steam consumption at no load=0.10*(P*0.7457/h_LP)in kg/s')\n",
"0.10*(P*0.7457/h_LP)\n",
"disp('The problem can be solved geometrically by drawing willans line as per scale on graph paper and finding out the HP stream requirement for getting 1000 hp if LP stream is available at 1.5 kg/s.')\n",
"disp('or,Analytically the equation for willans line can be obtained for above full load and no load conditions for HP and LP seperately.')\n",
"disp('Willians line for HP:y=m*x+C,here y=steam consumption,kg/s')\n",
"disp('x=load,hp')\n",
"disp('y_HP=m_HP*x+C_HP')\n",
"disp('0.254=m_HP*0+C_HP')\n",
"disp('so C_HP=0.254')\n",
"C_HP=0.254;\n",
"disp('2.54=m_HP*2500+C_HP')\n",
"disp('so m_HP=(2.54-C_HP)/2500')\n",
"m_HP=(2.54-C_HP)/2500\n",
"disp('so y_HP=9.144*10^-4*x_HP+0.254')\n",
"disp('Willans line for LP:y_LP=m_LP*x_LP+C_LP')\n",
"disp('0.481=m_LP*0+C_LP')\n",
"disp('so C_LP=0.481')\n",
"C_LP=0.481;\n",
"disp('4.81=m_LP*2500+C_LP')\n",
"disp('so m_LP=(4.81-C_LP)/2500')\n",
"m_LP=(4.81-C_LP)/2500\n",
"disp('so y_LP=1.732*10^-3*x_LP+0.481')\n",
"disp('Total output(load) from mixed turbine,x=x_HP+x_LP')\n",
"disp('For load of 1000 hp to be met by mixed turbine,let us find out the load shared by LP for steam flow rate of 1.5 kg/s')\n",
"y_LP=1.5;\n",
"disp('from y_LP=1.732*10^-3*x_LP+0.481,')\n",
"disp('x_LP=(y_LP-0.481)/1.732*10^-3')\n",
"x_LP=(y_LP-0.481)/(1.732*10^-3)\n",
"disp('since by 1.5 kg/s of LP steam only 588.34 hp output contribution is made so remaining(1000-588.34)=411.66 hp,411.66 hp should be contributed by HP steam.By willans line for HP turbine,')\n",
"x_HP=411.66;\n",
"disp('from y_HP=9.144*10^-4*x_HP+C_HP in kg/s')\n",
"y_HP=9.144*10^-4*x_HP+C_HP\n",
"disp('so HP steam requirement=0.63 kg/s')\n",
"\n",
""
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.15: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_15.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 15')\n",
"disp('Let us carry out analysis for 1 kg of steam generated in boiler.')\n",
"disp('Enthalpy at inlet to HPT,h2=2960.7 KJ/kg,s2=6.3615 KJ/kg K')\n",
"h2=2960.7;\n",
"s2=6.3615;\n",
"disp('state at 3  i.e. exit from HPT can be identified by s2=s3=6.3615 KJ/kg K')\n",
"s3=s2;\n",
"disp('Let dryness fraction be x3,s3=6.3615=sf at 2 bar+x3*sfg at 2 bar')\n",
"disp('so x3=(s3-sf)/sfg')\n",
"disp('from stem tables,at 2 bar,sf=1.5301 KJ/kg K,sfg=5.5970 KJ/kg K')\n",
"sf=1.5301;\n",
"sfg=5.5970;\n",
"x3=(s3-sf)/sfg\n",
"x3=0.863;//approx.\n",
"disp('h3=2404.94 KJ/kg')\n",
"h3=2404.94;\n",
"disp('If one kg of steam is generated in bolier then at exit of HPT,0.5 kg goes into process heater and 0.5 kg goes into separator')\n",
"disp('mass of moisture retained in separator(m)=(1-x3)*0.5 kg')\n",
"m=(1-x3)*0.5\n",
"disp('Therefore,mass of steam entering LPT(m_LP)=0.5-m kg')\n",
"m_LP=0.5-m\n",
"disp('Total mass of water entering hot well at 8(i.e. from process heater and drain from separator)=(0.5+0.685)=0.5685 kg')\n",
"disp('Let us assume the temperature of water leaving hotwell be T oc.Applying heat balance for mixing;')\n",
"disp('(0.5685*4.18*90)+(0.4315*4.18*40)=(1*4.18*T)')\n",
"disp('so T=((0.5685*4.18*90)+(0.4315*4.18*40))/4.18 in degree celcius')\n",
"T=((0.5685*4.18*90)+(0.4315*4.18*40))/4.18\n",
"disp('so temperature of water leaving hotwell=68.425 degree celcius')\n",
"disp('Applying heat balanced on trap')\n",
"disp('0.5*h7+0.0685*hf at 2 bar=(0.5685*4.18*90)')\n",
"disp('so h7=((0.5685*4.18*90)-(0.0685*hf))/0.5 in KJ/kg')\n",
"disp('from steam tables,at 2 bar,hf=504.70 KJ/kg')\n",
"hf=504.70;\n",
"h7=((0.5685*4.18*90)-(0.0685*hf))/0.5\n",
"disp('Therefore,heat transferred in process heater=0.5*(h3-h7)in KJ/kg steam generated')\n",
"0.5*(h3-h7)\n",
"disp('so heat transferred per kg steam generated=1023.175 KJ/kg steam generated')\n",
"disp('For state 10 at exit of LPT,s10=s3=s2=6.3615 KJ/kg K')\n",
"s10=s3;\n",
"disp('Let dryness fraction be x10')\n",
"disp('s10=6.3615=sf at 0.075 bar+x10*sfg at 0.075 bar')\n",
"disp('from steam tables,at 0.075 bar,sf=0.5764 KJ/kg K,sfg=7.6750 KJ/kg K')\n",
"sf=0.5764;\n",
"sfg=7.6750;\n",
"disp('so x10=(s10-sf)/sfg')\n",
"x10=(s10-sf)/sfg\n",
"x10=0.754;//approx.\n",
"disp('h10=hf at 0.075 bar+x10*hfg at 0.075 bar')\n",
"disp('from steam tables,at 0.075 bar,hf=168.79 KJ/kg,hfg=2406.0 KJ/kg')\n",
"hf=168.79;\n",
"hfg=2406.0;\n",
"disp('so h10=hf+x10*hfg in KJ/kg ')\n",
"h10=hf+x10*hfg \n",
"disp('net work output,neglecting pump work per kg of steam generated,')\n",
"disp('w_net=(h2-h3)*1+0.4315*(h3-h10) in KJ/kg steam generated')\n",
"w_net=(h2-h3)*1+0.4315*(h3-h10) \n",
"disp('Heat added in boiler per kg steam generated,q_add=(h2-h1)in KJ/kg')\n",
"q_add=(h2-4.18*68.425)\n",
"disp('thermal efficiency=w_net/q_add')\n",
"w_net/q_add\n",
"disp('in percentage')\n",
"w_net*100/q_add\n",
"disp('so Thermal efficiency=27.58%')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.16: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_16.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 16')\n",
"m=35;//mass flow rate in kg/s\n",
"disp('from steam tables,h1=3530.9 KJ/kg,s1=6.9486 KJ/kg K')\n",
"h1=3530.9;\n",
"s1=6.9486;\n",
"disp('Assuming isentropic expansion in nozzle,s1=s2=6.9486')\n",
"s2=s1;\n",
"disp('Letdryness fraction at state 2,x2=0.864')\n",
"disp('s2=sf at 0.2 bar+x2*sfg at 0.2 bar')\n",
"disp('from steam tables,sf=0.8320 KJ/kg K,sfg=7.0766 KJ/kg K')\n",
"sf=0.8320;\n",
"sfg=7.0766;\n",
"disp('so x2=(s2-sf)/sfg')\n",
"x2=(s2-sf)/sfg\n",
"x2=0.864;//approx.\n",
"disp('hence,h2=hf at 0.2 bar+x2*hfg at 0.2 bar in KJ/kg')\n",
"disp('from steam tables,hf at 0.2 bar=251.4 KJ/kg,hfg at 0.2 bar=2358.3 KJ/kg')\n",
"hf=251.4;\n",
"hfg=2358.3;\n",
"h2=hf+x2*hfg\n",
"disp('considering pump work to be of isentropic type,deltah_34=v3*deltap_34')\n",
"disp('from steam table,v3=vf at 0.2 bar=0.001017 m^3/kg')\n",
"v3=0.001017;\n",
"disp('or deltah_34=v3*(p3-p4)in KJ/kg')\n",
"p3=70;//;pressure of steam entering turbine in bar\n",
"p4=0.20;//condenser pressure in bar\n",
"deltah_34=v3*(p3-p4)*100\n",
"disp('pump work,Wp=deltah_34 in KJ/kg')\n",
"Wp=deltah_34\n",
"disp('turbine work,Wt=deltah_12=(h1-h2)in KJ/kg')\n",
"Wt=(h1-h2)\n",
"disp('net work(W_net)=Wt-Wp in KJ/kg')\n",
"W_net=Wt-Wp\n",
"disp('power produced(P)=mass flow rate*W_net in KJ/s')\n",
"P=m*W_net\n",
"disp('so net power=43.22 MW')\n",
"disp('heat supplied in boiler(Q)=(h1-h4) in KJ/kg')\n",
"disp('enthalpy at state 4,h4=h3+deltah_34 in KJ/kg')\n",
"h3=hf;\n",
"h4=h3+deltah_34 \n",
"disp('total heat supplied to boiler(Q)=m*(h1-h4)in KJ/s')\n",
"Q=m*(h1-h4)\n",
"disp('thermal efficiency=net work/heat supplied=W_net/Q')\n",
"P/Q\n",
"disp('in percentage')\n",
"P*100/Q\n",
"disp('so thermal efficiency=37.73%')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.17: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_17.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 17')\n",
"P=10*10^3;//output in KW\n",
"disp('from steam tables,h1=3625.3 KJ/s,s1=6.9029 KJ/kg K')\n",
"h1=3625.3;\n",
"s1=6.9029;\n",
"disp('due to isentropic expansion,s1=s2=s3=6.9029 KJ/kg K')\n",
"s2=s1;\n",
"s3=s2;\n",
"disp('at state 2,i.e at pressure of 2 MPa and entropy 6.9029 KJ/kg K')\n",
"disp('by interpolating state for s2 between 2 MPa,300 degree celcius and 2 MPa,350 degree celcius from steam tables,')\n",
"disp('h2=3105.08 KJ/kg ')\n",
"h2=3105.08;\n",
"disp('for state 3,i.e at pressure of 0.01 MPa entropy,s3 lies in wet region as s3<sg at 0.01 MPa.Let dryness fraction be x3 at this state')\n",
"disp('s3=sf at 0.01 MPa+x3*sfg at 0.01 MPa')\n",
"disp('from steam tables,sf at 0.01 MPa=0.6493 KJ/kg K,sfg at 0.01 MPa=7.5009 KJ/kg K')\n",
"sf=0.6493;\n",
"sfg=7.5009;\n",
"disp('so x3=(s3-sf)/sfg')\n",
"x3=(s3-sf)/sfg\n",
"x3=0.834;//approx.\n",
"disp('enthalpy at state 3,h3= hf at 0.01 MPa+x3*hfg at 0.01 MPa in KJ/kg')\n",
"disp('from steam tables,at 0.01 MPa,hf=191.83 KJ/kg,hfg=2392.8 KJ/kg')\n",
"hf=191.83;\n",
"hfg=2392.8;\n",
"h3=hf+x3*hfg\n",
"disp('let the mass of steam bled be mb per kg of steam from exit of HP for regenerative feed heating.')\n",
"disp('Considering state at exit from feed heater being saturated liquid the enthalpy at exit of feed heater will be,hf at 2 MPa')\n",
"disp('h6=hf at 2 MPa=908.79 KJ/kg')\n",
"h6=908.79;\n",
"disp('for adiabatic mixing in feed heater,for one kg of steam leaving boiler,the heat balance yields,')\n",
"disp('(1-mb)*h5+mb*h2=h6')\n",
"disp('while neglecting pump work,')\n",
"disp('h5=h4=hf at 0.01MPa=191.83 KJ/kg')\n",
"h4=191.83;\n",
"h5=h4;\n",
"disp('substituting in heat balance on the feed heater,')\n",
"disp('(1-mb)*h5+mb*h2=h6')\n",
"disp('so mb=(h6-h5)/(h2-h5)in kg per kg of steam entering HP turbine')\n",
"mb=(h6-h5)/(h2-h5)\n",
"mb=0.246;//approx.\n",
"disp('steam bled per kg of steam passing through HP stage=0.246 kg')\n",
"disp('let mass of steam leaving boiler be m kg/s')\n",
"disp('output(P)=m*(h1-h2)+m*(1-mb)*(h2-h3)')\n",
"disp('so m=P/((h1-h2)+(1-mb)*(h2-h3))in kg/s')\n",
"m=P/((h1-h2)+(1-mb)*(h2-h3))\n",
"m=8.25;//approx.\n",
"disp('neglecting pump work,h7=h6=908.79 KJ/kg')\n",
"h7=h6;\n",
"disp('heat supplied to boiler,Q_71=m*(h1-h7) in KJ/s')\n",
"Q_71=m*(h1-h7)\n",
"disp('so heat added=22411.21 KJ/s')\n",
"\n",
"\n",
""
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.18: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_18.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 18')\n",
"W_net=50*10^3;//net output of turbine in KW\n",
"disp('from steam tables,at inlet to first stage of turbine,h1=h at 100 bar,500oc=3373.7 KJ/kg,s1=s at 100 bar,500oc=6.5966 KJ/kg')\n",
"h1=3373.7;\n",
"s1=6.5966;\n",
"disp('Due to isentropic expansion,s1=s6=s2 and s3=s8=s4')\n",
"s2=s1;\n",
"s6=s2;\n",
"disp('State at 6 i.e bleed state from HP turbine,temperature by interpolation from steam table =261.6oc.')\n",
"disp('At inlet to second stage of turbine,h6=2930.572 KJ/kg')\n",
"h6=2930.572;\n",
"disp('h3=h at 10 bar,500oc=3478.5 KJ/kg,s3=s at 10 bar,500oc=7.7622 KJ/kg K')\n",
"h3=3478.5;\n",
"s3=7.7622;\n",
"s4=s3;\n",
"s8=s4;\n",
"disp('At exit from first stage of turbine i.e. at 10 bar and entropy of 6.5966 KJ/kg K,Temperature by interpolation from steam table at 10 bar and entropy of 6.5966 KJ/kg K')\n",
"disp('T2=181.8oc,h2=2782.8 KJ/kg')\n",
"T2=181.8;\n",
"h2=2782.8;\n",
"disp('state at 8,i.e bleed state from second stage of expansion,i.e at 4 bar and entropy of 7.7622 KJ/kg K,Temperature by interpolation from steam table,T8=358.98oc=359oc')\n",
"T8=359;\n",
"disp('h8=3188.7 KJ/kg')\n",
"h8=3188.7;\n",
"disp('state at 4 i.e. at condenser pressure of 0.1 bar and entropy of 7.7622 KJ/kg K,the state lies in wet region.So let the dryness fraction be x4.')\n",
"disp('s4=sf at 0.1 bar+x4*sfg at 0.1 bar')\n",
"disp('from steam tables,at 0.1 bar,sf=0.6493 KJ/kg K,sfg=7.5009 KJ/kg K')\n",
"sf=0.6493;\n",
"sfg=7.5009; \n",
"disp('so x4=(s4-sf)/sfg')\n",
"x4=(s4-sf)/sfg\n",
"x4=0.95;//approx.\n",
"disp('h4=hf at 0.1 bar+x4*hfg at 0.1 bar in KJ/kg ')\n",
"disp('from steam tables,at 0.1 bar,hf=191.83 KJ/kg,hfg=2392.8 KJ/kg')\n",
"hf=191.83;\n",
"hfg=2392.8;\n",
"h4=hf+x4*hfg\n",
"disp('given,h4=2464.99 KJ/kg,h11=856.8 KJ/kg,h9=hf at 4 bar=604.74 KJ/kg')\n",
"h4=2464.99;\n",
"h11=856.8;\n",
"h9=604.74;\n",
"disp('considering pump work,the net output can be given as,')\n",
"disp('W_net=W_HPT+W_LPT-(W_CEP+W_FP)')\n",
"disp('where,W_HPT={(h1-h6)+(1-m6)*(h6-h2)}per kg of steam from boiler.')\n",
"disp('W_LPT={(1-m6)+(h3-h8)*(1-m6-m8)*(h8-h4)}per kg of steam from boiler.')\n",
"disp('for closed feed water heater,energy balance yields;')\n",
"disp('m6*h6+h10=m6*h7+h11')\n",
"disp('assuming condensate leaving closed feed water heater to be saturated liquid,')\n",
"disp('h7=hf at 20 bar=908.79 KJ/kg')\n",
"h7=908.79; \n",
"disp('due to throttline,h7=h7_a=908.79 KJ/kg')\n",
"h7_a=h7;\n",
"disp('for open feed water heater,energy balance yields,')\n",
"disp('m6*h7_a+m8*h8+(1-m6-m8)*h5=h9')\n",
"disp('for condensate extraction pump,h5-h4_a=v4_a*deltap')\n",
"disp('h5-hf at 0.1 bar=vf at 0.1 bar*(4-0.1)*10^2 ')\n",
"disp('from steam tables,at 0.1 bar,hf=191.83 KJ/kg,vf=0.001010 m^3/kg')\n",
"hf=191.83;\n",
"vf=0.001010; \n",
"disp('so h5=hf+vf*(4-0.1)*10^2 in KJ/kg')\n",
"h5=hf+vf*(4-0.1)*10^2\n",
"disp('for feed pump,h10-h9=v9*deltap')\n",
"disp('h10=h9+vf at 4 bar*(100-4)*10^2 in KJ/kg')\n",
"disp('from steam tables,at 4 bar,hf=604.74 KJ/kg,vf=0.001084 m^3/kg ')\n",
"hf=604.74;\n",
"vf=0.001084;\n",
"h10=h9+vf*(100-4)*10^2\n",
"disp('substituting in energy balance upon closed feed water heater,')\n",
"disp('m6=(h11-h10)/(h6-h7)in kg per kg of steam from boiler')\n",
"m6=(h11-h10)/(h6-h7)\n",
"disp('substituting in energy balance upon feed water heater,')\n",
"disp('m8=(h9-m6*h7_a+m6*h5-h5)/(h8-h5)in kg per kg of steam from boiler')\n",
"m8=(h9-m6*h7_a+m6*h5-h5)/(h8-h5)\n",
"disp('Let the mass of steam entering first stage of turbine be m kg,then')\n",
"disp('W_HPT=m*{(h1-h6)+(1-m6)*(h6-h2)}')\n",
"disp('W_HPT/m=')\n",
"{(h1-h6)+(1-m6)*(h6-h2)}\n",
"disp('so W_HPT=m*573.24 KJ')\n",
"disp('also,W_LPT={(1-m6)*(h3-h8)+(1-m6-m8)*(h8-h4)}per kg of steam from boiler')\n",
"disp('W_LPT/m=')\n",
"{(1-m6)*(h3-h8)+(1-m6-m8)*(h8-h4)}\n",
"disp('so W_LPT=m*813.42 KJ')\n",
"disp('pump works(negative work)')\n",
"disp('W_CEP=m*(1-m6-m8)*(h5-h4_a)')\n",
"h4_a=191.83;//h4_a=hf at 0.1 bar\n",
"disp('W_CEP/m=')\n",
"(1-m6-m8)*(h5-h4_a)\n",
"disp('so W_CEP=m* 0.304')\n",
"disp('W_FP=m*(h10-h9)')\n",
"disp('W_FP/m=')\n",
"(h10-h9)\n",
"disp('so W_FP=m*10.41')\n",
"disp('net output,')\n",
"disp('W_net=W_HPT+W_LPT-W_CEP-W_FP ')\n",
"disp('so 50*10^3=(573.24*m+813.42*m-0.304*m-10.41*m)')\n",
"disp('so m=W_net/(573.24+813.42-0.304-10.41)in kg/s')\n",
"m=W_net/(573.24+813.42-0.304-10.41)\n",
"disp('heat supplied in boiler,Q_add=m*(h1-h11) in KJ/s')\n",
"Q_add=m*(h1-h11)\n",
"disp('Thermal efficenncy=W_net/Q_add')\n",
"W_net/Q_add\n",
"disp('in percentage')\n",
"W_net*100/Q_add\n",
"disp('so mass of steam bled at 20 bar=0.119 kg per kg of steam entering first stage')\n",
"disp('mass of steam bled at 4 bar=0.109 kg per kg of steam entering first stage')\n",
"disp('mass of steam entering first stage=36.33 kg/s')\n",
"disp('thermal efficiency=54.66%')\n",
"disp('NOTE=>In this question there is some caclulation mistake while calculating m6 in book,which is corrected above so some answers may vary.')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.1: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_1.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 1')\n",
"disp('T-S representation for carnot cycle operating between pressure of 7 MPa and 7KPa is shown in fig.')\n",
"disp('enthalpy at state 2,h2= hg at 7 MPa')\n",
"disp('from steam table,h=2772.1 KJ/kg')\n",
"h2=2772.1;\n",
"disp('entropy at state 2,s2=sg at 7MPa')\n",
"disp('from steam table,s2=5.8133 KJ/kg K')\n",
"s2=5.8133;\n",
"disp('enthalpy and entropy at state 3,')\n",
"disp('from steam table,h3=hf at 7 MPa =1267 KJ/kg and s3=sf at 7 MPa=3.1211 KJ/kg K')\n",
"h3=1267;\n",
"s3=3.1211;\n",
"disp('for process 2-1,s1=s2.Let dryness fraction at state 1 be x1 ')\n",
"s1=s2;\n",
"disp('from steam table, sf at 7 KPa=0.5564 KJ/kg K,sfg at 7 KPa=7.7237 KJ/kg K')\n",
"sf=0.5564;\n",
"sfg=7.7237;\n",
"disp('s1=s2=sf+x1*sfg')\n",
"disp('so x1=(s2-sf)/sfg ')\n",
"x1=(s2-sf)/sfg \n",
"x1=0.6806;//approx.\n",
"disp('from steam table,hf at 7 KPa=162.60 KJ/kg,hfg at 7 KPa=2409.54 KJ/kg')\n",
"hf=162.60;\n",
"hfg=2409.54;\n",
"disp('enthalpy at state 1,h1=hf+x1*hfg in KJ/kg')\n",
"h1=hf+x1*hfg\n",
"disp('let dryness fraction at state 4 be x4')\n",
"disp('for process 4-3,s4=s3=sf+x4*sfg')\n",
"s4=s3;\n",
"disp('so x4=(s4-sf)/sfg')\n",
"x4=(s4-sf)/sfg\n",
"x4=0.3321;//approx.\n",
"disp('enthalpy at state 4,h4=hf+x4*hfg in KJ/kg')\n",
"h4=hf+x4*hfg\n",
"disp('thermal efficiency=net work/heat added')\n",
"disp('expansion work per kg=(h2-h1) in KJ/kg')\n",
"(h2-h1)\n",
"disp('compression work per kg=(h3-h4) in KJ/kg(+ve)')\n",
"(h3-h4)\n",
"disp('heat added per kg=(h2-h3) in KJ/kg(-ve)')\n",
"(h2-h3)\n",
"disp('net work per kg=(h2-h1)-(h3-h4) in KJ/kg')\n",
"(h2-h1)-(h3-h4)\n",
"disp('thermal efficiency')\n",
"((h2-h1)-(h3-h4))/(h2-h3)\n",
"disp('in percentage')\n",
"(((h2-h1)-(h3-h4))/(h2-h3))*100\n",
"disp('so thermal efficiency=44.21%')\n",
"disp('turbine work=969.57 KJ/kg(+ve)')\n",
"disp('compression work=304.19 KJ/kg(-ve)')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.2: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_2.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 2')\n",
"disp('from steam tables,at 5 MPa,hf_5MPa=1154.23 KJ/kg,sf_5MPa=2.92 KJ/kg K')\n",
"disp('hg_5MPa=2794.3 KJ/kg,sg_5MPa=5.97 KJ/kg K')\n",
"hf_5MPa=1154.23;\n",
"sf_5MPa=2.92;\n",
"hg_5MPa=2794.3;\n",
"sg_5MPa=5.97;\n",
"disp('from steam tables,at 5 Kpa,hf_5KPa=137.82 KJ/kg,sf_5KPa=0.4764 KJ/kg K')\n",
"disp('hg_5KPa=2561.5 KJ/kg,sg_5KPa=8.3951 KJ/kg K,vf_5KPa=0.001005 m^3/kg')\n",
"hf_5KPa=137.82;\n",
"sf_5KPa=0.4764;\n",
"hg_5KPa=2561.5;\n",
"sg_5KPa=8.3951;\n",
"vf_5KPa=0.001005;\n",
"disp('as process 2-3 is isentropic,so s2=s3')\n",
"disp('and s3=sf_5KPa+x3*sfg_5KPa=s2=sg_5MPa')\n",
"s2=sg_5MPa;\n",
"s3=s2;\n",
"disp('so x3=(s3-sf_5KPa)/sfg_5KPa')\n",
"x3=(s3-sf_5KPa)/(sg_5KPa-sf_5KPa)\n",
"x3=0.694;//approx.\n",
"disp('hence enthalpy at 3,')\n",
"disp('h3=hf_5KPa+x3*hfg_5KPa in KJ/kg')\n",
"h3=hf_5KPa+x3*(hg_5KPa-hf_5KPa)\n",
"disp('enthalpy at 2,h2=hg_5KPa=2794.3 KJ/kg')\n",
"disp('process 1-4 is isentropic,so s1=s4')\n",
"s1=sf_5MPa;\n",
"disp('s1=sf_5KPa+x4*(sg_5KPa-sf_5KPa)')\n",
"disp('so x4=(s1-sf_5KPa)/(sg_5KPa-sf_5KPa)')\n",
"x4=(s1-sf_5KPa)/(sg_5KPa-sf_5KPa)\n",
"x4=0.308;//approx.\n",
"disp('enthalpy at 4,h4=hf_5KPa+x4*(hg_5KPa-hf_5KPa)in KJ/kg')\n",
"h4=hf_5KPa+x4*(hg_5KPa-hf_5KPa)\n",
"disp('enthalpy at 1,h1=hf_5MPa in KJ/kg')\n",
"h1=hf_5MPa\n",
"disp('carnot cycle(1-2-3-4-1) efficiency:')\n",
"disp('n_carnot=net work/heat added')\n",
"disp('n_carnot=((h2-h3)-(h1-h4))/(h2-h1)')\n",
"h2=hg_5MPa;\n",
"n_carnot=((h2-h3)-(h1-h4))/(h2-h1)\n",
"disp('in percentage')\n",
"n_carnot=n_carnot*100\n",
"disp('so n_carnot=42.95%')\n",
"disp('In rankine cycle,1-2-3-5-6-1,')\n",
"disp('pump work,h6-h5=vf_5KPa*(p6-p5)in KJ/kg')\n",
"p6=5000;//boiler pressure in KPa\n",
"p5=5;//condenser pressure in KPa\n",
"vf_5KPa*(p6-p5)\n",
"disp('h5=hf_5KPa=137.82 KJ/kg')\n",
"h5=hf_5KPa;\n",
"disp('hence h6 in KJ/kg')\n",
"h6=h5+(vf_5KPa*(p6-p5))\n",
"disp('net work in rankine cycle=(h2-h3)-(h6-h5)in KJ/kg')\n",
"(h2-h3)-(h6-h5)\n",
"disp('heat added=(h2-h6)in KJ/kg')\n",
"(h2-h6)\n",
"disp('rankine cycle efficiency(n_rankine)=')\n",
"n_rankine=((h2-h3)-(h6-h5))/(h2-h6)\n",
"disp('in percentage')\n",
"n_rankine=n_rankine*100\n",
"disp('so n_rankine=36.56%')\n",
""
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.3: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_3.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 3')\n",
"disp('from steam tables,h2=hg_40bar=3092.5 KJ/kg')\n",
"h2=3092.5;\n",
"disp('s2=sg_40bar=6.5821 KJ/kg K')\n",
"s2=6.5821;\n",
"disp('h4=hf_0.05bar=137.82 KJ/kg,hfg=2423.7 KJ/kg ')\n",
"h4=137.82;\n",
"hfg=2423.7;\n",
"disp('s4=sf_0.05bar=0.4764 KJ/kg K,sfg=7.9187 KJ/kg K')\n",
"s4=0.4764;\n",
"sfg=7.9187;\n",
"disp('v4=vf_0.05bar=0.001005 m^3/kg')\n",
"v4=0.001005;\n",
"disp('let the dryness fraction at state 3 be x3,')\n",
"disp('for ideal process,2-3,s2=s3')\n",
"s3=s2;\n",
"disp('s2=s3=6.5821=sf_0.05bar+x3*sfg_0.05bar')\n",
"disp('so x3=(s2-s4)/(sfg)')\n",
"x3=(s2-s4)/(sfg)\n",
"x3=0.7711;//approx.\n",
"disp('h3=hf_0.05bar+x3*hfg_0.05bar in KJ/kg')\n",
"h3=h4+x3*hfg\n",
"disp('for pumping process,')\n",
"disp('h1-h4=v4*deltap=v4*(p1-p4)')\n",
"disp('so h1=h4+v4*(p1-p4) in KJ/kg')\n",
"p1=40*100;//pressure of steam enter in turbine in mPa\n",
"p4=0.05*100;//pressure of steam leave turbine in mPa\n",
"h1=h4+v4*(p1-p4)\n",
"disp('pump work per kg of steam=(h1-h4) in KJ/kg')\n",
"(h1-h4)\n",
"disp('net work per kg of steam =(expansion work-pump work)per kg of steam')\n",
"disp('=(h2-h3)-(h1-h4) in KJ/kg)')\n",
"(h2-h3)-(h1-h4)\n",
"disp('cycle efficiency=net work/heat added')\n",
"((h2-h3)-(h1-h4))/(h2-h1)\n",
"disp('in percentage')\n",
"((h2-h3)-(h1-h4))*100/(h2-h1)\n",
"disp('so net work per kg of steam=1081.74 KJ/kg')\n",
"disp('cycle efficiency=36.67%')\n",
"disp('pump work per kg of steam=4.02 KJ/kg')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.4: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_4.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 4')\n",
"disp('Let us assume that the condensate leaves condenser as saturated liquid and the expansion in turbine and pumping processes are isentropic.')\n",
"disp('from steam tables,h2=h_20MPa=3238.2 KJ/kg')\n",
"h2=3238.2;\n",
"disp('s2=6.1401 KJ/kg K')\n",
"s2=6.1401;\n",
"disp('h5=h_0.005MPa in KJ/kg')\n",
"disp('from steam tables,at 0.005 MPa,hf=137.82 KJ/kg,hfg=2423.7 KJ/kg,sf=0.4764 KJ/kg K,sfg=7.9187 KJ/kg K')\n",
"hf=137.82;\n",
"hfg=2423.7;\n",
"sf=0.4764;\n",
"sfg=7.9187;\n",
"disp('h5=hf+0.9*hfg in KJ/kg')\n",
"h5=hf+0.9*hfg\n",
"disp('s5=sf+0.9*sfg in KJ/kg K')\n",
"s5=sf+0.9*sfg\n",
"disp('h6=hf=137.82 KJ/kg')\n",
"h6=137.82;\n",
"disp('it is given that temperature at state 4 is 500 degree celcius and due to isentropic processes s4=s5=7.6032 KJ/kg K.The state 4 can be conveniently located on mollier chart by the intersection of 500 degree celcius constant temperature line and entropy value of 7.6032 KJ/kg K and the pressure and enthalpy obtained.but these shall be approximate.')\n",
"disp('The state 4 can also be located by interpolation using steam table.The entropy value of 7.6032 KJ/kg K lies between the superheated steam states given under,p=1.20 MPa,s at 1.20 MPa=7.6027 KJ/kg K')\n",
"disp('p=1.40 MPa,s at 1.40 MPa=7.6027 KJ/kg K')\n",
"disp('by interpolation state 4 lies at pressure=')\n",
"1.20+((1.40-1.20)/(7.6027-7.6759))*(7.6032-7.6759)\n",
"disp('=1.399,approx.=1.40 MPa')\n",
"disp('thus,steam leaves HP turbine at 1.40 MPa')\n",
"disp('enthalpy at state 4,h4=3474.1 KJ/kg')\n",
"h4=3474.1;\n",
"disp('for process 2-33,s2=s3=6.1401 KJ/kg K.The state 3 thus lies in wet region as s3<sg at 1.40 MPa.Let dryness fraction at state 3 be x3.')\n",
"s3=s2;\n",
"disp('s3=sf+x3*sfg')\n",
"disp('from staem tables,at 1.4 MPa,sf=2.2842 KJ/kg K,sfg=4.1850 KJ/kg K')\n",
"sf=2.2842;\n",
"sfg=4.1850;\n",
"disp('so x3=(s3-sf)/sfg')\n",
"x3=(s3-sf)/sfg\n",
"disp('h3=hf+x3*hfg in KJ/kg')\n",
"disp('from steam tables,at 1.4 MPa,hf=830.3 KJ/kg,hfg=1959.7 KJ/kg')\n",
"hf=830.3;\n",
"hfg=1959.7;\n",
"h3=hf+x3*hfg\n",
"disp('enthalpy at 1,h1=h6+v6*(p1-p6) in KJ/kg')\n",
"disp('h1=hf at 0.005MPa+vf at 0.005MPa*(p1-p6)')\n",
"disp('from steam tables, at 0.005 MPa,h6=137.82 KJ/kg,v6=0.001005 m^3/kg')\n",
"h6=137.82;\n",
"v6=0.001005;\n",
"p1=20*1000;//steam entering HP turbine in KPa\n",
"p6=0.005*1000;//condensor pressure in KPa\n",
"h1=h6+v6*(p1-p6)\n",
"disp('net work per kg steam=(h2-h3)+(h4-h5)-(h1-h6)in KJ/kg')\n",
"(h2-h3)+(h4-h5)-(h1-h6)\n",
"disp('heat added per kg of steam=(h2-h1)in KJ/kg')\n",
"(h2-h1)\n",
"disp('thermal efficiency=net work/heat added')\n",
"((h2-h3)+(h4-h5)-(h1-h6))/(h2-h1)\n",
"disp('in percentage')\n",
"(((h2-h3)+(h4-h5)-(h1-h6))/(h2-h1))*100\n",
"disp('pressure of steam leaving HP turbine=1.40 MPa')\n",
"disp('thermal efficiency=56.39%')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.5: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_5.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 5')\n",
"P=50*10^3;//output of plant in KW\n",
"Cpw=4.18;//specific heat of water in KJ/kg K\n",
"Tw_in=15;//cooling water entering condenser temperature in degree celcius\n",
"Tw_out=30;//cooling water leaving condenser temperature in degree celcius\n",
"disp('from steam table,at inlet to turbine,')\n",
"disp('h2=h_10MPa,700oc')\n",
"disp('h2=3870.5 KJ/kg,s2=7.1687 KJ/kg K')\n",
"h2=3870.5;\n",
"s2=7.1687;\n",
"s3=s2;\n",
"disp('for process 2-3,s2=s3 and s3<sf at 0.005 MPa so state 3 lies in wet region .Let dryness fraction at state 3 be x3.')\n",
"disp('s3=7.1687=sf at 0.005 MPa+x3*sfg at 0.005 MPa')\n",
"disp('from steam tables,at 0.005 MPa,sf=0.4764 KJ/kg K,sfg=7.9187 KJ/kg')\n",
"sf=0.4764;\n",
"sfg=7.9187;\n",
"disp('so x3=(s3-sf)/sfg')\n",
"x3=(s3-sf)/sfg\n",
"x3=0.845;//approx.\n",
"disp('h3=hf at 0.005 MPa+x3*hfg at 0.005 MPa')\n",
"disp('from steam tables,at 0.005 MPa,hf=137.82 KJ/kg,hfg=2423.7 KJ/kg')\n",
"hf=137.82;\n",
"hfg=2423.7;\n",
"disp('so h3=hf+x3*hfg in KJ/kg')\n",
"h3=hf+x3*hfg\n",
"disp('h4=hf at 0.005 MPa')\n",
"h4=hf;\n",
"disp('for pumping process,(h1-h4)=v4*(p1-p4)')\n",
"disp('from steam tables,v4=vf at 0.005 MPa=0.001005 m^3/kg')\n",
"v4=0.001005;\n",
"disp('h1=h4+v4*(p1-p4)in KJ/kg')\n",
"p1=10;//pressure of steam leave boiler in MPa\n",
"p4=0.005;//pressure of steam leave turbine in MPa\n",
"h1=h4+v4*(p1-p4)*100\n",
"disp('net output per kg of steam,w_net=(h2-h3)-(h1-h4) in KJ/kg')\n",
"w_net=(h2-h3)-(h1-h4)\n",
"disp('mass flow rate of steam,ms=P/w_net in kg/s')\n",
"ms=P/w_net\n",
"ms=29.69;//approx.\n",
"disp('by heat balance on condenser,for mass flow rate of water being mw kg/s')\n",
"disp('(h3-h4)*ms=mw*Cpw*(Tw_out-Tw_in)')\n",
"disp('so mw=(h3-h4)*ms/(Cpw*(Tw_out-Tw_in)) in kg/s')\n",
"mw=(h3-h4)*ms/(Cpw*(Tw_out-Tw_in))\n",
"disp('the heat added per kg of steam (q_add)=(h2-h1)in KJ/kg')\n",
"q_add=(h2-h1)\n",
"disp('thermal efficiency=w_net/q_add')\n",
"w_net/q_add\n",
"disp('in percentage')\n",
"w_net*100/q_add\n",
"disp('ratio of heat supplied and rejected=(h2-h1)/(h3-h4)')\n",
"(h2-h1)/(h3-h4)\n",
"disp('mass of flow rate of steam=29.69 kg/s')\n",
"disp('mass flow rate of condenser cooling water=969.79 kg/s')\n",
"disp('thermal efficiency=45.12%')\n",
"disp('ratio of heat supplied and rejected=1.822')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.6: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_6.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 6')\n",
"disp('case (a) When there is no feed water heater')\n",
"disp('Thermal efficiency of cycle=((h2-h3)-(h1-h4))/(h2-h1)')\n",
"disp('from steam tables,h2=h at 200 bar,650oc=3675.3 KJ/kg,s2=s at 200 bar,650oc=6.6582 KJ/kg K,h4=hf at 0.05 bar=137.82 KJ/kg,v4=vf at 0.05 bar=0.001005 m^3/kg')\n",
"h2=3675.3;\n",
"s2=6.6582;\n",
"h4=137.82;\n",
"v4=0.001005;\n",
"disp('hf at 0.05 bar=137.82 KJ/kg,hfg at 0.05 bar=2423.7 KJ/kg,sf at 0.05 bar=0.4764 KJ/kg K,sfg at 0.05 bar=7.9187 KJ/kg K')\n",
"hf=137.82;\n",
"hfg=2423.7;\n",
"sf=0.4764;\n",
"sfg=7.9187;\n",
"disp('For process 2-3,s2=s3.Let dryness fraction at 3 be x3.')\n",
"s3=s2;\n",
"disp('s3=6.6582=sf at 0.05 bar+x3*sfg at 0.05 bar')\n",
"disp('so x3=(s3-sf)/sfg')\n",
"x3=(s3-sf)/sfg\n",
"x3=0.781;//approx.\n",
"disp('h3=hf at 0.05 bar+x3*hfg at 0.05 bar in KJ/kg')\n",
"h3=hf+x3*hfg\n",
"disp('For pumping process 4-1,')\n",
"disp('h1-h4=v4*deltap')\n",
"disp('h1=h4+v4*(200-0.5)*10^2 in KJ/kg')\n",
"h1=h4+v4*(200-0.5)*10^2\n",
"disp('Thermal efficiency of cycle=')\n",
"((h2-h3)-(h1-h4))/(h2-h1)\n",
"disp('in percentage')\n",
"((h2-h3)-(h1-h4))*100/(h2-h1)\n",
"disp('case (b) When there is only one feed water heater working at 8 bar')\n",
"disp('here,let mass of steam bled for feed heating be m kg')\n",
"disp('For process 2-6,s2=s6=6.6582 KJ/kg K')\n",
"s6=s2;\n",
"disp('Let dryness fraction at state 6 be x6')\n",
"disp('s6=sf at 8 bar+x6*sfg at 8 bar')\n",
"disp('from steam tables,hf at 8 bar=721.11 KJ/kg,vf at 8 bar=0.001115 m^3/kg,hfg at 8 bar=2048 KJ/kg,sf at 8 bar=2.0462 KJ/kg K,sfg at 8 bar=4.6166 KJ/kg K')\n",
"hf=721.11;\n",
"vf=0.001115;\n",
"hfg=2048;\n",
"sf=2.0462;\n",
"sfg=4.6166;\n",
"disp('substituting entropy values,x6=(s6-sf)/sfg')\n",
"x6=(s6-sf)/sfg\n",
"x6=0.999;//approx.\n",
"disp('h6=hf at at 8 bar+x6*hfg at 8 bar in KJ/kg')\n",
"h6=hf+x6*hfg\n",
"disp('Assuming the state of fluid leaving open feed water heater to be saturated liquid at 8 bar.h7=hf at 8 bar=721.11 KJ/kg')\n",
"h7=721.11;\n",
"disp('For process 4-5,h5=h4+v4*(8-.05)*10^2 in KJ/kg')\n",
"h5=h4+v4*(8-.05)*10^2\n",
"disp('Applying energy balance at open feed water heater,')\n",
"disp('m*h6+(1-m)*h5=1*h7')\n",
"disp('so m=(h7-h5)/(h6-h5) in kg')\n",
"m=(h7-h5)/(h6-h5)\n",
"disp('For process 7-1,h1=h7+v7*(200-8)*10^2 in KJ/kg')\n",
"disp('here h7=hf at 8 bar,v7=vf at 8 bar')\n",
"h7=hf;\n",
"v7=vf;\n",
"h1=h7+v7*(200-8)*10^2\n",
"disp('Thermal efficiency of cycle=((h2-h6)+(1-m)*(h6-h3)-{(1-m)*(h5-h4)+(h1-h7)})/(h2-h1)')\n",
"((h2-h6)+(1-m)*(h6-h3)-{(1-m)*(h5-h4)+(h1-h7)})/(h2-h1)\n",
"disp('in percentage')\n",
"((h2-h6)+(1-m)*(h6-h3)-{(1-m)*(h5-h4)+(h1-h7)})*100/(h2-h1)\n",
"disp('case (c) When there are two feed water heaters working at 40 bar and 4 bar')\n",
"disp('here, let us assume the mass of steam at 40 bar,4 bar to be m1 kg and m2 kg respectively.')\n",
"disp('2-10-9-3,s2=s10=s9=s3=6.6582 KJ/kg K')\n",
"s3=s2;\n",
"s9=s3;\n",
"s10=s9;\n",
"disp('At state 10.s10>sg at 40 bar(6.0701 KJ/kg K)so state 10 lies in superheated region at 40 bar pressure.')\n",
"disp('From steam table by interpolation,T10=370.6oc,so h10=3141.81 KJ/kg')\n",
"T10=370.6;\n",
"h10=3141.81;\n",
"disp('Let dryness fraction at state 9 be x9 so,') \n",
"disp('s9=6.6582=sf at 4 bar+x9*sfg at 4 bar')\n",
"disp('from steam tables,at 4 bar,sf=1.7766 KJ/kg K,sfg=5.1193 KJ/kg K')\n",
"sf=1.7766;\n",
"sfg=5.1193;\n",
"disp('x9=(s9-sf)/sfg')\n",
"x9=(s9-sf)/sfg\n",
"x9=0.9536;//approx.\n",
"disp('h9=hf at 4 bar+x9*hfg at 4 bar in KJ/kg')\n",
"disp('from steam tables,at 4 bar,hf=604.74 KJ/kg,hfg=2133.8 KJ/kg')\n",
"hf=604.74;\n",
"hfg=2133.8;\n",
"h9=hf+x9*hfg \n",
"disp('Assuming the state of fluid leaving open feed water heater to be saturated liquid at respective pressures i.e.')\n",
"disp('h11=hf at 4 bar=604.74 KJ/kg,v11=0.001084 m^3/kg=vf at 4 bar')\n",
"h11=604.74;\n",
"v11=0.001084;\n",
"disp('h13=hf at 40 bar=1087.31 KJ/kg,v13=0.001252 m^3/kg=vf at40 bar')\n",
"h13=1087.31;\n",
"v13=0.001252;\n",
"disp('For process 4-8,i.e in CEP.')\n",
"disp('h8=h4+v4*(4-0.05)*10^2 in KJ/kg')\n",
"h8=h4+v4*(4-0.05)*10^2\n",
"disp('For process 11-12,i.e in FP2,')\n",
"disp('h12=h11+v11*(40-4)*10^2 in KJ/kg')\n",
"h12=h11+v11*(40-4)*10^2\n",
"disp('For process 13-1_a i.e. in FP1,h1_a=h13+v13*(200-40)*10^2 in KJ/kg')\n",
"h1_a=h13+v13*(200-40)*10^2\n",
"disp('m1*3141.81+(1-m1)*608.64=1087.31')\n",
"disp('so m1=(1087.31-608.64)/(3141.81-608.64)in kg')\n",
"m1=(1087.31-608.64)/(3141.81-608.64)\n",
"disp('Applying energy balance  on open feed water heater 1 (OFWH1)')\n",
"disp('m1*h10+(1-m1)*h12)=1*h13')\n",
"disp('so m1=(h13-h12)/(h10-h12) in kg')\n",
"m1=(h13-h12)/(h10-h12)\n",
"disp('Applying energy balance on open feed water heater 2 (OFWH2)')\n",
"disp('m2*h9+(1-m1-m2)*h8=(1-m1)*h11')\n",
"disp('so m2=(1-m1)*(h11-h8)/(h9-h8) in kg')\n",
"m2=(1-m1)*(h11-h8)/(h9-h8)\n",
"disp('Thermal efficiency of cycle,n=[{(h2-h10)+(1-m1)*(h10-h9)+(1-m1-m2)*(h9-h3)}-{W_CEP+W_FP1+W_FP2}]/(h2-h1_a)')\n",
"disp('W_CEP=(1-m1-m2)*(h8-h4) in KJ/kg steam from boiler')\n",
"W_CEP=(1-m1-m2)*(h8-h4)\n",
"disp('W_FP1=(h1_a-h13)in KJ/kg of steam from boiler')\n",
"W_FP1=(h1_a-h13)\n",
"disp('W_FP2=(1-m1)*(h12-h11)in KJ/kg of steam from boiler')\n",
"W_FP2=(1-m1)*(h12-h11)\n",
"disp('W_CEP+W_FP1+W_FP2 in KJ/kg of steam from boiler')\n",
"W_CEP+W_FP1+W_FP2\n",
"disp('n=[{(h2-h10)+(1-m1)*(h10-h9)+(1-m1-m2)*(h9-h3)}-{W_CEP+W_FP1+W_FP2}]/(h2-h1_a)')\n",
"n=[{(h2-h10)+(1-m1)*(h10-h9)+(1-m1-m2)*(h9-h3)}-{W_CEP+W_FP1+W_FP2}]/(h2-h1_a)\n",
"disp('in percentage')\n",
"n=n*100\n",
"disp('so cycle thermal efficiency,na=46.18%')\n",
"disp('nb=49.76%')\n",
"disp('nc=51.37%')\n",
"disp('hence it is obvious that efficiency increases with increase in number of feed heaters.')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.7: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_7.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 7')\n",
"disp('from steam tables,')\n",
"disp('h2=h at 50 bar,500oc=3433.8 KJ/kg,s2=s at 50 bar,500oc=6.9759 KJ/kg K')\n",
"h2=3433.8;\n",
"s2=6.9759;\n",
"disp('s3=s2=6.9759 KJ/kg K')\n",
"s3=s2;\n",
"disp('by interpolation from steam tables,')\n",
"disp('T3=183.14oc at 5 bar,h3=2818.03 KJ/kg,h4= h at 5 bar,400oc=3271.9 KJ/kg,s4= s at 5 bar,400oc=7.7938 KJ/kg K')\n",
"T3=183.14;\n",
"h3=2818.03;\n",
"h4=3271.9;\n",
"s4=7.7938;\n",
"disp('for expansion process 4-5,s4=s5=7.7938 KJ/kg K')\n",
"s5=s4;\n",
"disp('let dryness fraction at state 5 be x5')\n",
"disp('s5=sf at 0.05 bar+x5*sfg at 0.05 bar')\n",
"disp('from steam tables,at 0.05 bar,sf=0.4764 KJ/kg K,sfg=7.9187 KJ/kg K')\n",
"sf=0.4764;\n",
"sfg=7.9187;\n",
"disp('so x5=(s5-sf)/sfg')\n",
"x5=(s5-sf)/sfg\n",
"x5=0.924;//approx.\n",
"disp('h5=hf at 0.05 bar+x5*hfg at 0.05 bar in KJ/kg')\n",
"disp('from steam tables,hf at 0.05 bar=137.82 KJ/kg,hfg at 0.05 bar=2423.7 KJ/kg')\n",
"hf=137.82;\n",
"hfg=2423.7;\n",
"h5=hf+x5*hfg \n",
"disp('h6=hf at 0.05 bar=137.82 KJ/kg')\n",
"h6=137.82;\n",
"disp('v6=vf at 0.05 bar=0.001005 m^3/kg')\n",
"v6=0.001005;\n",
"disp('for process 6-1 in feed pump,h1=h6+v6*(p1-p6) in KJ/kg')\n",
"p1=50;//steam generation pressure in bar\n",
"p6=0.05;//steam entering temperature in turbine in bar\n",
"h1=h6+v6*(p1-p6)*100\n",
"disp('cycle efficiency=W_net/Q_add')\n",
"disp('Wt=(h2-h3)+(h4-h5) in KJ/kg')\n",
"Wt=(h2-h3)+(h4-h5)\n",
"disp('W_pump=(h1-h6)in KJ/kg')\n",
"W_pump=(h1-h6)\n",
"disp('W_net=Wt-W_pump in KJ/kg')\n",
"W_net=Wt-W_pump\n",
"disp('Q_add=(h2-h1) in KJ/kg')\n",
"Q_add=(h2-h1)\n",
"disp('cycle efficiency=')\n",
"W_net/Q_add\n",
"disp('in percentage')\n",
"W_net*100/Q_add\n",
"disp('we know ,1 hp=0.7457 KW')\n",
"disp('specific steam consumption=0.7457*3600/W_net in kg/hp hr')\n",
"0.7457*3600/W_net\n",
"disp('work ratio=net work/positive work=W_net/Wt')\n",
"W_net/Wt\n",
"disp('so cycle efficiency=45.74%,specific steam consumption =1.78 kg/hp hr,work ratio=0.9967')\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
""
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.8: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_8.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 8')\n",
"T_cond=115;//condensate temperature in degree celcius\n",
"Cp=4.18;//specific heat at constant pressure in KJ/kg K\n",
"P=30*10^3;//actual alternator output in KW\n",
"n_boiler=0.9;//boiler efficiency\n",
"n_alternator=0.98;//alternator efficiency\n",
"disp('from steam tables,at state 2,h2=3301.8 KJ/kg,s2=6.7193 KJ/kg K')\n",
"h2=3301.8;\n",
"s2=6.7193;\n",
"disp('h5=hf at 0.05 bar=137.82 KJ/kg,v5= vf at 0.05 bar=0.001005 m^3/kg')\n",
"h5=137.82;\n",
"v5=0.001005;\n",
"disp('Let mass of steam bled for feed heating be m kg/kg of steam generated in boiler.Let us also assume that condensate leaves closed feed water heater as saturated liquid i.e')\n",
"disp('h8=hf at 3 bar=561.47 KJ/kg')\n",
"h8=561.47;\n",
"disp('for process 2-3-4,s2=s3=s4=6.7193 KJ/kg K')\n",
"s3=s2;\n",
"s4=s3;\n",
"disp('Let dryness fraction at state 3 and state 4 be x3 and x4 respectively.')\n",
"disp('s3=6.7193=sf at 3 bar+x3* sfg at 3 bar')\n",
"disp('from steam tables,sf=1.6718 KJ/kg K,sfg=5.3201 KJ/kg K')\n",
"sf_3bar=1.6718;\n",
"sfg_3bar=5.3201;\n",
"disp('so x3=(s3-sf_3bar)/sfg_3bar')\n",
"x3=(s3-sf_3bar)/sfg_3bar\n",
"x3=0.949;//approx.\n",
"disp('s4=6.7193=sf at 0.05 bar+x4* sfg at 0.05 bar')\n",
"disp('from steam tables,at 0.05 bar,sf=0.4764 KJ/kg K,sfg=7.9187 KJ/kg K')\n",
"sf=0.4764;\n",
"sfg=7.9187;\n",
"disp('so x4=(s4-sf)/sfg')\n",
"x4=(s4-sf)/sfg\n",
"x4=0.788;//approx.\n",
"disp('thus,h3=hf at 3 bar+x3* hfg at 3 bar in KJ/kg')\n",
"disp('here  from steam tables,at 3 bar,hf_3bar=561.47 KJ/kg,hfg_3bar=2163.8 KJ/kg K')\n",
"hf_3bar=561.47;\n",
"hfg_3bar=2163.8;\n",
"h3=hf_3bar+x3*hfg_3bar \n",
"disp('h4=hf at 0.05 bar+x4*hfg at 0.05 bar in KJ/kg')\n",
"disp('from steam tables,at 0.05 bar,hf=137.82 KJ/kg,hfg=2423.7 KJ/kg')\n",
"hf=137.82;\n",
"hfg=2423.7;\n",
"h4=hf+x4*hfg\n",
"disp('assuming process across trap to be of throttling type so,h8=h9=561.47 KJ/kg.Assuming v5=v6,')\n",
"h9=h8;\n",
"v6=v5;\n",
"disp('pumping work=(h7-h6)=v5*(p1-p5)in KJ/kg')\n",
"p1=60;//pressure of steam in high pressure turbine in bar\n",
"p5=0.05;//pressure of steam in low pressure turbine in bar\n",
"v5*(p1-p5)*100\n",
"disp('for mixing process between condenser and feed pump,')\n",
"disp('(1-m)*h5+m*h9=1*h6')\n",
"disp('h6=m(h9-h5)+h5')\n",
"disp('we get,h6=137.82+m*423.65')\n",
"disp('therefore h7=h6+6.02=143.84+m*423.65')\n",
"disp('Applying energy balance at closed feed water heater;')\n",
"disp('m*h3+(1-m)*h7=m*h8+(Cp*T_cond)')\n",
"disp('so (m*2614.92)+(1-m)*(143.84+m*423.65)=m*561.47+480.7')\n",
"disp('so m=0.144 kg')\n",
"m=0.144;\n",
"h6=137.82+m*423.65;\n",
"h7=143.84+m*423.65;\n",
"disp('steam bled for feed heating=0.144 kg/kg steam generated')\n",
"disp('The net power output,W_net=(h2-h3)+(1-m)*(h3-h4)-(1-m)*(h7-h6) in KJ/kg steam generated')\n",
"W_net=(h2-h3)+(1-m)*(h3-h4)-(1-m)*(h7-h6)\n",
"disp('mass of steam required to be generated=in kg/s')\n",
"P/(n_alternator*W_net)\n",
"disp('or in kg/hr')\n",
"26.23*3600\n",
"disp('so capacity of boiler required=94428 kg/hr')\n",
"disp('overall thermal efficiency=W_net/Q_add')\n",
"disp('here Q_add=(h2-h1)/n_boiler in KJ/kg')\n",
"Q_add=(h2-Cp*T_cond)/n_boiler \n",
"disp('overall thermal efficiency=')\n",
"W_net/Q_add\n",
"disp('in percentage')\n",
"W_net*100/Q_add\n",
"disp('so overall thermal efficiency=37.24%')"
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 8.9: Engineering_Thermodynamics_by_Onkar_Singh_Chapter_8_Example_9.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"// Display mode\n",
"mode(0);\n",
"// Display warning for floating point exception\n",
"ieee(1);\n",
"clear;\n",
"clc;\n",
"disp('Engineering Thermodynamics by Onkar Singh Chapter 8 Example 9')\n",
"P=15*10^3;//turbine output in KW\n",
"disp('At inlet to first turbine stage,h2=3230.9 KJ/kg,s2=6.9212 KJ/kg K')\n",
"h2=3230.9;\n",
"s2=6.9212;\n",
"disp('For ideal expansion process,s2=s3')\n",
"s3=s2;\n",
"disp('By interpolation,T3=190.97 degree celcius from superheated steam tables at 6 bar,h3=2829.63 KJ/kg')\n",
"T3=190.97;\n",
"h3=2829.63;\n",
"disp('actual stste at exit of first stage,h3_a=h2-0.8*(h2-h3) in KJ/kg')\n",
"h3_a=h2-0.8*(h2-h3)\n",
"disp('actual state 3_a shall be at 232.78 degree celcius,6 bar,so s3_a=7.1075 KJ/kg K')\n",
"s3_a=7.1075;\n",
"disp('for second stage,s3_a=s4;By interpolation,s4=7.1075=sf at 1 bar+x4*sfg at 1 bar')\n",
"s4=7.1075;\n",
"disp('from steam tables,at 1 bar,sf=1.3026 KJ/kg K,sfg=6.0568 KJ/kg K')\n",
"sf=1.3026;\n",
"sfg=6.0568;\n",
"disp('so x4=(s4-sf)/sfg')\n",
"x4=(s4-sf)/sfg\n",
"x4=0.958;//approx.\n",
"disp('h4=hf at 1 bar+x4*hfg at 1 bar in KJ/kg')\n",
"disp('from steam tables,at 1 bar,hf=417.46 KJ/kg,hfg=2258.0 KJ/kg')\n",
"hf=417.46;\n",
"hfg=2258.0;\n",
"h4=hf+x4*hfg\n",
"disp('actual enthalpy at exit from second stage,h4_a=h3_a-.8*(h3_a-h4) in KJ/kg')\n",
"h4_a=h3_a-.8*(h3_a-h4)\n",
"disp('actual dryness fraction,x4_a=>h4_a=hf at 1 bar+x4_a*hfg at 1 bar')\n",
"disp('so x4_a=(h4_a-hf)/hfg')\n",
"x4_a=(h4_a-hf)/hfg\n",
"disp('x4_a=0.987,actual entropy,s4_a=7.2806 KJ/kg K')\n",
"s4_a=7.2806;\n",
"disp('for third stage,s4_a=7.2806=sf at 0.075 bar+x5*sfg at 0.075 bar')\n",
"disp('from steam tables,at 0.075 bar,sf=0.5764 KJ/kg K,sfg=7.6750 KJ/kg K')\n",
"sf=0.5764;\n",
"sfg=7.6750;\n",
"disp('so x5=(s4_a-sf)/sfg')\n",
"x5=(s4_a-sf)/sfg\n",
"x5=0.8735;//approx.\n",
"disp('h5=2270.43 KJ/kg')\n",
"h5=2270.43;\n",
"disp('actual enthalpy at exit from third stage,h5_a=h4_a-0.8*(h4_a-h5)in KJ/kg')\n",
"h5_a=h4_a-0.8*(h4_a-h5)\n",
"disp('Let mass of steam bled out be m1 and m2 kg at 6 bar,1 bar respectively.')\n",
"disp('By heat balance on first closed feed water heater,(see schematic arrangement)')\n",
"disp('h11=hf at 6 bar=670.56 KJ')\n",
"h11=670.56;\n",
"disp('m1*h3_a+h10=m1*h11+4.18*150')\n",
"disp('(m1*2829.63)+h10=(m1*670.56)+627')\n",
"disp('h10+2159.07*m1=627')\n",
"disp('By heat balance on second closed feed water heater,(see schematic arrangement)')\n",
"disp('h7=hf at 1 bar=417.46 KJ/kg')\n",
"h7=417.46;\n",
"disp('m2*h4+(1-m1-m2)*4.18*38=(m1+m2)*h7+4.18*95*(1-m1-m2)')\n",
"disp('m2*2646.4+(1-m1-m2)*158.84=((m1+m2)*417.46)+(397.1*(1-m1-m2))')\n",
"disp('m2*2467.27-m1*179.2-238.26=0')\n",
"disp('heat balance at point of mixing,')\n",
"disp('h10=(m1+m2)*h8+(1-m1-m2)*4.18*95')\n",
"disp('neglecting pump work,h7=h8')\n",
"disp('h10=m2*417.46+(1-m1-m2)*397.1')\n",
"disp('substituting h10 and solving we get,m1=0.1293 kg and m2=0.1059 kg/kg of steam generated')\n",
"m1=0.1293;\n",
"m2=0.1059;\n",
"disp('Turbine output per kg of steam generated,Wt=(h2-h3_a)+(1-m1)*(h3_a-h4_a)+(1-m1-m2)*(h4_a-h5_a)in KJ/kg of steam generated')\n",
"Wt=(h2-h3_a)+(1-m1)*(h3_a-h4_a)+(1-m1-m2)*(h4_a-h5_a)\n",
"disp('Rate of steam generation required=P/Wt in kg/s')\n",
"P/Wt\n",
"disp('in kg/hr')\n",
"P*3600/Wt\n",
"disp('capacity of drain pump i.e. FP shown in layout=(m1+m2)*69192 in kg/hr')\n",
"(m1+m2)*69192\n",
"disp('so capacity of drain pump=16273.96 kg/hr')\n",
"\n",
"\n",
""
   ]
   }
],
"metadata": {
		  "kernelspec": {
		   "display_name": "Scilab",
		   "language": "scilab",
		   "name": "scilab"
		  },
		  "language_info": {
		   "file_extension": ".sce",
		   "help_links": [
			{
			 "text": "MetaKernel Magics",
			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
			}
		   ],
		   "mimetype": "text/x-octave",
		   "name": "scilab",
		   "version": "0.7.1"
		  }
		 },
		 "nbformat": 4,
		 "nbformat_minor": 0
}