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 [email protected]:
Send Nek5000-users mailing list submissions to [email protected]
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 [email protected]
You can reach the person managing the list at [email protected]
When replying, please edit your Subject line so it is more specific than "Re: Contents of Nek5000-users digest..."
Today's Topics:
1. Vorticity ([email protected]) 2. Re: Vorticity ([email protected]) 3. Re: Vorticity ([email protected])
----------------------------------------------------------------------
Message: 1 Date: Wed, 13 Jul 2011 17:55:35 -0400 From: [email protected] Subject: [Nek5000-users] Vorticity To: [email protected] Message-ID: <[email protected]> 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: [email protected] Subject: Re: [Nek5000-users] Vorticity To: [email protected] Message-ID: <[email protected]> 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, [email protected] 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 [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
------------------------------
Message: 3 Date: Wed, 13 Jul 2011 23:59:08 -0500 (CDT) From: [email protected] Subject: Re: [Nek5000-users] Vorticity To: [email protected] Message-ID: <[email protected]> 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: [email protected] To: [email protected] 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, [email protected] 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 [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
------------------------------
_______________________________________________ Nek5000-users mailing list [email protected] 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