Yes / Numerical simulation of wet granulation in a continuous granulator is carried out using Discrete Element Method (DEM) to discover the possibility of formation of seeded granules in a continuous process with the aim of reducing number of experimental trials and means of process control. Simple and scooped drum granulators are utilized to attain homogenous seeded granules in which the effects of drum rotational speed, particles surface energy, and particles size ratio are investigated. To reduce the simulation time a scale-up scheme is designed in which a dimensionless number (Cohesion number) is defined based on the work of cohesion and gravitational potential energy of the particles. Also a mathematical/numerical method along with a MATLAB code is developed by which the percentage of surface coverage of each granule is predicted precisely. The results show that use of continuous granulator is promising provided that a high level of shear is considered in the granulator design, e.g. it is observed that using baffles inside the drum granulators is essential for producing seeded granules. It is observed, moreover, that the optimum surface energy for scooped granulator with rotational speed of 50 rpm is 3 J/m2 which is close to the number predicted by Cohesion number. It is also shown that increasing the seed/fine size ratio enhances the seeded granulation both quantitatively (60% increase in seeds surface coverage) and qualitatively (more homogeneous granules).

Asphalt mixture compaction is an important procedure of asphalt mixture construction and can significantly affect the performance of asphalt pavement. Many laboratory compaction methods (or devices), have been developed to study the asphalt mixture compaction. Nevertheless, the whole process from the selection of aggregate to laboratory compaction is still time-consuming and requires significant human and material resources. In order to better understand asphalt mixture compaction, some researchers began to use finite element method (FEM) to study and analyze mixture compaction. However, FEM is a continuum approach and lacks the ability to take into account the slippage and interlocking of aggregates during compaction. Discrete Element Method (DEM) is a discontinuum analysis method, which can simulate the deformation process of joint systems or discrete particle assembly under quasi-static and dynamic condition. Therefore, it can overcome the shortcomings of FEM and is a more effective tool than FEM to simulate asphalt mixture compaction. In this study, an open source 3D DEM code implemented with the C++ programming language was modified and applied to simulate the compaction of hot-mix asphalt (HMA). A viscoelastic contact model was developed in the DEM code and was verified through comparison with well established analytical solutions. The input parameters of the newly developed contact model were obtained through nonlinear regression analysis of dynamic modulus test results. Two commonly used compaction methods (Superpave gyratory compaction and asphalt vibratory compaction) and one linear kneading compaction based on APA machine were simulated using the DEM code, and the DEM compaction models were verified through the comparison between the DEM predicted results and the laboratory measured test results. The air voids distribution within the asphalt specimens was also analyzed by post processing virtual DEM compaction digital specimens and the level of heterogeneity of the air void distribution within the specimens in the vertical and lateral directions was studied. The DEM simulation results in this study were in a relatively good agreement with the experimental data and previous research results, which demonstrates that the DEM is a feasible method to simulate asphalt mixture compaction under different loading conditions and, with further research, it could be a potentially helpful tool for asphalt mix design by reducing the number of physical compactions in the laboratory.

