L'érosion de contact est un type d'érosion interne qui se développe à l'interface entre deux couches de matériaux de granulométries différentes. Les particules d'un sol fin (sable, limon, argile…) sont détachées par l'écoulement et entraînées à travers les pores du sol grossier au contact (gravier…). Bien que l'on suspecte sa présence dans de nombreux ouvrages, ce processus d'érosion a été peu étudié jusqu'à présent. Aussi ce travail de thèse s'est attaché à mieux le comprendre pour parvenir à le modéliser. A l'échelle du pore, l'écoulement à l'interface entre deux milieux poreux a été caractérisé, grâce à un dispositif expérimental développé au Cemagref d'Aix-en-Provence qui combine fluorescence induite par laser, méthode PIV et milieu iso-indice. L'importance de la variabilité des sollicitations hydrauliques a ainsi été soulignée. A l'échelle de l'échantillon, des essais sur sols réels et sols reconstitués ont été menés au LTHE, afin d'identifier les phénomènes en jeu. Ont également été obtenues sur ce dispositif des mesures du taux d'érosion en fonction de l'intensité de l'écoulement, pour différents types de sols fins et de sols grossiers. Une modélisation stochastique de ces essais a ensuite été proposée. Enfin, des essais à grande échelle ont été conduits au laboratoire de la Compagnie Nationale du Rhône pour étudier l'éventuelle influence d'effets d'échelle ainsi que les conséquences de cette érosion de contact sur le comportement global et l'intégrité d'un ouvrage. Il a ainsi été mis en évidence la possibilité qu'une érosion de conduit soit initiée par érosion de contact. Cette thèse a été réalisée dans le cadre du projet ERINOH (ERosion INterne dans les Ouvrages Hydrauliques), en convention CIFRE avec le Centre d'Ingénierie Hydraulique d'EDF. / Contact erosion is a type of internal erosion which occurs at the interface between two soil layers of different particle sizes. Particles from the fine soil (sand, silt, clay…) are pulled and driven by the flow through the pores of the surrounding coarse soil (gravel ...). Although its presence is suspected in many earthfill embankments, this process has been little studied so far. The aim of this thesis was to better understand and, possibly, to model contact erosion. At pore scale, the flow at the interface between two porous media was characterized thanks to an experimental device, developed in Cemagref, combining Planar Laser Induced Fluorescence, PIV method and Refractive Index Matching. The spatial variability of flow shear stress has been emphasized. At sample scale, several tests on real and reconstituted soils were carried out in LTHE, in order to identify the phenomena involved in contact erosion. Erosion rates as a function of flow magnitude have also been measured for different types of fine and coarse soils. Then, a stochastic model was proposed to account for these experiments. Finally, large scale tests were conducted in the laboratory of Compagnie National du Rhône, to study the influence of scale effects as well as the consequences of contact erosion on global behavior and integrity of the embankment structure. The possibility of piping initiation by contact erosion was underlined. This thesis is part of the project ERINOH project and was funded by EDF-CIH.

