Ce travail porte sur la dynamique de la matière granulaire dans l'état fluide, et sa réponse à une excitation localisée dans deux limites : une faible perturbation suite à laquelle le système relaxe rapidement vers un état homogène, ou une agitation intense donnant lieu à une onde de choc telle qu'un souffle d'explosion. Cette réponse est affectée par deux caractéristiques des fluides granulaires : les particules macroscopiques qui les composent sont d'une part inélastiques, de sorte que leur dynamique est dissipative et ne possède pas d'état d'équilibre, et d'autre part polydisperses, c'est-à-dire hétérogènes en taille et en masse. Nous isolons d'abord un effet dynamique de la polydispersité en montrant qu'il existe un mélange optimal qui minimise le temps de relaxation du fluide vers son état asymptotique. Nous nous intéressons ensuite au cas où une seule des espèces est perturbée par l'application d'un champ extérieur, et étudions l'état stationnaire hors d'équilibre ainsi établi, dans la limite du traceur où les autres espèces constituent un bain stationnaire. Enfin, nous modélisons la croissance de formes autosimilaires dans ce bain suite à une intense libération ponctuelle d'énergie, que nous comparons au souffle d'une explosion dans un gaz moléculaire. / This work focuses on the dynamics of the fluid state of granular matter, and its response to a localized perturbation in two limiting cases : relaxation toward a homogeneous state or growth of a blast wave. This response is shaped by two distinctive features of granular media: their macroscopic constituent particles are both inelastic, entailing intrinsically non-equilibrium dynamics, and polydisperse or heterogeneous in their material properties. First, we isolate the effects of polydispersity in the return of a gas to its homogeneous asymptotic state, and evidence the existence of an optimal mixture for which the relaxation time is minimal. We then restrict the perturbation to accelerating a single species with an external field in order to study the induced non-equilibrium stationary state in the tracer limit, where other species are undisturbed by this process. Finally, we model the self-similar shock wave generated in such a dissipative bath by an intense yet localized release of energy, and contrast it with blast waves in molecular gases.

Afin de permettre la mise en place d’une filière de valorisation du CO2 issu de l’industrie sidérurgique, cette étude a pour but de proposer un nouveau procédé de dépoussiérage des fumées de hauts fourneaux qui serait placé à la suite des filières actuelles de traitement pour améliorer leur épuration. L’originalité de la technologie retenue, un lit granulaire arrosé, permet de maintenir la perte de charge constante en facilitant le réentraînement en continu des particules collectées à l’aide d’un écoulement de liquide (formant un film à la surface des collecteurs) pour ainsi pallier le problème de colmatage. Une étude hydrodynamique et des mesures d’efficacité vis-à-vis de particules de diamètre compris entre 0,5 et 10 µm ont montré que la perte de charge et l’efficacité initiale augmentaient avec la diminution du diamètre des collecteurs, ce qui est en accord avec les résultats de la littérature pour des lits granulaires secs. La présence d’un filet d’eau augmente l’efficacité initiale de collecte au détriment toutefois d’une perte de charge plus élevée. Ces observations sont liées à la diminution de la porosité du lit lorsque ce dernier est arrosé par un débit d’eau. Un bon compromis entre perte de charge et efficacité a été trouvé avec un lit composé de 50 cm de billes de 5 mm de diamètre et un débit d’eau de 12 L.min-1 pour un débit de gaz à traiter de 20 m3.h-1. Des essais de colmatage longue durée ont mis en évidence que l’efficacité du lit reste globalement constante au cours du temps, à une valeur proche de la valeur initiale. La perte de charge, bien que supérieure en début d’expérience à celle d’un lit granulaire sec, se stabilise rapidement à une valeur bien plus faible que celle mesurée pour un lit sec pour une même masse de particules collectées. Les expériences ayant été menées sur 30 h dans des conditions défavorables, nous pouvons nous attendre à ce que dans des conditions réelles de fonctionnement la perte de charge soit stable sur des temps beaucoup plus longs. Ainsi la présence de l’eau remplit bien le rôle prévu en limitant l’évolution du colmatage du lit granulaire par le réentrainement des particules collectées par le lit. Cette étude s’est achevée par le développement d’un modèle permettant de déterminer la perte de charge et l’efficacité initiale d’un lit granulaire arrosé, ce qui a abouti sur un exemple de pré-dimensionnement à échelle industrielle / In order to enable the establishment of a CO2 recovery chain sector from the steel industry, this study aims to propose a new filtration device for the dust removal of blast furnace fumes that would be placed after the current treatment chain to improve their purification. The originality of the chosen technology, a trickle bed, permits to maintain a constant pressure drop by re-entraining continuously the particles trapped on collectors thanks to a water film flow to overcome the clogging problem. A hydrodynamic study and efficiency measurements lead with particles (in a range of diameters between 0.5 and 10 μm) showed that the pressure drop and the initial efficiency increase with the decrease of the collector diameter, which is in agreement with the results of the literature for dry granular beds. The presence of water increases the initial collection efficiency at the expense of a higher pressure drop. These observations are related to the decrease of the bed porosity in presence of a water flow. A good compromise between pressure drop and efficiency was found for an air flow rate of 20 m3.h-1, a liquid flow rate of 12 L.min-1 and using a 50 cm bed composed of glass beads of 5 mm diameter. Long-term clogging tests showed that the efficiency of the bed remains globally constant over time at a value close to the initial one. The pressure drop, although higher at the beginning of the experiment than that of a dry granular bed, stabilizes rapidly at a much lower value than that measured for a dry granular bed for the same mass of particles collected. The experiments have been carried out over 30 h under disadvantageous conditions, and we can expect that under real operating conditions the pressure drop will be stable over much longer times. Thus, the presence of water plays its part by limiting the clogging evolution of the granular bed by the re-entrainment of particles collected by the bed. This study was completed by the development of a model to determine the pressure drop and initial efficiency of a trickle bed, which resulted in an example of pre-sizing of a process on an industrial scale

