Le travail concerne la simulation des ouvrages de stockage des déchets radioactifs en formation géologique profonde installés dans l’argilite du Callovo-Oxfordien.Il comporte 3 parties : La première présente le modèle SC2D (Swelling Creep Diagenesis Damage) utilisé pour évaluer les comportements de l’argilite ; La deuxième présente les applications du modèle SC2D pour décrire les zones EDZ des ouvrages excavés. Ces zones sont caractérisées par les domaines où les coefficients d’endommagement sont unitaires. Les coefficients d’endommagement sont évalués à partir d’une extension de la relation de Mazars en considérant les déformations d’extension dans les directions de diagenèse ; La troisième partie concerne les effets différés autour des ouvrages après l’excavation, ils sont représentés par les vecteurs vitesses de déformation dont les amplitudes sont gouvernées par les niveaux d’endommagement et les déviateurs des contraintes. La comparaison entre des mesures de convergences in situ et des résultats modélisés permet de valider le vecteur vitesse de déformation. / The study aims to evaluate the radioactive waste storage structures in argillite Callovo- Oxfordien : The first part presents the model SC2D (Swelling Creep Diagenesis Damage) used to simulate the behaviour of argilite, this model is based on special hydro-mechanical characteristics of the argillite Callovo-Oxfordien; The second part shows the applications of SC2D model to describe anisotropic damaged areas around the excavated structures. These EDZ zones are determined when the damage coefficients are unitary and the discontinuities activate when damage coefficient reaches unity, damaged coefficients are calculated from an extension of Mazars by considering the extension deformation in the directions of diagenesi; the final part discussess the delayed effects around the structures after excavation; they are represented by the deformation velocity vectors whose amplitudes are governed by the level of damage and the deviatoric stress. The comparison between in situ convergence measurements and modeled results validates the deformation velocity vectors.

