Btw.: the wrong triangle from the error message is part of a mesh produced by the NodeReordering tool?
Then there might be a bug in that tool. Could you open a bug report in gitlab and add your input/output meshes?
You are using msh2vtu from OGSTools. Which OGSTools version do you have? Maybe an update to the latest version changes something. There could have been some changes in msh2vtu, but I’m not entirely sure.
In the meantime a workaround could be the following (if you have linear geometries, i.e., no curved surfaces that you want to resolve as precisely as possible):
- generate a mesh with 1st order elements with gmsh
- convert that to vtu (
msh2vtu) - do NodeReordering on all meshes from step 2
- convert the meshes from step 3 to 2nd order: createQuadraticMesh
- “connect” boundary meshes and the bulk mesh via
identifySubdomainsfor each boundary mesh.