[Gmsh] Physical Entities & Tags problem in output file version 2.0 !

[Christophe Geuzaine]

thierry.gilles at tiscali.nl wrote:
> Bonjour (again :-)
> 
> Me revoici avec un nouveau problème GMSH..
> 
> Ce que j'aimerais faire, c'est associer aux éléments d'un maillage 
> (surfaces en prio, lignes si possible) des étiquettes permettant 
> d'identifier de quels objets physiques ces éléments sont la frontière. 
> Pour des surfaces, elles sont toujours la frontière de 2 objets (dont 
> l'un est parfois le vide), par contre pour des lignes, ça peut 
> allègrement dépasser les 2 (comme tu verras dans mon exemple simplissime 
> de deux pyramides accolées)
> 
> Ben oui, la solution évidente c'est de définir des "Physical Surface" et 
> puis les tags prévus dans la version fichier de sortie 2.0 !!
> 
> J'ai essayé toute une journée, en vain... :(
> Définir des Physical Points, Lines, Surfaces, Volumes, mailler le tout, 
> sauver aux deux formats 1.0 et 2.0.. jamais je ne suis parvenu à obtenir 
> dans le fichier de sortie, en regard des éléments du maillage généré, 
> les entiers correspondant aux diverses "Physical Entities" que j'ai 
> accolées aux éléments géométriques générateurs du maillage.. Non, au 
> lieu de ça, dans la version 2.0, je n'obtiens que 2 tags identiques 
> (jamais plus que 2..) correspondant à l'"Elementary", rien à propos des 
> Physical..

Did you check the "Save all elements" option in the dialog? You
shouldn't: checking this will discard all the physicals. (I've changed
the label in the cvs version to make it more explicit.)

I've attached a slightly modified .geo file as well as the resulting
meshes. Just load one of the meshes and use "Physical" in
Tools->Visibility to examine the physicals: they look fine to me.


> 
> Je t'envoie les très petits fichiers de test que j'ai utilisé, ainsi que 
> les questions en vrac qui me sont venues à l'esprit et dont je n'ai pas 
> trouvé la réponse dans le manuel pdf ni en faisant plein d'essais.
> 
> Question 1 :
> Comment forcer la sauvegarde de plus de 2 tags ? Et où et comment 
> définir ces tags ?

That's not how it's supposed to work: Gmsh always saves 2 tags
(elementary+physical), or 3 when the mesh is partitioned.

The "extra tag" facility is there for use by other applications. Gmsh
will simply ignore those tags.


> 
> Question 2 :
> Si c'est une meilleure alternative aux tags pour mon problème de 
> frontière de volumes, comment utiliser les partitions ? Comment 
> retrouver cette info dans le fichier .msh ?
> 
> Question 3 :
> Mon idée est de grouper en Entités Physiques les surfaces enfermant des 
> volumes de même propriété physique. Ce faisant, chaque surface 
> appartiendra à deux Entités Physiques. J'aimerais ensuite retrouver 
> cette info (les deux entités physiques) dans le fichier de sortie en 
> regard de chaque élément "surface" du maillage.. Est-ce selon toi la 
> meilleure solution ?
> 

I don't know... If a mesh element belongs to two physical entities, then 
this element will appear twice (with a different number) in the .msh file.

C.


>  
> 
> Merci encore une fois par avance pour tes commentaires !
> 
> Thierry GILLES
> 
> 


-- 
Christophe Geuzaine
Assistant Professor, Case Western Reserve University, Mathematics
http://www.case.edu/artsci/math/geuzaine
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: DoublePyramid2.msh
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051214/4963a19d/attachment.msh>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: DoublePyramid.msh
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051214/4963a19d/attachment-0001.msh>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: DoublePyramid.geo
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051214/4963a19d/attachment.geo>