c 02/06, àj 03/06/ et /03/09

Solution de convection naturelle dans une cavité chauffée

quelques éléments en rapports avec la "pal 06" de l'ENSTA. En fin de page, la version gerris du problème.


La cavité est refroidie à gauche et chauffée à droite.
Les murs inférieurs et supérieurs sont adiabatiques.
Le Grashof est de 106.
Calcul avec FreeFEM++, pour Gr =1000000
On observe les deux couches limites montantes et descendantes
à droite et à gauche, ainsi que la stratification de température induite.





Tracé de la solution numérique:
la température est stratifiée à peu près linéairement.
Profils de vitesse et de température
on voit bien les langues chaudes et froides
iso-psi:
gravité verticale
et iso-T

le code FreeFEM++ de la cavité chauffée



Tracé de la vitesse et de la température calculées ainsi que de la solution de "Gill 66" en exp fois cos et en exp fois sin! l'épaisseur de couche limite est en Gr-1/4

Utilisation de gerris

le code gerris de la cavité chauffée. et le code de GfsView .


en lançant ce code ( gerris2D KVITboussinesq.gfs ¦ gfsview2D KVITboussinesq.gfv ) on obtient

on observe la première phase de réchauffement symétrique t=0.2 0.4 0.6 0.8 1 et 2:

puis le basculement et la rotation... et par ex au temps t=34:



on pourrait reprendre les tracés ci dessus de température et de pression