La détection de peau consiste à détecter les pixels correspondant à une peau humaine dans une image couleur. Les visages constituent une catégorie de stimulus importante par la richesse des informations qu’ils véhiculent car avant de reconnaître n’importe quelle personne il est indispensable de localiser et reconnaître son visage. La plupart des applications liées à la sécurité et à la biométrie reposent sur la détection de régions de peau telles que la détection de visages, le filtrage d'objets 3D pour adultes et la reconnaissance de gestes. En outre, la détection de la saillance des mailles 3D est une phase de prétraitement importante pour de nombreuses applications de vision par ordinateur. La segmentation d'objets 3D basée sur des régions saillantes a été largement utilisée dans de nombreuses applications de vision par ordinateur telles que la correspondance de formes 3D, les alignements d'objets, le lissage de nuages de points 3D, la recherche des images sur le web, l’indexation des images par le contenu, la segmentation de la vidéo et la détection et la reconnaissance de visages. La détection de peau est une tâche très difficile pour différentes raisons liées en général à la variabilité de la forme et la couleur à détecter (teintes différentes d’une personne à une autre, orientation et tailles quelconques, conditions d’éclairage) et surtout pour les images issues du web capturées sous différentes conditions de lumière. Il existe plusieurs approches connues pour la détection de peau : les approches basées sur la géométrie et l’extraction de traits caractéristiques, les approches basées sur le mouvement (la soustraction de l’arrière-plan (SAP), différence entre deux images consécutives, calcul du flot optique) et les approches basées sur la couleur. Dans cette thèse, nous proposons des méthodes d'optimisation numérique pour la détection de régions de couleurs de peaux et de régions saillantes sur des maillages 3D et des nuages de points 3D en utilisant un graphe pondéré. En se basant sur ces méthodes, nous proposons des approches de détection de visage 3D à l'aide de la programmation linéaire et de fouille de données (Data Mining). En outre, nous avons adapté nos méthodes proposées pour résoudre le problème de la simplification des nuages de points 3D et de la correspondance des objets 3D. En plus, nous montrons la robustesse et l’efficacité de nos méthodes proposées à travers de différents résultats expérimentaux réalisés. Enfin, nous montrons la stabilité et la robustesse de nos méthodes par rapport au bruit. / Skin detection involves detecting pixels corresponding to human skin in a color image. The faces constitute a category of stimulus important by the wealth of information that they convey because before recognizing any person it is essential to locate and recognize his face. Most security and biometrics applications rely on the detection of skin regions such as face detection, 3D adult object filtering, and gesture recognition. In addition, saliency detection of 3D mesh is an important pretreatment phase for many computer vision applications. 3D segmentation based on salient regions has been widely used in many computer vision applications such as 3D shape matching, object alignments, 3D point-point smoothing, searching images on the web, image indexing by content, video segmentation and face detection and recognition. The detection of skin is a very difficult task for various reasons generally related to the variability of the shape and the color to be detected (different hues from one person to another, orientation and different sizes, lighting conditions) and especially for images from the web captured under different light conditions. There are several known approaches to skin detection: approaches based on geometry and feature extraction, motion-based approaches (background subtraction (SAP), difference between two consecutive images, optical flow calculation) and color-based approaches. In this thesis, we propose numerical optimization methods for the detection of skins color and salient regions on 3D meshes and 3D point clouds using a weighted graph. Based on these methods, we provide 3D face detection approaches using Linear Programming and Data Mining. In addition, we adapted our proposed methods to solve the problem of simplifying 3D point clouds and matching 3D objects. In addition, we show the robustness and efficiency of our proposed methods through different experimental results. Finally, we show the stability and robustness of our methods with respect to noise.

