[Nek5000-users] Vorticity

Thu Jul 14 18:00:29 CDT 2011

Hi Alireza,

I'm not certain of the context of your question, but you are
correct that you will have to repack sij() from a

    (lx1*ly1*lz1,6,lelt)

layout to a set of 6 vectors of the form:

    (lx1*ly1*lz1,lelt)

This will be best done with a loop of the form:

       nxyz=nx1*ny1*nz1

       do e=1,nelv
       do i=1,n
          v1(i,1,1,e) = sij(i,1,1,1,e)
          v2(i,1,1,e) = sij(i,1,1,2,e)
          v3(i,1,1,e) = sij(i,1,1,3,e)
          v4(i,1,1,e) = sij(i,1,1,4,e)
          v5(i,1,1,e) = sij(i,1,1,5,e)
          v6(i,1,1,e) = sij(i,1,1,6,e)
       enddo
       enddo

then call interp() with (v1,v2,v3...) and (v4,v5,v6).

You should also check the sij code to understand which parts of
the tensor are stored where and whether the computed sij includes
the 0.5 factor.  (I think in this case that it does not.)

hth,

Paul

On Thu, 14 Jul 2011, nek5000-users at lists.mcs.anl.gov wrote:

> Hi,
>
> Thanks for the code. I just have a problem regarding s_ij. Should I use 
> opcopy subroutine twice as follows?
>
> call opcopy(wrk(1,1),wrk(1,2),wrk(1,3),v1,v2,v3)
> call opcopy(wrk(1,4),wrk(1,5),wrk(1,6),v4,v5,v6)
>
> Or I should use another subroutine in this case in order to pack working 
> array?
>
> Regards,
> Alireza
>
>
>
>
>
> Quoting nek5000-users-request at lists.mcs.anl.gov:
>
>> Send Nek5000-users mailing list submissions to
>> 	nek5000-users at lists.mcs.anl.gov
>> 
>> To subscribe or unsubscribe via the World Wide Web, visit
>> 	https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>> or, via email, send a message with subject or body 'help' to
>> 	nek5000-users-request at lists.mcs.anl.gov
>> 
>> You can reach the person managing the list at
>> 	nek5000-users-owner at lists.mcs.anl.gov
>> 
>> When replying, please edit your Subject line so it is more specific
>> than "Re: Contents of Nek5000-users digest..."
>> 
>> 
>> Today's Topics:
>>
>>   1. Vorticity (nek5000-users at lists.mcs.anl.gov)
>>   2. Re: Vorticity (nek5000-users at lists.mcs.anl.gov)
>>   3. Re: Vorticity (nek5000-users at lists.mcs.anl.gov)
>> 
>> 
>> ----------------------------------------------------------------------
>> 
>> Message: 1
>> Date: Wed, 13 Jul 2011 17:55:35 -0400
>> From: nek5000-users at lists.mcs.anl.gov
>> Subject: [Nek5000-users] Vorticity
>> To: nek5000-users at lists.mcs.anl.gov
>> Message-ID: <20110713175535.1042430oct81vz2v at webmail.vt.edu>
>> Content-Type: text/plain; charset=ISO-8859-1; DelSp="Yes";
>> 	format="flowed"
>> 
>> 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
>> 
>> 
>> 
>> ------------------------------
>> 
>> Message: 2
>> Date: Wed, 13 Jul 2011 17:09:48 -0500 (CDT)
>> From: nek5000-users at lists.mcs.anl.gov
>> Subject: Re: [Nek5000-users] Vorticity
>> To: nek5000-users at lists.mcs.anl.gov
>> Message-ID: <Pine.LNX.4.64.1107131704140.26989 at jacobi.mcs.anl.gov>
>> Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed
>> 
>> 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
>> 
>> 
>> ------------------------------
>> 
>> Message: 3
>> Date: Wed, 13 Jul 2011 23:59:08 -0500 (CDT)
>> From: nek5000-users at lists.mcs.anl.gov
>> Subject: Re: [Nek5000-users] Vorticity
>> To: nek5000-users at lists.mcs.anl.gov
>> Message-ID:
>> 	<600285192.94128.1310619548370.JavaMail.root at zimbra.anl.gov>
>> Content-Type: text/plain; charset=utf-8
>> 
>> 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
>> 
>> 
>> ------------------------------
>> 
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>> 
>> 
>> End of Nek5000-users Digest, Vol 29, Issue 2
>> ********************************************
>> 
>
>
>
> -- 
> Alireza Karimi
> PhD Candidate
> Department of Engineering Science and Mechanics
> Virginia Polytechnic Institute and State University
> Blacksburg, VA 24061
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users