L’industrie aéronautique a exprimé un besoin en matériaux pour des zones fortement chargées thermiquement et mécaniquement. L’objectif est d’optimiser grâce à ces matériaux, notamment en termes de poids et de rendement, le coeur des turboréacteurs. Un procédé par voie liquide appelé «Reactive Melt infiltration» est industriellement envisagé pour fabriquer ces matériaux. La densification est tributaire de la compétition entre la montée capillaire et la réaction chimique entre le silicium liquide et la poudre préalablement introduite. Cette concurrence peut conduire à des phénomènes de «choking off» qui doivent être évités. Dans ce travail, l’approche numérique est à la fois macroscopique et microscopique. Des validations numériques, des applications dans le cadre de géométries modèles et des analyses physiques sont présentées pour les deux échelles. Les simulations d’imbibition réactive à l’échelle de la pièce sont effectués avec un outil développé au sein de SAFRAN alors que les simulations de montées capillaires dans des géométries modèles à l’échelle du pore sont réalisées avec le code de calcul Thétis (développé à l’I2M, Bordeaux). Une méthodologie expérimentale pour le suivi et l’analyse de l’imprégnation capillaire réactive d’un milieu granulaire constitué de matériaux modèles est également présentée. / The development of ceramic matrix composites (CMCs) to replace certain metal components ininternal hot parts of aircraft engines is an active research field for the aeronautical industry. Theseadvanced components may be manufactured with a fluid processing called ReactiveMelt Infiltration(RMI). The densification step is the capillary rise of a molten metal such as silicon within a wovenpreform. The molten metal can react with the previously introduced ceramic powder and lead to arefractory matrix. The competition between capillary rise and the reaction between molten siliconand introduced carbon powder has to be managed. The ambition is to prevent choking off effectsand closed pores created by preferential paths. Simulations of two-phase flows at macroscopic scaleare undertaken with a homemade code developed in SAFRAN whereas simulations of two-phaseflows atmicroscopic scale were undertaken with CFD code Thétis (developed at I2M, Bordeaux). Anexperimental methodology for analysis of reactive or not reactive capillary impregnation in granularmedia is described.

La combustion du noir de carbone (NC) dans le creuset d'une thermobalance est contrôlée à la fois par la réaction et par le transport de l'oxygène jusqu'à la surface du lit et à l'intérieur du lit poreux de NC.Les expériences menées en modifiant la masse de NC ont montré que la concentration en oxygène peut tomber à zéro avant d'atteindre le fond du lit. Ainsi, à un instant donné, seule une partie du lit est en combustion. Cette masse, appelée masse critique (mc) dépend de la température. Elle passe de 35 mg à 570°C à 17,5 mg à 650°C.Un gradient d'oxygène s'établi donc dans le lit. La modélisation du transport interne de l'oxygène a montré que la diffusion de Fick constitue une bonne approximation pour représenter ce transport.Des conseils pour extraire correctement une constante cinétique à partir d'expériences thermogravimétrique sont donnés. La procédure est adaptée en fonction de la précision souhaitée.Ainsi, l'utilisation de faibles masses afin de réduire au mieux l'effet de la masse et l'exothermicité de la réaction est fortement conseillée. L'influence de la diffusion externe du gaz peut être réduite en utilisant des creusets de très faibles hauteurs, ou en mettant l'échantillon le plus proche de la bouche du creuset en remplissant le fond du creuset avec un matériau inerte. / Combustion of carbon black (CB) in the crucible of a thermobalance is controlled by both carbon reactivity and oxygen transport from the oxidizing flux to the surface of the bed and within the porous bed.The experiments conducted by changing the mass of CB showed that the oxygen concentration can fall to zero before the bottom of the bed. Thus, at a given time, only a part of the bed is burning. This mass, called critical mass (mc), depends to temperature. It went from 35 mg at 570°C to 17.5 mg at 650°C.An oxygen gradient is thus established in the bed. The Modelling of the internal transport of oxygen showed that the Fick diffusion is a good approximation to represent the transport.Advices to correctly extract a kinetic constant using thermogravimetric experiments are given. The procedure is adjusted depending to the precision desired.Thus, the use of low masses to best reduce the mass and exothermic reaction effects is strongly recommended. The influence of stagnant gas can be reduced by using crucibles with very low height, or by placing the sample closest to the mouth of the crucible by filling the bottom of the crucible with an inert material.