In this thesis, the effects of size and stress on permeability, deformability and strength of fractured rock masses are investigated. A comparison study was carried out to examine the effects of considering, or not considering, the correlation between distributions of fracture apertures and fracture trace lengths on the hydro-mechanical behavior of fractured rocks. The basic concepts used are the fundamental principles of the general theory of elasticity, Representative Elementary Volume (REV), the tensor of equivalent permeability, and the strength criteria of the fractured rocks. Due to the size and stress dependence of the hydro-mechanical properties of rock fractures, the overall effective (or equivalent) hydro-mechanical properties of the fractured rocks are also size and stress-dependent. However, such dependence cannot be readily investigated in laboratory using small samples, and so numerical modeling becomes a necessary tool for estimating their impacts. In this study, a closed-form relation is established for representing the correlation between a truncated lognormal distribution of fracture apertures and a truncated power law distribution of trace lengths, as obtained from field mapping. Furthermore, a new nonlinear algorithm is developed for predicting the relationship between normal stress and normal displacement of fractures, based on the Bandis model and the correlation between aperture and length. A large number of stochastic Discrete Fracture Network (DFN) models of varying sizes were extracted from some generated large-sized parent realizations based on a realistic fracture system description from a site investigation programme at Sellafield, UK, for calculating the REV of hydro-mechanical properties of fractured rocks. Rotated DFN models were also generated and used for evaluation of the distributions of directional permeabilities, such that tensors of equivalent permeability could be established based on stochastically established REVs. The stress-dependence of the permeability and the stress-displacement behaviour were then investigated using models of REV sizes. The Discrete Element Method (DEM) was used for numerical simulation of the fluid flow, deformability properties and mechanical strength behavior of fractured rocks. The results show significant scale-dependency of rock permeability, deformability and strength, and its variation when the correlation between aperture and trace length of fractures are concerned, with the overall permeability and deformability more controlled by dominating fractures with larger apertures and higher transmissivity and deformability, compared with fracture network models having uniform aperture. As the second moment of aperture distribution increases, a fractured rock mass shows more discrete behavior and an REV is established in smaller value of second moment with much larger model size, compared with the models with uniform fracture aperture. When the fracture aperture pattern is more scattered, the overall permeability, Young’s modulus and mechanical strength change significantly. The effect of stress on permeability and fluid flow patterns in fractured rock is significant and can lead to the existence or non-existence of a permeability tensor. Stress changes the fluid flow patterns and can cause significant channeling and the permeability tensor, and REV may be destroyed or re-established at different applied stress conditions. With an increase in the confining stress on the DEM models, the strength is increased. Compared with the Hoek-Brown criterion, the Mohr-Coulomb strength envelope provides a better fit to the results of numerical biaxial compression tests, with significant changes of the strength characteristic parameters occurring when the second moment of the aperture distribution is increased. / QC 20100702

L'évaluation des risques naturels et du risque sismique en particulier a pris une importance croissante pour les autorités publiques ces dernières années, entraînant une adaptation importante des outils utilisés jusqu'à présent dans l'étude de la vulnérabilité du bâti. Il est nécessaire de mieux évaluer ces risques afin de mieux y faire face. De nouvelles approches doivent être proposées pour répondre à ces objectifs dont les approches numériques font partie. Parmi les méthodes numériques développées pour tenir compte de la spécificité du bâti maçonné figurent les approches par éléments discrets comme la méthode NSCD (« Non Smooth Contact Dynamics »). Les travaux présentés dans cette thèse visent à mieux appréhender le comportement mécanique et dynamique des ouvrages maçonnés en utilisant ce type de méthode implémenté dans le code calcul LMGC90, afin d'affiner la compréhension et la caractérisation de la vulnérabilité sismique des édifices en maçonnerie. La méthodologie que nous cherchons à mettre en place vise à traiter la vulnérabilité sismique de bâtiments spécifiques ou au contraire d'ensembles de bâtiments de même typologie. / The perception and the evaluation of natural risks and more particularly of the seismic hazard have never been so important and highlighted by the authorities. This phenomenon leads to a deep adjustment of tools used up to now in the assessment of the vulnerability of buildings. New approaches as numerical approaches have to be proposed so as to best estimate the risk. Among all the numeriacl methods developed to take into account specificities of masonry buildings, approaches by discrete elements methods as the NSCD method (“Non Smooth Contact Dynamics”) seem to be particularly well adapted. This PhD thesis aims to best understand the mechanical behaviour of masonry structures by using this discrete method implemented in the software LMGC90 so as to give an accurate estmation of the seismic vulnerability of buildings. The very purpose of this work is to develop a calculation tool able to deal with the vulnerabilityof both specific and common masonry buildings.

