Lift

on 12 Oct 2021

th=.5;

s=.5+i*.1;

% TRANS

r=1.2;

lm=r-s;

% CIRCULATION

k=2*sin(th);

G=k/(2*pi); %CIRCU

w=exp(i*th);

x = meshgrid(-5:.1:5); %mesh

y = x';

z = x+i*y;

f=w*(z)+(exp(-i*th)*r^2)./(z-s)+i*k*log(z);

J=z+lm^2./z;

a= 0:.1:2*pi;

zc=r*(cos(a)+i*sin(a))+s;

zf=zc+lm^2./zc;

L=300*1.2*G;

hold

Current plot held

contour(real(J),imag(J),imag(f),[-5:.1:5])

fill(real(zf),imag(zf),'y')

%axis equal

%axis([-5 5 -5 5])

%title(strcat('Flow Around the Corresponding Airfoil. Lift: ',L_str,' [N/m]'));