Re: [Nek5000-users] Vorticity
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
Hi Alireza, The best option is to call the vorticity/s_ij routines, then interpolate the fields that they generate. Attached is a routine that should do the interpolation --- I just hacked into another routine that was available and haven't tested, but it should work. Note that if you want to change your point list subsequent to the initial call, you should pass a handle so indicating. Interpolation requires two steps (done by the attached code): 1. Find the points (x,y,z)_i in the (r,s,t,e) parameter space (done once, and expensive) 2. Evaulate quantities of interest at (r,s,t,e)_i (potentially done often, and less expensive) Please let me know if you have any questions on usage, etc. Paul On Thu, 14 Jul 2011, [email protected] wrote:
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 _______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Hi, I need to monitor the vorticity in a single point of my domain, i've already add comp_vort3() and hpts() to userchk but only get velocity, pressure and temperature in hpts.out file. How i could modify hpts subroutine to get vorticity as output?. Best regards, Nicolás Thier
Hi Nicolas, You could look at the hemi.usr in the nekexamples directory. There you will see a routine: subroutine interp_v(uvw,xyz,n) that could be modified to get the vorticity at the point in question. If the point of interest is a grid point you could alternatively just print the value at that point --- but it's a bit of a pain to identify such points, particularly in parallel. Paul ________________________________________ From: [email protected] [[email protected]] on behalf of [email protected] [[email protected]] Sent: Tuesday, October 11, 2016 10:53 AM To: [email protected] Subject: Re: [Nek5000-users] Vorticity Hi, I need to monitor the vorticity in a single point of my domain, i've already add comp_vort3() and hpts() to userchk but only get velocity, pressure and temperature in hpts.out file. How i could modify hpts subroutine to get vorticity as output?. Best regards, Nicolás Thier _______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
participants (1)
-
nek5000-users＠lists.mcs.anl.gov