Les technologies de reconstruction multi-vues permettent de réaliser un clone virtuel d'un acteur à partir d'une simple acquisition vidéo réalisée par un ensemble de caméras à partir de multiples points de vue. Cette approche offre de nouvelles opportunités dans le domaine de la composition de scènes hybrides mélangeant les images réelles et virtuelles. Cette thèse a été réalisée dans le cadre du projet RECOVER 3D dont l'objectif est de développer une chaîne de production TV complète, de l'acquisition jusqu'à la diffusion, autour de la reconstruction multi-vues. Cependant la technologie utilisée dans ce contexte est mal adaptée à la reconstruction de scènes dynamiques. En effet, la performance d'un acteur est reproduite sous la forme d'une séquence d'objets 3D statiques qui correspondent aux poses successives du personnage au cours de la capture vidéo. L'objectif de cette thèse est de développer une méthode pour transformer ces séquences de poses en un modèle animé unique. Les travaux de recherches menés dans ce cadre sont répartis en deux étapes principales. La première a pour but de calculer un champ de déplacements qui décrit les mouvements de l'acteur entre deux poses consécutives. La seconde étape consiste à animer un maillage suivant les trajectoires décrites par le champ de mouvements, de manière à le déplacer vers la pose suivante. En répétant ce processus tout au long la séquence, nous parvenons ainsi à reproduire un maillage animé qui adopte les poses successives de l'acteur. Les résultats obtenus montrent que notre méthode peut générer un modèle temporellement cohérent à partir d'une séquence d'enveloppes visuelles. / 4D multi-view reconstruction technologies are more and more used in media production due to their abilities to produce a virtual clone of an actor from a simple video acquisition performed by a set of multi-viewpoint cameras. This approach is a major advance for the composition of animations which mix virtual and real images, and also offers new possibilities for the rendering of such complex hybrid scenes. The work described in this thesis takes parts in the RECOVER 3D project which aims at developing an innovative industrial framework for TV production, based on multi-view reconstruction, from studio acquisition to broadcasting. The major drawback of the methods used in this context is that they are not adapted to the reconstruction of dynamic scenes. The output are time series which describe the successive poses of the actor, figured as a sequence of static objects. The goal of this thesis is to transform these initial results into a dynamic 3D object where the actor is figured as an animated character. The research detailed in this manuscript presents two main contributions. The first one is centered on the computation of a motion flow which represents the displacements occurring in the reconstructed scene between two consecutive poses. The second one presents a mesh animation process that leads to the animation of a 3D model from one pose to another, following the motion flow. This two-step operation is repeated throughout the entire pose sequence to finally obtain a single animated mesh that matches the evolving shape of the reconstructed actor. Results show that our method is able to produce a temporally consistent mesh animation from various sequences of visual hulls.

Le but de cette thèse est de développer une méthode numérique simple, robuste, efficace et précise pour résoudre des problèmes d'ingénierie de grande taille à partir de la méthode Taylor Meshless (TMM) et fournir de nouvelles idées principales de TMM est d'utiliser comme fonctions de forme des polynômes d'ordre élevé qui sont des solutions approchées de l'EDP. Ainsi la discrétisation ne concerne que la frontière. Les coefficients de ces fonctions de forme sont obtenus en discrétisant les conditions aux limites par des procédures de collocation associées à la méthode des moindres carrés. TMM est alors une véritable méthode sans maillage sans processus d'intégration, les conditions aux limites étant obtenues par collocation. Les principales contributions de cette thèse sont les suivantes: 1) Basé sur TMM, un algorithme général et efficace a été développé pour résoudre des EDP elliptiques tridimensionnelles; 2) Trois techniques de couplage pour des résolutions par morceaux ont été discutées dans des cas de problèmes à grande échelle: la méthode de collocation par les moindres carrés et deux méthodes de couplage basées sur les multiplicateurs de Lagrange; 3) Une méthode numérique générale pour résoudre les EDP non-linéaires a été proposée en combinant la méthode de Newton, la TMM et la technique de différentiation automatique. 4) Pour résoudre des problèmes avec un bord non régulier, des solutions singulières satisfaisant l'équation de contrôle sont introduites comme des fonctions de forme complémentaires, ce qui fournit une base théorique pour la résolution de problèmes singuliers / Based on Taylor Meshless Method (TMM), the aim of this thesis is to develop a simple, robust, efficient and accurate numerical method which is capable of solving large scale engineering problems and to provide a new idea for the follow-up study on meshless methods. To this end, the influence of the key factors in TMM has been studied by solving three-dimensional and non-linear Partial Differential Equations (PDEs). The main idea of TMM is to use high order polynomials as shape functions which are approximated solutions of the PDE and the discretization concerns only the boundary. To solve the unknown coefficients, boundary conditions are accounted by collocation procedures associated with least-square method. TMM that needs only boundary collocation without integration process, is a true meshless method. The main contributions of this thesis are as following: 1) Based on TMM, a general and efficient algorithm has been developed for solving three-dimensional PDEs; 2) Three coupling techniques in piecewise resolutions have been discussed and tested in cases of large-scale problems, including least-square collocation method and two coupling methods based on Lagrange multipliers; 3) A general numerical method for solving non-linear PDEs has been proposed by combining Newton Method, TMM and Automatic Differentiation technique; 4) To apply TMM for solving problems with singularities, the singular solutions satisfying the control equation are introduced as complementary shape functions, which provides a theoretical basis for solving singular problems

