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

Re: [Nek5000-users] History points
by nek5000-users＠lists.mcs.anl.gov 06 May '10

06 May '10

Hi Jean-Christophe, This examples forces Nek5000 to output into the file blah.sch a time trace (H) of pressure (P), x & y velocity components (UV) for a point located in the element w/ global number IEL and collocation coordinates I,J K. Use 'grep History logfile' to get coordinates of the points. Aleks ----- Original Message ----- From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov Sent: Wednesday, May 5, 2010 12:22:03 PM GMT -06:00 US/Canada Central Subject: [Nek5000-users] History points Hello Nek's, I have a few questions on how to use the history points. I have seen on the wiki the following stuff: History points Specify grid points to monitor the specified variables in time. Example: ***** HISTORY AND INTEGRAL DATA ***** 1 POINTS. Var, Hcode, I, J, K, IEL UV P H 2 2 1 4 NOTE: • Make sure LHIS in SIZEu is large enough I bet Var stands for variable, but what do Hcode, I, J, K, and IEL stand for ? Once you have specified the location of your history point, let's say a point on the centerline of a stenotic flow at whatever streamwise station, how do you get the data concerning this history point? Do you have to change anything else in the output specification section ? I have tried to add the wiki lines up there directly in my blah.rea but once I ran the code I did not get any new files other than the expected blah.fld's. Regards, Jean-Christophe _______________________________________________ Nek5000-users mailing list Nek5000-users(a)lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

1 5

History points
by nek5000-users＠lists.mcs.anl.gov 05 May '10

05 May '10

Hello Nek's, I have a few questions on how to use the history points. I have seen on the wiki the following stuff: > > > History points > > Specify grid points to monitor the specified variables in time. Example: > > ***** HISTORY AND INTEGRAL DATA ***** > 1 POINTS. Var, Hcode, I, J, K, IEL > UV P H 2 2 1 4 > > > NOTE: > > * Make sure LHIS in SIZEu is large enough > I bet Var stands for variable, but what do Hcode, I, J, K, and IEL stand for ? Once you have specified the location of your history point, let's say a point on the centerline of a stenotic flow at whatever streamwise station, how do you get the data concerning this history point? Do you have to change anything else in the output specification section ? I have tried to add the wiki lines up there directly in my blah.rea but once I ran the code I did not get any new files other than the expected blah.fld's. Regards, Jean-Christophe

1 0

Generating curved sides with genbox?
by nek5000-users＠lists.mcs.anl.gov 04 May '10

04 May '10

Hello, Is it possible to generate curved sides with genbox? Or do I have to use prex for generating it? Is there a small example of generating a simple cylinder mesh? Thanks, Harish.

1 0

variable properties, uservp
by nek5000-users＠lists.mcs.anl.gov 04 May '10

04 May '10

Hello all, I am trying to understand how to set variable fluid properties depending on spatial location. I have two cases, for one "uservp" is never called, for the other it is. The difference in their *.rea files is as follows. I cannot see why in one case "uservp" is called as neither IFUSERVP=.true. or param(3)>0. Cheers, Frank [fmuldoo@tetra sim3]# diff *rea ../sim7/AR-.66-RE-1800-PR-04.0-dt-3.1e-5-time-006.0-grid-1-Intel/delet.rea | more 3c3 < 2 DIMENSIONAL RUN --- > 3 DIMENSIONAL RUN 6c6 < -1.00000 VISCOS --- > 1.00000 VISCOS 12c12 < 1.00000 CONDUCT --- > 0.25 CONDUCT 15,16c15,16 < 1000.0 NSTEPS < -5.e-6 DT --- > 193548 NSTEPS > -3.1e-05 DT 19c19 < 20.00000 IOSTEP --- > 1290 IOSTEP 25,26c25,26 < 1.000000E-07 DIVERGENCE < 1.000000E-07 HELMHOLTZ --- > 1.000000E-05 DIVERGENCE > 1.000000E-05 HELMHOLTZ 28,29c28,29 < 1.000000E-07 TOLREL < 1.000000E-07 TOLABS --- > 1.000000E-05 TOLREL > 1.000000E-05 TOLABS 39a40,45 > 0.000000E+00 XMAGNET > 0.000000E+00 NGRIDS > 0.000000E+00 NORDER2 > 0.000000E+00 NORDER3 > 0.000000E+00 NORDER4 > 0.000000E+00 NORDER5 46,51d51 < 1.000000E+00 p42=1 for std 2-level schwarz < 1.000000E+00 p43=1 etc. < 1.000000E+00 p44=1 < 0.000000E+00 < 0.000000E+00 < 0.3 Poisson ratio for mesh 107c107 < 0.05000 p103: filter strength, < 0 ==> no fil --- > 0.00000 p103: filter strength, < 0 ==> no fil 119c119 < T IFAXIS --- > F IFAXIS 125,126c125,126 < T IFMVBD < F IFCHAR --- > F IFMVBD > T IFCHAR -- Frank Herbert Muldoon, Ph.D. Mechanical Engineering Technische Universität Wien (Technical University of Vienna) Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid Mechanics and Heat Transfer) Resselgasse 3 1040 Wien Tel: +4315880132232 Fax: +4315880132299 Cell:+436765203470 fmuldoo (skype) http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.h…

1 1

Boundary condition "I"
by nek5000-users＠lists.mcs.anl.gov 04 May '10

04 May '10

Hello all, Could anyone tell me what the boundary condition "I" means? Cheers, Frank -- Frank Herbert Muldoon, Ph.D. Mechanical Engineering Technische Universität Wien (Technical University of Vienna) Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid Mechanics and Heat Transfer) Resselgasse 3 1040 Wien Tel: +4315880132232 Fax: +4315880132299 Cell:+436765203470 fmuldoo (skype) http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.h…

