pre-pos_processing_gmsh-herezh

maillage Caddie_Gmsh.msh: semble correct

Au niveau du fichier info:

1) Il faut penser à supprimer les noeuds non référencés par les
éléments. Ici il faut mettre après la définition du maillage:

1. -- def maillage
   < Maillage_nevez.her
   < Maillage_nevez.lis

1. --- demande de suppression des noeuds non references ----
   suppression_noeud_non_references_

ce qui conduit à la suppression de 4 noeuds qui n'appartiennent pas à des éléménts donc
qui n'ont aucune raideur affectée. Si ces noeuds n'étaient pas bloqués cela conduit systématiqueme
à une divergence.

2) Ensuite il faut penser à obtimiser la largeur de bande, car par défaut gmsh produit une
numérotation qui n'optimise pas la largeur de bande ce qui conduit à de très grosses
matrices et donc des temps d'exécution importants.
Pour l'optimisation, on indique dans le .info:

1. -- def maillage
   < Maillage_nevez.her
   < Maillage_nevez.lis

1. --- demande de suppression des noeuds non references ----
   suppression_noeud_non_references_

1. demande de renumerotation
   renumerotation_des_noeuds_

On observe que le calcul est beaucoup ... plus rapide !!

3) Si l'on veut une sortie des résultats pendant le calcul, il faut utiliser l'option ad hoc
(FREQUENCE_SORTIE_FIL_DU_CALCUL 1) et au bonne endroit !!

4) le problème fondamentale qui fait que le calcul ne converge pas est a priori le fait
que les normales aux éléments ne sont pas correctement orientées. La moitié bas du
maillage est orientée vers les z positifs et la moitié haute vers les z négatifs.
Les raideurs au niveau de la jonction s'annullent et donc on obtient une matrice qui
est incorrecte !!
Il faut donc réorienter les normales, soit directement via le mailleur ou le modeleur,
soit via un utilitaire d'herezh qui est prévu pour cela et qui travaille sur des facettes
triangulaire à la base.

Merci de me tenir au courant si ce n'est pas le bon diagnostique !

Bonjour,

La résolution du problème des normales aux éléments n'a pas résolu notre problème, on obtient ce type d'erreur : 
"warning: exception generee par un element mais dont la prise en compte  n'est pas prevu !, on ne fait rien et on continue le calcul erreur le pas de temps 9.34718e-17 est trop faible pour un calcul correcte de la vitesse moyenne".

De plus, lorsqu'on observe dans Caddie_Gmsh.msh les "elementary_tag" on a quelques éléments qui n'ont pas la même couleur que les autres (voir "Elementary_tag.png"). Peut être que l'erreur vient de là ?

Par contre, le jacobien n'est, à priori, plus négatif.

Merci d'avance.

Cordialement,

Adrien LAUNAY
Thomas JACOPIN

Je vous ai mis une version qui fonctionne. Le plus simple est que vous regardiez les paramètres que j'ai utilisés. J'ai également fait le calcul avec une précision plus importante (5.e-3 -> 1.e-3) et j'ai également multiplié le chargement par 10.