Dans cette thèse, nous nous intéressons à la construction d'éléments finis d'ordre élevé adaptés aux maillages hybrides, pour la résolution de systèmes hyperboliques linéaires en régimes harmonique et temporel. L'accent est plus particulièrement porté sur la construction d'éléments pyramidaux. On étudie trois formulations pour lesquelles on cherche des éléments finis "optimaux" au sens de la convergence dans la norme de l'espace considéré pour la formulation. Pour les formulations H^1 et H(rot), on construit des éléments finis "optimaux" nodaux et hp. Les matrices élémentaires sont évaluées grâce à des formules de quadrature adaptées et des estimations d'erreur sont effectuées pour vérifier la convergence des éléments optimaux construits. Pour la formulation discontinue LDG (Local Discontinuous Galerkin), on présente des éléments utilisant des fonctions de base orthogonales permettant de mettre au point une construction de la matrice de masse et un produit matrice-vecteur rapides. Dans le cas des trois formulations, on étudie les propriétés numériques des éléments construits, on vérifie que l'on retrouve bien numériquement la convergence théorique et on compare nos éléments avec d'autres éléments trouvés dans la littérature. Finalement, on présente des expériences numériques en 3D avec l'équation des ondes ou de Helmholtz, et les équations de Maxwell dans le cas des régimes temporels et harmoniques. On montre ainsi l'efficacité des maillages hybrides par rapport aux maillages purement tétraédriques ou aux maillages hexaédriques obtenus en découpant chaque tétraèdre d'un maillage purement tétraédrique en quatre hexaèdres.

Le sujet principal de cette thèse porte sur l'approximation polynômiale des fonctions aléatoires au moyen de la projection L2 aléatoire discrète, et son application aux problèmes inverses pour les équations aux dérivées partielles avec des données aléatoires. Les motivations proviennent de l'approximation paramétrique de la solution de modèles aux dérivées partielles. La thèse se compose de deux parties, avec un chapitre d'introduction qui résume les techniques modernes de l'approximation polynômiale des fonctions de variables aléatoires. La première partie, du chapitre 1 au chapitre 4, contient l'analyse théorique de la projection L2 aléatoire discrète pour résoudre le problème direct, par exemple, pour rapprocher les moments d'une fonction aléatoire à partir de ses observations, ou pour calculer la solution à un modèle numérique avec des coefficients stochastiques. La stabilité et l'optimalité de l'erreur d'approximation évaluée dans la norme L2 pondérée sont traités. Dans la dernière partie de la thèse, composé des chapitres 5 et 6, la méthodologie développée précédemment pour le problème direct est appliqué aux problèmes inverses pour les équations aux dérivées partielles à coefficients stochastiques. La méthode de factorisation est appliquée dans le cadre de la tomographie par impédance électrique, d'abord dans le cas de coefficient inhomogène, puis dans le cas de coefficient constante par morceaux, à valeurs dans chaque région affectée par l'incertitude. Enfin, dans le chapitre 6 les variantes de la méthode de factorisation proposées dans le chapitre précédent sont accélérés en utilisant les techniques qui ont été présentés dans la première partie de la thèse.

