Exploitation des résultats & Comparaison

A présent nous avons les résultats théoriques ainsi que ceux de simulation. Un dernier point à préciser est que pour les différentes grandeurs que nous allons comparer, nous n'avons pas accès directement par Parafoam. Pour chaque grandeur, il y en a des manipulations à a faire.

Le nombre de Nusselt :

Numériquement, nous avons les profils de température en tout point, pour calculer le nusselt, il suffit de calculer le gradient de température au niveau de la paroi. Pour cela, nous avons calculé la différence entre la température de la paroi et celle de la maille qui suit suivant la direction y. Pour que le calcul soit correct, il faut avoir des mailles suffisamment fines.

Remarque Pour tout calcul sur Parafoam, il ne faut surtout pas exporter les résultats pour des lignes données utilisant plotoverline, car ceci ne permet d'avoir que des valeurs pour un pas d'espace de 0.005 quelque soit le maillage, l'une des faiblesses de cet outil. Pour surmonter cette difficulté, il faut exporter les résultats de tout le domaine, et après faire un tri sur les zones qui intéressent l'utilisateur.

La comparaison des deux approches donne :

                                 

Nous remarquons que les deux profils ne collent pas du tout, quelque soit le maillage. L'origine de cet écart ne provient pas de la simulation ni du maillage, car nous avons bien vérifié que la simulation est correcte, les profils expérimentaux à l'entrée de la couche limite turbulente collent bien avec les profils à la sortie du premier bloc. Le problème vient plutôt de la définition du gradient de température que nous avons considéré pour calculer le Nusselt numériquement. En effet, théoriquement, on montre que le Nusselt augmente avec x, ce qui est bien physique. Cependant, Si pour notre définition du gradient, comme la température suivant x augmente avec x, donc le gradient diminue forcément, par suite Le Nusselt diminue avec x. Nous avons beau cherché pour résoudre ce problème, mais malheureusement nous nous sommes pas arrivés à le résoudre. Les idées que nous avons eu pour surmonter cette difficulté étaient de chercher si OpenFOAM permet de calculer automatiquement ces grandeurs, en les codant par exemple (comment ce qu'on peut faire pour les coefficients de portance et de traînée). L'autre idée était de se servir du Python pour les coder, car ce langage s'avère très puissant pour l'automatisation des calculs et l'obtention des résultats. Malheureusement, nous n'avons pas eu le temps pour faire ces manipulations.

Vitesse maximale :

Pour cette grandeur, le calcul est moins délicat. La comparaison des deux approches nous donne :

 

 

                                      

 

Nous constatons que les deux profils ont la même forme avec un "léger" écart.

Epaisseur de la couche limite dynamique :

Concernant cette grandeur, c'était la galère pour la déterminer numériquement. En effet, si nous faisons un plotoverline pour tous les points du maillage, le calcul de cette épaisseur est impossible avec Parafoam, vu que cet outil considère des pas d'espace de 0.005m et non pas les points du maillage, ce qui est quasiment pas faisable, puisque cette épaisseur est de l'ordre du millimètre. Nous avons essayé d'exporter les résultats de tout le domaine, cependant, vu le nombre de points du maillage, et les vitesses qui étaient très sensibles, le calcul de cette grandeur était imprécis.