Cette thèse est dédiée au développement et à l'analyse de schémas numériques de type volumes finis pour des équations de convection-diffusion, qui apparaissent notamment dans des modèles issus de la physique ou de la biologie. Nous nous intéressons plus particulièrement à la préservation de comportements asymptotiques au niveau discret. Ce travail s'articule en trois parties, composées chacune de deux chapitres. Dans la première partie, nous considérons la discrétisation du système de dérive-diffusion linéaire pour les semi-conducteurs par le schéma de Scharfetter-Gummel implicite en temps. Nous nous intéressons à la préservation par ce schéma de deux types d'asymptotiques : l'asymptotique en temps long et la limite quasi-neutre. Nous démontrons des estimations d'énergie--dissipation d'énergie discrètes qui permettent de prouver d'une part la convergence en temps long de la solution approchée vers une approximation de l'équilibre thermique, d'autre part la stabilité à la limite quasi-neutre du schéma. Dans la deuxième partie, nous nous intéressons à des schémas volumes finis préservant l'asymptotique en temps long dans un cadre plus général. Plus précisément, nous considérons des équations de type convection-diffusion non linéaires qui apparaissent dans plusieurs contextes physiques : équations des milieux poreux, système de dérive-diffusion pour les semi-conducteurs... Nous proposons deux discrétisations en espace permettant de préserver le comportement en temps long des solutions approchées. Dans un premier temps, nous étendons la définition du flux de Scharfetter-Gummel pour une diffusion non linéaire. Ce schéma fournit des résultats numériques satisfaisants si la diffusion ne dégénère pas. Dans un second temps, nous proposons une discrétisation dans laquelle nous prenons en compte ensemble les termes de convection et de diffusion, en réécrivant le flux sous la forme d'un flux d'advection. Le flux numérique est défini de telle sorte que les états d'équilibre soient préservés, et nous utilisons une méthode de limiteurs de pente pour obtenir un schéma précis à l'ordre deux en espace, même dans le cas dégénéré. Enfin, la troisième et dernière partie est consacrée à l'étude d'un schéma numérique pour un modèle de chimiotactisme avec diffusion croisée pour lequel les solutions n'explosent pas en temps fini, quelles que soient les données initiales. L'étude de la convergence du schéma repose sur une estimation d'entropie discrète nécessitant l'utilisation de versions discrètes d'inégalités fonctionnelles telles que les inégalités de Poincaré-Sobolev et de Gagliardo-Nirenberg-Sobolev. La démonstration de ces inégalités fait l'objet d'un chapitre indépendant dans lequel nous proposons leur étude dans un contexte assez général, incluant notamment le cas de conditions aux limites mixtes et une généralisation au cadre des schémas DDFV.

Ce travail de thèse porte sur l'extension et l'analyse d'un solveur volumes finis eulérien, co-localisé avec capture d'interfaces pour la simulation des écoulements multi-matériaux non miscibles. Les extensions proposées s'inscrivent dans la volonté d'élaborer un outil de simulation multi-physiques. Dans le cadre de ce mémoire, le caractère multi-physiques recouvre les champs que nous allons détailler. Nous traitons le cas des écoulements radiatifs modélisés par un système à deux températures qui couple les phénomènes purement hydrodynamiques aux phénomènes radiatifs. Nous proposons un solveur permettant la prise en compte des effets de tension superficielle à l'interface entre deux fluides. Nous développons un solveur implicite permettant la simulation précise d'écoulements faisant intervenir de faibles nombres de Mach par le biais d'une méthode de renormalisation de la diffusion numérique. Enfin, les effets tri-dimensionnels sont considérés ainsi que la possibilité d'étendre le schéma de base aux écoulements à un nombre quelconque de matériaux. A chaque étape, les solveurs développés sont validés sur des cas-tests.

Cette thèse est consacrée à l'étude de méthodes numériques pour la simulation des écoulements polyphasiques en milieu poreux, en vue de leur application à des problèmes d'ingénierie pétrolière ou environnementale. Nous présentons une formulation générique du modèle d'écoulements à nombre quelconque de composants présents dans un nombre quelconque de phases. Dans notre approche l'approximation des flux diffusifs (issus, par exemple, de la loi de Darcy) s'appuie sur de nouveaux schémas, appelés schémas gradient, qui ont plusieurs avantages sur les schémas industriels standard : ces derniers, qui sont des schémas volumes finis multi-points centrés aux mailles, ne sont généralement pas symétriques et convergent difficilement sur des cas à forts rapports d'anisotropie. Nous montrons en revanche que les schémas gradient conduisent naturellement à des approximations symétriques et convergentes. Parmi cette classe de schémas, nous étudions plus particulièrement le schéma "VAG" qui fait intervenir des inconnues au centre des mailles et aux sommets du maillage. Ce schéma conduit à la définition de flux entre le centre d'une maille et ses sommets, qui sont utilisés pour généraliser la méthode "VAG" au contexte polyphasique. Des tests numériques montrent alors que ce schéma est robuste, et conduit à un très bon compromis précision/coût, ce qui en fait un candidat idoine pour les applications industrielles. Nous présentons notamment un cas test, basé sur des observations de terrains, d'injection et de dissolution de CO2 dans la région proche d'un puits foré dans un aquifère salin. Nous montrons alors que le schéma numérique permet de simuler l'assèchement et la précipitation de minéral observée en pratique. Un chapitre de la thèse est enfin consacré à l'étude pratique et théorique d'une méthode numérique générique pour contrôler l'effet d'axe lors de l'utilisation de schémas industriels

Cette thèse est dédiée au développement et à l'analyse de schémas numériques de type volumes finis pour des équations de convection-diffusion, qui apparaissent notamment dans des modèles issus de la physique ou de la biologie. Nous nous intéressons plus particulièrement à la préservation de comportements asymptotiques au niveau discret. Ce travail s'articule en trois parties, composées chacune de deux chapitres. Dans la première partie, nous considérons la discrétisation du système de dérive diffusion linéaire pour les semi-conducteurs par le schéma de Scharfetter-Gummel implicite en temps. Nous nous intéressons à la préservation par ce schéma de deux types d'asymptotiques : l'asymptotique en temps long et la limite quasi-neutre. Nous démontrons des estimations d'énergie-dissipation d'énergie discrètes qui permettent de prouver d'une part la convergence en temps long de la solution approchée vers une approximation de l'équilibre thermique, d'autre part la stabilité à la limite quasi-neutre du schéma. Dans la deuxième partie, nous nous intéressons à des schémas volumes finis préservant l'asymptotique en temps long dans un cadre plus général. Plus précisément, nous considérons des équations de type convection-diffusion non linéaires qui apparaissent dans plusieurs contextes physiques : équations des milieux poreux, système de dérive-diffusion pour les semi-conducteurs... Nous proposons deux discrétisations en espace permettant de préserver le comportement en temps long des solutions approchées. Dans un premier temps, nous étendons la définition du flux de Scharfetter-Gummel pour une diffusion non linéaire. Ce schéma fournit des résultats numériques satisfaisants si la diffusion ne dégénère pas. Dans un second temps, nous proposons une discrétisation dans laquelle nous prenons en compte ensemble les termes de convection et de diffusion, en réécrivant le flux sous la forme d'un flux d'advection. Le flux numérique est défini de telle sorte que les états d'équilibre soient préservés, et nous utilisons une méthode de limiteurs de pente pour obtenir un schéma précis à l'ordre deux en espace, même dans le cas dégénéré. Enfin, la troisième et dernière partie est consacrée à l'étude d'un schéma numérique pour un modèle de chimiotactisme avec diffusion croisée pour lequel les solutions n'explosent pas en temps fini, quelles que soient les données initiales. L'étude de la convergence du schéma repose sur une estimation d'entropie discrète nécessitant l'utilisation de versions discrètes d'inégalités fonctionnelles telles que les inégalités de Poincaré-Sobolev et de Gagliardo-Nirenberg-Sobolev. La démonstration de ces inégalités fait l'objet d'un chapitre indépendant dans lequel nous proposons leur étude dans un contexte assez général, incluant notamment le cas de conditions aux limites mixtes et une généralisation au cadre des schémas DDFV.