Dans ce mémoire nous allons présenter une méthode de diffusion d’erreur originale qui peut reconstruire des images en demi-ton qui plaisent à l’œil. Cette méthode préserve des détails fins et des structures visuellement identifiables présentes dans l’image originale. Nous allons tout d’abord présenter et analyser quelques travaux précédents afin de montrer certains problèmes principaux du rendu en demi-ton, et nous allons expliquer pourquoi nous avons décidé d’utiliser un algorithme de diffusion d’erreur pour résoudre ces problèmes. Puis nous allons présenter la méthode proposée qui est conceptuellement simple et efficace. L’image originale est analysée, et son contenu fréquentiel est détecté. Les composantes principales du contenu fréquentiel (la fréquence, l’orientation et le contraste) sont utilisées comme des indices dans un tableau de recherche afin de modifier la méthode de diffusion d’erreur standard. Le tableau de recherche est établi dans un étape de pré-calcul et la modification est composée par la modulation de seuil et la variation des coefficients de diffusion. Ensuite le système en entier est calibré de façon à ce que ces images reconstruites soient visuellement proches d’images originales (des aplats d’intensité constante, des aplats contenant des ondes sinusoïdales avec des fréquences, des orientations et des constrastes différents). Finalement nous allons comparer et analyser des résultats obtenus par la méthode proposée et des travaux précédents, et démontrer que la méthode proposée est capable de reconstruire des images en demi-ton de haute qualité (qui préservent des structures) avec un traitement de temps très faible. / In this work we present an original error-diffusion method which produces visually pleasant halftone images while preserving fine details and visually identifiable structures present in original images. We first present and analyze the previous work to show the major problems in halftoning, and explain why we decided to use an error diffusion algorithm to solve the problems. Then we present our method which is conceptually simple and computationally efficient. The source image is analyzed, and its local frequency content is detected. The main components of the frequency content (main frequency, orientation, and contrast) serve as lookup table indices in a pre-computed database of modifications to a standard error diffusion. The modifications comprise threshold modulation and variation of error-diffusion coefficients. The whole system is calibrated in such a way that the produced halftone images are visually close to original images (patches of constant intensity, patches containing sinu- soidal waves of different frequencies/orientations/contrasts, as well as natural images of different origins). Finally, we compare and analyze the results obtained by our method and previous work, and show that our method can produre high-quality halftone image (which is struc- ture aware) within very short time.

Le CEA développe et utilise des logiciels de calcul, également appelés codes de calcul, dans différentes disciplines physiques pour optimiser les coûts de ses installations et de ses expérimentations. Lors d'une étude, plusieurs phénomènes physiques interagissent. Un couplage et des échanges de données entre plusieurs codes sont nécessaires.Chaque code réalise ses calculs sur une géométrie, généralement représentée sous forme d'un maillage contenant des milliers voire des millions de mailles. Cette thèse se focalise sur le transfert de déformations géométriques entre les maillages spécifiques de chacun des codes de calcul couplés. Pour cela, elle présente une méthode de couplage de plusieurs codes, dont le calcul des déformations est réalisé par l'un d'entre eux. Elle traite également de la mise en place d'un modèle commun aux différents codes de l'étude regroupant l'ensemble des données partagées. Enfin, elle porte sur les transferts de déformations entre des maillages représentant une même géométrie ou des géométries adjacentes. Les modifications géométriques sont de nature discrète car elles s'appuient sur un maillage. Afin de les rendre accessible à l'ensemble des codes de l'étude et pour permettre leur transfert, une représentation continue est calculée. Pour cela, deux fonctions sont développées : l'une à support global, l'autre à support local. Toutes deux combinent une méthode de simplification et un réseau de fonctions de base radiale. Un cas d'application complet est traité dans le cadre du réacteur Jules Horowitz. L'effet des dilatations différentielles sur le refroidissement d'un dispositif expérimental est étudié. / The CEA develops and uses scientific software, called physical codes, in various physical disciplines to optimize installation and experimentation costs. During a study, several physical phenomena interact, so a code coupling and some data exchanges between different physical codes are required.Each physical code computes on a particular geometry, usually represented by a mesh composed of thousands to millions of elements. This PhD Thesis focuses on the geometrical modification transfer between specific meshes of each coupled physical code. First, it presents a physical code coupling method where deformations are computed by one of these codes. Next, it discusses the establishment of a model, common to different physical codes, grouping all the shared data. Finally, it covers the deformation transfers between meshes of the same geometry or adjacent geometries. Geometrical modifications are discrete data because they are based on a mesh. In order to permit every code to access deformations and to transfer them, a continuous representation is computed. Two functions are developed, one with a global support, and the other with a local support. Both functions combine a simplification method and a radial basis function network. A whole use case is dedicated to the Jules Horowitz reactor. The effect of differential dilatations on experimental device cooling is studied.

