| Résumé :
|
Une approche multidimensionelle pour la résolution numérique des équations d'Euler et de Navier-Stokes sur maillages non-structurés est proposée. Dans une première partie, un exposé complet des schémas de distribution, dits de "fluctuation-splitting" ,est décrit, comprenant une étude comparative des schémas décentrés, positifs et de 2ème ordre, pour résoudre l'équation de convection à coefficients constants, ainsi qu'une étude théorique et numérique de la précision des schémas sur maillages réguliers et distordus. L'extension à des lois de conservation non-linéaires est aussi abordée, et une attention particulière est portée au problème de la linéarisation conservative. Dans une deuxième partie, diverses discrétisations des termes visqueux pour l'équation de convection-diffusion sont développées, avec pour but de déterminer l'approche qui offre le meilleur compromis entre précision et coût. L'extension de la méthode aux systèmes des lois de conservation, et en particulier à celui des équations d'Euler de la dynamique des gaz, représente le noyau principal de la thèse, et est abordée dans la troisième partie. Contrairement aux schémas de distribution classiques, qui reposent sur une extension formelle du cas scalaire, l'approche développée ici repose sur une décomposition du résidu par élément en équations scalaires, modélisant le transport de variables caracteristiques. La difficulté vient du fait que les équations d'Euler instationnaires ne se diagonalisent pas, et admettent une infinité de solutions élémentaires (ondes simples) se propageant dans toutes les directions d'espace. En régime stationnaire, en revanche, les équations se diagonalisent complètement dans le cas des écoulements supersoniques, et partiellement dans le cas des écoulements subsoniques. Ainsi, les équations sous forme conservative peuvent être remplacées par un système équivalent comprenant deux équations totalement découplées, exprimant l'invariance de l'entropie et de l'enthalpie totale le long des lignes de courant, et deux autres équations, modélisant les effets purement acoustiques. En régime supersonique, celles-ci se découplent aussi, et expriment la convection le long des lignes de Mach d'invariants de Riemann généralisés. La discrétisation de ces équations par des schémas scalaires décentrés permet de simuler des écoulements continus et discontinus avec une grande précision et sans oscillations. Finalement, dans une dernière partie, l'extension aux équations de Navier-Stokes est abordée, et la discrétisation des termes visqueux par une approche éléments finis est proposée. Les résultats numériques confirment la précision et la robustesse de la méthode. |
