[Nek5000-users] Problem in 2d pipe simulation
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jan 13 12:58:18 CST 2014
Hi,
When you create the .box file, it is important that the line containing the
boundary conditions follows the format as described on the nek wiki and
documentation.
It is a 3 character code that defines each boundary condition, and the
whitespace must be maintained. For example,
"v ,A ,W ,W , , "
would be a valid entry, whereas using less than 3 characters will result in
commas and whitespace being entered into your .rea file.
More info and examples can be found:
http://nek5000.mcs.anl.gov/index.php/Genbox
Best,
Katie
On Mon, Jan 13, 2014 at 8:49 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
> Hi,
>
> I created my .box file, with help of genbox I created .rea file.
> later I created a .usr and SIZE file.
>
> I have attached all the files I used with this mail
>
> I launched the simulation and I got this error
>
> *****************************************************************
> /----------------------------------------------------------\\
> | _ __ ______ __ __ ______ ____ ____ ____ |
> | / | / // ____// //_/ / ____/ / __ \\ / __ \\ / __ \\ |
> | / |/ // __/ / ,< /___ \\ / / / // / / // / / / |
> | / /| // /___ / /| | ____/ / / /_/ // /_/ // /_/ / |
> | /_/ |_//_____//_/ |_|/_____/ \\____/ \\____/ \\____/ |
> | |
> |----------------------------------------------------------|
> | |
> | NEK5000: Open Source Spectral Element Solver |
> | COPYRIGHT (c) 2008-2010 UCHICAGO ARGONNE, LLC |
> | Version: 1.0rc1 / SVN r982 |
> | Web: http://nek5000.mcs.anl.gov |
> | |
> \\----------------------------------------------------------/
>
>
> Number of processors: 1
> REAL wdsize : 8
> INTEGER wdsize : 4
>
>
> Beginning session:
> /home/kamal/neksamples/mine/box.rea
>
>
> timer accuracy: 0.0000000E+00 sec
>
> read .rea file
> nelgt/nelgv/lelt: 15 15 15
> lx1 /lx2 /lx3 : 8 6 8
>
> mapping elements to processors
> 0 15 15 15 15 NELV
> RANK 0 IEG 1 2 3 4 5 6 7 8
> 9 10 11 12 13 14 15
> element load imbalance: 0 15 15
> done :: mapping elements to processors
>
> 0 objects found
> done :: read .rea file 0.67091E-03 sec
>
> Reset the target Courant number to .5
> setup mesh topology
> Right-handed check complete for 15 elements. OK.
> setvert2d: 8 252 792 252 252
> call usrsetvert
> done :: usrsetvert
>
> gs_setup: 0 unique labels shared
> handle bytes (avg, min, max): 3972 3972 3972
> buffer bytes (avg, min, max): 0 0 0
> setupds time 1.9503E-04 seconds 0 8 252 15
> 4 max multiplicity
> done :: setup mesh topology
>
> call usrdat
> done :: usrdat
>
> generate geometry data
> vol_t,vol_v: 1.9999999999999920 1.9999999999999920
> done :: generate geometry data
>
> call usrdat2
> done :: usrdat2
>
> regenerate geometry data 1
> vol_t,vol_v: 1.9999999999999938 1.9999999999999938
> done :: regenerate geometry data 1
>
> verify mesh topology
> 0.0000000000000000 4.0000000000000000 Xrange
> 0.0000000000000000 0.50000000000000000 Yrange
> 0.0000000000000000 0.0000000000000000 Zrange
> done :: verify mesh topology
>
> 103 Parameters from file:/home/kamal/neksamples/mine/box.rea
> 1 1.00000 p1 DENSITY
> 2 -20.0000 p2 VISCOS
> 7 1.00000 p7 RHOCP
> 8 1.00000 p8 CONDUCT
> 11 1000.0 p11 NSTEPS
> 12 -1.000000E-04 p12 DT
> 15 100.000 p15 IOSTEP
> 18 0.250000E-01 p18 GRID
> 19 -1.00000 p19 INTYPE
> 20 12.0000 p20 NORDER
> 21 0.100000E-10 p21 DIVERGENCE
> 22 0.100000E-11 p22 HELMHOLTZ
> 24 0.100000E-09 p24 TOLREL
> 25 0.100000E-09 p25 TOLABS
> 26 1.00000 p26 COURANT/NTAU
> 27 3.00000 p27 TORDER
> 28 0. p28 TORDER: mesh velocity (0: p28=p27)
> 43 1.00000 p43 0=semg/1=schwarz
> 66 4.00000 p66 output : <0=ascii, else binary
> 67 4.00000 p67 restart: <0=ascii, else binary
> 93 20.0000 p93 Number of previous pressure solns saved
> 94 5.00000 p94 start projecting velocity after p94 step
> 95 5.00000 p95 start projecting pressure after p95 step
> 96 1.00000 p96 u0 translational velocity (in .usr)
> 97 0.30000 p97 v0 translational velocity (in .usr)
> 99 3.00000 p99 dealiasing: <0--> off/3--> old/4--> new
> 102 1.00000 p102 Dump out divergence at each time step
>
> IFTRAN = T
> IFFLOW = T
> IFHEAT = F
> IFSPLIT = F
> IFLOMACH = F
> IFUSERVP = F
> IFUSERMV = F
> IFSTRS = F
> IFCHAR = F
> IFCYCLIC = F
> IFAXIS = F
> IFMVBD = F
> IFMELT = F
> IFMODEL = F
> IFKEPS = F
> IFMOAB = F
> IFNEKNEK = F
> IFSYNC = T
>
> IFVCOR = T
> IFINTQ = F
> IFCWUZ = F
> IFSWALL = F
> IFGEOM = F
> IFSURT = F
> IFWCNO = F
>
> IFTMSH for field 1 = F
> IFADVC for field 1 = T
> IFNONL for field 1 = F
>
> Dealiasing enabled, lxd= 12
>
> Estimated eigenvalues
> EIGAA = 20.047633939712757
> EIGGA = 46987.971345031350
> EIGAE = 0.61685027506808487
> EIGAS = 3.8461538461538464E-003
> EIGGE = 46987.971345031350
> EIGGS = 2.0000000000000000
>
> verify mesh topology
> 0.0000000000000000 4.0000000000000000 Xrange
> 0.0000000000000000 0.50000000000000000 Yrange
> 0.0000000000000000 0.0000000000000000 Zrange
> done :: verify mesh topology
>
> E-solver strategy: 1 itr
> 0 15 BC FAIL
> EXIT: E INVALID BC FOUND in genfast 15
>
> call exitt: dying ...
>
> backtrace(): obtained 9 stack frames.
> ./nek5000(print_stack_+0x1a) [0x4fe90a]
> ./nek5000(exitt_+0x20f) [0x5bc0af]
> ./nek5000(exitti_+0x12a) [0x5be72a]
> ./nek5000(gen_fast_+0x581) [0x4fa701]
> ./nek5000(set_overlap_+0x125) [0x4eb465]
> ./nek5000(nek_init_+0x82d) [0x41affd]
> ./nek5000(main+0x24) [0x419e54]
> /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7fd9f1f33de5]
> ./nek5000() [0x419e99]
>
> total elapsed time : 3.74889E-03 sec
> total solver time incl. I/O : 0.00000E+00 sec
> time/timestep : 0.00000E+00 sec
> CPU seconds/timestep/gridpt : 0.00000E+00 sec
>
>
> real 0m0.030s
> user 0m0.020s
> sys 0m0.000s
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140113/2a0e92b4/attachment.html>
More information about the Nek5000-users
mailing list