Cette thèse est centrée autour de l’étude théorique et de l’analyse numérique des équations paraboliques non linéaires avec divers conditions aux limites. La première partie est consacrée aux équations paraboliques dégénérées mêlant des phénomènes non-linéaires de diffusion et de transport. Nous définissons des notions de solutions entropiques adaptées pour chacune des conditions aux limites (flux nul, Robin, Dirichlet). La difficulté principale dans l’étude de ces problèmes est due au manque de régularité du flux pariétal pour traiter les termes de bords. Ceci pose un problème pour la preuve d’unicité. Pour y remédier, nous tirons profit du fait que ces résultats de régularités sur le bord sont plus faciles à obtenir pour le problème stationnaire et particulièrement en dimension un d’espace. Ainsi par la méthode de comparaison "fort-faible" nous arrivons à déduire l’unicité avec le choix d’une fonction test non symétrique et en utilisant la théorie des semi-groupes non linéaires.L’existence de solution se démontre en deux étapes, combinant la méthode de régularisation parabolique et les approximations de Galerkin. Nous développons ensuite une approche directe en construisant des solutions approchées par un schéma de volumes finis implicite en temps. Dans les deux cas, on combine les estimations dans les espaces fonctionnels bien choisis avec des arguments de compacité faible ou forte et diverses astuces permettant de passer à la limite dans des termes non linéaires. Notamment, nous introduisons une nouvelle notion de solution appelée solution processus intégrale dont l’objectif, dans le cadre de notre étude, est de pallier à la difficulté de prouver la convergence vers une solution entropique d’un schéma volumes finis pour le problème de flux nul au bord.La deuxième partie de cette thèse traite d’un problème à frontière libre décrivant la propagation d’un front de combustion et l’évolution de la température dans un milieu hétérogène. Il s’agit d’un système d’équations couplées constitué de l’équation de la chaleur bidimensionnelle et d’une équation de type Hamilton-Jacobi. L’objectif de cette partie est de construire un schéma numérique pour ce problème en combinant des discrétisations du type éléments finis avec les différences finies. Ceci nous permet notamment de vérifier la convergence de la solution numérique vers une solution onde pour un temps long. Dans un premier temps, nous nous intéressons à l’étude d’un problème unidimensionnel. Très vite,nous nous heurtons à un problème de stabilité du schéma. Cela est dû au problème de prise en compte de la condition de Neumann au bord. Par une technique de changement d’inconnue et d’approximation nous remédions à ce problème. Ensuite, nous adaptons cette technique pour la résolution du problème bidimensionnel. A l’aide d’un changement de variables, nous obtenons un domaine fixe facile pour la discrétisation. La monotonie du schéma obtenu est prouvée sous une hypothèse supplémentaire de propagation monotone qui exige que la frontière libre se déplace dans les directions d’un cône prescrit à l’avance. / This thesis focuses on the theoretical study and numerical analysis of parabolic equations with boundary conditions.The first part is devoted to degenerate parabolic equation which combines features of a hyperbolic conser-vation law with those of a porous medium equation. We define suitable notions of entropy solutions foreach of the boundary conditions (zero-flux, Robin, Dirichlet). The main difficulty in these studies residesin the formulation of the adequate notion of entropy solution and in the proof of uniqueness. There isa technical difficulty due to the lack of regularity required to treat the boundaries terms. We take ad-vantage of the fact that boundary regularity results are easier to obtain for the stationary problem, inparticular in one space dimension. Thus, using strong-weak uniqueness approach we get the uniquenesswith the choice of a non-symmetric test function and using the nonlinear semigroup theory. The exis-tence of solution is proved in two steps, combining the method of parabolic regularization and Galerkinapproximations. Next, we develop a direct approach to construct approximate solutions by an implicitfinite volume scheme. In both cases, the estimates in the appropriately chosen functional spaces are com-bined with arguments of weak or strong compactness and various tricks to pass to the limit in nonlinearterms. In the appendix, we propose a result of existence of strong trace of a solution for the degenerateparabolic problem. In another appendix of independent interest, we introduce a new concept of solutioncalled integral process solution. We exploit it to overcome the difficulty of proving the convergence ofour finite volume scheme to an entropy solution for the zero-flux boundary problem.The second part of this thesis deals with a free boundary problem describing the propagation of a com-bustion front and the evolution of the temperature in a heterogeneous medium. So we have a coupledproblem consisting of the heat equation of bidimensional space and a Hamilton-Jacobi equation. The ob-jective is to construct a numerical scheme and to verify that the numerical solution converges to a wavesolution for a long time. Recall that an existence of wave solution for this problem was already proven inan analytical framework. At first, we focus on the study of a one-dimensional problem. Here, we face aproblem of stability of the scheme. This is due to a difficulty of taking into account the Neumann boun-dary condition. Through a technique of change of unknown, we can propose a monotone scheme. Wealso adapt this technique for solving two-dimensional problem. Using a change of variables, we obtaina fixed domain where the discretization becomes easy. The monotony of the scheme is proved under anadditional assumption of monotone propagation that requires the free boundary moves in the directionsof a cone given beforehand.

