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. / CO2 geological storage is a promising strategy to control the anthropogenic CO2 concentration in theatmosphere. Deep saline aquifers (DSA) were identified as viable options since they exhibit largestorage capacity. However, processes inherent to CO2 trapping suffer from a lack of fundamentalknowledge, since there are too few fast and reproducible experimental approaches able to work atgeological storage conditions. Therefore, to address these limitations, we have developed highpressure microreactors, so-called “geological labs on a chip” - GloCs – allowing mimicking porosityand permeability conditions of DSA for:(i) Measuring solubility of CO2 in water and brine through the combination of microsystemsand confocal Raman spectroscopy,(ii) Studying invasion mechanisms of CO2 in geological formations, including in particularbiphasic flows in porous media, phase separation and carbonates precipitation.

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.

Les premières centaines de mètres du sous-sol sont au coeur d'enjeux sociétaux importants liés aux fluides auxquels la société doit encore affronter, tels que la détéction et la surveillance des ressources en eau et de la pollution, ainsi que le rôle des fluides dans les différents types d'évaluation des risques ou en géotechnique. Cette gamme de profondeur est marquée par une forte complexité et une forte hétérogénéité aussi bien structurale et lithologique qu'en termes de mélanges de phases fluides. Parmi les méthodes géophysique disponibles, le phénomène de couplage sismo-électrique (et électro-sismique) a le potentiel d'émerger comme une nouvelle technique d'imagerie haute résolution, naturellement sensible aux contrastes des fluides. Deux signaux sismo-électrique ont été prédits et parfois observés: les signaux co-sismique éléctrique à forte amplitude et les faibles perturbations électromagnétiques générées en profondeur lorque des ondes sismiques traversent une interface. Ce deuxième phénomène est le plus prometteur, en termes d'imagerie en raison de sa sensibilité aux contrastes originaux par rapport à la réflexion sismique. L'émergence d'une méthode sismo-électrique a cependant été ralentie en raison des difficultés dans l'enregistrement de ces faibles signaux sismo-électrique. Cette thèse, basée sur la combinaison d'approches expérimentales et numériques, vise à aborder des questions pratiques concernant l'acquisition et le traitement de données sismo-électrique et électro-sismique.Dans la première partie, afin d'augmenter le rapport signal-sur-bruit (S/N) des signaux sismo-électriques, nous étudions théoriquement et expérimentalement leur sensibilité à différents arrangements d'électrodes. Pour cela, nous avons développé une théorie des filtres conçue pour évaluer l'influence d'ensemble d'électrodes pour l'acquisition des signaux sismo-électriques. Cette théorie a été confrontée avec succès à des données expérimentales sismo-électriques acquises en utilisant diverses dispositions d'électrodes (c'est-à-dire l'espacement entre les électrodes, le nombre des électrodes et la vitesse sismique apparente) et à des simulations numériques des formes d'onde complète. En particulier, cette approche explique la difficulté d'enregistrement des ondes converties sismo-électromagnétique générées en profondeur converties à une interface, lorsque la configuration dipôle classique est utilisé comme récepteur. Nous encourageons ensuite l'utilisation d'une configuration de multi-électrodes dans des mesures sismo-électriques, puisque: 1) il réduit le bruit électrique omniprésent causé par l'homme, et 2) ses propriétés de filtrage amplifient les ondes électromagnétiques provenant des interfaces en profondeur par rapport aux événements co-sismiques dominants.Dans la deuxième partie, nous nous concentrons sur les phénomènes réciproques électro-osmotique, c'est-à-dire la génération de signaux sismiques par une source électrique dans des milieux saturés poroélastiques. Nous avons d'abord effectué des études numériques pour évaluer la nature et les propriétés des différentes arrivées sismiques générées par une source électrique dans un matériau homogène et hétérogène, à l'échelle du terrain. Ensuite, nous avons essayé d'acquérir des données de laboratoire électro-sismique dans un matériau homogène, données que nous avons de nouveau confrontées à des calculs numériques de forme d'onde complète. Les résultats montrent que les ondes électro-sismiques peuvent être observées et qu'elles sont générées localement autour de la source électrique. Nous montrons enfin que ces signaux électro-sismiques sont influencés par la conductivité des fluides. / The first few hundred meters of the subsurface is the seat of important issues related to fluids which society still has to face, such as the detection and monitoring of water ressources and pollution as well as the role of fluids in various types of hazard assessment or in geotechnics. This depth scale is highly complex and heterogeneous in terms of lithology and fluid content. Among the available geophysical methods, seismo-electric (or electro-seismic) coupling phenomena have the potential to emerge as a new high-resolution imaging technique, naturally sensitive to fluid contrasts. Two seismo-electric signals have been predicted and sometimes observed: co-seismic electric signals (high amplitude) and weak electromagnetic (EM) disturbances generated at depth when seismic waves cross an interface. This second phenomenon is the most promising in terms of imaging due to its sensitivity to original contrasts compared to seismic reflection. The emergence of the seismo-electric method was however slowed down due to difficulties in recording these weak seismo-electric signals. Based on the combination of experimental and numerical approaches, this PhD thesis aims at addressing practical questions concerning the acquisition and processing of seismo-electric and electro-seismic data.In the first part, in order to increase the signal-to-noise (S/N) ratio of seismo-electric signals, we theoretically and experimentally study their sensitivity to various electrode arrangements. For this, we developed a filter theory designed to assess the influence of complex electrode arrays for seismo-electric signals acquisition. This theory was successfully confronted to experimental seismo-electric data acquired using various electrode arrangements (namely the spacing between the electrodes, the number of the electrodes and the apparent seismic velocity) and to full waveform numerical simulations. This combined analysis shows that electrode configuration properties strongly influence seismo-electric amplitudes and waveforms. In particular, this approach explains the difficulty in recording depth-generated seismo-electromagnetic waves converted at an interface, when the conventional dipole configuration is used as receivers. We then promote the use of a 3-electrode configuration in seismo-electric measurements, since: 1) it reduces ubiquitous man-made electric noise, and 2) its filtering properties do amplify the EM waves originating from interfaces at depth with respect to dominant coseismic events.In the second part, we focus on the reciprocal electro-osmotic phenomena, i.e. the generation of seismic signals by an electric source in poroelastic saturated media. We first performed numerical investigations to assess the nature and properties of the different seismic arrivals generated by an electrical source in a homogeneous and an heterogeneous material, at the field scale. Then we tried to acquire electro-seismic laboratory data in a homogeneous material, data that we again confronted with numerical full waveform computations. The results show that electro-seismic waves can be observed and that they are locally generated around the electric source. We finally show that these electro-seismic signals are influenced by the fluid conductivity.

