The implementation of the compute derivatives filter has been around a while so 
the author probably hasn't seen your email. I'm not super familiar with the 
code, but I believe the vorticity is computed from the first order derivative 
of the vector field.

It is not actually using a central difference scheme to compute the gradients. 
Instead, it is treating the field as a piecewise linear field (or piecewise 
quadratic or piecewise cubic depending on the cell type) based on the volumes 
of the cells. The derivatives are then computed in the middle of the piecewise 
C1 continuous parts of the cells from the known function.

You can take a look at the source for the algorithm here:

http://vtk.org/gitweb?p=VTK.git;a=blob;f=Filters/General/vtkCellDerivatives.cxx;h=4debdf45eeb7a11e98062ab2d9efa65e938433b6;hb=ab4dc579c340167b6b3ae6cd82b3fd645642796a

Vorticity is computed between lines 169 and 216.

-Ken

From: Jean Hertzberg <[email protected]<mailto:[email protected]>>
Date: Wednesday, February 11, 2015 at 12:01 PM
To: "[email protected]<mailto:[email protected]>" 
<[email protected]<mailto:[email protected]>>
Subject: [EXTERNAL] [Paraview] Vorticity calc algorithm?

Hi
I asked this last week, but got zero response. Please help me, or point me 
towards the right documentation:
Exactly what algorithm is used to calculate vorticity in the Compute 
Derivatives filter? Is it a first or a second order central difference scheme 
for gradients, or does it use the circulation method?

Thanks

Jean Hertzberg

Dept. Mechanical Engineering
University of Colorado, Boulder

_______________________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview

Reply via email to