Cette thèse porte sur la modélisation et la discrétisation d’écoulements Darcy dans les milieux poreux fracturés. Nous suivons l’approche des modèles, dits à dimensions hybrides, qui représentent les réseaux de fractures comme des surfaces de codimension 1 immergées dans la matrice. Les modèles considérés prennent en compte les interactions entre matrice et fractures et permettent de traiter des fractures agissant comme conduites ou comme barrières, ce que nécessite de prendre en compte les sauts de pression aux interfaces matrice-fracture. Dans le cas des écoulements diphasiques, nous proposons des modèles, qui prennent en compte les sauts de saturations aux interfaces matrice-fracture, dû à la capillarité. L’analyse numérique est menée dans le cadre général de la méthode de discrétisations gradients, qui est étendue aux modèles considérés. Deux familles de schémas numériques, le schéma Vertex Approximate Gradient et le schéma Volumes Finis Hybrides sont adaptées aux modèles à dimensions hybrides. On prouve via des résultats de densité que ce sont des schémas gradients, pour lesquels la convergence est établie. En diphasique, l’existence d’une solution est obtenue en passant. Plusieurs cas tests sont présentés. En monophasique, on observe la convergence sur des différents types de mailles pour une famille de solutions dans un milieux fracturé hétérogène et anisotrope. En diphasique, nous présentons une série de cas tests afin de comparer les modèles à dimensions hybrides au modèle de référence, dans lequel les fractures ont la même dimension que la matrice. A part quantifier le gain en performance de calcul, ces tests montrent la qualité des différents modèles réduits. / This thesis investigates the modelling of Darcy flow through fractured porous media and its discretization on general polyhedral meshes. We follow the approach of hybrid dimensional models, invoking a complex network of planar fractures. The models account for matrix-fracture interactions and fractures acting either as drains or as barriers, i.e. we have to deal with pressure discontinuities at matrix-fracture interfaces. In the case of two phase flow, we present two models, which permit to treat gravity dominated flow as well as discontinuous capillary pressure at the material interfaces. The numerical analysis is performed in the general framework of the Gradient Discretisation Method, which is extended to the models under consideration. Two families of schemes namely the Vertex Approximate Gradient scheme (VAG) and the Hybrid Finite Volume scheme (HFV) are detailed and shown to fit in the gradient scheme framework, which yields, in particular, convergence. For single phase flow, we obtain convergence of order 1 via density results. For two phase flow, the existence of a solution is obtained as a byproduct of the convergence analysis. Several test cases are presented. For single phase flow, we study the convergence on different types of meshes for a family of solutions. For two phase flow, we compare the hybrid-dimensional models to the reference equidimensional model, in which fractures have the same dimension as the matrix. This does not only provide quantitative evidence about computational gain, but also leads to deep insight about the quality of the proposed reduced models.

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.

Nous proposons une étude numérique bidimensionnelle des écoulements laminaires de convection mixte en conduite rectangulaire horizontale chauffée par le bas (encore appelés "écoulements de Poiseuille-Rayleigh-Bénard"). La méthode de résolution est de type volumes finis; elle est basée sur la méthode du Lagrangien augmenté. Nous présentons une revue bibliographique, historique et complète, des travaux effectués sur ce sujet. Nous rapportons de manière détaillée les résultats de trois thèses, analysant la stabilité linéaire de l'écoulement de Poiseuille et la stabilité faiblement non-linéaire des rouleaux thermoconvectifs transverses à l'axe de la conduite, et fournissant des mesures expérimentales très fines des champs de vitesse dans les domaines linéaires et non-linéaires. L'objectif de ce travail est d'étudier dans quelle mesure les simulations numériques s'accordent avec théories et expériences en effectuant des comparaisons quantitatives systématiques avec ces trois thèses. Nous étudions l'influence des conditions aux limites ouvertes et périodiques sur la stabilité et sur le développement spatial et temporel des rouleaux transversaux, en prenant en compte le caractère convectif ou absolu de l'instabilité. Cinq sortes de conditions aux limites ouvertes sont testées. On montre que l'amplitude de la perturbation en amont de la frontière de sortie peut considérablement varier d'une condition à l'autre, alors que la longueur de la zone perturbée varie très peu, et que la structure de l'écoulement en dehors de la zone perturbée se conserve. Avec les conditions aux limites périodiques, les rouleaux transversaux transitent vers l'écoulement de Poiseuille au passage du seuil entre instabilité convective et stabilité linéaire. Avec les conditions aux limites ouvertes, cette transition a lieu près du seuil entre instabilité absolue et instabilité convective. Tous les résultats de la stabilité faiblement non-linéaire sont reproduits par la simulation 2D, et certains résultats expérimentaux 3D sont retrouvés avec une grande précision. On montre que les points de désaccord qui sont mis en évidence résultent d'effets purement tridimensionnels et qu'ils ne pourront donc être traités que par la simulation numérique 3D.

