Vérification de code pour le modèle k-ϵ standard de OpenFOAM

Dans un contexte industriel toujours plus compétitif, les industries sont constamment comparées les unes aux autres. Elles ont donc intérêt à se moderniser et investir dans des outils de pointe afin d'améliorer la qualité de leur procédé et de diminuer leurs coûts d'opération. Les outils de CFD répondent précisément à ce besoin, ce qui explique la demande grandissante pour ces outils dans l'industrie. Toutefois, la crédibilité des résultats obtenus reste un obstacle majeur pour l'application plus répandue des outils de CFD. Nous avons donc comme but de contribuer à la crédibilité des résultats en confirmant que les méthodes numériques, telles qu'implémentées dans le code, se comportent conformément à la théorie. Dans le cadre de ce travail, nous avons choisi d'employer le code OpenFOAM, un logiciel à code source libre fréquemment employé dans l'industrie. Il offre une impressionnante collection d'outils entièrement gratuits, répondant en grande partie aux besoins industriels actuels. Cependant, ce logiciel est nettement plus complexe à employer qu'un logiciel de CFD commercial. De plus, on retrouve peu d'informations détaillant la vérification rigoureuse du code. Il nous semble donc approprié de porter plus d'attention à la vérification du code OpenFOAM pour les multiples solveurs, modèles et conditions limites disponibles. Dans ce mémoire nous effectuerons la vérification du code OpenFOAM pour un écoulement turbulent et incompressible. En effet, ce type d'écoulement est représentatif de la majorité des situations rencontrées dans un contexte industriel. On modélise la turbulence avec le modèle k-ϵ standard à cause de sa simplicité d'utilisation et son faible coût de calcul. Ce modèle est un compromis intéressant pour sa précision et son coût de calcul. Nous commencerons par vérifier le code pour un écoulement très simple, soit celui de Poiseuille. Nous appliquerons ensuite la méthode de solutions manufacturées imitant un écoulement laminaire 1D de type couche limite, un écoulement libre 2D laminaire, un écoulement libre 2D turbulent, et finalement un écoulement de type couche cisaillée. Afin de mieux comprendre les écarts entre le taux de convergence théorique et le taux de convergence observé, nous poursuivrons la vérification en y ajoutant l'étape de vérification de simulation à l'aide de l'extrapolation de Richardson. Ces cas nous permettront de confirmer si le taux de convergence théorique des schémas numériques est respecté ou non.