libmesh icon indicating copy to clipboard operation
libmesh copied to clipboard

Unable to read and set up second order mesh created using Salome (ERROR: UNV-element type 22 not supported)

Open raghwendra-78 opened this issue 3 years ago • 2 comments

I am trying to solve Laplace equation with curved boundary. Adapted example file introduction_ex4.C (the adapted file is attached below with .txt extension). Second order mesh was created using "Salome" software. After grouping faces to impose boundary conditions, mesh was exported in .unv format (compressed file of mesh exported is attached). When I try to execute the program, I get following error: UNV-element type 22 not supported (output at terminal is attached as log-second-order-mesh.txt).

introduction_ex4.C.txt test-second-order.unv.zip log-second-order-mesh.txt

Thanks, Raghwendra Kumar

raghwendra-78 avatar Oct 20 '21 09:10 raghwendra-78

Well, the error is correct, at least! There's no 22 in our switch statement there.

http://victorsndvg.github.io/FEconv/formats/unv.xhtml seems to imply that we could interpret that as a libMesh EDGE3? If so then you might try adding a case 22: block to unv_io.C (around line 670) and see if you can get it to work for you; if so then we'd love to get a patch! I'm not sure what the permutation would be for UNV->libMesh node ordering, but it's almost certainly 0,1,2 or 0,2,1, so you could just try both and see what works.

roystgnr avatar Oct 20 '21 21:10 roystgnr

Hmm... reading further at that link, the "Local order of nodes" section makes me think you'd want to try 0,2,1 first.

roystgnr avatar Oct 20 '21 21:10 roystgnr