Les modèles de turbulence de type Navier-Stokes en moyenne de Reynolds (RANS) au premier ordre sont étudiés dans cette thèse. Ils sont constitués des équations de Navier-Stokes, auxquelles on adjoint un système d'équations de bilan pour des échelles scalaires caractéristiques de la turbulence. L'évaluation de celles-ci permet, grâce à une relation algébrique, de calculer une viscosité additionnelle dite "turbulente", modélisant la contribution de l'agitation turbulente dans les équations de Navier-Stokes. Les problèmes d'analyse numérique abordés se placent dans le contexte d'un algorithme à pas fractionnaire constitué d'une approximation, sur un maillage régulier, des équations de Navier-Stokes par éléments finis non-conformes de Crouzeix-Raviart, ainsi que d'un ensemble d'équations de bilan de la turbulence de type convection-diffusion, discrétisées par la méthode de volumes finis standard. Un schéma numérique basé sur une discrétisation de volumes finis, permettant de préserver la positivité des échelles turbulentes telles que l'énergie cinétique turbulente (k) et son taux de dissipation (ε), est ainsi proposé dans le cas des modèles k − ε standard, k − ε RNG et leur extension k − ε − v2 − f. La convergence du schéma numérique proposé est ensuite étudiée sur un problème modèle constitué des équations de Stokes incompressibles et d'une équation de convection-diffusion stationnaires, couplées par les viscosités et le terme de production turbulente. Il permet d'aborder la difficulté principale de l'analyse d'un tel problème : l'expression du terme de production turbulente amène à considérer, pour les équations de bilan de la turbulence, un problème de convection-diffusion avec second membre appartenant à L1. Enfin, afin d'aborder le problème instationnaire, on montre la convergence du schéma de volumes finis pour une équation de convection-diffusion modèle avec second membre appartenant à L1. Les estimations a priori de la solution et de sa dérivée en temps sont obtenues dans des normes discrètes dont les espaces correspondants ne sont pas duaux. Un résultat de compacité plus général que le théorème de Kolmogorov usuel, qui se pose comme un équivalent discret du Lemme d'Aubin-Simon, est alors proposé et permet de conclure à la convergence dans L1 d'une suite de solutions discrètes.

Ce travail a pour objectif de développer et d'appliquer une méthode originale permettant de simuler l'écoulement dans un milieu poreux fracturé. Cette méthode repose sur une approche multicouches double continuum permettant de séparer le comportement des différents aquifères présents sur un site. La résolution des écoulements, basée sur la méthode des Eléments Finis de Crouzeix-Raviart, est associée à une méthode inverse (minimisation de type Quasi-Newton combinée à la méthode de l'état adjoint) et à une paramétrisation multi-échelle.La méthode est appliquée dans un premier temps sur l'aquifère fracturé du site expérimental de Poitiers. Les résultats montrent une bonne restitution du comportement de l'aquifère et aboutissent à des champs de transmissivité plus réguliers par rapport à ceux de l'approche simple continuum. L'application finale est réalisée sur le site de Cadarache (taille plus importante et données d'entrée moins denses). Le calage des deux aquifères présents sur le site est satisfaisant et montre que ceux-ci se comportent globalement de façon indépendante. Ce calage pourra être amélioré localement grâce à données de recharge plus fines.

