September 2011 - 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

Error in generalev
by nek5000-users＠lists.mcs.anl.gov 02 Jul '14

02 Jul '14

Hello Neks, I've recently come across an error I haven't seen before. It occurs in the routine "generalev" at the beginning of my run. It looks like the matrix "a" and "b" being fed to the routine have an Infinity or NaN at position (0,0). Are there any ideas of what could cause this type of error? I have attached the logfile to this email. Thanks, -- Josh Camp "All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke

1 11

problem with closed internal boundaries
by nek5000-users＠lists.mcs.anl.gov 23 Jul '13

23 Jul '13

Hello Neks, I'm trying to model a flow past an obstacle in a channel. As a first step I would like to do a 2D simulation. I'm encountering some problems during the geometry data generation in the nek5000 execution : it looks like the program can only generate NaN values although genmap seems to work well on the .rea file. I can also reopen the rea file with prenek without any problems. Do you have any idea what is going wrong here? Is it problematic to have a closed internal boundary condition? Thanks in advance for any help. Best, Michael B

1 9

opdiv and multd
by nek5000-users＠lists.mcs.anl.gov 10 Feb '13

10 Feb '13

Hi all, I am having some problems using the subroutines opdiv and multd. For some reasons they do not give me correct results. Do they need any special treatment or precautions? I was trying to compute the divergence of a 3D field, and i realized that in my case the multd subroutine used by opdiv does not gives good results, while the same derivative computed from gram1 gives no problems. For example, I overwrite an analytical field : ... parameter(lt=lx1*ly1*lz1*lelt) real Res(lt) ... vx(i,j,k,e)= 1.0*sin(x)*cos(y)*cos(z) ... call print_line(vx) ... iflg = 1 call multd (RES(1),vx,rym2,sym2,tym2,2,iflg) call print_line(Res(1) ) call gradm1(work1,Res(1),work2,vx) call print_line(Res(1) ) ... if (lx2.ne.lx1) then ... STOP endif and print output a (xr,:,zr) the result (performed by print_line) . In vxi.eps i compare the result (lines) with the analytical solution( symbols). The result of multd is not correct while all the others match perfectly the analytical solution. In vxy.eps i plot only the result of multd. Note that on y i have 12 elements the same number of peaks of the solution. Does anybody have an idea on how to fix this? Thank you very much for any help or suggestion francesco

1 14

Periodic B.C and Genmap Issue
by nek5000-users＠lists.mcs.anl.gov 09 Nov '12

09 Nov '12

Hi Paul, This post is in reference to the previous post by Michael (dated Jun 22). When I try to get a map file for a rea file that has periodic vertices, I get the following error. The rea file was generated the same way as before (Gambit -> Nek through Matlab). I checked on Matlab by plotting the periodic vertices and confirmed for a few faces that they are indeed periodic. Also, as you had earlier suggested (with the new genmap), I tried changing the tolerances fro 0.09 to 0.2 and I get the same periodic mismatch error. Could you provide any leads on where things are going wrong ? Please find below the error message. Thanks. <error> Input (.rea) file name: domain Input mesh tolerance (default 0.2): NOTE: smaller is better, but generous is more forgiving for bad meshes. 0.09 reading .rea file data ... start locglob_lexico: 8 3912 31296 8.99999999999999967E-002 locglob: 1 1 31296 locglob: 2 36 31296 locglob: 3 1224 31296 locglob: 1 4851 31296 locglob: 2 4851 31296 locglob: 3 4851 31296 locglob: 1 4851 31296 locglob: 2 4851 31296 locglob: 3 4851 31296 done locglob_lexico: 4851 4851 31296 10 start periodic vtx: 3912 4851 abort: PERIODIC MISMATCH 2: 1093 1 P ie 0 1 je 4.14470000000000029E-003 0.0000000000000000 bc 8 quit </error> Regards Shriram

1 10

Lagrangian Particle Tracking
by nek5000-users＠lists.mcs.anl.gov 22 Mar '12

22 Mar '12

Hi, I have couple of questions regarding the subroutine given for Lagrangian particle tracking at the following wiki page and I wonder if there are any hints. http://nek5000.mcs.anl.gov/index.php/Data_processing_example - Where should this subroutine be inserted? - What are the variables "rtl" and "mr" representing? - How to set the initial positions of the particles? - Should the constant parameters "npart" and "dt_p" be inserted in SIZE file? Thanks in advance, Regards, Alireza Karimi

1 3

