Dans de nombreux problèmes physiques, un modèle incertain peut être traduit par un ensemble d'équations aux dérivées partielles stochastiques. Nous nous intéressons ici à des problèmes présentant de nombreuses sources d'incertitudes localisées en espace. Dans le cadre des approches fonctionnelles pour la propagation d'incertitudes, ces problèmes présentent deux difficultés majeures. La première est que leurs solutions possèdent un caractère multi-échelle, ce qui nécessite des méthodes de réduction de modèle et des stratégies de calcul adaptées. La deuxième difficulté est associée à la représentation de fonctions de nombreux paramètres pour la prise en compte de nombreuses variabilités. Pour résoudre ces difficultés, nous proposons tout d'abord une méthode de décomposition de domaine multi-échelle qui exploite le caractère localisé des aléas. Un algorithme itératif est proposé, qui requiert une résolution alternée de problèmes globaux et de problèmes locaux, ces derniers étant définis sur des patchs contenant les variabilités localisées. Des méthodes d'approximation de tenseurs sont ensuite utilisées pour la gestion de la grande dimension paramétrique. La séparation multi-échelle améliore le conditionnement des problèmes à résoudre et la convergence des méthodes d'approximation de tenseurs qui est liée aux propriétés spectrales des fonctions à décomposer. Enfin, pour la prise en compte de variabilités géométriques localisées, des méthodes spécifiques basées sur les approches de domaines fictifs sont introduites.

Cette thèse est consacrée à l’étude de méthodes de domaines fictifs pour les éléments finis. Ces méthodes, initialement conçues pour l’approximation de problèmes d’interactions fluide/structure, consistent à prolonger un domaine réel par un domaine de géométrie simple appelé domaine fictif. On applique ces méthodes à un problème de contact unilatéral sans frottement en petite déformation entre deux corps élastiques séparés par une distance initiale non nulle et possédant par ailleurs des conditions aux bords de type Dirichlet et Neumann. Les deux premiers chapitres sont consacrés à l’introduction des méthodes de domaines fictifs et du problème unilatéral de contact de deux corps élastiques. Le chapitre 3 est consacré à l’analyse a priori et à l’étude numérique de ce problème de contact en domaine fictif avec les conditions aux bords de Dirichlet et de contact qui sont prises en compte à l’aide d’une méthode de type Nitsche. Des résultats théoriques de consistance de la méthode discrète, d’existence et d’unicité sont présentés. Afin d’obtenir une estimation d’erreur a priori optimale, une stabilisation de la méthode de domaine fictif est nécessaire. Ces résultats sont validés numériquement sur des cas tests en dimensions deux et trois. Le chapitre 4 est consacré à l’étude d’un estimateur d’erreur de type résidu d’un problème de contact sans domaine fictif entre un corps élastique et un corps rigide. Les résultats théoriques sont également validés sur deux cas tests numériques : un domaine rectangulaire avec seulement une partie de la zone de contact en contact effectif ainsi qu’un contact de type Hertz en dimensions deux et trois. Le chapitre 5 est une généralisation du chapitre 4 à l’approche domaine fictif et au cas de deux corps élastiques. / This thesis is dedicated to the study of the fictitious domain methods for the finite element methods. These methods, initially designed for the fluid-structure interaction, consist in immersing the real domain in a simply-shaped and a geometrically bigger domain called the fictitious domain. We apply these methods to a unilateral frictionless contact problem in small deformation of two deformable elastics bodies separated by an initial gap and satisfying boundary Dirichlet and Neumann conditions. The first two chapters are devoted to the introduction of these methods and to the unilateral contact problem. The chapter 3 is dedicated to a theoretical study for Dirichlet and contact boundary conditions taken into account with a Nitsche type method. Some theoretical results are presented: the consistency of the discrete method, existence and uniqueness results. To obtain an optimal a priori error estimate, a stabilized fictitious domain method is necessary. These results are numerically validated using Hertz contact in two and three dimensions. The chapter 4 is devoted to the study of a residual-based a posteriori error estimator, without the fictitious domain approach, between an elastic body and rigid obstacle. The numerical study of two tests cases will be performed: a rectangular domain with only a part of the potential zone of contact in effective contact as well as a Hertz contact in two and three dimensions. The chapter 5 is a generalization of the chapter 4 to the fictitious domain approach and the care of to two elastics bodies.

