[Nek5000-users] Trouble using cylindrical mesh

[nek5000-users at lists.mcs.anl.gov]

Jean-Christophe,

Please try the latest repo for the nek source - there was a formatting
issue with reading the input.

cd nek5_svn/trunk
mv nek nek487
svn update nek

should work.


Paul


On Tue, 20 Apr 2010, nek5000-users at lists.mcs.anl.gov wrote:

> Hello all,
>
> I have a few problems when using the blah.box template from the Nek 5000 
> Primer guide for cylindrical meshes.
> Here is the error message I get:
>
> ===================================================================================
>
> Number of processors:           1
> REAL    wdsize      :           8
> INTEGER wdsize      :           4
>
>
> Beginning session:
> /home/loiseau/BricoNek/Cylindre/Cylindre.rea 
>
>
> timer accuracy:   0.0000000E+00 sec
>
> read .rea file
> nelgt/nelgv/lelt:       72       72       80
> lx1  /lx2  /lx3 :        6        4        4
>
> mapping elements to processors
>          0          72          72          72          72  NELV
> RANK     0 IEG       1       2       3       4       5       6       7 
> 8
>                     9      10      11      12      13      14      15 
> 16
>                    17      18      19      20      21      22      23 
> 24
>                    25      26      27      28      29      30      31 
> 32
>                    33      34      35      36      37      38      39 
> 40
>                    41      42      43      44      45      46      47 
> 48
>                    49      50      51      52      53      54      55 
> 56
>                    57      58      59      60      61      62      63 
> 64
>                    65      66      67      68      69      70      71 
> 72
> done :: mapping elements to processors
>
> ERROR READING MESH DATA NEAR ELEMENT      1
> ABORTING IN ROUTINE RDMESH.
>
> call exitt: dying ...
>
> backtrace(): obtained 9 stack frames.
> ./nek5000 [0x551d88]
> ./nek5000 [0x653f02]
> ./nek5000 [0x44202e]
> ./nek5000 [0x43d215]
> ./nek5000 [0x403b32]
> ./nek5000 [0x403a8e]
> ./nek5000 [0x40389c]
> /lib64/libc.so.6(__libc_start_main+0xf4) [0x3393c1d974]
> ./nek5000 [0x4037a9]
>
> total elapsed time             :   0.00000E+00 sec
> total solver time incl. I/O    :   0.00000E+00 sec
> time/timestep                  :   0.00000E+00 sec
> CPU seconds/timestep/DOF       :   0.00000E+00 sec
>
> ===================================================================================
>
> Attached are all the files I've been using.
>
> I also have a question concerning the descriptors line in the .box file.
>
> ===================================================================================
>
> Cyl.rea
> 2                     Spatial dimension
> 1                       number of fields
> #
> #
> #
> Y                       cYlinder
> 3  -24  1               nelr, nel-theta, nelz
> .5  .3                  x0, y0
> ccbb                    descriptors: c-cyl, o-oct, b-box
> .5  .55  .7  .8         r0 ... r_nelr
> 0  1  1                  theta0/2pi theta1/2pi ratio
> v  ,W  ,E  ,E  ,      BC
>
> ===================================================================================
>
> What does the ccbb line stands for ?
> If I wanna mesh the cross section of a cylindrical pipe and then extrude it 
> to a 3D mesh using n2to3, what descriptors should I use ?
>
> Regards,
>
> Jean-Christophe
>