%% Exercise sur la détermination des courbes couple-vitesse (C-Ω) et couple puissance (P-Ω) d’une turbine éolienne

clear all;
close all;
clc;

% Une turbine éolienne à axe horizontal a un rayon R=38\ldotp 2m

R=38.2; %[m]
Ad=pi*R^2; %[m^s] area Ad de la turbine

% En connaissant la courbe Cp-\lambda nous souhaitons déterminer la caractéristique couple-vitesse de rotation de l’arbre de la turbine (C-Ω) dans une plage de variation de vitesse du vent  3\frac{m}{s}\le U_{\infty } \le 12\frac{m}{s}.

Uinf=3:1:12; %[m/s]
b=length(Uinf);

% La densité de l’air est \rho =1\ldotp 205\;\frac{\textrm{kg}}{m^3 }

rho=1.205 %[kg/m^3]

% La discrétisation de la courbe Cp-Ω  est la suivante

%Cp-Lamda discretization
lambda=[1 2 3 3.5 3.75 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 10 11 12 13 14 15];
Cp=[0.015 0.06 0.16 0.27 0.33 0.379 0.437 0.483 0.506 0.523 0.537 0.541 0.538 0.535 0.532 0.527 0.515 0.49 0.449 0.397 0.322 0.242];
a=length(lambda);

plot(lambda,Cp);
xlabel('\lambda')
ylabel('Cp')
title('Courbe CP-\lambda')

%
%% Questions
% Calculer la puissance et le couple de la turbine en correspondance une vitesse du vent U_{\infty } =5\frac{m}{s} et \lambda =8
% Déterminer les courbes couple-vitesse (C-Ω) et couple puissance (P-Ω) en correspondance de la plage de variation de vitesse du vent  3\frac{m}{s}\le U_{\infty } \le 12\frac{m}{s}.