Les glissements de terrain se produisent dans toutes les chaînes de montagnes où la résistance de massifs rocheux est insuffisante pour contrer l’action de la gravité. Les terrains volcaniques sont particulièrement susceptibles de s’effondrer car les édifices sont composés des lithologies diverses et variées qui peuvent être fortement fracturées. En plus, la croissance rapide des édifices volcaniques favorise leur instabilité et leur effondrement. L’activité magmatique est un facteur additionnel responsable de la déformation des édifices, tandis que l’activité hydrothermale réduit la résistance des roches volcaniques. Pour ces raisons, l’évaluation des aléas liés à l’effondrement des édifices et à la formation des avalanches des débris volcaniques mérite une attention particulière. Les caractéristiques physiques des composants des avalanches des débris ont une influence directe sur la dynamique de ce type d'écoulement. Les dépôts des avalanches de débris présentent une morphologie de surface composée des nombreuses collines (hummocks), qui montrent fréquemment les séquences volcaniques initiales, ce qui suggère un mécanisme de mise en place proche de celui des glissements de terrain. Cependant, d’autres dépôts présentent des crêtes allongées (rides) dont le mécanisme de formation est encore méconnu. Le volcan Tutucapa (sud du Pérou) a été affecté récemment par deux avalanches de débris. La plus ancienne, « Azufre », est d’âge Holocène et résulte de l’effondrement d’un complexe des dômes et d’une séquence volcanique altérée (hydrothermalisée) sous-jacente. La deuxième avalanche, « Paipatja », a eu lieu il y a seulement 200-230 ans BP et est associée à une grande éruption explosive du Tutupaca. Les dépôts de cette avalanche présentent notamment de nombreuses rides. Les deux dépôts d’avalanche montrent deux unités différentes : une unité inférieure, caractérisée par la présence des blocs altérés (hydrothermalisés) provenant de l’édifice basal, tandis que l’unité supérieure est constituée par des blocs du complexe de dômes actifs. Le travail de terrain montre que les rides de l’avalanche « Paipatja » présentent une forte variation de granulométrie entre leur partie centrale (enrichie en blocs grossiers) et leurs parties latérales, ce qui suggère un processus de ségrégation granulaire. Des expériences analogiques montrent que des écoulements de mélanges de particules des différentes tailles subissent un processus de ségrégation et de digitation granulaire qui engendre des rides par jonction de levées statiques qui délimitent un chenal d’écoulement. Le processus de formation des rides est facilité par de faibles différence de taille des particules dans des mélanges bidisperses. Ces résultats suggèrent que les rides observées au Tutupaca résultent d’un écoulement granulaire. Les principales caractéristiques morphologiques des structures formées lors de ces expériences de laboratoire ont été comparées qualitativement avec les structures observées dans les dépôts du Tutupaca. Les structures observées au Tutupaca montrent que deux mécanismes de mise en place peuvent coexister dans les avalanches de débris volcaniques : le glissement de blocs plus ou moins cohérents, et l’écoulement semblable à celui d’un matériau granulaire. Cela dépend probablement de la nature des différents matériaux à la source des avalanches. Cette information doit être prise en compte pour l’évaluation des aléas liés aux avalanches des débris car des mécanismes d’écoulement différents peuvent induire des fortes variations de la distance parcourue par ces avalanches. / Landslides occur in all mountainous terrain, where the rock strength is unable to support topographic loading. Volcanic rocks are particularly landslide prone, as they mix strong and weak lithologies and are highly pre-fractured. Also, volcanoes themselves, are peculiar mountains, as they grow, thus creating their own topographic instability. Magmatic activity also deforms the edifice, and hydrothermal activity reduces strength. For all these reasons, volcanoes need close consideration for hazards, especially for the landslide-derived rock avalanches. The characteristics and properties of different debris avalanche components influence their behavior during motion. Deposits are generally hummocky, preserving original layering, which indicates a slide-type emplacement. However, some deposits have ridged morphology for which the formation mechanisms are not well understood. Two recent debris avalanches occurred at the Tutupaca volcano (S Peru). The first one, “Azufre” is Holocene and involved the collapse of active domes and underlying older hydrothermally altered rocks. The second debris avalanche, “Paipatja” occurred 200-230 y BP and is associated with a large explosive event and this deposit is ridged. The excellent conservation state of the deposits and surface structures allows a comprehensive analysis of the ridges. Both deposits have two contrasting units: a lower basal edifice-derived hydrothermally-rich subunit and an upper dome-derived block-rich unit. Detailed fieldwork has shown that Paipatja ridges have coarser core material and are finer in troughs, suggesting grain size segregation. Using analog experiments, the process that allow ridge formation are explored. We find that the mixtures undergo granular segregation and differential flow that create fingering that forms ridges by junction of static léeves defining a channel flow. Granular segregation and fingering are favored by small particle size contrast during bi-dispersed flow. The results suggest that the ridges observed at Tutupaca are product of a granular flow We extract the morphological characteristics of the deposits of granular flows generated in the laboratory and make a qualitative comparison with the Tutupaca deposits. The description of the different landslide and debris avalanche features at Tutupaca shows that two types of debris avalanche motion can occur in volcanic debris avalanches: the sliding of blocks more or less coherent and a flow similar to a granular material. This probably depends on source materials and the conditions of different parts of the initial landslide. Such information should be taken into account when estimating hazards at other volcanic landslide sites, as the different behaviors may result in different run outs.