Cette thèse est dédiée au développement et à l’analyse de schémas numériques de type volumes finis pour des équations de convection-diffusion, qui apparaissent notamment dans des modèles issus de la physique ou de la biologie. Nous nous intéressons plus particulièrement à la préservation de comportements asymptotiques au niveau discret. Ce travail s’articule en trois parties, composées chacune de deux chapitres. Dans la première partie, nous considérons la discrétisation du système de dérive diffusion linéaire pour les semi-conducteurs par le schéma de Scharfetter-Gummel implicite en temps. Nous nous intéressons à la préservation par ce schéma de deux types d’asymptotiques : l’asymptotique en temps long et la limite quasi-neutre. Nous démontrons des estimations d’énergie–dissipation d’énergie discrètes qui permettent de prouver d’une part la convergence en temps long de la solution approchée vers une approximation de l’équilibre thermique, d’autre part la stabilité à la limite quasi-neutre du schéma. Dans la deuxième partie, nous nous intéressons à des schémas volumes finis préservant l’asymptotique en temps long dans un cadre plus général. Plus précisément, nous considérons des équations de type convection-diffusion non linéaires qui apparaissent dans plusieurs contextes physiques : équations des milieux poreux, système de dérive-diffusion pour les semi-conducteurs... Nous proposons deux discrétisations en espace permettant de préserver le comportement en temps long des solutions approchées. Dans un premier temps, nous étendons la définition du flux de Scharfetter-Gummel pour une diffusion non linéaire. Ce schéma fournit des résultats numériques satisfaisants si la diffusion ne dégénère pas. Dans un second temps, nous proposons une discrétisation dans laquelle nous prenons en compte ensemble les termes de convection et de diffusion, en réécrivant le flux sous la forme d’un flux d’advection. Le flux numérique est défini de telle sorte que les états d’équilibre soient préservés, et nous utilisons une méthode de limiteurs de pente pour obtenir un schéma précis à l’ordre deux en espace, même dans le cas dégénéré. Enfin, la troisième et dernière partie est consacrée à l’étude d’un schéma numérique pour un modèle de chimiotactisme avec diffusion croisée pour lequel les solutions n’explosent pas en temps fini, quelles que soient les données initiales. L’étude de la convergence du schéma repose sur une estimation d’entropie discrète nécessitant l’utilisation de versions discrètes d’inégalités fonctionnelles telles que les inégalités de Poincaré-Sobolev et de Gagliardo-Nirenberg-Sobolev. La démonstration de ces inégalités fait l’objet d’un chapitre indépendant dans lequel nous proposons leur étude dans un contexte assez général, incluant notamment le cas de conditions aux limites mixtes et une généralisation au cadre des schémas DDFV. / This dissertation is dedicated to the development and analysis of finite volume numericals chemes for convection-diffusion equations, which notably occur in models arising from physics and biology. We are more particularly interested in preserving asymptotic behavior at the discrete level. This dissertation is composed of three parts, each one including two chapters. In the first part, we consider the discretization of the linear drift-diffusion system for semiconductors with the implicit Scharfetter-Gummel scheme. We focus on preserving two kinds of asymptotics with this scheme : the long-time asymptotic and the quasineutral limit. We show discrete energy–energy dissipation estimates which constitute the main point to prove first the large time convergence of the approximate solution to an approximation of the thermal equilibrium, and then the stability at the quasineutral limit. In the second part, we are interested in designing finite volume schemes which preserve the long time behavior in a more general framework. More precisely, we consider nonlinear convection-diffusion equations arising in various physical models : porous media equation, drift-diffusion system for semiconductors... We propose two spatial discretizations which preserve the long time behavior of the approximate solutions. We first generalize the Scharfetter-Gummel flux for a nonlinear diffusion. This scheme provides satisfying numerical results if the diffusion term does not degenerate. Then we propose a discretization which takes into account together the convection and diffusion terms by rewriting the flux as an advective flux. The numerical flux is then defined in such a way that equilibrium states are preserved, and we use a slope limiters method so as to obtain second order space accuracy, even in the degenerate case. Finally, the third part is devoted to the study of a numerical scheme for a chemotaxis model with cross diffusion, for which the solutions do not blow up in finite time, even for large initial data. The proof of convergence is based on a discrete entropy estimate which requires the use of discrete functional inequalities such as Poincaré-Sobolev and Gagliardo-Nirenberg-Sobolev inequalities. The demonstration of these inequalities is the subject of an independent chapter in which we propose a study in quite a general framework, including mixed boundary conditions and generalization to DDFV schemes.

