[Gmsh] Gmsh Bug ? Report

[Christophe Geuzaine]

Hello Didier - Do you get the same duplicate nodes with the new version 
(Gmsh 2.0.0)?


Didier Vigneron wrote:
> 
> Bonjour,
> 
> J'utilise Gmsh pour créer des maillages non structurés pour la 
> résolution des équations de Navier-Stokes par la méthode des volumes finis.
> 
> Je crois avoir décelé un problème lors de l'extrusion d'un maillage.
> Je vous envoie le fichier Cylinder.geo que j'utilise. Il s'agit d'un 
> maillage autour d'un Cylindre. Le Cylindre est donc un trou.
> 
> La maillage 2D se compose de plusieurs parties.
> 
> 4 Ruled Surfaces transfinies circulaires autour du cylindre.  (partie // 
> BOUNDARY LAYER // dans le fichier Cylinder.geo )
> 1 Plane Surface   carrée autour du cylindre.                       
> (partie // SQUARE //                 dans le fichier Cylinder.geo )
> 1 Plane Surface   rectangulaire transfinie à droite du carré.  (partie 
> // RECTANGLE //           dans le fichier Cylinder.geo )
> 4 Plane Surfaces pour agrandir le maillage le plus possible. (partie // 
> EXTERIOR //              dans le fichier Cylinder.geo )
> 
> Ensuite le maillage est extrudé sur 4 couches sur une hauteur totale 
> z=0.1. (partie // EXTRUSION // dans le fichier Cylinder.geo)
> 
> Le problème est le suivant:
> 
> Dans le fichier Cylinder.msh les noeuds 25157 et 25931 ont exactement 
> les mêmes coordonnées (0.5125301401598015 0.4646092699136419 0.1)
> 
> Ce noeud unique mais dédoublé appartient à la surface supérieure créée 
> par extrusion.
> 
> De plus, lorsque j'inspecte le maillage, je remarque que ce noeud doit 
> appartenir à 6 volumes et à 6 faces frontières (les faces des 6 volumes).
> 
> Je constate que
> 
> -le noeud 25157 n'appartient qu'à une seule face (l'élément 16195 de la 
> surface physique 3 et surface elementaire 440) et à aucun volume !
> -le noeud 25931 appartient à 5 faces (les éléments 
> 16453,16821,16822,17477,17478 de la surface physique 3 et surface 
> elementaire 440)
>                                    et à 6 volumes (les éléments 
> 35378,36410,37882,37886,40506,40510 du volume physique 1 et volume 
> elementaire 5)
> 
> Si ce noeud n'en formait qu'un seul, il n'y aurait donc aucun problème.
> 
> Dans ce maillage je crois qu'il n'existe qu'un seul noeud dédoublé de ce 
> genre. Par contre lorsqu'on augmente le nombre de mailles en changeant 
> le paramètre N au début du fichier (en mettant N=1 plutôt que N=0.5 par 
> exemple) le nombre de noeuds dédoublés augmentent fortement !
> 
> Visiblement, le problème se pose aussi bien pour les prismes que pour 
> les hexaèdres. Je n'ai décelé ce problème que pour la noeuds qui se 
> trouvent sur la face supérieur créée par l'extrusion. Mais il se peut 
> que le problème arrive aussi pour des noeuds intérieurs , je n'en sais rien.
> 
> Est-ce un bug ou ai-je fait quelque chose d'incorrect ?
> 
> Tant que ce problème subsiste je n'arrive pas à importer ce maillage 
> dans mon code de calcul.
> 
> Merci de m'apporter votre aide.
> Je suis en tout cas très heureux d'utiliser ce mailleur gratuit ! de 
> qualité.
> Merci beaucoup d'avoir un jour eu l'idée de développer ce logiciel qui 
> doit, j'en suis convaincu, prendre beaucoup de votre temps !
> 
> Bonne journée,
> 
> 
>         Didier Vigneron
> 
> 
> --------------------------------------------------------
> Didier Vigneron - Assistant
> http://www.ulg.ac.be/aerodyn/fr/people/vigneron/index.html
> Université de Liège
> Département ASMA (Aérodynamique)
> Chemin des chevreuils, 1, Bât B52/3
> B-4000 Liège
> --------------------------------------------------------
> Phone : (+32) (0)4 366 94 39
> Fax     : (+32) (0)4 366 91 36
> www.ulg.ac.be/aerodyn
> <http://www.ulg.ac.be/aerodyn> www.ulg.ac.be/turbo
> <http://www.ulg.ac.be/turbo>--------------------------------------------------------
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh


-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine