[petsc-users] dmplex normal vector incorrect for periodic gmsh grids

[Praveen C]

Hello

In the attached test, I read a small grid made in gmsh with periodic bc.

This is a 2d mesh.

The cell numbers are shown in the figure.

All faces have length = 2.5

But using PetscFVFaceGeom I am getting length of 7.5 for some faces. E.g.,

face: 59, centroid = 3.750000, 2.500000, normal = 0.000000, -7.500000
===> Face length incorrect = 7.500000, should be 2.5
support[0] = 11, cent = 8.750000, 3.750000, area = 6.250000
support[1] = 15, cent = 8.750000, 1.250000, area = 6.250000

There are also errors in the orientation of normal.

If we disable periodicity in geo file, this error goes away.

Thanks
praveen


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221213/42de8afa/attachment-0004.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dmplex.c
Type: application/octet-stream
Size: 4821 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221213/42de8afa/attachment-0002.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221213/42de8afa/attachment-0005.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mesh.png
Type: image/png
Size: 10967 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221213/42de8afa/attachment-0001.png>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221213/42de8afa/attachment-0006.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ug_periodic.geo
Type: application/octet-stream
Size: 525 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221213/42de8afa/attachment-0003.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221213/42de8afa/attachment-0007.html>