On s'intéresse à l'écoulement d'un mélange d'eau et d'huile dans une matrice poreuse supposée hétérogène, et plus particulièrement apposition de différentes sous-matrices poreuses supposées homogènes. Si la modélisation et l'analyse des écoulements diphasiques dans des milieux poreux homogènes a fait l'objet de nombreuses études préalables, ce travail s'intéresse aux phénomènes liés aux forces provenant de la pression capillaire au niveau des interfaces entre des milieux différents.<br />Dans un premier temps, on suppose que l'on peut connecter les pressions au niveau des interfaces. Cela nécessite des hypothèses sur les profils de pression capillaire, afin que les raccords soient possibles. On démontre l'existence d'une solution faible du problème parabolique dégénéré obtenu par convergence d'une famille de solutions approchées obtenues à l'aide d'un schéma Volumes Finis. L'unicité est garantie, sous hypothèse sur les dégénérescence, par une méthode de dédoublement de variable aboutissant à un principe de contraction $L^1$.<br />La modélisation ne garantit pas forcément que le raccord des pressions capillaires aux interfaces soit possible. Dans le chapitre 3, on donne une condition de raccord graphique des pressions capillaires aux interfaces qui permet de traiter des cas beaucoup plus généraux. On montre que de le problème avec raccords graphiques admet une solution. Un résultat d'unicité et de contraction $L^1$ est donné dans le cas unidimensionnel.<br />Dans le chapitre 4, on montre la convergence d'une approximation Volumes Finis vers l'unique solution du problème unidimensionnel. Ce résultat utilise une borne uniforme sur les flux discrets, analogie discrète de la preuve dans le cas continue faite au chapitre précédent.<br />On étudie dans les chapitres 5 et 6 la limite des solutions lorsque la dépendance de la pression capillaire par rapport à l'inconnue saturation devient très faible, et que la pression capillaire ne dépend plus que du sous milieux poreux homogène. Il apparaît alors des phénomènes différents selon l'orientation des forces de gravité et de capillarité. Soit la solution su problème est la solution entropique d'une équation hyperbolique à flux discontinus, soit une solution faible, entropique à l'intérieur des sous-domaines homogènes, et laissant apparaître un choc non classique à l'interface.

Le circuit primaire d'une centrale nucléaire est composé d'un ensemble d'éléments très différents (cuve, coeur, réseau de conduite ...). A chacun de ces éléments correspond actuellement un ou plusieurs codes de calcul spécifiques basés sur des systèmes d'équations aux dérivées partielles spécifiques. Afin de permettre la simulation des écoulements diphasiques dans l'ensemble du circuit primaire, il faut envisager de coupler ces différents codes. L'approche proposée dans ce travail de thèse est de coupler les codes grâce à un échange d'information interfaciale instationnaire. Des flux numériques sont calculés au niveau des interfaces de couplage et servent de conditions aux limites à chacun des codes. Les méthodes permettant le calcul des flux de couplage sont dérivées du formalisme de Greenberg-Leroux proposé dans le cadre du décentrement des termes sources des systèmes hyperboliques non-homogènes stationnaires, et font intervenir un modèle d'interface. Trois cas de couplage ont été examinés : (i) le couplage du système des équations d'Euler en dimension un et deux ; (ii) le couplage de deux modèles diphasiques homogènes distincts ; (iii) le couplage d'un modèle homogène à quatre équations et du modèle bi-fluide standard.

Cette thèse de mathématiques appliquées traite de problèmes théoriques, numériques et asymptotiques en transport motivés par la physiologie rénale. Plus précisément, elle vise à comprendre et quantifier les échanges de solutés qui peuvent mener dans des cas pathologiques à des néphrocalcinoses, qui se caractérisent par des dépôts calciques dans le parenchyme rénal. Le manuscrit est constitué de deux parties. La première partie concerne le développement et l'analyse mathématique d'un modèle simplifié du rein. Il s'agit d'un système de 3 EDP hyperboliques à vitesses constantes, couplées par leur terme source non linéaire et assorti de conditions aux bords spécifiques. Le modèle rentre dans le cadre des modèles cinétiques avec un nombre fini de vitesses et des conditions aux bords de type réflexion. Nous montrons que ce système est bien posé, qu'il tend en temps grand vers un état stationnaire. On montre que le taux de convergence est exponentiel avec des éléments spectraux. Nous proposons l'étude du rôle de deux paramètres à travers une analyse asymptotique. L'une d'entre elles nous place dans le cadre de la relaxation hyperbolique vers une loi de conservation scalaire avec un flux hétérogène en espace sur un domaine borné. La deuxième partie concerne le développement et l'analyse numérique d'un modèle réaliste du rein. Il s'agit d'un système de 27 équations aux dérivées partielles de type hyperboliques dont les vitesses sont les solutions de 8 équations différentielles non linéaires, et toutes ces équations sont couplées par leur terme source. Les conditions aux bords sont là aussi spécifiques au modèle. Nous interprétons ensuite les résultats obtenus d'un point de vue physiologique, en proposant des prédictions de profils de concentration calciques dans le rein, dans le cas normal et dans certains cas pathologiques.

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.