L’érosion interne est la cause principale de ruptures des ouvrages hydrauliques en terre, tels que les digues et les barrages. Les conséquences de telles ruptures sont conséquentes et couteuses. L’objectif de cette thèse est de mieux comprendre une catégorie particulière de processus d’érosion, la suffusion, en en caractérisant l’initiation et l’évolution. Cette recherche, à dominante expérimentale, s'appuie sur une modélisation physique unidimensionnelle. Les essais sont réalisés dans une conduite en plexiglas de 18 cm de diamètre intérieur, équipée de capteurs de pression, d’un débitmètre et d’un turbidimètre. Un dispositif de collecte des particules érodées est installé à l’aval de la conduite pour quantifier les particules érodées à des intervalles de temps choisis. Dans cette étude, on a réalisé une série d’essais sur des matériaux pulvérulents. Dans ce type de matériaux, l’érosion se manifeste par suffusion. Elle correspond à l’arrachement et au transport des particules fines à travers l’espace poral des particules grossières. Dans un premier temps, un protocole de reconstitution de matériaux pulvérulents est élaboré. Le principe consiste à réaliser des mélanges humides avec une teneur en eau relative à la masse des particules fines. Pour surmonter les problèmes rencontrés dans la réalisation des essais, la conduite est disposée verticalement, et les matériaux testés sont reconstitués à partir de mélanges de particules sableuses. Une étude paramétrique est menée sur l’influence de certains paramètres, tels que la granulométrie, la pression d’écoulement, la nature et la teneur en particules argileuses ainsi que la densité de mise en place. Les résultats obtenus montrent que l’érodabilité des matériaux dépend de la forme des courbes granulométrique. En effet, l’érosion est plus accentuée pour les matériaux contenant moins de particules fines. Elle est également plus forte quand la courbe granulométrique est discontinue. L’ajout de particules argileuses augmente la résistance à la suffusion des matériaux testés. Trois catégories de matériaux argileux ou fins sont testées : une illite (argile verte de Velay), une kaolinite (Speswhite), et du sable broyé (C10). Les résultats montrent que les matériaux contenant de l’illite sont plus résistants à la suffusion que les matériaux contenant de la kaolinite, tandis que les matériaux contenant du sable broyé sont les moins résistants. Cette résistance à la suffusion augmente avec la teneur en particules argileuses. L’érosion des matériaux dépend aussi de la densité initiale des matériaux testés. En effet, pour la même pression d’entrée, la quantité des particules érodées diminue en augmentant la densité de mise en place / Internal erosion is the main cause of failure of hydraulic structures such dykes and dams. The consequences of such failures are substantial and costly. The objective of this thesis is to better understand one of the erosion phenomena, the suffusion, and to characterize the initiation and the evolution of this phenomenon. A physical modeling approach is used for this study. Tests are carried out in a Plexiglas pipe of 18 cm inner diameter, equipped with pressure sensors, flowmeter and turbidimeter. A device for collecting the eroded particles, at selected time intervals, is installed in the downstream part of the device. In this study we performed a series of tests on coarse cohesionless soil. In this type of material, erosion is manifested by suffusion. It corresponds to the detachment and transport of fine particles through the pore space of the coarse particles. First, a cohesionless soil reconstitution protocol is developed. The idea is to make wet mixtures with water content depending on the fine particle content. To overcome some problems encountered in carrying out the tests, the device is arranged vertically, and the materials tested are reconstituted from mixtures of sand particles. A parametric study on the influence of several parameters such as particle size distribution, flow pressure, nature and content of clay particles and initial density are conducted. The results show that the erodibility of the material depends on the shape of the particle size distribution. Indeed erosion is more pronounced for materials containing less fine particles. This erodability is more pronounced when the curve is gap-graded. Erodibility of the tested soils increases with the hydraulic load. The addition of clay particles increases the resistance to suffusion of the soils. Three categories of clayed or fine soils are tested: illite (Argile Verte de Velay), kaolinite (Speswhite), and crushed sand (C10). The results show that materials containing illite are more resistant to suffusion than materials containing kaolinite, whereas materials containing crushed sand are less resistant. This resistance to suffusion increases with clay particle content. The erosion of materials also depends on the initial density of the soil tested. Indeed, for the same applied pressure, the quantity of eroded particles decreases with increasing the initial density

Dams that comprise broadly and widely graded glacial materials, such as moraines, have been found to be susceptible to internal erosion, perhaps more than dams of other soil types. Internal erosion washes out fine-grained particles from the filling material; the erosion occurs within the material itself or at an interface to another dam zone, depending on the mode of initiation. Whether or not internal erosion proceeds depend on the adequacy of the filter material. If internal erosion is allowed, it may manifest itself as sinkholes on the crest, increased leakage and muddy seepage once it surfaces, which here is called surfacing internal erosion (i.e. internal erosion in the excessive erosion or continuation phase). In spite of significant developments since the 1980s in the field of internal erosion assessment, the validity of methods developed by others on broadly graded materials are still less clear because most available criteria are based on tests of narrowly graded granular soils. This thesis specifically addresses dams that are composed of broadly graded glacial soils and investigates typical indicators, signs and behaviors of internal erosion prone dams. Based on a review of 90+ existing moraine core dams, which are located mainly in Scandinavia as well as North America and Australia/New Zealand, this thesis will show that not only the filter’s coarseness needs to be reviewed when assessing the potential for internal erosion to surface (i.e., erosion in the excessive or continuing phase); in addition, the grading stability of the filter and the core material as well as non-homogeneities that are caused by filter segregation need to be studied. Cross-referencing between these aspects improves the assessment of potential for internal erosion in dams of broadly graded soils and furthermore it provides aid-to-judgment. / QC 20100715

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第20684号 / 工博第4381号 / 新制||工||1681(附属図書館) / 京都大学大学院工学研究科都市社会工学専攻 / (主査)教授 澤田 純男, 教授 清野 純史, 准教授 後藤 浩之 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

