[Nek5000-users] Problem converting 2d to 3d

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Tue Jan 14 10:28:07 CST 2014


Hi,

I created a 2d rectangular mesh using genbox,


then I tried to used n2to3 to convert this into a cylindrical pipe 
using  circular sweep but I get errors while loading the .rea file .

could some one help me fix this

Thank you

********************************
base.rea
2                      spatial dimension
1                      number of fields
#
#  Eddy case
#
Box
-50  -5                     nelx,nely,nelz for Box
0   10.     1.0             x0 x1 ratio
0   1      1.0                y0 y1 ratio
v  ,O  ,  ,W  ,              bc's
***************************************


Input for n2to3 script
*****************************************
Input old (source) file name:
eddy_uv
  Input new (output) file name:
fg
  For ASCII output only: 0; For .rea/.re2: 1
1
  input number of levels: (1, 2, 3,...; < 0 for circular sweep.):
-2
  input z min:
0
  input z max:
1
  input gain (0=custom,1=uniform,other=geometric spacing):
1
  This is for CEM: yes or no:
no
  Enter Z (5) boundary condition (P,v,O):
W
  Enter Z (6) boundary condition (v,O):
W
  this is cbz: W  W   <---
  F           0           1  ifheat

         2000 elements written to fg.rea

         2000 elements written to fg.re2


when I use genmap or try to launch the mesh in prex I get error
*****************************************************************
Input (.rea) file name:
fg
Input mesh tolerance (default 0.2):
NOTE: smaller is better, but generous is more forgiving for bad meshes.
0.2
  reading .rea file data ...
  reading .re2 file data ...
  start locglob_lexico:           8        2000       16000 
0.20000000000000001
  locglob:           1           1       16000
  locglob:           2         104       16000
  locglob:           3        1192       16000
  locglob:           1        3012       16000
  locglob:           2        3025       16000
  locglob:           3        3031       16000
  locglob:           1        3031       16000
  locglob:           2        3031       16000
  locglob:           3        3031       16000

      1     2     4  Matrix:  SELF!!
        1 SELF!!        1        2        1       12
        2 SELF!!       23       24       23       34
  cont: SELF!!           1   ??

  ABORT: SELF-CHK            1           5           1           0
  Try to tighten the mesh tolerance!

            0  quit
*****************************************************************************************




More information about the Nek5000-users mailing list