Ce travail porte sur l’étude de la propagation de décharges couronnes impulsionnelles à pression atmosphérique dans les milieux poreux et/ou alvéolaires. Face à la complexité des phénomènes mis en jeu, liés aux interactions entre la décharge et les surfaces du matériau qui la confine, nous proposons l’étude de décharges confinées par des structures élémentaires. L’étude du confinement radial des décharges, assuré par un large panel de capillaires, a été réalisée. Des diagnostics électriques et optiques de pointe permettent d’étudier la propagation de la décharge au sein des différents capillaires. La corrélation entre ces diagnostics a même permis des mesures de vitesse de propagation au sein de capillaires opaques. Les résultats montrent que la propagation de la décharge dépend grandement de la géométrie des capillaires et des paramètres électriques de génération de la décharge. Dans le cas de sections carrées ou rectangulaires, les arêtes induisent un renforcement local du champ qui attire la décharge. Dans le cas de capillaires cylindriques, le diamètre interne est le paramètre crucial qui détermine aussi bien la structure de la décharge que sa vitesse de propagation. Quelle que soit la nature du capillaire, la propagation présente alors une vitesse optimale à tout autre paramètre constant pour une valeur donnée du diamètre interne. Dans le cas du verre, la vitesse est maximale pour un diamètre interne de 200 µm. L’épaisseur et la permittivité diélectrique du capillaire possèdent également une influence sur la propagation de la décharge radialement confinée. Ainsi, diminuer l’épaisseur ou la permittivité diélectrique engendre une accélération de la décharge. Si l’épaisseur est très faible, la décharge peut même se déconfiner pour se propager à l’extérieur du capillaire. Une étude spectroscopique complémentaire montre que la réduction du diamètre de confinement implique une augmentation de la température du plasma, ce qui pourrait contribuer à l’obtention de ce profil de vitesse en fonction du diamètre de confinement. L’étude du confinement axial des décharges a ensuite été réalisée en insérant des membranes de différentes natures et caractéristiques, perpendiculairement à l’axe pointe plan. Les résultats montrent que la décharge présente une propagation en trois étapes : pointe/membrane, radialement au voisinage de la membrane, puis membrane/plan. Dans cette étude, nous avons mis en évidence l’importance du critère poreux ou non de la membrane. Dans le cas poreux, la propagation de la décharge dans l’ensemble du gap est continue, même pour des pores de l’ordre de la dizaine de µm. Dans le cas non poreux, la propagation est discontinue, et il est nécessaire pour assurer la propagation dans l’ensemble du gap qu’un ré-allumage ait lieu de l’autre côté de la membrane. Après l’instant de l’impact sur la membrane, la décharge marque un arrêt qui correspond à la réorganisation des charges et à la restructuration du champ électrique dans le gap. Elle se propage ensuite radialement au voisinage de la membrane en plusieurs fronts d’ionisation. Si les conditions de claquage sont réunies dans le volume membrane/plan, alors un ré-allumage apparaît à partir de la membrane pour atteindre le plan. L’étude de ces ré-allumages semble montrer l’importance de la position de la membrane au sein de l’espace inter-électrodes et de la dynamique des charges aux surfaces de la membrane. Plus on diminue la distance membrane/plan, plus il est facile d’en observer. Nous montrons également que la diminution de la permittivité diélectrique de la membrane ou l’augmentation de son épaisseur, semble augmenter la probabilité de ces ré-allumages. Dans le cas poreux, nous avons également mis en évidence l’influence de la taille des pores de la membrane sur l’ensemble des étapes de propagation. Lorsque la porosité est inférieure à 100 µm la propagation de la décharge est ralentie du fait de la difficulté de la décharge à traverser directement le matériau. / This study is an attempt to understand the mechanisms involved in the propagation of pulsed corona discharges at atmospheric pressure inside porous and/or alveolar media. Due to the complexity of these phenomena which hardly depends on plasma/surface interactions, the study was focused first on the propagation inside basic structures of confinement, before doing it in more complex media, such as monolithic cordierites or foams. Therefore, capillaries have been used to radially confine the discharge propagation. Thanks to highly resolved optical and current diagnostics, we succeed in describing precisely the propagation. A correlation of these diagnostics allowed us to measure propagation velocity inside opaque media. Results show that geometry is the key parameter which both governs the discharge structure and the propagation velocity. Electrical parameters of the discharge ignition have a great role in the propagation also. In case of square and rectangular capillaries, the local electric field enhancement due to edges attracts the discharge. In case of round capillaries, the internal diameter becomes the key parameter which governs the pattern of the discharge and its velocity. Whatever the kind of capillary used, the velocity of propagation shows an optimal value for a given internal diameter, at all others parameters constants. In case of glass capillaries, this maximum value is obtained for a 200 µm internal diameter. Parameters such as wall thickness or dielectric permittivity have also an influence on the propagation velocity. The decrease of the thickness or the dielectric permittivity implies an increase of the propagation velocity. If the thickness is small enough, we observe that the discharge is able to propagate outside the tube. A complementary spectroscopic study of that kind of discharge shows that the plasma temperature depends on the internal diameter, and warming observed for smallest internal diameters could contribute to the internal diameter/velocity of propagation relation. The study of the axial confinement of the propagation of the discharge has been made thanks to different kinds of dielectric membranes, perpendicularly inserted between the electrodes. Results show that the propagation of the discharge is a three step process: tip to membrane propagation, radial propagation near the membrane surfaces, and membrane to plane propagation under specific conditions. The porous feature of the membrane has a key role in the propagation. In case of porous membranes, the whole propagation in the gap is continuous, even for the smallest porosities under investigations (~10µm). In case of non porous films, the propagation by three step processes is no more continuous, and the propagation in the entire gap needs a re-ignition of the discharge on the other side of the membrane. After its impact on the membrane, the discharge stops on the dielectric while the charges and the electric field are self reorganized. Then a radial propagation of several ionization waves starts near the membrane surface. If the disruptive voltage is reached behind the obstacle, a re-ignition could occur in the membrane/plane gap and reach the cathode. Investigation on these reignition conditions shows that the key parameters seem to be the position of the membrane inside the gap and the dynamics of the surface charges on the membrane. Closer is the membrane to the plane, higher is the probability of seeing re-ignitions. Results also show that the decrease of the dielectric permittivity or the increase of the membrane thickness leads to more re-ignition events. In case of porous membrane, the pore size is the main parameter which will influence the propagation. When the pore size is below 100 µm, the discharge propagation is slower due to the discharge difficulties to penetrate totally inside the material.