Two new three-dimensional rock mass strength criteria are developed in this dissertation by extending an existing rock mass strength criterion. These criteria incorporate the effects of the intermediate principal stress, minimum principal stress and the anisotropy resulting from these stresses acting on the fracture system. In addition, these criteria have the capability of capturing the anisotropic and scale dependent behavior of the jointed rock mass strength by incorporating the effect of fracture geometry through the fracture tensor components. Another significant feature of the new rock mass strength criterion which has the exponential functions (equation 6.7) is having only four empirical coefficients compared to the existing strength criterion which has five empirical coefficients; if the joint sets have the same isotropic mechanical behavior, the number of the empirical coefficients reduces to two in this new strength criterion (equation 6.10). The new criteria were proposed after analyzing 452 numerical modeling results of the triaxial, polyaxial and biaxial compression tests conducted on the jointed rock blocks having one or two joint sets by the PFC3D software version 5. In this research to have several samples with the same properties a synthetic rock material that is made out of a mixture of gypsum, sand and water was used. In total, 20 joint systems were chosen and joint sets have different dip angles varying from 15 to 60 at an interval of 15 with dip directions of 30 and 75 for the two joint sets. Each joint set also has 3 persistent joints with the joint spacing of 42 mm in a cubic sample of size 160 mm and the joints have the same isotropic mechanical behavior. The confining stress combination values were chosen based on the uniaxial compressive strength (UCS) value of the modeled intact synthetic rock. The minimum principal stress values were chosen as 0, 20, 40 and 60 percent of the UCS. For each minimum principal stress value, the intermediate principal stress value varies starting at the minimum principal stress value and increasing at an interval of 20 percent of the UCS until it is lower than the strength of the sample under the biaxial loading condition with the same minimum principal stress value. The new rock mass failure criteria were developed from the PFC3D modeling data. However, since the joint sets having the dip angle of 60 intersect the top and bottom boundaries of the sample simultaneously, the joint systems with at least one of the joint sets having the dip angle of 60 were removed from the database. Thus, 284 data points from 12 joint systems were used to find the best values of the empirical coefficients for the new rock mass strength criteria. λ, p and q were found to be 0.675, 3.16 and 0.6, respectively, through a conducted grid analysis with a high R2 (coefficient of determination) value of 0.94 for the new criterion given by equation 6.9 and a and b were found to be 0.404 and 0.972, respectively, through a conducted grid analysis with a high R2 value of 0.92 for the new criterion given by equation 6.10. The research results clearly illustrate how increase of the minimum and intermediate principal stresses and decrease of the joint dip angle, increase the jointed rock block strength. This dissertation also illustrates how different confining stress combinations and joint set dip angles result in different jointed rock mass failure modes such as sliding on the joints, failure through the intact rock and a combination of the intact rock and joint failures. To express the new rock mass strength failure criteria, it was necessary to determine the intact rock strengths under the same confining stress combinations mentioned earlier. Therefore, the intact rock was also modeled for all three compression tests and the intact rock strengths were found for 33 different confining stress combinations. Suitability of six major intact rock failure criteria: Mohr-Coulomb, Hoek-Brown, Modified Lade, Modified Wiebols and Cook, Mogi and Drucker-Prager in representing the intact rock strength was examined through fitting them using the aforementioned 33 PFC3D data points. Among these criteria, Modified Lade, Modified Mogi with power function and Modified Wiebols and Cook were found to be the best failure criteria producing lower Root Mean Square Error (RMSE) values of 0.272, 0.301 and 0.307, respectively. Thus, these three failure criteria are recommended for the prediction of the intact rock strength under the polyaxial stress condition. In PFC unlike the other methods, macro mechanical parameters are not directly used in the model and micro mechanical parameter values applicable between the particles should be calibrated using the macro mechanical properties. Accurate calibration is a difficult or challenging task. This dissertation emphasized the importance of studying the effects of all micro parameter values on the macro mechanical properties before one goes through calibration of the micro parameters in PFC modeling. Important effects of two micro parameters, which have received very little attention, the particle size distribution and the cov of the normal and shear strengths, on the macro properties are clearly illustrated before conducting the said calibration. The intact rock macro mechanical parameter values for the Young’s modulus, uniaxial compression strength (UCS), internal friction angle, cohesion and Poisson's ratio were found by performing 3 uniaxial tests, 3 triaxial tests and 5 Brazilian tests on a synthetic material made out of a mixture of gypsum, sand and water and the joint macro mechanical parameter values were found by conducting 4 uniaxial compression tests and 4 direct shear tests on jointed synthetic rocks with a horizontal joint. Then the micro mechanical properties of the Linear Parallel Bond Model (LPMB) and Modified Smooth Joint Contact Model (MSJCM) were calibrated to represent the intact rock and joints respectively, through the specific procedures explained in this research. The similar results obtained between the 2 polyaxial experiments tests of the intact rock and 11 polyaxial experimental tests of the jointed rock blocks having one joint set and the numerical modeling verified the calibrated micro mechanical properties and further modification of these properties was not necessary. This dissertation also proposes a modification to the Smooth Joint Contact Model (SJCM) to overcome the shortcoming of the SJCM to capture the non-linear behavior of the joint closure varying with the joint normal stress. Modified Smooth Joint Contact Model (MSJCM) uses a linear relation between the joint normal stiffness and the normal contact stress to model the non-linear relation between the joint normal deformation and the joint normal stress observed in the compression joint normal stiffness test. A good agreement obtained between the results from the experimental tests and the numerical modeling of the compression joint normal test shows the accuracy of this new model. Moreover, another shortcoming associated with the SJCM application known as the interlocking problem was solved through this research by proposing a new joint contact implementation algorithm called joint sides checking (JSC) approach. The interlocking problem occurs due to a shortcoming of the updating procedure in the PFC software related to the contact conditions of the particles that lie around the intended joint plane during high shear displacements. This problem increases the joint strength and dilation angle and creates unwanted fractures around the intended joint plane.

