Le procédé Resin Transfer Molding est une technologie de fabrication de pièces en matériaux composites. Une compréhension et une modélisation des phénomènes physiques qui interviennent dans certaines étapes de la fabrication apportent une contribution importante à la maîtrise de ce procédé. Les travaux abordés dans cette étude se situent dans le cadre d'un programme plus général visant le contrôle et le pilotage du RTM. Tout d'abord, une caractérisation de la perméabilité des renforts dans la direction perpendiculaire aux plis va compléter les résultats obtenus dans des études préalables sur les perméabilités planes. Ensuite, nous étudierons de près les transferts thermiques durant l'injection et la cuisson. Pour cela, plusieurs étapes sont nécessaires: caractérisation des propriétés thermiques (avec notamment le développement d'un montage destiné à mesurer les conductivités thermiques dans le plan des fibres), étude théorique et étude expérimentale à l'aide d'un moule instrumenté. Des simulations numériques permettront alors de valider les modèles développés. Ces travaux devront ensuite servir de base à l'amélioration du logiciel de simulation LCMFLOT. Enfin, une première étape vers l'optimisation du remplissage a été réalisée afin de pouvoir piloter les machines d'injection en imposant une vitesse d'imprégnation constante pendant le remplissage. Cette optimisation est basée sur l'exploitation des résultats de simulation.

Le couplage entre l'endommagement et la perméabilité des roches suscite un intérêt grandissant pour des applications industrielles. C'est le cas par exemple du problème de la fissuration autour des galeries (EDZ) de stockage de déchets radioactifs qui peut augmenter considérablement la perméabilité et constituer, de ce fait, un risque sérieux pour l'étanchéité et la sûreté de ces ouvrages. Le même phénomène peut se produire dans le béton des enceintes des centrales nucléaires ainsi que dans les roches réservoirs pétrolifères. Des recherches à caractère expérimental sur ce sujet, particulèrement au laboratoire G.3S-LMS de l'Ecole Polytechnique, ont permis de déterminer les ordres de grandeur de l'effet de l'endommagement sur la perméabilité. Mais une modélisation numérique de ce phénomène permettant d'établir des lois rhéologiques restait à faire, et c'était l'objectif principal de ce travail.

La caractérisation de deux propriétés de transport de matière a l'intérieur de milieux poreux consolides modèle (réalisés par le frittage partiel d'un empilement de billes de verre) a nécessité l'analyse de ces textures poreuses et ce par porosimétrie au mercure et/ou par analyse d'image. Un dispositif expérimental approprie a été mis au point pour étudier chacune de ces deux propriétés à savoir la conductivité électrique d'un électrolyte et l'écoulement d'un liquide. Le facteur de formation, grandeur caractérisant la conductivité électrique, est indépendant de la taille des particules composant le milieu et de l'électrolyte utilise. Il vérifie la loi d'Archie pour des porosités comprises entre 0,20 et 0,40. Il permet d'évaluer la tortuosité t du milieu, en assimilant l'espace poreux à un réseau de tubes tortueux. Deux modèles ont été utilisés pour expliciter les lois expérimentales de variation de la chute de pression subie par le liquide entre les deux extrémités du milieu en fonction de sa vitesse. Le modèle d'Ergun est inadapté à notre étude. Par contre, le modèle de Comiti lui est préférable.

