Define side of a surface


I’d like to create a surface triangular mesh of a non-manifold geometry and pass to the exodus file the information that describes for each triangle that lies on an interface of volumes, to which volume of the interface does it point to. What is the best strategy to do this task?

What I’ve tried to do so far is: after creating the geometry and merging all the surfaces, then I scan each pair of volumes and create the corresponding groups (Let’s say that we have 5 volumes. The command: group 5 add surface in volume 1 except surface in volume 2 3 4 with is_merged. This command will add to group 5 the surfaces that lie on the boundary of volumes 1 and 5). Consequently, I define the orientation of the group (for example: surface in group 5 normal volume 1) and eventually I create a block for each group. Therefore I can finally write to the exodus file as an attribute the two volumes that each mesh triangle lies on the interface of as well as its orientation.

However, the resulting number of blocks is huge (obviously).

Is there a more formal way/strategy to complete this task?

Thanks in advance for any kind of help!

Keep up the good work!

So just to clarify, you would like:

  1. Each surface that is shared between volumes to be its own block.
  2. Each block needs to also have a attribute indicating which volume the elements in the block reference in a forward sense? (element normals and surface normals point in the same direction)

Yes that’s exactly what I want to do.

Up to recent! Any ideas anyone?

Yes anybody able to do this. I would love to see a well defined solution.

GSA Search Engine Ranker Download | GSA Captcha Breaker download | GSA Search Engine Ranker Coupon

The attached script should do what you need. You may want sidesets in addition to blocks. Sidesets reference the ‘sides’ of elements and in the case of surfaces on two volumes, you have the options of using the ‘wrt’ option:

SideSet <sideset_id> [ADD|Remove]{Curve|Surface|Edge|Face|Tri|Group} SideSet <sideset_id> [Add] Surface [wrt {{Volume|Surface} |all}] [FORWARD|Reverse|Both]
merged_surf_blocks.jou (996 Bytes)