Wind Analysis

clc
clear
M=30;
g=9.81;
MAX_VEL=1;
for i=1:M;
Lamda1(i)=2+90*rand(1,1)*exp(-0.001*M);
a(i)=1;
Time_Scale(i)=a(i)*Lamda1(i);
end

for i=1:M;
DT(i)=0.0098+randn(1,1)*exp(-0.01*M);
TV(i)=10;
DZ(i)=5;
TVV(i)=DT(i)/(TV(i)*DZ(i));
end
for i=1:M;
NB(i)=(g/TVV(i))^0.5;
Lamda(i)=2*pi*MAX_VEL*0.7/NB(i);
end

for i=1:M;
NB(i)=(g/TVV(i))^0.5;
Lamda(i)=2*pi*MAX_VEL*0.7/NB(i);
end
R1=625.856;
R2=181.415;
R_AV=(R1+R2)/2;
L=2*pi*R_AV;
MID_L=L/2;
DX=MID_L/M;
DTHETA=pi/M;
THETA(1)=0;
x(1)=R_AV*cos(THETA(1));
for i=1:M-1;
    THETA(i+1)=THETA(i)+DTHETA;
    x(i+1)=R_AV*cos(THETA(i+1));
end
b=0.01;  %damping factor
for i=1:M;
z1(i)=5;
z(i)=z1(i)*exp(-x(i)/(b*Lamda(i)))*cos((2*pi*x(i))/Lamda(i));
end

for i=1:M;
w(i)=10;
FR(i)=Lamda(i)/w(i);
end
figure(1)
plot(x,z)
xlabel('x')
ylabel('z')
grid on
figure(2)
plot(x,FR)
xlabel('x')
ylabel('FR')
grid on
hold on
plot(x,DT,'r')
figure(3)
plot((1:1:M),Lamda1)
xlabel('x')
ylabel('Lamda1')
grid on
pause(10)
close all
alpha=2.02;
DM=M/M;
M0=1;
for i=1:M;
    M(i)=2*rand(1,1)*exp(-i);
    P(i)=(alpha*DM*((M(i))^(alpha-1))/((M0)^(alpha)))*exp((M(i)/M0)^(alpha));
end
bar(P)
grid on

---

Unless otherwise noted, all content on this site is @Copyright by Ahmed Al Makky 2012-2013 - http://cfd2012.com