Curl, Divergence calculation from velocity data

8 views (last 30 days)
Hi,
I have calculate the Curl, Divergence from the velocity data file available in 4650 x 4 array. Here 1, 2, 3, and 4 columns reprsents x, y, u, v. I have tried with below script, but I am getting some error. Please help..
X = data (:,1);
Y = data (:,2);
UU = data (:,3);
VV = data (:,4);
[X,Y] = meshgrid(X,Y);
[UU, VV] = meshgrid(UU,VV);
f = curl(UU,VV);
  2 Comments
Sudheer Bhimireddy
Sudheer Bhimireddy on 6 Aug 2020
"some error" doesn't help much. Post your error so that we can help.
Turbulence Analysis
Turbulence Analysis on 7 Aug 2020
Hi, The error msg is as below
Index in position 1 is invalid. Array indices must be positive integers or logical values.

Sign in to comment.

Answers (3)

Bruno Luong
Bruno Luong on 7 Aug 2020
Edited: Bruno Luong on 7 Aug 2020
Your data is already gridded, no need to gridded on top of it
load('matlab.mat')
UU=reshape(data(:,3),[75 62]);
VV=reshape(data(:,4),[75 62]);
X=reshape(data(:,1),[75 62]);
Y=reshape(data(:,2),[75 62]);
close all
figure
quiver(X,Y,UU,VV);
f = curl(UU,VV);
figure
imagesc(f');
  1 Comment
Turbulence Analysis
Turbulence Analysis on 7 Aug 2020
Hi, Bruno,
Many thanks for your quick reply..
Actually, I am getting only vector plot. Stiil, getting below error while exceuting 'curl'
Index in position 1 is invalid. Array indices must be positive integers or logical values.

Sign in to comment.


Turbulence Analysis
Turbulence Analysis on 7 Aug 2020
Hi, Bruno,
Now, I am able to execute curl function. However, it seems, some mistake in curl computtaion. Actually, the calculated vorticity should look something similar to the attached fig.
Furthermore, I would like to get the vector plot and vorticty plot on the same fig, I tried with 'hold on', but I am not gettting correct figure..
I would be great if you could help me with this...
  5 Comments
Turbulence Analysis
Turbulence Analysis on 7 Aug 2020
Thnaks..
For instnace if I need to get only dVx/dy , how to implement this ..

Sign in to comment.


Turbulence Analysis
Turbulence Analysis on 7 Aug 2020
Hi,
Many thanks for your support..
For gradieant I have used as follows
[dx, dy]= gradient (UU');
I am getting good results
I am getting good results...

Categories

Find more on Visual Exploration in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!