%% Minicurso E2A-Aprimoramento aerodinâmico.
%
% implementação da geometria parametrizada NACA 4-digitos. Além de gerar a
% geometria baseada nos parametros da família de aerofólios, aplica-se a
% rotina LumpedVortexVec para calcular a curva de Cl Vs alpha do aerofólio.
%
% A implementação da geometria não está restrita à parametros inteiros. Os
% digitos do aerofólio podem aqui ser definidos como números reais
% positivos.
% 

close all
clear all

%% Definição da geometria
%
% O perfil Naca 4 digitos define a linha de arqueamento por dois
% parametros (M,P):
%
%  * M é o valor máximo do arqueamento em unidades de percentual
%  * P é a coordenadas em x onde ocorre o máximo arqueamento em desenas de
%    percentual
%

n = 50;
digit = [2.5 3.7 12];
[xy,x,yt,yc]=NACA4Digits(digit,n,'Plota');

alpha = pi/180*[-5:5];
xy = [x(:) yc(:)];
[Cl,Cm0,xv,gamma,dCp]=LumpedVortexVec(xy,alpha);


figure
plot(alpha*180/pi,Cl)
xlabel('alpha')
ylabel('Cl')
grid minor

% calculo do dCl/dalpha
dCldalpha = 180/pi*(Cl(3:end)-Cl(1:end-2))/2;

figure
plot(alpha(2:end-1)*180/pi,dCldalpha)
xlabel('\alpha')
ylabel('dCl/d\alpha')
grid on
axis([alpha(2) alpha(end-1) 0 2*pi+0.1])