In soil, where trees and crops grow, heavy vehicles shear and compact the soil, leading to reduced plant growth and diminished nutrient recycling. Computer simulations offer the possibility to improve the understanding of these undesired phenomena. In this thesis, soils were modelled as large collections of contacting spherical particles using the Discrete Element Method (DEM) and the physics engine AGX Dynamics, and these entities were analyzed. In the first part of the thesis, soils, which were considered to be continua, were subjected to various controlled deformations and fields for quantities such as stress and strain were visualized using coarse graining (CG). These fields were then compared against analytical solutions. The main goal of the thesis was to evaluate the usefulness, accuracy, and precision of this plotting technique when applied to DEM-soils. The general behaviour of most fields agreed well with analytical or expected behaviour. Moreover, the fields presented valuable information about phenomena in the soils. Relative errors varied from 1.2 to 27 %. The errors were believed to arise chiefly from non-uniform displacement (due to the inherent granularity in the technique), and unintended uneven particle distribution. The most prominent drawback with the technique was found to be the unreliability of the plots near the boundaries. This is significant, since the behaviour of a soil at the surface where it is in contact with e.g. a vehicle tyre is of interest. In the second part of the thesis, a vehicle traversed a soil and fields were visualized using the same technique. Following a limited analysis, it was found that the stress in the soil can be crudely approximated as the stress in a linear elastic solid.

This work aims to develop and implement a linear elastic grain-level micromechanical model based on the discrete element method using bonded contacts and an improved fracture criteria to capture both intergranular and transgranular microcrack initiation and evolution in polycrystalline ceramics materials. Gaining a better understanding of the underlying mechanics and micromechanics of the fracture process of brittle polycrystalline materials will aid in high performance material design. Continuum mechanics approaches cannot accurately simulate the crack propagation during fracture due to the discontinuous nature of the problem. In this work we distinguish between predominately intergranular failure (along the grain boundaries) versus predominately transgranular failure (across the grains) based on grain orientation and microstructural parameters to describe the contact interfaces and present the first approach at fracturing discrete elements. Specifically, the influence of grain boundary strength and stiffness on the fracture behavior of an idealized ceramic material is studied under three different loading conditions: uniaxial compression, brazilian, and four-point bending. Digital representations of the sample microstructures for the test cases are composed of hexagonal, prismatic, honeycomb-packed grains represented by rigid, discrete elements. The principle of virtual work is used to develop a microscale fracture criteria for brittle polycrystalline materials for tensile, shear, torsional and rolling modes of intergranular motion. The interactions between discrete elements within each grain are governed by traction displacement relationships.

Micromechanics of solidluid interactions can play a key role controlling macro-scale engineering behavior of granular media. The main objective of this study is to numerically investigate the micromechanics involved in solidluid mixtures to develop a better understanding of the macroscopic behavior of granular media for different applications. This is accomplished by developing a numerical model coupling the Discrete Element Method (DEM) and the Lattice Boltzmann Method (LBM) and employing it to study three distinct yet interrelated applications throughout the course of this research. In the first application, the DEM model is used to provide a clear relationship between energy dissipated by micro-scale mechanisms versus the traditional engineering definition based on macro-scale (continuum) parameters to develop a better understanding for the frictional behavior of granular media. Macroscopic frictional behavior of granular materials is of great importance for studying several complex problems such as fault slip and landslides. In the second application, the DEM-LBM model is employed for studying the undrained condition of dense granular media. While the majority of previous modeling approaches did not realistically represent non-uniform strain conditions that exist in geomechanical problems, including the LBM in the proposed model offers a realistic approach to simulate the undrained condition since the fluid can locally conserve the system volume. For the third application, the DEM-LBM model is used to study discontinuous shear thickening in a dense solidluid suspension. Shear thickening in a fluid occurs when the viscosity of the fluid increases with increasing applied strain rate. The DEM-LBM results for discontinuous shear thickening were compared to experimental data and proved to be an accurate approach at reproducing this phenomenon. The validated DEM-LBM model is then used to develop a physics-based constitutive model for discontinuous shear thickening-shear thinning in granular medialuid suspension. A closedorm model is then calibrated using the DEM-LBM model and validated against existing experimental test results reported in the literature. Findings of this research demonstrate how micromechanical modeling can be employed to address challenging problems in granular media involving solidluid interaction.

