<div dir="ltr"><div><div><div>Hi Samuel,<br><br></div>The issue with the symmetric tensor magnitude calculation is that ParaView doesn't know it's a symmetric tensor so it doesn't know that it needs to sum the off-diagonal values twice. All that is know for these arrays are how many components there are. A three-component array could be representing either a 3D vector or a 3D anti-symmetrix, second order tensor so it's the user's responsibility to compute the magnitude properly where "shortened" representations of tensors are used.<br><br></div>Regards,<br></div>Andy<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, May 23, 2015 at 12:03 PM, Samuel Key <span dir="ltr"><<a href="mailto:samuelkey@bresnan.net" target="_blank">samuelkey@bresnan.net</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
Reference: ParaView Version 4.3.1<br>
<br>
I have two technical questions:<br>
<br>
(Q1) In the PV filter 'Gradient Of Unstructured Data Set' is the
"QCriterion" function, the Q-Criterion defined by J. C. R. Hunt, A.
Wray, and P. Moin (1988) "Eddies, stream, and convergence zones in
turbulent flows," Center for Turbulence Research Report CTR-S88, and
referenced by G. Haller in "<i>An objective definition of a vortex</i>,"
J. Fluid Mech. (2005), vol. 525, pp.1–26?<br>
<br>
(Answer to Q1) Courtesy of Don Simons of Egility Corp., Yes, the
formula being<br>
Q=-0.5*(Vx,x*Vx,x + Vy,y*Vy,y + Vz,z*Vz,z) - (Vx,y*Vy,x + Vx,z*Vz,x
+ Vy,z*Vz,y)<br>
<br>
(Q2) In the PV filter 'Calculator' what is the formula used by the
"mag" function to calculate mag(<b>S</b>) where <b>S</b> is a
symmetric 2nd-order tensor? Using matrix notation, the tensor <b>S</b>
would have the following components:<br>
<br>
---------------<br>
| Sxx Sxy Sxz |<br>
| Sxy Syy Syz |<br>
| Sxz Syz Szz |<br>
---------------<br>
<br>
(Answer to Q2) The PV formula for mag(<b>S</b>) is the square root
of the sum of the squares of the 6-tuple (Sxx, Syy, Szz, Sxy, Sxz,
Syz). <br>
<br>
For consistency with a three-component vector, I would have expected
it to be the 6-tuple (Sxx, Syy, Szz, 2Sxy, 2Sxz, 2Syz)?<br>
<br>
<br>
Samuel W Key<br>
FMA Development, LLC<br>
1005 39th Ave NE<br>
Great Falls, Montana 59404<br>
USA
</div>
<br>_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" target="_blank">http://markmail.org/search/?q=ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/paraview" target="_blank">http://public.kitware.com/mailman/listinfo/paraview</a><br>
<br></blockquote></div><br></div>