Cette thèse est consacrée è la modélisation et a la simulation de l'ensemble PISTON /SEGMENT/ CHEMISE, d'un moteur thermique. Dans un moteur, le piston muni de segments est en mouvement relatif dans la chemise. Ce système est lubrifié par projection d'huile, les segments et la surface de la chemise "contrôlent" la quantité de lubrifiant restant sur la paroi. Pour modéliser ce phénomène, le modèle conservatif parabolique-hyperbolique d'Elrod-Adams (P − θ), prenant en compte la cavitation (présence de gaz dans le fluide, aspect multiphasique de l'écoulement), est couplé a un modèle de contact, celui de Greenwood-Tripp (basé sur une approche statistique de la surface). En se restreignant a l'écoulement et au contact, nous négligeons beaucoup de phénomènes physiques pouvant intervenir comme le grippage ou les effets thermiques. Dans une première partie, les phénomènes mécaniques sont étudiés, de la réalisation de la surface des chemises è la modélisation de la physique de la lubrification du contact segment/chemise. Cette étude expose les simplifications effectuées et le choix des conditions aux limites pour la prise en compte de la cavitation en tant que problème a frontière libre, notamment au niveau du débit d'entrée pour simuler un fonctionnement "normal". Dans une deuxième partie, la mise en oeuvre et la modification de l'algorithme "classique" pour des surfaces mesurées sont décrites en détail. Nous présentons les performances simulées de différentes surfaces, et nous les comparons è des essais effectués par Renault. Dans une troisième partie, des modifications sont introduites au modèle (P − θ) utilisé précédemment afin de rendre l'écoulement cohérent avec le modèle de Navier-Stokes. Bien que ce nouveau problème puisse être posé en dimension 2, l'étude mathématique ne concerne que la dimension 1. On se ramène a un système dynamique dont les inconnues sont les frontières libres. L'unicité de la solution impose l'introduction d'un paramètre artificiel. Le comportement local et/ou global en temps est déterminé par les paramètres géométriques et les conditions aux limites. L'introduction de deux paramètres supplémentaires permet de démontrer l'existence et l'unicité de la solution. Dans une dernière partie, nous comparons le modèle de cavitation proposé avec le modèle P − θ.

Cette thèse présente une étude expérimentale visant à caractériser les propriétés hydrauliquesd'écoulements à surface libre de fluides complexes: fluides viscoplastiques et mélangesfluides viscoplastiques-grains. Ce travail est motivé par la nécessité d'améliorer lesconnaissances sur les propriétés des écoulements de laves torrentielles, dans un contexte deprotection contre les risques naturels en montagne. Nous avons développé un dispositif delaboratoire original, un canal à fond mobile, permettant de générer des coulées gravitaires quisont globalement stationnaires dans le référentiel du laboratoire. Les coulées présentent unfront très raide, suivi à l'amont d'une zone d'épaisseur uniforme. Les fluides viscoplastiquesutilisés (carbopol et kaolin) ont été choisis de manière à modéliser la matrice des couléesboueuses naturelles en prenant en compte les critères de similitude adaptés. Ces fluidessuivent une loi de comportement de Herschel-Bulkley, et leurs propriétés rhéologiques ont étédéterminées au moyen d'un rhéomètre de laboratoire en accordant un soin particulier àl'établissement des incertitudes associées. Nous présentons des mesures précises de l'évolutionde l'épaisseur des coulées dans la zone uniforme en fonction de la vitesse du fond du canal etde la pente. Ces résultats sont comparés à des prédictions théoriques correspondant à unécoulement permanent et uniforme d'un fluide de Herschel-Bulkley, et nous montrons quel'accord est satisfaisant moyennant les incertitudes sur les mesures rhéologiques. Nous avonségalement analysé la forme du front des coulées et, là-aussi, l'accord avec des modèleshydrauliques (modèles de type Saint-Venant) est bon. Ces deux études permettent de validerles prédictions hydrauliques obtenues à partir de lois de comportement mesurées aurhéomètre. Enfin, dans une dernière partie, nous présentons une étude préliminaire de ladynamique de particules rigides isolées placées au sein de la coulée, en nous intéressant auxvitesses longitudinales et aux vitesses de chute pour différentes tailles et différentes densitésde particules.

L'objectif de cette thèse est l'accélération des solveurs de Gradient Conjugué avec déflation utilisés pour la résolution de l'équation de Poisson pour la pression, dans le cas de la simulation d'écoulements à faible nombre de Mach sur des maillages non structurés. Une méthode de redémarrage basée sur une estimation de l'effet des erreurs numériques a été mise en œuvre et validée. Par la suite, une méthode à trois niveaux de maillage a été créée, et deux techniques ont dû être développées pour réduire le nombre d'itérations sur les niveaux grossiers : l'une permet la création de solutions initiales grâce à une méthode de projection adaptée, l'autre consiste en une adaptation du critère de convergence sur les niveaux grossiers. Les résultats numériques sur des simulations massivement parallèles montrent entre autres une réduction considérable du temps de calcul global. D'autres pistes de recherche sont introduites, notamment concernant l'équilibrage dynamiques de charge de calcul.

Ce mémoire présente quelques contributions à la prise en compte de diverses conditions limites ou d'interface lors de la résolution de problèmes par la méthode des éléments finis. Diverses techniques sont passées en revue, avec un focus sur celle de Nitsche. Les problèmes traités proviennent de la mécanique des solides et des fluides, comme par exemple l'interaction fluide-structure ou le contact.

Une capsule est une goutte de liquide enveloppée par une membrane fine et déformable. Les propriétés mécaniques de la membrane sont essentielles pour le mouvement de la capsule. L'analyse de l'écoulement d'une suspension de capsules dans un canal microfluidique au moyen d'un modèle mécanique est une technique permettant de déterminer les propriétés élastiques de la membrane. Un modèle numérique tridimensionnel a été développé pour résoudre ce problème d'interaction fluide-structure en écoulement confiné. Il couple une méthode des intégrales de frontières pour les écoulements des fluides et une méthode éléments finis pour la déformation de la membrane. Le modèle est utilisé pour étudier l'écoulement d'une capsule initialement sphérique dans des canaux de différentes sections. Dans un canal cylindrique, on montre que l'effet de confinement du canal conduit à la compression de la capsule. Cela engendre la formation de plis sur la membrane autour de l'axe de l'écoulement, phénomène également observé expérimentalement. Dans un canal de section carrée, les effets de la loi constitutive de la membrane, du rapport de taille et du débit d'écoulement sur la déformation de la capsule sont systématiquement étudiés. La comparaison entre les résultats expérimentaux et numériques nous permet de déduire les propriétés mécaniques de la membrane d'une population de capsules artificielles. Ce travail démontre la faisabilité de la mesure de propriétés mécaniques d'une membrane en utilisant une technique microfluidique en canal carré. Il pourrait être étendu par l'étude d'écoulements instationnaires dans un canal de section variable ou avec bifurcations.

Les caractéristiques d’un nuage de poussières avant l’activation de la source d’inflammation ont une grande influence sur la sensibilité et la sévérité de son explosion. Pour cette raison, le procédé de dispersion d’une poussière combustible qui se développe à l’intérieur des chambres d’explosion standardisés telles que le tube de Hartmann modifié et la sphère 20 L a été décrit par deux approches complémentaires. Dans un premier temps, une caractérisation expérimentale de l’évolution du nuage de poussière a identifié les niveaux de turbulence de l’écoulement de gaz, ainsi que les principales variations de la distribution de la taille des particules combustibles. Ces résultats ont été complétés par l’identification des niveaux de ségrégation de la poudre dispersée à l’intérieur des chambres d’explosion. Par la suite, une série de simulations CFD basée sur une approche Euler-Lagrange a été développée pour prédire le comportement du nuage. Cette étude a été réalisée grâce à la caractérisation des principaux mécanismes d’interaction tels que le transfert de quantité de mouvement (couplage bidirectionnel) et le phénomène de fragmentation des poudres. Deux types de particules microniques ont été spécifiquement étudiés : une poussière métallique, l’aluminium et un composé organique, l’amidon de blé. Ces résultats ont démontré qu’il est possible de définir ab initio, par simulation numérique, les conditions de tests les plus pertinentes (les plus réalistes ou les plus pénalisantes) en vue de la quantification de ces risques majeurs / The pre-ignition stage of an explosibility test determines the ignitability of a dust cloud as well as the main characteristics of the flame propagation. For this reason, the dispersion process of a combustible dust that develops inside the explosion chambers of the modified Hartmann tube and the 20 L sphere has been described by two complementary approaches. Initially, an experimental characterization of the evolution of the dust cloud identified the evolution of the turbulence levels of the gas flow along with the main variations of the particle size distribution of the combustible dust. These results were complemented by the identification of the segregation levels of the dispersed powder inside the explosion chambers. Thereafter, a set of CFD simulations based on the Euler-Lagrange formulation was developed to predict the behavior of the combustible dust cloud. This study was accomplished through the characterization of the main interaction mechanisms such as the momentum exchange (two-way coupling) and the dust fragmentation phenomenon. In this manner, the research work constituted for this thesis allowed determining the most appropriate conditions to ignite a dust cloud formed by a metallic (Aluminum micro-Al 42) or organic powders (wheat starch) in a case study. Thus, this research work presents a methodology that can be extended for the analyses of combustible dusts and the further development of the standard test methods

Cette thèse traite de la modélisation des réacteurs chimiques par la « modélisation compartimentale », qui consiste à diviser le système en un réseau d’une dizaine à quelques centaines de volumes interconnectés, appelés compartiments. La structure du réseau est déduite à partir d’informations provenant d’expériences de traçage, d’informations techniques sur le réacteur chimique, de simulations de mécanique des fluides numérique et des objectifs de la modélisation. Cette méthode procure un bon compromis entre temps de calcul et finesse des résultats. Quand ils sont correctement menés, les modèles à compartiments donnent des prédictions similaires, en termes de réactions chimiques, à ceux issus des simulations de mécanique des fluides numérique réactive avec un temps de calcul plus court et une représentation physique plus concrète du comportement du réacteur. Chaque étude issue de la littérature est consacrée à un réacteur spécifique avec une approche particulière qui ne peut pas être directement transposée sur un autre réacteur. L’objectif de cette thèse est d’apporter une contribution au développement d’une méthodologie la plus générale possible et de développer un outil de génération automatique et de résolution du système d’équations différentielles qui doit être résolu. Dans le premier chapitre, un état de l’art est réalisé, définissant le champ d’application de notre méthode, dans le but d’identifier les méthodes de découpage les plus pertinentes et les différentes méthodes pour calculer les échanges entre les compartiments. Dans un second chapitre, une méthode générale pour de la modélisation compartimentale est développée. Une approche polyvalente est proposée, consistant à découper le réacteur en tranches identiques. Le calcul des échanges entre compartiments, dus à la convection et la turbulence, est présenté en détail, avec la description des trois méthodes de calcul des échanges turbulents. Une interface a été développée permettant de construire n’importe quel réseau de compartiments. À partir de cette interface, les équations sont écrites et automatiquement résolues. La méthode est appliquée dans un troisième chapitre sur un cas défavorable au découpage en tranches. Cela a permis de tester les limites de cette approche. En particulier, deux points ont été étudiés : (1) l’applicabilité du découpage en tranches identiques et (2) la comparaison entre les méthodes de calcul des échanges turbulents. Le premier test a prouvé la robustesse de l’approche par division mais le second test n’a pas permis d’établir si une méthode de calcul est meilleure qu’une autre. Finalement, la méthode a été valorisée et transférée en implémentant les algorithmes développés dans un logiciel commercial. Ce logiciel permet de simuler la dispersion d’espèces réactives et non réactives (traceurs), dans un modèle contenant plusieurs centaines de compartiments organisés en tranches identiques / This PhD deals with modelling of chemical reactors with the “compartmental modelling” approach, which consists in dividing the system into a network from a dozen to several hundreds of interconnected volumes, called compartments. The structure of the network is deduced from tracer experiments, technical information about the chemical reactor and computational fluid dynamics flow simulations. This method provides a good compromise between computation time and results accuracy. When they are properly set-up, compartmental models give similar predictions, in terms of chemical reactions, as those of CFD simulations with a shorter calculation time and a more concrete representation of the reactor behavior. Every study from the literature is devoted to a specific reactor with a particular approach that cannot be straightforwardly transposed to other reactors. The aim of this PhD is to provide a contribution to the development of the most general possible methodology and to develop an automatic tool of generation and resolution of the differential equations system which must be solved. In the first chapter, a state of the art is proposed, defining the field of application of our method, in order to identify the most relevant division methods and the different methods to calculate the exchange between compartments. In the second chapter, a general methodology for compartmental modelling is developed. A versatile approach is proposed, consisting in dividing the reactor in identical slices. The calculation of exchange between compartments, both due to convection and turbulence, is presented in detail, with the description of three calculation methods for turbulent exchange. An interface has been developed, allowing to build any network of compartments. From this interface, the equations are written and solved automatically. The methodology is applied in the third chapter to an unfavorable case for slice cutting. This has allowed to test the limit of this approach. In particular, two points have been studied: (1) the applicability of division into identical slices and (2) the comparison between the turbulent exchange calculation methods. The first test has proved the robustness of the division approach but the second test has not allowed to establish whether one calculation method is better than another. Finally, the methodology has been promoted and transferred by implementing the developed algorithms within a commercial software. This software allows to simulate the dispersion of reactive and non-reactive (tracers) species, in model containing hundreds of compartments organized in identical slices

Ce travail porte sur l'étude du comportement rhéologique de fluide en milieu confiné. Pour cela le levier d'un microscope à force atomique (AFM) est utilisé pour sonder les propriétés rhéologiques d'un fluide confiné entre deux surfaces : la surface d'une sphère collée à l'extrémité du levier et une surface plane sur lequel le fluide est déposé. Le dispositif expérimental est constitué du système de mesure d'un AFM et d'un piézoélectrique permettant d'approcher ou d'éloigner de la sphère la surface plane. Un modèle analytique permet d'extraire les propriétés rhéologiques du fluide confiné à partir de la déflexion du levier induite par le pincement du fluide. Cette méthode a été validée pour les fluides newtoniens. Par contre pour les fluides non-newtoniens comme par exemple la solution de polyacrylami de nous avons trouvé que la viscosité dépend de la distance D et que le cisaillement n'est pas le seul paramètre pertinent pour interpréter les propriétés rhéologiques. / This thesis focuses on the study of the rheological behavior of confined fluids. For this purpose, the microcantilever of an atomic force microscope (AFM) is used to probe the rheological properties of a fluid confined between two surfaces, the surface of a sphere glued to the free-end of the AFM microcantilever and a flat solid surface on which the fluid is deposited. The set-up consists of an AFM, an electrical system for the deflection measurement and a piezoelectric device to move the solid surface (approach, oscillation, etc.). An analytical model allows to determine the rheological properties of the confined fluid from the measurement of the microcantilever deflection due to the hydrodynamic force exerted by the fluid on the sphere.This method has been validated for Newtonian fluids. For non-Newtonian fluids, such as polyacrylamide solution, we found that the viscosity depends on the distance D between the sphere and the plane surface and the shear rate is not the only relevant parameter for interpreting the rheological properties.

La vitrification des déchets radioactifs de haute activité et vie longue permet d'assurer leur confinement à long terme grâce à l'intégration à l'échelle atomique des radioisotopes au sein d'une matrice vitreuse. Parmi les différents procédés de vitrification, celui basé sur l'utilisation d'un creuset froid à induction électromagnétique directe présente de nombreux avantages. Il est à ce titre en fonctionnement depuis 2010 à l'usine française de retraitement de la Hague. Affiner la compréhension et la modélisation des phénomènes en jeu dans un creuset froid lors de l'élaboration de verre constitue un enjeu décisif pour le soutien industriel et l'appui des démarches d'optimisation du procédé. Dans cette optique, nos travaux proposent un modèle multiphysique ainsi qu'une simulation numérique prédictive prenant en compte l'ensemble des phénomènes couplés régissant le comportement de la charge de verre lors de son élaboration en creuset froid, en termes de convection, de flux thermiques et d'interactions électromagnétiques. Dans un premier temps sont formulées les hypothèses de travail permettant la modélisation mathématique du système magnéto-thermo-hydraulique avec un degré d'approximation contrôlé. Nous mettons notamment en place une approche multi-échelle pour la prise en compte de la nature diphasique du verre nucléaire, constitué d'une matrice homogène présentant des particules métalliques micrométriques de platinoïdes en suspension. La résolution numérique du modèle couplé est alors assurée par le développement et l'utilisation de codes de simulation précis, stables et optimisés, en proposant des schémas numériques et des méthodes de couplage adaptés. Une grande part de nos travaux est dédiée à l'analyse des phénomènes physiques caractéristiques de la fusion du verre, liés de manière générale aux propriétés physico-chimiques variables du matériau (viscosité, masse volumique, conductivité électrique...). Sont notamment étudiés deux grands axes qui constituent des apports applicatifs importants : le transport, la réponse à l'induction et la sédimentation des particules de platinoïdes ainsi que les échanges thermiques entre le verre en fusion et les parois refroidies du creuset. Des lois d'échelle pour l'évaluation de l'intensité respective de ces phénomènes sont dérivées de l'analyse du modèle couplé. Enfin, des comparaisons expérimentales précises, à la fois à petite échelle et à l'échelle industrielle, permettent de valider l'approche choisie et les résultats obtenus. / The vitrification of high-level nuclear waste ensures their long-term confinement through atomic integration of radioisotopes within a vitreous matrix. Among the various vitrification processes, the one based on the use of a cold crucible with direct electromagnetic induction shows many advantages. Thus it has been in operation since 2010 at the French reprocessing plant at La Hague. Refining the understanding and the modeling of the phenomena involved during glass melting inside a cold crucible constitutes a decisive stake for industrial support and process optimization. In this context, our work presents a multiphysics model as well as a predictive numerical simulation taking into account all the coupled phenomena governing the behavior of the glass melt, in terms of convection, heat fluxes and electromagnetic interactions. The working hypotheses are first formulated, allowing the mathematical modeling of the magneto-thermo-hydraulic system with a controlled degree of approximation. In particular, we design a multi-scale approach to take into account the diphasic nature of nuclear glass, consisting of a homogeneous matrix seeded with micrometric particles mostly made of platinum-group metals. The numerical resolution of the coupled model is then assured by the development and use of precise, stable and optimized simulation codes, by proposing numerical schemes and adapted coupling methods. A large part of our work is dedicated to the analysis of the physical phenomena characteristic of glass melting, generally related to the variable physicochemical properties of the material (viscosity, density, electrical conductivity...). In particular, two main axes are studied here with important applications: transport, response to the induction and settling of the platinum-group metal particles and heat transfers between the molten glass and the cooled walls of the crucible. Scaling laws for the evaluation of the respective intensity of these phenomena are derived from the coupled model analysis. Finally, precise experimental comparisons, both on a small scale and on an industrial scale, make it possible to validate the chosen approach and the outcoming results.

Le mélange de solutés par les écoulements en milieux poreux contrôle les réactions chimiques dans un grand nombre d'applications souterraines, dont le transport et la remédiation des contaminants, le stockage et l'extraction souterrains d'énergie, et la séquestration du CO2. Nous étudions les mécanismes du mélange à l'échelle du pore et plus précisément comment la topologie de l'écoulement est reliée à la dynamique du mélange d'espèces conservatives; en particulier, l'émergence d'un mélange chaotique est-elle possible dans un milieu poreux tridimensionnel (3D) ? Nous calculons donc numériquement ou mesurons expérimentalement les vitesses d'écoulement et l'évolution temporelle des champs de concentration afin de caractériser la déformation et le mélange à l'échelle du pore. Une première étude, expérimentale, permet de caractériser le mélange dans un fluide s’écoulant à travers un milieu poreux bidimensionnel (2D). Nous mesurons les vitesses par suivi de microparticules solides (''PTV''). L’évolution temporelle de la distance séparant deux particules permet de caractériser la dynamique de la déformation lagrangienne. Des mesures de transport conservatif dans le même milieu fournissent l'évolution temporelle du gradient de concentration moyen (une mesure du mélange). À partir de ces résultats expérimentaux nous proposons la première validation expérimentale à l'échelle du pore de la théorie lamellaire du mélange, reliant les propriétés de la déformation du fluide à la dynamique du mélange. Dans une deuxième étude nous examinons les conditions d'apparition du mélange chaotique dans l’écoulement dans des milieux poreux 3D granulaires ordonnés. Nous effectuons des calculs numériques hautement résolus de d'écoulement de Stokes entre des sphères empilées selon une structure cristalline (cubique simple ou cubique centrée), périodique. La déformation lagrangienne, obtenue à partir des champs de vitesse à l'aide d'outils numériques développés spécifiquement, met en lumière une large variété de dynamiques de la déformation dans ces milieux 3D, selon l'orientation de l'écoulement. Quand la direction de l'écoulement n'est pas normale à l'un des plans de symétrie de réflection du cristal, l'évolution temporelle de la déformation est exponentielle, traduisant une advection chaotique. L’émergence (ou non) du chaos est contrôlée par un mécanisme similaire à la ''transformation du boulanger'': les particules fluides se déplaçant autour d'un grain solide se retrouvent séparées par une surface virtuelle (appelée “variété”) qui émerge de la surface du grain. De multiples variétés existent dans l’écoulement, et la façon dont elles s'intersectent contrôle la nature - chaotique ou non - du mélange, et l'intensité du chaos. En particulier, l'exposant de Lyapunov (une mesure du chaos), est contrôlé par la fréquence spatiale des intersections appropriées à la génération du chaos, nommées ''connections hétéroclinines'' entre variétés. L'image conventionnelle, 2D, des mécanismes du mélange, impose des contraintes topologiques qui ne permettent pas le développement de ces mécanismes 3D. Elle pourrait donc être inadaptée aux milieux poreux naturels. La troisième étude a deux objectifs: (i) fournir une preuve expérimentale de la nature chaotique de l'advection, par la visualisation des variétés et par l'obtention d'une mesure de l'exposant de Lyapunov; et (ii), évaluer si nos résultats numériques obtenus pour des milieux granulaires ordonnés peuvent être généralisés à des milieux désordonnés, plus proches des milieux naturels. L’expérience est fondée sur un empilement désordonné de sphères rendu transparent par l'ajustement optique du liquide avec les sphères. La fluorescence induite par laser (''LIF'') permet de détecter les variétés au sein de l'écoulement, et des techniques PTV de mesurer les vitesses d'écoulement et quantifier l'exposant de Lyapunov. Les premiers résultats expérimentaux sont prometteurs. / Solute mixing in porous media flows plays a central role in driving chemical reactions in a number of subsurface applications, including contaminant transport and remediation, subsurface energy storage and extraction, and CO2 sequestration. We study the mechanisms of solute mixing, in particular how the pore scale flow topology is related to the mixing dynamics of conservative solutes, with a particular emphasis on the possible emergence of chaotic mixing processes in three-dimensional (3D) porous media. To do so, we perform numerical computations or experimental measurements of the flow velocities and temporal evolution of the concentration fields, and characterize fluid deformation and mixing at the pore scale. This PhD work consists of three main studies. In the first study, we experimentally characterize mixing in a fluid flowing through a two-dimensional (2D) porous medium built by lithography. We measure the velocity distributions from Particle Tracking Velocimetry (PTV). The time evolution of the separation distance between two particles is analyzed to characterize the Lagrangian deformation dynamics. In parallel we perform conservative transport experiments with the same porous media, and quantify the temporal evolution of the mean concentration gradient, which is a measure of the mixing rate. From these experimental results we obtain the first experimental pore scale validation of the lamella mixing theory, which relates the fluid deformation properties to the mixing dynamics. In the second study, we investigate the conditions of emergence of chaotic mixing in the flow through 3D ordered granular porous media. In these periodic cubic crystalline packings (Simple Cubic - SC - and Body-Centered Cubic - BCC) of spheres, we are able to perform highly resolved computations of the 3D Stokes flow. Using custom-developed numerical tools to measure the Lagrangian deformation from the computed velocity fields, we uncover the existence of a rich array of Lagrangian deformation dynamics in these 3D media, depending on the flow orientation. When the flow direction is not normal to one of the reflection symmetry planes of the crystalline lattice, we find that the Lagrangian deformation dynamics follow an exponential law, which indicates chaotic advection. This chaotic behavior is controlled by a mechanism akin to the baker's transformation: fluid particles traveling around a solid grain along different paths end up either separated by, or on the same side of, a virtual surface projecting from the grain surface and called a manifold. Multiple such manifolds exist within the flow, and the way they intersect controls the nature of mixing (that is, either non-chaotic or chaotic), and the strength of chaos. We show in particular that the magnitude of the Lyapunov exponent (a measure of the vigor of chaos) is controlled by the spatial frequency of transverse connections between the manifolds (called heteroclinic intersections). We thus demonstrate that the conventional 2D picture of the mechanisms of mixing may not be adapted for natural porous media because that picture imposes topological constraints which cannot account for these important 3D mechanisms. The third study has two objectives: (i) provide experimental evidence of the chaotic nature of pore scale advection/mixing, both by visualizing the manifolds and by obtaining a quantitative estimate of the Lyapunov exponent; and (ii) assess if the results obtained numerically in ordered packings of spheres extend to random packings, which are closer to natural porous media. The experiment features a random packing of glass beads rendered transparent by optical index-matching between the fluid and solid grains. We use Laser Induced Fluorescence (LIF) to detect the manifolds, and PTV techniques to measure flow velocities and subsequently quantify Lyapunov exponent. The first experimental results are promising.