Dans ce travail de thèse nous nous intéressons aux propriétés interfaciales du CO2 en conditions de réservoir avec une attention particulière pour le stockage du CO2 dans des aquifères salins profonds. Chaque partie de cette thèse représente une échelle physique différente et même si elles ont été réalisées avec des dispositifs expérimentaux et des outils de modélisation très différents, elles sont très fortement liées. Les résultats obtenus la première partie de cette thèse constituent un jeu de données complet de tension interfaciale saumure-CO2 et représentatif de la gamme de pression, température et salinité des opérations de stockage actuellement en cours d'injection. Une équation semi analytique basée sur ces résultats est proposée avec le but de faciliter le travail des ingénieurs réservoir. La deuxième partie aborde l'étude de ces propriétés à l'échelle du pore en utilisant un micromodèle en verre transparent (milieu poreux 2D) pour différentes conditions de mouillabilité. Cette partie montre le comportement mouillant du CO2 sur des surface hydrophobes ou mouillables à l'huile. Un modèle de type réseau a été utilisé pour l'interprétation et l'exploitation de ces résultats. Le troisième volet correspond à deux phases expérimentales sur carotte pour des mouillabilités différentes, associées à une modélisation à l'échelle de Darcy. Cette partie de notre étude constitue le coeur de la validation du simulateur compositionnel de réservoirs COORES, développé par l'IFP. Elle nous a permis également de déterminer des propriétés polyphasiques, telles que les Pc et kr, pour des systèmes saumure-CO2. L'ensemble de cette thèse présente les trois échelles nécessaires pour modéliser du stockage du CO2 dans des aquifères salins profonds

La modélisation numérique de la migration des hydrocarbures dans les bassins sédimentaires permet de déterminer les accumulations d’hydrocarbures au sein des formations géologiques. A partir de cela on peut prévoir la hauteur d’hydrocarbure piégé. Cette détermination est essentielle dans l’industrie du pétrole. Cependant grâce à certaines études, on a pu s’apercevoir que des erreurs numériques pouvaient apparaître lors de l’utilisation de pression capillaires polynomiales. Dans cette thèse, nous travaillons principalement sur les modèles dits de Darcy et nous évoquons sommairement les modèles de type percolation. L’objectif de ce travail est de fournir de nouveaux modèles des pressions capillaires, donnant de meilleurs résultats que ceux actuellement utilisés. Dans un premier temps, nous décrivons les mécanismes de la migration des hydrocarbures dans les couches. Ensuite nous étudions plus attentivement les lois des pressions capillaires permettant l’écoulement des fluides. Puis nous établissons la discrétisation, suivant la méthode des volumes finis, du problème. Dans la partie suivante nous testons en une dimension de nouveaux modèles de pressions capillaires affines par morceaux. Puis dans une autre partie, nous faisons les tests en deux dimensions de ces modèles auxquels nous rajoutons un autre modèle, bâti à partir des deux précédents. En conclusion, nous synthétisons l’ensemble des résultats et évoquons certaines perspectives concernant l’amélioration des modèles étudiés. / Numerical modeling of hydrocarbon migration in sedimentary basins makes it possible to determine hydrocarbon accumulations within geological formations. From this it is possible to predict the trapped hydrocarbon height. This determination is essential in the petroleum industry. However, thanks to some studies, it has been found that numerical errors can occur when using polynomial capillary pressure. In this thesis, we work mainly on the so-called models of Darcy and we briefly discuss percolation-type models. The objective of this work is to provide new models of capillary pressures, giving better results than those currently used. First, we describe the mechanisms of hydrocarbon migration in the layers. Then we study more closely the laws of capillary pressures allowing the flow of fluids. Then we establish the discretization, according to the finite volume method, of the problem. In the next part we test in one dimension new models of capillary pressures affine in pieces. Then in another part, we do the two-dimensional tests of these models to which we add another model, built from the two previous ones. In conclusion, we summarize all the results and discuss some perspectives concerning the improvement of the studied models.

