Plot the streamlines and pressure for flow around a rotating cylinder

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');

回答(1 个)

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 中查找有关 2-D and 3-D Plots 的更多信息

提问:

erq
2013-12-13

回答:

2013-12-13

Community Treasure Hunt

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

Start Hunting!

Translated by