<div>Thermal management in next generation power electronic devices, radar applications and semiconductor packaging architectures is becoming increasingly challenging due to the need to reject localized high heat fluxes as well as large total powers. Air cooling has been considered as a simple and reliable method for thermal management compared to architectures that incorporate liquid cooling. However, air-cooled heat sinks typically require effective heat spreading to provide the requisite level of area enhancement to dissipate high heat fluxes. Compared to solid metallic heat spreaders, advanced heat sinks that incorporate two-phase heat transfer devices such as vapor chambers can significantly enhance the power dissipation capabilities in such configurations. Vapor chambers are devices that utilize evaporation/boiling processes within a sealed cavity to achieve efficient heat spreading. In high-heat-flux applications, boiling can occur within the internal wick structure of the vapor chamber at the location of the heat input (i.e., the evaporator). The maximum dryout heat flux and thermal resistance of the device is dictated by the resulting two-phase flow and heat transfer in the porous evaporator due to boiling. While various works in the literature have introduced new evaporator wick designs to improve the dryout heat flux during boiling, the enhancement is limited to small, millimeter scale hotspots or at a very high thermal resistance. In additixon, the effective design of such evaporator systems requires mechanistic models that can accurately predict the dryout limit and thermal performance. </div><div> This thesis first explores the usage of a novel ‘two-layer’ evaporator wick for passive high heat flux dissipation over large heater areas at a low thermal resistance. Moreover, a new mechanistic (first principles based) model framework is introduced for dryout limit and thermal performance prediction during boiling in capillary fed evaporators, by considering the resulting simultaneous flow of two phases (liquid and vapor) within the microscale porous media.</div><div> The novel two-layer wick concept uses a thick ‘cap’ layer of porous material to feed liquid to a thin ‘base’ layer through an array of vertical liquid-feeding ‘posts’. Vapor ‘vents’ in the cap layer allow for vapor formed during the boiling process (which is constrained to the base layer) to escape out of the wick. This two-layer structure decouples the functions of liquid resupply and capillary-fed boiling heat transfer, making the design realize high heat flux dissipation greater than 500 W/cm2 over large heat input areas of ~1 cm2. A reduced-order model is first developed to demonstrate the performance of a vapor chamber incorporating such a two-layer evaporator wick design. The model comprises simplified hydraulic and thermal resistance networks for predicting the capillary-limited maximum heat flux and the overall thermal resistance, respectively. The reduced-order model is validated against a higher fidelity numerical model and then used to analyze the performance of the vapor chamber with varying two-layer wick geometric feature sizes. The fabrication of the proposed two-layer wick is then presented. The thermal performance of the fabricated wicks is characterized using a boiling test facility that utilizes high speed visualization to identify the characteristic regimes of boiling operation in the wicks. The performance is also benchmarked to conventional single-layer wicks. </div><div> It is observed that single-layer wicks exhibit an unfavorable boiling regime where the center of the heater area dries out locally, leading to a high value of thermal resistance. The two-layer wicks avoid local dryout due to the distributed feeding provided by the posts and enhance the dryout heat flux significantly compared to single-layer wicks. A two-layer design that consists of a 10 × 10 array of liquid feeding posts provided a 400% improvement in the dryout heat flux. Following a parametric analysis of the effect of particle size, two-layer wicks composed of 180 – 212 µm particles and a 15 × 15 array of liquid feeding posts yielded a maximum heat flux dissipation of 485 W/cm2 over a 1 cm2 heat input area while also maintaining a low thermal resistance of only ~0.052 K/W. The effect of vapor venting and liquid-feeding areas is also experimentally studied. By understanding these effects, a parametrically optimized design is fabricated and shown to demonstrate an extremely high dryout limit of 512 W/cm2. We identify that the unique area-scalability of the two-layer wick design allows it to achieve an unprecedented combination of high total power and low-thermal-resistance heat dissipation over larger areas than was previously possible in the literature.</div><div> The results from the characterization of two-layer wicks revealed that the overall performance of the design was limited by the boiling process in the thin base wick layer. A fundamental model-based understanding of the resulting two-phase flow and heat transfer process in such thin capillary-fed porous media was still lacking. This lack of a mechanistic model precluded the accurate prediction of dryout heat flux and thermal performance of the two-layer wick. Moreover, such an understanding is needed for the optimal design of advanced hybrid evaporator wicks that leverage capillary-fed boiling. Despite the existence of various experimental works, there are currently no mechanistic approaches that model this behavior. To fill this unmet need, this thesis presents a new semi-empirical model for prediction of dryout and thermal resistance of capillary-fed evaporator systems. Thermal conduction across the solid and volumetric evaporation within the pores are solved to obtain the temperature distribution in the porous structure. Capillary-driven lateral liquid flow from the outer periphery of the evaporator to its center, with vapor flow across the thickness, is considered to obtain the local liquid and vapor pressures. Experiments are conducted on sintered copper particle evaporators of different particle sizes and heater areas to collect data for model calibration. To demonstrate the wider applicability of the model for other types of porous evaporators, the model is further calibrated against a variety of dryout limit and thermal resistance data collected from the literature. The model is shown to predict the experimentally observed trends in the dryout limit with mean particle/pore size, heater size, and evaporator thicknesses. This physics–based modeling approach is then implemented into a vapor chamber model to predict the thermal performance limits of air-cooled heat sinks with embedded vapor chambers. The governing energy and momentum equations of a low-cost analytical vapor chamber modeling approach is coupled with the evaporator model to capture the effect of boiling in the evaporator wick. An example case study illustrating the usage of the model is demonstrated and compared to a purely evaporation-based modeling approach, for quantifying the differences in dryout limit prediction, signifying the need to account for boiling in the evaporator wick. </div><div> The understanding gained from this thesis can be utilized for the prediction of dryout and thermal performance during boiling in capillary limited evaporator systems. The work also suggests the usage of a universal relative permeability correlation for the two-phase flow configuration studied herein for capillary-fed boiling, based on a wide calibration to experimental data. The modeling framework can also be readily leveraged to find novel and unexplored designs of advanced evaporator wicks. From an application standpoint, the new vapor chamber model developed here can be used for the improved estimation of performance limits specifically when high heat fluxes are encountered by the device. This will enable better and informed design of air-cooled heat sink architectures with embedded vapor chambers for high performance applications. </div><div><br></div>

