L17_1005

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

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

View Full Document Right Arrow Icon
Background image of page 1

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

View Full DocumentRight Arrow Icon
Background image of page 2
Background image of page 3

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

View Full DocumentRight Arrow Icon
Background image of page 4
Background image of page 5

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

View Full DocumentRight Arrow Icon
Background image of page 6
Background image of page 7

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

View Full DocumentRight Arrow Icon
Background image of page 8
Background image of page 9

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

View Full DocumentRight Arrow Icon
Background image of page 10
Background image of page 11

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

View Full DocumentRight Arrow Icon
Background image of page 12
Background image of page 13
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);ygfilinspace(-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

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 Right Arrow Icon
Ask a homework question - tutors are online