Streamlines over an airfoil

I need help to figure out how to plot streamlines over an airfoil (only the top half, figure attached). I think i have to use the contour function or maybe the streamslice function but not quite sure. The streamlines should be going from left to right, and over the hump of the airfoil at the bottom. The freestream speed (Vinf) is 1 m/s.

Vinf = 1;
c = 1;
L = 11*c;
H = 4*c;
N = 220;
deltax = L/N;
M = 40;
deltay = H/M;

yi=linspace(0,L,M);
xi=linspace(0,H,N);
z=zeros(M,N);

t = 0.12;
p0= 0.2969;
p1=-0.1260;
p2=-0.3516;
p3= 0.2843;
p4=-0.1015;
X0 = 0:0.001:1;
Y0 = (t/0.2).*(p0*sqrt(X0)+p1*(X0)+p2*((X0).^2)+p3*((X0).^3)+p4*((X0).^4));
 
figure(1)
plot (X0+5,Y0,'r','LineWidth',1.5)
xlim([0 L])
ylim([0 H])

[darova]

Do you have a vector field of a stream?

Answer 1

[Sulaymon Eshkabilov]

Hi,

Here is a good script to simulate Joukowski airfoil:

https://www.mathworks.com/matlabcentral/fileexchange/8870-joukowski-airfoil-transformation