Ce travail a porté principalement sur le développement et l'étude de méthodes numériques de type éléments finis multi-échelles pour un problème d'advection diffusion multi-échelles dominé par l'advection. Deux types d'approches sont envisagées: prendre en compte l'advection dans la construction de l'espace d'approximation, ou appliquer une méthode de stabilisation. On commence par l'étude d'un problème d'advection diffusion, dominé par l'advection, dans un milieu hétérogène. On poursuit sur des problèmes d'advection-diffusion, sous le régime où l'advection domine, posés dans un domaine perforé. On se focalise ici sur la condition aux bords de type Crouzeix Raviart pour la construction des éléments finis multi-échelles. On considère deux situations différentes selon la condition prescrite au bord des perforations: la condition de Dirichlet homogène ou la condition de Neumann homogène. Cette étude repose sur une hypothèse de coercivité.Pour finir, on se place dans un cadre général où l'opérateur d'advection-diffusion est non coercif, possiblement dominé par l'advection. On propose une approche éléments finis basée sur une mesure invariante associée à l'opérateur adjoint. Cette approche est bien posée inconditionnellement en la taille du maillage. On la compare numériquement à une méthode standard de stabilisation / This work essentially deals with the development and the study of multiscale finite element methods for multiscale advection-diffusion problems in the advection-dominated regime. Two types of approaches are investigated: Take into account the advection in the construction of the approximation space, or apply a stabilization method. We begin with advection-dominated advection-diffusion problems in heterogeneous media. We carry on with advection-dominated advection-diffusion problems posed in perforated domains.Here, we focus on the Crouzeix-Raviart type boundary condition for the construction of the multiscale finite elements. We consider two different situations depending on the condition prescribed on the boundary of the perforations: the homogeneous Dirichlet condition or the homogeneous Neumann condition. This study relies on a coercivity assumption.Lastly, we consider a general framework where the advection-diffusion operator is not coercive, possibly in the advection-dominated regime. We propose a Finite Element approach based on the use of an invariant measure associated to the adjoint operator. This approach is unconditionally well-posed in the mesh size. We compare it numerically to a standard stabilization method

Cette thèse est consacrée à la simulation d'écoulements turbulents, incompressibles ou à faible nombre de Mach pour des applications touchant à la sûreté nucléaire. En particulier, nous nous concentrons sur le développement et l'analyse mathématique de schémas numériques performants pour la méthode dite de Simulation des Grandes Echelles. Ces schémas sont basés sur des méthodes à pas fractionnaires de type correction de pression et des éléments finis non conformes de bas degré. Deux arguments semblent essentiels à la construction de tels schémas: le contrôle de l'énergie cinétique et la précision pour des écoulements à convection dominante. Concernant la discrétisation en temps, nous proposons un schéma de type Crank-Nicolson et nous montrons qu'il satisfait un contrôle de l'énergie cinétique. Ce schéma présente de plus l'avantage d'être peu dissipatif numériquement (résidu d'ordre deux en temps). Concernant le défaut de précision de la discrétisation par l'élément fini de Rannacher-Turek, nous envisageons deux approches. La première consiste à construire un schéma pénalisé contraignant les degrés de liberté tangents aux faces des cellules à s'écrire comme combinaison linéaire des degrés de liberté normaux alentour. La deuxième approche repose sur l'enrichissement de l'espace discret d'approximation pour la pression. Enfin, différents tests numériques sont présentés en dimensions deux et trois et pour des maillages généraux, afin d'illustrer les capacités des schémas étudiés et de confronter les résultats théoriques et expérimentaux.

Cette thèse est consacrée à la simulation d'écoulements turbulents, incompressibles ou à faible nombre de Mach pour des applications touchant à la sûreté nucléaire. En particulier, nous nous concentrons sur le développement et l'analyse mathématique de schémas numériques performants pour la méthode dite de Simulation des Grandes Echelles. Ces schémas sont basés sur des méthodes à pas fractionnaires de type correction de pression et des éléments finis non conformes de bas degré. Deux arguments semblent essentiels à la construction de tels schémas: le contrôle de l'énergie cinétique et la précision pour des écoulements à convection dominante. Concernant la discrétisation en temps, nous proposons un schéma de type Crank-Nicolson et nous montrons qu'il satisfait un contrôle de l'énergie cinétique. Ce schéma présente de plus l'avantage d'être peu dissipatif numériquement (résidu d'ordre deux en temps). Concernant le défaut de précision de la discrétisation par l'élément fini de Rannacher-Turek, nous envisageons deux approches. La première consiste à construire un schéma pénalisé contraignant les degrés de liberté tangents aux faces des cellules à s'écrire comme combinaison linéaire des degrés de liberté normaux alentour. La deuxième approche repose sur l'enrichissement de l'espace discret d'approximation pour la pression. Enfin, différents tests numériques sont présentés en dimensions deux et trois et pour des maillages généraux, afin d'illustrer les capacités des schémas étudiés et de confronter les résultats théoriques et expérimentaux. / This thesis is devoted to the simulation of incompressible or low Mach turbulent flows, for nuclear safety applications. In particular, we focus on the development and analysis of performing numerical schemes for the Large Eddy Simulation technique. These schemes are based on fractional step methods of pressure correction type and on nonconforming low degree finite elements. Two requirements seems essential to build such schemes, namely a control of kinetic energy and the accuracy for convection dominated flows. Concerning the time marching algorithm, we propose a Crank-Nicolson like scheme for which we prove a kinetic energy control. This scheme has the advantage to be numerically low dissipative (numerical dissipation residual is second order in time). Concerning the low accracy of the Rannacher-Turek discretization, two approaches are investigated in this work. The first one consists in building a penalized scheme constraining the velocity degrees of freedom tangent to the faces to be written as a linear combination of the normal ones. The second approach relies on the enrichment of the pressure approximation discrete space. Finally, various numerical tests are presented in both two and three dimensions and for general meshes, to illustrate the capacity of the schemes and compare theoretical and experimental results.

En mécanique des fluides numérique, un enjeu est le développement de méthodes d'approximation d'ordre élevé, comme celles de Galerkin Discontinues (GD). Si ces méthodes permettent d'envisager la simulation d'écoulements complexes en alternative aux méthodes usuelles d'ordre deux, elles souffrent cependant d'une forte restriction sur le pas de temps lorsqu'elles sont associées à une discrétisation explicite en temps. Ce travail de thèse consiste à mettre en œuvre une stratégie d'intégration temporelle explicite-implicite efficace, associée à une discrétisation spatiale GD d'ordre élevé, pour les écoulements instationnaires à convection dominante de fluides visqueux compressibles modélisés par le système des équations de Navier-Stokes. La discrétisation spatiale de la méthode GD est associée à des flux numériques de fluides parfaits et visqueux à stencil compact. En présence de frontières matérielles courbes, l'ordre élevé est garanti par la discrétisation du domaine de calcul à l'aide d'une représentation iso-paramétrique. La stratégie d'intégration temporelle repose sur une décomposition d'opérateurs de Strang, où les termes de convection sont résolus explicitement et ceux de diffusion implicitement. Son efficacité résulte d'une simplification du schéma implicite, où le calcul de la matrice implicite est approché avec une méthode sans jacobienne et où les degrés de liberté du schéma sont découplés. De fait, la taille du système linéaire à résoudre et le temps de calcul de la résolution sont significativement réduits. Enfin, la validation et l'évaluation des performances du schéma numérique sont réalisées à travers cinq cas tests bien référencés en deux dimensions d'espace.