Le but de cette étude est d'améliorer les performances, en termes d'espace mémoire et de temps de calcul, des simulations actuelles de l'Interaction mécanique Pastille-Gaine (IPG), phénomène complexe pouvant avoir lieu lors de fortes montées en puissance dans les réacteurs à eau sous pression. Parmi les méthodes de raffinement de maillage, méthodes permettant de simuler efficacement des singularités locales, une approche multi-grille locale a été choisie car elle présente l'intérêt de pouvoir utiliser le solveur en boîte noire tout en ayant un faible nombre de degrés de liberté à traiter par niveau. La méthode Local Defect Correction (LDC), adaptée à une discrétisation de type éléments finis, a tout d'abord été analysée et vérifiée en élasticité linéaire, sur des configurations issues de l'IPG, car son utilisation en mécanique des solides est peu répandue. Différentes stratégies concernant la mise en oeuvre pratique de l'algorithme multi-niveaux ont également été comparées. La combinaison de la méthode LDC et de l'estimateur d'erreur a posteriori de Zienkiewicz-Zhu, permettant d'automatiser la détection des zones à raffiner, a ensuite été testée. Les performances obtenues sur des cas bidimensionnels et tridimensionnels sont très satisfaisantes, l'algorithme proposé se montrant plus performant que des méthodes de raffinement h-adaptatives. Enfin, l'algorithme a été étendu à des problèmes mécaniques non linéaires. Les questions d'un raffinement espace/temps mais aussi de la transmission des conditions initiales lors du remaillage ont entre autres été abordées. Les premiers résultats obtenus sont encourageants et démontrent l'intérêt de la méthode LDC pour des calculs d'IPG. / The aim of this study is to improve the performances, in terms of memory space and computational time, of the current modelling of the Pellet-Cladding mechanical Interaction (PCI),complex phenomenon which may occurs during high power rises in pressurised water reactors. Among the mesh refinement methods - methods dedicated to efficiently treat local singularities - a local multi-grid approach was selected because it enables the use of a black-box solver while dealing few degrees of freedom at each level. The Local Defect Correction (LDC) method, well suited to a finite element discretisation, was first analysed and checked in linear elasticity, on configurations resulting from the PCI, since its use in solid mechanics is little widespread. Various strategies concerning the implementation of the multilevel algorithm were also compared. Coupling the LDC method with the Zienkiewicz-Zhu a posteriori error estimator in orderto automatically detect the zones to be refined, was then tested. Performances obtained on two-dimensional and three-dimensional cases are very satisfactory, since the algorithm proposed is more efficient than h-adaptive refinement methods. Lastly, the LDC algorithm was extended to nonlinear mechanics. Space/time refinement as well as transmission of the initial conditions during the remeshing step were looked at. The first results obtained are encouraging and show the interest of using the LDC method for PCI modelling.