<p>Most of the food, pharmaceutical, and chemical industries rely heavily on the supply of free-flowing powders that finds their application in raw materials, additives, and manufactured products. Improper storage conditions combined with environmental factors affect the free-flowing ability of powders. An undesirable transformation of these free-flowing powders into a coherent mass that resists flow is called caking.  </p> <p>Given the difficulty in quantifying the interparticle forces, both experimentally and numerically, most studies have considered only the humidity effect in powder caking. In this study, the interparticle forces in caked powders were quantified using the simplified Johnson-Kendall-Roberts (JKR) model to account for the material and environmental factors that influence powder caking. The cohesion energy density, which is the ratio of cohesive energy to volume of the particle, was used as the indicator of caking in powders. Simulated force chain network was used to track the relay of interparticle forces under compression. The model was validated experimentally by using caked isomalt powder. The results of the simulations demonstrated that an initial interparticle force of less than 0.01 N did not result in a caked mass. The cohesion energy density was found to be more sensitive to moisture content than consolidation pressure. A 33% increase in moisture at the same consolidation pressure increased the cohesion energy density by 42.45% while a 50% increase in consolidation pressure at the same moisture content increased the cohesion energy density only by 12.23%. </p> <p>In similar, to understand the progression of caking at the bulk level, the development of tensile strength in isomalt with changes in temperature, relative humidity, and consolidation pressures was modeled and validated using the finite element method. In this model, Darcy's equation and species transport equation was used to model the continuity and momentum transfer in porous media. The heat transfer equation was used to solve the energy and the solid bridging model was used to the tensile strength. This study revealed that storing isomalt above 25 ˚C and 85±0.1% RH could initiate caking or increase in tensile strength. An increase in RH from 85% to 86% increased the tensile strength magnitude by 42.7%. Additionally, the study recommends lowering the consolidation pressures during storage to less than 3 kPa.</p> <p>To mitigate caking, a powder flow aid device that could transmit vibration energy to powders through direct contact was developed. The device could be controlled remotely using an android application. The portable flow aid device was then tested under static and dynamic conditions and thereby the evolution of stresses during the operation of the device was mathematically analyzed. The decrease in static angle of repose of isomalt using the developed flow aid device for moisture contents of 3.84, 4.84, and 5.92 % was 45, 42.5, and 33 %. The dynamic analysis revealed that the developed device improved the flow rate of isomalt at 3.82% moisture by about 17.64%. On the other hand, a flow obstruction was observed in isomalt at moisture contents of 4.79% and 5.88%. The device was found to aid the flow of isomalt at 4.79% moisture. These observations were mathematically explained using the stress evolution model which predicted a flow obstruction for isomalt at 4.79 and 5.88% moisture contents.</p> <p> </p> <p><br></p>

