May 2016 - Nek5000-users - mailman3.cels.anl.gov

Re: [Nek5000-users] Vorticity
by nek5000-users＠lists.mcs.anl.gov 11 Oct '16

11 Oct '16

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(a)lists.mcs.anl.gov: > Send Nek5000-users mailing list submissions to > nek5000-users(a)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(a)lists.mcs.anl.gov > > You can reach the person managing the list at > nek5000-users-owner(a)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(a)lists.mcs.anl.gov) > 2. Re: Vorticity (nek5000-users(a)lists.mcs.anl.gov) > 3. Re: Vorticity (nek5000-users(a)lists.mcs.anl.gov) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 13 Jul 2011 17:55:35 -0400 > From: nek5000-users(a)lists.mcs.anl.gov > Subject: [Nek5000-users] Vorticity > To: nek5000-users(a)lists.mcs.anl.gov > Message-ID: <20110713175535.1042430oct81vz2v(a)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(a)lists.mcs.anl.gov > Subject: Re: [Nek5000-users] Vorticity > To: nek5000-users(a)lists.mcs.anl.gov > Message-ID: <Pine.LNX.4.64.1107131704140.26989(a)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(a)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(a)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(a)lists.mcs.anl.gov > Subject: Re: [Nek5000-users] Vorticity > To: nek5000-users(a)lists.mcs.anl.gov > Message-ID: > <600285192.94128.1310619548370.JavaMail.root(a)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(a)lists.mcs.anl.gov > To: nek5000-users(a)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(a)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(a)lists.mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > _______________________________________________ > Nek5000-users mailing list > Nek5000-users(a)lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > > > ------------------------------ > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users(a)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

1 3

Spectral coefficients for velocity
by nek5000-users＠lists.mcs.anl.gov 01 Jun '16

01 Jun '16

Hi Neks, Priviously, the parameters p101 and p103 have been discussed: "In each direction, on each element, the solution is represented as a polynomial of degree N in the reference element (i.e., on the interval [-1,1]). For any Nth-order polynomial there is a modal expansion of the form: u(x) = \sum_k=0^N uh_k phi_k(x) where phi_0 = 1, phi_1 = x, and phi_k(x)=L_k(x) - L_{k-1} (x) for k > 1, where L_k(x)=kth Legendre polynomial." How can I obtain these coefficients "uh_k" on the several elements? Thanks in advance for any suggestion, Andrey

1 2

laplacian in axisymmetric mode (IFAXIS=T)
by nek5000-users＠lists.mcs.anl.gov 31 May '16

31 May '16

Hi neks, I have a question about how to compute the laplacian of a scalar quantity (e. g. temperature), if IFAXIS is set to true, such as ( \nabla^2(TEMP) ). Is there an easy way to do it, or someone has already faced this problem? Thanks Mirko Farano

1 0

Temperature field values in mesh 1 only
by nek5000-users＠lists.mcs.anl.gov 19 May '16

19 May '16

Hi Neks! I am running a Conjugate Heat Transfer simulation with a fluid box (velocity mesh) and two solid boxes (all boxes- temperature mesh). I am interested in computing the heat flux budget for this simulation. While doing this, on multiple occasions, I need to multiply the temperature or passive scalars throughout the fluid domain field with the velocity component values at the respective locations. I also need to calculate the gradients of these variables throughout the fluid domain. For this I use the sub-routine "gradm1" present in navier5.f . But I believe gradm1 calculates the gradients for field variables in the velocity mesh only. For all these functions, I think it's best if I could isolate an array with the values of the temperature and multiple passive scalars in the fluid domain (or velocity mesh) only, which should be a subset of the whole temperature array. Thus my question is, is there any way to get this? I am using multiple processors which complicates things. I want to create an array that contains the values of the temperature and the passive scalars at all the velocity mesh elements and this I have to do across all processors. In other words, from the standard NEK array t(lx2,ly2,lz2,lelt,ldimt) which contains the temperature and the passive scalars, I have to create an array like, say, tm1(lx1,ly1,lz1,lelv,ldimt) which would have these values. I am not sure if I have been able to explain my query clearly enough, but I would greatly appreciate if anybody could help me with this or if you could think of an easier way to calculate the gradient of the temperature and passive scalars throughout the velocity field. Let me know if you want more clarity with the question, I will try to explain better. Thanks in advance and cheers! Pradyumna

1 2

Two-way coupling and mpi calls to userf()
by nek5000-users＠lists.mcs.anl.gov 13 May '16

13 May '16