La compréhension des phénomènes thermiques et mécaniques mis en jeu lors de la mise en forme des matériaux est généralement réalisée avec l’aide de simulations numériques. Ces simulations montrent leurs limites pour les procédés qui conduisent à de très grandes déformations de la matière. Dans ce cas, de très fortes distorsions du maillage se produisent pendant le calcul, entrainant une augmentation de l’erreur, voire l’arrêt prématuré de la simulation. Cette étude porte sur le développement d’une stratégie de remaillage adaptative afin d’éviter les distorsions des éléments pendant les simulations en grandes transformations. La méthode proposée a été intégrée dans un environnement de calcul utilisant le solveur ABAQUS/Explicit, un mailleur 3D et un algorithme de transfert de champ.La méthode h-adaptative en combinaison avec un critère de contrôle basé sur l’endommagement et un estimateur d’erreur de type Zienkiewicz-Zhu Z2 (SPR-amélioré) ont été implantés. Le maillage initial est remplacé par un nouveau maillage avec le niveau de qualité désiré par l’utilisateur, tout en minimisant le nombre des degrés de liberté. Cette technique s’est montrée robuste et entièrement automatique pour déterminer la taille optimale des nouveaux éléments. Une fois le nouveau maillage généré, toutes les variables doivent être soigneusement transférées. Plusieurs techniques de transfert sont décrites et comparées. Des améliorations permettant d’augmenter leurs efficacités en termes de diffusion de l’information et de stabilité numérique ont été proposées. Une attention particulière est portée à la restauration de l'équilibre mécanique local du système. Les différentes techniques développées ont permis de modéliser différents procédés entrainant de grandes déformations élastoplastiques avec endommagement. Dans toutes les applications testées, il a été montré une amélioration de la précision et de la qualité des résultats numériques obtenus. Pour des opérations d’usinage, des mesures de champs cinématiques à travers la technique de corrélation d’images ont été réalisées afin de déterminer les champs de déformation en pointe d’outil. Ces mesures ont servi à la validation de la simulation numérique à l’échelle locale. La comparaison des champs cinématiques expérimentaux avec ceux issus du calcul éléments finis met en évidence la robustesse du processus d’adaptation du maillage proposée pour retranscrire les phénomènes locaux observés expérimentalement. En effet, la reproduction de l’écoulement de la matière sur les bords et la géométrie du copeau sont en très bonne corrélation avec les résultats expérimentaux. Ce développement a permis de proposer une description nouvelle du processus de formation des bandes de cisaillement. / In this work, a fully automated adaptive remeshing strategy, based on a tetrahedral element to simulate various 3D metal forming processes, was proposed. The aim of this work is to solve problems associated with the severe mesh distortion that occurs during the computation and which may be incompatible with the evolution of the physical behavior of the FE solution. Indeed, the quality of the mesh conditions affects the accuracy of the calculations. The proposed strategy is integrated in a computational platform which integrates a finite element solver (Abaqus/Explicit), 3D mesh generation and a field transfer algorithm.The base idea is to use the h-adaptive methodology in the combination with a damage-criterion error and Zienkiewicz-Zhu Z2 type error estimator (SPR-improved) to locally control the mesh modification-as-needed. Once a new mesh is generated, all history-dependent variables need to be carefully transferred between subsequent meshes. Therefore, different transfer techniques are described and compared. An important part of this work concerns the presentation of the proposed modification of the field transfer operator and a special attention is given to restore the local mechanical equilibrium of the system. During the large elasto-plastic deformation simulation with damage, the necessary steps for remeshing the mechanical structure are presented. The several types of applications are also given. For all studied applications, the above strategy can improve the accuracy and quality of numerical results. It also has benefits to decide how refined a mesh needs to be to reach a particular level of accuracy, or how coarse the mesh can be without unacceptably impacting solution accuracy.For the machining processes, kinematic field measurements using Digital image Correlation were performed to validate the numerical simulation at the local level. The comparison of the experimental kinematic fields and those resulting from the FE calculation highlights the robustness of the proposed mesh adaptation process which can transcribe the experimental local phenomena. Also, the reproduction of the material flow at the edges and the chip are correlated with the experimental results accurately. Finally, the physical study of the numerical results can be allowed to propose an innovative description of ASB formation.

