Cette thèse présente la construction de schémas distribuant le résidu (RD) d'ordre très élevés, pour la discrétisation d'équations d'advection-diffusion multidimensionnelles et stationnaires sur maillages non structurés. Des schémas linéaires ainsi que des schémas non linéaires sont considérés. Une approximation de la solution polynomiale par morceaux et continue sur chaque élément est adoptée, de plus une procédure de reconstruction du gradient que celle de la solution numérique est utilisée afin d'avoir une représentation continue de la solution numérique et de son gradient. Il est montré que le gradient doit être reconstruit avec la même précision de la solution, sans quoi la précision formel du schéma numérique est perdue dans les cas où les effets de diffusion prévalent sur les effets d'advection, et aussi quand l'advection et la diffusion sont également importants. Ensuite, la méthode est étendue à des systèmes d'équations, en particulier aux équations de Navier-Stokes et aux équations RANS. La précision, l'efficacité et la robustesse du solveur RD implicite sont démontrées sur plusieurs cas tests.

Le bon fonctionnement et les performances des turbines à vapeur sont liés à l’état de la vapeur et notamment au taux d’humidité qu’elle contient. EDF souhaite pouvoir maîtriser les phénomènes spécifiques à ces problématiques afin d’améliorer l’utilisation et l’évolution de ses turbines. Le sujet de recherche concerne la modélisation de la formation de l’humidité dans un corps de turbine et l’étude des couplages entre la phase liquide et les instationnarités. Dans ce contexte, la démarche adoptée est la suivante : la présence d’humidité est prise en compte à l’aide d’un modèle homogène, couplé à des modèles de condensation permettant de prendre en compte les phénomènes hors-équilibre thermodynamique : le grossissement et la nucléation des gouttes d’eau dans la vapeur. Pour mener à bien les calculs, des méthodes numériques adaptées aux gaz réels ont été utilisées et testées à l’aide d’un code monodimensionnel avant d’être intégrées dans le code 3D elsA. Deux types de modèles de condensation ont été mis en œuvre, considérant ou non la polydispersion des gouttes dans la vapeur. Les couplages instationnaires entre la condensation et l’écoulement principal ont été étudiés à différents niveaux d’observations (1D, 1D − 3D, 3D). Il a été montré que la méthode des moments apporte une richesse supplémentaire par rapport à un modèle mono-dispersé, et permet de mieux capter les couplages instationnaires entre l’humidité et le champ principal. / In addition to conventional turbomachinery problems, both the behavior and performances of steam turbines are highly dependent on the vapour thermodynamic state and the presence of a liquid phase. EDF, the main French electricity producer, is interested in further developing its’ modelling capabilities and expertise in this area to allow for operational studies and long-term planning. This PhD thesis explores the modelling of wetness formation and growth in a steam turbine and an analysis of the coupling between the liquid phase and the main flow unsteadiness. To this end, the work in this thesis took the following approach. Wetness was accounted for using a homogeneous model coupled with transport equations to take into account the effects of non-equilibrium phenomena, such as the growth of the liquid phase and nucleation. The real gas attributes of the problem demanded adapted numerical methods. Before their implementation in the 3D elsA solver, the accuracy of the chosen models was tested using a developed one-dimensional nozzle code. In this manner, various condensation models were considered, including both polydispersed and monodispersed behaviours of the steam. Finally, unsteady coupling effects were observed from several perspectives (1D, 1D − 3D, 3D), demonstrating the ability of the method of moments to sustain unsteady phenomena which were not apparent in a simple monodispersed model.

Cette thèse contribue à la simulation numérique des écoulements d'un fluide compressible modélisé par les équations de Euler et Navier-Stokes: étude d'un schéma d'ordre élévé basé sur une matrice de masse, modélisation des écoulement turbulents compressibles à très haut Reynolds, développement des estimateurs et correcteurs d'erreurs a posteriori et a priori, et adaptation de maillage anisotrope pour les fonctionnelles d'observation. Concernant la prédiction des écoulements turbulents, on s'est intéressé aux modèles hybrides de type RANS/LES comportant les nouveautés suivantes: traitement des tourbillons de grande échelle utilisant la formulation VMS (Variational Multi-Scale) et du RANS employé sur la paroi sur une distance imposée via une zone de protection conçue afin d'éviter le phénomène assez commun apelé "grid induced model depletion". Le niveau de viscosité du modèle VMS-LES est de plus controllé par un procédé de double filtre dynamique. La seconde partie concerne l'adaptation de maillage anisotrope pour mieux observer une fontionnelle d'observation. Les estimations a priori sont réalises pour le modèle des équations d'Euler et Navier-Stokes en instationnaire en 2D et 3D. A partir de ces estimations on sait définir les maillages optimaux au cours du calcul instationnaire, en fonction de l'état et de l'état adjoint. Le système d'optimalité est discrétisé et résolu à l'aide d'une méthode de point fixe instationnaire global, comportant une stratégie de stockage/recalcul pour le couplage état/ état adjoint. Des applications à la propagation d'ondes de choc et d'ondes acoustiques sont présentées.