Le calcul d'écoulements dans les générateurs de vapeur des réacteurs à eau pressurisée est un problème complexe, faisant intervenir différents régimes d'écoulement et plusieurs échelles de temps et d'espace. Un scénario accidentel peut être caractérisé par des variations très rapides pour un nombre de Mach de l'ordre de l'unité. A l'inverse en régime nominal l'écoulement peut être stationnaire, à bas nombre de Mach. De plus quelque soit le régime considéré, la complexité de la géométrie d'un générateur de vapeur conduit à modéliser le faisceau de tubes par un milieu poreux, d'où le problème de couplage à l'interface avec le milieu fluide.Un schéma de correction de pression tout-Mach en volumes finis colocalisés a été introduit pour les équations d'Euler et de Navier-Stokes. L'existence d'une solution discrète, la consistance du schéma au sens de Lax et la positivité de l'énergie interne ont été démontrées. Le schéma a été ensuite étendu aux modèles diphasiques homogènes du code GENEPI développé au CEA. Enfin un algorithme Multigrille-AMR a été adaptée pour permettre de mettre en oeuvre notre schéma sur des maillages adaptatifs.Concernant la seconde problématique, une extension de la loi de Beavers-Joseph a été proposée pour le régime convectif. En introduisant un saut d'énergie cinétique à l'interface, on retrouve une loi de type Beavers-Joseph mais avec un coefficient de glissement non-linéaire, qui dépend de la vitesse fluide à l'interface et de la vitesse Darcy. La validité de cette nouvelle condition d'interface a été évaluée en réalisant des calculs de simulation numérique directe à différents nombres de Reynolds. / The numerical simulation of steam generators of pressurized water reactors is a complex problem, involving different flow regimes and a wide range of length and time scales. An accidental scenario may be associated with very fast variations of the flow with an important Mach number. In contrast in the nominal regime the flow may be stationary, at low Mach number. Moreover whatever the regime under consideration, the array of U-tubes is modelled by a porous medium in order to avoid taking into account the complex geometry of the steam generator, which entails the issue of the coupling conditions at the interface with the free-fluid.We propose a new pressure-correction scheme for cell-centered finite volumes for solving the compressible Navier-Stokes and Euler equations at all Mach number. The existence of a discrete solution, the consistency of the scheme in the Lax sense and the positivity of the internal energy were proved. Then the scheme was extended to the homogeneous two-phase flow models of the GENEPI code developed at CEA. Lastly a multigrid-AMR algorithm was adapted for using our pressure-correction scheme on adaptive grids.Regarding the second issue addressed in this work, an extension to the Beavers-Joseph law was proposed for the convective regime. By introducing a jump in the kinetic energy at the interface, we recover an interface condition close to the Beavers-Joseph law but with a non-linear slip coefficient, which depends on the free-fluid velocity at the interface and on the Darcy velocity. The validity of this new transmission condition was assessed with direct numerical simulations at different Reynolds numbers.