La pyrolyse est la décomposition thermochimique de la biomasse en gaz de synthèse valorisables en biocarburants. Cette technologie, propre et renouvelable, nécessite aujourd’hui des efforts de recherche et de développement afin de prouver sa compétitivité par rapport aux autres sources d’énergie. L’objectif de cette thèse est de développer un modèle macroscopique de la pyrolyse en utilisant la méthode de prise de moyenne volumique. Le modèle sera ensuite utilisé pour faire des études numériques afin de caractériser le procédé et améliorer les performances des réacteurs. Une approche en deux temps a été établie afin d’atteindre notre objectif. D’abord, des modèles macroscopiques d’écoulements à masse volumique variable en milieu poreux ont été développés. Ce type d’écoulements est similaire à celui mis en jeu en pyrolyse pour deux deux raisons: la masse volumique varie sous l’effet de gradients forts de température et le réacteur de pyrolyse peut être considéré comme un milieu poreux à double porosité (porosité à l’échelle du lit et porosité à l’échelle de la particule). Les résultats théoriques ont montré que les équations de conservation macroscopiques (continuité, quantité de mouvement et énergie) et les propriétés effectives (masse volumique, perméabilité et diffusivité thermique) font apparaitre de nouveaux termes résultants de la variation de densité. La forme explicite de ces termes a été établie et validée par simulations numériques. Les résultats obtenus ont été utilisés dans un deuxième temps afin de développer un modèle macroscopique de la pyrolyse. / Pyrolysis is a thermo-chemical conversion of biomass into bio-fuels. This technology has not been fully developed and its competitiveness against other sources of energy is yet to be proven. The aim of this work is to derive a macroscopic model of pyrolysis by means of volume averaging method. The obtained macroscopic model can then be used to conduct fast and low-cost numerical simulations to characterize the process and improve the reactor efficiency. To achieve our objective, a two-steps methodology has been established. First, the fundamental problem of variable density flow in porous media has been investigated. The physical phenomena in this kind of problem are very similar to those involved in pyrolysis for two reasons: the fluid density varies due to high temperature gradients and the pyrolysis reactor can be considered as a double porosity medium (porosity at the reactor scale and porosity at the biomass particle scale). The obtained macroscopic conservation equations (continuity, momentum and energy) and the effective properties (density, permeability and thermal diffusivity) contain additional terms resulting from the fluid density variation. The explicit form of these terms has been established and their components have been computed. The resulting models of the first step have then been used to develop a macroscopic model of the pyrolysis in the second part of our study.

