3 files changed, 84 insertions, 0 deletions
diff --git a/3434/CH7/EX7.1/Ex7_1.sce b/3434/CH7/EX7.1/Ex7_1.sce
new file mode 100644
index 000000000..ac89e4d55
--- /dev/null
+++ b/3434/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,42 @@
+// given data
+clear 
+clc
+rho=1.226 // air density in kG/m^3
+alpha =0.14
+H=10.0 // height at which wind speed is given in m
+uH=12.0 // speed in m/s
+z=100.0 // tower height in m
+D=80.0 // diameter  in m
+effigen=0.85 // efficiency og generator
+
+A=%pi*(D**2)/4 // area in m^3
+u0=uH*(z/H)**alpha // velocity at 100 m in m/s
+u1=0.8*u0 // exit velocity in m/s 
+Po=(A*rho*u0**3)/2 // Total Power in Wind
+// Part 1 
+printf("Total Power in Wind is %0.2f MW \n",Po/1000000)
+
+// Part 2
+a=(u0-u1)/u0 // interference factor
+Cp=4*a*(1-a)**2 // Power Coefficient
+PT=Cp*Po/1000000 // power to turbine in MW
+
+printf("The power extracted by turbine is %0.2f MW \n",PT)
+
+// Part 3
+Pelec=effigen*PT // electrical power generated in MW
+
+printf("The Electrical power generated is %0.2f MW \n",Pelec)
+
+// Part 4
+FA=4*a*(1-a)*(A*rho*u0**2)/2 // axial thrust in N
+
+printf("The axial thrust is %0.2f N \n",FA)
+
+// Part 5
+
+Fmax=(A*rho*u0**2)/2 // maximum thrust in N
+printf("Maximum axial thrust is %0.2f N \n ",Fmax)
+
+
+
diff --git a/3434/CH7/EX7.2/Ex7_2.sce b/3434/CH7/EX7.2/Ex7_2.sce
new file mode 100644
index 000000000..dc0d2768e
--- /dev/null
+++ b/3434/CH7/EX7.2/Ex7_2.sce
@@ -0,0 +1,27 @@
+// given data
+clear 
+clc
+
+u0=20.0 // wind speed in m/s
+T=273+27.0 // temp in kelvin
+P=1.01325e5 // pressure in Pa
+R=287.0 // gas constant
+r=80/2.0 // radius of rotor in m
+w=2*%pi*40/60.0 // rotor speed in rad/s
+A=%pi*r**2 // area of rotor in m^2
+
+// soln:
+rho=P/(R*T) // density in Kg/m^3
+a=1/3.0 // condition for maximum output
+Cpmax=4*a*(1-a)**2 // Power Coefficient
+Lambda=r*w/u0 //tip speed ratio
+
+Po=(A*rho*u0**3)/2000000 // Total Power in Wind in W
+
+Tm=Po*r/u0 // Torque in N
+
+Ctmax=Cpmax/Lambda // torque coefficient 
+
+Tshmax=Tm*Ctmax // torque at shaft
+
+printf("The torque at shaft for maximum output is %0.2f N",Tshmax)
diff --git a/3434/CH7/EX7.3/Ex7_3.sce b/3434/CH7/EX7.3/Ex7_3.sce
new file mode 100644
index 000000000..505c50461
--- /dev/null
+++ b/3434/CH7/EX7.3/Ex7_3.sce
@@ -0,0 +1,15 @@
+// Given Data
+// given data
+clear();
+clc();
+u0=15.0 // wind speed in m/s
+R=80/2.0 // radius of rotor in m
+n=3 // number of blades
+
+Lambda=4*%pi/n // condition of tip ratio for maximum output
+
+w=Lambda*u0/R // using Eq 7.21 rotor speed in rad/s
+
+N=w*60/(2*%pi) // rotor speed in RPM
+
+printf( "For optimum energy the rotor speed should be %.1f rpm",N)