January 2018 - Nek5000-users - mailman3.cels.anl.gov

Re: [Nek5000-users] How to change the x y z coordinates to be polar
by nek5000-users＠lists.mcs.anl.gov 18 Jan '18

18 Jan '18

Hi Philipp, Thanks for your reply. Currently I have the mean statistics. Can I change the averaged velocities and velocities square to be based on polar coordinates, then have the rms in radial and azimuthal? Or I need to change the coordinate firstly before the averaging. Kind regards, Jian

1 2

Advection for passive scalars in par file
by nek5000-users＠lists.mcs.anl.gov 18 Jan '18

18 Jan '18

Dear Neks, when including passive scalars in my simulation, I set "advection = yes" for velocity, temperature, and passive scalars. However, in PARDICT the key pardictkey(93) / 'SCALAR%%:ADVECTION' / does not exist so I have added it myself to the list. Is that correct? I am using the official release (tarball version 17.0) -- Karlsruhe Institute of Technology (KIT) Institute of Fluid Mechanics M.Sc. Steffen Straub Doctoral Researcher Kaiserstraße 10 Building 10.23 76131 Karlsruhe, Germany Phone: +49 721 608-43027 E-mail: steffen.straub(a)kit.edu Web: http://www.istm.kit.edu Registered office: Kaiserstraße 12, 76131 Karlsruhe, Germany KIT – The Research University in the Helmholtz Association

1 1

How to change the x y z coordinates to be polar
by nek5000-users＠lists.mcs.anl.gov 17 Jan '18

17 Jan '18

Hi Nek, Do you know how to change the result of rms to be based on polar coordinates? Currently I have x, y and z directions for the circular pipe. I have to plot the rms in axial,radial and azimuthal. Cheers, Jian

1 2

Re: [Nek5000-users] Nek5000-users Digest, Vol 107, Issue 24
by nek5000-users＠lists.mcs.anl.gov 16 Jan '18

16 Jan '18