Le Laboratoire de Thermique, Energétique et Procédés de Pau travaille en parallèle depuis plusieurs années sur les deux thématiques suivantes, la modélisation des réacteurs multifonctionnels et la modélisation des transports en milieux poreux. Dans ce travail, un modèle général de transports en milieu poreux multiphasique multiconstituant réactif est développé puis appliqué au cas particulier de l’absorption réactive afin de décrire les phénomènes se produisant au sein du garnissage. La modélisation débute par une description classique de chacune des phases continues grâce aux équations de conservation de la masse, de la quantité de mouvement et de l’énergie. Les phases fluides sont multiconstituantes et le siège de réactions homogènes. Bien qu’à cette échelle les mécanismes soient parfaitement décrits, le passage à la simulation impose une étape d’homogénéisation par prise de moyenne. Ce changement d’échelle conduit à un système d’équations à l’échelle locale. Le modèle est appliqué à l’absorption du dioxyde de carbone dans une solution aqueuse d’hydroxyde de sodium au sein d’une colonne garnie. Les simulations bidimensionnelles instationnaires fournissent des informations locales sur le procédé. Les résultats obtenus montrent une bonne adéquation avec la connaissance classique des mécanismes d’absorption réactive. / The “Laboratoire de Thermique, Energétique et Procédés de Pau” works on the two following subjects, modeling of multifunctional reactor and modeling of transport in porous media. In this work, a general model of transport in multiphase multi component reactive porous media is developed and then applied to the specific case of reactive absorption in order to describe the phenomena which occur within the packing. Modeling begins with a classic description of each continuous phase through the equations of conservation of mass, momentum and energy. The fluid phases are multi component and reactive. Although mechanisms are perfectly described at this scale, the simulation imposes a homogenization by volume averaging. This change of scale leads to a set of equations at the local scale. The model is applied to the absorption of carbon dioxide in an aqueous solution of sodium hydroxide within a packed column. Unsteady two-dimensional simulations provide local information about the process. The results obtained show a good agreement with the classical knowledge of reactive absorption mechanisms.

Ce travail s'intéresse a` l'identification des mécanismes physiques en jeu lors de la déshydratation assistée thermiquement. Une cellule de filtration/compression a été instrumentée avec des capteurs de pression liquide. Ces mesures permettent de suivre la formation du gâteau de filtration, sa consolidation puis les évolutions thermomécaniques induites par l’apport de chaleur. Sous l'effet du chauffage, elles montrent l'apparition d’un gradient de pression liquide, qui induit un écoulement additionnel. Un modèle mathématique est ensuite proposé. Des caractérisations en cellule de compression/perméabilité sont utilisées pour de´terminer les lois d'évolution des propriétés du milieu poreux. Le modèle permet de retrouver les évolutions des grandeurs macroscopiques et locales mais pas les cinétiques. Ceci a été attribué au modèle de déformation du milieu poreux choisi. / In the scope of thermally assisted mechanical dewatering process, this work focuses on the identification of physical mechanisms involved in the thermally assisted mechanical dewatering. Operating conditions ensure that the water is expelled in liquid phase. An experimental study on talc and cellulose saturated suspensions highlights the gains in terms of final dry solid contend involved by the thermal intensification. The filtrationcompression test cell is instrumented with pore liquid sensors along the cake thickness that allows to measures the pressure of the interstitial liquid phase. Obtained data highlight an increase of the liquid pressure in the heated side of the cake during thermally assisted mechanical dewatering. Indeed, as the temperature increases the water density decrease. As the cake consolidation restricts the flow a liquid pressure gradient reappears inside the cake. As the temperatures on the filter side of the cake increase, the thermally induced liquid pressure gradient vanishes with an additional filtrate outflow. In order to measure the porous media properties a compression-permeability cell has been build. This apparatus allows permeability measurement to be conducted for different temperature and loading on the porous media. A physical model including the thermal pressurization has been implemented in COMSOL Multiphysics in order to simulate the process. The model confirms the thermal pressurization occurring during thermally assisted mechanical dewatering, but because of deviation in the calculated mechanical behavior and temperature of the porous media, the model can’t depict the experimental additional filtrate outflow.

