How do I extract all the elements on the grain boundaries?

[smallrain1996]

Is there a way to output the ids of all elements located on the grain boundaries? If you can tell me what command to use to create an Elset composed of all elements in the model that are located at the grain boundary, or how to output a statistical file that shows whether each element is an in-grain element or a grain boundary element?

[mkasemer]

Try: neper -M … -statelt elsetbody.

To be clear: this will contain 0 for elements which are on the grain boundary, 1 for elements which are separated from the grain boundary by a single layer of elements, etc.

To a question below, use neper -M ... -statelt body to find the same information relative to the domain, rather than the grain boundaries.

[smallrain1996]

I'm using the latest version of Neper, so you can try executing my code. It's true that all the elements are 0, and I checked with the tool.

[mkasemer]

Try refining the mesh further. Otherwise, @rquey will have to comment.

[smallrain1996]

Maybe it was a Bug, I used a mesh like 100 100 100, and so on.

[rquey]

It won't work with hex meshing...

[smallrain1996]

Crying instantly, this is really a very sad thing.

[charanjo]

Is there a way to extract elements on the surface of the domain?

[mkasemer]

Try: neper -M … -statelt body