On considère le problème d'advection-diffusion stationnaire v(∇u, ∇v)+( β*∇u, v) = (f, v) et non stationnaire d/dt (u(t), v) + v(∇u, ∇v)+( β*∇u, v) = (g(t), v), ainsi que le problème d'advection (β*∇u, v) = (f, v) sur un domaine polygonal borné du plan. Le terme de diffusion est approché par des éléments de Crouzeix Raviart et le terme de convection par une méthode upwind sur des volumes barycentriques finis avec un maillage triangulaire. Pour le problème stationnaire d'advection-diffusion, la L²-stabilité (c'est-à-dire indépendante du coefficient de diffusion v) est démontrée pour la solution du problème approché obtenue par cette méthode d'éléments finis et de volumes finis. Pour cela une condition sur la géométrie doit être satisfaite. Des exemples de maillages sont donnés. Toujours avec cette condition géométrique sur le maillage, une inégalité de stabilité (où la discrétisation en temps n'est pas couplée à une condition sur la finesse du maillage) est obtenue pour le cas non-stationnaire. La discrétisation en temps y est faite par un schéma d'Euler implicite. Une majoration de l'erreur, proportionnelle au pas en temps et à la finesse du maillage, est ensuite proposée et exprimée explicitement en fonction des données du problème. Pour le problème d'advection, une approche utilisant la théorie des graphes est utilisée pour obtenir l'existence et l'unicité de la solution, ainsi que le résultat de stabilité. Comme pour la stabilité du problème d'advection-diffusion, une condition géométrique - qui est équivalente pour les points intérieurs du maillage à celle du problème d'advection-diffusion - est nécessaire.

Ce travail a pour objectif de développer et d’appliquer une méthode originale permettant de simuler l’écoulement dans un milieu poreux fracturé. Cette méthode repose sur une approche multicouches double continuum permettant de séparer le comportement des différents aquifères présents sur un site. La résolution des écoulements, basée sur la méthode des Eléments Finis de Crouzeix-Raviart, est associée à une méthode inverse (minimisation de type Quasi-Newton combinée à la méthode de l’état adjoint) et à une paramétrisation multi-échelle.La méthode est appliquée dans un premier temps sur l’aquifère fracturé du site expérimental de Poitiers. Les résultats montrent une bonne restitution du comportement de l’aquifère et aboutissent à des champs de transmissivité plus réguliers par rapport à ceux de l’approche simple continuum. L’application finale est réalisée sur le site de Cadarache (taille plus importante et données d’entrée moins denses). Le calage des deux aquifères présents sur le site est satisfaisant et montre que ceux-ci se comportent globalement de façon indépendante. Ce calage pourra être amélioré localement grâce à données de recharge plus fines. / The aim of this study is to develop and validate a new method for the simulation of flow in fractured porous media. This method is based on a multi-layered and dual continuum approach allowing to discriminate the behavior of different aquifers present on a site. The flow equations are solved using a Crouzeix-Raviart Finite Element method, in association with an inverse method (Quasi-Newton minimization combined with the adjoint state method) and a multi-scale parameterization.The method is first applied and validated on the fractured aquifer of the Hydrogeological Experimental Site of Poitiers. The results closely reproduce the flow behavior of the aquifer and lead to a transmissivity field much more homogeneous than the one obtained with a simple continuum approach. The final application is performed on the site of Cadarache (large scale problem with heterogeneously distributed input data). The model calibration of both aquifers is rather satisfactory and shows that their behavior is globally independent. It could locally be improved if more accurate groundwater recharge data is made available.

