Dans cette thèse nous nous intéressons à l'analyse d'erreur a priori et a posteriori de méthodes d'éléments finis mixtes et non-conformes. Nous considérons en particulier les équations de Darcy à perméabilité fortement variable et les équations de convection-diffusion-réaction en régime de convection dominante. Nous discrétisons les équations de Darcy par une méthode d'éléments finis mixtes non-conformes de type Petrov-Galerkin appelée schéma boîte. Les techniques d'estimations d'erreur a posteriori par résidu et hiérarchique conduisent à des estimateurs d'erreur a posteriori fiables et optimaux indépendamment des fluctuations de la perméabilité. Les résultats théoriques sont validés numériquement sur différents cas tests présentant de forts contrastes de perméabilité. Enfin, nous montrons comment les indicateurs d'erreur obtenus permettent de générer des maillages adaptatifs. Nous discrétisons les équations de convection-diffusion-réaction par des éléments finis nonconformes. Deux méthodes de stabilisation sont étudiées: la stabilisation par viscosité de sous-maille, conduisant à un schéma boîte et la méthode de pénalisation sur les faces. Nous montrons que les deux schémas ainsi obtenus ont les mêmes propriétés de convergence que les approximations par éléments finis conformes. Grâce aux techniques d'estimations d'erreur par résidu nous obtenons des estimateurs d'erreur a posteriori fiables et optimaux. Certains des indicateurs d'erreur sont robustes au sens de Verfürth, c'est à dire que le rapport des constantes intervenant dans les inégalités de fiabilité et d'optimalité explose en au plus l'inverse du nombre de Péclet. Les résultats théoriques sont validés numériquement et les indicateurs d'erreur a posteriori obtenus permettent de générer des maillages adaptatifs sur des problèmes présentant des couches intérieures.

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.

Dans ce mémoire on s'intéresse à la modélisation du changement de phase solide-liquide en convection naturelle et en convection forcée en utilisant l'équation de diffusion de la chaleur et une forme modifiée de la conductivité. Comme on ne cherche pas à résoudre le champ des vitesses, on intègre à l'intérieur de la conductivité modifiée les effets thermiques des mouvements convectifs, puis on résout l'équation de conduction en utilisant ces conductivités modifiées. L'objectif du mémoire est de prédire la position de l'interface solide-liquide en fonction du temps dans un processus de fusion ou de solidification en présence de convection, et ce avec un modèle conductif. Dans les deux cas types que nous avons traités, la valeur de la conductivité modifiée à utiliser dans le modèle conductif simplifié est d'abord estimée par une étude d'ordre de grandeur. Le premier cas type étudié est celui du changement de phase en présence de convection forcée se déroulant dans une conduite à paroi froide. On a réussi à développer une relation permettant d'obtenir la valeur de la conductivité modifiée directionnelle à utiliser dans le modèle conductif en fonction du nombre de Reynolds de l'écoulement ainsi que du facteur de forme de la conduite. Cette relation a montré de bons résultats en comparaison avec ceux issus d'un modèle numérique complet plus classique (CFD). Le deuxième cas type étudié est celui de la fusion d'un matériau sous l'effet de la convection naturelle à l'intérieur d'une enceinte. Nous avons réussi à paramétrer l'évolution des valeurs de conductivité modifiée à utiliser par notre modèle conductif, celles-ci sont fonction de l'avancement du front de fusion et du nombre de Nusselt. Cette modélisation est comparée à la fois avec les résultats fournis par un modèle numérique complet (CFD) mais aussi avec des résultats expérimentaux issus de la littérature. Ces comparaisons ont montré une bonne concordance entre notre modélisation et la réalité. / In this work we have developed an enhanced conduction model for predicting solid-liquid interface positions in convection-dominated phase-change processus. The flow field is not calculated and the effects of convection are taken into account via the modification of the material conductivity. Our objective is to obtain a good approximation of the solid-liquid interface evolution. It is shown that the enhanced thermal conductivity of the melt may be formulated in terms of directional thermal conductivity components and that their value may be correlated in terms of dimensional numbers obtained from an order of magnitude analysis. The proposed approach is then tested for two different cases. The first test case is devoted to forced convection dominated solidification in a duct while the second test case is concerned with buoyancy dominated melting in an enclosure. The results of the simulations using our conduction model are then compared with a full CFD model and in the case of melting in an enclosure experimental data, they show good agreements.

Les questions de sûreté et d'incertitudes sont au centre des études de faisabilité pour un site de stockage souterrain de déchets nucléaires, en particulier l'évaluation des incertitudes sur les indicateurs de sûreté qui sont dues aux incertitudes sur les propriétés du sous-sol et des contaminants. L'approche globale par les méthodes probabilistes de type Monte Carlo fournit de bons résultats, mais elle demande un grand nombre de simulations. La méthode déterministe étudiée ici est complémentaire. Reposant sur la décomposition en valeurs singulières de la dérivée du modèle, elle ne donne qu'une information locale, mais elle est beaucoup moins coûteuse en temps de calcul. Le modèle d'écoulement suit la loi de Darcy et le transport des radionucléides autour du site de stockage est modélisé par une équation de diffusion-convection linéaire. Différentiation à la main et différentiation automatique sont comparées sur ces modèles en mode direct et en mode adjoint. Une étude comparée des deux approches probabiliste et déterministe pour l'analyse de la sensibilité des flux de contaminants aux exutoires par rapport aux variations des paramètres d'entrée est menée sur des données réalistes fournies par l'ANDRA. Des outils génériques d'analyse de sensibilité et de couplage de code sont développés en langage Caml. Ils permettent à l'utilisateur de ces plates-formes génériques de ne fournir que la partie spécifique de l'application dans le langage de son choix. Une étude sur les écoulements diphasiques eau/air partiellement saturés en hydrogéologie porte sur les limitations des approximations de Richards et de la formulation en pression globale issue du domaine pétrolier.