[Nek5000-users] help for Element mesh mismatch

stefanke stefanke at lav.mavt.ethz.ch
Tue Sep 22 07:31:41 CDT 2009


Hi,

it looks like your mesh coordinates and connectivity do not match.
The given elemental boundary point (see err msg) is a shared D.O.F  
(the same physical point is represented multiple times) but there  
seems to be a difference in coordinates which is larger than the  
tolerance (1e-3).

hth,
Stefan


On Sep 21, 2009, at 1:07 PM, wangzhicheng wrote:

> Hi developers
>
> The problems with cavity flow in previous email has been solved.
> However, I meet another trouble today.
> I attached a document  that contain .rea file for the case of  flow  
> past a cylinder. The mesh is transformed  from the format used in  
> code Semtex to Nek5000, by  tool sem2nek with some change.
> But when I try to run this case, nek5000 give me errors like this:
>
> ****************************************************************************************************************************************************************
>
>  Number of processors:           1
>  REAL    wdsize      :           8
>  INTEGER wdsize      :           4
>
>
>   Beginning session:
>   /home/wzc/Code/nek5_svn/examples/2Dcylinder/cylinder01.rea
>
>
>  timer accuracy:   1.0967255E-06 sec
>
>  read .rea file
>  nelgt/nelgv/lelt:       54       54      300
>  lx1  /lx2  /lx3 :        8        6        8
>
>  mapping elements to processors
>            0          54          54          54          54  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
>  done :: mapping elements to processors
>
>            0  objects found
>  readat time  3.292083740234375E-003  sec
>  done :: read .rea file
>
>  setup domain topology
>    Right-handed check complete for      54 elements. OK.
>    setvert2d:           8           8           1
>    setupds2d:          8          0         69         69         
> 321        807
>    gs_setup: 0 unique labels shared
>    gs_init time 3.2997E-04 seconds   8        54      3456        
> 321         0
>            5  max multiplicity
>  done :: setup domain topology
>
>  call usrdat
>  done :: usrdat
>
>  generate geomerty data
>  vol_t,vol_v:   134.648547351613        134.648547351613
>  done :: generate geomerty data
>
>  call usrdat2
>  done :: usrdat2
>
>  regenerate geomerty data           1
>  vol_t,vol_v:   134.648547351613        134.648547351613
>  done :: regenerate geomerty data           1
>
>  verify mesh topology
>   -5.00000000000000        10.0000000000000       Xrange
>   -5.00000000000000        5.00000000000000       Yrange
>    0.00000000000000        0.00000000000000       Zrange
>  WARNING1 Element mesh mismatch at:
>  i,j,k,ie:    8    2    1        2
>  Near X =  0.43169461      0.47521817       0.0000000    , d:    
> 0.0000000     -0.15500086E-02   2.1757104
>
> call exitt: dying ...
>
> backtrace(): obtained 8 stack frames.
> ./nek5000 [0x4d040a]
> ./nek5000 [0x51f1a6]
> ./nek5000 [0x42af2d]
> ./nek5000 [0x403b86]
> ./nek5000 [0x40334c]
> ./nek5000 [0x61f8be]
> /lib/libc.so.6(__libc_start_main+0xf4) [0x7f8bddd421c4]
> ./nek5000 [0x403299]
>
> total elapsed time         :   0.12290E-01 sec
> total solve time incl. I/O :    0.0000     sec
> time/timestep              :    0.0000     sec
> ********************************************************************************************************************************************************************
>
> With my best wishes
>
>
> Zhicheng
>
>
>
>
>
>
> <cylinder01.tar>_______________________________________________
> 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