A l'heure actuelle, des sédiments fluviaux et lacustres se trouvent pollués par les polychlorobiphényles (PCB), à des degrés de gravité divers. Les PCB se présentent comme des micropolluants persistants dans l'environnement et nuisibles à la santé. Un programme de recherche " PCB-Axelera " est conduit afin de pencher sur cette problématique et explorer les solutions de diminution de cette pollution. En France, l'ensemble des arrêtés préfectoraux (30/10/2002 et 14/06/2000) et le décret du 18/04/2002 considèrent les sédiments pollués comme déchets et autorisent leur stockage dans les installations de stockage de déchets si leurs teneurs en PCB sont comprises entre 0,68 ppm et 50 ppm. L'objectif principal de cette étude est de concevoir une méthodologie adaptée pour une installation de déchets non-dangereux ISDND. Les géomatériaux dans ce stockage seront utilisés comme barrière d'étanchéité pour les sédiments pollués. L'étude est composée de trois parties : Dans la première partie la caractérisation hydromécanique chimique des sédiments pollués par les PCB est présentée. La classification des sédiments est effectuée. Le comportement mécanique du sédiment est déterminé à l'aide de l'œdomètre. Quant au comportement hydromécanique il est étudié à l'aide du filtre presse et du perméamètre à paroi rigide (PPR). Cette caractérisation est indispensable pour proposer la meilleure mise en œuvre des sédiments dans le centre de stockage. Dans la deuxième partie, les performances hydraulique (perméabilité) et mécanique (gonflement) des différents géomatériaux en contact avec un fluide extrait du sédiment et d'autres fluides synthétiques contenant des PCB sont explorées. L'étude du gonflement des géomatériaux est réalisée à l'échelle du laboratoire tandis que, l'étude de sa performance hydrique est réalisée sur deux échelles centimétrique et métrique. L'étude menée sur un pilote de 1,2 m3 de volume comprenant un système de pluviation automatique a permis de suivre les tassements à l'aide d'un capteur ultrason et de mesurer la perméabilité des sédiments et des géomatériaux sous conditions réelles de stockage. Dans la troisième partie, un modèle hydraulique 1D est développé afin de suivre le front d'eau au sein d'un échantillon de sédiment soumis à une charge hydraulique représentant l'expérimentation dans un PPR. La méthode VER est utilisée. Les résultats expérimentaux montrent que le degré de contamination en PCB de ces sédiments n'influence pas le comportement hydromécanique et donc la qualité d'étanchéité des géomatériaux étudiés. Les résultats numériques de simulation du comportement hydraulique du sédiment sont en concordance avec les résultats expérimentaux. La solution de stockage hors site des sédiments pollués est alors envisageable puisqu'elle confine la pollution en éliminant le risque de contamination du milieu environnant.

L'imbibition est d'une importance fondamentale dans de nombreuses applications technologiques et intervient par ailleurs dans de nombreux phénomènes naturels (industrie textile, industrie pharmaceutique, érosion des sols ....). Malgré l'importance de ce phénomène, la description et la modélisation des mécanismes d'imbibition sont encore sujettes à discussion dans la littérature, en particulier, la prédiction des cinétiques d'imbibition à partir de la connaissance de la topographie du milieu poreux. L'objectif de ce travail de thèse a été de relier les cinétiques d'imbibition à la structure du milieu poreux imbibé. Pour cela, nous avons étudié l'imbibition capillaire (ou spontanée) en suivant une démarche expérimentale dans laquelle nous avons utilisé des substrats poreux modèles dont nous contrôlons le volume et la structure de pores. La configuration d'imbibition choisie dans ce travail est celle d'une goutte sessile (ou posée).Ces travaux ont permis de mettre en œuvre une technique de construction de pastilles macroscopiques, autosupportées et cohésives, par assemblage de microbilles de polymère. Ces systèmes modèles ont été utilisés pour étudier les cinétiques d'imbibition capillaire de liquides dans des supports tridimensionnels à structure de pores complexe en fonction de la taille des microbilles constituant la pastille poreuse, de la reconstruction thermique, de la perméabilité pour une structure poreuse bicouches et de la force capillaire. Ces études ont permis entre autres de mettre en évidence des régimes cinétiques et des transitions d'imbibition inattendues dans ces systèmes (régimes visqueux et inertiel).

