<html>
<body>
<br>
Bonjour,<br><br>
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.<br><br>
Je crois avoir décelé un problème lors de l'extrusion d'un maillage.<br>
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.<br><br>
La maillage 2D se compose de plusieurs parties.<br><br>
4 Ruled Surfaces transfinies circulaires autour du cylindre.
(partie // BOUNDARY LAYER // dans le fichier Cylinder.geo )<br>
1 Plane Surface carrée autour du
cylindre.
(partie // SQUARE
//
dans le fichier Cylinder.geo )<br>
1 Plane Surface rectangulaire transfinie à droite du
carré. (partie // RECTANGLE
// dans le
fichier Cylinder.geo )<br>
4 Plane Surfaces pour agrandir le maillage le plus possible. (partie //
EXTERIOR
//
dans le fichier Cylinder.geo )<br><br>
Ensuite le maillage est extrudé sur 4 couches sur une hauteur totale
z=0.1. (partie // EXTRUSION // dans le fichier Cylinder.geo)<br><br>
Le problème est le suivant:<br><br>
Dans le fichier Cylinder.msh les noeuds 25157 et 25931 ont exactement les
mêmes coordonnées (<font face="Courier New, Courier">0.5125301401598015
0.4646092699136419 0.1</font>)<br><br>
Ce noeud unique mais dédoublé appartient à la surface supérieure créée
par extrusion.<br><br>
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).<br><br>
Je constate que<br><br>
-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 !<br>
-le noeud 25931 appartient à 5 faces (les éléments
16453,16821,16822,17477,17478 de la surface physique 3 et surface
elementaire 440) <br>
<x-tab> </x-tab>
et à 6 volumes (les éléments 35378,36410,37882,37886,40506,40510 du
volume physique 1 et volume elementaire 5)<br><br>
Si ce noeud n'en formait qu'un seul, il n'y aurait donc aucun
problème.<br><br>
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 !<br><br>
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.<br><br>
Est-ce un bug ou ai-je fait quelque chose d'incorrect ?<br><br>
Tant que ce problème subsiste je n'arrive pas à importer ce maillage dans
mon code de calcul.<br><br>
Merci de m'apporter votre aide.<br>
Je suis en tout cas très heureux d'utiliser ce mailleur gratuit ! de
qualité. <br>
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 !<br><br>
Bonne journée,<br><br>
<br>
<x-tab> </x-tab>Didier
Vigneron<br><br>
<br>
<x-sigsep><p></x-sigsep>
--------------------------------------------------------<br>
Didier Vigneron - Assistant<br>
<a href="http://www.ulg.ac.be/aerodyn/fr/people/vigneron/index.html" eudora="autourl">
http://www.ulg.ac.be/aerodyn/fr/people/vigneron/index.html<br>
</a>Université de Liège<br>
Département ASMA (Aérodynamique)<br>
Chemin des chevreuils, 1, Bât B52/3<br>
B-4000 Liège<br>
--------------------------------------------------------<br>
Phone : (+32) (0)4 366 94 39<br>
Fax : (+32) (0)4 366 91 36<br>
<a href="http://www.ulg.ac.be/aerodyn" eudora="autourl">
www.ulg.ac.be/aerodyn<br>
</a><a href="http://www.ulg.ac.be/turbo" eudora="autourl">
www.ulg.ac.be/turbo<br>
</a>--------------------------------------------------------<br>
</body>
</html>