Les matériaux granulaires présentent une large gamme de lois de comportement lorsqu'ils sont soumis à différents chemins de chargement. Le développement de modèles constitutifs permettant de rendre compte de ces caractéristiques a été une préoccupation constante de nombreux chercheurs depuis des décennies. Parmi les différentes options possibles, les approches par changement d’échelle semblent prometteuses. Dans ces approches, le modèle constitutif est formulé en reliant les propriétés macroscopiques du matériau aux propriétés micro-structurelles correspondantes.Cette thèse propose un modèle micromécanique tridimensionnel (le modèle H-3D) prenant en compte une échelle intermédiaire (méso-échelle). Il permet ainsi de décrire de manière naturelle un grand nombre de caractéristiques constitutives des matériaux granulaires non cohésifs. La comparaison entre essais expérimentaux et simulations numériques révèle la capacité prédictive de ce modèle. En particulier, des simulations réalisées avec différentes pressions de confinement et différents rapports de vide initiaux ont permis de démontrer la capacité du modèle à rendre compte quantitativement de l'état critique sans nécessiter d’équation spécifique et de paramètre d'état critique. Le modèle est également analysé à l’échelle microscopique, où l'évolution de certains paramètres microscopiques clés est présentée.Une approche multi-échelle 3D est ensuite présentée afin d’étudier le comportement mécanique d'un échantillon macroscopique constitué d'un assemblage granulaire, en tant que problème aux conditions limites. Le cœur de cette approche est un couplage multi-échelle, où la méthode des éléments finis est utilisée pour résoudre le problème aux conditions limites et le modèle H-3D est utilisé pour calculer la loi de comportement à l’échelle d’un volume élémentaire représentatif. Cette approche fournit un moyen pratique de relier les observations macroscopiques avec les mécanismes microscopiques intrinsèques. Des conditions de chargement biaxiaux en déformations planes sont appliquées pour simuler le phénomène de localisation des déformations. Une série de tests est effectuée, où différents motifs de rupture sont observés et analysés. Un système de bande de cisaillement apparaît naturellement dans un spécimen initialement homogène. En définissant la zone de la bande de cisaillement, les mécanismes microstructuraux sont étudiés séparément à l'intérieur et à l'extérieur de celle-ci. En outre, une analyse directionnelle de travail du second ordre est effectuée en appliquant des petits incréments de contrainte à différents états de contrainte-déformation sur des chemins de chargement biaxiaux drainés. Le travail de second ordre normalisé, introduit comme un indicateur d’instabilité du système, est analysé non seulement à l’échelle macroscopique mais aussi à l’échelle microscopique.Enfin, une analyse du travail de second ordre appliquée à des problèmes géotechniques et utilisant l'approche multi-échelle développée dans cette thèse est présentée. L'approche multi-échelle est utilisée afin de simuler des problèmes aux conditions limites homogènes et non homogènes, offrant ainsi la possibilité d’interpréter et de comprendre les micro-mécanismes qui à l’origine des phénomènes de rupture dans les problèmes géotechniques. Cette approche multi-échelle utilise un schéma numérique d’intégration dynamique-explicite afin de pouvoir étudier la rupture post-pic sans avoir à recourir à des outils mathématiques trop sophistiqués. Ainsi, en changeant le type de condition de chargement de déplacement à contrainte lorsque le système atteint son état limite, son effondrement se traduit par une augmentation soudaine de l'énergie cinétique découlant de la différence entre les travaux internes et externes du second ordre. / Granular materials exhibit a wide spectrum of constitutive features when submitted under various loading paths. Developing constitutive models which succeed in accounting for these features has been challenged by scientists for decades. A promising direction for achieving this can be the multi-scale approach. Through this approach, the constitutive model is formulated by relating material’s macroscopic properties to their corresponding microstructure properties.This thesis proposes a three-dimensional micro-mechanical model (the so-called 3D-H model) taking into account an intermediate scale (meso-scale) which makes it possible to describe a variety of constitutive features in a natural way. The comparison between experimental tests and numerical simulations reveals the predictive capability of this model. Particularly, several simulations are carried out with different confining pressures and initial void ratios, based on the fact that the critical state is quantitatively described without requiring any critical state formulations and parameter. The model is also analyzed from a microscopic view, wherein the evolution of some key microscopic parameters is investigated.Then, a 3D multi-scale approach is presented to investigate the mechanical behavior of a macroscopic specimen consisting of a granular assembly, as a boundary value problem. The core of this approach is a multiscale coupling, wherein the finite element method is used to solve a boundary value problem and the 3D-H model is employed to build the micro constitutive relationship used at a representative volume element scale. This approach provides a convenient way to link the macroscopic observations with intrinsic microscopic mechanisms. Plane-strain biaxial loading conditions are selected to simulate the occurrence of strain localization. A series of tests are performed, wherein distinct failure patterns are observed and analyzed. A system of shear band naturally appears in a homogeneous setting specimen. By defining the shear band area, microstructural mechanisms are separately investigated inside and outside the shear band. Moreover, a second-order work directional analysis is performed by applying strain probes at different stress-strain states along drained biaxial loading paths. The normalized second order work introduced as an indicator of an unstable trend of the system is analyzed not only on the macroscale but also on the microscale.Finally, a second order work analysis in application to geotechnical problems by using the aforementioned multiscale approach is presented. The multiscale approach is used to simulate a homogeneous and a non-homogeneous BVP, opening a road to interpret and understand the micro mechanisms hiding behind the occurrence of failure in geotechnical issues. This multiscale approach utilizes an explicit-dynamic integral method so that the post-peak failure can be investigated without requiring over-sophisticated mathematical ingredients. Thus, by switching the loading method from a strain control to a stress control at the limit state, the collapse of the system can be reflected in an abrupt increase of kinetic energy, stemming from the difference between both internal and external second-order works.

