{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

L17_1005

# L17_1005 - a;¢.4.‘:ha‘nm 1w.” g s i(9 I [email protected]_i:1i...

This preview shows pages 1–13. Sign up to view the full content.

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: _ a;¢.4.‘.., :ha‘nm 1w.” * g s i. (9%) _ __ I. iggs[@s_i:1i M‘ H N;\Pursonalx_LTLAB\Ll % visualize the flow due to a superposition of a point source and a uniform flow. do The source is at (0,0); All distances are in units of the distance from the stagnation point to the source All velocities are in units of free stream velocity o\° e\o' o\° o‘xO % calculate and plot the dividing streamline: th=linspace{.l, 2*pi—.l); r=—(th—pi)./sin(th); x=r.*cos(th); y:r.*sin(th); Plot (x! Y) xlabel{'x');ylabel('y') axis equal axis([—2 6 —4 4]} % get pressure distribution on the dividing streamline: u=l + x./r.“2; V=y./r.“2; Cp=l—(u.“2+v.“2); figure plot(x,cp) Xlabel('x'); ylabel('cp') %conventionally, it's plotted with cp increasing down % note that cp goes negative with min{cp)=—.59. % This means the max speed is 26% above free stream. % streamlines: figure 'xg=linspace(-2,6);ygﬁlinspace(-4,4); [yg,xg]=meshgrid(yg,xg); psi=yg+atan2(yg,xg); contour(xg,yg,psi,21); hold on plot(x,y,'b') Ititle('streamlines‘) %pressure contours figure title('isobars') rg=sqrt(xg.“2+yg.”2); u=l + xg./rg.“2; v=yg-/rg-“2: . .1 cp=l—(u.A2+v.“2); Cp=l-(u.A2+v.“2); contour(xg,yg,cp;[-.6:.02:l]); hold on plot(x,y,‘b') axis equal axis([-2 6 -4 4]) fill(x,y,"); %this masks what's inside the dividing streamline title('isobars') CD 0.8 0.6 0.4 0.2 20 25 30 35 streamlines —4 . _ 2 1 o 1 2 3 4 5 6 ? / K W ...
View Full Document

{[ snackBarMessage ]}

### Page1 / 13

L17_1005 - a;¢.4.‘:ha‘nm 1w.” g s i(9 I [email protected]_i:1i...

This preview shows document pages 1 - 13. Sign up to view the full document.

View Full Document
Ask a homework question - tutors are online