August 2014 - 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 when creating own mesh
by nek5000-users＠lists.mcs.anl.gov 01 Mar '16

01 Mar '16

Dear Neks, I am have attempted to use nek5000, firstly by following some of the examples etc. Now I am trying to create my own simulation, which is based upon the ocyl examples. My efforts simply include creating my own grid points in the .rea file to create a slightly different dimension... I believe I have understood how to define the problem however I continually receive errors when creating the genmap it reads as follows.... Input (.rea) file name: ocyl Input mesh tolerance (default 0.2): NOTE: smaller is better, but generous is more forgiving for bad meshes. 0.2 reading .rea file data ... ERROR: error reading 1 11 24 aborting 510 in routine rdbdry. 1 quit All I have done again is change the dimensions of the elements in the .rea file so I do not see why it would now not work? Thank you for any help in advance, Friedrich.

1 7

Magnetoconvection setup
by nek5000-users＠lists.mcs.anl.gov 09 Sep '14

09 Sep '14

Dear NEK developers, We are trying to set up a magneto-convection case with NEK. Can you kindly give us pointers on how to do that? Also we would like to know the form of equations that are implemented in case of magneto-convection. with best wishes and regards, Supriyo Paul *********************************************************** *Dr. Supriyo Paul* CFD/CAE Centre for Development of Advanced Computing (C-DAC) Pune University Campus, Ganeshkhind, Pune - 411007 INDIA Phone:+91-202-570-4175 Mobile:+91-860-596-5538

1 1

Outpost question
by nek5000-users＠lists.mcs.anl.gov 04 Sep '14

04 Sep '14

Howdy Nek's, I have a couple questions regarding the ouput of files: 1. How can force, within the code (say in usrchk), some files to be output in 64 bits for high-accuracy checkpointing and other in regular 8 bits for simple visualisation? 2. Say, during the first 1000 iterations, the code outpost a file every single iteration so that I have blah0.f00001 to blah0.f01000. How can I then force the code to outpost my new files (from istep = 1001 to 2000 say) rewriting on blah0.f00001 to blah0.f01000? Thanks a lot. Regards, JC -- Jean-Christophe Loiseau Homepage <https://sites.google.com/site/loiseaujc/>

1 5

Constant mass flow rate in simulation
by nek5000-users＠lists.mcs.anl.gov 31 Aug '14

31 Aug '14

Dear Neks, Recently I've been struggling with constant mass flow rate in my turbulent channel flow simulation. Based on the userchk subroutine in 'turbChannel' tutorial case, I also used torque_calc(1.0,x0,.false.,.false.) to calculate the wall shear and then got the mean wall shear. After that I set mean pressure gradient 'ffx' equal to the mean wall shear in order to keep constant mass flow rate. In this case, the flow rate started from 1 (nondimensionalised case, bulk mean equal to 1) and kept increasing. Although the increment was very small at each time step, it never stopped and the bulk mean velocity has increased to 1.002. Is it appropriate to set constant mass flow rate by updating mean pressure gradient based on mean wall shear in Nek5000 or am I using the wrong calculation? Since the bulk mean is still increasing, the simulation is like acceleration channel flow. Does anyone has any clue about how to set constant mass flow rate for turbulent channel flow in Nek5000? Any idea is appreciated and many thanks in advance. Thanks and regards, Tony

1 1

Linking other libs to Nek, e.g. ESSL in Blue Gene
by nek5000-users＠lists.mcs.anl.gov 29 Aug '14

29 Aug '14

Dear Neks, I have a short technical question. For our post processing analysis of turbulent convection we want to use routines from ESSL Library. Where is the best place to include the link command to this library into the makefile or makenek? Thanks in advance, Joerg. ==================================================== Joerg Schumacher Professor for Fluid Mechanics Institute of Thermodynamics and Fluid Mechanics Technische Universität Ilmenau P.O. Box 100565 D-98684 Ilmenau Germany E-mail: joerg.schumacher(a)tu-ilmenau.de<mailto:[email protected]> http://www.tu-ilmenau.de/tsm Phone: +49-3677-69-2428 Fax: +49-3677-69-2411 ====================================================

1 3

boyd filtering in dyn. Smagorinsky model
by nek5000-users＠lists.mcs.anl.gov 28 Aug '14

28 Aug '14

Hello Neks, In *turbChannel* example, explicit modal filtering has been used in the dynamaic Smagorinsky case in the eddy_visc(eddif,e) routines. This routine calls the setup filtering routine *set_ds_filt(fh,fht,nt,diag,nx)*, where coefficients for diag() are defined manually. Cases for nt = nx1-1 and nt = nx1-2 are mentioned in the routines. I am wondering if anybody could tell me how to set up the routines for nt = nx1-3 and nt = nx1-4, where nt is the filtered grids. Thanks in advance, Sincerely, Tanmoy

1 0

Re: [Nek5000-users] Lagrangian Particle Tracking
by nek5000-users＠lists.mcs.anl.gov 28 Aug '14

28 Aug '14

Hi Arash, I guess you are looking for an LPT module to track inertial particles with wall-collision and so on. At the moment we are working on a more unified implementation of such routine to be added to the nek trunk. A version of that will be finalised soon. Cheers Azad 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. Lagrangian Particle Tracking (nek5000-users(a)lists.mcs.anl.gov) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 27 Aug 2014 15:10:28 +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] Lagrangian Particle Tracking > Message-ID: > <mailman.4524.1409152234.3424.nek5000-users(a)lists.mcs.anl.gov> > Content-Type: text/plain; charset="us-ascii" > > Hi, > > I searched about Lagrangian Particle Tracking subject in history > texts and I found that "hemi" example relates to it. > After a review of this example, some questions are: > - How can I have a constant number of particles in the domain? It > needs that the number of particles which come to the domain should > be equal to the number of particles exit the domain. > - How about the particles which contact to the wall? How do they come back? > - How can I attach some information (for example Temperature) to any > particle and change it each time step according to the flow > properties and then have ensemble average of the magnitude of these > information particles? > > Thank you > Arash > > > ------------------------------ > > _______________________________________________ > 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 66, Issue 11 > *********************************************

1 1

About governing equations
by nek5000-users＠lists.mcs.anl.gov 27 Aug '14

27 Aug '14

Hi, I read in website that NEK5000 can't solve fully compressible flows. How it changes density in low mach number flows (in which we have thermal dilatation effects)? Thanks

1 1

emergency exit because of number of processors
by nek5000-users＠lists.mcs.anl.gov 20 Aug '14

20 Aug '14

Hello Neks, I am running pipe flow with nek5000, however, once I increase the number of processors, after few iterations CFL dramatically goes up immediately and it breaks out my simulation (Emergency Exit). Does any one has any idea what is wrong with that?! Thanks, Ami

1 0