[Nek5000-users] Problem in 2d pipe simulation
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jan 13 13:03:54 CST 2014
Thanks a lot.
Could you please tell me where I can find more information of defining the .usr file.
for example if I have boundary condition like this on a pipe
" v,O,W,W "
and If I want to define a parabolic profile on the inlet say a value of 10 . How do I specific the edge ?
Thank you.
Kamal
On Jan 13, 2014, at 7:58 PM, nek5000-users at lists.mcs.anl.gov wrote:
> 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
>
>
> _______________________________________________
> 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/055bb78e/attachment-0001.html>
More information about the Nek5000-users
mailing list