[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:30:29 CST 2014
Hi Kamal,
Go to
https://nek5000.mcs.anl.gov/index.php/UG#Documentation
and click on the Nek5000 manual. On pg. 35 there are
examples of how to specify functions.
Several of the cases in /examples also have prescribed functions
for boundary conditions --- for example, the Kovasznay case.
Best,
Paul
On Mon, 13 Jan 2014, nek5000-users at lists.mcs.anl.gov wrote:
> 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
>
>
More information about the Nek5000-users
mailing list