Les travaux de cette thèse portent sur des méthodes de volumes finis sur maillages quelconque pour la discrétisation de problèmes d'évolution non linéaires modélisant le transport de contaminants en milieu poreux et les écoulements diphasiques.Au Chapitre 1, nous étudions une famille de schémas numériques pour la discrétisation d'une équation parabolique dégénérée de convection-reaction-diffusion modélisant le transport de contaminants dans un milieu poreux qui peut être hétérogène et anisotrope. La discrétisation du terme de diffusion est basée sur une famille de méthodes qui regroupe les schémas de volumes finis hybrides, de différences finies mimétiques et de volumes finis mixtes. Le terme de convection est traité à l'aide d'une famille de méthodes qui s'appuient sur les inconnues hybrides associées aux interfaces du maillage. Cette famille contient à la fois les schémas centré et amont. Les schémas que nous étudions permettent une discrétisation localement conservative des termes d'ordre un et d'ordre deux sur des maillages arbitraires en dimensions d'espace deux et trois. Nous démontrons qu'il existe une solution unique du problème discret qui converge vers la solution du problème continu et nous présentons des résultats numériques en dimensions d'espace deux et trois, en nous appuyant sur des maillages adaptatifs.Au Chapitre 2, nous proposons un schéma de volumes finis hybrides pour la discrétisation d'un problème d'écoulement diphasique incompressible et immiscible en milieu poreux. On suppose que ce problème a la forme d'une équation parabolique dégénérée de convection-diffusion en saturation couplée à une équation uniformément elliptique en pression. On considère un schéma implicite en temps, où les flux diffusifs sont discrétisés par la méthode des volumes finis hybride, ce qui permet de pouvoir traiter le cas d'un tenseur de perméabilité anisotrope et hétérogène sur un maillage très général, et l'on s'appuie sur un schéma de Godunov pour la discrétisation des flux convectifs, qui peuvent être non monotones et discontinus par rapport aux variables spatiales. On démontre l'existence d'une solution discrète, dont une sous-suite converge vers une solution faible du problème continu. On présente finalement des cas test bidimensionnels.Le Chapitre 3 porte sur un problème d'écoulement diphasique, dans lequel la courbe de pression capillaire admet des discontinuité spatiales. Plus précisément on suppose que l'écoulement prend place dans deux régions du sol aux propriétés très différentes, et l'on suppose que la loi de pression capillaire est discontinue en espace à la frontière entre les deux régions, si bien que la saturation de l'huile et la pression globale sont discontinues à travers cette frontière avec des conditions de raccord non linéaires à l'interface. On discrétise le problème à l'aide d'un schéma, qui coïncide avec un schéma de volumes finis standard dans chacune des deux régions, et on démontre la convergence d'une solution approchée vers une solution faible du problème continu. Les test numériques présentés à la fin du chapitre montrent que le schéma permet de reproduire le phénomène de piégeage de la phase huile.