1 8

genbox error
by nek5000-users＠lists.mcs.anl.gov 03 May '10

03 May '10

Hello all, I have come across what may be an error in genbox. It concerns reading in data. If instead of 8 elements, 7 are used or if the trailing zeros are removed, all is fine. It almost seems as if genbox has a limitation on the length of the lines defining the element edges. Cheers, Frank Reading 32 = 8 4 4 elements for box 1. Reading 2 = error: 8 9 error: 0 5 The following is the input file: openboat2D.rea -2 spatial dimension 2 number of fields #======================================================== # # Open boat with the Marangoni effect at the top surface. # #======================================================== # Box 8 4 -4.000000 -3.946236 -3.806717 -3.497438 -3.000000 -2.500000 -2.000000 -1.500000 -1.000000 -5.000000 -3.347128 0.000000 3.754690 5.000000 W ,W ,W ,shl, V bc's ! NB: 3 characters each ! t ,t , , , T bc's ! You must have 2 spaces!! -- Frank Herbert Muldoon, Ph.D. Mechanical Engineering Technische Universität Wien (Technical University of Vienna) Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid Mechanics and Heat Transfer) Resselgasse 3 1040 Wien Tel: +4315880132232 Fax: +4315880132299 Cell:+436765203470 fmuldoo (skype) http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.h…

1 6

User feedback on Wiki
by nek5000-users＠lists.mcs.anl.gov 03 May '10

03 May '10

Howdy Nek users, bundling user feedback from Dr. Duggleby's research group here at A&M, I started a "Not-only-for-Christmas-Nek-wishlist" on the wiki: https://nek5000.mcs.anl.gov/index.php/USER_feedback Please add to it any items you would like to see implemented/improved in the nek5000 package or other feedback. Markus

1 1

Calculating forces on sphere
by nek5000-users＠lists.mcs.anl.gov 03 May '10

03 May '10

Dear Neks, for my case 'periodic channel flow around a sphere' I wanted to calculate the forces acting on the sphere. I've seen that there exists a 'drag_calc' subroutine in navier5.f. But as I understood I have to define my sphere as an object, that this works. Is this right?? Because until now the surface of the sphere is defined as 'wall'. And if I have to define it as an object, I wanted to ask if there is an easy and fast doing this? Is there maybe a way doing this similar to defining the boundary conditions with 'set entire level', or do I have to add all the faces of the sphere by hand in prenek?? Thanks a lot for you time. Best, Fred

1 1

Restarting simulation from instantenous & average velocity field files
by nek5000-users＠lists.mcs.anl.gov 02 May '10

02 May '10

Hello Paul & Stefan, I am running a DNS simulation of backward facing simulation (BFS) at Re=3000. The Reynolds # is based on step height (h=1) & free stream velocity (ux = 1.0). The BFS geometry has 5184 elements and lx1=ly1=lz1=8. After running for about 10 flow throughs, we started computing averages. I want to restart simulation after 2 flow throughs, but I have 2 field files one for instantaneous (blah.fld01) and the other for average (avgblah.fld01). To continue computing averages I need to read in both the velocity field files. Is there an option to restart from 2 field files? Or how else can I do this? I will also require this option later for computing "rms" velocities. Thanks, Harish.

1 5

Periodic BC
by nek5000-users＠lists.mcs.anl.gov 30 Apr '10

30 Apr '10

Hello all, I have a real simple mesh, one element. Running genmap with periodic BC's fails. Replace the periodic BC's with walls and then genmap works. Anybody have any ideas? Cheers, Frank ELEMENT 1 [ 1a] GROUP 0 -0.400000E+01 0.400000E+01 0.400000E+01 -0.400000E+01 -0.500000E+00 -0.500000E+00 0.500000E+00 0.500000E+00 -0.500000E+00 -0.500000E+00 -0.500000E+00 -0.500000E+00 -0.400000E+01 0.400000E+01 0.400000E+01 -0.400000E+01 -0.500000E+00 -0.500000E+00 0.500000E+00 0.500000E+00 0.500000E+00 0.500000E+00 0.500000E+00 0.500000E+00 ***** CURVED SIDE DATA ***** 0 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE ***** BOUNDARY CONDITIONS ***** ***** FLUID BOUNDARY CONDITIONS ***** W 1 1 0.000000 0.000000 0.000000 0.000000 0.000000 W 1 2 0.000000 0.000000 0.000000 0.000000 0.000000 shl 1 3 0.000000 0.000000 0.000000 0.000000 0.000000 W 1 4 0.000000 0.000000 0.000000 0.000000 0.000000 P 1 5 1.000000 6.000000 0.000000 0.000000 0.000000 P 1 6 1.000000 5.000000 0.000000 0.000000 0.000000 ***** THERMAL BOUNDARY CONDITIONS ***** 1 1 0.000000 0.000000 0.000000 0.000000 0.000000 t 1 2 0.000000 0.000000 0.000000 0.000000 0.000000 1 3 0.000000 0.000000 0.000000 0.000000 0.000000 t 1 4 0.000000 0.000000 0.000000 0.000000 0.000000 P 1 5 1.000000 6.000000 0.000000 0.000000 0.000000 P 1 6 1.000000 5.000000 0.000000 0.000000 0.000000 -- Frank Herbert Muldoon, Ph.D. Mechanical Engineering Technische Universität Wien (Technical University of Vienna) Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid Mechanics and Heat Transfer) Resselgasse 3 1040 Wien Tel: +4315880132232 Fax: +4315880132299 Cell:+436765203470 fmuldoo (skype) http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.h…

1 1