La lubrification eX-Poro-HydroDynamique (XPHD) est un mécanisme de lubrification d’inspiration biomimétique. Il s’agit principalement d’un écoulement dans un milieu poreux, dont la phase solide représentée par des fibres, induit des forces élastiques de compression considérées comme négligeables par rapport aux forces hydrodynamiques générées à l’intérieur du milieu poreux. L’idée essentielle de la lubrification XPHD consiste en remplacement du matériau antifriction et du film de fluide mince, traditionnellement utilisé dans les solutions classiques, par une couche poreuse imbibée d’un fluide qui fournit plus grande capacité de charge. Ce type de lubrification représente une solution technologique totalement nouvelle (en rupture avec la solution classique) qui peut permettre de remplacer les lubrifiants pétroliers, de créer des systèmes tribologiques autolubrifiants et donc plus écologiques et moins coûteux. Ce travail de recherche est donc focalisé sur l’évolution des performances de la lubrification XPHD dans le cadre d’un mouvement tangentiel, adapté à l’étude de butées pour faible et moyenne vitesses de rotation. Dans ce contexte scientifique, une étude approfondie d’un matériau poreux présélectionné (mousse en polyuréthane) fut réalisée avec le but de déterminer les caractéristiques physiques et les paramètres cruciaux pour la lubrification XPHD: la porosité et la perméabilité du matériau poreux. Les modèles théorique et numérique de lubrification XPHD proposés se basent sur l’équation de Darcy-Brinkman et les hypothèses de la lubrification classique, ainsi que l’écoulement à l’intérieur du milieu poreux sont prédit avec une nouvelle forme de l’équation de Reynolds. Le banc d’essais spécialement développé permet d’étudier expérimentalement le mécanisme de lubrification XPHD pour des différents types de dislocateur en combinaison avec des liquides newtoniens et nonnewtoniens. Une description détaillée du banc d’essais et de tous les dispositifs expérimentaux utilisés ainsi que la comparaison des résultats de modélisation et des résultats expérimentaux sont présentés. / The eX-Poro-HydroDynamic (XPHD) lubrication is a new biomimetic inspired lubrication mechanism. It consists of self-sustained fluid films generated within highly compressible porous layers imbibed with liquids, whose solid phase represented by fibers, induces compressive elastic forces considered negligible compared to the hydrodynamic forces generated inside the porous medium. The essential idea of XPHD lubrication is to replace the antifriction material and the thin fluid film, traditionally used for classical sliding motion, with a porous layer imbibed with a fluid that provides a greater load capacity. This type of lubrication represents a completely new technological solution (in rupture with the classical one) that can replace petroleum lubricants, create self-lubricating and therefore more ecological and less expensive tribological systems. This research work is focused on the evolution of XPHD lubrication performances in the context of a tangential movement, adapted to the study of thrust bearing for low and medium rotation speeds. In this scientific context, a thorough study of a preselected porous material (polyurethane foam) was carried out with the aim of determining the physical characteristics and the crucial parameters for XPHD lubrication: the porosity and the permeability of the porous material. The theoretical and numerical models of the XPHD lubrication proposed are based on the Darcy-Brinkman equation and the classical lubrication hypothesis, as well as the flow within the porous media is predicted with a new form of the Reynolds equation. The specially developed test rig is used to investigate experimentally the mechanism of the XPHD lubrication for different types of thrust bearing in combination with Newtonian and non-Newtonian liquids. A detailed description of the test rig and all used experimental devices, as well the comparison between the experimental and numerical results are presented.