No description available.

Cette thèse est une étude expérimentale et numérique du comportement au cisaillement de milieu granulaire humide sous l’effet de la quantité de liquide introduite et la contrainte normale appliquée. Les expériences ont été faites sur une cellule de cisaillement annulaire, pour une large gamme de contraintes appliquées allant de presque 0.3 kPa à 12 kPa. Les résultats donnent la variation de la contrainte de cisaillement en régime stationnaire en fonction de la contrainte normale pour une large variation de la quantité de liquide. Le liquide dans le milieu granulaire va de ponts liquides formés au point de contact jusqu’au remplissage totale de l'espace entre les grains. L’effet de liquide sur la résistance au cisaillement et la porosité de milieu granulaire a été analysé. Différents régimes du comportement de milieu granulaire humide ont été identifiés. Afin d’acquérir une compréhension microscopique du comportement au cisaillement de milieu granulaire sec et partiellement humide, la méthode des éléments discrets (DEM) a été utilisée. Des billes de verre de grande taille (2 mm de diamètre) ont été utilisées pour réduire le temps de simulation et faciliter la caractérisation à l’échelle de particule. Une première partie a été consacrée à l’étude de l’effet des propriétés microscopiques de particule (Module de Young et la friction de glissement) sur les propriétés macroscopiques de milieu granulaire sec et humide (le nombre de coordination, la porosité, le ratio de contraintes et la vitesse de particules). Une deuxième partie a été concernée par l’étude du comportement au cisaillement de milieu granulaire humide pour différentes fractions de liquide et différentes contraintes normales appliquées. En particulier, les forces capillaires et le nombre de ponts liquide ont été quantitativement analysés. / We study experimentally and numerically the shear behaviour of wet granular material. We investigate the effect of the liquid content and the applied normal stresses to this behaviour. An annular shear cell was used to carry out the experiments, for a large range of applied normal stress from about 0.3 kPa to 12 kPa. The results give the variation of the shear stress at steady-state as a function of the normal stress for a wide range of liquid fraction. The incorporated liquid goes from forming bridges at the contact point to completely filling the space between grains. The shear resistance and the voidage fraction variations with the liquid fraction were analysed. Depending on the applied normal stress and the liquid fraction, different regimes of the shear resistance were identified. The discrete element method (DEM) was used to gain a microscopic understanding of the shear behaviour of dry and partially wet granular material in the shear cell. Large size glass beads were used to speed up the computational time and to facilitate characterisation at the particle scale. First, the influence of the microscopic properties of the particle (The Young’s modulus and the sliding friction) on the macroscopic properties of dry and wet granular materials (the coordination number, the voidage fraction, the shear ratio and the velocity of particles) was investigated. Secondly, the shear behaviour of the partially wet granular material for different liquid fractions and normal stresses was studied. The capillary forces and the number of liquid bridges were quantitatively analysed.