In this study, experiments and simulations were conducted to investigate ice and snow. The ice sintering force as a function of temperature, pressing force (contact load), contact duration, and particle size during the primary stage of sintering was formulated using experimental methods along with an approximate, semi-analytic, close-form solution. It was shown that the ice sintering force increases nearly linear with increasing external pressing force but best approximated as a power law for dependency on both contact duration and particle size. Moreover, the exponent of the power law for size dependence is around the value predicted by general sintering theory. The temperature dependence of the sintering force is highly nonlinear and follows the Arrhenius equation. It was observed that at temperatures closer to the melting point, a liquid bridge is observed upon these paration of the contacted ice particles. The ratio of ultimate tensile strength of ice to the axial stress concentration factor in tension is found as an important factor in determining the sintering force, and a value of nearly 1.1 MPa was estimated to best catch the sintering force of ice in different conditions. From the temperature dependency, the activation energy is calculated to be around 41.4 kJ/mol, which is close to the previously reported value. Also, the results for the sintering force suggest that smaller particles are “stickier” than larger particles. Moreover, cavitation and surface cracking is observed during the formation of the ice particles and these can be one of the sources for the variations observed in the measured ice sintering force values. The presence of a capillary bridge in contact between an ice particle and a "smooth" (or rough) Aluminum surface at relative humidity around 50% and temperatures below the melting point was experimentally demonstrated. Experiments were conducted under controlled temperature conditions and the mechanical instability of the bridge upon separation of the ice particle from the Aluminum surface with a constant speed was considered. It was observed that a liquid bridge with a more pronounced volume at temperatures near the melting point is formed. It was showen that the separation distance is proportional to the cube root of the volume of the bridge. The volume of the liquidbridge is used to estimate the thickness of the liquid layer on the ice particle and the estimated value was shown to be within the range reported in the literature. The thickness of the liquid layer decreases from nearly 56 nm at -1.7◦C to 0.2 nm at -12.7◦C. The dependence can be approximated with a power law, proportional to (TM − T)−β, where β &lt; 2.6. We further observe that for a rough surface, the capillary bridge formation in the considered experimental conditions vanishes. The Discrete Element Method (DEM) was employed to simulate the filling behavior of dry snow. Snow as a heterogeneous, hot material which is constituted from spherical ice particles which can form bonds. The bonding behavior of ice particles is important in determining the macroscopic behavior of snow. The bond diameter of ice-ice contacts as a function of time, compressive load, and strain rate is used and a DEM for dry snow was developed and programmed in MATLAB. A beam element with implemented damage model was used in the simulation. The simulated parameters were macroscopic angle of repose, packing density, and surface conditions as a function of temperature and fillingrate. The DEM results were able to verify the existing published experimental data. The simulation results showed that angle of repose of snow decreased with decreasing the temperature, the surface became irregular due to particles rotation and re-arrangement for lower falling speeds of particles, and density increased with depth of deposition.

Granular soils exhibit a generally beneficial change in engineering properties with time after deposition or densification, during a process commonly known as aging. Soil properties reported to change during aging include the small strain modulus and stiffness, penetration resistance, liquefaction resistance, and pile setup. Different hypotheses have been proposed to explain the occurrence of aging in granular soils, including cementation induced by dissolution of silica and precipitation at the particle contacts, cementation due to microbiological activity, internal stress redistribution caused by particle crushing, and buckling of the load chains due to particle slippage. New evidence points out that internal and time-dependent changes in the soil structure caused by particle slippage and rearrangement as the source of the time-dependent variations in soil properties. This investigation is focused on the study of aging at the particle scale to determine its main driving mechanism and define the factors which affect it. Results from an extensive laboratory testing program and computer simulations based on the discrete element method provide insight into the causes of aging and its effects on the macroscopic properties of sands based on the analysis of the interaction between sand grains. / Ph. D.