Dans une démarche de valorisation des déchets, les granulats recyclés sont introduits dans la formulation des bétons pour donner naissance à de nouveaux bétons recyclés qui représentent l’objet du projet national RECYBETON et du projet ANR ECOREB. Cette thèse se concentre sur l’étude du béton frais et en particulier la maîtrise du retrait plastique et l’effet du degré de saturation initial des gravillons recyclés sur le comportement des bétons recyclés. Deux familles de bétons avec des rapports eau/ciment respectifs de 0,6 et 0,45 ont été testés en conditions endogènes ainsi qu’en dessiccation (Vvent = 8 m/s). Chaque famille de bétons est constituée d’un béton naturel de référence et de deux bétons recyclés différenciés pas le degré initial de saturation des gravillons recyclés (50 et 120 % de l’absorption nominale). Les résultats expérimentaux soulignent la capacité des gravillons recyclés initialement partiellement saturés à capter rapidement l’eau contenue dans la pâte de ciment, modifiant ainsi le rapport E/C, les propriétés rhéologiques du béton frais et les résistances mécaniques du béton durci. Après saturation en eau, si les conditions de séchage conduisent à un manque d’eau dans le béton, les gravillons recyclés peuvent fournir de l’eau à la pâte. Ils constituent donc un potentiel de cure interne. Le retrait plastique sous vent est explicitement lié au ressuage, au développement de la pression capillaire et à la fissuration. Le temps d’initiation de la fissuration dépend de la quantité d’eau totale dans le béton et de sa capacité à ressuer tandis que l’ouverture de fissure varie avec la valeur de retrait plastique mesurée. Dans cette étude, le développement de la pression capillaire est la cause de la fissuration qui apparaît dès l’entrée d’air dans le matériau poreux. Les différences de comportements les plus importantes sont observées entre bétons ayant une quantité d’eau totale différente plutôt qu’entre bétons naturel et recyclé. La compilation des résultats expérimentaux a permis de mettre sur pied des modélisations qui illustrent les comportements observés. Les pores concernés par l’entrée d’air dans les bétons recyclés et naturels au moment de la fissuration sont les plus gros pores de la pâte. Enfin, un couplage hygrothermique séchage-température du béton peut influer sur le démarrage de l’hydratation. / In the context of sustainable development, the reuse of construction and demolition waste is necessary to conserve nonrenewable natural aggregate resources, so recycled aggregates are introduced in concrete mix design. This is the aim of the national projet RECYBETON and the research project ECOREB. This study deals with the fresh concrete and more specifically with shrinkage control and the effects of the initial saturation degrees of recycled coarse aggregates on concrete behavior.Two concrete families, with two different water/cement ratios 0,60 and 0,45, are tested under endogenous and drying (wind speed equal to 8 m/s) conditions. Each concrete family contains a reference natural concrete and two recycled concretes. The initial saturation degree is the difference between them (recycled coarse aggregates saturated or semi saturated).Experimental results underline the capacity of non-saturated aggregates to quickly absorb water from cement paste, modifying the W/C ratio, rheological properties of the fresh concrete and the mechanical strength (at 28 days) of recycled concretes. After saturation in water, recycled aggregates can release water into the cement paste if the undergone drying conditions lead to a lack of water in the cement matrix. The recycled coarse aggregates can be seen as an internal curing potential.Experimental plastic shrinkage studies carried out under drying conditions highlight a link between bleeding, capillary pressure, plastic shrinkage and cracking. It should be pointed out that the initial cracking is dependent on the total quantity of water in the concrete and on its bleeding capacity. The opening cracks vary with the plastic shrinkage values measured during the test. The analysis of the results emphasize that the capillary pressure is the determining parameter and that the air entry value matches the cracks. The major behavior differences are found between concretes with different volumes of water rather than between natural and recycled concretes.Finally, the analysis of all the experimental results have allowed concrete modelling and understanding why concretes do not behave in the same way. When it cracks, the air come in the biggest pores of the concrete paste. Moreover, a hygrothermal coupling exists between the drying and the temperature in concrete. It can affect hydration start up.