La pollution d’aquifères par des hydrocarbures est très persistante et difficile à traiter, devenant un enjeu majeur pour l’environnement en raison de l’effet négatif sur la santé humaine. La génération de mousse in situ combinée au lavage de sols est une technique de dépollution innovante, permettant un meilleur contrôle de la mobilité des fluides dans des formations hétérogènes. Le but de cette méthode est de bloquer l’écoulement dans les couches à forte perméabilité pour améliorer l’efficacité de balayage des strates peu perméables.Le mémoire se divise en deux parties principales décrivant une étude expérimentale et une étude numérique.Une sélection préalable de tensioactifs respectueux de l’environnement à été réalisée sur la base de leur moussabilité. Le sucrose laurate (émulsifiant alimentaire) a été retenu et comparé avec un tensioactif dit classique. Des tests sur colonnes 1-D ont permis d’étudier l’influence de la perméabilité, du type de tensioactif et de la qualité de la mousse sur le facteur de résistance à la mobilité du gaz. Les résultats, interprétés à l’aide de la pression capillaire, révèlent deux régimes d’écoulement, correspondants aux deux textures de mousse : « weak » et « strong ». Des expériences sur pilotes 2-D hétérogènes ont permis de visualiser les différentes phases du processus de remédiation et de comparer l’efficacité de balayage avec et sans une région bloquée par de la mousse.La partie modélisation comprend une simulation numérique des tests sur colonne 1-D en utilisant le code UTCHEM et une étude de sensibilité des paramètres les plus importants du modèle, ainsi qu’une modélisation des tests de traçage sur pilote 2-D avec le code MODFLOW. / Aquifer pollution with hydrocarbons is very persistent and difficult to treat, becoming a major issue for the environment because of the negative effect on human health. In situ foam generation combined with soil washing is an innovative remediation technology that allows a better fluid mobility control in heterogeneous formations. The purpose of this method is to block the flow through the high permeability layers in order to improve the sweep efficiency of low permeability strata.The thesis is divided into two main parts describing an experimental study and a numerical study.A preliminary selection of environmentally-friendly surfactants was carried out on the basis of their foamability. Sucrose laurate (emulsifier from food industry) was selected and compared with a conventional surfactant. 1-D column tests were performed to study the influence of permeability, surfactant type and foam quality on the resistance factor. The results, interpreted using the capillary pressure, show two flow regimes, corresponding to the two foam textures: "weak" and "strong". Experiments on 2-D heterogeneous pilot allowed visualizing the different steps of the remediation process and comparing the sweep efficiency with and without a region blocked by foam.The numerical part includes modeling 1-D column tests using UTCHEM, a sensitivity study of the most important parameters of the model, and a simulation of tracer tests on 2-D pilot using MODFLOW.

La production de méthane à partir des champs hydratifères des fonds océaniques est un procédé promis à se développer et à atteindre l'échelle industrielle au cours des dix prochaines années. Cependant, les premiers essais d'extraction de méthane se sont révélés infructueux et difficiles à mener. Le phénomène de dissociation des hydrates de méthane au sein d'une matrice poreuse reste mal connu et maîtrisé: la fusion des hydrates de méthane engendre des écoulements de fluide et des transferts thermiques au sein d'un milieu poreux. Les propriétés de ce milieu poreux évoluent donc avec la disparition de la phase hydrate et la présence éventuelle d'une phase glace au cours de la dissociation. Les transferts de masse et de chaleur doivent donc être couplés, afin de tenir compte premièrement du changement de la perméabilité et deuxièmement de l'évolution de la conductivité thermique du système.<br />Dans cette étude, les transferts de masse et de chaleur ont été étudiés numériquement et expérimentalement. Un modèle numérique 2D est proposé dans lequel les transferts de chaleur et de masse gouvernent la dissociation des hydrates de méthane. Les résultats numériques montrent la présence de gradients de pression et de température au sein du milieu poreux et l'évolution de la frontière de dissociation selon le type de sédiment utilisé. Ce modèle est utilisé afin de dimensionner un dispositif expérimental de dissociation de carottes sédimentaires partiellement saturées en hydrates de méthane qui permet un suivi précis de la cinétique de dissociation.<br />Le montage expérimental est composé de cinq zones cylindriques de même diamètre (1/2 inch) mais de tailles différentes (pour une longueur totale de 2,6 m). Chaque zone est contrôlée en pression et en température. Chaque expérience consiste en une cristallisation d'hydrates de méthane au sein d'un milieu poreux, puis en une dissociation de ces mêmes hydrates par une méthode de contrôle de la pression à une extrémité du dispositif. La cinétique de dissociation est étudiée par le suivi de la pression dans un ballast (situé en aval du dispositif). Les résultats obtenus, via une étude paramétrique, permettent de cibler les paramètres clés de la dissociation de ces hydrates sédimentaires et d'observer leur impact sur la cinétique de dissociation. Deux régimes bien distincts de dissociation sont mis en évidence et caractérisés selon les propriétés du sédiment partiellement saturé en hydrates de méthane. Le rôle de la glace au cours de la dissociation est également étudié pour ces deux types de dissociation. <br />Enfin, la correspondance des résultats numériques et expérimentaux est mise en évidence par comparaison de courbes de cinétique de dissociation et de courbes d'évolution de la pression au sein du milieu poreux.