Optimizing proppant pack conductivity and proppant-transport and -deposition patterns in a hydraulic fracture is of critical importance to sustain effective and economical production of petroleum hydrocarbons. In this research, a numerical modeling approach, combining the discrete element method (DEM) with the lattice Boltzmann (LB) simulation, was developed to provide fundamental insights into the factors regulating the interactions between reservoir depletion, proppant-particle compaction and movement, single-/multiphase flows and non-Darcy flows in a hydraulic fracture, and fracture conductivity evolution from a partial-monolayer proppant concentration to a multilayer proppant concentration. The potential effects of mixed proppants of different sizes and types on the fracture conductivity were also investigated. The simulation results demonstrate that a proppant pack with a smaller diameter coefficient of variation (COV), defined as the ratio of standard deviation of diameter to mean diameter, provides better support to the fracture; the relative permeability of oil was more sensitive to changes in geometry and stress; when effective stress increased continuously, oil relative permeability increased nonmonotonically; the combination of high diameter COV and high effective stress leads to a larger pressure drop and consequently a stronger non-Darcy flow effect. The study of proppant mixtures shows that mixing of similar proppant sizes (mesh-size-20/40) has less influence on the overall fracture conductivity than mixing a very fine mesh size (mesh-size-100); selection of proppant type is more important than proppant size selection when a proppant mixture is used. Increasing larger-size proppant composition in the proppant mixture helps maintain fracture conductivity when the mixture contains lower-strength proppants. These findings have important implications to the optimization of proppant placement, completion design, and well production. In the hydraulic-mechanical rock-proppant system, a fundamental understanding of multiphase flow in the formation rock is critical in achieving sustainable long-term productivity within a reservoir. Specifically, the interactions between the critical dimensionless numbers associated with multiphase flow, including contact angle, viscosity ratio, and capillary number (Ca), were investigated using X-ray micro computed tomography (micro-CT) scanning and LB modeling. The primary novel finding of this study is that the viscosity ratio affects the rate of change of the relative permeability curves for both phases when the contact angle changes continuously. Simulation results also indicate that the change in non-wetting fluid relative permeability was larger when the flow direction was switched from vertical to horizontal, which indicated that there was stronger anisotropy in larger pore networks that were primarily occupied by the non-wetting fluid. This study advances the fundamental understanding of the multiphysics processes associated with multiphase flow in geologic materials and provides insight into upscaling methodologies that account for the influence of pore-scale processes in core- and larger-scale modeling frameworks. During reservoir depletion processes, reservoir formation damage is an issue that will affect the reservoir productivity and various phases in fluid recovery. Invasion of formation fine particles into the proppant pack can affect the proppant pack permeability, leading to potential conductivity loss. The combined DEM-LB numerical framework was used to evaluate the role of proppant particle size heterogeneity (variation in proppant particle diameter) and effective stress on the migration of detached fine particles in a proppant supported fracture. Simulation results demonstrate that a critical fine particle size exists: when a particle diameter is larger or smaller than this size, the deposition rate increases; the transport of smaller fines is dominated by Brownian motion, whereas the migration of larger fines is dominated by interception and gravitational settling; this study also indicates that proppant packs with a more heterogeneous particle-diameter distribution provide better fines control. The findings of this study shed lights on the relationship between changing pore geometries, fluid flow, and fine particle migration through a propped hydraulic fracture during the reservoir depletion process. / Doctor of Philosophy / Hydraulic fracturing stimulation design is required for unconventional hydrocarbon energy (e.g., shale oil and gas) extraction due to the low permeability and complex petrophysical properties of unconventional reservoirs. During hydrocarbon production, fractures close after pumping due to the reduced fluid pressure and increased effective stress in rock formations. In the oil and gas industry, proppant particles, which are granular materials, typically sand, treated sand, or man-made ceramic materials, are pumped along with fracturing fluids to prevent hydraulic fractures from closing during hydrocarbon extraction. In order to relate the geomechanical (effective stress), geometric (pore structure and connectivity), and transport (absolute permeability, relative permeability, and conductivity) properties of a proppant assembly sandwiched in a rock fracture, a geomechanics-fluid mechanics framework using both experiment and simulation methods, was developed to study the interaction and coupling between them. The outcome of this research will advance the fundamental understanding of the coupled, multiphysics processes with respect to hydraulic fracturing and benefit the optimization of proppant placement, completion design, and well production.