Lors de la rentrée atmosphérique nous sommes amenés à modéliser trois phénomènes physiques différents. Tout d'abord, l'écoulement autour du véhicule entrant dans l'atmosphère est hypersonique, il est caractérisé par la présence d'un choc fort et provoque un fort échauffement du véhicule. Nous modélisons l'écoulement par les équations de Navier-Stokes compressibles et l'échauffement du véhicule au moyen de la thermique anisotrope. De plus le véhicule est protégé par un bouclier thermique siège de réactions chimiques que l'on nomme communément ablation.<br /><br /> Dans le premier chapitre de cette thèse nous présentons le schéma numérique de diffusion CCLAD (Cell-Centered LAgrangian Diffusion) que nous utilisons pour résoudre la thermique anisotrope. Nous présentons l'extension en trois dimensions de ce schéma ainsi que sa parallélisation.<br /> Nous continuons le manuscrit en abordant l'extension de ce schéma à une équation de diffusion tensorielle. Cette équation est obtenue en supprimant les termes convectifs de l'équation de quantité de mouvement des équations de Navier-Stokes. Nous verrons qu'une pénalisation doit être introduite afin de pouvoir inverser la loi constitutive et ainsi appliquer la méthodologie CCLAD. Nous présentons les propriétés numériques du schéma ainsi obtenu et effectuons des validations numériques.<br /> Dans le dernier chapitre, nous présentons un schéma numérique de type Volumes Finis permettant de résoudre les équations de Navier-Stokes sur des maillages non-structurés obtenu en réutilisant les deux schémas de diffusion présentés précédemment.

L'objet de la thèse est de contribuer à l'étude numérique des lois de conservation hyperboliques avec termes sources, ce qui est motivé par les applications aux équations de Saint-Venant pour les eaux peu profondes. La première partie traite des questions habituelles de l'analyse des approximations numériques des lois de conservation scalaires. On se concentre sur des schémas aux volumes finis semi-discrets, dans le cas général d'un maillage non-uniforme. Pour définir des discrétisations appropriées du terme source, on introduit le formalisme spécifique de la méthode "Upwind Interface Source" et on établit des conditions sur les fonctions numériques telles que le solveur discret préserve les solutions stationnaires. Une définition rigoureuse de consistance est ensuite formulée, adaptée aux "schémas équilibres", pour laquelle on est capable de prouver un théorème de convergence faible de type Lax-Wendroff. La méthode considérée dans un premier temps est essentiellement d'ordre un en espace. Pour améliorer la précision, on développe des approches à haute résolution pour la méthode "Upwind Interface Source" et on montre que celles-ci sont un moyen efficace de dériver des schémas d'ordre plus élevé avec des propriétés convenables. On prouve une estimation d'erreur dans $L^p$, $1\le p < +\infty$, qui est un résultat optimal dans le cas d'un maillage uniforme. On conclut alors que les mêmes taux de convergence $O(h)$ et $O(h^2)$ que pour les systèmes homogènes correspondants sont valables. La deuxième partie présente un schéma numérique pour approcher les équations de Saint-Venant, avec un terme source géométrique, qui vérifie les propriétés théoriques suivantes: il préserve les états stationnaires de l'eau au repos, vérifie une inégalité d'entropie discrète, préserve la positivité de la hauteur de l'eau et reste stable avec des profiles du fond discontinus. Cela est obtenu grâce à une approche cinétique au système; dans ce contexte, on utilise une description formelle du comportement microscopique du système pour définir les flux numériques aux interfaces d'un maillage non-structuré. On utilise aussi le concept de variables conservatives centrées (typique de la méthode des volumes finis) et des termes sources décentrés aux interfaces. Finalement, on présente des simulations numériques du système des équations de Saint-Venant modifiées pour prendre en compte le frottement et la viscosité, afin de retrouver les résultats de certaines études expérimentales. Une application à la modélisation des termes de frottement pour les avalanches de neige est discutée dans l'Appendice.