소개글
[유체역학] Potential Flow로 나타낸 Air Foil 주변 유동 해석에 대한 자료입니다.
목차
Ⅰ. Air Foil
Ⅱ. Programming Algorithm
2.1 Flow Chant
Ⅲ. Program Coding
3.1 소스와측정위치의좌표
3.2 측정위치사이의단위벡터
3.4 내적을통해소스구하기
3.5 Grid를나누어속도구하기
Ⅳ. Velocity Fields
4.1 각도변화에 따른 속도장
Ⅴ. Conclusion
본문내용
Ⅲ. Program Coding
3.1 소스와측정위치의좌표
grid=100;
N=28;
Sx=[0.05, 0.058518543,...,0.05];
Sy=[0, 0.064704761, 0.125, ... , 0];
Nx=[0, 0.010222252, ... , 0];
Ny=[0, 0.077645714, ... , 0];
3.2 측정위치사이의단위벡터
for i=1:N
dx(i)=Nx(i)-Nx(i+1);
dy(i)=Ny(i)-Ny(i+1);
end
r=sqrt(dx.^2+dy.^2);
unitvector_x=dx./r;unitvector_y=dy./r;
normal_vector_x=unitvector_y;
normal_vector_y=-unitvector_x;
3.3 각소스에서중점으로의단위벡터
for i=1:N
x_m(i)=(Nx(i)+Nx(i+1))/2;
y_m(i)=(Ny(i)+Ny(i+1))/2;
end
for i=1:N
for j=1:N
dSmx(i,j)=x_m(i)-Sx(j);
dSmy(i,j)=y_m(i)-Sy(j);
end
end
Sr=sqrt(dSmx.^2+dSmy.^2);
unitvector_Sx=dSmx./Sr;unitvector_Sy=dSmy./Sr;
3.4 내적을통해소스구하기
for i=1:N
for j=1:N
A(i,j)=(normal_vector_x(i)*unitvector_Sx(i,j)+normal_vector_y(i)*unitvector_Sy(i,j))/(2*pi*Sr(i,j));
end
end
pB=-normal_vector_x;B=pB';
Q=inv(A)*B;