[MOAB-dev] Neumann boundary condition

[kanaev at ibrae.ac.ru]

 So, as i can understand, any value of NEUMANN_SET tag corresponds to
some certain rule implemented in solver  
 I've got one more question - having only nodes and hexagonals in
.h5m mesh file, what is the best way to find and store the opposite
face from the neighbouring cell for each hexagonal's face?  
 Thanks   
 On Wed, 11 Dec 2013 15:25:00 -0600 (CST), Iulian Grindeanu  wrote:  
p { margin: 0; }  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
           Range 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
-------------------------
  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/94e73811/attachment.html>