Dans cette thèse, nous nous sommes intéressés à la mesure et à la modélisation des propriétés d'écoulements diphasiques (perméabilités relatives) dans des milieux poreux carbonatés représentatifs de la zone de transition de réservoirs pétroliers carbonatés. Notre étude expérimentale a porté sur deux roches carbonatées de même minéralogie mais de propriétés pétrophysiques différentes (unimodale et bimodale). Les fluides utilisés sont une saumure, une huile minérale pour les expériences en mouillabilité franche à l'eau et une huile brute pour les expériences en mouillabilité altérée (obtenue par macération). Durant les expériences en « écoulements permanents », la perte de charge est mesurée continûment et le champ de saturation locale est mesuré en fin d'expérience par atténuation gamma. Les principaux résultats expérimentaux de cette étude sont les suivantes :- La mouillabilité dépend de la saturation initiale en huile et donc de la hauteur de la zone de transition. Elle a une influence sur les perméabilités relatives à l'huile ET à l'eau.- L'hystérésis des perméabilités relatives à l'eau est d'autant plus prononcée que la saturation initiale en huile augmente. Une observation similaire est faite sur la courbure des perméabilités relatives à l'huile.- Les évolutions des Sorw en fonction des Soi ne vérifient pas toutes, une loi monotone de type Land. Leur évolution dépend aussi de la structure poreuse.- Les deux types de roches carbonatés présentent une altération de mouillabilité différente, la structure vacuolaire favorisant l'altération de mouillabilité. Nos données expérimentales ont été interprétées en utilisant différents modèles d'hystérésis de la littérature. Il apparait clairement des différences entre les prédictions numériques et nos données expérimentales car ces modèles ne prennent pas en compte l'évolution de la mouillabilité. Nous avons proposé un modèle de calcul des perméabilités relatives qui décrit nos données expérimentales de manière satisfaisante. / In this work, we take an interest in measuring and modeling the multiphase flow properties (relative permeabilities) in different carbonate rock types that are representative of the transition zone of carbonate reservoirs. Our experimental study focused on two different carbonate rock types, with almost the same mineralogy, but different petrophysical properties (bimodal and unimodal). The system of fluids used is composed of brine, a mineral oil (for experiments at water-wet conditions), and a crude oil (for experiments including wettability alteration after an ageing process at high temperature). During the experiments performed by the steady-state core flooding method, the pressure drop is continuously measured along the core, and the saturation profiles are obtained at the end of the flooding, by gamma ray measurements. The main results of this study are as follows:- Wettability depends on the initial oil saturation and hence the heights above the oil-water contact. It affects the oil relative permeabilities AND water relative permeabilities- The hysteresis observed on water relative permeability is even more pronounced that the initial oil saturation increases. A similar observation is made on the curvatures of the oil relative permeability.- The evolution of Sorw as function of Soi does not always follow a classical Land correlation. Their evolution also depends on the pore structure.- The two rock types exhibit different responses to wettability alteration. The vuggy structure might promote oil-wetness.Our experimental data were interpreted using different hysteretic models of the literature. We observed significant differences between the results predicted by the models and our experimental observations which are attributed mainly to the failures to take into account the change in wettability. We proposed an improved hysteresis model, which described our experimental data satisfactorily.

