[Nek5000-users] Vorticity

[nek5000-users at lists.mcs.anl.gov]

Hi Alireza,

For vorticity computation, you can call

       call comp_vort3(vort,work1,work2,vx,vy,vz)

from userchk with the appropriate definitions of vorticity and work arrays 
like

       common /scren/   vort (lx1,ly1,lz1,lelv,3) ! x y z components
      $               , work1(lx1,ly1,lz1,lelv)
      $               , work2(lx1,ly1,lz1,lelv)

Let me check what we have for the calculation of strain tensor.
Best,
Aleks


On Wed, 13 Jul 2011, nek5000-users at lists.mcs.anl.gov wrote:

> Hi,
>
> I need to calculate the vorticity vector and the rate of strain tensor at any 
> arbitrary point in the domain. I wonder if these values are computed 
> somewhere in the code or if I should interpolate them in the usr file. I 
> appreciate any kind of help.
>
> Thanks in advance,
>
>
> Regards,
> Alireza Karimi
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users