Scribd
Upload
142 views3 pages

Scilab Codes for Wind Turbine Analysis

1) The document presents Scilab code to solve problems related to wind turbines. 2) For the first problem, the code calculates the total power available in wind, power extracted by the turbine, maximum torque on the turbine shaft, and maximum axial thrust on the turbine given parameters like wind speed, rotor diameter, and turbine efficiency. 3) For the second problem, the code calculates total power available in wind, power extracted by the turbine, electrical power generated, axial thrust on the turbine, and maximum axial thrust on the turbine given additional parameters like tower height and wind speed reduction.

Uploaded by

Akhil Kotadiya
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
142 views3 pages

Scilab Codes for Wind Turbine Analysis

1) The document presents Scilab code to solve problems related to wind turbines. 2) For the first problem, the code calculates the total power available in wind, power extracted by the turbine, maximum torque on the turbine shaft, and maximum axial thrust on the turbine given parameters like wind speed, rotor diameter, and turbine efficiency. 3) For the second problem, the code calculates total power available in wind, power extracted by the turbine, electrical power generated, axial thrust on the turbine, and maximum axial thrust on the turbine given additional parameters like tower height and wind speed reduction.

Uploaded by

Akhil Kotadiya
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 3

Scilab Codes

Name : Akhil Kotadiya


Q1 Atmospheric air having wind velocity of 20 m/ s has density of 1.22 kg/m3 at high
altitude where the windmill is installed. Rotor diameter is 20 m and it runs at
60 RPM. Actual efficiency of windmill is 35 %, find actual power produced, the
maximum torque and axial thrust on the turbine shaft.

Solution:

Scilab code
clc
clear all

roh = 1.22 // air density_kg/m^3

u0 = 12.0 // Wind speed_m/s


D = 20.0 // Diameter of rotor_m
N = 60 // RPM
effgen = 35 // Efficiency of the generator_%

A = (%pi/4)*(D^2) // Area swept by the rotor_m^2

// (1) Total power available in wind

P0 = (1/2)*A*roh*(u0^3)
mprintf('Total power available in the wind is = %f Watts\n',P0)

// (2) Actual power extracted by the turbine

Pt = (effgen/100)*P0
mprintf('Power extracted by the turbine is = %f Watts\n',Pt)

// (3) Torque at maximum efficieny

Tmax = (2/27)*((roh*D*(u0^3))/(N/60))
mprintf('Torque at maximum efficieny is = %f Newton-meter\n',Tmax)

// (4) Maximum axial thrust on the turbine

Fxmax = (%pi/9)*(D^2)*roh*(u0^2)
mprintf('Maximum axial thrust on the turbine is = %f Newton\n',Fxmax)
Scilab output:
Total power available in the wind is = 331148.998430 Watts
Power extracted by the turbine is = 115902.149450 Watts
Torque at maximum efficieny is = 3123.200000 Newton-meter
Maximum axial thrust on the turbine is = 24529.555439 Newton

Q2 A propeller type wind turbine has the following data:


Speed of free wind at a height of 10 m = 12 m/s
Air density = 1.226 kg /m3
α = 0.14
Height of tower = 100 m/ s
Diameter of the rotor = 80 m
Wind velocity at the turbine reduces by 20 %
Generator efficiency = 85 %

Find:

1. Total power available in the wind


2. Power extracted by the turbine
3. Electrical power generated
4. Axial thrust on the turbine
5. Maximum axial thrust on the turbine

Solution:

Scilab code
clc
clear

roh = 1.226 // air density_kg/m^3


alpha = 0.14
H = 10.0 // Height at which wind speed is given_m/s
z = 100.0 // Tower height_m
uH = 12.0 // Wind speed at given height_m/s
D = 80.0 // Diameter of rotor_m
effgen = 85 // Efficiency of the generator_%
reduce = 20 // Reduction of speed by_%

u0 = uH*((z/H)^alpha) // Wind speed at the given height_m/s


u1 = (1 - (reduce/100))*u0 // Wind speed at the rotor_m/s
A = (%pi/4)*(D^2) // Area swept by the rotor_m^2

// (1) Total power available in wind


P0 = (1/2)*A*roh*(u0^3)
mprintf('Total power available in the wind is = %f Watts\n',P0)

// (2) Power extracted by the turbine

a = (u0 - u1)/u0 // Interference factor


Cp = 4*a*(1 - a)^2 // Power coefficient

Pt = Cp*P0 // Power extracted by the turbine


mprintf('Power extracted by the turbine is = %f Watts\n',Pt)

// (3) Electrical power generated

Pe = (effgen/100)*Pt
mprintf('Electrical power generated is = %f Watts\n',Pe)

// (4) Axial thrust on the turbine

Fa = 4*a*(1 - a)*(A*roh*((u0^2)/2))
mprintf('Axial thrust on the turbine is = %f Newton\n',Fa)

// (5) Maximum axial thrust on the turbine

Famax = A*roh*((u0^2)/2)
mprintf('Maximum axial thrust on the turbine is = %f Newton\n',Famax)

Scilab output:
Total power available in the wind is = 14004708.317772 Watts
Power extracted by the turbine is = 7170410.658699 Watts
Electrical power generated is = 6094849.059895 Watts
Axial thrust on the turbine is = 541094.096840 Newton
Maximum axial thrust on the turbine is = 845459.526313 Newton

You might also like