Afin de prévoir et réglementer l'impact environnemental des activités humaines telle que l'exploitation minière de l'uranium et des déchets radioactifs, il est nécessaire de comprendre le comportement des actinides, la migration des radionucléides dans l'environnement et analyser leur interaction avec l'argile minérale. Le comportement des actinides dans le sol concerne principalement les interactions et les adsorptions dans les surfaces qui changent les formes des éléments radioactifs et réduisent la mobilité des actinides dans les systèmes naturels. Par conséquent, il est important de comprendre comment les actinides interagissent avec l'argile minérale et étudier le processus fondamental de précipitation de surface.L'uranium est la teneur en métal lourd la plus prépondérante des déchets ultimes dans le cycle du combustible nucléaire ( > 95 % de UO2 ). En outre, l'uranium est un contaminant majeur dans le sol, le sous-sol et les eaux souterraines en raison de l'activité humaine. Dans des conditions environnementales standards, la forme chimique la plus stable de U ( VI ) est l' ion uranyle UO22+, qui est potentiellement très mobile et se combine facilement avec la matière organique et inorganique.D’un autre côté, le dioxyde de carbone est un important gaz à effet de serre, réchauffant la surface de la terre à une température plus élevée en réduisant vers l'extérieur le rayonnement. Toutefois, des problèmes peuvent se produire lorsque la concentration atmosphérique des gaz à effet de serre augmente. Des quantités importantes d'émissions de dioxyde de carbone ont été produites depuis la révolution industrielle, ce qui est derrière l’important réchauffement climatique et l’augmentation du niveau de la mer. Les minéraux argileux sont d'une grande importance pratique, particulièrement dans le stockage du dioxyde de carbone en raison de sa perméabilité hydraulique et sa capacité de conserver les particules mobiles. Nous avons choisi la kaolinite et la montmorillonite comme des prototypes de minéraux argileux de type 1:1 et 2:1. Les méthodes de Monte Carlo (MC) et la dynamique moléculaire (MD) ont été utilisées dans ce travail afin de comprendre le comportement d'adsorption des radionucléides et de dioxyde de carbone à la surface des argiles. Dans cette thèse, nous étudierons d'abord l'adsorption de l’uranyle à la surface de la kaolinite par le biais de la dynamique moléculaire et de la Monte Carlo. Plusieurs sites d'adsorption ont été modélisés en considérant des défauts de surface dans le but d'avoir des complexes inter ou externe-sphère. Ensuite, les adsorptions des particules d'uranyle sur des surfaces de la montmorillonite en présence des différents ions seront effectuées. L'énergie d'interaction entre des feuilles de montmorillonite et le travail d'adhérence entre le radionucléide et surface MMT seront également discutés. Enfin, nous présenterons l’étude du comportement d'adsorption de dioxyde de carbone dans le MMT, en détaillant au passage les propriétés thermodynamiques, structurales et dynamiques. / In order to predict and regulate the environmental impact of human activities such as uranium mining and radioactive waste disposal, it is necessary to understand the behavior of actinides in the environment because their interaction with clay mineral is an important factor to control the migration of radionuclide in the environment. The behavior of actinides in the soil is mainly the surface adsorption interactions, which change the forms of radioactive elements and reduces the mobility of actinides in the natural systems. Therefore, it is important to search how the actinides interact with clay mineral such as the fundamental process of surface precipitation. Uranium is the predominant heavy metal content of the final waste in the nuclear fuel cycle (>95% UO2). In addition, uranium is a major contaminant in the soil, subsurface and groundwater as a result of human activity. Under standard environmental conditions, the most stable chemical form of U(VI) is the uranyl ion UO22+, which is potentially very mobile and readily complexates with organic and inorganic matter. On the other hand, Carbon dioxide is an important greenhouse gas, warming the earth’s surface to a higher temperature by reducing outward radiation. However, problems may occur when the atmospheric concentration of greenhouse gases increases. Amounts of carbon dioxide were produced since the industrial revolution, which is behind the significant global warming and rising sea level. Clay minerals are of great practical importance here, in storage of carbon dioxide due to its hydraulic permeability and ability to retain mobile species. We have chosen kaolinite and montmorillonite as prototypes of clay minerals of 1:1 and 2:1. Classical Monte Carlo (MC) and molecular dynamics (MD) methods have been used in this work in order to understand the adsorption behaviour of radionuclide and carbon dioxide in clays surface. In this thesis, we will investigate –first- the adsorption of uranyl on kaolinite surface by means of Monte Carlo and Molecular Dynamics simulation methods. Several adsorption sites have been modeled by considering surface defects in order to have inner or outer-sphere complexes. Then, the adsorptions of uranyl species onto Montmorillonite surfaces in the presence of different counterions will be performed. Interaction energy between Montmorillonite sheets and work of adhesion between the radionuclide and MMT surface will be discussed as well. Finally, we will study the adsorption behavior of carbon dioxide in MMT, and investigate at the same time thermodynamics, structural and dynamical properties.

Cette thèse s’inscrit dans le cadre des recherches et vise à contribuer à la compérension du comportement mécanique d’une roche tendre (i.e. cataclastique grès). La roche étudiée se trouve dans une fondation d’un barrage, situé dans le sud-ouest de la chine. L'accent de des travaux réalisés dans cette thèse a concerné le comportement à long terme du matériau étudié.Pour caractériser le comportement mécanique et identifier les principaux paramètres qui le gouvernent, il était nécessaire d’aborder en première lieu avec une approche expérimentale. De nombreux essais ont été réalisé, tel que des essais de compression, des essais de fluage triaxiale et des mesures de perméabilité. On constate que la réponse mécanique du matériau dépends fortement la microstructure et les propriétés mécaniques ainsi que les contraintes appliquées. En se basant sur les résultats expérimentaux, un modèle unifiée est proposé pour décrire le comportement à court terme et à long terme du matériau étudié. Dans le modèle proposé, les principales caractéristiques observées dans les expériences sont prises en compte, en particulier des déformations irréversibles, la dépendance de la pression, la transition volumétrique entre le compaction et la dilatation ainsi que le fluage. La performance du modèle proposé est tout d’abords validé sur les essais homogènes. A la fin, une validation plus poussée de ce modèle a été réalisée dans le cadre d’applications à étudier la performance de la fondation. / The short term and long term mechanical behaviors of cataclastic sandstone are investigated in this dissertation. The studied material is cored from the fault zone in the dam foundation of a hydropower station in southwestern China. The emphasis is put on the long-term behavior of studied material. For this, conventional compression tests, multi-stage triaxial creep tests as well as the permeability measurements are performed. The mechanical behavior of studied rock is derived based on the experimental investigation. The mechanical responses of studied material are not only functions of the material properties, but are also related to applied stress. Based on the experimental investigation, a unified creep model is proposed to describe the instantaneous and time-dependent behavior of cataclastic sandstone. Main features observed in experiments are taken into account, in particular irreversible deformations, pressure dependency, volumetric transition between compaction and dilatancy and creep. The performance of the model is examined by comparing numerical simulations with test data in representative load paths. Finally, the model is applied to analyze the long term stability of a foundation of hydropower station.

