[Nek5000-users] Vorticity

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

Hi,

Thanks for your quick reply. Just I wanna know if I can use the  
interpolated values of velocity at some arbitrary points in the domain  
(not necessarily the grid points) as the input arguments to theses  
subroutines.

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