Define side of a surface

Topics related to boundary conditions.

Define side of a surface

Postby bouzias » Fri Dec 20, 2013 1:34 pm

Hello!

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!
bouzias
 
Posts: 5
Joined: Fri Dec 20, 2013 7:17 am

Re: Define side of a surface

Postby corey_ernst » Tue Jan 07, 2014 4:52 pm

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)
corey_ernst
 
Posts: 12
Joined: Tue Jan 07, 2014 4:27 pm

Re: Define side of a surface

Postby bouzias » Mon Jan 13, 2014 3:19 am

Yes that's exactly what I want to do.
bouzias
 
Posts: 5
Joined: Fri Dec 20, 2013 7:17 am

Re: Define side of a surface

Postby bouzias » Mon Mar 10, 2014 7:08 am

Up to recent! Any ideas anyone?
bouzias
 
Posts: 5
Joined: Fri Dec 20, 2013 7:17 am

Re: Define side of a surface

Postby statistics » Fri May 23, 2014 12:46 am

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
statistics
 
Posts: 1
Joined: Thu May 01, 2014 2:34 am

Re: Define side of a surface

Postby csimsoft Moderator » Mon Aug 31, 2015 3:33 pm

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} <range>SideSet <sideset_id> [Add] Surface <range> [wrt {{Volume|Surface} <range>|all}] [FORWARD|Reverse|Both]
You do not have the required permissions to view the files attached to this post.
csimsoft
csimsoft Moderator
 
Posts: 249
Joined: Tue Feb 14, 2012 5:59 pm
Location: American Fork, UT


Return to Boundary Conditions

Who is online

Users browsing this forum: No registered users and 1 guest

cron

Search / Mark

User Control Panel

Login

Who is online

In total there is 1 user online :: 0 registered, 0 hidden and 1 guest (based on users active over the past 5 minutes)
Most users ever online was 58 on Mon Jun 26, 2017 4:19 pm

Users browsing this forum: No registered users and 1 guest