Le codage réseau est apparu comme une technique alternative au routage au niveau de la couche réseau permettant d'améliorer le débit et d'optimiser l'utilisation de la capacité du réseau. Récemment, le codage réseau a été appliqué au niveau de la couche physique des réseaux sans-fil pour profiter de la superposition naturelle des signaux effectuée par le lien radio. Le codage réseau peut être vue comme un traitement interne du réseau pour lequel différentes techniques de relayage peuvent être utilisées. Cette thèse étudie un ensemble de traitements ayant des compromis variés en terme de performance et complexité. Nous considérons le canal bidirectionnel à relais, un modèle de canal de communication typique dans les réseaux coopératifs, où deux terminaux s'échangent mutuellement des messages par l'intermédiaire d'un relais. La communication se déroule en deux phases, une phase à accès multiple et une phase de broadcast. Pour ce scénario, nous analysons, dans une première partie, une stratégie de "decode-and-forward". Nous considérons, pour cette étude, des alphabets de taille finie et nous calculons les probabilités moyennes d'erreur de bout-en-bout en se basant sur la métrique d'exposant d'erreur du codage aléatoire. Puis, nous dérivons les régions des débits atteignables par rapport à une probabilité d'erreur maximale tolérable au niveau de chaque nœud. Dans une deuxième partie de la thèse, nous proposons deux schémas de codage réseau pratiques, avec complexité réduite, qui se basent sur la stratégie de relayage "compress-and-forward" (CF). Le premier schéma utilise un codage en réseau de points imbriqués (nested lattices). Le deuxième schéma est une version améliorée qui permet d'atteindre des débits de données supérieurs pour l'utilisateur qui a les meilleures conditions canal. Nous construisons les régions des débits atteignables par les deux schémas proposés tout en optimisant la répartition du temps alloué à chacune des deux phases de transmission. Après l'étude du régime asymptotique, nous analysons le schéma de codage CF avec des réseaux de points de dimension finie. Nous nous concentrons sur le problème de la transmission analogique où la distorsion est optimisée. Enfin, nous étudions l'application d'un schéma de codage, basé sur la stratégie CF avec des réseaux de points imbriqués, pour le canal bidirectionnel à canaux parallèles. Ainsi, nous présentons deux régions de débits atteignables selon la technique de traitement, conjoint ou séparé, des sous-canaux par le relais. / Network coding has emerged as an alternative technique to routing that enhances the throughput at the network layer. Recently, network coding has been applied at the physical layer to take advantage of the natural signal superposition that occurs in the radio link. In this context, the physical-layer network coding can be seen as an in-network processing strategy for which multiple forwarding schemes can be proposed. This thesis investigates a set of processing schemes tailored to the network coding at the physical layer with various compromises between performance and complexity. We consider a two-way relay channel, a typical communication system in cooperative networks, where two terminals communicate with each other via a relay node. This communication occurs during two transmission phases, namely a multiple-access phase and a broadcast phase. For TWRC scenario, we first analyze a decode-and-forward strategy with finite size alphabets. We calculate the end-to-end average error probabilities based on random coding error exponents. Then, we derive the achievable rate regions with respect to a maximal probability of error allowed at each terminal. Next, we propose two low-complexity and practical schemes based on compress-and-forward relaying strategy. The first scheme employs nested lattice coding. The second is an improved version which enables higher data rates for the user experiencing the best channel conditions. We present an information-theoretic framework to reconstruct the achievable rate regions of both schemes by considering optimal time division between both transmission phases. After the asymptotic regime analysis, we study single-layer lattice coding scheme with finite dimension lattices. We focus on the analog transmission problem where the distortion is optimized. Finally, we investigate single-layer lattice coding scheme for parallel Gaussian two-way relay channel. We present two achievable rate regions based on whether the relay processes all the sub-channels jointly or separately.

