greetings of the day
I am unable to extract the velocity values around the vortex eye inside a closed contour, so a to get the circulation. Dose anyone know how to isolate the velocity contour containg the vortex eye and find the circulation it. Ihave a .mat file of the velocity vector feild. I have calculated the curl but the i think its for the whole vector feild not at the particular region of recirculation . heres my code;
clear all
close all
load('B00001_X0_70_AvgV.mat')
x=x1';
y=y1';
vy=vy';
figure
contour(x,y,vy,15,'linewidth',1,'linecolor','b')
hold on
sr=surf(x,y,vy);
sr.EdgeColor='none';
sr.FaceColor='interp';
sr.FaceAlpha=0.8;
% [c,h]=contourf(x,y,vy,[0,-30.46],'linewidth',2);
contour(x,y,vy,[0,-30.46],'linewidth',1,'linecolor','k');
qv=quiver(x,y,vx',vy,'Color','k');
qv.LineWidth=2;
qv.AutoScaleFactor=1.5;
grid off
figure
surf(x,y,curl(vx',vy))
view(2)
% [curlx,curly,cav]=curl(x,y,vy);
%%It would be of great help if i can get some leads on it.
thanks