我试图在MATLAB中使用Sullivan涡模型在圆柱坐标系下求解轴对称不可压缩Navier Stokes方程。
以下是示例代码：
drive.google.com/open?id=0B-rqPMNkwWbNV2NBOURhUzZEcWM
如果r表示径向距离，z表示轴向距离。
你好，v＆amp; w分别是r，theta和z方向的速度，仅取决于'r'和'z'，不依赖于θ，因此是2-D矩阵。p是压力，仅取决于'r'和'z'，不依赖于θ，因此是2-D矩阵。
U，V，W＆amp; P分别是r，theta和z方向的校正数据，仅取决于'r'和'z'，但是也包括3-D矩阵，因为θ方向也包括在内，但参数在θ方向上不变。
我尝试了以下但没有一个有用。
[th,R]=meshgrid(theta,r);
[x1,y1] = pol2cart(R,th);
contour(x1,y1,P(:,:,34),500)
%xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
[x1,y1] = pol2cart(th,R);
contour(x1,y1,P(:,:,34),500)
quiver(x1,y1,U(:,:,34),V(:,:,34))
%xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
contour(R.*cos(th),R.*sin(th),P(:,:,34),500)
现在，我正在尝试根据这些数据绘制轮廓图和箭头图，但这些图也会随着θ而变化。在其他一些情况下，我的圈子不完整。
也许，我无法为速度矢量指定正确的方向。
请告诉我一些合适的方法。