L'étude expérimentale présentée vise à étudier l'effet de la température sur le comportement hydraulique et poromécanique d'un calcaire. de nombreux essais (mesures de porosité et de perméabilité, essais de compressions uniaxiales et hydrostatiques) ont été menés afin d'évaluer l'influence des cycles thermiques et donc de la microfissuration thermique sur le comportement de la roche. En parrallèle, un dispositif de mesure de la perméabilité sous températures élevées (jusqu'à 200°C) a été conçu, permettant ainsi d'évaluer directement la variation de la perméabilité du calcaire sous sollicitations thermiques. Par la suite, notre étude a été étendue au cas des matériaux cimentaires, l'objectif étant d'analyser l'effet de la température sur la perméabilité d'un mortier normalisé.

Le travail actuel traite le problème de l'écoulement gaz-liquide compositionnel pour la représentation d'un puits dans des simulations de réservoir. L'objectif est de développer les rapports analytiques qui pourraient lier la pression de puits, la saturation et les concentrations de composant à leurs valeurs moyennes dans chaque zone de l'influence de puits. Nous avons montre que des N-2 équations décrivant le transport des concentrations de phase peuvent être transformées en équations ordinaires (différentiation en ce qui concerne la pression) indépendantes du temps et de l'espace examinant le long des lignes de courant. Ces équations transformées représentent des relations thermodynamiques additionnelles réduisant le degré de liberté thermodynamique. En raison de ceci le variance thermodynamique du modèle compositionnel limite s'avère égal à 1 pour tout nombre de composants chimiques. Cette transformation assurent se découplage total du modèle compositionnel limite dans le nouveau modèle thermodynamique et le modèle hydrodynamique, qui peut être resoved inedpendently d'un un autre. Le modèle thermodynamique décompose est totalement indépendant sur l'hydrodynamique, et décrit le comportement d'équilibre d'un système gaz-liquide ouvert. Ce modèle contient les équations d'équilibre et la EOS classiques, aussi bien que les N-2 nouvelles équations appelées la "Delta-loi", qui déterminent la variation de composition d'un système ouvert dans lequel la masse de chaque composant n'est pas conservée. Le modèle hydrodynamique décompose a été utiliser pour développer les solutions asymptotiques des problèmes d'écoulement de gaz-condensat. Le problème a été montré perturbé singulièrement avec la formation d'une couche limite à voisinage du puits. Dans cette couche la propriété basique de contraste des mobilities de gaz et de liquide est perturbée. On développe une technique spéciale qui permet de construire des expansions asymptotiques sous forme de deux diverses séries: le primer est valide loin du puits (l'expansion extérieure), alors que le second dans valide à voisinage du puits (la couche limite ou l'expansion intérieure). En appliquant la méthode asymptotique suggérée, nous avons développé les solutions asymptotiques pour le problème de l'écoulement multicompositionnel de gaz-condensat àu puits dans un domaine borné à un débit variable. En plusieurs cas la solution peut être obtenue sous la forme analytique, alors qu'en cas général de l'écoulement la méthode mène à une solution de semi-analytical, présentée comme problème initial pour une équation. Cette solution, même étant présenté en forme non-analytique, est beaucoup plus simple que le modèle compositionnel original, car l'équation pour la saturation ne dépend pas de la pression locale, mais dépend de la pression de bord seulement. Dans le dernier chapitre nous avons prolongé cette approche au cas quand la pression capillaire n'est pas négligée. Nous avons supposé cependant que les forces capillaires sont inférieures à la différence de pression entre le puits et la bord de réservoir, dus à ce que nous avons appliqué la méthode de perturbation pour petit nombre capillaire inverse. On obtient les solutions asymptotiques améliorées qui tiennent compte de l'effet capillaire. Simulations numériques montrées que ces effets sont maximaux àu voisinage du puits. Le cas d'une exploitation à long terme du réservoir. Tout d'abord, la simulation traditionnelle du comportement de réservoir peut être effectuée avec l'ECLIPSE en ajoutant la méthode de Peaceman de représentation bonne, qui est une relation analytique pour la pression de puits par l'intermédiaire du débit de production. Cette relation inclut une saturation condensat qui peut être évaluée comme saturation moyenne de réservoir. Une telle simulation fournit un bon résultat pour la pression de puits (ou le débit de production), et un bon résultat pour la saturation de bord, mais des données faibles pour la saturation de puits. Cette valeur peut être calculée en utilisant les solutions asymptotiques suggérées dans le présent projet. Le cas d'un puits de production à court terme (un essai de puits). Il est suffisant de simuler le comportement de réservoir dans le domaine de l'influence de puits, en supposant que l'état de frontière demeure invariable (et connu a priori). Dans ce cas-ci les solutions asymptotiques suggérées dans le travail de presnet peuvent être directement employées pour simuler le problème (sans employer l'ECLIPSE). Le problème de l'écoulement de gaz-condensat à une fracture. Nous avons construit un champ plutôt arbitraire avec des lignes de courant orientées à la fracture, en supposant que la fracture joue le rôle d'une décharge, et les lignes de courant sont stationnaire. Pour une ligne de courant arbitraire nous avons reformulé le modèle d'écoulement de gaz-condensat dans des coordonnées cartésiennes. Pour ce problème nous avons développé les expansions asymptotiques / The present work deals with the problem of the compositional gas-liquid flow for the well representation in reservoir simulations. The objective is to develop analytical relationships which would be able to link the wellbore pressure, saturation and component concentrations to their mean values within each zone of the well influence. It is shown that N-2 equations describing the transport of phase concentrations can be transformed into the space- and time-independent ordinary differential equations (differentiation with respect to pressure) when examined along flow streamlines. These transformed equations represent additional thermodynamic relations reducing the thermodynamic degree of freedom. Due to this the thermodynamic variance of the limit compositional model is shown to be equal to 1 for any number of chemical components. This transformation ensure a total splitting of the limit compositional model into the new thermodynamic model and a hydrodynamic model, which may be resoved inedpendently of one another. The split thermodynamic model is totally independent on the hydrodynamic one, and describes the equilibrium behaviour of an open gas-liquid system. This model contains the classic equilibrium equations and EOS, as well as N-2 new differential equations called the "delta-law" which determine the composition variation in an open system, in which the mass of each component is not conserved. The split hydrodynamic model consists of two equations for pressure and saturation. The split hydrodynamic model was used to develop asymptotic solutions of gas-condensate flow problems. The problem was shown to be singularly perturbed with formation of a boundary layer in the vicinity of the well. In this layer the basic contrast property of gas and liquid mobilities is perturbed. A special technique is developed which enables to construct asymptotic expansions in the form of two various series, one of them is valid far from the well (the exterior expansion), while the second one in valid in the vicinity of the well (the boundary-layer or interior expansion). By applying the suggested asymptotic method, we have developed the asymptotic solutions for the problem of multi-component gas-condensate flow to a well in a bounded domain at a variable flow rate. In several cases the solution may be obtained in the analytical form, while in general case of flow the method leads to a semi-analytical solution presented as an initial problem for a differential equation. This solution, even being presented in non-analytical form, is much simpler than the original compositional model, as the equation for saturation does not depend on the local pressure, but on the boundary pressure only. In the last chapter we extended this approach to the case when the capillary pressure is not neglected. We assumed however that the capillary forces are lower than the pressure difference between the wellbore and reservoir boundary, due to which we applied the perturbation method over the small inverse capillary number. The improved asymptotic solutions are obtained which take into account the capillary effect. Numerical simulations shown that these effects are maximal in the vicinity of the well. For the practice, the obtained asymptotic solutions may be used in the following way to resolve the problem of gas-condensate well representation. The case of a long-term exploitation of the reservoir}. First of all, the traditional simulation of the reservoir behaviour can be performed with ECLIPSE by adding the Peaceman method of well representation, which is an analytical relation for the wellbore pressure via the production rate. This relation includes a condensate saturation which can be evaluated as a mean reservoir saturation. Such a simulation provides a good result for the wellbore pressure (or the production rate), and a good result for the boundary saturation, but poor data for the wellbore saturation. This value can be calculated next by using the asymptotic solutions suggested in the presented project. The case of a short-term well production (a well test). It is sufficient to simulate the reservoir behaviour in the domain of the well influence, by assuming that the boundary state remains invariable (and known a priori). In this case the asymptotic solutions suggested in the presnet work can be directly used to simulate the problem (without using ECLIPSE)

Lorsqu'un fluide complexe s'écoule à travers un milieu poreux, à la non-linéarité de l'écoulement s'ajoute la spécificité de la structure poreuse, qui est souvent multi-échelle. Il émerge alors un grand nombre de problématiques fondamentales liées à l'interaction entre le fluide et la structure poreuse. L'interprétation et la modélisation de la grande variété des phénomènes physiques à petite échelle ainsi que leurs répercussions à grande échelle soulèvent de nombreuses questions. Dans cette thèse, les fluides étudiés sont des solutions de polymères, et les milieux poreux sont, entre autre, des roches issues de réservoirs pétroliers. Dans le contexte des méthodes de récupération améliorée pour les gisements pétroliers, l’injection d’eau avec polymères fait en effet partie des méthodes couramment utilisées, permettant d’augmenter l’efficacité du balayage et donc la production d’huile sur différents types de réservoirs. De part la rhéologie non-Newtonienne ainsi que les phénomènes particuliers proches de la paroi que développent les molécules de polymères, les fluides modélisés dans cette thèse peuvent être qualifiés de complexes. L'objectif de cette thèse est d'étudier la rhéologie non-Newtonienne ainsi que le comportement des molécules de polymère proches de la paroi. On relie alors ces phénomènes aux propriétés effectives macroscopiques. Pour cela, on simule numériquement les écoulements à travers desimages tomographiques de milieux poreux. Ainsi, on souhaite mieux cerner la physique qui est en jeu et également proposer des pistes d'amélioration des modèles actuellement implémentés dans les simulateurs de réservoirs. Dans une première partie de ce travail de thèse, on s'intéresse à la transition du régime d'écoulement macroscopique, de Newtonien à non-Newtonien, induite parune solution de polymère. Par des simulations numériques de l'écoulement à travers un large panel de milieux poreux, on étudie la transition entre ces deux régimes. Une analyse de la mécanique de l'écoulement permet de proposer un modèle simple et d'évaluer en ordre de grandeur la vitesse moyenne de transition. Ensuite, on étudie le glissement apparent induit par un mécanisme de répulsion des chaînes de polymère à la paroi liquide/solide. On propose un modèle à l'échelle des pores de ce phénomène et, par comparaison avec des données expérimentales,on montre que ce modèle permet de retrouver avec une précision acceptable les comportements macroscopiques observés. Enfin, avec des simulations directes sur des milieux périodiques, on relie les phénoménologies micro- et macroscopiques d'écoulements non-Newtoniens. D'un point de vue fondamental, on étudie notamment la compétition entre la non-linéarité induite par la rhéologie non-Newtonienne et le désordre inhérent à la structure poreuse. Les modèles actuellement utilisés dans les simulateurs de réservoirs sont reconsidérés au vu des résultats. / When a complex fluid is flowing through a porous medium, in addition to the fluid intrinsic physics, the multi-scale properties of the porous structure play a significant role. From the interaction between these features arise a great number of complex physical phenomena. The understandingand the modelling of the variety of these phenomena involved at the small scales and their impact on the large scales is the subject of intense work. In this thesis, the fluids we consider are polymer solutions and the porous media are typical sandstones met in petroleum applications. In petroleumengineering, the injection of a polymer slug into the oil-bearing reservoir is indeed a method commonly used in enhanced oil recovery. This method allows to increase the sweep efficiency, hence to improve the oil production of the reservoir. Due to the non-Newtonian rheology induced by the polymer molecules as well as specific mechanisms occurring at the liquid/solid interface, a polymer solution may be qualified as a complex fluid. The goal of this thesis is to investigate the non-Newtonian rheology and the behaviour of the polymer molecules near the liquid/solid interface. These phenomena are related to macro-scale effective properties. To thoroughly address this goal, we perform numerical simulations of flows through porous media. The goal is to obtain a better understanding of the underlying physics and, furthermore, we wish to propose possible improvements of the models that are currently used in reservoir simulators. Primarily, we are interested in the macro-scale transition from a Newtonian to a non-Newtonian flow. This transition is induced by the non-Newtonian rheology. By simulating flows through a wide panel of porous media, we study the macro-scale transition. An analysis of the fluid mechanics involved allows us to propose a simple model for the critical average velocity at which the transition occurs. In addition, we study the apparent slip induced by a repulsive mechanism of the polymer chains from the liquid/solid interface. We propose a pore-scale model to this mechanism. By performing comparisons with experimental datasets, we show that the model allows for a good description of the observed macro-scale behaviour. Finally, with direct simulations over periodic media, we link the phenomenology at the micro- and macro-scale for the flow of non-Newtonian fluids. On a fundamental level, we study the competition that emerges between the nonlinearity induced by the rheology and the disorder inherent to the porous structure. The results are related to modelscommonly used in reservoir simulators.

Les matelas fibreux à base de fibres de verre produits par voie humide ont des applications industrielles en tant que séparateur de batterie et matériaux d’isolation (cœur de panneaux d’isolation sous vide). Ces matériaux utilisent principalement des fibres submicroniques, relativement chères et pouvant présenter des risques pour la santé. Ce projet est une contribution à l'élaboration d'un procédé de fabrication par voie humide de matelas de fibres de verre visant à valoriser des fibres plus grossières, le produit final devant respecter un cahier des charges précis. Nous avons été amenés à étudier le comportement des fibres de verre dans les différentes étapes du procédé et à caractériser les matelas résultants. Nous nous sommes intéressés en particulier au comportement physico-chimique des suspensions aqueuses de fibres de verre. Nous avons caractérisé les propriétés de contexture, la résistance mécanique en traction, la compressibilité et la conductivité thermique des matelas fibreux. Les travaux expérimentaux nous ont permis de proposer une formulation de la composition optimisée et des conditions opératoires du procédé afin que le matelas final soit conforme au cahier des charges. Cette nouvelle composition intègre des fibres de renfort en faible quantité. Elle permet d’améliorer les caractéristiques mécaniques sans affecter les autres propriétés.Enfin, nous avons quantifié les coûts de production et les avons comparés à ceux du procédé actuellement utilisé avec des fibres grossières (production par voie sèche).Mots-clés : Génie papetier, physico-chimie, milieu poreux, fibres de verre, caractérisation. / Glass fibre - based mats produced by a wetlaid process have industrial applications as a battery separator and insulation materials (core of vacuum insulation panels). These materials are mainly made with sub-micron fibres which relatively expensive and can present a risk to health. This project is a contribution to the production of glass fibre-based mats by a wet-laid process to add value to coarser fibres, the final product should respect precise specifications. We have been led to study the behavior of glass fibers in different stages of the process and to characterize the resulting mats. We investigated especially the physico-chemical behavior of aqueous suspensions of glass fiber. We have characterized structure properties, the mechanical resistance to traction, the compressibility and the thermal conductivity of fibrous mats. The experimental work has allowed us to give a formulation of the optimized composition and operational conditions of the process so that the final mattress conforms to the specifications. This new composition includes reinforcement fibres in small quantities. It allows improving the mechanical characteristics without affecting the other properties. Finally, we quantified the production costs and compared them to those of the process currently used with coarse fiber (dry-laid production).Keywords: Paper engineering, physical chemistry, porous media, glass fibre, characterization.

La cristallisation d'un sel initialement dissous dans l'eau contenue dans des matériaux poreux tels que les ciments, béton, mortiers, briques, ..., est une cause majeure de dégradations intervenant dans les constructions et dans certains éléments de notre patrimoine culturel (monuments, murs, sculptures, fresques, ...). Les sels cristallisés à la surface du milieu poreux (efflorescence) n'affectent pas généralement la durabilité de l'ouvrage mais pose un réel problème esthétique pour les façades des bâtiments. Dans le cas d'une cristallisation interne au milieu, les contraintes exercées sur la matrice solide peuvent endommager l'ouvrage allant jusqu'à l'apparition de fissures. Dans ce contexte, le présent travail porte sur l'évaporation avec cristallisation de sel en milieu granulaire et plus précisément sur l'étude de la localisation des cristaux et des déplacements mécaniques induits. La première partie de cette thèse porte sur l'étude des facteurs contrôlant la localisation des cristaux à la surface d'un milieu hétérogène et s'appuie sur l'identification de deux situations de base. Dans la première situation dite de mèche, le milieu poreux est en contact permanent avec une solution saline qui l'alimente à sa base. Cette situation peut correspondre à des parties basses de murs ou d'ouvrages, lorsque le matériau est proche de l'eau des sols. La deuxième situation dite de séchage peut correspondre à des matériaux situés suffisamment en hauteur dans un mur ou un ouvrage. Contrairement au cas de la situation de mèche, le milieu poreux n'est pas alimenté en solution quand l'évaporation se produit. Ces deux situations conduisent à des localisations exactement opposées: à la surface du milieu grossier en séchage, à la surface du milieu fin pour la situation de mèche. Dans le cas de la situation de mèche, la modélisation proposée montre que la structure de l'écoulement joue un rôle clé dans la localisation de la cristallisation alors que pour la situation de séchage, la désaturation préférentielle du milieu grossier est un effet dominant. Concernant la colonisation de la totalité de la surface du milieu poreux par l'efflorescence, l'étude suggère qu'il existe une vitesse d'évaporation critique au regard de laquelle ce phénomène apparait. La deuxième partie de la thèse est dédiée au phénomène de surrection de surface (" surface heave ") à l'aide d'expériences au sein de cellules quasi bidimensionnelles de type Hele Shaw complétés par des expériences en milieu granulaire confiné au sein de tubes cylindriques. Un modèle de croissance de la subflorescence est proposé dans lequel la loi de croissance de la subflorescence est contrôlée par l'évaporation. Ce modèle prédit une très légère sursaturation en haut de la subflorescence compatible avec l'existence d'une pression de cristallisation. Ce modèle permet d'identifier deux régimes principaux de déplacement: un régime dit de désaturation de la subflorescence et un régime de colmatage. Des expressions simples sont proposées pour estimer la croissance de la subflorescence et donc le déplacement induit. Ces relations permettent de déterminer la sursaturation au sommet de la subflorescence et donc la pression de cristallisation. / The crystallization of salt, initially dissolved in water, from porous materials (like cements, concrete, mortar, bricks) is a major cause of degradation occurring in buildings and in certain elements of our cultural landscape (monuments, walls, sculptures, frescoes,). Salt crystallization at the surface of a porous medium (efflorescence) don't affect structure durability but create an aesthetic problem for building construction. However, when salt crystallization occurs inside the porous media (subflorescence), pressure on the solid matrix can damage structures by cracking. In this context, the present work focuses on evaporation with salt crystallization in a granular medium and in particular on the locus of crystals and induced mechanical displacements. The first part of this thesis focuses on the study of the factors controlling the localization of crystals atthe surface of an heterogeneous medium. This study is based on the identification of two basic situations of evaporation. In the first so-called evaporation wicking situation, the porous medium is in contact at its bottom with an aqueous solution and the medium remains fully saturated by the solution during evaporation. This situation may correspond to lower parts of walls or structures, when the material is close to the ground water. The second so-called drying situation can correspond to materials located sufficiently high in a wall or a structure. Contrary to the case of wicking situation, the porous medium isn't supplied with a saline solution when evaporation occurs. These two situations lead to a markedly different locus of the efflorescence formation: on the surface of the coarse medium in drying situation and on the surface of the fine medium for the wicking situation. The study emphasizes the key-role of the velocity field induced in the porous domain in the case of the evaporation-wicking situation. In the case of the drying situation, a key aspect lies in the local increase in the ion mass fraction due to the local desaturation, i.e. the local shrinking of the liquid volume containing the ions. Regarding the colonization of the entire surface of the porous medium by efflorescence, the study suggests that there is a critical evaporation rate against which this phenomenon appears. The second part of the thesis deals with the study of surface heave phenomena using experiments in quasi-two-dimensional Hele-Shaw cell supplemented by experiments in a confined granular medium of cylindrical tubes. A growth model of subflorescence is proposed in which the growing subflorescence is controlled by evaporation and not by the precipitation of salt. This model predicts a very slight supersaturation on the top of the subflorescence compatible with the existence of a crystallization pressure. This model makes it possible to identify two main regimes of displacement: a regime known as desaturation of the subflorescence and a clogging regime. Simple expressions are proposed to estimate the growth of subflorescence and thus the induced displacement. These relationships make it possible to determine the supersaturation at the top of the subflorescence and therefore the crystallization pressure.

Le dimensionnement au séisme, du cœur d'un réacteur à eau pressurisée, est une préoccupation majeure de l'industrie du nucléaire.<br />Nous proposons, dans ce mémoire de thèse, d'établir les équations globales, du comportement du cœur, par une approche milieu poreux. Les équations locales, du fluide et de la structure, sont moyennées sur un volume de contrôle, nous définissons ainsi un fluide équivalent et une structure équivalente, dont les inconnues sont définies sur tout le domaine spatial. Le caractère non linéaire des assemblages combustibles est modélisé par une loi de comportement visco-élastique quadratique. Le couplage fluide-structure est pris en compte par une force volumique dont l'expression est issue de formules empiriques des forces fluides s'exerçant sur un tube soumis à un écoulement axial. Les équations ainsi obtenues sont résolues à l'aide d'une méthode éléments finis.<br />Une validation du modèle est proposée sur trois séries d'essais. La première présente deux assemblages combustibles soumis à un écoulement axial. L'un des deux assemblages est écarté de sa position d'équilibre et lâché, tandis que l'autre est laissé au repos. La deuxième met en œuvre six assemblages en ligne, immergés dans une eau stagnante, posés sur une table vibrante pouvant simuler un séisme. Enfin, la dernière propose neuf assemblages, disposés en un réseau trois par trois, soumis à un écoulement axial. Le déplacement de l'assemblage central est imposé. Les simulations sont en accord avec les expériences ; le modèle reproduit l'influence de la vitesse d'écoulement du fluide sur la dynamique et le couplage des assemblages.

Une technique de renforcement des structures en béton consiste à lier rigidement à ces structures un matériau plus résistant à la traction, par exemple un tissu de fibres de carbone dans une matrice plastique. Le lien mécanique entre le matériau renforçant et la structure se fait actuellement par une colle de type réactive (époxy). Les différentes applications ont montré qu'un facteur limitant de ce procédé était la qualité de l'interface béton/colle à travers laquelle transitent les efforts mécaniques sous forme de contraintes de cisaillement. <br />On propose dans ce travail une modélisation de la phase de pénétration d'une colle fluide (colle à solvant) dans un matériau cimentaire (pâte de ciment). Une première partie consiste à définir les différents phénomènes physiques associés au processus de pénétration puis à établir le système d'équations différentielles décrivant analytiquement ces phénomènes. Une seconde partie décrit les moyens expérimentaux et les sources bibliographiques permettant de disposer des caractéristiques nécessaires à la résolution des équations. Le problème de la colle déposée en surface du matériau cimentaire est enfin étudié par une méthode numérique utilisant les concepts de volumes finis et d'éléments finis. Les résultats obtenus sont présentés et sont comparés avec différentes méthodes de validation expérimentales développées au cours de ce travail.

Le travail présenté dans ce mémoire est consacré à l'application des méthodes de changement d'échelle pour la modélisation du comportement des milieux poreux non saturés dont la phase solide est constituée d'un matériau élastique linéaire. Pour aborder le cas non sturé, il est essentiel de connaître la répartition des phases fluides au sein de l'espace poreux. Les outils permettant de déterminer les positions d'équilibre d'une interface capillaire au sein de d'un milieu poreux et d'étudier leur stabilité sont présentés ici, en particulier dans le cas où il existe des échanges de matière entre les phases fluides. On utilise ensuite les outils de l'homogénéisation linéaire "classique" pour aborder le comportement mécanique d'un milieu poreux saturé par deux phases fluides. Les résultats obtenus sont appliqués pour estimer les déformations de séchage à contrainte imposée constante en condition isotherme en tenant compte de l'influence des caractéristiques morphologiques de l'espace poreux et des domaines occupés par les phases fluides. On propose également une modélisation du comportement d'un milieu mésofissuré non saturé prenant en compte le couplage entre les efforts capillaires et les changements de géométrie des fissures. On montre sur l'exemple du séchage que, dans le cas où les fissures sont très aplaties, la prise en compte de cette non-linéarité modifie radicalement le comportement par rapport à celui prévu par un modèle linéaire. Le mémoire s'achève avec l'examen du comportement d'un matériau granulaire non saturé. On utilise pour cela une approche numérique reposant sur l'hypothèse de périodicité du milieu étudié. L'ensemble des résultats obtenus au cours de ce travail démontre clairement l'intérêt qu'il y a à se référer aux phénomènes se déroulant à l'échelle des pores ainsi qu'aux propriétés morphologiques des domaines occupés par chacune des phases pour construire des modèles de comportement pertinents pour les milieux poreux non saturés.

HOMOGENEISATION NUMERIQUE DE PARAMETRES PETROPHYSIQUES POUR DES MAILLAGES DESTRUCTURES EN SIMULATION DE RESERVOIR : Dans cette thèse on s?intéresse aux méthodes d?homogénéisations qui nécessitent la résolution de problèmes locaux soumis à des conditions aux limites. Dans la première partie on étudie une méthode avec des conditions aux limites linéaires applicables sur des maillages déstructurés. On montre qu?en utilisant une méthode d?éléments finis mixtes [resp. conformes] on obtient une approximation numérique inférieure [resp. supérieure] du tenseur équivalent. On démontre également un résultat de stabilité pour la G-convergence parabolique. Des approximations numériques du tenseur équivalent calculées sur des domaines 2D et 3D déstructurés sont données ainsi que des simulations d?écoulements monophasiques et diphasiques permettant de valider la méthode. Dans la deuxième partie on définit une nouvelle méthode de mise à l?échelle permettant de prendre en compte des conditions aux limites quelconques appliquées aux problèmes locaux. Sur des maillages déstructurés le tenseur équivalent est déterminé en minimisant l?écart des énergies dissipées (ou. des vitesses moyennes) locales et globales. On obtient, en utilisant les techniques du contrôle optimal, un algorithme de calcul effectif qui permet de retrouver, pour les conditions aux limites classiques, les résultats bien connus. Des résultats de convergence et des estimations d?erreurs sont établis. On montre enfin que cette méthode est stable pour la G-convergence et quelques essais numériques 2D sont présentés.

L'humidité dans les murs aggrave les pertes thermiques dans les bâtiments. La présence d'eau est due soit à des remontées capillaires du sol soit à une infiltration de façade en présence d'une pluie battante.Nous tentons dans ce travail une application aux matériaux de construction des méthodes classiques développées pour l'étude des sols. C'est-à-dire une validation des hypothèses et une adaptation des moyens de mesure dans une autre gamme de matériaux a priori plus compacts, de porosité plus faible et plus fine. L'étude sera menée sur un matériau type, a priori bien défini par des normes et facilement reproductible, homogène macroscopiquement, c'est-à-dire au delà d'un volume représentatif relativement faible. Elle portera sur un mortier de ciment. Celui-ci est en effet présent dans un béton où il enrobe les gros granulats; d'autre part, son utilisation est envisagée dans l'isolation thermique des parois sous forme de crépi extérieur ou au sein de revêtements multicouches plus complexes, en particulier dans la rénovation de bâtiments anciens. Dans la définition d'un éventail de moyens de mesure adaptés à ce genre de matériaux, nous serons amenés à présenter la mise au point d'un appareil original de mesure des propriétés conductrices, son application dans divers milieux poreux et son développement vers une nouvelle méthode de mesure de la teneur en eau