[Gmsh] Gmsh Bug ? Report
Christophe Geuzaine
cgeuzaine at ulg.ac.be
Fri Feb 9 10:23:30 CET 2007
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