清空筛选条件
清空筛选条件

How to calculate vorticity and rotation of the vector field

70 次查看（过去 30 天）
Petr Michalek
Petr Michalek 2021-9-2
评论： Bjorn Gustavsson 2021-9-3
Hello,
I need to calculate the vorticity and rotation of the vector field with the curl function, but I get only Infs and NaNs results. I have 4000 snapshots of a 2D flow field, each snapshot is 159x99 vectors, containts x and y coordinates in mm and U and V components in m/s. The x and y variables are 159x99 double, the Udatar and Vdatar variables are 4000x159x99 double. The included MAT file contains the first snapshot of the flow field.
For the beginning I want to calculate the curl for the first snapshot. Thanks for your help.
My code:
close all; clear;
input_folder=uigetdir;
cd(input_folder);
load('PIVkolmo_example.mat');
I=159; J=99; IJ=[159 99];
% Ures1=reshape(Udatar(1,:,:),I,J);
% Vres1=reshape(Vdatar(1,:,:),I,J);
[curlz1,cav1] = curl(x,y,Ures1,Vres1);
  8 个评论
显示 6更早的评论
Bjorn Gustavsson
Bjorn Gustavsson 2021-9-3
From the help of curl:
[CURLZ, CAV] = curl(U,V) assumes
[X Y] = meshgrid(1:N, 1:M) where [M,N]=SIZE(U).
Since curl is a differential operator you simply have to divide by dx and dy, in the correct places...

请先登录，再进行评论。

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 MATLAB 的更多信息

标签

产品


版本

R2020b

Community Treasure Hunt

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

Start Hunting!

Translated by