[pt] Em uma barragem de terra homogênea responsável pelo abastecimento de água de um empreendimento privado, observou-se um processo inicial de carreamento de partículas finas próximo ao pé do talude de jusante, situação tida como um dos principais modos de falha de barragens. Embora esta anomalia tenha sido devidamente controlada na estrutura em estudo, este trabalho buscou encontrar as causas para a ocorrência detectada, visando aprofundar o conhecimento sobre o assunto. Para isso, investigou-se o histórico da estrutura, coletou-se amostras indeformadas próximo a região de interesse com a função de realizar ensaios de laboratório para caracterizar o solo no local, determinar os parâmetros físicos, aferir coeficientes de permeabilidade sob diferentes tensões efetivas e encontrar a sua curva característica (SWCC) através de ensaios de papelfiltro. Com os resultados dos ensaios, foram elaborados modelos computacionais sob diferentes cenários para investigar o comportamento dos gradientes hidráulicos. Foi empregado o software Plaxis LE para esta função. Avaliou-se se a adoção de parâmetros não saturados ou a consideração de análises tridimensionais promoveria alguma alteração no padrão de fluxo na estrutura em estudo, o que não ocorreu. Também foi considerado se a interação solo-atmosfera poderia gerar gradientes hidráulicos críticos através de ciclos de chuva e evaporação, mas também não houve resultados críticos. Ao variar o coeficiente de permeabilidade no corpo da barragem, observou-se que existe uma relação inversamente proporcional com o gradiente hidráulico. Desta forma, entende-se que a presença de solos orgânicos no pé do talude de jusante pode ter acarretado o processo de carreamento observado. / [en] In a homogeneous earthfill dam responsible for the water supply of a private company, an initial process of internal erosion was observed near the base of the downstream slope, which is listed as one of the main dam failures causes in the world. Although this problem has been correctly solved at the studied structure, this thesis aimed to find the causes for the detected anomaly in order to increase knowledge about this subject. To do so, the history of the structure was investigated, undisturbed samples were collected near the region of interest to carry out laboratory tests to characterize the soil, determine the physical parameters, measure permeability coefficients under different effective stresses and find its characteristic curve (SWCC) through filter paper tests. With the results of the tests, computational models were elaborated under different scenarios to investigate the behavior of hydraulic gradients. The software Plaxis LE was used in this work. It was evaluated whether the adoption of unsaturated parameters or 3D analysis would promote any change in the flow pattern in the studied structure, which did not occur. It was also considered whether the soil-atmosphere interaction could generate critical hydraulic gradients through rainfall and evaporation cycles, but there were no critical results either. Although, permeability seems to have a an inversely proportional relationship with hydraulic gradient. Thus, it is understood that the formation of organic soils at the toe of the downstream slope may have led to the observed internal erosion process.

La plupart des digues sont constituées de matériaux granulaires compactés. Elles sont ainsi perméables et constamment soumises à des écoulements d’eau dans leur volume. Dans certaines conditions, ces écoulements peuvent altérer leur microstructure par érosion interne et générer des instabilités mécaniques responsables de ruptures inopinées lors de crues. Cette thèse s’intéresse à l’analyse multi-échelle des instabilités mécaniques dans les matériaux granulaires soumis à l'érosion interne. Dans ce travail, le comportement mécanique de ces matériaux est simulé en 3D à l’échelle de volumes élémentaires représentatifs, et ce, pour différents états de contraintes et gradients hydrauliques. Grâce à l’utilisation du critère du travail du second ordre et d’outils micromécaniques, leur stabilité est analysée avant et après l’application d’un écoulement interne. Il est établi que l’origine micro-inertielle des instabilités observées provient du déconfinement et de la flexion des chaînes de force ainsi que des déformations plastiques importantes résultant de leur effondrement. Par leur capacité à enrayer rapidement le développement de telles déformations plastiques, il est montré que les particules libres contribuent à assurer la stabilité mécanique des matériaux granulaires. Ce résultat est fondamental pour analyser les conséquences de l’érosion interne en termes de stabilité mécanique car les particules libres sont facilement transportables sous l’action d’un écoulement interne. Selon si elles sont colmatées ou érodées, un écoulement interne aura un effet stabilisateur ou déstabilisateur vis-à-vis du comportement mécanique des matériaux granulaires soumis à l’érosion interne / Dikes are most of the time built of compacted granular materials that are permeable and continuously subjected to internal fluid flows. In some cases, microstructure modifications resulting from internal erosion generate mechanical instability that will lead to unexpected failures in case of serious flooding. This thesis focuses on multi-scale analysis of mechanical instability in granular materials subjected to internal erosion. In this work, the mechanical behavior of such materials is simulated in three dimensions at the scale of representative elementary volumes subjected to different stress states and hydraulic gradients. Thanks to the use of the second order work criterion and micromechanical tools, the mechanical stability of these materials is tested before and after internal erosion. It is established that the micro-inertial origin of the observed instabilities is linked to force chain deconfinement and bending as well as to the development of large plastic strains resulting from force chain collapse. By preventing the development of such plastic strains, it is shown that rattlers contribute to ensure the mechanical stability of granular materials. This key finding is of a particular significance in relation with internal erosion as rattlers can be easily transported under the action of an internal fluid flow. Depending on whether they get clogged or eroded, an internal fluid flow has thus either a stabilizing or a destabilizing effect on the mechanical behavior of granular materials subjected to internal erosion