La France et l’Espagne connaissent depuis plus de vingt ans des processus de recompositions territoriales, dont l’un des enjeux est la création de nouveaux territoires intermédiaires à l’échelle locale (structures intercommunales, “Pays” en France, comarcas en Espagne). Entre commune/municipio et département/provincia, ces derniers sont censés favoriser une refonte de l’organisation territoriale, mais aussi le développement de nouveaux modes de gestion territoriale, plus intégrés et plus horizontaux. À partir des exemples des régions Rhône-Alpes et Catalogne, parmi les premières à développer ces nouveaux espaces locaux, cette thèse interroge leur caractère intermédiaire et le type d’intermédiarité qui les caractérisent. L’analyse révèle la fragilité des mailles intermédiaires en construction, liée à l’échec de leurs potentielles fonctions d’intermédiarité et à la faiblesse née de leur caractère intermédiaire : elles connaissent une neutralisation interne et externe, et peinent donc à trouver leur légitimité par rapport aux niveaux existants et à d’autres nouveaux découpages locaux. La multiplication de ces derniers, qui révèle autant qu’elle explique l’échec des mailles étudiées, montre cependant que c’est moins l’échelle ou le niveau intermédiaire qui sont en cause que les modalités de leur construction. En s’appuyant sur un vaste corpus législatif et sur des enquêtes de terrain (Pays Roannais, comarques du Bages et d’Osona), en croisant de multiples échelles spatio-temporelles dans l’analyse de ces découpages locaux, ce travail invite à systématiser, dans le cadre d’une possible convergence européenne, les études géographiques des systèmes d’administration territoriale. / For more than twenty years, France and Spain have undergone various processes of territorial reshaping, among which the creation of new intermediate authorities at the local scale (inter-city structures, Pays in France, comarcas in Spain). Halfway between commune/municipio and département/provincia, these government units are intended to favor the reform of the territorial organization as well as the design of better integrated and more horizontal territorial management policies. Based on the case studies of pioneer regions Rhône-Alpes and Catalunya, among the first to develop these new local spaces, this work questions the level of intermediateness that characterizes these local units under construction. It shows that the failure of their potential intermediate functions is due to their intermediate nature which itself makes them very fragile. They are subject to both internal and external neutralization and struggle to assert their legitimacy vis à vis existing government levels and other new local territories. The dramatic increase of the latter shows that it is the way the intermediate local units under study are built rather than their mere intermediate scale that explains their failure. Based on a large corpus of legislative material, on field studies (Pays Roannais, comarques of Bages and Osona) and on a multi-level cross-analysis in both space and time, this work aims to systematise the geographical approach of local government systems, a methodology particularly relevant in the perspective of converging local policies at the European level.

Aujourd'hui, avec l'amélioration des puissances de calcul informatique, la simulation numérique est devenue un outil essentiel pour la prédiction des phénomènes physiques et l'optimisation des procédés industriels. La modélisation de ces phénomènes pose des difficultés scientifiques car leur résolution implique des temps de calcul très longs malgré l'utilisation d'importantes ressources informatiques.Dans cette thèse, on s'intéresse à la résolution de problèmes complexes couplant écoulements et transferts thermiques. Les problèmes physiques étant fortement anisotropes, il est nécessaire d'avoir un maillage avec une résolution très élevée pour obtenir un bon niveau de précision. Cela implique de longs temps de calcul. Ainsi il faut trouver un compromis entre précision et efficacité. Le développement de méthodes d'adaptation en temps et en espace est motivé par la volonté de faire des applications réelles et de limiter les inconvénients inhérents aux méthodes de résolution non adaptatives en terme de précision et d'efficacité. La résolution de problèmes multi-échelles instationnaires sur un maillage uniforme avec un nombre de degrés de liberté limité est souvent incapable de capturer les petites échelles, nécessite des temps de calcul longs et peut aboutir à des résultats incorrects. Ces difficultés ont motivé le développement de méthodes de raffinement local avec une meilleure précision aux endroits adéquats. L'adaptation en temps et en espace peut donc être considérée comme une composante essentielle de ces méthodes.L'approche choisie dans cette thèse consiste en l'utilisation de méthodes éléments finis stabilisées et le développement d'outils d'adaptation espace-temps pour améliorer la précision et l'efficacité des simulations numériques.Le développement de la méthode adaptative est basé sur un estimateur d'erreur sur les arrêtes du maillage afin de localiser les régions du domaine de calcul présentant de forts gradients ainsi que les couches limites. Ensuite une métrique décrivant la taille de maille en chaque noeud dans les différentes directions est calculée. Afin d'améliorer l'efficacité des calculs la construction de cette métrique prend en compte un nombre fixe de noeuds et aboutit à une répartition et une orientation optimale des éléments du maillage. Cette approche est étendue à une formulation espace-temps où les maillages et les pas de temps optimaux sont prédits sur des intervalles de temps en vue de contrôler l'erreur d'interpolation sur la domaine de calcul. / Nowadays, with the increase in computational power, numerical modeling has become an intrinsic tool for predicting physical phenomena and developing engineering designs. The modeling of these phenomena poses scientific complexities the resolution of which requires considerable computational resources and long lasting calculations.In this thesis, we are interested in the resolution of complex long time and large scale heat transfer and fluid flow problems. When the physical phenomena exhibit sharp anisotropic features, a good level of accuracy requires a high mesh resolution, hence hindering the efficiency of the simulation. Therefore a compromise between accuracy and efficiency shall be adopted. The development of space and time adaptive adaptation techniques was motivated by the desire to devise realistic configurations and to limit the shortcomings of the traditional non-adaptive resolutions in terms of lack of solution's accuracy and computational efficiency. Indeed, the resolution of unsteady problems with multi-scale features on a prescribed uniform mesh with a limited number of degrees of freedom often fails to capture the fine scale physical features, have excessive computational cost and might produce incorrect results. These difficulties brought forth investigations towards generating meshes with local refinements where higher resolution was needed. Space and time adaptations can thus be regarded as essential ingredients in this recipe.The approach followed in this work consists in applying stabilized finite element methods and the development of space and time adaptive tools to enhance the accuracy and efficiency of the numerical simulations.The derivation process starts with an edge-based error estimation for locating the regions, in the computational domain, presenting sharp gradients, inner and boundary layers. This is followed by the construction of nodal metric tensors that prescribe, at each node in the spatial mesh, mesh sizes and the directions along which these sizes are to be imposed. In order to improve the efficiency of computations, this construction takes into account a fixed number of nodes and generates an optimal distribution and orientation of the mesh elements. The approach is extended to a space-time adaptation framework, whereby optimal meshes and time-step sizes for slabs of time are constructed in the view of controlling the global interpolation error over the computation domain.