Cette thèse présente la construction d’un schéma aux Résidus Distribués p-adaptatif pour la discrétisation des équations d’Euler ainsi qu’un schéma aux Résidus Distribués hp-adaptatif pour les équations de Navier- Stokes pénalisées. On rappelle tout d’abord les équations d’Euler et de Navier-Stokes ainsi que leurs versions non dimensionnelles. Les définitions et propriétés de base des schémas aux Résidus Distribués sont ensuite présentées. On décrit alors la construction d’un schéma aux Résidus Distribués p-adaptatif pour les équations d’Euler. La construction du schéma p-adaptatif est basée sur la possibilité d’exprimer le résidu total d’un élément K de degré k (au sens où l’élément fini (K; P; Sigma ) est un élément fini de degré k) comme une somme pondérée des résidus totaux de ses sous-éléments de degré 1. La solution discrète ainsi obtenue est en général discontinue à l’interface entre un élément subdivisé et un élément non subdivisé. Ceci contredit l’hypothèse de continuité de la solution qui est utilisée pour démontrer le théorème de Lax-Wendroff discret pour les schémas aux Résidus Distribués. Cependant, on montre que cette hypothèse peut être assouplie. La conséquence pratique est que si l’on emploie des quadratures particulières dans l’implémentation numérique, on peut quand même démontrer le théorème de Lax-Wendroff discret, ce qui garantit la convergence du schéma numérique vers une solution faible des équations d’origine. Les formules qui permettent d’exprimer le résidu total comme une somme pondérée des résidus totaux des sous-éléments sont à la base de la méthode de p-adaptation présentée ici. Dans le cas quadratique, la formule est obtenue avec les classiques fonctions de base de Lagrange en dimension deux et avec des fonctions de base de Bézier en dimension trois. Ces deux formules sont ensuite généralisées à des degrés polynomiaux quelconques en dimension deux et trois avec des fonctions de base de Bézier. Dans la deuxième partie de la thèse, on présente l’application du schéma p-adaptatif aux équations pénalisées de Navier-Stokes avec adaptation de maillage anisotrope. . En pratique, on combine le schéma p-adaptatif avec la méthode IBM-LS-AUM (Immersed Boundary Method with Level Sets and Adapted Unstructured Meshes). La méthode IBM-LS-AUM permet d’imposer les conditions aux bords grâce à la méthode de pénalisation et l’adaptation anisotrope du maillage à la solution numérique et à la level-set augmente la précision de la solution et de la représentation de la surface. Une fois la méthode IBM-LS-AUM combinée avec le schéma p-adaptatif, il est alors possible d’utiliser des éléments d’ordre élevés en-dehors de la zone où la pénalisation est appliquée. La méthode est robuste comme le montrent les diverses expérimentations numériques à des vitesses faibles à élevées et à différents nombres de Reynolds. / This thesis presents the construction of a p-adaptive Residual Distribution scheme for the steady Euler equations and a hp-adaptive Residual Distribution scheme for the steady penalized Navier-Stokes equations in dimension two and three. The Euler and Navier-Stokes equations are recalled along with their non dimensional versions. The basis definitions and properties of the steady Residual Distribution schemes are presented. Then, the construction of a p-adaptive Residual Distribution scheme for the Euler equations is considered. The construction of the p-adaptive scheme is based upon the expression of the total residual of an element of a given degree k (in the Finite Element sense) into the total residuals of its linear sub-elements. The discrete solution obtained with the p-adaptive scheme is then a one degree polynomial in the divided elements and a k-th degree polynomial in the undivided ones. Therefore, the discrete solution is in general discontinuous at the interface between a divided element and an undivided one. This is in apparent contradiction with the continuity assumption used in general to demonstrate the discrete Lax-Wendroff theorem for Residual Distribution schemes. However, as we show in this work, this constrain can be relaxed. The consequence is that if special quadrature formulas are employed in the numerical implementation, the discrete Lax-Wendroff theorem can still be proved, which guaranties the convergence of the p-adaptive scheme to a weak solution of the governing equations. The formulas that express the total residual into the combination of the total residuals of the sub-elements are central to the method. In dimension two, the formula is obtained with the classical Lagrange basis in the quadratic case and with the Bézier basis in dimension three. These two formulas are then generalized to arbitrary polynomial degrees in dimension two and three with a Bézier basis. In the second part of the thesis the application of the p-adaptive scheme to the penalized Navier-Stokes equations with anisotropic mesh adaptation is presented. In practice, the p-adaptive scheme is used with the IBM-LS-AUM (Immersed Boundary Method with Level Sets and Adapted Unstructured Meshes) method. The IBM-LS-AUM allows to impose the boundary conditions with the penalization method and the mesh adaptation to the solution and to the level-set increases the accuracy of the representation of the surface and the solution around walls. When the IBM-LSAUM is combined with the p-adaptive scheme, it is possible to use high-order elements outside the zone where the penalization is applied. The method is robust as shown by the numerical applications at low to large Mach numbers and at different Reynolds in dimension two and three.