Nous nous intéressons aux méthodes permettant d'approcher les solutions de systèmes d'équations aux dérivées partielles conservatives. Dans les cas où l'écoulement est très<br />complexe - lorsqu'il y a plusieurs modèles physiques à calculer sur des zones difficiles à délimiter, on utilise des méthodes de couplage par recouvrement de domaine. <br />Nous présentons ici un algorithme, nouveau et performant, calculé grâce à une superposition de deux maillages correspondant à deux schémas différents. On utilise des projections conservatives de la solution d'un maillage vers l'autre.<br />Cette méthode de décomposition de domaine ne fait<br />pas intervenir de conditions aux limites artificielles. Elle est basée sur une régularisation de la fonction de Heaviside sur la zone de couplage. Elle est parfaitement conservative et donc bien indiquée pour l'étude des lois de conservation. <br />L'analyse mathématique est réalisée pour les problèmes hyperboliques, dans le cas scalaire multidimensionnel. Elle est basée sur le convergence des schémas volumes finis. Tout d'abord, on obtient la convergence de la solution mesure grâce aux travaux de Diperna, puis on estime l'erreur de convergence en $h^(^1/_4)$. Une nouvelle estimation de type $H^1$ faible permet d'estimer les erreurs induites par le couplage.<br />De nombreuses applications numériques en mécanique des fluides avec les tubes à chocs et de détente montrent que la méthode est très stable et conservative. Nous utilisons aussi la méthode sans grille appelée Smooth Particule Hydrodynamics - plus précisément sa nouvelle variante renormalisée - pour calculer la création d'un jet en couplant la méthode volumes finis à la méthode SPH. On montre ainsi la robustesse de l'algorithme de couplage et sa souplesse pour le calcul des écoulement complexes. <br />Cette étude à fait l'objet d'une collaboration avec l'équipe du Pr. D. Kröner de l'Institut des Mathématiques Appliquées à l'Université de Frieburg (Allemagne).

Cette thèse porte sur le couplage chimie-transport pour la modélisation et la simulation du stockage géologique de CO2. Nous présentons un modèle d'écoulement multiphasique et un modèle géochimique permettant de décrire un modèle couplé d'écoulement multiphasique réactif. Nous proposons ensuite deux méthodes de résolution, l'une est une méthode globale, l'autre est une méthode de splitting utilisée à l'IFP dans le logiciel COORES. Le splitting effectué pour cette méthode repose sur des hypothèses physiques. La méthode de couplage utilisée est une méthode de couplage non itérative dans laquelle l'erreur de splitting est corrigée à l'aide d'un terme de pénalisation. Une étude de convergence sur un cas simplifié permet de montrer que le schéma pénalisé est convergent vers la même solution que le schéma global. Une partie de cette thèse est consacrée à l'étude des phénomènes de diffusion-dispersion. On s'intéresse particulièrement à ce terme car il ne peut être intégré facilement dans un schéma de splitting si l'on souhaite résoudre le modèle de transport réactif de façon locale (nécessaire pour l'utilisation de sous-pas de temps locaux). Après avoir mis en évidence l'importance de ce terme sur un cas test représentatif, nous montrons la difficulté de l'intégrer dans le schéma de splitting. Finalement, on étudie un problème d'écoulement miscible en 1D d'un point de vue mathématique. Les difficultés proviennent de la non linéarité due à la solubilité non nulle du gaz dans l'eau. Nous proposons une définition d'une solution faible pour ce problème dont l'existence est montrée à l'aide de la convergence d'un schéma volumes finis de type Godunov.

Cette thèse concerne l'étude des interactions entre le plasma et la paroi d'un réacteur à fusion nucléaire de type tokamak. L'objectif est de proposer des méthodes de résolution des systèmes d'équations issus de modèles de plasma de bord. Nous nous sommes intéressés au traitement de deux difficultés qui apparaissent lors de la résolution numérique de ces modèles. La première difficulté est liée à la forme complexe de la paroi du tokamak. Pour cela, il a été choisi d'utiliser des méthodes de pénalisation volumique. Des tests numériques de plusieurs méthodes de pénalisation ont été réalisés sur un problème hyperbolique non linéaire avec un domaine 1D. Une de ces méthodes a été étendue à un système hyperbolique quasilinéaire avec bord non caractéristique et conditions aux limites maximales strictement dissipatives sur un domaine multidimensionnel : il est alors démontré que cette méthode de pénalisation ne génère pas de couche limite. La deuxième difficulté provient de la forte anisotropie du plasma, entre la direction parallèle aux lignes de champ magnétique et la direction radiale. Pour le potentiel électrique, cela se traduit par une résistivité parallèle très faible. Afin d'éviter les difficultés liées au fait que le problème devient mal posé quand la résistivité parallèle tend vers 0, nous avons utilisé des méthodes de type asymptotic-preserving (AP). Pour les problèmes non linéaires modélisant le potentiel électrique avec un domaine 1D et 2D, nous avons fait l'analyse théorique ainsi que des tests numériques pour deux méthodes AP. Des tests numériques sur le cas 1D ont permis une étude préliminaire du couplage entre les méthodes de pénalisation volumique et AP. / This thesis deals with the study of wall plasma interactions in a nuclear fusion reactor such as a tokamak. The goal is to propose methods to solve partial differential equations issued from edge plasma models. We focus on two difficulties for the numerical resolution of these models. The first issue concerns the complex shape of the tokamak wall: we choose volume penalty methods. Numerical tests on several penalization methods have been performed on a nonlinear hyperbolic problem. One of these methods has been extended to a quasilinear hyperbolic system with a non characteristic boundary and maximally strictly dissipative boundary conditions on a multidimensional domain: it is proven that this penalty method does not generate any boundary layer. The second question comes from the strong plasma anisotropy between the direction parallel to the magnetic field lines and the radial one. Concerning the electrical potential, this results in a very low parallel resistivity. In order to avoid the troubles due to the ill-posedness of the equations when the parallel resistivity tends to 0, we study asymptotic preserving (AP) methods. For 1D and 2D nonlinear models of the electrical potential, we performed the theoretical analysis and numerical simulations for two AP methods. A preliminary study of the coupling between volume penalty and AP methods is also presented.