Dans cette thèse, nous étudions l'équation de la chaleur couplée avec la loi de Darcy à travers de la viscosité non-linéaire qui dépend de la température pour les dimensions d=2,3 (Hooman et Gurgenci ou Rashad). Nous analysons ce problème en introduisant la formulation variationnelle équivalente et en la réduisant à une simple équation de diffusion-convection pour la température où la vitesse dépend implicitement de la température.Nous démontrons l'existence de la solution sans la restriction sur les données par la méthode de Galerkin et du point fixe de Brouwer. L'unicité globale est établie une fois la solution est légèrement régulière et les données se restreignent convenablement. Nous introduisons aussi une formulation variationnelle alternative équivalente. Toutes les deux formulations variationnelles sont discrétisées par quatre schémas d'éléments finis pour un domaine polygonal ou polyédrique. Nous dérivons l'existence, l'unicité conditionnée, la convergence et l'estimation d'erreur a priori optimale pour les solutions des trois schémas. Par la suite, ces schémas sont linéarisés par des algorithmes d'approximation successifs et convergentes. Nous présentons quelques expériences numériques pour un problème modèle qui confirme les résultats théoriques de convergence développées dans ce travail. L'estimation d'erreur a posteriori est établie avec deux types d'indicateurs d'erreur de linéarisation et de discrétisation. Enfin, nous montrons des résultats numériques de validation. / In this thesis, we study the heat equation coupled with Darcy's law by a nonlinear viscosity depending on the temperature in dimension d=2,3 (Hooman and Gurgenci or Rashad). We analyse this problem by setting it in an equivalent variational formulation and reducing it to an diffusion-convection equation for the temperature where the velocity depends implicitly on the temperature.Existence of a solution is derived without restriction on the data by Galerkin's method and Brouwer's Fixed Point. Global uniqueness is established when the solution is slightly smoother and the dataare suitably restricted. We also introduce an alternative equivalent variational formulation. Both variational formulations are discretized by four finite element schemes in a polygonal or polyhedral domain. We derive existence, conditional uniqueness, convergence, and optimal a priori error estimates for the solutions of the three schemes. Next, these schemes are linearized by suitable convergent successive approximation algorithms. We present some numerical experiments for a model problem that confirm the theoretical rates of convergence developed in this work. A posteriori error estimates are established with two types of errors indicators related to the linearisation and discretization. Finally, we show numerical results of validation.

Les simulations en mécanique du solide présentent des difficultés comme le traitement de l'incompressibilité ou les non-linéarités dues aux grandes déformations, aux lois de comportement et de contact. L'objectif principal de ce travail est de proposer des méthodes éléments finis capables de gérer l'incompressibilité en grandes déformations en utilisant des éléments de faible ordre. Parmi les approches de la littérature, les formulations mixtes offrent un cadre théorique intéressant. Dans ce travail, une formulation mixte à trois champs (déplacements, pression, gonflement) est introduite. Dans certains cas, cette formulation peut être condensée en formulation à deux champs. Cependant, il est connu que le problème discret obtenu par une approche éléments finis de type Galerkin n'hérite pas automatiquement de la condition de stabilité “inf-sup” du problème continu : les éléments finis utilisés, et notamment les ordres d'interpolation doivent être choisis de sorte à vérifier cette condition de stabilité. Cependant, il est possible de s'affranchir de cette contrainte en ajoutant des termes de stabilisation à la formulation EF Galerkin. Cette approche permet entre autres d'utiliser des ordres d'interpolation égaux. Dans ce travail, des éléments finis stables de type P2/P1 sont utilisés comme référence, et comparés à une formulation P1/P1, stabilisée soit avec une fonction bulle, soit avec une méthode VMS (Variational Multi-Scale) basée sur un espace sous-grille orthogonal à l'espace EF. Combinées à un modèle grandes déformations basé sur des déformations logarithmiques, ces approches sont d'abord validées sur des cas académiques puis sur des cas industriels. / Simulations in solid mechanics exhibit difficulties as dealing with incompressibility or nonlinearities due to finite strains, constitutive laws and contact. The basic motivation of our work is to propose efficient finite element methods capable of dealing with incompressibility in finite strain context, and using low order elements. Among the approaches in the literature, mixed formulations offer an interesting theoretical framework. In this work, a three-field mixed formulation (displacement, pressure, volumetric strain) is investigated. In some cases, this formulation can be condensed in a two-field formulation. However, it is well-known that the discrete problem given by the Galerkin finite element technique, does not inherit the “inf-sup” stability condition from the continuous problem: the finite elements used, and in particular the interpolation orders must be chosen so as to satisfy this stability condition. However, it is possible to circumvent it, by adding terms stabilizing the FE Galerkin formulation. The latter approach allows the use of equal order interpolation. In this work, stable finite elements of type P2/P1 are used as reference, and compared to a P1/P1 formulation, stabilized with a bubble function, or with a VMS method (Variational Multi-Scale) based on a sub-grid-space orthogonal to the FE space. Combined to a finite strain model based on logarithmic strain, these approaches are first validated on academic cases and then on industrial cases.