En simulation de réservoir, les variations spatiales des propriétés des roches sont données par un modèle géologique défini sur une grille pouvant comporter plusieurs centaines de millions de mailles. Simuler des écoulements polyphasiques sur ce type de grille peut être très coûteux en temps de calculs. Les ingénieurs de gisement procèdent souvent à une mise à l'échelle de ces propriétés sur une grille plus grossière, utilisée ensuite pour la simulation. Cette façon de procéder ne permet pas toujours de bien prendre en compte l'influence sur l'écoulement des hétérogénéités présentes à l'échelle fine. Plusieurs méthodes multi-échelles ont donc été proposées ces dernières années afin d'obtenir un meilleur compromis entre les temps de calculs et la précision des solutions obtenues. Actuellement, la plupart de ces méthodes permettent de résoudre plus efficacement l'équation en pression. Les équations de transport sont, quant à elles, toujours résolues sur la grille fine. La première partie de cette thèse présente la mise en œuvre et les résultats de tests réalisés sur la méthode des éléments finis mixtes multi-échelles dans un environnement de calcul parallèle. Le second et principal objectif de ce travail est de proposer une nouvelle méthode multi-échelle pour la simulation du transport d'un traceur dans un milieu poreux. Cette méthode est construite à partir de résultats rigoureux d'homogénéisation en milieux périodiques. Cette thèse fournit une vitesse de convergence explicite pour ce procédé d'homogénéisation ainsi qu'une estimation d'erreur pour la méthode numérique multi-échelle associée. Des exemples numériques sont ensuite présentés.

Dans cette thèse nous mettons en oeuvre une méthode vortex hybride pénalisée (HVP) afin desimuler des écoulements incompressibles autour de corps non profilés dans des milieux complexessolides-fluides-poreux. Avec cette approche particulaire hybride, le phénomène de convection estmodélisé à l’aide d’une méthode vortex afin de bénéficier du caractère peu diffusif et naturel desméthodes particulaires. Un remaillage des particules est alors réalisé systématiquement sur unegrille cartésienne sous-jacente afin d’éviter les phénomènes de distorsion. D’autre part, les effetsdiffusifs et d’étirement ainsi que le calcul de la vitesse sont traités sur la grille cartésienne, àl’aide de schémas eulériens. Le traitement des conditions de bords aux parois de l’obstacle esteffectué à l’aide d’une technique de pénalisation, particulièrement bien adaptée au traitementde milieux solides-fluides-poreux.Cette méthode HVP est appliquée au contrôle passif d’écoulement. Cette étude de contrôleest effectuée respectivement en 2D et en 3D autour d’un demi-cylindre et d’un hémisphère parl’ajout d’un revêtement poreux à la surface de l’obstacle. La présence de cette couche poreusemodifiant la nature des conditions aux interfaces, permet de régulariser l’écoulement global etde diminuer la traînée aérodynamique de l’obstacle contrôlé. A travers des études paramétriquessur la perméabilité, l’épaisseur et la position du revêtement poreux, ce travail vise à identifier desdispositifs de contrôles efficaces pour des écoulements autour d’obstacles comme des rétroviseursautomobiles. / In this work we use a hybrid vortex penalization method (HVP) to simulate incompressibleflows past bluff bodies in complex solid-fluid-porous media. In this hybrid particle approach,the advection phenomenon is modeled through a vortex method in order to benefit from thenatural description of the flow supplied by particle methods and their low numerical diffusionfeatures. A particle remeshing is performed systematically on an underlying Cartesian grid inorder to prevent distortion phenomena. On the other hand, the viscous and stretching effects aswell as the velocity calculation are discretized on the mesh through Eulerian schemes. Finally,the treatment of boundary conditions is handled with a penalization method that is well suitedfor the treatment of solid-fluid-porous media.The HVP method is applied to passive flow control. This flow control study is realized pasta 2D semi-circular cylinder and a 3D hemisphere by adding a porous layer on the surface of thebody. The presence of such porous layer modifies the characteristics of the conditions at theinterfaces and leads to a regularization of the wake and to a decrease of the aerodynamic dragof the controlled obstacle. Through parametric studies on the permeability, the thickness andthe position of the porous coating, this works aims to identify efficient control devices for flowsaround obstacles like the rear-view mirrors of a ground vehicle.