La modélisation et la simulation numérique des interfaces sont au coeur de nombreuses applications en mécanique des fluides et des solides, telles que la biologie cellulaire (déformation des globules rouges dans le sang), l'ingénierie pétrolière et la sismique (modélisation de réservoirs, présence de failles, propagation des ondes), l'aérospatiale (problème de rupture, de chocs) ou encore le génie civil. Cette thèse porte sur l'approximation des problèmes d'interface et d'interphase en mécanique des fluides par la méthode NXFEM, qui permet de prendre en compte de façon précise une discontinuité non alignée avec le maillage. Nous nous sommes d'abord intéressés au développement de la méthode NXFEM pour des éléments finis non-conformes pour prendre en compte une interface séparant deux milieux. Nous avons proposé deux approches pour les équations de Darcy et de Stokes. La première consiste à modifier les fonctions de base de Crouzeix-Raviart sur les cellules coupées et la deuxième consiste à rajouter des termes de stabilisation sur les arêtes coupées. Les résultats théoriques obtenus ont été ensuite validés numériquement. Par la suite, nous avons étudié la modélisation asymptotique et l'approximation numérique des problèmes d'interphase, faisant apparaître une couche mince. Nous avons considéré d'abord les équations de Darcy en présence d'une faille et, en passant à la limite dans la formulation faible, nous avons obtenu un modèle asymptotique où la faille est décrite par une interface, avec des conditions de transmission adéquates. Pour ce problème limite, nous avons développé une méthode numérique basée sur NXFEM avec éléments finis conformes, consistante et stable. Des tests numériques, incluant une comparaison avec la littérature, ont été réalisés. La modélisation asymptotique a été étendue aux équations de Stokes, pour lesquelles nous avons justifié le modèle limite obtenu. Enfin, nous nous sommes intéressés à la modélisation de la membrane d'un globule rouge par un fluide non-newtonien viscoélastique de Giesekus, afin d'appréhender la rhéologie du sang. Pour un problème d'interphase composé de deux fluides newtoniens (l'extérieur et l'intérieur du globule) et d'un liquide de Giesekus (la membrane du globule), nous avons dérivé formellement le problème limite, dans lequel les équations dans la membrane sont remplacées par des conditions de transmission sur une interface. / Numerical modelling and simulation of interfaces in fluid and solid mechanics are at the heart of many applications, such as cell biology (deformation of red blood cells), petroleum engineering and seismic (reservoir modelling, presence of faults, wave propagation), aerospace and civil engineering etc. This thesis focuses on the approximation of interface and interphase problems in fluid mechanics by means of the NXFEM method, which takes into account discontinuities on non-aligned meshes.We have first focused on the development of NXFEM for nonconforming finite elements in order to take into account the interface between two media. Two approaches have been proposed, for Darcy and Stokes equations. The first approach consists in modifying the basis functions of Crouzeix-Raviart on the cut cells and the second approach consists in adding some stabilization terms on each part of a cut edge. We have studied them from a theoretical and a numerical point of view. Then we have studied the asymptotic modelling and numerical approximation of interphase problems, involving a thin layer between two media. We have first considered the Darcy equations in the presence of a highly permeable fracture. By passing to the limit in the weak formulation, we have obtained an asymptotic model where the 2D fracture is described by an interface with adequate transmission conditions. A numerical method based on NXFEM with conforming finite elements has been developed for this limit problem, and its consistency and uniform stability have been proved. Numerical tests including a comparison with the literature have been presented. The asymptotic modelling has been finally extended to Stokes equations, for which we have justified the limit problem. Finally, we have considered the mechanical behaviour of red blood cells in order to better understand blood rheology. The last part of the thesis is devoted to the modelling of the membrane of a red blood cell by a non-Newtonian viscoelastic liquid, described by the Giesekus model. For an interphase problem composed of two Newtonian fluids (the exterior and the interior of the red blood cell) and a Giesekus liquid (the membrane), we formally derived the limit problem where the equations in the membrane are replaced by transmission conditions on an interface.