Dans le but d’améliorer le comportement thermo-mécanique des aubes de turbine présentes dans les moteurs d’avion développés par SAFRAN, il est nécessaire de mettre au point un nouveau matériau composite. Le procédé de fabrication de ce matériau est complexe et requiert une densification par voie liquide divisée en deux étapes. Cette thèse s’intéresse à la modélisation numérique de la première étape appelée Slurry Cast/APS.Celle-ci consiste en l’injection et le confinement, dans la préforme fibreuse, de particules préalablement mises en suspension. Pour cela, nous avons développé à l’échelle des fibres,un modèle qui utilise les équations de Navier-Stokes incompressibles et monophasiques ,l’équation de Phillips [Phillips et al., 1992] et une loi rhéologique [Krieger, 1972]. Après validation des résultats numériques par comparaison avec des résultats expérimentaux [Hampton et al., 1997] et théoriques [Belfort et al., 1994], le modèle est utilisé pour simuler l’écoulement autour de géométries de tissage proches du matériau étudié. / In order to improve thermo-mechanical behavior of tubine blades in SAFRAN engines plane, a new composite material is necessary. The manufacturing process to obtain this composite is intricate and requires a two steps fluid densification process. This thesis focuses on numerical simulation of the first one called Slurry Cats/APS. In this step, suspended particles are introducted and captured in the reinforcement. For that purpose,we carry out a model at fiber scale, using Navier-Stokes equations in incompressible and monophasic formulation, Phillips equations [Phillips et al., 1992] and a rheological law [Krieger, 1972]. After validation step consisting in a comparison of computational results with experiments [Hampton et al., 1997] and theorical law [Belfort et al., 1994], this model has been used to simulate flow around geometries similar to those encountered in our composite material.