Cette thèse est consacrée au développement d'une méthode numérique visant à modéliser des écoulements dans des milieux poreux fissurés, ainsi qu'à déterminer leur perméabilité effective à partir des avancements théoriques récents. En parallèle, elle a été aussi l'occasion de continuer sur la voie théorique et d'obtenir de nouveaux résultats sur ce plan.Les équations générales gouvernant l'écoulement dans de tels matériaux sont rappelées, et plus particulièrement, la conservation de la masse à l'intersection entre fissures est établie explicitement [132, 139]. Des solutions générales du potentiel sont proposées sous la forme d'une équation intégrale singulière décrivant l'écoulement dans et autour des fissures entourées par une matrice infinie soumise à un champ lointain [136, 139]. Ces solutions représentent le champ de pression dans le milieu infini en fonction de l'infiltration dans les fissures qui prennent en compte complètement l'interaction et l'intersection entre fissures. En considérant le problème d'une fissure super-conductrice, des solutions analytiques de l'écoulement ont été développées et ces solutions sont utilisées comme cas de référence pour valider la solution numérique. De plus, la solution théorique de ce problème dans le cas 3D permet de comparer le modèle d'écoulement de Poiseuille dans une fissure elliptique d'épaisseur nulle et le modèle d'inclusion ellipsoïdale aplatie soumise à l'écoulement de Darcy [140]. Des outils numériques ont été développés en se basant sur la méthode des équations intégrales singulières afin de résoudre les équations générales du potentiel [132, 180]. Cela permet, d'une part, de modéliser l'écoulement stationnaire dans un domaine poreux contenant un grand nombre de fissures et, d'autre part, de proposer une solution semi-analytique de l'infiltration dans une fissure isolée dépendant de la perméabilité de la matrice, de la conductivité de la fissure et de la variable géométrique de la fissure. Cette dernière est l'élément essentiel pour déterminer de la perméabilité effective d'un milieu poreux fissuré en utilisant des schémas d'homogénéisation. Ce modèle auto-cohérent révèle un seuil de percolation qui est alors applicable pour l'estimation de la perméabilité effective d'un matériau contenant un grand nombre de fissures. L'approche par sous-structuration permet l'extension de la solution générale du potentiel, écrite pour un domaine infini, à celle d'un domaine fini [181]. Une solution analytique de l'écoulement dans et autour d'une fissure partiellement saturée est établie, fondée sur la solution analytique pour la fissure super-conductrice. Celle-ci est alors utilisée pour estimer la perméabilité effective d'un milieu poreux fissuré non-saturé [141]. Le modèle de la perméabilité effective est appliqué dans le contexte du stockage géologique du CO2 en vue d'étudier le comportement d'une zone de faille constituée par un noyau argileux et des zones fissurées. La pression d'injection provoque l'augmentation de la pression interstitielle dans le réservoir. Cette surpression engendrée dans le réservoir peut affecter la perméabilité de zones fissurées ce qui conduit à des phénomènes hydromécaniques couplés. Les résultats de simulations numériques permettent d'évaluer, d'une part, le risque de la remontée de la saumure à l'aquifère supérieur, et d'autre part, le risque de l'initiation d'une rupture sur le plan de la faille