According to the Mine Health and Safety Administration (MSHA), between 2006 and 2016, the underground stone mining industry had the highest fatality rate in 4 out of 10 years, compared to any other type of mining in the United States. Additionally, the National Institute for Occupational Safety and Health (NIOSH) stated that structurally controlled instability is a predominant failure mechanism in underground limestone mines. This type of instability occurs when the different discontinuity sets intercept with each other forming rock blocks that displace inwards the tunnel as the excavation takes place, posing a great hazard for miners and overall mine planning. In recent years, Terrestrial laser scanning (TLS) has been used for mapping and characterizing fractures present in a rock mass. TLS is a technology that allows to generate a three-dimensional multimillion point cloud of a scanned area. In addition to this, the advances in computing power throughout the past years, have allowed numerical modeling codes to represent more realistically the behavior of a fractured rock masses. This work presents and implements a methodology that integrates laser scanning technology along with Discrete Element Modeling as tools for characterizing, preventing, and managing structurally controlled instability that may affect large-opening underground mines. The stability of an underground limestone mine that extracts a dipping ore body with a room and pillar (and eventual stoping) mining method is analyzed with this approach. While this methodology is proposed based on a specific case study that does not meet the requirements to be designed with current NIOSH published guidelines, this process proposes a general methodology that can be applied in any mine experiencing similar failure mechanisms, considering site-specific conditions. The aim of this study is to ensure the safety of mine workers and to reduce accidents that arise from ground control issues. The results obtained from this methodology allowed us to generate Probability Density Functions to estimate the probability of rock fall in the excavations. These models were also validated by comparing the numerical model results with those obtained from the laser scans. / M.S. / According to the Mine Health and Safety Administration (MSHA), between 2006 and 2016, the underground stone mining industry had the highest fatality rate in 4 out of 10 years, compared to any other type of mining in the United States. Additionally, the National Institute for Occupational Safety and Health (NIOSH) stated that structurally controlled instability is one of the main causes of rock falls in underground limestone mines. This type of instability occurs when the fractures present in the rock mass intercept each other forming rock blocks that displace into the tunnel as the excavation takes place and poses a great hazard for miners. In recent years, Terrestrial laser scanning (TLS) has been used for mapping and characterizing fractures present in a rock mass. TLS is a technology that allows to generate a three-dimensional multimillion point cloud of a scanned area. In addition to this, the advances in computing power throughout the past years, have allowed simulation softwares such as the Discrete Element Model (DEM) to represent more realistically the behavior of a fractured rock mass under excavation. The aim of this work was to develop and evaluate a methodology that could complement already exisiting design guidelines that may not apply to all kind of underground mines. The presented methodology evaluates rock failure due to presence of discontinuites, through the integration of TLS with DEM and considers site specific conditions. An area of a case study mine was assessed with this methodology, where several laser scans were performed. Information extracted from this laser scans was used to simulate the response of the rock mass under excavation by running Discrete Element Numerical Models. Results from these models allowed us to estimate the probability of rock failure in the analized areas. These, rock block failure probability estimations provide engineers a tool for characterizing, preventing, and managing structurally controlled instability, and ultimately improving workers safety.

Granular matter is defined as a collection of particle grains, such as sand.This type of matter have different characteristics (solid, liquid and gas) depending on the energy level per grain. There are several approaches to modeling and numerical simulations of granular matter. They are used by different groups for different purposes, and the choice between the approaches is based on knowledge and tradition rather than what might be best for the purpose. The key questions are when to use what method and what physical quality is lost depending on the choice.Two regimes of discrete element granular simulations emerge: smooth and non-smooth. To compare the efficiency and physical quality of the two approaches, four physics softwares are examined including Bullet Physics, LMGC90, AgX and LIGGGHTS. Test scenes are setup in each software and the results are compared to each other or to the results of other work.The thesis is performed at UMIT Research Lab at Umeå University.

Yes / Twin-screw granulation has emerged as a key process in powder processing industries and in the pharmaceutical sector to produce granules with controlled properties. This comprehensive review provides an overview of the simulation techniques and approaches that have been employed in the study of twin-screw granulation processes. This review discusses the major aspects of the twin-screw granulation process which include the fundamental principles of twin-screw granulation, equipment design, process parameters, and simulation methodologies. It highlights the importance of operating conditions and formulation designs in powder flow dynamics, mixing behaviour, and particle interactions within the twin-screw granulator for enhancing product quality and process efficiency. Simulation techniques such as the population balance model (PBM), computational fluid dynamics (CFD), the discrete element method (DEM), process modelling software (PMS), and other coupled techniques are critically discussed with a focus on simulating twin-screw granulation processes. This paper examines the challenges and limitations associated with each simulation approach and provides insights into future research directions. Overall, this article serves as a valuable resource for researchers who intend to develop their understanding of twin-screw granulation and provides insights into the various techniques and approaches available for simulating the twin-screw granulation process.