La connaissance des réservoirs géologiques et de leur évolution est cruciale pour toutes les applications les mettant en jeu (prospection/extraction de pétrole ou de gaz, stockage de gaz ou de déchets, réhabilitation de sites industriels, dépollution des sols et des nappes phréatiques). La complexité des phénomènes couplés de réaction chimique et transport, la non-linéarité des équations en découlant et la diversité des échelles d'observation (en espace et en temps) amènent à avoir recours à la modélisation. Lors des injections de gaz acides dans des structures géologiques potentiellement capables de les stocker, les minéraux présents primaires) les plus réactifs se dissolvent, du fait de l'acidification de l'eau en place. En conséquence, le fluide se charge en éléments issus de ces réactions jusqu'à atteindre un état de sursaturation par rapport à d'autres minéraux primaires, mais aussi par rapport à des minéraux jusque là absents de l'assemblage (secondaires) qui peut les amener à précipiter. Dans ces configurations, la justesse des prédictions du comportement du système, qui permettent de s'assurer de la pérennité du stockage du gaz et de se prémunir contre les risques mécaniques éventuels, va dépendre de l'aptitude du programme à bien choisir les solides qui précipitent et à connaître précisément à la fois les cinétiques de croissance des minéraux primaires et les cinétiques de formation des minéraux secondaires.<br />Dans le code ARCHIMEDE (développé à l'E.N.S.M-S.E, Ecole Nationale Supérieure des Mines de Saint-Etienne, en collaboration avec l'I.F.P, Institut Français du Pétrole), seule la partie géochimie est concernée. Outre quelques difficultés qui empêchaient de traiter les variations de volume occasionnées par les réactions et que nous nous sommes efforcés de résoudre en reprenant l'ensemble du programme, l'objectif principal de ce travail de thèse a été la conception et l'implémentation d'un modèle de germination/mûrissement/croissance pour les minéraux secondaires.<br />Notre démarche a consisté à analyser précisément l'apparition de nouveaux minéraux dans des assemblages naturels d'où ils sont absents initialement. Quelques simulations numériques ont permis de montrer les limites inhérentes à la représentation, inadaptée, de la précipitation de ce genre de minéraux par croissance cristalline. La mise en avant de ces faiblesses a défini alors le cadre pour un nouveau modèle de précipitation propre aux minéraux secondaires. Une étude détaillée de la phase de germination, à savoir l'apparition des premiers cristaux d'un minéral, a conduit à la construction d'un modèle pour sa cinétique. Par suite, la mise en compétition du processus de germination avec celui de croissance, destiné à prendre plus tard le relais de la production de volume du nouveau solide, a nécessité la prise en charge de cristaux de tailles variées et le recours au mûrissement d'Ostwald, via un calcul parallèle. Il en a résulté un algorithme élaboré pour gérer, pour chaque minéral secondaire, la phase initiale de germination/mûrissement et ensuite le passage à la phase de croissance cristalline. La sensibilité de ce nouveau modèle complet de germination/mûrissement/croissance vis-à-vis des différents paramètres qui le composent a pu être analysée moyennant de nouvelles simulations numériques. Son comportement, en terme d'acuité des prédictions, a également pu être mis à l'épreuve et finalement jugé satisfaisant.