[Nek5000-users] Periodicity with curvature

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

Dear Nek Users,

I am trying to get periodicity to work in the lastest version of Nek. I've
made an extremely simple mesh within prenek (prex) in polar coordinates and
cannot seem to run this case. When I run genmap, I get the error below. Did
I do something wrong in setting up this mesh?



>> genmap
 reading a.rea
Input mesh tolerance (default 0.2):
NOTE: smaller is better, but generous is more forgiving for bad meshes.
0.01
 reading mesh data ...
 start locglob_lexico:           4           2           8
 1.0000000000000000E-002
 locglob:           1           1           8
 locglob:           2           5           8
 locglob:           1           6           8
 locglob:           2           6           8
 done locglob_lexico:           6           6           8           2
 start periodic vtx:           2           6
     1     4     4  Matrix:    x0
  1   x0    1.0000E+00  1.0000E+00  5.4347-323  4.9407-324
  2   x0   -3.5355E-01  3.5355E-01  6.9423-310  6.9423-310
  3   x0    3.5355E-01 -3.5355E-01  6.9529-310  6.9529-310
  4   x0    6.9423-310  6.9529-310  2.3477-313  3.1632-317
     2     4     4  Matrix:    x1
  1   x1    1.0000E+00  1.0000E+00  6.3240-322  0.0000E+00
  2   x1   -3.5355E-01  3.5355E-01  0.0000E+00  6.9423-310
  3   x1   -3.5355E-01  3.5355E-01  6.9423-310  6.9529-310
  4   x1    6.9423-310  6.9423-310  6.9423-310  6.9529-310
     1     4     4  Matrix:    z0
  1   z0    1.0000E+00  1.0000E+00  5.4347-323  4.9407-324
  2   z0   -2.1213E+00 -1.4142E+00  6.9423-310  6.9423-310
  3   z0    2.1213E+00  1.4142E+00  6.9529-310  6.9529-310
  4   z0    6.9423-310  6.9529-310  2.3477-313  3.1632-317
     2     4     4  Matrix:    z1
  1   z1    1.0000E+00  1.0000E+00  6.3240-322  0.0000E+00
  2   z1    1.4142E+00  2.1213E+00  0.0000E+00  6.9423-310
  3   z1    1.4142E+00  2.1213E+00  6.9423-310  6.9529-310
  4   z1    6.9423-310  6.9423-310  6.9423-310  6.9529-310
           1 3  0  1.00000000E-03  2.12132000E+00 abort: FACE MATCH FAIL
           2 4  3  2.12132000E-03  1.00000031E+00 abort: FACE MATCH FAIL

           0  quit





The mesh within my .rea file is as follows:
   20.0000       20.0000      -8.00000      -4.00000
 XFAC,YFAC,XZERO,YZERO
 **MESH DATA** 1st line is X of corner 1,2,3,4. 2nd line is Y.
         2         2         2 NEL,NDIM,NELV
      ELEMENT          1 [    1A]    GROUP     0
  -2.121320      -1.414213     -0.8742276E-07 -0.1311341E-06
   2.121320       1.414213       2.000000       3.000000
      ELEMENT          2 [    1B]    GROUP     0
 -0.1311341E-06 -0.8742276E-07   1.414213       2.121320
   3.000000       2.000000       1.414213       2.121320
  ***** CURVED SIDE DATA *****
         4 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE
  2  1  -2.00000       0.00000       0.00000       0.00000       0.00000
 C
  4  1   3.00000       0.00000       0.00000       0.00000       0.00000
 C
  2  2  -2.00000       0.00000       0.00000       0.00000       0.00000
 C
  4  2   3.00000       0.00000       0.00000       0.00000       0.00000
 C
  ***** BOUNDARY CONDITIONS *****
  ***** FLUID   BOUNDARY CONDITIONS *****
 P    1  1   2.00000       3.00000       0.00000       0.00000       0.00000
 W    1  2   0.00000       0.00000       0.00000       0.00000       0.00000
 E    1  3   2.00000       1.00000       0.00000       0.00000       0.00000
 W    1  4   0.00000       0.00000       0.00000       0.00000       0.00000
 E    2  1   1.00000       3.00000       0.00000       0.00000       0.00000
 W    2  2   0.00000       0.00000       0.00000       0.00000       0.00000
 P    2  3   1.00000       1.00000       0.00000       0.00000       0.00000
 W    2  4   0.00000       0.00000       0.00000       0.00000       0.00000
  ***** NO THERMAL BOUNDARY CONDITIONS *****

[image: Inline image 2]
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180113/158333a2/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 7486 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180113/158333a2/attachment.png>