2D Vorticity with any stencil

version 1.0.0 (2.16 KB) by Fernando Zigunov
Computes curl() with arbitrary number of points in the stencil


Updated 21 May 2019

View License

curl(X,Y,U,V) in Matlab will give you the angular velocity or vorticity of a 2D vector field (U,V) with a 3-point stencil. This implementation will compute the curl of a **UNIFORM*** grid with a 'npts' point stencil (only supports central differences). Since it uses convolution to accomplish that, the corners and edges of the grid will be incorrect.

Using for PIV post-processing, where noise can be an issue when computing derivatives.

Shoutout to Greg von Winckel for providing the stencil generator code (https://www.mathworks.com/matlabcentral/fileexchange/5269-finite-difference-weights)!

Cite As

Fernando Zigunov (2022). 2D Vorticity with any stencil (https://www.mathworks.com/matlabcentral/fileexchange/71612-2d-vorticity-with-any-stencil), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2019a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

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

Start Hunting!