Parmi les options en cours d’investigation, le stockage souterrain de l'hydrogène dans les formations sédimentaires comme les grès pourrait offrir un potentiel unique pour stocker de grandes quantités d'énergie. L'évaluation des modalités de stockage souterrain de l'hydrogène nécessite donc à la fois une connaissance précise des transformations minéralogiques dues à la présence de l'hydrogène et l’acquisition de données sur le comportement hydrodynamique des fluides. Par conséquent, cette étude se composera de trois parties : 1- Etude des interactions géochimiques de l’hydrogène dans des formations sédimentaires gréseuses : Les produits expérimentaux portent la marque d'une réaction très limitée entre les minéraux du grès et l'hydrogène. Si les résultats expérimentaux sont combinés aux résultats numériques, l’étude démontre que l'hydrogène, une fois injecté, peut être considéré comme relativement inerte. De façon globale, nos résultats renforcent la faisabilité du confinement de l'hydrogène dans des réservoirs géologiques comme les grès. 2- Etude de la migration de l'hydrogène dans les grès : détermination de la perméabilité relative et de la pression capillaire du système hydrogène-eau : Afin de fournir des données quantitatives pour le développement du stockage souterrain de l'hydrogène, la pression capillaire et la perméabilité relative ont été mesurées pour le système hydrogène-eau en deux conditions potentielles. Les résultats indiquent que les données obtenues sont applicables à l’ensemble des conditions de stockage de l'hydrogène. 3- Modélisation numérique d’un site de stockage géologique d’hydrogène : La simulation numérique a été effectuée pour caractériser l'évolution dynamique d’un site de stockage d'hydrogène pur. Une fluctuation saisonnière du fonctionnement du réservoir et l'effet des fuites d'hydrogène dus aux réactions ont été pris en compte. / Underground hydrogen storage has been introduced as storage solution for renewable energy systems as it offers a unique potential to store large amounts of energy, especially in sedimentary formations such as sandstones. However, evaluating the underground hydrogen storage requires a precise knowledge of the hydrodynamic behavior of the fluids and of mineralogical transformations due to the presence of hydrogen that may affect the storage properties. Therefore, this study is consists in three parts: 1- Study of geochemical reactivity of hydrogen in sandstone sedimentary formations: The experimental products bear the mark of only very limited reaction between sandstone minerals and hydrogen. Taken together with the numerical results, this study demonstrates that hydrogen, once injected, can be considered as relatively inert. Overall, our results support the feasibility of hydrogen confinement in geological reservoirs such as sandstones. 2- Study of the migration of hydrogen in sandstone: determination of relative permeability and capillary pressure of hydrogen-water system: To provide quantitative data for the development of underground hydrogen storage, capillary pressures and relative permeabilities of hydrogen-water system have been measured at two potential conditions. The interpretation of the results would suggest that the obtained data are applicable for the entire range of hydrogen storage conditions. Interfacial tensions and contact angles for the hydrogen-water system have been also derived. 3- Numerical simulation of a geological hydrogen storage site: The numerical simulation was performed to characterize the evolution of pure hydrogen storage, by considering the seasonal fluctuation of renewable energy and the effect of hydrogen loses due to the biotic reactions.

Ce mémoire est consacré à l'analyse des effets du confinement de l'air et du caractère multidimensionnel de l'écoulement lors des essais de caractérisation hydrodynamique d'un sol par la méthode du drainage interne. L'étude expérimentale et numérique de la compression de l'air au cours de l'infiltration montre que les caractéristiques du sol ne sont pas affectées par l'écoulement de l'air, dès lors que la pression capillaire est mesurée et non celle de l'eau uniquement. L'étude numérique du caractère axisymétrique de l'écoulement lors d'un essai à l'infiltromètre double anneaux montre que malgré l'importance du transfert radial, son influence sur la caractérisation d'un milieu homogène reste faible. En revanche, la caractérisation de milieux stratifiés horizontalement est très délicate.

Une solution prometteuse pour diminuer les émissions anthropogéniques de gaz à effet de serre consiste à injecter une partie des rejets industriels de CO2 dans des formations souterraines. Celles-ci comportent un réservoir entouré de roches de couverture, qui constituent la première barrière à la migration des fluides. La caractérisation de leurs propriétés de confinement et de leur évolution en présence de CO2 est donc un élément clé de la sécurité d’un site de stockage. Le travail présenté propose une méthodologie, appliquée ici à des roches carbonatées du bassin parisien, permettant de mesurer les paramètres de transport de roches de couverture et les conséquences d’un vieillissement en conditions représentatives de celles d’un stockage en aquifère profond. La pression de percée, le coefficient de diffusion des produits de dissolution du CO2, et la perméabilité, paramètres contrôlant les principaux mécanismes de fuite, ont été mesurés avant et après altération des matériaux par réaction avec une saumure saturée en CO2 dans des conditions thermodynamiques typiques d’un réservoir (environ 80°C et 100 bar). Les résultats obtenus ont révélé un bon comportement global des roches, mais également une forte diminution du potentiel de confinement en présence de défauts structurels initiaux (fractures rebouchées, pores de grand diamètre...). Une simulation numérique décrivant les évolutions de la formation rocheuse non-fissurée sur une durée de 1000 ans a été réalisée en s’appuyant sur des paramètres mesurés directement ou obtenus par modélisation des essais d’altération. Elle a montré que les transformations engendrées par le stockage de CO2 sous une roche de couverture homogène restent très limitées spatialement / A promising solution to reduce anthropogenic emissions of greenhouse effect gases consists in the injection and long-term storage of a part of the industrial carbon dioxide discharges in underground formations. These formations must be composed of a reservoir surrounded by tight caprocks, which represent the first barrier preventing fluids migration. The characterization of their confining properties and of their evolution in presence of CO2 is thus a key element regarding a storage site security. This work presents a methodology allowing the measurement of caprocks transport parameters and the consequences of an alteration under representative conditions of deep aquifers storage. This methodology was applied to carbonate rocks from the Paris basin. The breakthrough pressure, the diffusion coefficient of CO2 dissolution products, and the permeability, controlling parameters of leakage mechanisms, were measured before and after alteration of the materials by reaction with a CO2-saturated brine under reservoir thermodynamic conditions (about 80°C and 100 bar). Results revealed a satisfactory global behaviour under these aggressive conditions, but also a strong diminution of the confinement potential in presence of initial structural faults (sealed fractures, large-diameter pores…) forming higher-permeability zones. A numeric simulation describing the evolution of a homogeneous rock formation during 1000 years was also realized based on parameters directly measured or obtained by modelling of the alteration experiments. It showed that the transformations brought by the CO2 storage under a rock formation with no initial faults remain very localized spatially