Des argilites sont étudiées comme barrières géologiques potentielles pour le stockage des déchets radioactifs. Ce travail est consacré à l’étude expérimentale et à la modélisation du comportement mécanique des argilites de Meuse- Haute Marne avec prise en compte des influences de la saturation, de la composition minéralogique et de la température.Un grand nombre d’études expérimentales ont été effectuées sur le comportement mécanique. Ce travail propose de réaliser de nouvelles expériences pour compléter les données existantes. L’accent est mis sur les influences de la saturation, et la composition minéralogique et de la température. Des essais de micro-indentation et de mini compression sont effectués sur des échantillons de petites tailles sous différentes humidités relatives et températures. Ces essais permettent d’obtenir une meilleure distribution de saturation au sein des échantillons. Les variations des propriétés élastiques et de la résistance mécaniques des argilites sont étudiées en fonction de saturation et de température. Des échantillons provenant de différentes profondeurs sont étudiées et les influences de la composition minéralogique sur les propriétés mécaniques sont mises en évidence.Ce travail propose d’utiliser un modèle micromécanique basé sur une méthode d’homogénéisation non linéaire, ceci afin de prendre en compte des influences de la composition minéralogiques comparant avec des modèles phénoménologiques. L’argilite est schématisée par un matériau à trois phases : la matrice argileuse élastoplastique, des grains de calcite et de quartz élastiques linéaires. Des comparaisons entre les simulations et les données expérimentales sont présentées. / Clayey rocks are largely investigated for the feasibility study of geological storage of radioactive wastes. In that context, the clayey rocks are subjected to solicitations of different origin: stress variation due to cavity excavation, change of water saturation degree (drying and resaturation), variation of temperature due to exothermal wastes, and chemical reaction processes. Therefore, it is necessary to develop coupled THM-C modelling for the short and long term feasibility analysis of geological storage facilities. The present work is devoted to experimental investigation and modelling of mechanical behaviour of Meuse- Haute Marne (MHM) argillites by taking into account influences of saturation, mineralogical composition and temperature.A number of laboratory investigations have been performed for studying mechanical behaviour of MHM argillites. We propose to realize new investigations to extend and complete existing data. The emphasis is put on the influences of saturation and temperature. Micro-indentation and mini compression tests are conducted on small size samples under different values of relative humidity and temperature. Argillite samples from different depth are then investigated and the influences of mineralogical composition on mechanical properties are shown.In the present work, we propose to apply a homogenization-based micromechanical model. Comparing with the macroscopic models, it can describe influences of mineralogical compositions on macroscopic behaviours. The argillites are considered as a three phase composite: clay matrix, calcite and quartz grains. Comparisons between numerical modelling and experimental data will be presented.

L'objectif de la thèse est d'explorer les comportements mécaniques effectifs des matériaux poreux impliqués dans les problèmes de pores ou d'inclusion, avec pour objectif principal d'homogénéiser un matériau aussi hétérogène que possible, plus d'informations de géométrie physique à une échelle donnée. Nous examinons séparément ici quatre microstructures typiques de matériaux ressemblant à des roches couvrant différents constituants: nous avons cherché à déterminer comment la présence de pores et une inclusion de morphologies différentes influent sur les comportements macroscopiques élastiques, débit plastique et dépendant du temps. Un modèle d'homogénéisation par étapes pour les comportements élastoplastiques et viscoplastiques est proposé, basé sur une méthode d'homogénéisation par transformée de Fourier rapide, afin d'améliorer le comportement local depuis la transition à la micro-échelle vers la méso-échelle, puis vers la macroscopie. Une série de simulations est calculée et comparée. Nous portons une attention particulière au comportement anisotrope des matériaux à la fois pore et inclusion configurée à la même échelle et à différentes échelles Nous fournissons une solution de référence pour ce type de matériaux couvrant l’effet des géométries de pores et d’inclusion en raison de la limitation à l’obtention d’un critère analytique. En plus nous proposons un modèle d’endommagement plastique multi-échelle. les comparaisons entre les résultats numériques et les données expérimentales montrent que ce modèle peut bien caractériser son comportement en défaillance avec des micro-structures complexes. / The objective of this thesis aims to explore the effective mechanical behaviors of porous materials involved with pore or inclusion problem. The key point to this problem is to homogenize such a highly heterogeneous material incorporating more physical geometry information at a given scale. In this work, four typical microstructure of rock-like materials are respectively considered here covering different constituents. We sought to determine how the presence of pores and inclusion with different morphologies influence the macroscopic elastic, plastic flow and time-dependent behaviors. For this propose, a two-step homogenization model for elastoplastic and viscoplastic behaviors is proposed based on a Fast Fourier Transform homogenization method to upscale the local behavior from the micro-scale transition to meso-scale and then to macroscale. A series of simulations are computed and compared. Particular attentions are focused on the anisotropic behavior for the materials with both pore and inclusion configured at same and different scales. We provide a reference solution for this kinds of materials covering the effect of pore and inclusion geometries due to a limitation to obtain an analytical criterion. Moreover, we proposed a damage plastic model with multi-scale characters. Comparisons between numerical results and experimental data show that this model can well characterize its failure behavior with complex micro-structures.

Dans ce travail, un nouveau modèle de liaison avec critère de rupture de cisaillement non linéaire est d'abord proposé et mis en œuvre dans PFC pour décrire le comportement mécanique de matériaux granulaires cohésifs isotropes tels que le grès. Un grand nombre de tests de compression ont été effectués sur des échantillons 3D. Les résultats montrent que l'effet de la pression de confinement sur la résistance à la compression et le modèle de rupture est bien décrit par ce modèle de liaison. Les effets du chemin de chargement et de la contrainte principale intermédiaire sur la déformation et la rupture ont également été étudiés. Après avoir étendu les paramètres d'élasticité et de résistance, deux modèles de liaison sont couplés pour caractériser l'anisotropie de la résistance et de la déformation dans les matériaux cohésifs anisotropes tels que les roches sédimentaires. Une série de test de compression avec différents chemin de chargement ont été effectués et les résultats numériques sont en bon accord avec les données expérimentales, en particulier pour les propriétés de réponse élastique et de résistance. La transformation du mode de défaillance entre la couche de faiblesse et la matrice de roche sous différentes pressions de confinement a été bien décrite. De plus, le modèle de liaison proposé a été étendu pour étudier le processus de fracturation hydraulique dans les matériaux cohésifs. Un processus de propagation de fracture hydraulique représentatif a été présenté. Enfin, nous avons développé une méthode de reconstruction dimensionnelle de l'irrégularité de la forme des blocs et étudié ses effets sur les impacts des blocs en utilisant une approche basée sur l'énergie. / In this work, a new bond model with nonlinear shear failure criterion is first proposed and implemented in PFC for describing mechanical behavior of isotropic cohesive granular materials such as sandstone. A large number of compression tests have been performed on 3D samples. The results show that the effect of confining pressure on compressive strength and failure pattern is well described by the proposed bond model. Effects of loading path and the intermediate principal stress on deformation and failure have been also investigated. After further extending parameters on elastic and strength, two bond models: the proposed bond model and the improved smooth joint model, are coupled to character the anisotropy of strength and deformation in anisotropic cohesive materials such as sedimentary rocks. A series of conventional triaxial compression tests with different loading paths have been performed and numerical results are in good agreement with experimental data, especially for elastic response and strength properties. The failure mode transformation between weakness layer and rock matrix under different confining pressures has been well described. Furthermore, the proposed bond model has been extended to study hydraulic fracturing process in cohesive materials. A representative hydraulic fracture propagation process has been presented. Influences of different factors, such as confining pressure, fluid viscosity and fluid injection rate, on hydraulic fracture extension have been investigated. Finally, we have developed a dimensional reconstruction method of block shape irregularity and studied its effects on block impacts using an energy based approach.

- / -

L'objectif de cette thèse est de développer une approche micromécanique pour modéliser l'endommagement et le frottement dans des géomatériaux initialement anisotropes. Dans ce but, deux modèles d'endommagement micromécaniques en 3D ont été proposés. L'estimation du tenseur de Hill pour les microfissures sphéroïdales 3D dans une matrice isotrope transverse est le cœur de ces deux modèles. Les propriétés élastiques effectives du milieu fissuré sont déterminées par trois schémas d'homogénéisation (schéma Dilué, schéma MT, schéma PCW) et analysées par une série d'études de sensibilité. Le premier modèle d'endommagement formulé dans le cadre de la thermodynamique irréversible est proposé pour prendre en compte des interactions entre l'anisotropie initiale et la fissuration induite. Ensuite, le deuxième modèle micromécanique est développé pour modéliser l'endommagement et le glissement plastique par frottement dans des fissures fermées. Les évolutions des déformations inélastiques sont provoquées par deux phénomènes couplés: le glissement avec frottement et l'endommagement induit par la propagation des microfissures. Elles sont décrites par des critères spécifiques locaux de l'endommagement et du frottement en fonction des forces thermodynamiques conjuguées. Les principales caractéristiques du comportement mécanique des géomatériaux initialement anisotropes (la croissance des microfissures induites, les déformations irréversibles, boucles d'hystérésis causées par les effets du frottement qui apparaissent sur les lèvres des microfissures fermées, le couplage entre l'anisotropie inhérente et celle induite, les effets de la pression de confinement et de la dilatance volumétrique) sont prises en compte. / The objective of this thesis is to develop a micromechanical method for modeling induced damage and friction sliding for a class of initially anisotropic quasi-brittle materials. To this end, two three-dimensional micromechanical damage models have been developed respectively for open and closed microcracks. The key step of these two models is an efficient numerical method to estimate Hill's tensor for three dimensional spheroidal microcracks in a transversely isotropic solid matrix. Then, the effective elastic properties of cracked material are determined through a rigorous up-scaling procedure based on three homogenization methods (Dilute scheme, MT scheme, PCW scheme) and analyzed by a series of sensitivity studies. Further, a damage model based on the irreversible thermodynamics framework is proposed to account for interactions between initial anisotropy and induced growth of cracks. Furthermore, the second micromechanical model for frictional closed cracks is developed for the description of induced damage and frictional sliding in initially anisotropic materials. The evolutions of inelastic strains due to frictional sliding and induced damage by the propagation of microcracks are coupled and described by specific local friction and damage criteria as functions of conjugate thermodynamic forces. The main features of mechanical behaviors of initially anisotropic quasi-brittle materials are taken into account, such as induced growth of micro-cracks, irreversible deformation and hysteretic loops due to friction effects of crack surfaces, coupling between inherent anisotropy and induced anisotropy, effects of confining pressure and volumetric dilatancy.

L'objectif de ce travail est d'étudier les effets de microstructures sur les comportements mécaniques macroscopiques des géomatériaux hétérogènes. Profitant de la méthode numérique basée sur la Transformée de Fourier Rapide (TFR), les microstructures complexes de géomatériaux peuvent être simulées aussi proche que possible de la vraie microstructure. Avec des calculs en champ complet, les contraintes et les déformations locales sont fournies à l’échelle microscopique. Cette caractéristique permet de surmonter les hypothèses fortes sur la microstructure par des approches traditionnelles d'homogénéisation. Matrix-inclusion composites avec microstructures différentes sont d'abord étudiés. Les influences de la forme, la taille, la distribution et l'orientation des inclusions sur le comportement macroscopique sont prises en compte. Ensuite, ce modèle numérique basé sur TFR est appliqué à l'argilite du Callovo-Oxfordien qui est traitée comme une matrice élastoplastique renforcé par des grains élastiques de quartz et de calcite. Avec une règle d'écoulement non associé dans la matrice, les effets de grains sphériques de minéraux sont explicitement pris en compte. Ce modèle est en outre étendu par la considération du processus de détérioration progressive due à la croissance de microfissures. Après cela, les effets de pores (formes, tailles, orientations et distributions) sur le comportement effectif de matériaux poreux sont considérés. A titre d'exemple, un géomatériau poreux typique le grès a été étudié en détail. En outre, la détérioration autour des pores et des régions interagies sont simulées avec un critère simple d’endommagement pour l'évolution des pores et le phénomène de dégradation est entièrement exposé. Les comparaisons entre les résultats numériques et les données expérimentales vérifient l'efficacité et la précision de cette méthode numérique basée sur TFR pour les géomatériaux hétérogènes. / The objective of this work is to study the effects of microstructures on the macroscopic mechanical behaviors of heterogeneous geomaterials. Taking advantage of the Fast Fourier Transforms (FFT) based numerical method, the complex microstructures of geomaterials can be considered as close as possible to the real microstructure and a full field micromechanical analysis of strain and stress can be provided at the microscopic scale. This feature overcomes the strong assumptions on microstructure by traditional homogenization approaches. Matrix-inclusion type composites with different microstructure are firstly investigated. The influences of inclusion shape, size, distribution and orientation on the effective behavior are fully studied. Then it is applied to the Callovo-Oxfordian argillite which is modeled as an elastoplastic matrix reinforced by linear elastic quartz and calcite spherical grains. With a non-associated plastic flow rule in the matrix, the effects of mineral grains are explicitly taken into account. The model is further extended by introducing the progressive damage process due to the growth of microcracks. After that, the pore effects (shapes, sizes, orientations and distributions) on the effective behavior of porous materials are considered with a specific application to the typical porous geomaterial sandstone. Furthermore, the deterioration around the pores and the interacted regions are simulated with a simple damage criterion for pore evolution and the degradation phenomenon is fully exhibited. The comparisons between the numerical results and experimental data verify the efficiency and accuracy of the FFT based method for heterogeneous geomaterials.

Ce travail porte sur des modélisations, par techniques d’homogénéisation, de la plasticité des matériaux poreux ductiles. Différentes classes de milieux poreux (avec ou sans compressibilité plastique de la matrice) sont considérées. Pour une matrice de von Mises, nous avons obtenu une expression du critère macroscopique de plasticité des milieux contenant des cavités sphéroïdales aléatoirement orientées. Une extension du critère de Gurson est proposée pour des matériaux dont la matrice est plastiquement compressible et obéit à un critère elliptique. Les résultats obtenus sont appliqués et illustrés dans le cas des milieux à double porosité. Un accent particulier a été ensuite mis sur les géomatériaux poreux. L’argilite du Callovo Oxfordien, géomatériau auquel on s’est intéressé, est modélisé comme un milieu hétérogène, constitué à l’échelle mésoscopique d’une matrice élastoplastique poreuse et d’inclusions élastiques. La transition micro-méso aboutissant à la loi de la matrice argileuse est réalisée pour une phase solide de type Drucker-Prager, à écoulement associé ou non, et pour des cavités sphériques. Tandis que le passage méso-macro s’inscrit dans une démarche d'homogénéisation incrémentale de Hill. Les capacités prédictives du modèle basé sur une phase solide à loi d’écoulement non associée ont été démontrées à travers une large validation expérimentale. Enfin, par souci de simplification du modèle micro-macro proposé pour l’argilite, une version simplifiée du modèle est proposée en s’appuyant sur l’hypothèse d'inclusions rigides. L'introduction d'une règle d'écoulement macroscopique non associée a conduit à des prédictions d’excellente qualité pour l’argilite. / This work focused on the plasticity modeling of ductile porous materials. Different classes of porous media (with or without plastic compressibility of the matrix) are considered. Concerning the case of incompressible matrix obeying to the von Mises criterion, by using limit analysis combined with homogenization techniques, we obtained, a closed form expression of the macroscopic yield function for materials containing randomly oriented spheroidal cavities. An extension of the Gurson criterion is proposed for materials whose matrix is plastically compressible and obeys to an elliptic criterion. The obtained results are applied and illustrated in the case of double porous media. Then, an emphasis is put on porous geomaterials. The Callovo Oxfordian argilite, a geomaterial which is of a primary interest in this study, is modeled as a heterogeneous medium made, at the mesoscale, of a porous elastoplastic matrix and elastic inclusions. The micro-meso transition leading to the constitutive law of the clay matrix is performed for an associated or a non-associated Drucker-Prager solid phase and spherical cavities. The meso-macro transition is performed by means of a Hill incremental homogenization method. The predictive capabilities of the micro-macro model based on a non-associated solid phase have been demonstrated via a wide experimental validation. Finally, a simplified version of the micro-macro model is formulated by assuming rigid inclusions. The introduction of a macroscopic non-associated flow rule allows to obtain excellent predictions for the argillite.