In recent years, the increasing impacts of landslide hazards on human lives and lifeline facilities worldwide has advanced the necessity to find out both economically acceptable and useful techniques to predict the occurrence and destructive power of landslides. Though many projects exist to attain this goal, the current investigation set out to establish an understanding of the initiation and propagation mechanisms of landslides via numerical simulations, so that mitigation strategies to reduce the long-term losses from landslide hazards can be made. In this research, the Discrete Element Method (DEM) and Computational Fluid Dynamics (CFD) have been used to investigate the mechanical and hydraulic behaviour of granular materials involved in landslides. The main challenge is to provide rational analyses of large scale landslides via small scale numerical simulations. To solve this problem, dimensional analyses have been performed on a simple granular column collapse model. The influence of governing dimensionless groups on the debris runout distance and deposit height has been studied for the terrestrial and submerged granular flows. 3D DEM investigations of granular flows in plane strain conditions have been performed in this research. The input parameters of the DEM model have been calibrated by the numerical triaxial tests, based on which, the relationships between the microscopic variables and the macroscopic soil strength properties are analysed. Using the simple granular column collapse model, the influences of column aspect ratio, characteristic strain, model size ratio and material internal friction angle on the runout distance and deposit height of granular materials have been examined. Additionally, the deformation and energy evolution of dry granular materials are also discussed. The DEM-CFD coupling model has been employed to study the mechanical and hydraulic behaviour of highly mobilized terrestrial / submarine landslides. This model has been validated via numerical simulations of fluid flow through a porous soil sample and grain batch sedimentations. The simulations of granular flows in the submerged environment have led to some meaningful insights into the flow mechanisms, such as the mobilization of sediments, the generation and dissipation of excess pore water pressures and the evolution of effective stresses. Overall, this study shows that the proposed numerical tools are capable of modelling the mechanical and hydraulic behaviour of terrestrial and submarine landslides.

La modélisation multi-échelle multi-physique de la solidification d'alliages métalliques demande de combiner des phénomènes à l'échelle macroscopique du produit et microscopiques à l'échelle des structures de solidification. Dans cette thèse, l'empilement aléatoire des grains équiaxes avec des morphologies typiques de solidification est étudié. Nous mettons tout d'abord en évidence les paramètres hydrodynamiques adimensionnels qui régissent l'empilement de grains équiaxes : le nombre de Stokes, St, le nombre d'Archimède, Ar, et le rapport entre le temps caractéristique de la croissance et le temps caractéristique du mouvement, Γ. Un dispositif expérimental a été conçu par similitude hydrodynamique avec le phénomène réel de l'empilement de la solidification afin d'étudier l'influence de la géométrie des grains équiaxes et l'influence des conditions hydrodynamiques sur la fraction d'empilement. En outre, un outil numérique basé sur le méthode des éléments discrets a été développé pour compléter le travail expérimental de détermination de : la fraction d'empilement locale, le nombre de particules voisines en contact et l'orientation des particules. Des fractions d'empilement entre environ 0,53 et 0,67 ont été mesurées et calculées pour les grains sphériques non-cohésifs, alors que des valeurs allant jusqu'à environ 0,30 sont trouvées pour les grains dendritiques non-cohésifs. Enfin, nous étudions la dynamique de l'empilement, qui est la transition d'un régime de sédimentation à l'équilibre mécanique. L'évolution des variables comme la fraction locale de solide, le nombre de particules voisines en contact et l'orientation du grain en fonction du temps est présentée / Solidification multiphase multiscale modeling of metal alloys is based on the combination of the phenomena at the macroscopic scale of the product and at the microscopic scale of the solidification structures. In this thesis, the random packing of the typical equiaxed grain morphologies in metal alloy solidification is investigated. Firstly, we highlight the hydrodynamic dimensionless parameters governing the grain packing in the melt: the Stokes number, St, the Archimedes number, Ar, and the growth-to-motion ratio, Γ. Subsequently, an experimental setup is designed by hydrodynamic similarity with the actual solidification packing phenomenon in order to investigate the influence of the equiaxed grain geometry and the hydrodynamic conditions on the average solid packing fraction. Additionally, a numerical Discrete Element Method tool is developed to complement the experimental work by accessing to those granular variables which result difficult to be experimentally obtained such as the local packing fraction, the contacting neighbors and the particle orientation. Packing fractions between approximately 0.53 and 0.67 are measured and computed for the spherical noncohesive grains, for different hydrodynamic, frictional and polydispersity conditions, whereas values down to approximately 0.30 are found for noncohesive dendrite envelopes. Finally, we investigate the packing dynamics, which is the transition from a sedimentation regime to the mechanical equilibrium (packing). The evolution of the local solid fraction, contacting neighbors, mechanical contacts and grain orientation are given