L'analyse de la pression exercée par un écoulement granulaire sur un obstacle nécessite la connaissance du coefficient de traînée. Néanmoins, les études bibliographiques expérimentales ne fournissent pas les données internes de l'écoulement intervenant dans son calcul. C'est pourquoi, nous proposons la détermination du coefficient de traînée par l'application d'outils numériques tridimensionnels, basés sur la méthode des éléments discrets (MED), et validés par des études expérimentales. L'expérience est un canal d'écoulement de laboratoire, intégralement modélisé par la MED. Les lois locales de contact intègrent un comportement normal élastique, hystérétique associé à un critère de glissement tangentiel. La validation du modèle intégral repose sur la similarité entre les résultats expérimentaux et numériques des paramètres physiques de l'écoulement et de l'effort d'impact sur l'obstacle. Cette validation est obtenue sans calibration importante des paramètres numériques. Le modèle intégral est ensuite optimisé, en termes de temps de calcul et de possibilités d'études, sous la forme d'un canon granulaire, dans lequel n'est simulée que la zone d'interaction entre l'écoulement et l'obstacle. Les propriétés de l'écoulement : vitesse, épaisseur, densité, sont contrôlées car intégrées en tant que paramètres initiaux numériques. L'influence de caractéristiques associées à l'écoulement, comme le nombre de Froude, à la taille des grains, ou à la forme et la taille de l'obstacle, sur le coefficient de traînée, est alors analysée. Enfin, une extension micromécanique est proposée dans l'analyse de l'influence d'une adhésion et d'une cohésion intégrées aux lois locales de contact.