Il est connu depuis des décennies que la propagation des fissures de fatigue dans les matériaux élastoplastiques est très sensible à l’histoire du chargement car le comportement non-linéaire du matériau peut avoir une grande influence sur les vitesses de propagation. Cependant, le calcul brut de millions de cycles de fatigue avec des comportements matériaux non-linéaires sur des structures tridimensionnelles réalistes conduirait à des temps de calcul prohibitifs. Ainsi, nous proposons de coupler deux approches de réduction de modèle a priori et a posteriori, afin de diminuer considérablement le coût de calcul de ce type de problèmes. Tout d’abord, considérant l’hypothèse de plasticité confinée, une stratégie de réduction de modèle a posteriori du comportement plastique de la structure fissurée est proposée. Le modèle réduit ainsi obtenu fournit incrémentalement l’état plastique autour du front de fissure, duquel est déduite la vitesse instantanée de la fissure. De plus, une seconde approche de réduction de modèle, a priori cette fois, est aussi mise en place afin d’accélérer encore plus les temps de résolution du problème global. Cette approche a priori consiste à construire incrémentalement —et sans calculs préalables— une base réduite spécifique à chaque cas-test, en extrayant de l’information des champs de déplacement de la structure au cours du temps et pendant la propagation éventuelle de la fissure. Ainsi, les champs de déplacement solutions de la géométrie fissurée réactualisée sont vus comme une combinaison linéaire de cette base réduite de vecteurs. La méthode numérique considérée ici est la méthode des éléments finis. De fait, pendant la propagation de la fissure, la discrétisation spatiale du modèle doit être réactualisée afin d’être conforme avec le front de la fissure. Dans ce but, une technique spécifique de déformation de maillage est utilisée, et permet de discrétiser la géométrie variable du modèle avec des maillages de même topologie. Cette technique de déformation de maillage apparaît comme une étape clé de la stratégie de réduction de modèle. Finalement, une approche adaptative est construite autour de cette stratégie. Elle permet de garantir la qualité des résultats obtenus par rapport à un critère de précision donné. La précision et l’efficacité de cette stratégie globale sont démontrées à travers de nombreux exemples bidimensionnels et tridimensionnels dans le cadre de propagation de fissure en model, de même que pour un exemple industriel d’une pièce fissurée d’hélicoptère. / It has been known for decades that fatigue crack propagation in elastic-plastic media is very sensitive to load history since the nonlinear behavior of the material can have a great influence on propagation rates. However, the raw computation of millions of fatigue cycles with nonlinear material behavior on tridimensional structures would lead to prohibitive calculation times. In this respect, we propose a global model reduction strategy, mixing both the a posteriori and a priori approaches in order to drastically decrease the computational cost of these types of problems. First, the small scale yielding hypothesis is assumed, and an a posteriori model reduction of the plastic behavior of the cracked structure is performed. This reduced model provides incrementally the plastic state in the vicinity of the crack front, from which the instantaneous crack growth rate is inferred. Then an additional a priori model reduction technique is used to accelerate even more the time to solution of the whole problem. This a priori approach consists in building incrementally and without any previous calculations a reduced basis specific to the considered test-case, by extracting information from the evolving displacement field of the structure. Then the displacement solutions of the updated crack geometries are sought as linear combinations of those few basis vectors. The numerical method chosen for this work is the finite element method. Hence, during the propagation the spatial discretization of the model has to be updated to be consistent with the evolving crack front. For this purpose, a specific mesh morphing technique is used, that enables to discretize the evolving model geometry with meshes of the same topology. This morphing method appears to be a key component of the model reduction strategy. Finally, the whole strategy introduced above is embedded inside an adaptive approach, in order to ensure the quality of the results with respect to a given accuracy. The accuracy and the efficiency of this global strategy have been shown through several examples; either in bidimensional and tridimensional cases for model crack propagation, including the industrial example of a helicopter structure.