Ces travaux de thèse sont consacrés à l’examen du rôle du rayonnement thermique dans le processus de propagation des flammes issues de la combustion des particules d’aluminium dans l’air. Le sujet étant complexe et d’un intérêt industriel, il nécessite de prendre en compte le couplage de nombreux phénomènes physico-chimiques afin de prédire finement les conséquences des explosions de poussières. Une analyse bibliographique approfondie est proposée, concernant les mécanismes d’inflammation et de combustion des particules d’aluminium et aussi concernant les connaissances relatives à la propagation des flammes de poussières. La question spécifique de la nature des échanges thermiques et de l’influence du rayonnement thermique est étudiée. La revue bibliographique souligne les approximations et les hypothèses simplificatrices utilisées dans la littérature permettant donc de définir les pistes d’améliorations. Compte tenu des limitations importantes concernant la physique de ces flammes, un outil de simulation de physique numérique nommé « RADIAN », proche de la simulation numérique directe, a été développé proposant un couplage fin entre les différents modes d’échanges thermiques et la combustion pour modéliser la propagation de la flamme dans un nuage de poussières. La méthode des éléments discrets (MED) est utilisée pour modéliser numériquement les échanges radiatifs entre les particules et les échanges conductifs entre gaz et particules. La méthode des différences finies est utilisée pour modéliser numériquement la conduction thermique dans la phase gazeuse et la combustion. Un modèle radiatif est proposé se basant sur la théorie de Mie sur les interactions rayonnement-particules. Les résultats des simulations sont comparés avec des solutions analytiques et des données expérimentales de la littérature. Mais en plus, une étude expérimentale est aussi conduite afin de mesurer la distribution du flux radiatif devant la flamme et la vitesse de combustion laminaire pour des flammes Méthane-Sic, Méthane-Alumine et Al-air. Un bon accord entre les simulations et les expériences est démontré. La loi de Beer-Lambert relative au transfert radiatif devant le front de flamme s’avère inapplicable et une nouvelle solution analytique est proposée. La présence de particules absorbantes du rayonnement promeut la propagation de la flamme. En particulier, il a été montré expérimentalement et confirmé numériquement que les mélanges riches d’AL-air sont susceptibles d’accélérer rapidement. / In this thesis, the role of thermal radiation in aluminum-air flames propagation is studied. The subject being complex and of industrial interest, it requires the coupling of many physiochemical phenomena to accurately predict the consequences of dust explosions. A thorough literature review is proposed about the ignition and the combustion of aluminum particles and about the available theoretical models of dust flames propagation. The specific question of the nature of thermal exchanges and the influence of thermal radiation is studied. The bibliographic review underlines the simplifying assumptions and hypotheses used in the literature making possible the definition of improvement areas. Because of the limited amount of knowledge available to address these questions, a numerical tool “RADIAN” is developed enabling an accurate coupling between the different modes of heat exchange and combustion. The Discrete Element Method (DEM) is used to numerically model the radiative exchanges between particles and the gas-particle thermal conduction. The Finite Difference method is used to numerically model the thermal conduction through the gas phase and combustion. A radiative model based on Mie theory for radiation-particles interactions is incorporated. The results of the simulations are compared with available analytical solutions and experimental data. An original experimental study is also conducted to measure the distribution of irradiance ahead of the flame front and the laminar burning velocity for methane-air-Sic, methane-air-alumina and Al-air flames. A good agreement between numerical simulations and experiments is demonstrated. The Beer-Lambert’s law for radiative transfer in front of the flame front is found to be inapplicable and a new analytical solution is proposed. The presence of absorbing particles may promote the flame propagation. In particular, it is shown experimentally and confirmed theoretically/numerically that Al-air rich mixtures are likely to rapidly accelerate.