[MOAB-dev] Neumann boundary condition
    Iulian Grindeanu 
    iulian at mcs.anl.gov
       
    Wed Dec 11 15:25:00 CST 2013
    
    
  
programmatically? 
a neumann set in moab is usually a set of boundary faces (quads in a hexahedral mesh, or triangles in a tetrahedral mesh) 
what kind of solver do you have? 
if you know the tag value, you can obtain the set that has that tag value, with something like 
R ange neumann_sets; 
Tag neu_tag; 
int value_neumann_tag= 100; // your value , set it accordingly 
mb->tag_get_handle(NEUMANN_SET_TAG_NAME, neu_tag, 1, MB_TYPE_INTEGER); 
mb->get_entities_by_type_and_tag(0, MBENTITYSET, &neu_tag, &value_neumann_tag, 1, neumann_sets); 
neumann_sets[0] should be your set. 
then you can get entities from it with 
Range entities; 
mb->get_entities_by_handle(neumann_sets[0], entities); 
Entities should be your quads or triangles. it is up to the specific solver syntax to set pressure on them (or other BCs) 
Maybe I didn't understand your question 
Iulian 
----- Original Message -----
| Hello MOAB-dev,
| Can you please explain how can i set a certain Neumann boundary
| condition
| for selected set of faces via NEUMANN_SET tag?
| For example the normal derivative of pressure
| Thanks
| Anton
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20131211/84db3382/attachment.html>
    
    
More information about the moab-dev
mailing list