Plot the streamlines and pressure for flow around a rotating cylinder

2013 12 13

1 个回答

3 次查看（30 天）

0 个投票

r0=1,U=1,p=1,0<theta<2pi,r0<r<3r0, generate gamma=0, gamma=-2pir0U.

LAMBDA = 1.; Pinf = 1.; U = 1.;

th = linspace(0,2*pi,100); r = linspace(0.1,1,100); [TH,R] = meshgrid(th,r); [X,Y] = pol2cart(TH,R);

PSI = -LAMBDA.*sin(TH)./R; VR = -LAMBDA./(R.^2).*cos(TH); VTH = -LAMBDA./(R.^2).*sin(TH); P = Pinf + U^2/2 - (VR.^2 + VTH.^2)/2;

figure contourf(X,Y,P,20,'LineColor','none'); hold on; contour(X,Y,PSI,20,'-k.','LineWidth',1); axis equal; xlabel('X'); ylabel('Y');

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

Simon 2013-12-13

0 个投票

Hi!

Please use the code formatting for pasting code, thanks!

And what is the question? What did you do? Where does the code come from? What does it calculate?

在帮助中心和 File Exchange 中查找有关 Vector Fields 的更多信息

Find the treasures in MATLAB Central and discover how the community can help you!

Translated by