nek tools bug ?
by nek5000-users＠lists.mcs.anl.gov 30 Sep '11

30 Sep '11

hi all, I had problems merging grids with the nek tools. When the final merged grid was having more than 100000 elements genmap was failing. I think the problem was due to an inconsistency between the writing formats for the BC section in n2to3,nekmerge and genmap (at least in my version, see below). In my case i have solved the problem by changing the formats in nekmerge.f. just my 2 cents francesco ------------------------------ n2to3.f lines 441-... ------------------------------ if (neln.lt.1000) then write(11,20) cbc(k,e),id,k,(bc(j,k,e),j=1,5) elseif (neln.lt.100000) then write(11,21) cbc(k,e),id,k,(bc(j,k,e),j=1,5) else write(11,22) cbc(k,e),id,k,(bc(j,k,e),j=1,5) endif .... 20 FORMAT(1x,A3,2I3,5G14.7) 21 FORMAT(1x,A3,i5,i1,5G14.7) 22 FORMAT(1x,A3,i6,i1,5G14.7) ---------------------------------- nekmerge.f lines 664-... ---------------------------------- if (nel.lt.1000) then write(11,20) cbc(f,e,fld),e,f,(bc(j,f,e,fld),j=1,5) elseif (nel.lt.100000) then write(11,21) cbc(f,e,fld),e,f,(bc(j,f,e,fld),j=1,5) else write(11,22) cbc(f,e,fld),e,(bc(j,f,e,fld),j=1,5) endif ... 20 format(1x,a3,2i3,5g14.6) 21 format(1x,a3,i5,i1,5g14.6) 22 format(1x,a3,i10,5g14.6) -------------------------------- genmap.f lines 507-... -------------------------------- if (nel.lt.1000) then read(io,50,err=510,end=600) $ chtemp, $ cbc(f,e),id1,id2, $ (bc(ii,f,e),ii=1,nbcrea) 50 format(a1,a3,2i3,5g14.7) elseif (nel.lt.100 000) then read(io,51,err=520,end=600) $ chtemp, $ cbc(f,e),id1,id2, $ (bc(ii,f,e),ii=1,nbcrea) 51 format(a1,a3,i5,i1,5g14.7) elseif (nel.lt.1 000 000) then read(io,52,err=530,end=600) $ cbc(f,e),id1,(bc(ii,f,e),ii=1,nbcrea) 52 format(1x,a3,i7,5g14.6) elseif (nel.lt.10 000 000) then read(io,53,err=540,end=600) $ cbc(f,e),id1,(bc(ii,f,e),ii=1,nbcrea) 53 format(1x,a3,i7,5e20.12) else read(io,*,err=550,end=600) $ cbc(f,e),id1,(bc(ii,f,e),ii=1,nbcrea) endif

1 0

Re: [Nek5000-users] Running Simulations in Parallel with nekbmpi
by nek5000-users＠lists.mcs.anl.gov 27 Sep '11

27 Sep '11

Hi, I have lp=512 in SIZE, and I did compile with IFMPI=true. I recompiled yesterday because I thought this may be the case, but I was still getting the same error. When I first start a simulation and run top, it shows four processors running nek5000, but this changes to one processor after a few seconds. Thanks, John

1 2

CUBIT/Moab
by nek5000-users＠lists.mcs.anl.gov 21 Sep '11

21 Sep '11

Neks, We are interested in using CUBIT meshes in our Nek runs. I took a look at the "moab" example in the source, but I'm still unsure of the proper process of getting a CUBIT-generated mesh into Nek. I have created a CUBIT mesh (.cub), and have also downloaded, compiled, and installed MOAB. I see in the example that a .h5m file is used (the native format of MOAB), but I don't know how to generate it through MOAB (from my .cub file). I realize that this is probably more of a MOAB question than a Nek question, but are there some simple instructions/examples to help me in understanding how to use a CUBIT mesh in Nek? Any help would be appreciated! Thanks, -- Josh Camp "All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke

1 3

Problems with grid. Nek gets stuck.
by nek5000-users＠lists.mcs.anl.gov 20 Sep '11

20 Sep '11

Hi Paul, I'm having trouble with a run involving a box with dimensions 5:5:1 (length:breadth:height), for which I would like to go to a Rayleigh number of 10^8. I generated the grid (.rea and .map) using the attached .box file, but the run gets stuck as shown in the log. I think that there is a problem with the grid and the .map file. Could you please look into it? Thanks, Mani

1 0