Extremely small pore size, low porosity, and ultra-low permeability are among the characteristics of shale rocks. In tight shale reservoirs, the nano-confinement effects that include large gas-oil capillary pressure and critical property shifts could alter the phase behaviors, thereby affecting the oil or gas production. In this research, two in-house simulation models, i.e., a compositionally extended black-oil model and a fully composition model are developed to examine the nano-pore confinement effects on convective and diffusive multicomponent multiphase fluid transport. Meanwhile, the effect of nano-confinement and rock intrinsic properties (porosity and tortuosity factor) on predicting effective diffusion coefficient are investigated. First, a previously developed compositionally extended black-oil simulation approach is modified, and extended, to include the effect of large gas-oil capillary pressure for modeling first contact miscible (FCM), and immiscible gas injection. The simulation methodology is applied to gas flooding in both high and very low permeability reservoirs. For a high permeability conventional reservoir, simulations use a five-spot pattern with different reservoir pressures to mimic both FCM and immiscible displacements. For a tight oil-rich reservoir, primary depletion and huff-n-puff gas injection are simulated including the effect of large gas-oil capillary pressure in flow and in flash calculation on recovery estimations. A dynamic gas-oil relative permeability correlation that accounts for the compositional changes owing to the produced gas injection is introduced and applied to correct for changes in interfacial tension (IFT), and its effect on oil recovery is examined. The results show that the simple modified black-oil approach can model well both immiscible and miscible floods, as long as the minimum miscibility pressure (MMP) is matched. It provides a fast and robust alternative for large-scale reservoir simulation with the purpose of flaring/venting reduction through reinjecting the produced gas into the reservoir for EOR. Molecular diffusion plays an important role in oil and gas migration in tight shale formations. However, there are insufficient reference data in the literature to specify the diffusion coefficients within porous media. Another objective of this research is to estimate the diffusion coefficients of shale gas, shale condensate, and shale oil at reservoir conditions with CO2 injection for EOR/EGR. The large nano-confinement effects including large gas-oil capillary pressure and critical property shifts could alter the phase behaviors. This study estimates the diffusivities of shale fluids in nanometer-scale shale rock from two perspectives: 1) examining the shift of diffusivity caused by nanopore confinement effects from phase change (phase composition and fluid property) perspective, and 2) calculating the effective diffusion coefficient in porous media by incorporating rock intrinsic properties (porosity and tortuosity factor). The tortuosity is obtained by using tortuosity-porosity relations as well as the measured tortuosity of shale from 3D imaging techniques. The results indicated that nano-confinement effects could affect the diffusion coefficient through altering the phase properties, such as phase compositions and densities. Compared to bulk phase diffusivity, the effective diffusion coefficient in porous shale rock is reduced by 102 to 104 times as porosity decreases from 0.1 to 0.03. Finally, a fully compositional model is developed, which enables us to process multi-component multi-phase fluid flow in shale nano-porous media. The validation results for primary depletion, water injection, and gas injection show a good match with the results of a commercial software (CMG, GEM). The nano-confinement effects (capillary pressure effect and critical property shifts) are incorporated in the flash calculation and flow equations, and their effects on Bakken oil production and Marcellus shale gas production are examined. The results show that including oil-gas capillary pressure effect could increase the oil production but decrease the gas production. Inclusion of critical property shift could increase the oil production but decrease the gas production very slightly. The effect of molecular diffusion on Bakken oil and Marcellus shale gas production are also examined. The effect of diffusion coefficient calculated by using Sigmund correlation is negligible on the production from both Bakken oil and Marcellus shale gas huff-n-puff. Noticeable increase in oil and gas production happens only after the diffusion coefficient is multiplied by 10 or 100 times. / Doctor of Philosophy / Shale reservoir is one type of unconventional reservoir and it has extremely small pore size, low porosity, and ultra-low permeability. In tight shale reservoirs, the pore size is in nanometer scale and the oil-gas capillary pressure reaches hundreds of psi. In addition, the critical properties (such as critical pressure and critical temperature) of hydrocarbon components will be altered in those nano-sized pores. In this research, two in-house reservoir simulation models, i.e., a compositionally extended black-oil model and a fully composition model are developed to examine the nano-pore confinement effects on convective and diffusive multicomponent multiphase fluid transport. The large nano-confinement effects (large gas-oil capillary pressure and critical property shifts) on oil or gas production behaviors will be investigated. Meanwhile, the nano-confinement effects and rock intrinsic properties (porosity and tortuosity factor) on predicting effective diffusion coefficient are also studied.

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.