L'érosion interne est le résultat de l'arrachement des particules d'un sol sous l'action d'un écoulement hydraulique. Ce phénomène est à l'origine des ruptures hydrauliques des barrages et des digues. Sur le plan mondial, 46% des désordres observés, sur les ouvrages en terre, ont pour origine l'érosion interne. En France, 70 cas critiques ont déjà été détectés. Malheureusement, à cause du phénomène d'érosion lui-même, la prévision du risque reste délicate et les autorités ont du mal à mettre en place des plans d'urgence. Compte tenu du risque potentiel, la compréhension des phénomènes à l'origine de l'érosion interne apparaît comme un enjeu scientifique majeur. L'objectif du travail de thèse est de mettre au point une méthode de diagnostic de l'érosion interne dans les digues et autres ouvrages en terre. Nous cherchons à établir un protocole qui nous guidera sur les essais à faire et qui nous aidera à interpréter les résultats. Cette méthode permettra d'expertiser les risques d'érosion interne. (i) Dans un premier temps, les principaux résultats des reconnaissances géomécaniques permettent d'identifier des couches géotechniques de sol de caractéristiques homogènes dans lesquelles sont réalisés les essais pressiométriques. L'analyse cluster des résultats des essais pressiométriques permet de déterminer le type, la nature de sol, ainsi que les caractéristiques géotechniques moyennes. (ii) Les résultats obtenus servent ensuite à déterminer si le sol est instable et quel est son seuil de sensibilité à l'aide du logiciel expert. Ce logiciel permet de calculer les critères d'érosion interne et de qualifier la sensibilité du sol à l'érosion en déterminant son seuil érosif. (iii) Un nouveau dispositif expérimental appelé l'Essai d'Erosion Transverse (Cross Erosion Test - CET) permet de caractériser expérimentalement le sol par rapport au risque d'érosion interne. Cet essai consiste à injecter de l'eau dans un forage, puis à pomper cette eau dans un forage parallèle au premier. L'eau extraite est chargée en particules produites par l'érosion interne du sol. L'érosion est estimée en se basant sur la mesure de l'évolution de la masse des particules recueillies au cours du temps. Les résultats montrent que cette expérience permet de caractériser l'érosion interne dans un sol spécifique. L'avantage de cette technique c'est qu'elle est transposable in situ et qu'elle peut être utilisée pour prédire le risque de suffusion dans les barrages et les digues. (iv) La validation de l'expérience est faite par un modèle numérique 3D réalisé avec le logiciel Comsol Multiphysics 4.3b. Ce modèle permet de montrer que sous les conditions de perméabilité observées, le gradient hydraulique se concentre autour des points d'injection et de pompage. Une modélisation de base en 2D est développée. Cette première approche permet de décrire le phénomène d'érosion et le transport des particules fines dans le milieu poreux du sol. / Internal erosion is the displacement of the fine particles of a soil under the action of an internal flow. This mechanism could be the origin of the damage on embankments and earth dams. In the word, 46% of failure observed on earthwork, are caused by internal erosion. In France, 70 critical cases have been detected. Unfortunately, due to internal erosion, the prediction of risk remains difficult. The understanding of this phenomenon of internal erosion appears as a major scientific challenge. The aim of this thesis is to develop a diagnostic method of internal erosion in dams and other earthworks. We seek to establish a protocol that will guide us to the test and help us to interpret the results. This method allows concluding about the risk of internal erosion appearance. (i) Initially, the main results of geotechnical survey allow identifying the layers of homogeneous characteristics of soil which are measured by the pressuremeter test. The results of cluster analysis of pressuremeter tests allow determining the localisation, the nature and the mean geotechnical characteristics of the soil. (ii) The results are used to determine the threshold of sensibility to internal erosion with the help of expert software. This software allows to calculate the criteria of internal erosion and to describe the sensitivity of the soil to internal erosion. (iii) A new experimental device is developed in our laboratory. This experiment called “Cross Erosion Test" (CET) allows determining the experimental resistance of the soil to the risk of internal erosion. The test consists of the injection, in a first drilling, of clear water and of the recovery, in another drilling, of water charged with eroded particles. For different initial state of the soils, it is possible to measure and to characterize the internal erosion by visualisation of the water flow and the measurement of the weight of the extracted eroded particles. The results show that this experience allows characterizing the internal erosion in a specific soil. The advantage of this technique is that it can be used in situ to predict the risk of suffusion in dams and dikes. (iv) A validation of experiences, with a 3D finite element method is carried out with the help of the Comsol Multiphysics 4.3b software. This model shows that under experiment hydraulic conditions, the hydraulic gradient is concentrated around the injection and the pumping. A finite element 2D model is developed to simulate erosion process. This approach describes the phenomenon of internal erosion and transport of fines particles in the porous medium of soils tested.