On considère le problème d’advection-diffusion stationnaire v(∇u, ∇v)+( β•∇u, v) = (f, v) et non stationnaire d/dt (u(t), v) + v(∇u, ∇v)+( β•∇u, v) = (g(t), v), ainsi que le problème d’advection (β•∇u, v) = (f, v) sur un domaine polygonal borné du plan. Le terme de diffusion est approché par des éléments de Crouzeix Raviart et le terme de convection par une méthode upwind sur des volumes barycentriques finis avec un maillage triangulaire. Pour le problème stationnaire d’advection-diffusion, la L²-stabilité (c’est-à-dire indépendante du coefficient de diffusion v) est démontrée pour la solution du problème approché obtenue par cette méthode d’éléments finis et de volumes finis. Pour cela une condition sur la géométrie doit être satisfaite. Des exemples de maillages sont donnés. Toujours avec cette condition géométrique sur le maillage, une inégalité de stabilité (où la discrétisation en temps n’est pas couplée à une condition sur la finesse du maillage) est obtenue pour le cas non-stationnaire. La discrétisation en temps y est faite par un schéma d’Euler implicite. Une majoration de l’erreur, proportionnelle au pas en temps et à la finesse du maillage, est ensuite proposée et exprimée explicitement en fonction des données du problème. Pour le problème d’advection, une approche utilisant la théorie des graphes est utilisée pour obtenir l’existence et l’unicité de la solution, ainsi que le résultat de stabilité. Comme pour la stabilité du problème d’advection-diffusion, une condition géométrique - qui est équivalente pour les points intérieurs du maillage à celle du problème d’advection-diffusion - est nécessaire. / We consider the stationary linear convection-diffusion equation v(∇u, ∇v)+( β•∇u, v) = (f, v), the time dependent d/dt (u(t), v) + v(∇u,∇v)+( β•∇u, v)= (g(t), v) equation and the linear advection equation (β•∇u, v) = (f, v) on a two dimensional bounded polygonal domain. The diffusion term is discretized by Crouzeix-Raviart piecewise linear finite elements, and the convection term by upwind barycentric finite volumes on a triangular grid. For the stationary convection-diffusion problem, L²-stability (i.e. independent of the diffusion coefficient v) is proven for the approximate solution obtained by this combined finite-element finite-volume method. This result holds if the underlying grid satisfies a condition that is fulfilled, for example, by some structured meshes. Using again this condition on the grid, stability is shown for the time dependent convection-diffusion equation (without any link between mesh size and time step). An implicit Euler approach is used for the time discretization. It is shown that the error associated with this scheme decays linearly with the mesh size and the time step. This result holds without any link between mesh size and time step. The dependence of the corresponding error bound on the diffusion coefficient is completely explicit. For the stationary advection equation, an approach using graph theory is used to obtain existence, uniqueness and stability. As in the stationary linear convection-diffusion equation, the underlying grid must satisfy some geometric condition. / Gegenstand der Arbeit ist die zweidimensionale stationäre Konvektion-Diffusionsgleichung v(∇u, ∇v)+( β•∇u, v) = (f, v), die zeitabhängige Konvektion-Diffusionsgleichung d/dt (u(t), v) + v(∇u,∇v)+( β•∇u, v)= (g(t), v), sowie die Konvektionsgleichung (β•∇u, v) = (f, v). Der Diffusionsterm ist diskretisiert mittels Crouzeix-Raviart stückweise lineare Finite Elemente. Das Gebiet ist in Dreiecke unterteilt und der Konvektionsterm ist mittels einer upwind Methode auf Baryzentrische Finite Volumenelemente definiert. Für die stationäre Konvektion-Diffusionsgleichung, wird (d.h. von v unabhängige) L²-Stabilität der numerischen Lösung bewiesen. Voraussetzung dafür, ist die Erfüllung gewisser geometrischer Bedingungen an die Unterteilung des Gebiets. Beispiele von Unterteilungen die diese Bedingungen erfüllen, werden gegeben. Wieder an dieser geometrischen Bedingung geknüpft, wird Stabilität (d.h. die Zeitdiskretisierung ist entkoppelt von der Netzweite) für die zeitabhängige Konvektion-Diffusionsgleichung, bewiesen. Für die Zeitableitung wird dabei eine Implizite Euler Diskretisierung verwendet. Eine obere Schranke für den Diskretisierungsfehler, proportional zum Zeitdiskretisierungsparameter und zur Netzfeinheit, ausgedrückt als Funktion der Daten der Differenzialgleichung, wird gezeigt. Für die Konvektionsgleichung wird ein graphentheoretischer Zugang verwendet, der es ermöglicht Existenz, Eindeutigkeit und Stabilität, zu bekommen. Für die Stabilität, werden ähnliche geometrische Bedingungen an die Unterteilung des Gebiets gestellt, wie beim stationären Konvektion-Diffusionsproblem.