[Nek5000-users] Vorticity

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

Hi Alireza,


For a strain stress computation you can call subroutine comp_sij

       call comp_sij(sij,nij,vx,vy,vz
      $                     ,ur,us,ut,vr,vs,vt,wr,ws,wt)


(navier5.f:2225) with nij=6 for 3D somewhere from userchk once you declare the stress component array sij and arrays with local coordinate r-s-t derivatives of velocity components

       common /csij/
      $     sij(lx1*ly1*lz1,6,lelv) ! 6 components
      $    ,ur (lx1*ly1*lz1)
...
      $    ,wt (lx1*ly1*lz1)


Best,
Aleks






----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Wednesday, July 13, 2011 5:09:48 PM
Subject: Re: [Nek5000-users] Vorticity

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
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users