Le phénomène de suffusion correspond à la migration des particules fines érodées dans un milieu poreux sous l'action d’un flux hydraulique. Ce mécanisme pourrait être la cause principale des ruptures des remblais et des barrages en terre. De nombreuses études expérimentales ont été menées pour comprendre la suffusion et mettre en évidence les phénomènes couplés associés. Au niveau numérique, on trouve beaucoup de modèles analytiques et numériques mais rares sont ceux qui prennent en compte l’arrachement des particules, leur transport et leur dépôt éventuel (par un processus de colmatage ou pour des vitesses d’écoulement suffisamment faibles).En se basant sur une loi d’écoulement de type Darcy, sur la loi d’érosion de Papamichos et al (2001) modifiée par Bendahmane (2005)et sur une loi d’évolution de la viscosité de fluide (relation d’Einstein), nous développons une nouvelle approche hydromécanique de la suffusion. La résolution numérique du problème est réalisée par la méthode des éléments finis. Deux aspects sont présentés: le cadre général du modèle proposé (réalisé en 2D) avec une étude paramétrique et la validation du modèle avec une modélisation 3D réalisée avec le modèle géométrique de l’essai (Cross Erosion Test) (CET), ce qui permet de décrire à la fois les phénomènes d’érosion, de transport et de dépôt des particules érodées dans un milieu poreux saturé.Pour cela, nous avons choisi un outil de modélisation approprié à la représentation de phénomènes couplés tels que ceux mobilisés par l’érosion interne. Le logiciel Comsol-Multiphysics 3.4b, basé sur la méthode des éléments finis, a été retenu car il permet d’implanter relativement facilement de nouvelles équations constitutives des phénomènes. Dans un deuxième temps, le modèle a été validé sur des expériences de référence reflétant les phénomènes étudiés. Nous avons utilisé les résultats expérimentaux de Nguyen (2013) qui a mis au point un nouveau dispositif expérimental au laboratoire 3SR appelé " l'Essai d’Érosion Transverse " ou Cross Érosion Test. Deux types de sol ont été considérés : un sol gravelo-sableux(noté S3) et un sol limoneux de l’Isère (noté S4 et S4-a).Afin de discuter des aspects transitoire et dynamique du phénomène de suffusion, les vitesses d’écoulement, les gradients hydrauliques, la distribution de la porosité et de la perméabilité, la distribution de la concentration des particules érodées en suspension, et enfin la masse érodée sortantes ont mesurés au cours du temps et comparés aux résultats du modèle physique. / The phenomenon of Suffusion corresponds to the migration of fine eroded particles in a porous medium under the action of an internal water flow. This mechanism could be the main origin of damage in embankments and earth dams. Many experimental studies have been done to analyse the suffusion but actually understanding this phenomenon appears as a major scientific challenge. On the other hand, numerically, only a few studies were found which associate the mechanisms of detachment, transport and deposition of eroded particles.Based on the Darcy’s law model, the erosion model of Papamichos et al (2001) modified by Bendahmane (2005) and the fluid viscosity evolution law (Einstein relation), we developed a new hydro-mechanical approach of suffusion. The numerical solution of the problem is carried out by using a finite element method. Two aspects are presented: first, the general framework of the proposed model (implemented in 2D) with a parametric study and then the validation of the numerical model realized with a 3D study of the geometric model of the experimental study of Cross Erosion Test (CET). This model describes the three phenomena of suffusion: erosion, transport and deposition of eroded particles in a saturated soil.Furthermore, the numerical modeling of coupled phenomena of internal erosion is performed using the commercial software "Comsol Multiphysics-3.4b". However, to validate the model, experimental results on two types of soil: sandy gravel soil (S3) and the silt of Isère (S4 and S4-a) are used. These results are found in the PhD thesis of Nguyen(2013) who has developed a new apparatus called Cross Erosion Test (CET).Moreover, to discuss the transient and dynamic aspects of the phenomenon, the evaluation of the mass conservation with different figures of flow velocity, hydraulic gradient, hydraulic head, porosity and permeability distribution, concentration distribution, and of outgoing eroded particles, are measured over time and compared with experimental results.

