[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