Le but de cette thèse est de développer une méthode de volumes finis qui s'applique à une classe de maillages beaucoup plus grande que celle des méthodes classiques, limitées par des conditions d'orthogonalité très restrictives. On construit des opérateurs différentiels discrets agissant sur les trois maillages décalés nécessaires à la construction de la méthode. Ces opérateurs vérifient des propriétés discrètes analogues à celles des opérateurs continus. La méthode est tout d'abord appliquée au problème divergence-rotationnel qui peut etre considéré comme une brique du problème de Stokes. Ensuite, le problème de Stokes est discrétisé avec diverses conditions aux limites. Par ailleurs, il est bien connu que lorsque le domaine est polygonal et non-convexe, l'ordre de convergence des méthodes numériques se dégrade. Par conséquent, nous avons étudié sous quelles conditions un raffinement local approprié permet de restaurer l'ordre de convergence optimal. Enfin, nous avons discrétisé le problème non-linéaire de Navier-Stokes, en utilisant la formulation rotationnelle du terme de convection, associée à la pression de Bernoulli. Par un algorithme itératif, nous sommes amenés à résoudre un problème de point-selle à chaque itération, pour lequel nous testons quelques préconditionneurs issus des éléments finis, que l'on adapte (quand c'est possible) à la méthode. Chaque problème est illustré par des cas tests numériques sur des maillages "arbitraires", tels que des maillages fortement non-conformes.

L'étude des interactions entre un plasma et un champ magnétique joue un rôle important dans différents domaines tels que la fusion thermonucléaire par confinement magnétique, les plasmas astrophysiques. En évolution, ces interactions sont décrites par les équations de la magnétohydrodynamique (MHD). A l'équilibre, les équations de la MHD se réduisent à celles de la magnétostatique.<br />Les équations de la magnétostatique forment un système d'équations aux dérivées partielles non linéaires en dimension 3 faisant intervenir le champ magnétique et la pression cinétique du plasma. Quand on néglige la pression, le champ magnétique est alors dit de Beltrami. Nous proposons de résoudre numériquement les équations régissant les champs de Beltrami par un algorithme itératif de type point fixe associé à des méthodes d'éléments finis. Cette stratégie itérative est étendue au cas des configurations d'équilibres avec pression.<br />On s'intéresse ensuite à l'approximation des équations de la MHD idéale instationnaires. Il s'agit d'un système de loi de conservation hyperbolique non linéaire. Nous proposons une approche de type volumes finis dans laquelle les flux sont calculés par une méthode de Roe sur un maillage tétraédrique et où les flux du champ magnétique sont modifiés afin de satisfaire la contrainte de divergence nulle qui lui est imposée. <br />Les méthodes proposées ont été implantées dans deux nouveaux codes tridimensionnels TETRAFFF pour les équilibres, et TETRAMHD pour la MHD. Les résulats numériques obtenus par ces codes montrent la performance des méthodes employées.

Reproduction de : Thèse doctorat : Mathématiques appliquées : Pau : 2008.

Reproduction de : Thèse de doctorat : Mathématiques appliquées : Lille 1 : 2005. / N° d'ordre (Lille 1) : 3608. Résumé. Titre provenant de la page de titre du document numérisé. Bibliogr. p. 245-250.

La présente thèse comporte trois parties indépendantes.<br>La première partie présente une preuve d'existence de solutions faibles globales pour un modèle de sprays de type Vlasov-Navier-Stokes-incompressible avec densité variable. Ce modèle est obtenu par une limite formelle à partir d'un modèle Vlasov-Navier-Stokes-incompressible avec fragmentation, où seules deux valeurs de rayons de particules sont considérées : un rayon r1 pour les particules avant fragmentation, et un rayon r2<

Cette thèse concerne l'étude mathématique et numérique d'équations aux dérivées partielles hyperboliques non-linéaires. Une première partie traite d'une problématique émergente: le couplage d'équations hyperboliques. Les applications poursuivies relèvent du couplage mathématique de plateformes de calcul, en vue d'une simulation adaptative de phénomènes multi-échelles. Nous proposons et analysons un nouveau formalisme de couplage construit sur des systèmes EDP augmentés permettant de s'affranchir de la description géométrique des frontières. Ce nouveau formalisme permet de poser le problème en plusieurs variables d'espace en autorisant l'éventuel recouvrement des modèles à coupler. Ce formalisme autorise notamment à munir la procédure de couplage de mécanismes de régularisation visqueuse utiles à la sélection de solutions discontinues naturelles. Nous analysons alors les questions d'existence et d'unicité dans le cadre d'une régularisation parabolique autosemblable. L'existence est acquise sous des conditions très générales mais de multiples solutions sont susceptibles d'apparaître dès que le phénomène de résonance survient. Ensuite, nous montrons que notre formalisme de couplage à l'aide de modèles EDP augmentés autorise une autre stratégie de régularisation basée sur l'épaississement des interfaces. Nous établissons dans ce cadre l'existence et l'unicité des solutions au problème de Cauchy pour des données initiales $L^\infty$. À cette fin, nous développons une technique de volumes finis sur des triangulations générales que nous analysons dans la classe des solutions à valeurs mesures entropiques de DiPerna. La seconde partie est consacrée à la définition d'un schéma de volumes finis pour l'approximation des solutions non classiques d'une loi de conservation scalaire basée sur une relation cinétique. Ce schéma présente la particularité d'être stricto sensu conservatif contrairement à une approche à la Glimm qui ne l'est que statistiquement. Des illustrations numériques étayent le bien-fondé de notre approche.