summaryrefslogtreecommitdiff
path: root/3785/CH5/EX5.13/Ex5_13.sce
blob: e77b1d298e7a0213d8b457418f87abc7fa05b57c (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
// Example 5_13
clc;funcprot(0);
// Given data
h=100;// Height in m
A_n=1.0;// The area of the turbine jet stream in in^2
alpha=20;// The blade angle in degree
g=9.807;// The acceleration due to gravity in m/s^2
rho=1*10^3;// The density of water in kg/m^3

// Calculation
// (a)
V_n=sqrt(2*g*h);// The nozzle velocity V_n in m/s
printf("\n(a)The nozzle velocity V_n=%2.2f m/s",V_n);
// (b)
maxP_b=((1+cosd(alpha))/2)*(rho*(A_n*2.54*10^-2)^2*(V_n^3/2))/1000;// The maximum power P, of the turbine in kW
printf("\n(b)The maximum power P_t of the turbine is %2.2f kW.",maxP_b);
// (c)
V_b=V_n/2;// The blade speed in m/s
F_b=rho*(A_n*2.54*10^-2)^2*(V_n-V_b)^2*(1+cosd(alpha));// The force in N
printf("\n(c)The blade speed,V_b=%2.2f m/s \n   The force when maximum power is being produced,F_b=%3.1f N",V_b,F_b);