[Nek5000-users] Problem in 2d pipe simulation
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jan 13 08:49:56 CST 2014
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
-------------- next part --------------
15 24 3 8 60 24 0
7 8 7 19 18
7 7 16 18 17
6 16 24 17 23
3 24 6 23 11
3 6 5 11 10
5 19 18 15 14
5 18 17 14 22
6 17 23 22 21
2 23 11 21 13
2 11 10 13 12
4 15 14 4 3
4 14 22 3 20
1 22 21 20 2
1 21 13 2 9
0 13 12 9 1
-------------- next part --------------
****** PARAMETERS *****
2.610000 NEKTON VERSION
2 DIMENSIONAL RUN
103 PARAMETERS FOLLOW
1.00000 p1 DENSITY
-20.0000 p2 VISCOS
0.
0.
0.
0.
1.00000 p7 RHOCP
1.00000 p8 CONDUCT
0.
0. p10 FINTIME
1000.0 p11 NSTEPS
-1.000000E-04 p12 DT
0. p13 IOCOMM
0. p14 IOTIME
100.000 p15 IOSTEP
0. p16 PSSOLVER
0.
0.250000E-01 p18 GRID
-1.00000 p19 INTYPE
12.0000 p20 NORDER
0.100000E-10 p21 DIVERGENCE
0.100000E-11 p22 HELMHOLTZ
0 p23 NPSCAL
0.100000E-09 p24 TOLREL
0.100000E-09 p25 TOLABS
1.00000 p26 COURANT/NTAU
3.00000 p27 TORDER
0. p28 TORDER: mesh velocity (0: p28=p27)
0. p29 magnetic visc if > 0, = -1/Rm if < 0
0. p30 > 0 ==> properties set in uservp()
0. p31 NPERT: #perturbation modes
0. p32 #BCs in re2 file, if > 0
0.
0.
0.
0.
0.
0.
0.
0.
0. p41 1-->multiplicative SEMG
0. p42 0=gmres/1=pcg
1.00000 p43 0=semg/1=schwarz
0. p44 0=E-based/1=A-based prec.
0. p45 Relaxation factor for DTFS
0. p46 reserved
0. p47 vnu: mesh matieral prop
0.
0.
0.
0.
0. p52 IOHIS
0.
0. p54 1,2,3-->fixed flow rate dir=x,y,z
0. p55 vol.flow rate (p54>0) or Ubar (p54<0)
0.
0.
0.
0. p59 !=0 --> full Jac. eval. for each el.
0. p60 !=0 --> init. velocity to small nonzero
0.
0. p62 >0 --> force byte_swap for output
0. p63 =8 --> force 8-byte output
0. p64 =1 --> perturbation restart
0. p65 #iofiles (eg, 0 or 64); <0 --> sep. dirs
4.00000 p66 output : <0=ascii, else binary
4.00000 p67 restart: <0=ascii, else binary
0. p68 iastep: freq for avg_all
0.
0.
0.
0.
0.
0. p74 verbose Helmholtz
0.
0.
0.
0.
0.
0.
0.
0.
0.
0. p84 !=0 --> sets initial timestep if p12>0
0. p85 dt ratio if p84 !=0, for timesteps>0
0. p86 reserved
0.
0.
0.
0.
0.
0.
20.0000 p93 Number of previous pressure solns saved
5.00000 p94 start projecting velocity after p94 step
5.00000 p95 start projecting pressure after p95 step
1.00000 p96 u0 translational velocity (in .usr)
0.30000 p97 v0 translational velocity (in .usr)
0.
3.00000 p99 dealiasing: <0--> off/3--> old/4--> new
0.
0. p101 No. of additional filter modes
1.00000 p102 Dump out divergence at each time step
.000 p103 weight of stabilizing filter (.01)
4 Lines of passive scalar data follows2 CONDUCT; 2RHOCP
1.00000 1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000
13 LOGICAL SWITCHES FOLLOW
T IFFLOW
F IFHEAT
T IFTRAN
T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)
F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)
F IFAXIS
F IFSTRS
F IFSPLIT
F IFMGRID
F IFMODEL
F IFKEPS
F IFMVBD
F IFCHAR
4.00000 0.50000 0.00 0.00 XFAC,YFAC,XZERO,YZERO
**MESH DATA** 6 lines are X,Y,Z;X,Y,Z. Columns corners 1-4;5-8
15 2 15 NEL,NDIM,NELV
ELEMENT 1 [ 1a] GROUP 0
0.00000 0.800000 0.800000 0.00000
0.00000 0.00000 0.166667 0.166667
ELEMENT 2 [ 1b] GROUP 0
0.800000 1.60000 1.60000 0.800000
0.00000 0.00000 0.166667 0.166667
ELEMENT 3 [ 1c] GROUP 0
1.60000 2.40000 2.40000 1.60000
0.00000 0.00000 0.166667 0.166667
ELEMENT 4 [ 1d] GROUP 0
2.40000 3.20000 3.20000 2.40000
0.00000 0.00000 0.166667 0.166667
ELEMENT 5 [ 1e] GROUP 0
3.20000 4.00000 4.00000 3.20000
0.00000 0.00000 0.166667 0.166667
ELEMENT 6 [ 1f] GROUP 0
0.00000 0.800000 0.800000 0.00000
0.166667 0.166667 0.333333 0.333333
ELEMENT 7 [ 1g] GROUP 0
0.800000 1.60000 1.60000 0.800000
0.166667 0.166667 0.333333 0.333333
ELEMENT 8 [ 1h] GROUP 0
1.60000 2.40000 2.40000 1.60000
0.166667 0.166667 0.333333 0.333333
ELEMENT 9 [ 1i] GROUP 0
2.40000 3.20000 3.20000 2.40000
0.166667 0.166667 0.333333 0.333333
ELEMENT 10 [ 1j] GROUP 0
3.20000 4.00000 4.00000 3.20000
0.166667 0.166667 0.333333 0.333333
ELEMENT 11 [ 1k] GROUP 0
0.00000 0.800000 0.800000 0.00000
0.333333 0.333333 0.500000 0.500000
ELEMENT 12 [ 1l] GROUP 0
0.800000 1.60000 1.60000 0.800000
0.333333 0.333333 0.500000 0.500000
ELEMENT 13 [ 1m] GROUP 0
1.60000 2.40000 2.40000 1.60000
0.333333 0.333333 0.500000 0.500000
ELEMENT 14 [ 1n] GROUP 0
2.40000 3.20000 3.20000 2.40000
0.333333 0.333333 0.500000 0.500000
ELEMENT 15 [ 1o] GROUP 0
3.20000 4.00000 4.00000 3.20000
0.333333 0.333333 0.500000 0.500000
***** CURVED SIDE DATA *****
0 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE
***** BOUNDARY CONDITIONS *****
***** FLUID BOUNDARY CONDITIONS *****
,W 1 1 0.00000 0.00000 0.00000 0.00000 0.00000
E 1 2 2.00000 4.00000 0.00000 0.00000 0.00000
E 1 3 6.00000 1.00000 0.00000 0.00000 0.00000
v , 1 4 0.00000 0.00000 0.00000 0.00000 0.00000
,W 2 1 0.00000 0.00000 0.00000 0.00000 0.00000
E 2 2 3.00000 4.00000 0.00000 0.00000 0.00000
E 2 3 7.00000 1.00000 0.00000 0.00000 0.00000
E 2 4 1.00000 2.00000 0.00000 0.00000 0.00000
,W 3 1 0.00000 0.00000 0.00000 0.00000 0.00000
E 3 2 4.00000 4.00000 0.00000 0.00000 0.00000
E 3 3 8.00000 1.00000 0.00000 0.00000 0.00000
E 3 4 2.00000 2.00000 0.00000 0.00000 0.00000
,W 4 1 0.00000 0.00000 0.00000 0.00000 0.00000
E 4 2 5.00000 4.00000 0.00000 0.00000 0.00000
E 4 3 9.00000 1.00000 0.00000 0.00000 0.00000
E 4 4 3.00000 2.00000 0.00000 0.00000 0.00000
,W 5 1 0.00000 0.00000 0.00000 0.00000 0.00000
,A 5 2 0.00000 0.00000 0.00000 0.00000 0.00000
E 5 3 10.0000 1.00000 0.00000 0.00000 0.00000
E 5 4 4.00000 2.00000 0.00000 0.00000 0.00000
E 6 1 1.00000 3.00000 0.00000 0.00000 0.00000
E 6 2 7.00000 4.00000 0.00000 0.00000 0.00000
E 6 3 11.0000 1.00000 0.00000 0.00000 0.00000
v , 6 4 0.00000 0.00000 0.00000 0.00000 0.00000
E 7 1 2.00000 3.00000 0.00000 0.00000 0.00000
E 7 2 8.00000 4.00000 0.00000 0.00000 0.00000
E 7 3 12.0000 1.00000 0.00000 0.00000 0.00000
E 7 4 6.00000 2.00000 0.00000 0.00000 0.00000
E 8 1 3.00000 3.00000 0.00000 0.00000 0.00000
E 8 2 9.00000 4.00000 0.00000 0.00000 0.00000
E 8 3 13.0000 1.00000 0.00000 0.00000 0.00000
E 8 4 7.00000 2.00000 0.00000 0.00000 0.00000
E 9 1 4.00000 3.00000 0.00000 0.00000 0.00000
E 9 2 10.0000 4.00000 0.00000 0.00000 0.00000
E 9 3 14.0000 1.00000 0.00000 0.00000 0.00000
E 9 4 8.00000 2.00000 0.00000 0.00000 0.00000
E 10 1 5.00000 3.00000 0.00000 0.00000 0.00000
,A 10 2 0.00000 0.00000 0.00000 0.00000 0.00000
E 10 3 15.0000 1.00000 0.00000 0.00000 0.00000
E 10 4 9.00000 2.00000 0.00000 0.00000 0.00000
E 11 1 6.00000 3.00000 0.00000 0.00000 0.00000
E 11 2 12.0000 4.00000 0.00000 0.00000 0.00000
, 11 3 0.00000 0.00000 0.00000 0.00000 0.00000
v , 11 4 0.00000 0.00000 0.00000 0.00000 0.00000
E 12 1 7.00000 3.00000 0.00000 0.00000 0.00000
E 12 2 13.0000 4.00000 0.00000 0.00000 0.00000
, 12 3 0.00000 0.00000 0.00000 0.00000 0.00000
E 12 4 11.0000 2.00000 0.00000 0.00000 0.00000
E 13 1 8.00000 3.00000 0.00000 0.00000 0.00000
E 13 2 14.0000 4.00000 0.00000 0.00000 0.00000
, 13 3 0.00000 0.00000 0.00000 0.00000 0.00000
E 13 4 12.0000 2.00000 0.00000 0.00000 0.00000
E 14 1 9.00000 3.00000 0.00000 0.00000 0.00000
E 14 2 15.0000 4.00000 0.00000 0.00000 0.00000
, 14 3 0.00000 0.00000 0.00000 0.00000 0.00000
E 14 4 13.0000 2.00000 0.00000 0.00000 0.00000
E 15 1 10.0000 3.00000 0.00000 0.00000 0.00000
,A 15 2 0.00000 0.00000 0.00000 0.00000 0.00000
, 15 3 0.00000 0.00000 0.00000 0.00000 0.00000
E 15 4 14.0000 2.00000 0.00000 0.00000 0.00000
***** NO THERMAL BOUNDARY CONDITIONS *****
0 PRESOLVE/RESTART OPTIONS *****
7 INITIAL CONDITIONS *****
C Default
C Default
C Default
C Default
C Default
C Default
C Default
***** DRIVE FORCE DATA ***** BODY FORCE, FLOW, Q
4 Lines of Drive force data follow
C
C
C
C
***** Variable Property Data ***** Overrrides Parameter data.
1 Lines follow.
0 PACKETS OF DATA FOLLOW
***** HISTORY AND INTEGRAL DATA *****
0 POINTS. Hcode, I,J,H,IEL
***** OUTPUT FIELD SPECIFICATION *****
6 SPECIFICATIONS FOLLOW
T COORDINATES
T VELOCITY
T PRESSURE
F TEMPERATURE
F TEMPERATURE GRADIENT
0 PASSIVE SCALARS
***** OBJECT SPECIFICATION *****
0 Surface Objects
0 Volume Objects
0 Edge Objects
0 Point Objects
-------------- next part --------------
c-----------------------------------------------------------------------
subroutine uservp (ix,iy,iz,ieg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'
C
udiff =0.
utrans=0.
return
end
c-----------------------------------------------------------------------
subroutine userf (ix,iy,iz,ieg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'
c
ffx = 0.0
ffy = 0.0
ffz = 0.0
return
end
c-----------------------------------------------------------------------
subroutine userq (ix,iy,iz,ieg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'
C
qvol = 0.0
source = 0.0
return
end
c-----------------------------------------------------------------------
subroutine userchk ! called once per step
include 'SIZE'
include 'TOTAL'
c
n = nx1*ny1*nz1*nelv
visc = param(2)
u0 = param(96)
v0 = param(97)
if (istep.eq.0 ) call outpost(ue,ve,vx,pr,t,' ')
call sub3 (ud,ue,vx,n)
call sub3 (vd,ve,vy,n)
if (istep.eq.nsteps) call outpost(ud,vd,vx,pr,t,' ')
umx = glamax(vx,n)
vmx = glamax(vy,n)
uex = glamax(ue,n)
vex = glamax(ve,n)
udx = glamax(ud,n)
vdx = glamax(vd,n)
if (nid.eq.0) then
write(6,11) istep,time,udx,umx,uex,u0,' X err'
write(6,11) istep,time,vdx,vmx,vex,v0,' Y err'
11 format(i5,1p5e14.6,a7)
endif
if (istep.le.5) then ! Reset velocity to eliminate
call copy (vx,ue,n) ! start-up contributions to
call copy (vy,ve,n) ! temporal-accuracy behavior.
endif
return
end
c-----------------------------------------------------------------------
subroutine userbc (ix,iy,iz,iside,ieg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'
ux = 1.0
uy = 0.0
temp = 0
return
end
c-----------------------------------------------------------------------
subroutine useric (ix,iy,iz,ieg)
include 'SIZE'
include 'TOTAL'
include 'NEKUSE'
ux = 1.0
uy = 0.0
temp = 0
return
end
c-----------------------------------------------------------------------
subroutine usrdat
return
end
c-----------------------------------------------------------------------
subroutine usrdat2
return
end
c-----------------------------------------------------------------------
subroutine usrdat3
return
end
c-----------------------------------------------------------------------
-------------- next part --------------
C Dimension file to be included
C
C HCUBE array dimensions
C
parameter (ldim=2)
parameter (lx1=8,ly1=lx1,lz1=1,lelt=15,lelv=lelt)
parameter (lxd=12,lyd=lxd,lzd=1)
parameter (lelx=20,lely=20,lelz=1)
parameter (lzl=3 + 2*(ldim-3))
parameter (lx2=lx1-2)
parameter (ly2=ly1-2)
parameter (lz2=lz1 )
parameter (lx3=lx1)
parameter (ly3=ly1)
parameter (lz3=lz1)
parameter (lp = 512)
parameter (lelg = 4100)
c
c parameter (lpelv=lelv,lpelt=lelt,lpert=3) ! perturbation
c parameter (lpx1=lx1,lpy1=ly1,lpz1=lz1) ! array sizes
c parameter (lpx2=lx2,lpy2=ly2,lpz2=lz2)
c
parameter (lpelv=1,lpelt=1,lpert=1) ! perturbation
parameter (lpx1=1,lpy1=1,lpz1=1) ! array sizes
parameter (lpx2=1,lpy2=1,lpz2=1)
c
c parameter (lbelv=lelv,lbelt=lelt) ! MHD
c parameter (lbx1=lx1,lby1=ly1,lbz1=lz1) ! array sizes
c parameter (lbx2=lx2,lby2=ly2,lbz2=lz2)
c
parameter (lbelv=1,lbelt=1) ! MHD
parameter (lbx1=1,lby1=1,lbz1=1) ! array sizes
parameter (lbx2=1,lby2=1,lbz2=1)
C LX1M=LX1 when there are moving meshes; =1 otherwise
parameter (lx1m=1,ly1m=1,lz1m=1)
parameter (ldimt= 2) ! 2 passive scalars + T
parameter (ldimt1=ldimt+1)
parameter (ldimt3=ldimt+3)
c
c Note: In the new code, LELGEC should be about sqrt(LELG)
c
PARAMETER (LELGEC = 1)
PARAMETER (LXYZ2 = 1)
PARAMETER (LXZ21 = 1)
PARAMETER (LMAXV=LX1*LY1*LZ1*LELV)
PARAMETER (LMAXT=LX1*LY1*LZ1*LELT)
PARAMETER (LMAXP=LX2*LY2*LZ2*LELV)
PARAMETER (LXZ=LX1*LZ1)
PARAMETER (LORDER=3)
PARAMETER (MAXOBJ=4,MAXMBR=LELT*6)
PARAMETER (lhis=100) ! # of pts a proc reads from hpts.in
! Note: lhis*np > npoints in hpts.in
C
C Common Block Dimensions
C
PARAMETER (LCTMP0 =2*LX1*LY1*LZ1*LELT)
PARAMETER (LCTMP1 =4*LX1*LY1*LZ1*LELT)
C
C The parameter LVEC controls whether an additional 42 field arrays
C are required for Steady State Solutions. If you are not using
C Steady State, it is recommended that LVEC=1.
C
PARAMETER (LVEC=1)
C
C Uzawa projection array dimensions
C
parameter (mxprev = 20)
parameter (lgmres = 30)
C
C Split projection array dimensions
C
parameter(lmvec = 1)
parameter(lsvec = 1)
parameter(lstore=lmvec*lsvec)
c
c NONCONFORMING STUFF
c
parameter (maxmor = lelt)
C
C Array dimensions
C
COMMON/DIMN/NELV,NELT,NX1,NY1,NZ1,NX2,NY2,NZ2
$,NX3,NY3,NZ3,NDIM,NFIELD,NPERT,NID
$,NXD,NYD,NZD
c automatically added by makenek
parameter(lxo = lx1) ! max output grid size (lxo>=lx1)
c automatically added by makenek
parameter(lpart = 1 ) ! max number of particles
c automatically added by makenek
integer ax1,ay1,az1,ax2,ay2,az2
parameter (ax1=lx1,ay1=ly1,az1=lz1,ax2=lx2,ay2=ly2,az2=lz2) ! running averages
c automatically added by makenek
parameter (lxs=1,lys=lxs,lzs=(lxs-1)*(ldim-2)+1) !New Pressure Preconditioner
c automatically added by makenek
parameter (lfdm=0) ! == 1 for fast diagonalization method
-------------- next part --------------
****** PARAMETERS *****
2.610000 NEKTON VERSION
2 DIMENSIONAL RUN
103 PARAMETERS FOLLOW
1.00000 p1 DENSITY
-20.0000 p2 VISCOS
0.
0.
0.
0.
1.00000 p7 RHOCP
1.00000 p8 CONDUCT
0.
0. p10 FINTIME
1000.0 p11 NSTEPS
-1.000000E-04 p12 DT
0. p13 IOCOMM
0. p14 IOTIME
100.000 p15 IOSTEP
0. p16 PSSOLVER
0.
0.250000E-01 p18 GRID
-1.00000 p19 INTYPE
12.0000 p20 NORDER
0.100000E-10 p21 DIVERGENCE
0.100000E-11 p22 HELMHOLTZ
0. p23 NPSCAL
0.100000E-09 p24 TOLREL
0.100000E-09 p25 TOLABS
1.00000 p26 COURANT/NTAU
3.00000 p27 TORDER
0. p28 TORDER: mesh velocity (0: p28=p27)
0. p29 magnetic visc if > 0, = -1/Rm if < 0
0. p30 > 0 ==> properties set in uservp()
0. p31 NPERT: #perturbation modes
0. p32 #BCs in re2 file, if > 0
0.
0.
0.
0.
0.
0.
0.
0.
0. p41 1-->multiplicative SEMG
0. p42 0=gmres/1=pcg
1.00000 p43 0=semg/1=schwarz
0. p44 0=E-based/1=A-based prec.
0. p45 Relaxation factor for DTFS
0. p46 reserved
0. p47 vnu: mesh matieral prop
0.
0.
0.
0.
0. p52 IOHIS
0.
0. p54 1,2,3-->fixed flow rate dir=x,y,z
0. p55 vol.flow rate (p54>0) or Ubar (p54<0)
0.
0.
0.
0. p59 !=0 --> full Jac. eval. for each el.
0. p60 !=0 --> init. velocity to small nonzero
0.
0. p62 >0 --> force byte_swap for output
0. p63 =8 --> force 8-byte output
0. p64 =1 --> perturbation restart
0. p65 #iofiles (eg, 0 or 64); <0 --> sep. dirs
4.00000 p66 output : <0=ascii, else binary
4.00000 p67 restart: <0=ascii, else binary
0. p68 iastep: freq for avg_all
0.
0.
0.
0.
0.
0. p74 verbose Helmholtz
0.
0.
0.
0.
0.
0.
0.
0.
0.
0. p84 !=0 --> sets initial timestep if p12>0
0. p85 dt ratio if p84 !=0, for timesteps>0
0. p86 reserved
0.
0.
0.
0.
0.
0.
20.0000 p93 Number of previous pressure solns saved
5.00000 p94 start projecting velocity after p94 step
5.00000 p95 start projecting pressure after p95 step
1.00000 p96 u0 translational velocity (in .usr)
0.30000 p97 v0 translational velocity (in .usr)
0.
3.00000 p99 dealiasing: <0--> off/3--> old/4--> new
0.
0. p101 No. of additional filter modes
1.00000 p102 Dump out divergence at each time step
.000 p103 weight of stabilizing filter (.01)
4 Lines of passive scalar data follows2 CONDUCT; 2RHOCP
1.00000 1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000
13 LOGICAL SWITCHES FOLLOW
T IFFLOW
F IFHEAT
T IFTRAN
T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)
F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)
F IFAXIS
F IFSTRS
F IFSPLIT
F IFMGRID
F IFMODEL
F IFKEPS
F IFMVBD
F IFCHAR
More information about the Nek5000-users
mailing list