Στην παρούσα διατριβή, ερευνάται αν οι διακυμάνσεις της τριχοειδούς πίεσης είναι δυνατό να παράσχουν ποσοτική πληροφορία που αφορά τη χωρική κατανομή της διαβροχής σε πορώδη μέσα κλασματικής διαβροχής. Για το σκοπό αυτό πραγματοποιήθηκαν πειράματα εκτόπισης της ελαϊκής φάσης από υδατική υπό σταθερή παροχή σε πορώδη μέσα κλασματικής διαβροχής. Τα πορώδη μέσα κλασματικής διαβροχής κατασκευάστηκαν με τυχαία ανάμιξη γυάλινων και πλαστικών (PTFE) σφαιρών, απαράλλαχτης διαμέτρου, έτσι ώστε η γεωμετρία των πόρων να διατηρηθεί αναλλοίωτη. Η χωρική κατανομή των διαφορετικών ειδών σφαιρών διέφερε ακόμη και όταν τα σχετικά τους κλάσματα ήταν ίδια. Κατά τη διάρκεια των πειραμάτων λαμβάνονταν μετρήσεις της μεταβατικής απόκρισης της τριχοειδούς πίεσης και στιγμιότυπα της εκτόπισης ανά τακτά χρονικά διαστήματα. Αρχικά πραγματοποιήθηκαν πειράματα εκτόπισης σε πορώδη μέσα που ήταν ενδιάμεσα υδατοδιαβρεκτά και ισχυρά ελαιοδιαβρεκτά. Τα ληφθέντα σήματα τριχοειδούς πίεσης, αναλύθηκαν με τη μέθοδο των κυματιδίων (wavelet) και συσχετίσθηκαν με τα γεγονότα ροής που ελάμβαναν χώρα. Η διαβροχή των περιοχών των πορωδών μέσων που εισδύονταν από υδατική φάση περιγράφηκε μέσω ενός εμπειρικού δείκτη διαβροχής ο οποίος συσχετίστηκε με τα φάσματα λεπτομέρειας που προέκυψαν από την ανάλυση των σημάτων. Στην περίπτωση ομοιόμορφης και ενδιάμεσης υδατοδιαβροχής ευνοούνται γεγονότα ροής που συνοδεύονται από σημαντικές διακυμάνσεις της τριχοειδούς πίεσης ώστε να καθίσταται δύσκολη η παραπάνω συσχέτιση. Για το λόγο αυτό πραγματοποιήθηκαν πειράματα εκτόπισης με χρήση ενός δεύτερου ζεύγους ρευστών. Τα πορώδη μέσα κλασματικής διαβροχής ως προς το δεύτερο ζεύγος ρευστών ήταν υδατοδιαβρεκτά και ισχυρά ελαιοδιαβρεκτά. Στην περίπτωση του δεύτερου ζεύγους ρευστών, τα σήματα τριχοειδούς πίεσης χαρακτηρίζονταν από αργές και μικρές σε ύψος διακυμάνσεις. Αυτή η διαφορετικότητα επέβαλλε τη βελτίωση της μεθόδου ανάλυσης των σημάτων ώστε να μπορεί να εφαρμοστεί σε όλα τα σήματα που ελήφθησαν στα πειράματα εκτόπισης της παρούσας εργασίας. Υιοθετήθηκαν πλέον δύο εμπειρικοί δείκτες διαβροχής: ο πρώτος δείκτης, Λ, δηλώνει τη διαβροχή της ενεργούς διεπιφάνειας και ο δεύτερος δηλώνει τη διαβροχή της επιφάνειας του πορώδους μέσου που γεμίζει με υδατική φάση όταν κινηθεί η προαναφερθείσα διεπιφάνεια. Οι δύο δείκτες διαβροχής συσχετίστηκαν με την «ενέργεια» των σημάτων, τη «συχνότητα» των φασμάτων λεπτομερειών του βέλτιστου επιπέδου ανάλυσης των σημάτων και τις χρονομέσες τιμές τριχοειδούς πίεσης. / A transparent porous medium of controlled fractional wettability is fabricated by mixing water-wet glass microspheres with strongly oil-wet polytetrafluouroethylene microspheres and packing them between two transparent glass plates. The displacement experiments performed for two pairs of fluids: a) silicon oil and water and b) paraffin oil and water. The growth pattern is video recorded and the transient response of the pressure drop across the pore network is measured for various fractions of oil-wet particles. The measured global capillary pressure fluctuates as the result of the variation of the equilibrium curvature of menisci between local maxima and local minima. By using wavelet transform on the capillary pressure signal, the best level detail spectrum arises and its peaks indict the parts of the original signal with the most important fluctuations, each of those is analyzed again. The best level detail coefficients give the “energy” of the analyzed part of the capillary pressure signal while the number of the most important fluctuations of the best level detail spectrum to the time interval they took place gives the “frequency” of the best level detail spectrum. The time averaged capillary pressure, the “energy” and the “frequency” are associated with the frontal wettability of the interface separating the two fluids before the invasion step and one number which denotes the regional wettability of the invaded area.

