Many powders and particulate solids are stored and handled in large quantities across various industries. These solids often encounter handling and storage difficulties that are caused by the material cohesion. The cohesive strength of a bulk material is a function of its past consolidation stress. For example, high material cohesive strength as a result from high storage stresses in a silo can cause ratholing problems during discharge. Therefore, it is essential to consider the stress-history dependence when evaluating such handling behaviour. In recent years the Discrete Element Method (DEM) has been used extensively to study the complex behaviour of granular materials. Whilst extensive DEM studies have been performed on cohesionless solids, much less work exists on modelling of cohesive solids. The commonly used DEM models to model adhesion such as the JKR, DMT and linear cohesion models have been shown to have difficulty in predicting the stress-history dependent behaviour for cohesive solids. DEM modelling of cohesive solid at individual particle level is very challenging. To apply the model at single particle level accurately would require one to determine the model parameters at particle level and consider the enormous complexity of interfacial interaction. Additionally it is computationally prohibitive to model each and every individual particle and cohesion arising from several different phenomena. In this study an adhesive elasto-plastic contact model for the mesoscopic discrete element method (DEM) with three dimensional non-spherical particles is proposed with the aim of achieving quantitative predictions of cohesive powder flowability. Simulations have been performed for uniaxial consolidation followed by unconfined compression to failure using this model. Additionally, the scaling laws necessary to produce scale independent predictions for cohesionless and cohesive solids was also investigated. The influence of DEM input parameters and model implementation have been explored to study the effect of particle (meso-scale) properties on the bulk behaviour in uniaxial test simulation. The DEM model calibration was achieved using the Edinburgh Powder Tester (EPT) – an extended uniaxial tester to measure flowability of bulk solids. The EPT produced highly repeatable flowability measurements and was shown to be a good candidate for DEM model calibration. The implemented contact model has been shown to be capable of predicting the experimental flow function (unconfined compressive strength versus the prior consolidation stress) for a limestone powder which has been selected as a reference solid in the Europe wide PARDEM research network. Contact plasticity in the model is shown to affect the flowability significantly and is thus essential for producing satisfactory computations of the behaviour of a cohesive granular material. The model predicted a linear relationship between a normalized unconfined compressive strength and the product of coordination number and solid fraction. Significantly, it has been found that contribution of adhesive force to the limiting friction has a significant effect on bulk unconfined strength. Failure to include the adhesive contribution in the calculation of the frictional resistance may lead to under-prediction of unconfined strength and incorrect failure mode. The results provide new insights and propose a micromechanical based measure for characterising the strength and flowability of cohesive granular materials. Scaling of DEM input parameters in a 3D simulation of the loading regimes in a uniaxial test indicated that whilst both normal and tangential contact stiffness (loading, unloading, and load dependent) scales linearly with radius of the particle, the adhesive forces scales with the square of the radius of the particles. This is a first step towards a mesoscopic representation of a cohesive powder that is phenomenological based to produce the key bulk characteristics of a granular solid and the results indicate that it has potential to gain considerable computational advantage for large scale DEM simulations. The contact model parameters explored include particle contact normal loading stiffness, tangential stiffness, and contact friction coefficient. The DEM model implementation parameters included numerical time step, strain rate, and boundary condition. Many useful observations have been made with significant implications for the relative importance of the DEM input parameters. Finally the calibration procedure was applied to a spray dried detergent powder and the simulation results are compared to whole spectrum of loading regime in a uniaxial experiment. The experimental and simulation results were found to be in reasonable agreement for the flow function and compression behaviour.

Dans le cadre du projet ASTRID, le procédé de fabrication du combustible nucléaire par métallurgie des poudres, pour les RNR, est revisité dans le but d’être simplifié. En particulier, nous cherchons à supprimer l’étape de granulation mécanique de la poudre qui lui confère un bon comportement à l’écoulement. La problématique, est d’arriver à faire s’écouler spontanément et rapidement une poudre, à travers un orifice par lequel elle ne s’écoule pas naturellement. Par ailleurs, la poudre alterne entre des phases d’écoulement et de non-écoulement.Pour garantir l’écoulement de la poudre à travers le moule de presse, nous appliquons des vibrations horizontales. Les vibrations permettent de faire s’écouler la poudre et d’atteindre le débit requis par les cadences de fabrication. Toutefois, elles ont également le défaut, lors des phases de non-écoulement, de compacter la poudre et de retarder les futurs écoulements. Tout l’art du remplissage des moules de presse assisté par vibrations consiste à maîtriser le caractère ambivalent des vibrations.L’évolution singulière de la compacité des poudres d’actinides, lors de la phase de non écoulement, nous a amené à définir un modèle stochastique unidimensionnel de compaction. La confrontation des résultats de simulation, avec les lois empiriques de compaction, a permis d’identifier le sens physique des paramètres d’ajustement des lois empiriques. Nous avons aussi proposé une nouvelle loi d’évolution de la compacité à deux exponentielles étirées. Cette nouvelle loi rend compte, non seulement de la cinétique de compaction des poudres d’actinides mais également de l’ensemble des résultats que nous avons trouvé dans la littérature. / In the framework of the ASTRID project, the nuclear fuel production process by powder metallurgy, for Fast Neutron Reactors, is revisited in order to be simplified. In particular, we seek to remove the mechanical granulation step of the powder which gives a good flow behavior during the filling of press molds. The aim is to reach a spontaneous and quick powder flow through a hole in which the powder does not flow without external energy. Furthermore, the powder alternates between flow phases during the filling of press molds, and non-flow phases during the compaction and ejection of the pellet.We hence apply horizontal vibrations to ensure the flow of the powder through the press mold. The vibrations help the powder to flow and increase the production rates. However, they have the disadvantage to compact the powder and delay the future flows, during the non-flow phases. The art of filling the press mold assisted by vibrations is to master / control the ambivalent nature of the vibrations.The remarkable packing fraction evolution of actinides powders, during the non-flow phases, allows us to define a simple 1D stochastic model to understand the compaction kinetics. The comparison of the stochastic model with the empirical compaction laws found in the literature helps us to identify the physical meaning of fitting parameters proposed by the empirical models. Furthermore, we have also proposed a new compaction law with two-stretched exponentials. This new law not only reflects the compaction kinetics of actinides powders, but also of all the compaction data we found in the literature.