Hi Neks, I trust you are all well, I'm having some trouble with my two-way LPT code, in particular using information from the usrchk() subroutine (where all my LPTing is done) and userf() (where I want to apply a force to all meshpoints based on the number of particles in the surrounding cell). My issue is that the code doesn't count the number of particles in each cell properly, even though the xmc, xxc etc. arrays store the correct minimum and maximum values for the cell boundaries as in the usrchk() code I use these to calculate a full volume of the mesh and it is correct. I feel like my understanding on how the mpi works is lacking, so any advice on that would also be helpful - I understand that subroutine userf (ix,iy,iz,ieg) is called once per local element? So far my 2W code in the userf() looks like this: c !=== TWO WAY COUPLING ===! if(TWOWAYON .EQ. 1 .AND. istep .GT. 1) then real twowayforcex(lx1,ly1,lz1,lelv) real twowayforcey(lx1,ly1,lz1,lelv) real twowayforcez(lx1,ly1,lz1,lelv) real xmc(lx1,ly1,lz1,lelv) real xxc(lx1,ly1,lz1,lelv) real ymc(lx1,ly1,lz1,lelv) real yxc(lx1,ly1,lz1,lelv) real zmc(lx1,ly1,lz1,lelv) real zxc(lx1,ly1,lz1,lelv) real vol(lx1,ly1,lz1,lelv) real partpostw(ldim,lpart) real partftw(ldim,lpart) integer npartcell(lx1,ly1,lz1,lelv) integer TWOWAYON integer iel common /twowayr/ twowayforcex $ , twowayforcey $ , twowayforcez $ , xmc, xxc $ , ymc, yxc $ , zmc, zxc, vol $ , partpostw, partftw common /twowayi/ npartcell $ , TWOWAYON iel = gllel(ieg) if(istep .GT. 2) then do jp = 1,300000 if(partpostw(1,jp) .GE. xmc(ix,iy,iz,iel) .AND. $ partpostw(1,jp) .LE. xxc(ix,iy,iz,iel) .AND. $ partpostw(2,jp) .GE. ymc(ix,iy,iz,iel) .AND. $ partpostw(2,jp) .LE. yxc(ix,iy,iz,iel) .AND. $ partpostw(3,jp) .GE. zmc(ix,iy,iz,iel) .AND. $ partpostw(3,jp) .LE. zxc(ix,iy,iz,iel)) then npartcell(ix,iy,iz,iel) = npartcell(ix,iy,iz,iel) + 1 ffx = ffx + partftw(1,jp) ! Still need to average these ffy = ffy + partftw(2,jp) ffz = ffz + partftw(3,jp) endif enddo else ffx = 0.0 ffz = 0.0 ffy = 0.0 endif endif Thankyou for reading this, any help would be greatly appreciated!

1 0

Re: [Nek5000-users] Nek5000-users Digest, Vol 87, Issue 1
by nek5000-users＠lists.mcs.anl.gov 12 May '16

12 May '16

Thanks a lot Paul and Max the answers were very helpful since I am new to Linux and coding. ________________________________________ From: nek5000-users-bounces(a)lists.mcs.anl.gov <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf of nek5000-users-request(a)lists.mcs.anl.gov <nek5000-users-request(a)lists.mcs.anl.gov> Sent: Wednesday, May 11, 2016 1:00 PM To: nek5000-users(a)lists.mcs.anl.gov Subject: Nek5000-users Digest, Vol 87, Issue 1 Send Nek5000-users mailing list submissions to nek5000-users(a)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(a)lists.mcs.anl.gov You can reach the person managing the list at nek5000-users-owner(a)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. problem in running a case (nek5000-users(a)lists.mcs.anl.gov) 2. Re: problem in running a case (nek5000-users(a)lists.mcs.anl.gov) 3. Re: problem in running a case (nek5000-users(a)lists.mcs.anl.gov) ---------------------------------------------------------------------- Message: 1 Date: Tue, 10 May 2016 19:57:12 +0000 From: nek5000-users(a)lists.mcs.anl.gov To: "nek5000-users(a)lists.mcs.anl.gov" <nek5000-users(a)lists.mcs.anl.gov> Subject: [Nek5000-users] problem in running a case Message-ID: <mailman.21233.1462912411.24727.nek5000-users(a)lists.mcs.anl.gov> Content-Type: text/plain; charset="iso-8859-1" Hello, I am mansi patel and I am downloading nek5000 I was able to compile makenek eddy_uv but after that I am not able to run a case nekb eddy_uv Do I need to stay in eddy directory? even if I change directory I am getting same error command not found . what should I do?

1 0

problem in running a case
by nek5000-users＠lists.mcs.anl.gov 10 May '16

10 May '16

Hello, I am mansi patel and I am downloading nek5000 I was able to compile makenek eddy_uv but after that I am not able to run a case nekb eddy_uv Do I need to stay in eddy directory? even if I change directory I am getting same error command not found . what should I do?

1 2