L’accident de Fukushima a démontré qu’aujourd’hui encore, malgré les progrès déjà réalisésdans le domaine de la sûreté nucléaire, une interruption prolongée du circuit de refroidissement primaire du réacteur (plusieurs semaines ici) était possible. La France s’est par conséquent vue contrainte de réévaluer le niveau de sûreté de ses centrales. Plus spécifiquement, le cas le plus défavorable qu’EDF considérait jusqu’alors, qui consistait en un arrêt total du système de refroidissement primaire de 24h, a été réévalué à deux semaines. Ce changement d’échelle temporelle a introduit des problématiques de fluage, d’évolution hygrométrique du béton ainsi que de flux de vapeur – celle-ci transportant les radionucléides – jusque-là laissées de côté. De ces considérations est né le projet ANR/RSNR MACENA (MAîtrise du Confinement d’une ENceinte en Accident), qui vise à évaluer le confinement d’une enceinte soumise à une température de 180°C et à une pression de 5 bar pendant deux semaines. Ce travail s’inscrit dans la lignée de celui entamé par Khaddour [2014] dans le but de parvenir à une meilleure prédiction des propriétés de transport des milieux poreux en se basant uniquement sur une caractérisation indirecte, la distribution des tailles de pore obtenue grâce à la technique de porosimétrie par intrusion de mercure. Si l’aptitude du modèle initial à estimer correctement les perméabilités intrinsèques a été confirmée sur des matériaux représentant une large gamme de perméabilités, il n’était pas, en l’état, capable de prédire les perméabilités relatives au gaz et au liquide. C’est pourquoi une approche dite redistributive a été développée, laquelle a montré de bons résultats sur diverses PSD. Enfin, un nouveau modèle a été développé pour essayer de s’affranchir de cette approche tout en améliorant la capacité à déterminer les perméabilités relatives. Bien qu’efficace sur des PSD monomodales, étroites ou larges, il a montré ses limites sur une PSD réelle bimodale. En parallèle, une étude de faisabilité a exhibé la possibilité de suivre un front de saturation à l’intérieur d’un matériau poreux, ce qui permettrait à terme de valider l’inclusion d’autres phénomènes dans le modèle, liés notamment à la pression de perçage. / The Fukushima catastrophe that struck Japan in 2011 demonstrated that despite significant progress in the field of nuclear safety a prolonged reactor primary cooling circuit breakdown was possible (several weeks in this specific case). With 4 nuclear power plants located on its shoreline, France therefore needed to reassess the safety level of its facilities. More specifically, the worst case scenario considered up until that point by EdF — which consisted in a 24h breakdown of the primary cooling system — was revised up to two weeks. This time-scale shift induced creep, drying and vapour flow problems previously left aside. Thus came to be the ANR/RSNR MACENA (MAîtrise du Confinement d’une ENceinte en Accident) project, which aims at bettering the tightness assessment of a nuclear containment vessel submitted to a temperature of 180°C and to a pressure of 5 bar for two weeks. This work falls in line with that initiated by Khaddour [2014], who set out to better predict porous materials’ transport properties based solely on an indirect characterisation of their topology, namely Pore Size Distribution (PSD) which is obtained via an experimental technique called Mercury Intrusion Porosimetry (MIP). The initial model’s ability to correctly estimate intrinsic permeabilities was confirmed on several different materials whose intrinsic permeabilities span several orders of magnitude. However, it was not fit to accurately account for gas and liquid relative permeabilities. This led to the introduction of a so-called redistributive approach, which yielded better results on various PSD and corresponding experimental datasets. Finally, a new model was developed to try and avoid said approach while bettering relative permeability predictions. Although it behaved well with monomodal pore size distributions, be they wide or narrow, it fell short when applied to a real bimodal PSD. In parallel, an experimental feasibility study demonstrated the possibility of tracking a saturation front within a porous material, which should ultimately allow for the implementation of several other phenomena into the model, linked among others to breakthrough capillary pressure.