Hey Paul, As suggested by you, I did look into the .usr of Rayleigh Benard Convection example and it did help me. But I face the following error while compiling using makenek. obj/plan4.o: In function `plan4_': plan4.f:(.text+0x1781): additional relocation overflows omitted from the output collect2: error: ld returned 1 exit status makefile:114: recipe for target 'nek5000' failed make: *** [nek5000] Error 1 I have attached The SIZE file, .usr file,.par file and the .box file used for the simulation. Really appreciate your help. On Mon, Jan 15, 2018 at 11:30 PM, <nek5000-users-request(a)lists.mcs.anl.gov> wrote: > 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. Re: Low-Density jet using low -Mach number formulation > (nek5000-users(a)lists.mcs.anl.gov) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 14 Jan 2018 21:00:59 +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: Re: [Nek5000-users] Low-Density jet using low -Mach number > formulation > Message-ID: > <mailman.6277.1515963679.19102.nek5000-users(a)lists.mcs.anl.gov> > Content-Type: text/plain; charset="us-ascii" > > > Dear Nitish, > > > There are several examples of Rayleigh-Benard convection that show how to > do buoyancy driven flows. Looking at those usr files would be a first > step. > > > hth, > > > Paul > > > ________________________________ > From: Nek5000-users <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf > of nek5000-users(a)lists.mcs.anl.gov <nek5000-users(a)lists.mcs.anl.gov> > Sent: Sunday, January 14, 2018 8:08:20 AM > To: nek5000-users(a)lists.mcs.anl.gov > Subject: [Nek5000-users] Low-Density jet using low -Mach number formulation > > Hello Nek Users, > > I want to simulate low density jet (heated air whose density is lesser > than ambient air) expanding out from a pipe into stagnant air. > I want to run a 2D planar case. Im planning to export the mesh from exodus > II to nek. Can somebody guide me what I have to do from here? What Stress > formulation to use? Thanks for the help. > > -- > Nitish Kovalam > Aerospace Engineering > Indian Institute of Space Science and Technology >

1 1

conjugate heat transfer; NELT NELV; genmap
by nek5000-users＠lists.mcs.anl.gov 16 Jan '18

16 Jan '18

Hello Neks, I am using conjugate heat transfer module. Elements in velocity domain = 3072 = NELV Elements in temperature domain = 3072 Elements in PS1 domain = 7424 = NEL Elements in PS2 domain = 7424 = NEL When I generate .map file, genmap throws the following error " reading .rea file data ... ERROR: error reading 2 3073 7424 aborting 520 in routine rdbdry. " When I modify Elements in temperature domain = 7424 = NEL It runs fine. I am able to generate .map and run a simulation and obtained expected results. When the number of elements in atleast one scalar is not equal to NEL, there is a problem. Is it possible to have the number of elements in temperature field = NELV and the number of elements in other scalars = NEL? Thank you advance. Cheers, Sandeep

1 0

Low-Density jet using low -Mach number formulation
by nek5000-users＠lists.mcs.anl.gov 14 Jan '18

14 Jan '18

Hello Nek Users, I want to simulate low density jet (heated air whose density is lesser than ambient air) expanding out from a pipe into stagnant air. I want to run a 2D planar case. Im planning to export the mesh from exodus II to nek. Can somebody guide me what I have to do from here? What Stress formulation to use? Thanks for the help. -- Nitish Kovalam Aerospace Engineering Indian Institute of Space Science and Technology

1 1

Periodicity with curvature
by nek5000-users＠lists.mcs.anl.gov 13 Jan '18

13 Jan '18

Dear Nek Users, I am trying to get periodicity to work in the lastest version of Nek. I've made an extremely simple mesh within prenek (prex) in polar coordinates and cannot seem to run this case. When I run genmap, I get the error below. Did I do something wrong in setting up this mesh? >> genmap reading a.rea Input mesh tolerance (default 0.2): NOTE: smaller is better, but generous is more forgiving for bad meshes. 0.01 reading mesh data ... start locglob_lexico: 4 2 8 1.0000000000000000E-002 locglob: 1 1 8 locglob: 2 5 8 locglob: 1 6 8 locglob: 2 6 8 done locglob_lexico: 6 6 8 2 start periodic vtx: 2 6 1 4 4 Matrix: x0 1 x0 1.0000E+00 1.0000E+00 5.4347-323 4.9407-324 2 x0 -3.5355E-01 3.5355E-01 6.9423-310 6.9423-310 3 x0 3.5355E-01 -3.5355E-01 6.9529-310 6.9529-310 4 x0 6.9423-310 6.9529-310 2.3477-313 3.1632-317 2 4 4 Matrix: x1 1 x1 1.0000E+00 1.0000E+00 6.3240-322 0.0000E+00 2 x1 -3.5355E-01 3.5355E-01 0.0000E+00 6.9423-310 3 x1 -3.5355E-01 3.5355E-01 6.9423-310 6.9529-310 4 x1 6.9423-310 6.9423-310 6.9423-310 6.9529-310 1 4 4 Matrix: z0 1 z0 1.0000E+00 1.0000E+00 5.4347-323 4.9407-324 2 z0 -2.1213E+00 -1.4142E+00 6.9423-310 6.9423-310 3 z0 2.1213E+00 1.4142E+00 6.9529-310 6.9529-310 4 z0 6.9423-310 6.9529-310 2.3477-313 3.1632-317 2 4 4 Matrix: z1 1 z1 1.0000E+00 1.0000E+00 6.3240-322 0.0000E+00 2 z1 1.4142E+00 2.1213E+00 0.0000E+00 6.9423-310 3 z1 1.4142E+00 2.1213E+00 6.9423-310 6.9529-310 4 z1 6.9423-310 6.9423-310 6.9423-310 6.9529-310 1 3 0 1.00000000E-03 2.12132000E+00 abort: FACE MATCH FAIL 2 4 3 2.12132000E-03 1.00000031E+00 abort: FACE MATCH FAIL 0 quit The mesh within my .rea file is as follows: 20.0000 20.0000 -8.00000 -4.00000 XFAC,YFAC,XZERO,YZERO **MESH DATA** 1st line is X of corner 1,2,3,4. 2nd line is Y. 2 2 2 NEL,NDIM,NELV ELEMENT 1 [ 1A] GROUP 0 -2.121320 -1.414213 -0.8742276E-07 -0.1311341E-06 2.121320 1.414213 2.000000 3.000000 ELEMENT 2 [ 1B] GROUP 0 -0.1311341E-06 -0.8742276E-07 1.414213 2.121320 3.000000 2.000000 1.414213 2.121320 ***** CURVED SIDE DATA ***** 4 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE 2 1 -2.00000 0.00000 0.00000 0.00000 0.00000 C 4 1 3.00000 0.00000 0.00000 0.00000 0.00000 C 2 2 -2.00000 0.00000 0.00000 0.00000 0.00000 C 4 2 3.00000 0.00000 0.00000 0.00000 0.00000 C ***** BOUNDARY CONDITIONS ***** ***** FLUID BOUNDARY CONDITIONS ***** P 1 1 2.00000 3.00000 0.00000 0.00000 0.00000 W 1 2 0.00000 0.00000 0.00000 0.00000 0.00000 E 1 3 2.00000 1.00000 0.00000 0.00000 0.00000 W 1 4 0.00000 0.00000 0.00000 0.00000 0.00000 E 2 1 1.00000 3.00000 0.00000 0.00000 0.00000 W 2 2 0.00000 0.00000 0.00000 0.00000 0.00000 P 2 3 1.00000 1.00000 0.00000 0.00000 0.00000 W 2 4 0.00000 0.00000 0.00000 0.00000 0.00000 ***** NO THERMAL BOUNDARY CONDITIONS ***** [image: Inline image 2]

1 9

Error using [PROBLEMTYPE] in .par file
by nek5000-users＠lists.mcs.anl.gov 13 Jan '18

13 Jan '18

Greetings, I have been trying to execute the periodic hill example as described on the Nek Website. But there is a problem with the statement "equation = incompNS" in hillp.par under the [PROBLEMTYPE] category.It doesnt seem to recognise the statement. On the other hand while running the 'low Mach test' problem from the examples which has the statement "lowMachnumber = yes" in the .par file, Im able to get the solution up and running. Can somebody please tell what i'm missing? Nitish -- Nitish Kovalam Aerospace Engineering Indian Institute of Space Science and Technology

1 2

Error using [PROBLEMTYPE] in .par file
by nek5000-users＠lists.mcs.anl.gov 13 Jan '18

13 Jan '18

Sir, Yes. Im using the latest version of Nek. The problem still persists. Nitish -- Nitish Kovalam Aerospace Engineering Indian Institute of Space Science and Technology

1 0

A postprocessing problem and a query about using LES
by nek5000-users＠lists.mcs.anl.gov 12 Jan '18

12 Jan '18

Dear experts in Nek5000 : Thank you for reading my email. I am a freshman in using Nek5000 and I am willing to simulate unsteady flow phenomenon with LES method through parallel processes. And Nek5000 is the best solver for resolving my problem and I hope to reach my goal with the outstanding SEM-CFD tools. And I have been wondering 2 questions for a certain period of time: First of all, could you please introduce an appropriate sample for me to learn how to run Nek5000 within the LES formulation? And for the 2nd problem, when I dump the final result in .nek5000 format into the PARAVIEW postprocessing software, I could only see some principle fuild variaties fields to be selected in the GUI such as Vx Vy pressure and so on, and could you please give a proper example to show me how to add another paramters like fluctuating velocity in turbulent flow calculating by LES through the list in order to give more details, to illustrate my flow problem better. Thank you all the time for your supporting me and I am waiting for your earlier reply . Best regards Haoran JU

1 1