Dans un contexte d'adaptabilité ciblée, l'erreur commise sur une quantité d'intérêt peut être représentée grâce aux erreurs globales des problèmes direct et adjoint. Cette représentation de l'erreur est majorée par la somme des indicateurs d'erreurs élémentaires. Ces derniers sont alors utilisés pour produire des raffinements de maillage optimaux. Dans ces travaux, nous proposons de représenter l’erreur du problème adjoint via un opérateur alternatif. L’avantage principal de notre approche est que lorsque l'on choisit correctement l'opérateur alternatif, la majoration correspondante de l'erreur à la quantité d'intérêt devient plus précise, pour autant l'adaptabilité issue de l'utilisation de ces nouveaux indicateurs s'en trouve améliorée. Ces représentations peuvent être employées pour concevoir des algorithmes adaptatifs en espace (h), en ordre d’approximation (p) ou les deux (hp), basés sur la norme d’énergie ou bien ciblés sur une quantité d'intérêt. Bien que la méthode puisse être appliquée à une large gamme de problèmes, nous nous concentrons tout d’abord sur des problèmes unidimensionnels (1D), comme le problème d’Helmholtz et le problème de convection-diffusion stationnaire à convection dominante. Les résultats numériques en 1D montrent que, pour les problèmes de propagation d'ondes, les avantages de notre méthode sont notoires lorsque l'on considère l'opérateur de Laplace pour la représentation de l'erreur. Plus précisément, les majorations issues de la nouvelle représentation sont plus précises que celles provenant de la méthode classique et ce si l'on considère l'énergie globale ou bien une quantité d'intérêt particulière. Le phénomène est d’autant plus notable lorsque l'erreur de dispersion (pollution) est significative. Le problème 1D de convection-diffusion stationnaire à convection dominante avec des conditions limites de Dirichlet homogènes présente une couche limite qui produit une perte de stabilité numérique. La nouvelle représentation d'erreur délivre des majorations plus précises. Lorsqu’appliquée à une p-adaptabilité ciblée, la représentation d'erreur alternative permet une capture plus efficace la couche limite, malgré les oscillations numériques parasites existantes. Devant ces résultats encourageants, nous nous penchons sur l'équation d'Helmholtz à deux et trois dimensions (2D et 3D). Nous montrons, au travers de multiples simulations numériques, que les majorations fournies par les représentations d'erreur alternatives sont plus précises que celle de la représentation classique. Lorsque l'on utilise les indicateurs d'erreur alternatifs, un processus naïf de p-adaptabilité ciblée converge, tandis que dans les mêmes conditions, la méthode classique échoue et requiert l'utilisation d'un opérateur de projection ou d'autre techniques pour récupérer la convergence. Dans ce travail, nous fournissons également des directives pour déterminer les opérateurs qui fournissent des représentations d’erreur induisant de majorations précises. Des résultats similaires sont aussi établis tant pour un problème 2D de convection-diffusion stationnaire à convection dominante que pour des problèmes 2D ayant des coefficients de matériaux discontinus. Nous considérons un problème de diagraphie ultra-sonique en cours de forage pour illustrer l'applicabilité de la méthode proposée. / In Goal-Oriented Adaptivity (GOA), the error in a Quantity of Interest (QoI) is represented using global error functions of the direct and adjoint problems. This error representation is subsequently bounded above by element-wise error indicators that are used to drive optimal refinements. In this work, we propose to replace, in the error representation, the adjoint problem by an alternative operator. The main advantage of the proposed approach is that, when judiciously selecting such alternative operator, the corresponding upper bound of the error representation becomes sharper, leading to a more efficient GOA. These representations can be employed to design novel h, p, and hp energy-norm and goal-oriented adaptive algorithms. While the method can be applied to a variety of problems, in this Dissertation we first focus on one-dimensional (1D) problems, including Helmholtz and steady state convection-dominated diffusion problems. Numerical results in 1D show that for the Helmholtz problem, it is advantageous to select the Laplace operator for the alternative error representation. Specifically, the upper bounds of the new error representation are sharper than the classical ones used in both energy-norm and goal-oriented adaptive methods, especially when the dispersion (pollution) error is significant. The 1D steady state convection-dominated diffusion problem with homogeneous Dirichlet boundary conditions exhibits a boundary layer that produces a loss of numerical stability. The new error representation based on the Laplace operator delivers sharper error upper bounds. When applied to a p-GOA, the alternative error representation captures earlier the boundary layer, despite the existing spurious numerical oscillations. We then focus on the two- and three-dimensional (2D and 3D) Helmholtz equation. We show via extensive numerical experimentation that the upper bounds provided by the alternative error representations are sharper than the classical ones. When using the alternative error indicators, a naive p-adaptive process converges, whereas under the same conditions, the classical method fails and requires the use of the so-called Projection Based Interpolation (PBI) operator or some other technique to regain convergence. We also provide guidelines for finding operators delivering sharp error representation upper bounds. / En un contexto de adaptatividad orientada a un objetivo, el error en una cantidad de interés está representado a través de los errores globales de los problemas directo y adjunto. Esta representación del error se acota superiormente por una suma de indicadores de error de cada elemento. Estos se utilizan para producir refinamientos óptimos. En este trabajo, proponemos representar el error del problema adjunto utilizando un operador alternativo. La principal ventaja de nuestro enfoque es que cuando se elige correctamente dicho operador alternativo, la correspondiente cota superior se vuelve más cercana al error en la cantidad de interés, lo que permite una adaptatividad más eficiente. Estas representaciones pueden ser utilizadas para diseñar algoritmos adaptativos en h, p o hp, basados en la norma de la energía o para aproximar una cantidad de interés específica. Aunque el método propuesto se puede aplicar a una amplia gama de problemas, en esta tesis doctoral nos centramos primero en problemas unidimensionales (1D), tales como el problema de Helmholtz y el problema estacionario de convección-difusión con convección dominante. Los resultados numéricos en 1D muestran que, para los problemas de propagación de ondas, las ventajas de este método son notorias cuando se considera el operador de Laplace para la representación del error. Específicamente, las cotas superiores derivadas de la nueva representación son más cercanas a la cantidad de interés que las del método convencional. Esto es cierto tanto para la norma de la energía global como para una cantidad de interés particular, especialmente cuando el error de dispersión es significativo. El problema estacionario 1D de convección-difusión con convección dominante y con condiciones de Dirichlet homogéneas tiene una capa límite que produce una pérdida de estabilidad numérica. La nueva representación del error proporciona cotas superiores más cercanas a la cantidad de interés. Cuando se aplica a un algoritmo adaptativo en p orientado a un objetivo, la representación alternativa del error captura antes la capa límite, a pesar de las existentes oscilaciones numéricas no físicas. En esta tesis doctoral, también nos centramos en la ecuación de Helmholtz en dos y tres dimensiones (2D y 3D). Mostramos a través de múltiples experimentos numéricos que las cotas superiores proporcionadas por las representaciones alternativas del error son más cercanas a la cantidad de interés que cuando uno considera la representación clásica. Al utilizar los indicadores alternativos del error, un algoritmo adaptativo en p sencillo converge, mientras que en las mismas condiciones, el método convencional falla y requiere el uso de operadores de proyección o de otras técnicas para recuperar la convergencia. En este trabajo, también determinamos operadores que proporcionan representaciones del error que inducen cotas superiores más ajustadas. Establecemos resultados similares tanto para el problema estacionario de convección-difusión con convección dominante en 2D como para problemas 2D con materiales discontinuos. Finalmente, se considera un problema sónico en pozos petrolíferos para ilustrar la aplicabilidad del método propuesto.