No description available.

$aL’objectif de cette thèse est de comprendre le phénomène de suffusion (érosion interne) dans ses processus et ses conséquences sur le comportement mécanique des sols à faible cohésion par l’établissement d’une relation directe ou indirecte entre la stabilité interne de ces matériaux et leur critère de rupture. L’érosion interne vue essentiellement sous l’angle de la suffusion a été analysée par des modélisations numériques et des essais expérimentaux. Les principaux essais utilisés sont : - Ecoulements d’eau à travers un matériau granulaire dans une colonne. - Essais de cisaillement triaxial non drainés. Pour la modélisation numérique, PFC3D a été utilisé pour les simulations d’essais triaxiaux et la prise en charge du couplage fluide solide a été possible grâce à une option de type CFD (Coupling Fluid Dynamique). Ce code DEM a permis le développement d’une nouvelle méthode de génération d’échantillon granulaire à partir de sa courbe granulométrique et de sa porosité cibles. L’évolution des porosités, celle des courbes granulométriques et de l’état de stabilité interne avant et après écoulement ont été identifiées comme indicateur du niveau de suffusion. Les critères de stabilité de Kenney & Lau et de Kezdi ont permis d’apprécier la stabilité interne des échantillons. Les conséquences de la suffusion se sont manifestées par : Une diminution de l’angle de frottement interne dans le cas d’une suffusion interne; une augmentation de l’angle de frottement dans le cas d’une suffusion externe. Les comparaisons des résultats expérimentaux et numériques tant pour les écoulements que pour les essais triaxiaux ont présenté de bonnes concordances. Enfin, une proposition d’un nouveau critère de stabilité interne est établie, elle est déduite de l’analyse des relations entre le comportement mécanique et la stabilité interne, notamment la contractance du matériau. / The objective of this thesis is to understand the phenomenon of suffusion (internal erosion) in its process and its effects on the mechanical behavior of loose soil by establishing a direct or indirect relationship between the internal stability of these materials and failure criterion. Numerical modeling and experimental tests analyzed the internal erosion in terms of the suffusion. The main tests used are : - Water flows through a granular material in a column. - Undrained shear triaxial tests. For numerical modeling, PFC3D was used for simulations of triaxial tests and the treatment of solid fluid coupling was possible through a CFD type of option (Coupling Fluid Dynamics). This DEM code has allowed the development of a new method of generating granular sample from its grading curve and its target porosity. The evolutions of porosity, grading curves and the state of internal stability before and after discharge were identified as indicator of suffusion. Kenney & Lau and Kezdi stability criteria were used to assess the internal stability of samples. The consequences of suffusion were expressed by: A reduction of the internal friction angle in the case of an internal suffusion; an increase of the friction angle in the case of an external suffusion. Comparisons of experimental and numerical results for both flows and triaxial tests have shown good agreements. Finally, a proposal of a new criterion of internal stability is established, it is deducted from the analysis of the relationship between the mechanical behavior and internal stability, including the material contractance.