Les travaux de cette thèse portent sur des méthodes de volumes finis sur maillages quelconque pour la discrétisation de problèmes d'évolution non linéaires modélisant le transport de contaminants en milieu poreux et les écoulements diphasiques.Au Chapitre 1, nous étudions une famille de schémas numériques pour la discrétisation d'une équation parabolique dégénérée de convection-reaction-diffusion modélisant le transport de contaminants dans un milieu poreux qui peut être hétérogène et anisotrope. La discrétisation du terme de diffusion est basée sur une famille de méthodes qui regroupe les schémas de volumes finis hybrides, de différences finies mimétiques et de volumes finis mixtes. Le terme de convection est traité à l'aide d'une famille de méthodes qui s'appuient sur les inconnues hybrides associées aux interfaces du maillage. Cette famille contient à la fois les schémas centré et amont. Les schémas que nous étudions permettent une discrétisation localement conservative des termes d'ordre un et d'ordre deux sur des maillages arbitraires en dimensions d'espace deux et trois. Nous démontrons qu'il existe une solution unique du problème discret qui converge vers la solution du problème continu et nous présentons des résultats numériques en dimensions d'espace deux et trois, en nous appuyant sur des maillages adaptatifs.Au Chapitre 2, nous proposons un schéma de volumes finis hybrides pour la discrétisation d'un problème d'écoulement diphasique incompressible et immiscible en milieu poreux. On suppose que ce problème a la forme d'une équation parabolique dégénérée de convection-diffusion en saturation couplée à une équation uniformément elliptique en pression. On considère un schéma implicite en temps, où les flux diffusifs sont discrétisés par la méthode des volumes finis hybride, ce qui permet de pouvoir traiter le cas d'un tenseur de perméabilité anisotrope et hétérogène sur un maillage très général, et l'on s'appuie sur un schéma de Godunov pour la discrétisation des flux convectifs, qui peuvent être non monotones et discontinus par rapport aux variables spatiales. On démontre l'existence d'une solution discrète, dont une sous-suite converge vers une solution faible du problème continu. On présente finalement des cas test bidimensionnels.Le Chapitre 3 porte sur un problème d'écoulement diphasique, dans lequel la courbe de pression capillaire admet des discontinuité spatiales. Plus précisément on suppose que l'écoulement prend place dans deux régions du sol aux propriétés très différentes, et l'on suppose que la loi de pression capillaire est discontinue en espace à la frontière entre les deux régions, si bien que la saturation de l'huile et la pression globale sont discontinues à travers cette frontière avec des conditions de raccord non linéaires à l'interface. On discrétise le problème à l'aide d'un schéma, qui coïncide avec un schéma de volumes finis standard dans chacune des deux régions, et on démontre la convergence d'une solution approchée vers une solution faible du problème continu. Les test numériques présentés à la fin du chapitre montrent que le schéma permet de reproduire le phénomène de piégeage de la phase huile. / In Chapter 1 we study a family of finite volume schemes for the numerical solution of degenerate parabolic convection-reaction-diffusion equations modeling contaminant transport in porous media. The discretization of possibly anisotropic and heterogeneous diffusion terms is based upon a family of numerical schemes, which include the hybrid finite volume scheme, the mimetic finite difference scheme and the mixed finite volume scheme. One discretizes the convection term by means of a family of schemes which makes use of the discrete unknowns associated to the mesh interfaces, and contains as special cases an upwind scheme and a centered scheme. The numerical schemes which we study are locally conservative and allow computations on general multi-dimensional meshes. We prove that the unique discrete solution converges to the unique weak solution of the continuous problem. We also investigate the solvability of the linearized problem obtained during Newton iterations. Finally we present a number of numerical results in space dimensions two and three using nonconforming adaptive meshes and show experimental orders of convergence for upwind and centered discretizations of the convection term.In Chapter 2 we propose a finite volume method on general meshes for the numerical simulation of an incompressible and immiscible two-phase flow in porous media. We consider the case that it can be written as a coupled system involving a degenerate parabolic convection-diffusion equation for the saturation together with a uniformly elliptic equation for the global pressure. The numerical scheme, which is implicit in time, allows computations in the case of a heterogeneous and anisotropic permeability tensor. The convective fluxes, which are non monotone with respect to the unknown saturation and discontinuous with respect to the space variables, are discretized by means of a special Godunov scheme. We prove the existence of a discrete solution which converges, along a subsequence, to a solution of the continuous problem. We present a number of numerical results in space dimension two, which confirm the efficiency of the numerical method.Chapter 3 is devoted to the study of a two-phase flow problem in the case that the capillary pressure curve is discontinuous with respect to the space variable. More precisely we assume that the porous medium is composed of two different rocks, so that the capillary pressure is discontinuous across the interface between the rocks. As a consequence the oil saturation and the global pressure are discontinuous across the interface with nonlinear transmission conditions. We discretize the problem by means of a numerical scheme which reduces to a standard finite volume scheme in each sub-domain and prove the convergence of a sequence of approximate solutions towards a weak solution of the continuous problem. The numerical tests show that the scheme can reproduce the oil trapping phenomenon.