La transition métal-isolant (TMI) induite par le dopage et la supraconductivité ont été étudiées dans le diamant monocristallin fortement dopé au bore (propriétés structurales, approche en loi d'échelle pour la TMI, substitution isotopique, structure de bande vs effet Hall, mode de vibration des dimères de bore) et des matériaux voisins. La supraconductivité a été découverte pour le silicium monocristallin fortement dopé au bore (préparé par la méthode GILD), la TMI dans le cas du carbure de silicium 4H fortement dopé à l'aluminium (préparé par la méthode VLS). L'étude de la TMI et du magnétotransport dans le diamant nanocristallin fortement dopé au bore a montré l'importance de la granularité de ce système.

Nous étudions les phénomènes de ségrégation par taille dans des suspensions concentrées<br />de particules en écoulement. Les particules sont plus lourdes que le fluide, comme dans la<br />plupart des situations pratiques (boues, bétons, propergols solides. . .). Le matériau modèle<br />est transparent, et nous pouvons mesurer les champs de vitesse et de concentration des<br />espèces sous cisaillement.<br /><br />Dans les deux géométries étudiées (Couette et plan-plan), on observe que l'écoulement<br />principal (azimutal) est en général localisé. De plus, nous montrons l'existence systématique<br />d'une composante non azimutale des trajectoires.<br /><br />Expérimentalement, nous établissons les relations entre figure de ségrégation, écoulements<br />primaire et secondaire. Sur cette base, nous discutons les mécanismes possibles de<br />la ségrégation.

La propagation d'ondes sonores dans des empilements granulaires présente des comportements originaux et mal-élucidés liés à l'aspect divisé de ces matériaux. Ceux-ci se comportent notamment comme des milieux désordonnés et non-linéaires. Notre étude est originale car nous employons des grains photoélastiques. Couplés à une caméra rapide, nous visualisons le champ de déformation associé à la propagation de l'onde sonore. Nous abordons d'abord le cas de l'empilement le plus simple : une chaîne unidimensionnelle de grains cylindriques en contact. Nous étudions la propagation d'impulsions sonores dans les deux cas où leurs amplitudes sont petites (régime linéaire) ou grandes (régime non-linéaire) devant la force de confinement statique appliquée à la chaîne. Dans le premier cas, nous observons l'effet des imperfections de surface des grains sur la vitesse des ondes. Nous constatons aussi que la dissipation est essentiellement due au frottement solide. Pour les grandes amplitudes, l'impulsion initiale se décompose en un train de pics d'amplitudes décroissantes, de largeurs comprises entre 3 et 4 grains et de vitesses supersoniques. La vitesse du pic principal, adimensionnée par la vitesse du son linéaire, ne dépend que du rapport entre l'amplitude de l'onde et la force statique. Ces observations sont interprétées en généralisant les résultats de Nesterenko. Nous présentons ensuite le cas des empilements bidimensionnels. Les ondes linéaires se propagent dans les chaînes de forces et leurs vitesses croient avec la force statique. Elles restent cependant toujours inférieures à celles mesurées à 1D. Cela suggère que la dynamique de l'onde le long d'une chaîne de force est influencée par les contacts latéraux à cette chaîne de force.

L'idée que nous avons développé dans ce travail est de faire usage d'une commande électrique unique pour modifier à la fois les propriétés diélectrique et magnétique de nanomatériaux contenus dans le dispositif hyperfréquence. Pour cela, nous cherchons à exploiter le concept de propriété produit, dont un exemple est donné par l'effet magnétoélectrique (ME). Si l'on considère un mélange de deux phases condensées en contact, la première étant magnétostrictive et la seconde piézoélectrique, alors l'application d'un champ magnétique externe sur la phase magnétostrictive génère un champ de contrainte mécanique locale activant la piézoélectricité et ainsi fait apparaître un gradient de potentiel. En d'autres termes, la magnétoélectricité conduit à une variation de la permittivité par l'application d'un champ magnétique externe, et par commutativité, à une dépendance de la perméabilité magnétique sous l'action d'un gradient de potentiel. Les nanostructures granulaires présentent de nombreux avantages par rapport aux céramiques, aux couches minces et aux matériaux monophasés. La faible conductivité, l'absence d'interdiffusion et de réactions chimiques parasites dues au frittage, l'absence de substrat, la conservation de l'échelle nanométrique, la facilité et le faible coût du protocole d'élaboration de ces NCs en font des matériaux de choix pour la réalisation d'un effet ME en vue de proposer des dispositifs hyperfréquences originaux. Le fait saillant de ce travail a été de mettre en évidence un couplage magnétoélectrique (ME) dans des nanocomposites (NCs) granulaires biphasiques dans le domaine micro-onde à température ambiante. L'usage de la résonance ferromagnétique (RFM-LMB) sur une large gamme de fréquences (6-28 GHz) et de spectroscopie micro-onde (100 MHz-6 GHz) combiné à une étude du transport électronique et à l'utilisation de plusieurs sondes morphologiques nous a permis de corréler les propriétés dynamiques de l'aimantation, les paramètres électromagnétiques des matériaux, et les informations microstructurales au coefficient ME.

La réponse à des vibrations et l'étude de la stabilité des pentes ont une place particulière dans l'étude des milieux granulaires. Ces études peuvent d'une part nous renseigner sur le comportement des milieux granulaires présents dans le milieu naturel et souvent soumis à ces deux types de sollicitations (avalanches de grains, éboulements de roches, dynamique des dunes, tremblements de terre, ...). D'autre part, elles permettent d'appréhender ces milieux comme des systèmes modèles hors équilibre et d'analyser les résultats à la lumière de ceux obtenus sur la dynamique des systèmes vitreux. Nous avons tout d'abord étudié expérimentalement par des méthodes optiques et acoustiques la réponse d'un milieu granulaire lentement incliné. La dynamique présente deux régimes : un régime de petits réarrangements indépendants et localisés est suivi d'un régime intermittent caractérisé par de grands événements de type "précurseurs" au voisinage de l'angle d'avalanche. Ce régime dit de "précurseurs" correspond à un mouvement collectif de l'ensemble des grains situés à la surface de l'empilement. Nous mettons en évidence l'importance du détail microscopique de l'interaction entre les grains et de la condition à la limite à la base de l'empilement sur la dynamique du système dans ce régime. La dynamique globale du système a été étudiée et comparée à un modèle simple de type automate cellulaire. Dans un second temps, nous avons étudié numériquement la réponse d'un milieu granulaire soumis à des vibrations. Le comportement du système présente des caractéristiques qui sont communes à tous les systèmes hors équilibre (distributions non-gaussiennes des déplacements, ralentissement et hétérogénéités de la dynamique,...). L'influence de la condition à la limite au fond de l'empilement sur la structuration de celui-ci et sur le ralentissement de la dynamique a également été discutée. Nous avons, en particulier, mis en évidence une zone d'épaisseur environ égale à 5 diamètres de grains où des hétérogénéités dynamiques et de fortes corrélations entre les mouvements des grains apparaissent.