Ce travail de thèse a pour objectif de simuler numériquement des écoulements de gaz dans des matrices poreuses dont les pores sont de taille micrométrique. On étudie l'influence des phénomènes de glissement hydrodynamique qui apparaissent lorsque la dimension caractéristique de micro- conduites est caractérisée par des nombres de Knudsen compris entre Kn = 0,01 et Kn = 0,1.Le mémoire de thèse est composé de cinq chapitres suivis d'une conclusion dans laquelle nous présentons quelques perspectives pour une suite de ce travail. Le chapitre I constitue le travail préliminaire de thèse qui s'est ensuite orienté vers des approches complémentaires. Le principe des méthodes d'homogénéisation périodique est d'abord exposé. Suit une présentation de deux méthodes de résolution dans l'espace de Fourier : l'approche en déformation et l'approche en contrainte. L'extension de ces méthodes à la résolution d'écoulements régis par l'équation de Stokes est ensuite décrite. Des applications aux cas d'écoulements à travers des réseaux de cylindres, avec condition d'adhérence ou avec condition de glissement, sont ensuite discutées. Deux techniques de modélisation des phénomènes de transport dans des milieux poreux saturés par un fluide monoconstituant sont présentées dans le second chapitre. La première est basée sur la méthode des développements asymptotiques, appelée aussi méthode d'homogénéisation. On explique que le processus consiste en trois étapes : description locale, localisation et description macroscopique. La seconde technique s'appuie sur la méthode de calcul de moyennes à l'échelle d'un VER. Le point de départ de cette méthode est basé sur des théorèmes donnant les expressions des moyennes de tous les opérateurs intervenant dans une équation de transport. Après une brève présentation du logiciel commercialisé que nous avons utilisé, nous exposons les études de convergence spatiale que nous avons effectuées et nous comparons nos solutions avec des résultats de la littérature dans le chapitre III. Diverses géométries sont considérées (allant de géométries planes à des empilements 3D de cubes ou de sphères).L'effet du glissement sur la perméabilité de milieux microporeux est abordé dans le chapitre IV. Le formalisme résultant de l'homogénéisation de structures périodiques est utilisé pour simuler numériquement des écoulements isothermes de gaz dans divers empilements de complexités croissantes. Les perméabilités sont déterminées en calculant les moyennes spatiales des champs de vitesses, solutions des équations de Stokes. Les valeurs obtenues en imposant des conditions d'adhérence sont comparées à celles obtenues avec des conditions de glissement du premier ordre. Dans le chapitre V, nous présentons des solutions pour des écoulements anisothermes et étudions l'effet du glissement sur la conductivité effective de milieux microporeux 2D et 3D. Dans ce chapitre, nous résolvons les équations de Navier-Stokes et de l'énergie en imposant des conditions de symétries dans une ou deux directions. A partir des solutions locales, sont calculées les moyennes intrinsèques des champs de vitesse et de température. Nous considérons des cas pour lesquels la condition d'équilibre thermique local peut être considérée comme satisfaite et d'autres correspondant à un non-équilibre thermique (NTLE). On détermine les conductivités de dispersion en fonction du nombre du Péclet et on montre l'influence du glissement sur les composantes longitudinales et transverses pour différentes porosités et longueur de glissement. Dans les cas NLTE, le coefficient macroscopique de transfert fluide-solide est aussi calculé

Le stockage géologique du CO2 est une stratégie prometteuse pour limiter la concentration de CO2anthropique dans l'atmosphère. Les aquifères salins (AS) ont été identifiés comme des optionsviables car ils possèdent de grandes capacités potentielles de stockage. Toutefois, les processusrelatifs au piégeage du CO2 souffrent d'un manque de connaissances fondamentales car il existe peude méthodes d'expérimentation rapides et reproductibles, travaillant dans les conditions du stockagegéologique. Ainsi, nous avons développé des microréacteurs haute pression, véritables laboratoiresgéologiques sur puce (GLoCs), recréant les conditions de porosité et de perméabilité des AS pour :(i) Mesurer la solubilité du CO2 dans l'eau et les saumures via un couplage microsystèmes /spectroscopie Raman ;(ii) Etudier les mécanismes d'invasion du CO2 dans les formations géologiques, incluantnotamment les écoulements diphasiques en milieux poreux, les séparations de phases etla précipitation des carbonates.