[harshabose]

I am currently facing an issue with my mesh generation process. I have developed a C++ code that generates polyMesh files (points, faces, cells, owners etc) to create a purely hexahedral mesh. The code functions correctly, and the polyMesh files are created as expected. However, when I open these polyMesh files in paraView, I notice the presence of additional points and cells.

According to my C++ code, the mesh should consist of exactly 36,080 points and 17,550 cells. However, upon inspection in paraView, I observe 36,158 points and 17,940 cells. This discrepancy is perplexing, as the generated polyMesh files clearly indicate the correct number of points and cells. Furthermore, the additional points seem to be created at coordinates (X, Y, Z) = (..., ..., 0.00025), while my intended coordinates only include (X, Y, Z) values of (..., ..., 0) or (..., ..., 0.0005).

I have also checked the results from the "checkMesh" command, which report the correct number of points and cells. However, even in the "checkMesh" output, I notice the presence of polyhedra cells and points at coordinates (X, Y, Z) = (..., ..., 0.00025).

Mesh stats 
    points:           36080
    internal points:  0
    faces:            70689
    internal faces:   34611
    cells:            17550
    faces per cell:   6
    boundary patches: 10
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     17472
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     78
    Breakdown of polyhedra by number of faces:
        faces   number of cells
            5   39
            7   39

I checked my code, even manually checked the polyMesh files. Was not able to find the problem. I need to know why these extra points and cells are being created. I am including few images and my files below.


Please find attached a few images and my mesh files for reference.

I would greatly appreciate any assistance or insights into resolving this issue.

Thank you in advance for your help

Extra Points.jpg
Extra Points

Extra Cells.jpg
Extra Cells
Case File Link (Google Drive)

[Yann]

Hello Lyly,

Everything seems fine to me! (cf screenshot)
In paraView, just uncheck "Decompose polyhedra" when loading your case.

Cheers,
Yann

[harshabose]

Hello Yann,

That works. Thank you for the help. Really appreciate it.😊
There was still the problem of polyhedral cells appearing for some reason.

Face Count.jpg
Face Count

I was able to convert these polyhedra to hexahedra using the "Convert Polyhedral Cells" filter. Seems to work fine.