Ce travail de thèse porte sur l'extension et l'analyse d'un solveur volumes finis eulérien, co-localisé avec capture d'interfaces pour la simulation des écoulements multi-matériaux non miscibles. Les extensions proposées s'inscrivent dans la volonté d'élaborer un outil de simulation multi-physiques. Dans le cadre de ce mémoire, le caractère multi-physiques recouvre les champs que nous allons détailler. Nous traitons le cas des écoulements radiatifs modélisés par un système à deux températures qui couple les phénomènes purement hydrodynamiques aux phénomènes radiatifs. Nous proposons un solveur permettant la prise en compte des effets de tension superficielle à l'interface entre deux fluides. Nous développons un solveur implicite permettant la simulation précise d'écoulements faisant intervenir de faibles nombres de Mach par le biais d'une méthode de renormalisation de la diffusion numérique. Enfin, les effets tri-dimensionnels sont considérés ainsi que la possibilité d'étendre le schéma de base aux écoulements à un nombre quelconque de matériaux. A chaque étape, les solveurs développés sont validés sur des cas-tests. / This work is devoted to the extension of a eulerian cell-centered finite volume scheme with interfaces capturing for the simulation of multimaterial fluid flows. Our purpose is to develop a simulation tool which could be able to handle multi-physics problems in the following sense. We address the case of radiating flows, modeled by a two temperature system of equations where the hydrodynamics are coupled to radiation transport. We address a numerical scheme for taking surface tension forces into account. An implicit scheme is proposed to handle low Mach number fluid flows by means of a renormalization of the numerical diffusion. Eventually, the scheme is extended to three-dimensional flows and to multimaterial flows, that is with an arbitrary number of materials. At each step, numerical simulations validate our schemes.

Cette thèse est consacrée à l'étude de méthodes numériques pour la simulation des écoulements polyphasiques en milieu poreux, en vue de leur application à des problèmes d'ingénierie pétrolière ou environnementale. Nous présentons une formulation générique du modèle d'écoulements à nombre quelconque de composants présents dans un nombre quelconque de phases. Dans notre approche l'approximation des flux diffusifs (issus, par exemple, de la loi de Darcy) s'appuie sur de nouveaux schémas, appelés schémas gradient, qui ont plusieurs avantages sur les schémas industriels standard : ces derniers, qui sont des schémas volumes finis multi-points centrés aux mailles, ne sont généralement pas symétriques et convergent difficilement sur des cas à forts rapports d'anisotropie. Nous montrons en revanche que les schémas gradient conduisent naturellement à des approximations symétriques et convergentes. Parmi cette classe de schémas, nous étudions plus particulièrement le schéma "VAG" qui fait intervenir des inconnues au centre des mailles et aux sommets du maillage. Ce schéma conduit à la définition de flux entre le centre d'une maille et ses sommets, qui sont utilisés pour généraliser la méthode "VAG" au contexte polyphasique. Des tests numériques montrent alors que ce schéma est robuste, et conduit à un très bon compromis précision/coût, ce qui en fait un candidat idoine pour les applications industrielles. Nous présentons notamment un cas test, basé sur des observations de terrains, d'injection et de dissolution de CO2 dans la région proche d'un puits foré dans un aquifère salin. Nous montrons alors que le schéma numérique permet de simuler l'assèchement et la précipitation de minéral observée en pratique. Un chapitre de la thèse est enfin consacré à l'étude pratique et théorique d'une méthode numérique générique pour contrôler l'effet d'axe lors de l'utilisation de schémas industriels / This thesis is focused on numerical methods dedicated to the simulation of multiphase flow in porous media, involved in petroleum or environmental engineering. We present a generic formulation of the flow model which is able to take into account any number of components within any number of phases. In our approach the approximation of the diffusive fluxes (mainly resulting from Darcy's law) is based on new schemes, called gradient schemes, which show several advantages over the standard industrial numerical schemes : these schemes, which belong to the class of the cell-centred MultiPoint Flux Approximation finite volume schemes, are not symmetric and may lead to difficulties of convergence in the case of high anisotropy ratios. We indeed show that gradient schemes are naturally providing symmetric and convergent approximations. We particularly study the Vertex Approximate Gradient scheme (called the VAG scheme), whose discrete unknowns are the values at the cell centres and at the vertices of the mesh. This scheme implies the definition of fluxes between the centre of a given cell and its vertices, used for the extension of the scheme to multiphase flow. Numerical tests show the robustness and the accuracy of the method for a low computational cost, which enables the use of the VAG scheme in an industrial framework. A test case, based on experimental data of injection and dissolution of CO2 in the near-well region within a saline aquifer, shows the aptitude of the scheme for reproducing drying and salt precipitation, which are practically observed. Finally, a chapter of the thesis is devoted to the theoretical and practical study of a general numerical method for controlling Grid Orientation Effect in industrial simulators