Les effets capillaires sont très courant dans la Nature. Dans le contexte du séchage de milieux poreux dont la taille de pore est dans la gamme micromètre-millimètre, ils jouent un effet dominant en contrôlant la répartition des phases (liquide-vapeur) dans l’espace poral, au fur et a` mesure que le séchage se produit. L’idée du présent travail est d’étudier le séchage d’un fluide pur et mouillant dans des micromodèles, c’est-à-dire des milieux poreux modèles quasi-2D et micro-fabriqués. Nous présentons des résultats obtenus pour différentes géométries. Typiquement, les micromodèles utilisés sont constitués de réseaux de cylindres pris en sandwich entre deux plaques. La distribution des phases et le taux d’évaporation dans de tels micromodèles peuvent être aisément mesurés par visualisation directe puis traitement d’images.En jouant sur l’arrangement spatial des cylindres, on obtient dans un premier temps des micromodèles pour lesquels le taux de séchage est quasi-constant, depuis le début de l’expérience de séchage jusqu’à l’évaporation totale du liquide saturant initialement le système. Typiquement, cette situation est obtenue quand la taille des pores décroît en allant du centre du micromodèle vers sa périphérie (les micromodèles sont axisymmétriques). Au contraire, quand la taille des pores croît du centre vers la périphérie, l’invasion d’un front de séchage stable est observée, d’ou` un temps de séchage total bien supérieur.Nous avons aussi réalisé un autre type de microsystèmes, au sein duquel les cylindres sont arrangés en spirale de Fibonacci, en nous inspirant de motifs observés en phillotaxie. Dans de tels systèmes, des films liquides épais se développent le long des spirales, au cours du séchage, et jouent un rôle crucial dans la cinétique d’évaporation. Cette situation rappelle celle déjà étudiée par Chauvet dans des tubes capillaires de section carrée. Cependant, elle est plus complexe, de par la nature poreuse du micromodèle (alors qu’un tube capillaire, tel qu’étudié par Chauvet, peut être vu comme un pore unique) et parce que les films liquides y ont une forme plus complexe. Pour de tels systèmes, nous présentons des résultats expérimentaux quantifiant l’effet des films liquides sur la cinétique de séchage, en lien avec des prédictions théoriques issues d’un modèle de séchage visco-capillaire. Un tel modèle nécessite l’utilisation du logiciel Surface Evolver pour modéliser la forme des films liquides, couplée avec des simulations directes de l’écoulement de Stokes dans les films liquides, pour y calculer la résistance visqueuse a` l’écoulement induit par l’évaporation.Enfin, dans un dernier chapitre, plusieurs expériences d’évaporation sont conduites sur des micromodèles déformables. Des effets élasto-capillaires peuvent en effet induire des changements de géométrie de l’espace poral en cours d’évaporation, ce qui, comme vu précédemment, peut affecter la distribution des phases et la cinétique de séchage. / Capillarity is a common phenomenon encountered in Nature. In the context of the drying of porous media with pore size in the micrometer-millimeter size range, capillary effects play a dominant role in controlling the phases (liquid or vapor) distribution in the pore space as drying occurs. The basic idea of the present work is to study the drying of pure, wetting fluids in micro-fabricated, quasi-2D, model porous media (hereafter called micromodels). We present results obtained for different micromodel geometries. Typically, the micromodels used consist of arrangements of cylinders sandwiched between a top and bottom plate. Phases distribution and evaporation rates in such micromodels can easily be measured by direct visualizations and subsequent image processing.By tuning the cylinders pattern, one can first obtain micromodels for which the drying rate is almost constant, from the beginning of the drying experiment to the total evaporation of the liquid initially filling the system. Typically, this situation is obtained when the pores size decreases from the micromodel center to the periphery (the micromodels are axisymmetric). On the contrary, when the pores size increases from the center to the periphery, invasion of a stable drying front is observed, resulting in a much longer total drying time.We also designed another type of micromodel where the cylinders are arranged in a Fibonacci spiral pattern, a design inspired by phyllotaxic structure. In such systems, thick liquid films develop along the spirals during drying and play a key role in the drying kinetics. This situation is reminiscent of that already studied by Chauvet in capillary tubes with square cross-sections. However, it is more complex because of the porous nature of the micromodel (whereas a single capillary tube, as studied by Chauvet, can be viewed as a unique pore), and because of the much more complex liquid films shapes. For such systems, we present some experimental results on the liquid films effects on the drying kinetics, together with theoretical prediction, based on a visco-capillary drying model. Such a modelling requires the use of the Surface Evolver software to model the film shape, coupled with DNS simulations of the Stokes flow within the liquid films to compute the viscous resistance to the evaporation-induced flow.Finally, as a last part of this thesis, several evaporation experiments performed on deformable micromodels are presented. This preliminary work aims at reaching a situation where elasto-capillary effects modify the pore space geometry during evaporation. This, as seen above, should in turn alter the phase distribution during evaporation and the drying kinetics.

Cette thèse a pour objectif de modéliser le couplage fluide-structure pouvant survenir dans les faisceaux tubulaires des réacteurs nucléaires. Leurs simulations numériques directes étant hors de portée, on met en œuvre une approche multi-échelle: il s'agit de tirer profit du coût modeste d'une description macroscopique, et à la fois, de la précision des informations microscopiques. Vis-à-vis des modèles existants, le travail de développement se focalise sur la prise en compte de la convection dans le calcul des champs hydrodynamiques, mais surtout, sur la possibilité de restituer des réponses vibratoires variées au sein d'un même faisceau. L'homogénéisation aboutit à un système d'équations gouvernant les Interactions Fluide-Solide à une échelle macroscopique. Ces équations sont couplées par une source en quantité de mouvement, traduisant les charges hydrodynamiques exercées sur une structure donnée. Cette force à modéliser représente une loi de fermeture du problème homogénéisé, mettant en jeu des coefficients a priori inconnus. Une méthode d'estimation est proposée à partir des champs microscopiques obtenus par simulation directe sur un domaine réduit et représentatif du large système de référence. Les capacités prédictives du modèle homogénéisé sont évaluées en comparaison avec des données de référence, issues de calculs numériques directs microscopiques. Chaque système considéré présente une variété de réponses en déplacement que le modèle homogénéisé restitue avec un accord satisfaisant. Cette approche multi-échelle semble être un bon compromis entre le coût des réalisations numériques et la précision attendue des données vibratoires et hydrodynamiques. / Vibration of tubes arrays is a matter of safety assessments of nuclear reactor cores or steam generators. Such systems count up thousands of slender-bodies immersed in viscous flow, involving multi-physics mechanisms caused by nonlinear dynamic interactions between the fluid and the solid materials. Direct numerical simulations for predicting these phenomena could derive from continuum mechanics, but require expensive computing resources. Therefore, one alternative to the costly micro-scale simulations consists in describing the interstitial fluid dynamics at the same scale as the structures one. Such approach rely on homogenization techniques intended to model mechanics of multi-phase systems. Homogenization results in coupled governing equations for the fluid and solid dynamics, whose solution provides individual tubes displacements and average fluid fields for each periodic unit cell. An hydrodynamic force term arises from the formulation within this set of homogenized equations: it depends on the micro-scale flow in the vicinity of a given tube-wall, but needs to be estimated as a function of the macro-scale fields in order to close the homogenized problem. The fluid force estimation relies on numerical micro-scale solutions of fluid-solid interactions over a tube array of small size. The multi-scale model is assessed for arrays made up of hundreds tubes, and is compared with solutions coming from the numerical micro-scale simulations. The macro-scale solution reproduces with good agreement the averaged solution of the micro-scale simulation, indicating that the homogenization method and the hydrodynamic force closure are suitable for such tube array configurations.

Les réactions de catalyse hétérogène impliquant un gaz et un liquide sont mises en oeuvre dans des réacteurs à lit fixes. Ces réacteurs peuvent être assimilés à un milieu poreux. La nature complexe de ce milieu rend la compréhension des interactions entre phases difficile, et nécessite une étude exhaustive à l’échelle globale et locale afin d’identifier les paramètres clés de l’hydrodynamique, des transferts de chaleur et de matière. Nous avons donc développé une cellule miniaturisée bidimensionnelle transparente, qui permet l’observation directe des écoulements avec une très bonne résolution spatiale et temporelle. En faisant varier le débit total et le rapport des débits gaz/liquide, nous avons mis en évidence l’apparition des régimes ruisselant et pulsé, observés par ailleurs dans des systèmes tridimensionnels. Grâce à une méthode d’analyse d’image, nous sommes capables de quantifier et cartographier la saturation liquide locale apparente et la morphologie des phases. L’analyse des variances a permis d’étudier les transitions de régimes pour différentes propriétés de la phase liquide. Cette approche a permis de comparer avec la littérature, mais aussi de s’intéresser aux mécanismes de génération et propagation des instabilités lors des transitions. Il a été mis en évidence que l’apparition des instabilités responsables de la déstabilisation du régime ruisselant intervient pour un nombre deWeber liquide fixé, indiquant que le régime pulsé apparait suite à la déstabilisation des interfaces gaz /liquides par les forces inertielles. Enfin, une étude préliminaire des transferts thermiques dans la cellule a été réalisée. La cellule est utilisée pour réaliser la réaction exothermique d’hydrogénation de l’alpha-methylstyrène. Un modèle simple de transferts thermique a été utilisé pour caractériser l’augmentation de température dans la cellule. Bien que cette approche ne permette pas des mesures quantitatives, elle ouvre à la perspective de suivi de transferts thermiques par caméra infra-rouge. / Heterogeneously catalyzed reactions involving a gas and a liquid phase are frequently achieved in fixed bed reactors. These reactors can be described as a porous medium. The complex nature of this medium makes the understanding of the interplay between phases difficult, and requires a thorough study at the global andlocal scale to identify the key parameters of hydrodynamics, mass and heat transfers. Therefore, we have developed a miniaturized two-dimensional system that is transparent to allow the direct observation of the flow with very high spatial and temporal resolution. While varying the total flow rate and the gas/liquid flow rate ratio, we observe the appearance of the trickle and the pulse flow regime, which can be observed in threedimensional beds. Thanks to some image analysis techniques, we are able to quantify and to map the local apparent liquid saturation and the morphology of the phases. Variances analysis allowed the study of the transition for different liquid properties. This approach allowed the comparison with the existing state of art, but also the study of the onset and propagation mechanisms of the instabilities during the transition. We report that the onset of instabilities responsible for the destabilization of the trickle flow regime occurs at a fixed Weber number. This indicates that the pulsed regime is due to the destabilization of the gas/liquid interface by inertial forces. Finally, a preliminary study of thermal transfers in the device was realized. The device was used to perform the exothermic hydrogenation of alpha-methylstyrene. A simple model is used to characterize the temperature increase in the device. Although this approach does not allow quantitative measurements, it opens up the perspective of monitoring thermal transfers with an infra-red camera.

Le coefficient d'électrofiltration est simulé par méthode Lattice Boltzmann dans un chenal 2-D sur une grande gamme de salinité. L'influence de la permittivité et de la viscosité est discutée. La validité de l'équation d'Helmholtz Smoluchowski à de forts potentiels zeta est évaluée. Un modèle de conductivité intrinsèque est développé en prenant en compte les variations locales de conductivité, qui ont un impact significatif en la présence d'espèces polyvalentes. Étendu aux conditions non saturées, l'algorithme montre que la densité de charge électrique associée à l'interface eau-air est une composante clé. Le coefficient présente une attitude non monotone, avec une augmentation par rapport à l'état saturé. L'amplitude de cette augmentation dépend de l'état dynamique des bulles, mobiles ou piégées. L'aspect transport multiphasique est complété par une étude numérique de l'impact de la forme des échantillons sur la mesure des lois reliant saturation et pression capillaire en hydrologie. / The electrokinetic coefficient is simulated in a large range of salinities using the Lattice Boltzmann method in a 2-D channel. The effect of permittivity and viscosity is discussed. The validity of the Helmholtz Smoluchowski equation using strong zeta potentials is assessed. A model of bulk fluid conductivity is derived, taking into account the local variations of conductivity which have a significant impact in the presence of polyvalent counterions. Extended to unsaturated conditions, the model shows that the electrical charge density associated to the air-water interface is a key component. The coefficient shows a non monotonous behaviour, with an enhancement compared to the saturated state. The magnitude of this enhancement depends on the dynamic state of the bubbles, moving or entrapped. The multiphase transport aspect is associated to a numerical study of the influence of the sample geometry on the measurement of the capillary pressure / saturation relationships used in hydrology.

La liquéfaction des sols s'observe lors de forts séismes dans des zones saturées en eau et peut causer un enfoncemenl des bâtiments dans le sol. On propose ici un nouveau modèle qui explique de nombreux cas de liquéfaction incompris jusque là. Un modèle analytique de milieu granulaire a été conçu avec une sphère modélisant un bâtiment. Le milieu est soumis à une oscillation horizontale et on caractérise l'état liquéfié par l'amplitude d'enfoncement de la sphère. On suppose que la liquéfaction a lieu lorsque la secousse permet aux grains du milieu de glisser les uns sur les autres. La fenêtre d'accélération permettant ce glissement dépend de la hauteur d'eau dans le milieu, du coefficient de frottement et de la densité du matériau. Les expériences et les simulations numériques réalisées confirment nos prédictions précédentes, caractérisent la micromécanique des milieux liquéfiés, expliquent la vitesse de pénétration de la sphère et ouvrent des pistes pour l'étude des sables mouvants. / Sail liquefaction happens during important earthquakes in water-saturated zones and can make the buildings sink in the soil. A new model explaining many occurences of liquefaction we cannot understand yet is presented in this paper. A new analytical model of granular medium, with a sphere representing a building, has been developed. A horizontal oscillation is applied on the medium. The penetration rate of the sphere then caracterizes the liquefaction state of the modeled soil. We suppose that liquefaction happens when the shaking makes the grains slide against one another. The interval of acceleration allowing this sliding depends on the water height in the medium, on the friction coefficient and on the material density. Experiments and numerical simulations have been carried out. They confirmed our previous predictions, they caracterize the micromecanics of liquefied media, they explain the penetration speed of the sphere and they suggest perspectives to study quicksands.

Afin de mieux comprendre les spécificités de l'écoulement des fluides en seuil en géométries confinées, nous avons opté pour une approche multi-échelle expérimentale et/ou numérique dans des milieux poreux complexes et modèles. Nous montrons qu'il est possible d'utiliser la RMN pour visualiser des écoulements de fluides à seuil en géométrie complexe. Dans un milieu poreux, il est également possible de mesurer la distribution statistique des vitesses, ceci sans problème de résolution spatiale, grâce à la méthodologie de réglage d'une expérience d'injection sous IRM que nous avons mise en place. A l'aide de ces techniques, nous montrons que l'écoulement d'un fluide à seuil dans un pore modèle (une expansion-contraction axisymétrique) se localise dans la partie centrale du pore, dans le prolongement du tube d'entrée, tandis que les régions extérieures restent dans le régime solide. Des simulations numériques confirment ces résultats et montrent que la localisation de l'écoulement provient du confinement engendré par la géométrie. A l'inverse, nous montrons que pour un fluide à seuil s'écoulant dans un milieu poreux réel (en trois dimensions), il n'existe pas de zones au repos. De plus, la distribution de vitesse est identique à celle d'un fluide newtonien. Une analyse de ces résultats nous permet de prédire la forme de la loi de Darcy pour les fluides à seuil et de comprendre l'origine physique des paramètres déterminés par des expériences d'injection « macroscopiques » / To better understand the specifics of the flow of yield stress fluids in confined geometries, we opted for a multi-scale experimental and / or numerical approach in complex and model porous media. We show the usefulness of NMR for the study of yield stress fluid's flows in complex geometry. In a porous medium, we can also measure the true probability density function of fluid velocities without spatial resolution problem thanks to a complete optimisation of the design process of a NMR-PGSE experiment. Using these measurement technics, we find that the flow of a yield stress fluid in a model pore (an axisymetric expansion-contraction) is localised in the central part of the pore, i.e. in the continuity of the entry duct, and the external region stay at rest in the solid regime. Numerical simulations confirm those results and point out that the flow localisation is due to the confinement caused by the geometry. On the contrary, no region at rest exists for a yield stress fluid flowing through a real porous media (in 3D). Furthermore, the velocity distribution is the same as a newtonian fluid. The analysis of the results makes it possible to deduce the form of the Darcy's law for yield stress fluids and provides an insight in the physical origin of the coefficients found by “macroscopical” injection experiments