L'objectif de cette recherche est l'étude du transfert d'eau en phase liquide dans les parois du bâtiment et de définir une nouvelle métrologie pour déterminer le coefficient de perméabilité liquide des matériaux du bâtiment. Une partie théorique fait l'état des connaissances en matière de modélisation du transfert d'eau en phase liquide. Elle est suivie par une partie qui relate les méthodes utilisées pour mesurer ou imposer des pressions capillaires (ou succions) avec un intérêt particulier à la méthode osmotique qui est utilisée plus tard pour imposer des conditions aux limites en terme de succion. Cette dernière partie s'ouvre sur une revue des méthodes de mesure de la perméabilité liquide et sera conclue par la définition et la mise au point d'une nouvelle méthode de mesure de la perméabilité liquide, qui a fait l'objet d'un dépôt de brevet d'invention, s'appuyant sur la méthode osmotique et qui fonctionne en régime permanent, mesurant les débits entrant et sortant de l'échantillon. Une campagne d'essais est menée sur du béton cellulaire autoclave avec différents couples de conditions aux limites. Il s'ensuit une recherche d'identification d'une relation entre la perméabilité liquide et la succion, démarche qui mérite plus d'approfondissement. Finalement cette nouvelle méthode de mesure mérite d'être diffusée, et son caractère opérationnel est manifeste.