Two-phase flow in porous media is an important subsurface process that has significant impacts on the global economy and environments. To study two-phase system in porous media, capillary pressure (Pc ), relative permeability (Kr), bulk electrical conductivity (σb) and bulk relative permittivity (εb) are often employed as characterization parameters. Interestingly, all of these parameters are functions of water saturation (S). However, the non-uniqueness in the Pc -S, Kr-S,σb-S and εb-S relationships pose considerable challenges in employing them for effective monitoring and control of the two-phase flow processes. In this work, laboratory scale experiments and numerical simulations were conducted to investigate the factors and conditions contributing to the non-uniqueness in the above relationships for silicone oil-water and supercritical CO2-water flow in porous media, with a special emphasis on geological carbon sequestration. Specifically, the dynamic capillary pressure effect, which indicates the dependence of the Pc - S relationship on the rate of change of saturation (αS/αt) during two-phase flow in porous media was investigated. Using a silicone oil-water system, the dynamic capillary pressure effect was quantified in term of the parameter named the dynamic coefficient, τ , and it was found to be dependent on the domain scale and the viscosity ratio of the two fluids. It was found that τ increases with the domain scale and the viscosity ratio. It is inversely affected by αS αt , which is related to the degree of resistance to the fluid motion, namely, viscosity. In almost all cases, τ was found to decrease monotonically with an increase in water saturation, S. An order increase in magnitude of τ was observed as the domain scale increases from 4cm scale to 8cm in height. A similar order of increase in τ was observed in the 12cm high domain scale. There is an order increase in the value of τ for the silicone oilwater system as the viscosity ratio increases from 200 to 500. For the supercritical CO2 (scCO2) and water system in porous media, the experiments and numerical simulations showed that τ increases with rising system temperature and decreasing porous media permeability. Dimensionless analysis of the silicone oil-water experimental results showed that by constructing non-dimensional groups of quantities expressing a relationship among different variables on which τ depends, it is possible to summarise the experimental results and determine their functional relationship. A generalised scaling relationship for τ was derived from the dimensionless analysis which was then validated against independent literature data. The exercise showed that the τ-S relationship obtained from the literature and the ii scaling relationship match reasonably well. This work also demonstrated the applicability of an artificial neural network (ANN) as an alternative computational platform for the prediction of the domain scale dependence of τ . The dependence of the Kr-S relationship on αS/αt was also investigated. The results showed that the Kr-S curve under dynamic flow condition is different from that under the quasi-static condition. Kr for water (Krw) increases with increasing water saturation and decreases with the increase in viscosity ratio while Kr for silicone oil (Krnw) increases with decreasing water saturation as well as with the increase in viscosity ratio. Also, Krw decreases while Krnw increases with the increasing boundary pressure. However, the εb-S and σb-S relationships were found to be independent of αS/αt for the scCO2-water system in carbonate and silicate porous media. Nevertheless, the εb and σb values decrease as the water saturation decreases in the two porous media samples. While εb decreases with increase in temperature in silica sand, the trend in the limestone showed a slight increase with temperature, especially at high water saturation. Also, the εb-S relationship is shown to be affected by pressure in silica sand increasing with the pressure of the domain. On the contrary, the σb-S relationship increases as the temperature increases with more significance at higher water saturation in the silica sand sample. This work further demonstrated the application of a membrane in the monitoring of the CO2 in geological sites used for carbon sequestration. Commercial silicone rubber coupled with a pressure transducer showed potential in the detection of CO2 leakage from geological sites. The response of the device in terms of the mass of permeated gas, permeability and gas flux were investigated for both CO2 and N2. In addition, the monitoring of potable water contamination in a shallow aquifer by the migrating or leaking of CO2 is demonstrated with the combination of the pH analysis, geoelectrical measurement techniques and the membrane-sensor system. Overall, the work in this PhD research demonstrated robust applications of two-phase systems'characterization parameters under different scenarios in the porous media. Implications of the findings in this work to the monitoring and control of two-phase systems in porous media are expatiated.

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