L'étude du transfert de chaleur et de masse dans les milieux poreux saturés et insaturés fortement chauffés à leur surface possèdent de nombreuses applications, notamment en archéologie, en agriculture et en géothermie. La première partie de ce travail concerne l'amélioration de la méthode AHC (Accumulation de chaleur latente) qui permet de traiter le changement de phase, dans un milieu homogène : l'intervalle de changement de température au moment du changement de phase apparaît comme un paramètre important, et il doit être choisi proportionnel à la taille des mailles. Des résultats à la fois précis et lisses sont obtenus grâce à un raffinement du maillage localisé près de l'interface de changement de phase. La deuxième partie se rapporte à l'estimation des propriétés thermophysiques du sol par problème inverse à l'aide de données à la fois synthétiques et expérimentales. La méthode de Gauss-Newton avec relaxation et l'algorithme de Levenberg-Marquardt sont utilisés pour résoudre le problème inverse. Le choix de l'intervalle de température de la méthode AHC apparaît crucial : la convergence n'est obtenue parfois qu'au prix d'un enchaînement de plusieurs problèmes inverses. La troisième partie présente un modèle simple pour calculer la conductivité thermique effective d'un milieu granulaire contenant une faible quantité d'eau liquide. La forme exacte de ces ménisques est calculée à l'équilibre. Les résultats montrent un phénomène très net d'hystérésis quand on étudie la variation de la conductivité thermique effective en fonction de la quantité d'eau liquide ; un futur travail concernant un nouveau modèle insaturé, limité au cas du régime pendulaire et présenté à la fin de cette thèse, devrait pouvoir utiliser ces résultats. / In this work, we are interested in studying heat and mass transfer in water saturated and unsaturated porous medium with a strong heating at the surface. Applications concerned are archaeology, agriculture and geothermal engineering. The first part of this work concerns the improvement of the AHC (Apparent Heat Capacity) method used in the numerical resolution of phase change problem in a homogeneous medium: the phase change temperature interval, over which the heat capacity varies, appears as a key parameter which must be chosen proportional to the mesh size. Accurate and smooth results are obtained thanks to a local refinement of the mesh near the phase change interface. The second part is about the estimation of the thermophysical properties of the soil by inverse problem using both synthetic and experimental data. The Damped Gauss-Newton and the Levenberg-Marquardt algorithms are used to solve the problem. In relation with the AHC method, the choice of the phase change temperature interval caused convergence problems which have been fixed by chaining many inverse problems. The obtained results show good convergence to the desired solution. The third part presents a simple model to calculate the effective thermal conductivity of a granular medium which contains a small quantity of liquid water. The exact shape of the liquid menisci between the grains is calculated at equilibrium. The effective thermal conductivity experiences a hysteresis behavior with respect to the liquid volume. A future work that concerns a new unsaturated model, restricted to the pendular regime and detailed at the end of this thesis, should be able to use this result.

Ce travail de thèse s’inscrit dans le cadre de l’étude de l’Interaction mécanique Pastille-Gaine (IPG) se produisant dans les crayons combustibles des réacteurs à eau pressurisée. Ce mémoire porte sur le développement de méthodes de raffinement de maillage permettant de simuler plus précisément le phénomène d’IPG tout en conservant des temps de calcul et un espace mémoire acceptables pour des études industrielles. Une stratégie de raffinement automatique basée sur la combinaison de la méthode multi-grilles Local Defect Correction (LDC) et l’estimateur d’erreur a posteriori de type Zienkiewicz et Zhu est proposée. Cette stratégie s’appuie sur l’erreur fournie par l’estimateur pour détecter les zones à raffiner constituant alors les sous-grilles locales de la méthode LDC. Plusieurs critères d’arrêt sont étudiés afin de permettre de stopper le raffinement quand la solution est suffisamment précise ou lorsque le raffinement n’apporte plus d’amélioration à la solution globale.Les résultats numériques obtenus sur des cas tests 2D élastiques avec discontinuité de chargement permettent d’apprécier l’efficacité de la stratégie proposée.Le raffinement automatique de maillage dans le cas de problèmes de contact unilatéral est ensuite abordé. La stratégie proposée dans ce travail s’étend aisément au raffinement multi-corps à condition d’appliquer l’estimateur d’erreur sur chacun des corps séparément. Un post-traitement est cependant souvent nécessaire pour garantir la conformité des zones de raffinement vis-à-vis des frontières de contact. Une variété de tests numériques de contact entre solides élastiques confirme l’efficacité et la généricité de la stratégie proposée. / This Ph.D. work takes place within the framework of studies on Pellet-Cladding mechanical Interaction (PCI) which occurs in the fuel rods of pressurized water reactor. This manuscript focuses on automatic mesh refinement to simulate more accurately this phenomena while maintaining acceptable computational time and memory space for industrial calculations. An automatic mesh refinement strategy based on the combination of the Local Defect Correction multigrid method (LDC) with the Zienkiewicz and Zhu a posteriori error estimator is proposed. The estimated error is used to detect the zones to be refined, where the local subgrids of the LDC method are generated. Several stopping criteria are studied to end the refinement process when the solution is accurate enough or when the refinement does not improve the global solution accuracy anymore.Numerical results for elastic 2D test cases with pressure discontinuity shows the efficiency of the proposed strategy.The automatic mesh refinement in case of unilateral contact problems is then considered. The strategy previously introduced can be easily adapted to the multibody refinement by estimating solution error on each body separately. Post-processing is often necessary to ensure the conformity of the refined areas regarding the contact boundaries. A variety of numerical experiments with elastic contact (with or without friction, with or without an initial gap) confirms the efficiency and adaptability of the proposed strategy.