La modélisation des problèmes de contact pose de sérieuses difficultés qu'elles soient conceptuelles, mathématiques ou informatiques. Motivés par le rôle fondamental que jouent les phénomènes de contact, nous nous intéressons à la modélisation, l'analyse et la simulation de problèmes de contact intervenant en mécanique des solides et des fluides. Dans une première partie théorique, on étudie le comportement asymptotique de solutions de problèmes variationnels dépendant du temps issus de la mécanique du contact frottant. La deuxième partie est consacrée au contrôle de la qualité des calculs en mécanique des solides. Guidés par la recherche de la formulation et l'étude du contact dans la méthode des éléments finis étendus (XFEM), nous étudions notamment les estimateurs d'erreur par résidu pour la méthode XFEM dans le cas linéaire, ceux pour le problème de contact unilatéral avec frottement de Coulomb approchés par une méthode d'éléments finis standard et l'extension au cas de méthodes mixtes stabilisées (i.e., ne nécessitant pas de condition inf-sup). Cette partie s'achève par la définition du problème de contact avec XFEM suivie d'une estimation a priori de l'erreur. La troisième partie concerne la simulation numérique en mécanique des fluides, plus précisément du problème de contact de la dynamique des globules rouges évoluant dans un fluide régi par les équations de Navier-Stokes en dimension deux.