Ce travail s'inscrit dans l'étude mathématique d'écoulements de fluides environnementaux. Nous en abordons deux aspects, à travers deux contextes distincts d'application.En lien avec la simulation des écoulements en milieux poreux, on s'intéresse dans une première partie à la discrétisation d'opérateurs de diffusion anisotropes hétérogènes par des méthodes de volumes finis sur des maillages généraux. Dans le but d'obtenir des solutions approchées qui respectent les bornes physiques des modèles, notre attention se porte sur la conservation du principe du maximum pour les opérateurs elliptiques. Nous présentons des mécanismes généraux permettant de corriger tout schéma volumes finis afin de garantir un principe du maximum discret tout en préservant certaines de ses propriétés principales. On étudie en particulier les propriétés de coercivité et de convergence des schémas corrigés.La deuxième partie est consacrée à la construction de modèles approchés pour la propagation des vagues en eaux peu profondes et sur des topographies irrégulières. A cet effet, nous proposons tout d'abord une adaptation de la démarche d'étude classique à des écoulements bidimensionnels sur des topographies polygonales. Dans un cadre plus général, nous développons ensuite une démarche formelle qui débouche sur des alternatives non locales à quelques modèles classiques (équations de Saint-Venant, équations de Serre, système de Boussinesq). Ces nouveaux modèles contiennent des termes régularisants pour les contributions du fond. / This work investigates two research questions associated with environmental flows and their mathematical modeling.The first part is devoted to the development of finite volume methods for anisotropic and heterogeneous diffusion operators arising in models of porous media flows. To ensure that the approximate solutions lie within physical bounds, we aim at maintaining a discrete analogous of the maximum principle for elliptic operators. Starting from any given cell-centered finite volume scheme, we present a general approach to devise non-linear corrections providing a discrete maximum principle while retaining some main properties of the scheme. In particular, we study the coercivity and convergence properties of the modified schemes.The second part of this work focuses on the derivation of approximate models for shallow water wave propagation over rough topographies. In the particular case of one-dimensional polygonal bottom profiles, we first propose an adaptation of the usual derivation method using complex analysis tools. We then develop a formal approach to account for more general topographies. We propose nonlocal alternatives to some classical models (namely Saint-Venant equations, Serre equations and Boussinesq system). All these alternative models only involve smoothing contributions of the bottom.

Ce travail s'inscrit dans l'étude mathématique d'écoulements de fluides environnementaux. Nous en abordons deux aspects, à travers deux contextes distincts d'application.En lien avec la simulation des écoulements en milieux poreux, on s'intéresse dans une première partie à la discrétisation d'opérateurs de diffusion anisotropes hétérogènes par des méthodes de volumes finis sur des maillages généraux. Dans le but d'obtenir des solutions approchées qui respectent les bornes physiques des modèles, notre attention se porte sur la conservation du principe du maximum pour les opérateurs elliptiques. Nous présentons des mécanismes généraux permettant de corriger tout schéma volumes finis afin de garantir un principe du maximum discret tout en préservant certaines de ses propriétés principales. On étudie en particulier les propriétés de coercivité et de convergence des schémas corrigés.La deuxième partie est consacrée à la construction de modèles approchés pour la propagation des vagues en eaux peu profondes et sur des topographies irrégulières. A cet effet, nous proposons tout d'abord une adaptation de la démarche d'étude classique à des écoulements bidimensionnels sur des topographies polygonales. Dans un cadre plus général, nous développons ensuite une démarche formelle qui débouche sur des alternatives non locales à quelques modèles classiques (équations de Saint-Venant, équations de Serre, système de Boussinesq). Ces nouveaux modèles contiennent des termes régularisants pour les contributions du fond.

Dans cette thèse on étudie des systèmes d’EDP non linéaires modélisant la propagation électromagnétique dans des milieux de type Kerr. On considère deux modèles. Le premier dit de Kerr-Debye, suppose un temps de réponse non nul du matériau à l’onde électromagnétique. Le second, dit de Kerr, suppose une réponse instantanée. On est ainsi confronté à des systèmes de relaxation tels que définis par Chen-Levermore-Liu (CPAM 1994). Nous établissons ici des résultats d’existence globale de solutions fortes à données petites en 3D pour le problème de Cauchy et un problème mixte. Puis nous construisons des schémas volumes finis asymptotic preserving et nous étudions leurs performances sur des cas physiques. / In this thesis, we study non-linear PDE systems modeling the electromagnetic propagation in Kerr media. We consider two models. The first one is the Kerr-Debye model, it assumes a finite response time of the medium. The second one is the Kerr model, it assumes an instantaneous response. We deal with relaxation systems as defined by Chen-Levermore-Liu (CPAM 1994). For small data, we establish results of global existence of smooth solutions in 3D for the Cauchy problem and the IBVP. Then we investigate asymptotic preserving finite volume schemes and we study their performance on physical cases.