diff options
Diffstat (limited to 'sample_notebooks/KavinkumarD')
| -rw-r--r-- | sample_notebooks/KavinkumarD/Chapter_11__Impulse_and_Reaction_Turbines.ipynb | 345 |
1 files changed, 345 insertions, 0 deletions
diff --git a/sample_notebooks/KavinkumarD/Chapter_11__Impulse_and_Reaction_Turbines.ipynb b/sample_notebooks/KavinkumarD/Chapter_11__Impulse_and_Reaction_Turbines.ipynb new file mode 100644 index 00000000..5415ad01 --- /dev/null +++ b/sample_notebooks/KavinkumarD/Chapter_11__Impulse_and_Reaction_Turbines.ipynb @@ -0,0 +1,345 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:38f9fe4fd8a5c174c9e1dd9b5dc21976f4cdd814f7eb8fcfe0c266e278f9a77b"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter 11 : Impulse and Reaction Turbines"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.1 and Page No:454"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "p02=6; # Inlet pressure in bar\n",
+ "T02=900; # Inlet temperature in kelvin\n",
+ "p0fs=1; # Outlet pressure in bar\n",
+ "eff_isenT=0.85; # insentropic efficiency of turbine\n",
+ "alpha_2=math.radians(75); # Nozzle outlet angle in degree and conversion to radians\n",
+ "u=250; # Mean blade velocity in m/s\n",
+ "Cp=1.15*10**3; # Specific heat in J/ kg K\n",
+ "r=1.333; # Specific heat ratio\n",
+ "\n",
+ "#Calculations\n",
+ "T0fs=T02/(p02/p0fs)**((r-1)/r); # Isentropic temperature at the exit of the final stage\n",
+ "Del_Toverall=eff_isenT*(T02-T0fs); # Actual overall temperature drop\n",
+ "c2=2*u/math.sin (alpha_2); # absolute velocity\n",
+ "c3= c2*math.cos (alpha_2);# absolute velocity\n",
+ "c1=c3; # From velocity triangles\n",
+ "Del_Tstage=(c2**2-c1**2)/(2*Cp); # Stage temperature drop\n",
+ "n=Del_Toverall/Del_Tstage; # Number of stages\n",
+ "\n",
+ "#Results\n",
+ "print \"Number of stages n =\",round (n,0);\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Number of stages n = 3.0\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.2 and Page No:455"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "N=10000; # Speed of gas turbine in rpm\n",
+ "T01=700+273.15; # Total head temperature at nozzle entry in kelvin\n",
+ "P01=4.5; #Total head pressure at nozzle entry in bar\n",
+ "P02=2.6; # Outlet pressure from nozzle in bar\n",
+ "p3=1.5;# Pressure at trbine outlet annulus in bar\n",
+ "M=0.5; # Mach number at outlet\n",
+ "alpha_2=math.radians(70); # outlet nozzle angle in degrees and conversion to radians\n",
+ "D=64; # Blade mean diameter in cm\n",
+ "m=22.5; # Mass flow rate in kg/s\n",
+ "eff_T=0.99; # turbine mechanical efficiency\n",
+ "Cp=1.147; # Specific heat in kJ/kg K\n",
+ "r=1.33; # Specific heat ratio\n",
+ "fl=0.03; # frictional loss\n",
+ "R=284.6; # characteristic gas constant in J/kg K\n",
+ "\n",
+ "#Calculations\n",
+ "eff_N=1-fl; # Nozzle efficiency\n",
+ "T_02=(P02/P01)**((r-1)/r)*T01; # Isentropic temperature after expansion\n",
+ "T02=T01-eff_N*(T01-T_02); # Actual temperature after expansion\n",
+ "c2=math.sqrt (2*Cp*10**3*(T01-T02)); # Absolute velocity\n",
+ "u=(3.14*D*10**-2*N)/60; # Mean blade velocity\n",
+ "# From velocity triangles\n",
+ "wt2=c2*math.sin( (alpha_2))-u;\n",
+ "ca=c2*math.cos( (alpha_2));\n",
+ "beta_2=(math.atan((wt2)/ca));\n",
+ "T3=T02/(P02/p3)**((r-1)/r); # Assuming rotor losses are negligible\n",
+ "c3=M*math.sqrt (r*R*T3); # Absolute velocity\n",
+ "beta_3=(math.atan(u/c3));\n",
+ "ct2=c2*math.sin((alpha_2));\n",
+ "P=eff_T*m*(ct2)*u/1000; # Power developed\n",
+ "\n",
+ "#Results\n",
+ "print \"(i).\"\n",
+ "print \"\\tGas angle at entry = \",round (math.degrees(beta_2),3),\"degree\"\n",
+ "print \"\\tGas angle at exit = \",round (math.degrees(beta_3),3),\"degree\"\n",
+ "print \"(ii).\"\n",
+ "print \"\\tPower developed = \",round(P,3),\"kW (roundoff error)\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "(i).\n",
+ "\tGas angle at entry = 41.411 degree\n",
+ "\tGas angle at exit = 51.609 degree\n",
+ "(ii).\n",
+ "\tPower developed = 3680.184 kW (roundoff error)\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.3 and Page No:457"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "alpha_2=math.radians(65); # Nozzle discharge angle in degree and conversion to radians\n",
+ "c3=300; # Absolute velocity in m/s\n",
+ "alpha_3=math.radians(30); # in degrees and conversion to radians\n",
+ "\n",
+ "#Calculations\n",
+ "ca2=c3*math.cos (alpha_3); # Axial velocity\n",
+ "c2=ca2/math.cos(alpha_2); # Absolute velocity\n",
+ "# ca3=ca2=ca and equal blade angles then\n",
+ "ca=ca2;\n",
+ "beta_2=math.atan((c2*math.sin(alpha_2)+c3*math.sin(alpha_3))/(2*ca)); # Blade angle\n",
+ "beta_3=beta_2; # equal blade angles\n",
+ "u=c2*math.sin(alpha_2)-ca2*math.tan(beta_2); # Mean blade velocity\n",
+ "# From velocity triangles\n",
+ "ct2=c2*math.sin(alpha_2);\n",
+ "ct3=c3*math.sin(alpha_3);\n",
+ "WT=u*(ct2+ct3)/1000; # Work done\n",
+ "sigma=u/c2; # optimum speed ratio\n",
+ "eff_B=4*(sigma*math.sin(alpha_2)-sigma**2);\n",
+ "\n",
+ "#Results\n",
+ "print \"Blade angle = beta_2= beta_3 = \",round (math.degrees(beta_2),3),\"degree\"\n",
+ "print \"Power Produced = \",round(WT,3),\"kJ/kg (roundoff error)\"\n",
+ "print \"Blade efficiency = \",round(eff_B*100,2),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Blade angle = beta_2= beta_3 = 53.692 degree\n",
+ "Power Produced = 143.963 kJ/kg (roundoff error)\n",
+ "Blade efficiency = 76.19 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.4 and Page No:458"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "P01=7; # Pressure at inlet in bar\n",
+ "T01=300+273.15; # Temperature at inlet in kelvin\n",
+ "P02=3; # Pressure at outlet in bar\n",
+ "alpha_2=math.radians(70); # Nozzle angle in degree and conversion to radians\n",
+ "eff_N=0.9; # Isentropic efficiency of nozzle\n",
+ "WT=75; # Power Produced in kW\n",
+ "Cp=1.15; # Specific heat in kJ/kg K\n",
+ "r=1.33; # Specific heat ratio\n",
+ "\n",
+ "#Calculations\n",
+ "T_02=T01*(P02/P01)**((r-1)/r); # Isentropic temperature after expansion\n",
+ "T02=T01-eff_N*(T01-T_02); # Actual temperature after expansion\n",
+ "c2=math.sqrt (2*Cp*10**3*(T01-T02)); # Absolute velocity\n",
+ "# For optimum blade speed ratio\n",
+ "u=(c2*math.sin (alpha_2)/2); # Mean blade velocity\n",
+ "beta_2=math.atan((c2*math.sin(alpha_2)-u)/(c2*math.cos(alpha_2))); # Blade angle\n",
+ "# From velocity triangles\n",
+ "ct2=c2*math.sin(alpha_2);\n",
+ "w2=c2*math.cos(alpha_2)/math.cos(beta_2);\n",
+ "w3=w2; # Equal inlet and outlet angles\n",
+ "beta_3=54; # in degrees\n",
+ "ct3=w3*math.sin(beta_3)-u;\n",
+ "m=(WT*10**3)/(u*(ct2+ct3)); # Gas mass flow rate\n",
+ "\n",
+ "#Results\n",
+ "print \"Blade angle = \",round(math.degrees(beta_2),3),\"degree\"\n",
+ "print \"Gas Mass Flow Rate = \",round(m,3),\"kg/s\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Blade angle = 53.948 degree\n",
+ "Gas Mass Flow Rate = 4.89 kg/s\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.5 and Page No:460"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "P01=4.6; # Total head inlet pressure in bar\n",
+ "T01=700+273.15; # Total head inlet temperature in kelvin\n",
+ "P2=1.6; # Static head pressure at mean radius in bar\n",
+ "Dm_h=10; # Mean blade diameter/blade height\n",
+ "lc=0.1; # Nozzle losses coefficient\n",
+ "alpha_2=math.radians(60); # Nozzle outlet angle in degree and conversion to radians\n",
+ "Cp=1.147; # Specific heat in kJ/kg K\n",
+ "r=1.33; # Specific heat ratio\n",
+ "m=20; # Mass flow rate in kg/s\n",
+ "R=284.6; # characteristic gas constant in J/kg K\n",
+ "\n",
+ "#Calculations\n",
+ "T_2=T01*(P2/P01)**((r-1)/r); # Isentropic temperature after expansion\n",
+ "T2=(lc*T01+T_2)/(1+lc); # Actual temperature after expansion\n",
+ "c2=math.sqrt(2*Cp*10**3*(T01-T2)); # Absolute velocity\n",
+ "# From velocity triangles\n",
+ "ca=c2*math.cos(alpha_2);\n",
+ "row=P2*10**5/(R*T2); # Density of gas\n",
+ "A=m/(ca*row); # Area\n",
+ "Dm=math.sqrt (A*Dm_h/3.14); # Mean Diameter\n",
+ "h=Dm/10; # Blade height\n",
+ "rm=Dm/2; # Mean radius\n",
+ "# At root\n",
+ "r_root=(Dm-h)/2;\n",
+ "#At the tip\n",
+ "r_tip=(Dm+h)/2;\n",
+ "# Free vorte flow\n",
+ "ct_mean=c2*math.sin (alpha_2);\n",
+ "# At the root\n",
+ "ct2_root=(ct_mean*rm)/r_root;\n",
+ "alpha2_root=math.atan(ct2_root/ca);\n",
+ "c2_root=ct2_root/math.sin (alpha2_root);\n",
+ "T2_root=T01-c2_root**2/(2*Cp*10**3);\n",
+ "# At the tip\n",
+ "ct2_tip=ct_mean*rm/r_tip;\n",
+ "alpha2_tip = math.atan (ct2_tip/ca);\n",
+ "c2_tip=ct2_tip/math.sin(alpha2_tip);\n",
+ "T2_tip=T01-c2_tip**2/(2*Cp*10**3);\n",
+ "\n",
+ "#Results\n",
+ "print \"A the Root\"\n",
+ "print \"\\tGas Temperature at the root = \",round(T2_root,3),\"K\"\n",
+ "print \"\\tGas velocity at the root = \",round(c2_root,3),\"m/s\"\n",
+ "print \"\\tDischarge angle at the root = \",round(math.degrees(alpha2_root),3),\"degree\"\n",
+ "print \"\\nA the Tip\"\n",
+ "print \"\\tGas Temperature at the tip = \",round(T2_tip,3),\"K\"\n",
+ "print \"\\tGas velocity at the tip = \",round(c2_tip,3),\"m/s\"\n",
+ "print \"\\tDischarge angle at the tip = \",round(math.degrees(alpha2_tip),3),\"degree\"\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "A the Root\n",
+ "\tGas Temperature at the root = 733.345 K\n",
+ "\tGas velocity at the root = 741.696 m/s\n",
+ "\tDischarge angle at the root = 62.543 degree\n",
+ "\n",
+ "A the Tip\n",
+ "\tGas Temperature at the tip = 795.766 K\n",
+ "\tGas velocity at the tip = 637.902 m/s\n",
+ "\tDischarge angle at the tip = 57.581 degree\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |