De nombreux matériaux du génie civil, comme le béton ou les enrobés bitumineux, sont des pâtes granulaires, c'est-à-dire des suspensions dont la concentration volumique en inclusions solides est très élevée. Ce travail porte sur la modélisation et la simulation numérique discrète du comportement et du compactage de pâtes granulaires 2D non browniennes et non colloïdales constituées de grains élastiques rugueux de forme circulaire, baignant dans un liquide newtonien de très grande viscosité. Après avoir introduit un modèle de pâte granulaire dans lequel le seul effet du liquide est d'induire des interactions de lubrification entre paires de grains proches, nous présentons une méthode numérique quasi-statique permettant de simuler le comportement de ces systèmes lorsqu'ils sont soumis à un état de contrainte donné. L'étude de l'état d'équilibre adopté par une pâte granulaire lubrifiée lorsqu'on lui applique une pression isotrope montre que, dans la limite où la rugosité des grains est nulle, un tel système se comporte comme un assemblage sec de grains non frottants. C'est la raison pour laquelle nous étudions ensuite le comportement de ces assemblages en 3D lorsqu'on les soumet à des chargements statiques anisotropes ainsi qu'à des essais de cisaillement simple dans la limite quasi-statique. Nous montrons en particulier que les coefficients de frottement interne statique et dynamique sont non nuls et coïncident, que ces assemblages ne présentent aucune dilatance et qu'ils satisfont un critère de rupture de Lade-Duncan. Le comportement d'une pâte granulaire lubrifiée 2D faite de grains de rugosité finie soumise à un cisaillement maintenu est ensuite analysé et nous comparons nos résultats avec ceux obtenus en cisaillant des assemblages granulaires secs 2D. Finalement, nous étudions de manière détaillée les paramètres influençant le compactage d'une pâte granulaire lubrifiée et mettons en évidence la réduction de la mobilisation du frottement intergranulaire accompagnant la densification du matériau.

Cette thèse porte sur un système dissipatif modèle, les gaz granulaires. Au moyen de méthodes issues de la théorie cinétique et des processus stochastiques, nous avons cherché à déterminer les propriétés statistiques d'observables globales, mimant ainsi la démarche qui prévaut l'équilibre. Parmi celles-ci, l'énergie qu'il faut fournir au gaz granulaire pour le maintenir dans un état stationnaire joue un rôle central. Elle est d'intérêt expérimental, mais elle est aussi pertinente comme mesure de la distance à l'équilibre, et c'est à ce titre qu'elle est récemment intervenue dans le contexte des relations de fluctuation. Nous avons complètement caractérisé la distribution de cette énergie injectée dans le système. Nous avons montré qu'elle ne pouvait être invoquée dans les relations de fluctuation, mais aussi que la distribution de cette grandeur macroscopique permettait d'obtenir des informations sur la dynamique microscopique. Plusieurs questions se sont développées autour de ce travail, dont la pertinence de la notion d'entropie dans un système dissipatif ou les limitations intrinsèques des relations de fluctuations.

La raréfaction des combustibles fossiles et la prise en compte de plus en plus importante des enjeux environnementaux incite à rechercher de nouvelles pistes de production, stockage et transport de l’énergie. Le développement de la filière hydrogène dépend des améliorations de chacun de ces principaux maillons qui sont la production, le stockage, la distribution et l’utilisation de l’hydrogène.Cette thèse se consacre au stockage de l’hydrogène par absorption dans des matériaux solides pulvérulents composés d’hydrures intermétalliques. Les performances de ce type de stockage dépendent grandement de la gestion des échanges de chaleurs imposés par le caractère exothermique de la réaction d’hydruration (absorption d’hydrogène par le matériau) et respectivement le caractère endothermique de la réaction de déshydruration (désorption d’hydrogène par le matériau). Cette nécessité a guidé la conception des réservoirs vers une architecture interne cellulaire. Si la thermique des réservoirs est très étudiée, les mécanismes de gonflement-dégonflement cycliques des hydrures lors de l’absorption et la désorption d’hydrogène sont assez peu abordés dans la littérature. Pourtant, les interactions mécaniques de la poudre hydrure sur les parois des cellules peuvent remettre en cause l’intégrité du réservoir. La mécanique des milieux granulaires constitue ainsi un apport essentiel à l’analyse et à termes à la prédiction du comportement mécanique du lit de poudre hydrure au sein du réservoir.Cette thèse présente en premier lieu la caractérisation expérimentale du comportement mécanique d’un hydrure intermétallique Ti-Cr-V à différents états de cyclage sous hydrogène en cellule instrumentée en incluant l’analyse des mécanismes de gonflement-dégonflement. L’analyse du comportement mécanique est abordée par une approche de type mécanique des milieux granulaires basée principalement sur la méthode des éléments discrets (DEM). Ces travaux ont permis d’identifier les caractéristiques et le comportement mécanique de l’hydrure (granulométrie, compression en matrice, écoulement) ainsi que l’évolution de paramètres liés au cyclage. Les effets du mécanisme de réarrangement opérant lors du cyclage à l’échelle des grains d’hydrure ont pu être identifiés expérimentalement. En particulier, la variation de volume cyclique du matériau entraine un tassement progressif du lit de poudre, une diminution de la porosité et une augmentation des niveaux de contraintes. Ce comportement a été reproduit numériquement par des simulations impliquant des particules discrètes de formes sphériques et de formes agrégées complexes. / It is harder and harder to find fossil fuels and it becomes therefore more and more expensive. Furthermore the environmental impacts of energy are more and more taken into account. For those reasons new solutions to produce, store, transport energy are currently researched and/or developed. Among those solutions, hydrogen system could be a good solution depending and how its main parts (production, storage, supply and use) are improved.This PhD is relative to hydrogen storage in intermetallic hydride compouds. Hydrogen tanks using these kind of materials have to allow good thermal control of the hydride powder especially because of the exothermicity of the hydriding and endothermicity of the unhydriding reaction. That's why the internal architecture of these vessels are often cellular, the cell walls playing the role of a thermal conductor. If the thermal aspects relative to these tanks are often studied, it is not the case of the mechanical phenomenon induced by the swelling and shrinking of the grains during absorption and desorption of hydrogène by the material. However, the mechanical interaction between the powder and the cell walls could endanger the tank.This PhD consists of two main parts. The first part is a mechanical study of the behavior of a Ti-Cr-V compound while the second is a mechanical modelling and analysis mainly by the Discretes Elements Method (DEM). Thanks to this work the main features and mechanical properties of the hydride (granulometry, matrix compaction, granular flow, powder density …) and their evolution due to cycling were measured. The influence of hydride grains rearrangement induced by hydride breathing were analysed experimentally. It lead to a progressive densification of the powder bed in instrumented cells that resulted in a decrease of porosity and an increase of stresses on cell walls. This effect was reproduced in discrete elements simulations of spherical and clump particles.

Nous étudions la rhéologie d'un ensemble de grains sphériques et frottants, par la simulation numérique, à l'échelle des grains, d'écoulements de cisaillement sous une contrainte normale P contrôlée, en présence d'un liquide interstitiel. En faible teneur, ce liquide se présente sous forme de ménisques intergranulaires qui transmettent des forces capillaires attractives ; s'il sature l'espace intergranulaire, on s'intéresse alors à l'écoulement de Stokes de la suspension dense ainsi constituée, où dominent les forces visqueuses. Les assemblages de grains secs constituent un système de référence aux propriétés mécaniques bien connues, en particulier l'approche de l'état critique de la mécanique des sols dans la limite quasi-statique. L'effet des ménisques capillaires qui joignent les grains en présence d'un liquide en faible saturation (régime pendulaire) est étudié pour les taux de cisaillement allant du régime quasi statique au régime inertiel. La rhéologie est caractérisée par le frottement interne apparent, la compacité de l'assemblage, les différences de contraintes normales et diverses variables internes, fonctions de deux paramètres de contrôle adimensionnés : le nombre inertiel I et la pression réduite P*, qui compare les forces de confinement à l'adhésion dans les contacts. Notre étude concerne les états homogènes, ce qui exclut les états de cisaillement localisés observés à faible P*, de l'ordre de 0,1. Le coefficient de frottement interne augmente de 0.35 (cas sec) à 0.9 environ pour P*=0.4, tandis que la compacité décroît de 0.59 à 0.52. L'important effet des forces capillaires sur la rhéologie, sensible pour des P* de plusieurs unités, est relié à la texture anisotrope des contacts et des ponts liquides. Lorsque P* décroît, nombre de contacts cohésifs sont maintenus pour des intervalles de déformation de plusieurs unités, survivant aux effets de rotation et de cisaillement de l'écoulement, et forment des amas percolants dans le système entier. Les résultats sont modérément sensibles à la saturation dans le régime pendulaire, mais fortement affectés par l'hystérèse de la conformation des ménisques. En présence de forces visqueuses et non plus capillaires, une version simplifiée de la dynamique stokésienne est adoptée dans laquelle les forces de lubrification entre proches voisins, supposées dominantes, sont les seules interactions hydrodynamiques. La rhéologie est fortement influencée par les contacts intergranulaires directes, qu'autorise la coupure à courte distance de la singulérité de lubrification du fait de la rugosité de surface des particules. Le même état critique que celui des grains secs est approché dans la limite quasi-statique. Nous discutons de lois rhéologiques exprimées en fonction du nombre visqueux qui remplace alors le nombre inertiel, et de la divergence de la viscosité effective à l'approche de la compacité critique en écoulement permanent, ou de la compacité maximale des assemblages aléatoires pour les configurations isotropes désordonnées / We numerically simulate the shear flow of dense assemblies of 3D frictional spherical grains under a fixed normal stress P in steady-state, either in the presence of a small amount of an interstitial liquid, which gives rise to capillary menisci and attractive forces, or in the fully saturated state, when the mechanical properties of suspensions in Stokes flow are controlled by hydrodynamic and contact forces. Dry grain assemblies are used as a reference system for which the rheological properties - in particular the approach to the critical state – are rather well known and can be measured with good accuracy. A non-saturating wetting fluid creates capillary attractive intergranular forces, the effects of which on the rheology are investigated in the pendular state, with shear rates ranging from quasistatic to inertial regimes. The system behavior is characterized by the dependence of internal friction coefficient, solid fraction, normal stress differences and internal state parameters on two dimensionless control parameters: the inertial number, I and the reduced pressure, P*, comparing confining forces to contact tensile strength. We focus on steady homogeneous flows, excluding localized flow patterns which we observe to occur for low P* (of order 0.1). The apparent internal friction coefficient increases to 0.9 in the quasistatic limit for P*=0.4, from its dry value 0.35, while solid fraction decreases from 0.59 to 0.52. We relate the significant effect of capillary forces on the macroscopic behavior of the system, up to P* values of several unities, to fabric anisotropy parameters of contact and distant interactions. As P* decreases, many cohesive contacts are observed to survive the tumbling motion associated to the shear flow, and their average age exceeds the reciprocal shear rate. Corresponding clusters of grains with enduring capillary bonds gather a large proportion of grains and percolate through the sample. The results are shown to be moderately sensitive to saturation within the pendular range, yet rather strongly affected by the hysteretic nature of liquid bridges. In the presence of viscous forces, assuming lubrication effects to dominate the hydrodynamic interactions, we adopt a simplified version of the (overdamped) Stokesian dynamics approach, in which hydrodynamic interactions only couple close neighbours. Rheological properties are strongly influenced by direct intergranular contacts and friction, which are permitted due to a very small distance lubrication cutoff modeling surface asperities. The same critical state as in the dry case is approached in the quasistatic limit. We discuss expressions of rheological laws involving the viscous number instead of the inertial number, and the divergence of effective viscosities in steady flow and in isotropic random suspensions as either the critical state or the random close packing solid fraction are approached

La détermination de la vitesse de propagation des ondes de cisaillement (Vs) dans un sol est fondamentale dans tous les problèmes en dynamique des sols. Plusieurs techniques permettent la mesure de Vs dans les sols. En effet, le paramètre Vs permet la caractérisation et l'évaluation du potentiel de liquéfaction des sols granulaires. Cependant, indépendamment de la technique de mesures utilisées, plusieurs valeurs de Vs peuvent être obtenues. Cette diversification peut se justifier par plusieurs facteurs, notamment la granulométrie. D'où cette étude qui porte sur l'influence de la granulométrie sur la vitesse des ondes de cisaillement des sols granulaires. Dans ce mémoire de maîtrise, l'effet des éléments caractéristiques des sols granulaires (D[indice inférieur 50], e, Cu, [sigma][indice inférieur m]) sur la vitesse des ondes de cisaillement est étudié. Quatre types de sols granulaires ayant des propriétés granulométriques différentes sont considérés. Il s'agit de : (1) Le sable d'Ottawa (C-190) (D[indice inférieur 50] = 0,37 mm et Cu=1,96), (2) Le sable de Milby qui est un sable très fin ((D[indice inférieur 50] = 0,247 mm et Cu= 3,2), (3) Le sable Eastman (D[indice inférieur 50] = 0,6 mm et Cu = 11,05) et (4) Le Till Eastman (D[indice inférieur 50] = 0,1 mm et Cu = 22). Les essais en laboratoire ont été réalisés sur une cellule oedométrique équipée d'un dispositif d'anneau piézoélectrique qui permet de déterminer à la fois la valeur le module oedométique (Eoed) à partir de la cellule et la mesure de Vs à partir de l'essai ring actuator développé à l'université de Sherbrooke. Il ressort de cette étude que la vitesse des ondes de cisaillement augmente avec l'augmentation du D[indice inférieur 50]. De plus, cette étude révèle qu'il existe une corrélation directe entre le module oedométrique et le module de cisaillement Gmax. Cette corrélation peut se traduire par une fonction de la forme de la forme Eoed / Gmax = f([sigma]).

La motivation principale de la présente étude est la compréhension des phénomènes d’agglomération particulaire des systèmes granulaires pharmaceutiques secs dans l’étape de mélange. Une meilleure compréhension des phénomènes d’agglomération est nécessaire puisqu’ils peuvent créer un problème lors de la production pharmaceutique. La traduction de ces connaissances en méthodes peut minimiser les phénomènes d’agglomération non désirés pendant les procédés de production de mélanges particulaires pharmaceutiques secs. Cette étude pourrait contribuer à une meilleure qualité du médicament au patient tout en réduisant les délais et les coûts causés par les mesures correctives en production. Ces travaux de recherche visent à comprendre quelle est la cause fondamentale et quelle est l’influence des différents paramètres sur l’agglomération dans les systèmes granulaires pharmaceutiques secs. L’état de l’art résume les éléments-clés tels que les facteurs intrinsèques, la théorie du mélange particulaire, la ségrégation, la théorie de l’agglomération particulaire, les mécanismes d’agglomération, les forces d’adhésion, les facteurs extrinsèques, les techniques d’analyse de laboratoire et l’analyse multivariée. Plusieurs paramètres pouvant avoir un impact sur l’agglomération lors d’un mélange complexe pharmaceutique ont été étudiés par imagerie dans un tambour rotatif. L’humidité de la poudre favorise l’agglomération, une petite distribution des tailles de particules tend à rendre le mélange cohésif et une valeur élevée de temps de mélange ainsi que de vitesse de rotation du mélangeur tend à désagglomérer. Certaines matières premières ont été identifiées comme problématiques lors de l’ajout dans un mélange. Une forte humidité de 90 %HR du milieu ambiant protège contre la désagglomération. Cependant, tous ces paramètres ont montré des tendances peu significatives. Le temps statique (chargement) est le paramètre dominant ayant eu le plus d’influence sur l’agglomération lors de l’étude du mélange. Étant donné la découverte du fort comportement agglomératif du chlorure de potassium comme résultat, une autre étude paramétrique a permis de comprendre davantage son mécanisme d’agglomération. Les paramètres significatifs de l’étude sont l’humidité de conditionnement, le séchage, le temps de conditionnement et la distribution de taille de particules. La pression exercée à l’échantillon n’est pas un paramètre significatif pour l’agglomération. Le présent mémoire permettra d’enrichir la connaissance dans cette science particulaire et celle de l’agglomération qui sont omniprésentes. Ce mémoire propose deux méthodes pour étudier l’agglomération ; par imagerie en milieu dynamique et par des tests de rupture d’agglomérat en milieu statique. Ce projet démontre le désir et l’initiative canadienne pharmaceutique de poursuivre l’avancement des technologies d’analyse des procédés (PAT) dans les domaines novateurs reliés à l’assurance qualité et à la production. L’industrie pharmaceutique arrivera à mieux maîtriser leur procédé de mélange et à éviter les problèmes engendrés par l’agglomération.

Un bon dimensionnement d’une fondation superficielle nécessite une bonne caractérisation du sol sur lequel elle reposera. Pour le cas des sols granulaires, les essais de reconnaissances destructifs (SPT, CPT, etc.) impliquent la conversion des profils obtenus du sol exploré en un profil de rigidité afin de pouvoir estimer sa portance. Une telle tâche peut induire des erreurs au niveau de l’estimation étant donné qu’elle nécessite le passage d’un nombre de paramètres (N, qc) n’illustrant pas l’état réel du sol (rigidité) à des paramètres caractérisant sa performance (module élastique, coefficient de Poisson). La caractérisation des sols granulaires en utilisant la vitesse des ondes de cisaillement (non destructive) nous permet de contourner ce problème étant donné qu’elle nous donne directement la variation d’un paramètre intrinsèque au sol (Vs) à partir duquel il est possible d’effectuer le dimensionnement. La disponibilité d’une technique de mesure de vitesses de cisaillement (P-RAT) sur des échantillons de sols au laboratoire géotechnique de l’UdeS, nous a permis de mettre en œuvre l’application de Vs dans l’estimation de la portance des sols granulaires. Au cours de ce mémoire, une procédure servant à estimer le tassement des fondations superficielles reposant sur des dépôts granulaires est présenté. Les paramètres de cisaillement (Vs, G) d’un sable déterminés à l’aide des essais de caractérisation (P-RAT) ont été utilisés pour construire un modèle numérique de fondation. La consistance des résultats numériques du modèle a été confirmée à travers des essais de chargement de plaque. Ce travail montre donc l’importance des paramètres de cisaillement dans la modélisation d’une application fréquente en géotechnique. Il constitue un départ considérable pour la construction de relations constitutives basées sur le paramètre vitesse de cisaillement.

This doctoral research aims to deepen knowledge on the impact of electrostatic forces on powder rheology. The initial focus of this thesis was to answer the following question:"Does the presence of electrostatic charges on a pharmaceutical powder carrier have an influence on the spatial lubricant distribution over the carrier particles and thus on the granular rheology?". After the development of an electrification setup and many experiments, the answer to this question turns out to be negative. However, this question led us to discover several phenomena, such as deformation mechanisms of powder columns during shear tests with the FT4 rheometer; and the formation of self-assemblies after the imposition of electrostatic charges on pharmaceutical excipients. These phenomena are explored in this thesis. Therefore, we demonstrate the inability of a FT4 rheometer to conduct shear tests when pharmaceutical formulations are lubricated with magnesium stearate. Indeed,"monolithic cake" behaviour appears. Two numerical models (method of characteristics and analysis of Janssen) are used to assess the stress state and explain this particular case. In addition, a characterization of the shape and particle size distribution of the electrostatic self-assemblies is presented in this thesis. Finally, the previous parts of this thesis showed the importance of wall friction effects during shear tests. Thus, a final part of this thesis focuses on the mechanism of wall lubrication.This latest effort brings us back to the study of the influence of electrostatic forces on powder rheology. Indeed, a high variability with the wall friction results is obtained with unlubricated powders at a relative humidity of 20%. The presence of magnesium stearate and/or higher humidity conditions reduces this variability. However, some tests indicate that electrostatic phenomena are part of the explanation.

La fragmentation des matériaux est un phénomène qui entre en jeu dans de nombreux systèmes naturels et industriels, et à différentes échelles. D'un point de vue industriel, la fragmentation est d'une grande importance dans la production de granulats de carrière, qui sont soumis à des critères stricts en termes de tailles et de formes. Néanmoins, les phénomènes à l'origine de la fracture des matériaux restent encore aujourd'hui mal maîtrisés. Les méthodes aux éléments discrets permettent une modélisation des interactions à l'échelle de la microstructure du matériau. Un des avantages de ces méthodes est que la fissure apparaît et se propage dans l'empilement de manière naturelle, sans qu'il soit nécessaire de la décrire par modèles de fissuration.Un modèle aux éléments discrets permettant de représenter la cohésion au sein d'un matériau composé de particules de formes quelconques a été développé durant cette thèse. Des simulations numériques ont été effectuées pour deux types d'empilements représentant soit des matériaux cimentés, soit des matériaux “pleins”. Les résultats obtenus ont permis de confirmer que ce modèle est capable de reproduire le comportement macroscopique de la rupture de matériau fragile. Le modèle a été ensuite appliqué à l'étude de l'influence du positionnement de points de contact externes sur la résistance d'une particule cylindrique. Une étude expérimentale d'impact a été menée et a permis de confirmer ces résultats. Enfin, nous avons appliqué notre approche à la reconstruction de grains issus de données tomographiques. / The fragmentation of materials is a phenomenon which arises in several natural and industrial systems, and for a wide range of scales. From an industrial point of view, the crushing process is very important in the production of aggregates, which are often required to meet high criteria in terms of size and shape. However, the phenomena behind fracture of materials are still not completely understood. The discrete element methods allow to model interactions at the scale of the material micro-structure, by means of simple models. One of the advantages of this kind of methods is the natural way the crack initiates and propagates in the sample, without any need of a crack model. A discrete element model allowing to describe cohesion inside the material composed of particles of arbitrary shapes has been developed in this work. Numerical simulations have been conducted for two kinds of samples describing both cemented and plain materials. The results obtained have shown the ability of the numerical model to reproduce the macroscopic behavior of the fracture of a brittle material. The developed model has been applied to study the influence of the positioning of external contact points on the effective strength of a cylindrical shaped grain. The results have demonstrated an increase of the required force to break the grain depending on the position of the contacts. An experimental study has confirmed the numerical results. Finally, the model has been applied to the building of grains based on tomographic data.

Cette thèse a pour objet le développement d’un modèle de transport sédimentaire avec la méthode SPH (Smoothed Particle Hydrodynamics). Si les modèles couramment proposés dans la littérature reposent sur un solveur hydrodynamique couplé à des lois semi-empiriques qui modélisent le transport sédimentaire, une approche différente est proposée ici. Dans le modèle proposé dans ce travail, la dynamique du sédiment est également résolue. Celui-ci est assimilé à un milieu continu dont la loi de comportement rend compte de la nature granulaire.Pour ce faire, le modèle élastique-viscoplastique d’Ulrich (2013) a été implémenté dans un code SPH préexistant programmé en langage Cuda, et amélioré du point de vue physico-numérique. Le comportement mécanique du sédiment dépend donc d’une contrainte de rupture déterminée conformément au critère de Drucker-Prager. Dans les zones du matériau où la rupture n’a pas eu lieu, les contraintes de cisaillement sont calculées selon la loi de Hooke généralisée. Dans les zones où la contrainte de rupture a été dépassée, le matériau est assimilé à un fluide rhéofluidifiant. Numériquement, la transition entre les deux états est opérée à l’aide d’une fonction de raccord qui dépend notamment du l’amplitude du taux de déformation et des propriétés granulaires du sédiment.L’eau et le sédiment sont modélisées comme deux phases immiscibles, dans le cadre d’une formulation SPH multi-phase. Pour ce faire, le modèle de Hu et Adams (2006) a été adapté aux modèles de conditions limites semi-analytiques (Ferrand, 2013). Enfin, un schéma d’intégration implicite des forces visqueuses a été développé dans ce contexte, afin d’améliorer les performances du solveur lors de modélisation d’écoulement à bas Reynolds.Plusieurs cas tests sont proposés pour valider le modèle multiphasique, le schéma implicite et le modèle élastique-viscoplastique. De manière générale, les résultats sont en bon accord avec les données expérimentales et analytiques. Le modèle permet de représenter des écoulements multi-fluide avec une bonne précision, même en présence de grand rapport de densité entre les phases. Il en va de même pour les écoulements de fluide non-newtonien et les écoulements à bas Reynolds, pour lesquels le schéma implicite conduit à des résultats très satisfaisants. Enfin, le modèle élastique-viscoplastique a été appliqué à divers cas d’écoulements granulaires, dans le cas d’un matériau sec et saturé, ainsi qu’à des cas d’érosion et d’affouillement. Là encore, les résultats sont globalement en bon accord avec l’expérience / This thesis presents the development and application of a Smoothed Particle Hydrodynamics (SPH) model to bed-load transport. While state of the art simulation methods commonly rely on a fluid dynamics solver coupled to semi-empirical relationships to model the sediment transport, a completely different approach is investigated in this work. The sediment is treated as a continuum whose behaviour law takes account for its granular nature. citepos{ulrich2013smoothed} elastic-viscoplastic model is thus implemented in an in-house code based on the Cuda language, and improved on physical and numerical aspects. The sediment behaviour depends on a yield stress determined according to Drucker-Prager's criterion. In unyielded regions, the shear stresses are calculated in line with the linear elastic theory. In yielded regions, a shear thinning rheological law is used and the transitions between solid and liquid states are ensured by a blending function driven by the strain rate magnitude and sediment granular properties. Water and sediment are modelled as two immiscible phases in the frame of a multi-phase SPH model with semi-analytical wall boundary conditions cite{ferrand2013unified}. An implicit viscous forces integration scheme is also developed to improve the code performance as for low-Reynolds flows.The multi-phase model, as well as the implicit viscous forces integration scheme, were validated on analytical test cases and good agreement was obtained. The multi-phase formulation has also proven its capability to handle flows involving high density ratio, while the implicit viscous forces integration scheme was successfully applied to the simulation of a non-Newtonian flow. The elastic-viscoplastic model was tested on dry and submerged granular flow problems. The model was able to correctly capture the liquid and solid states of the granular material, as well as the failure and the regime transitions. It was also applied to bed-load transport problems for which a good agreement with the experiment was generally found