[arunsmec]

Hi

I am trying to find the curl of a 3D vector field using the Matlab function curl. However, I get only -Inf and Inf in the output.
The data I have is the velocity field in 3 different 3d arrays U, V and W and the grid points as one-dimensional monotonic arrays x, y and z. I do the following operation to find the curl.

[Y,X,Z]=meshgrid(y,x,z)
[Wx,Wy,Wz,omega]=curl(X,Y,Z,U,V,W)

I am unable to figure out what causes division by zero.

[FMDenaro]

curl requires first the velocity field, then the grid, see: https://it.mathworks.com/help/symbol...id=gn_loc_drop

[arunsmec]

I was referring to this function.
https://in.mathworks.com/help/matlab/ref/curl.html

Originally Posted by FMDenaro

curl requires first the velocity field, then the grid, see: https://it.mathworks.com/help/symbol...id=gn_loc_drop

The function which you have mentioned works with symbolic variables.

[FMDenaro]

Originally Posted by arunsmec

I was referring to this function.
https://in.mathworks.com/help/matlab/ref/curl.html



The function which you have mentioned works with symbolic variables.

I suggest to start first a check by using successive 2D planes. For example for several Z values you can first check if the vector plot in each plane is correct and then compute the 2D curl.

[arunsmec]

Thank you for your suggestion.
I have found out that the problem is associated with employing the meshgrid function. It has been sorted out by resizing the velocity arrays.