Cette thèse porte sur la modélisation et la discrétisation d’écoulements Darcy dans les milieux poreux fracturés. Nous suivons l’approche des modèles, dits à dimensions hybrides, qui représentent les réseaux de fractures comme des surfaces de codimension 1 immergées dans la matrice. Les modèles considérés prennent en compte les interactions entre matrice et fractures et permettent de traiter des fractures agissant comme conduites ou comme barrières, ce que nécessite de prendre en compte les sauts de pression aux interfaces matrice-fracture. Dans le cas des écoulements diphasiques, nous proposons des modèles, qui prennent en compte les sauts de saturations aux interfaces matrice-fracture, dû à la capillarité. L’analyse numérique est menée dans le cadre général de la méthode de discrétisations gradients, qui est étendue aux modèles considérés. Deux familles de schémas numériques, le schéma Vertex Approximate Gradient et le schéma Volumes Finis Hybrides sont adaptées aux modèles à dimensions hybrides. On prouve via des résultats de densité que ce sont des schémas gradients, pour lesquels la convergence est établie. En diphasique, l’existence d’une solution est obtenue en passant. Plusieurs cas tests sont présentés. En monophasique, on observe la convergence sur des différents types de mailles pour une famille de solutions dans un milieux fracturé hétérogène et anisotrope. En diphasique, nous présentons une série de cas tests afin de comparer les modèles à dimensions hybrides au modèle de référence, dans lequel les fractures ont la même dimension que la matrice. A part quantifier le gain en performance de calcul, ces tests montrent la qualité des différents modèles réduits. / This thesis investigates the modelling of Darcy flow through fractured porous media and its discretization on general polyhedral meshes. We follow the approach of hybrid dimensional models, invoking a complex network of planar fractures. The models account for matrix-fracture interactions and fractures acting either as drains or as barriers, i.e. we have to deal with pressure discontinuities at matrix-fracture interfaces. In the case of two phase flow, we present two models, which permit to treat gravity dominated flow as well as discontinuous capillary pressure at the material interfaces. The numerical analysis is performed in the general framework of the Gradient Discretisation Method, which is extended to the models under consideration. Two families of schemes namely the Vertex Approximate Gradient scheme (VAG) and the Hybrid Finite Volume scheme (HFV) are detailed and shown to fit in the gradient scheme framework, which yields, in particular, convergence. For single phase flow, we obtain convergence of order 1 via density results. For two phase flow, the existence of a solution is obtained as a byproduct of the convergence analysis. Several test cases are presented. For single phase flow, we study the convergence on different types of meshes for a family of solutions. For two phase flow, we compare the hybrid-dimensional models to the reference equidimensional model, in which fractures have the same dimension as the matrix. This does not only provide quantitative evidence about computational gain, but also leads to deep insight about the quality of the proposed reduced models.