Notre travail est consacré à la définition, l'analyse et l'implémentation de nouveaux algorithmes numériques pour l'approximation de la solution de problèmes à 2 dimensions du type problème de Stefan. Dans ce type de problèmes une équation aux dérivée partielle parabolique posée sur un ouvert omega quelconque est couplée avec une autre équation qui contrôle la frontière gamma du domaine lui même. Les difficultés classiquement associés à ce type de problèmes sont: la formulation en particulier de l'équation pour le bord du domaine, l'approximation de la solution liées à la forme quelconque du domaine, les difficultés associées à l'implication des opérateurs de trace (approximation, conditionnement), les difficultés liées aux de régularité fonds du domaine.De plus, de nombreuse situations d'intérêt physique par exemple demandent des approximations de haut degré. Notre travail s'appuie sur une formulation de type espaces de niveaux (level set) pour l'équation du domaine, et une formulation de type domaine fictif (Omega) pour l'équation initiale.Le contrôle des conditions aux limites est effectué à partir de multiplicateurs de Lagrange agissant sur une frontière (Gamma) dite de contrôle différente de frontière(gamma) du domaine (omega). L'approximation est faite à partir d'un schéma aux différences finies pour les dérivées temporelle et une discrétisation à l'aide d'ondelettes bi-dimensionelles pour l'équation initiale et une dimensionnelle pour les multiplicateurs de Lagrange. Des opérateurs de prolongement de omega à Omega sont également construits à partir d'analyse multiéchelle sur l'intervalle. Nous obtenons aussi: une formulation pour laquelle existence de la solution est démontrées, un algorithme convergent pour laquelle une estimation globale d'erreur (sur Omega) est établie, une estimation intérieure prouvant sur l'erreur à un domaine omega, overline omega subset Xi, des estimations sur les conditionnement associés a l'opérateur de trace, des algorithmes de prolongement régulier. Différentes expériences numériques en 1D ou 2D sont effectuées. Le manuscrit est organisé comme suit: Le premier chapitre rappelle la construction des analyses multirésolutions, les propriétés importantes des ondelettes et des algorithmes numériques liées à l'application d'opérateurs aux dérivées partielles. Le second chapitre donne un aperçu des méthodes de domaine fictif classiques, approchées par la méthode de Galerkin ou de Petrov-Galerkin. Nous y découvrons les limites de ces méthodes ce qui donne la direction de notre travail. Le chapitre trois présente notre nouvelle méthode de domaine fictif que l'on appelle méthode de domaine fictif lisse.L'approximation est grâce à une méthode d'ondelettes de type Petrov-Galerkin. Cette section contient l'analyse théorique et décrit la mise en œuvre numérique. Différents avantages de cette méthode sont démontrés. Le chapitre quatre introduit une technique de prolongement régulier. Nous l'appliquons à des problèmes elliptiques en 1D ou 2D.\par Le cinquième chapitre décrit quelques simulations numériques de problème de Stefan. Nous testons l'efficacité de notre méthode sur différents exemples dont le problème de Stefan à 2 phases avec conditions aux limites de Gibbs-Thomson. / Our work is devoted to the definition, analysis and implementation of a new algorithms for numerical approximation of the solution of 2 dimensional Stefan problem. In this type of problem a parabolic partial differential equation defined on an openset Omega is coupled with another equation which controls the boundary gamma of the domain itself. The difficulties traditionally associated with this type of problems are: the particular formulation of equation on the boundary of domain, the approximation of the solution defined on general domain, the difficulties associated with the involvement of trace operation (approximation, conditioning), the difficulties associated with the regularity of domain. Addition, many situations of physical interest, for example,require approximations of high degree. Our work is based on aformulation of type level set for the equation on the domain, and aformulation of type fictitious domain (Omega) for the initialequation. The control of boundary conditions is carried out throughLagrange multipliers on boundary (Gamma), called control boundary, which is different with boundary (gamma) of the domain (omega). The approximation is done by a finite difference scheme for time derivative and the discretization by bi-dimensional wave letfor the initial equation and one-dimensional wave let for the Lagrange multipliers. The extension operators from omega to Omega are also constructed from multiresolution analysis on theinterval. We also obtain: a formulation for which the existence of solution is demonstrated, a convergent algorithm for which a global estimate error (on Omega) is established, interior error estimate on domain omega, overline omega subset estimates on the conditioning related to the trace operator, algorithms of smooth extension. Different numerical experiments in 1D or 2D are implemented. The work is organized as follows:The first chapter recalls theconstruction of multiresolution analysis, important properties of wavelet and numerical algorithms. The second chapter gives an outline of classical fictitious domain method, using Galerkin or Petrov-Galerkin method. We also describe the limitation of this method and point out the direction of our work.\par The third chapter presents a smooth fictitious domain method. It is coupled with Petrov-Galerkin wavelet method for elliptic equations. This section contains the theoretical analysis and numerical implementation to embody the advantages of this new method. The fourth chapter introduces a smooth extension technique. We apply it to elliptic problem with smooth fictitious domain method in 1D and 2D. The fifth chapter is the numerical simulation of the Stefan problem. The property of B-spline render us to exactly calculate the curvature on the moving boundary. We use two examples to test the efficiency of our new method. Then it is used to resolve the two-phase Stefan problem with Gibbs-Thomson boundary condition as an experimental case.