Search Documents

By C. Eswaraiah

A drop weight tester (DWT) is designed and fabricated for the purpose of analyzing single particle impact breakage characteristics of different ores. The test results were evaluated through the breakage distributions of different size fractions at five different input energy levels. The breakage parameter, t10 (Narayanan, 1986) is used to represent the degree of size reduction, which is assumed to be representative of the breakage product size distribution obtained from single particle breakage. Representation of breakage distribution with this method consist one more breakage dependent parameter (a). To evaluate the impact breakage characteristic of different ores, a mathematical relationship was established to analyze the dependency of breakage parameters, t10 and a on specific input energy (EIS) or fracture energy (Tavares et al, 1998, Fandrich, 1998), specimen size (?) and Bond?s operating work index (WI0). Mathematical results developed gives adequate model for all ball mill modeling. It can be concluded that the drop weight test method is a useful and practical way of evaluating the impact strengths of various ores and results of which can be used in modeling mills.

Sep 1, 2012

With the ambitious target of adding another 2,00,000 MW of power generation by the year 2017 in India the demand for coal is bound to grow. The non-coking coal deposits in India are characterized as high ash and difficult-to-wash coals. Quality improvement by beneficiation is almost mandatory due to MoEF notification to use <34% ash coal for all power plants except those located at pithead. Most of the coal washeries in India are based on wet process, where the coal fines are allowed to settle in thickeners. The coarser particles settle faster leaving the fine particles suspended for longer period. The dispersed fine coal particles contain some surface charge, which leads to the development of electrical repulsion between the particles, and as a result of this, the suspension becomes more stable. In order to destabilize the colloidal suspension, flocculants are usually added, so that the settling rate of the fine particles also increases. These flocculants are of cationic, anionic and non-anionic depending upon the type of charge released when dispersed in water. The flocculants may act partly for charge neutralization and partly for inter particle bridging and electrostatic patch formation. The effectiveness of any type of flocculants depends upon the molecular weight and charge density of the polymer. The paper discusses the settling behavior of non-coking coal fines suspension by using different flocculants. Keywords: floculent, coal fines, sedimentation rate, effluent

Sep 1, 2012

By D. Kosior

Influence of frothers (a-terpineol and n-octanol) concentration on kinetics of a three phase contact (TPC) formation by a bubble colliding with a hydrophobic surface (Teflon) of different roughness was studied using a high-speed camera (1040 frames per second). It was found that the Teflon surface roughness and the frother concentration strongly affected time of the TPC formation (tTPC) and time of drainage (tD) of the liquid film separating the colliding bubble from the solid surfaces. With increasing Teflon surface roughness, from ca 1 to 80 µm, the time of the TPC formation was shortened by over the order of magnitude, from ca 100 to a few milliseconds only. At low a-terpineol and n-octanol concentrations, thetTPC values were shorter than in distilled water due to lowering the bubble velocity in the frother presence and diminished bouncing of the colliding bubble. However, at high frothers concentrations the tTPC values were prolonged, despite that the bubble velocity was similar. The prolongation of the tTPC at high frother concentrations indicates that air (micro- and/or nano-bubbles) was present at the Teflon surfaces of different roughness. The mechanism of the tTPCprolongation at high frother dosage is described.

Sep 1, 2012

By Mohammad Kor, David Beattie, Agnieszka Mierczynksa-Vasilev

The selectivity of polymer adsorption in mixed mineral systems is the basis for polymer depression of gangue minerals in flotation. Effective polymer depressants not only reduce gangue mineral recovery, they also have minimal effect on the recovery of the value phase, presumably through targeted adsorption on the gangue phase, and minimal uptake on the value phase. Traditionally, polymer depressant design has been guided by single mineral phase adsorption tests, with affinity and maximum adsorbed amount determined for gangue and value phases separately. However, there is a need to determine the actual extent of adsorption on the gangue phase while in the presence of the value phase. In this work, we report on the development of in situ atomic force microscopy (AFM) imaging for determination of polymer depressant coverage and morphology on mixed mineral surfaces. Our group has previously used AFM imaging on single minerals to connect adsorbed layer characteristics with bubble-particle attachment and flotation outcomes. Two sample types were used to test the methodology: (i) a mixture of two crushed/ground minerals encased in resin; and (ii) a composite mineral sample containing macroscopic regions of three different mineral phases. The samples have been used to probe the selective uptake of two modified dextrin polymers (octenyl succinate anhydride dextrin). AFM images acquired in situ (while all surfaces are immersed in the same polymer solution) of each mineral phase have revealed the natural selectivity of the polymer for the mineral surfaces. In all cases, the polymers adsorb on the all mineral surfaces present, but to a varying extent.

Jan 1, 2014

By G. Bulut, Ü. Yenial

"Heavy metals contamination in aqueous systems can be found in a variety of industries, in particularly mining, metal processing, finishing and plating. Water containing more than 0.5 ppm heavy metal considered as a toxic due to harmful effects on health. The major techniques employed for heavy metal removal from industrial wastewaters include precipitation, reduction, ion-exchange, and adsorption. Adsorption is a heavy metal removal method for aqueous systems. As a common adsorbent, activated carbon has been used in treatment but it is ineffective on some elements. Thus, there has been a research for developing highly effective and selective adsorbents. In this study the adsorption of lead and arsenic ions on to coal waste and red mud waste was studied. The experimental parameters were determined as pH, equilibrium time and the initial concentrations, activation by heat treatment and the blend of adsorbents. According to their adsorption capabilities, red mud was found to adsorb lead and arsenic ion with high efficiency but coal waste did not adsorb arsenic ion. However the binary adsorptions of these ions were enhanced the adsorption efficiencies. In addition the blend of coal waste and red mud affected the adsorption efficiency. INTRODUCTION As a result of rapidly developing industry and excessive use of water resources, the need of clean water and wastewater management appear as a problem. Wastewater recycling activities and researches come to the fore for an effective use of wastewater. Recycling of wastewater, which contains heavy metals harmful to the environment and human, is a difficult process. Heavy metals are generally divided into three categories: toxic metals, precious metals, radionuclides. Wastewaters from various industries, containing heavy metals, are the main polluters of water and soil. Although their toxic structures, heavy metals are used extensively in industry as they have technological significance. Usually mining industry, metal industry and industrial organizations produce heavy metal containing wastewaters (Volesky 1990, Ting and Lawson 1991, Bishop 2002)."

Jan 1, 2018

By Thenepalli Thriveni, Seong-Young Nam, Nam-Il Um

New modern water treatment technologies required for arsenic removal from waste waters contaminated extensively by various background sources from mining, industrial waste, agricultural waste, painting industry and natural geological bed rock system etc. Arsenic is a widely dispersed element in the earth&apos;s crust and exists at an average concentration of approximately 5 mg/kg. Arsenic (As) and heavy metals may be released from the mine wastes to the ground and surface water systems, as well as the geological environment due to their solubility and mobility. This research focuses on the effect of accelerated carbonation process for successful removal of heavy metal arsenic anions from artificial waste water and mining waste waters. This accelerated carbonation process is a simple and novel method with 85 ~ 88% removal efficiency. Carbonation has been studied as a process to further improve immobilization of solidified/stabilized waste and granular waste. After successful removal of arsenic anions from solid mining waste and its potential benefits for producing construction materials. Besides, the industrial and urban management systems are generating bulk solid mining wastes, and most often dumping them in open fields. To prevent the landfilling with mining waste, many efforts are being made for the recycling of different types of mining solid wastes with a view to utilizing them in the production of various construction materials.

Jan 1, 2014

By Lucy van de Ruit, Megan Becker, Aubrey Mainza, Jenny Wiese

"UG2 chromitite ore is a complex, highly disseminated platinum ore from the Bushveld Complex in South Africa, with typical platinum group mineral (PGM) grain sizes below 15µm. The introduction of stirred mills to the platinum industry has enabled grinding to finer product sizes, which has improved liberation and recovery of the finely grained PGMs. However, reduction in particle size leads to an increase in chromite entrainment in the flotation concentrate, which has a negative impact on the subsequent smelting process. Factors such as particle size, froth structure, solids concentration and water recovery have all been identified as impacting entrainment, but whether particle shape also has a significant effect is unclear. The purpose of this study was to determine whether the recovery of chromite particles by entrainment during the flotation of UG2 ore is affected by particle morphology. A laboratory vertical stirred mill with inert ceramic media was used for grinding UG2 ore which was then floated in a batch flotation cell. Flotation feed and concentrate samples were analyzed using QEMSCAN to quantify the angularity and aspect ratio distributions of the particles. The stirred mill product was found to have relatively broad aspect ratio and angularity distributions, with the major mineral components (chromite and silicates) having very similar shape characteristics. Analysis of flotation feed and concentrate samples suggested that particle shape does not have a significant effect on chromite entrainment."

Jan 1, 2014

By N. Schreithofer

Varying flow regimes, particle populations and consequently different types of bubble-particle interactions characterize the different zones of a flotation cell. The understanding of the role of the different zones of turbulent intensity on the particle-bubble interactions is of primary importance to optimization and scale-up of flotation processes and machinery. In order to study particle-bubble interaction in turbulent flow a completely new approach was taken. A laboratory scale stationary turbulence device was developed and built. The set-up allows the observation of particle-bubble interactions in various turbulence intensities, particle and bubble sizes and surfactant concentrations. The experimental set-up consists of a transparent, specially designed tank in which liquid flows between carefully designed and placed obstacles that create different turbulent flow fields. Bubble-particle aggregates are introduced into this turbulent zone. The behaviour of the aggregates, such as attachment /detachment of particles to/from air bubbles is recorded by a High Speed Imaging (HSI) system (5400 fps). The images are then analysed and the dependence of attachment/detachment process is described as a function of turbulent intensity in different conditions for different surfactant concentrations. This paper introduces this method to study particle-bubble interactions in turbulent flow and correlation to Computational Fluid Dynamics (CFD) simulation. CFD modelling is validated for two and three phase systems; liquid-gas with liquid and solid tracers respectively. The CFD simulations using the RANS k-S turbulence model showed good correlation to the experimental data and the experimental method is a promising tool for investigation of particle-bubble interaction in turbulent flow.

Jan 1, 2014

By P. K. Sen, Mrinal K. Ghose

Tailing slurry is usually abrasive and possesses high viscosity. This paper examines the existing iron ore beneficiation process and focuses on some of the key environmental issues related to the disposal of tailings. It identifies that almost 32% of the iron ore extracted ends up as tailings and they are usually disposed in slurry form in tailing ponds. By the year 2000, 27 Mt of iron ore tailings were generated in India, which has become one of the biggest waste handling problems to the industry. Characteristics of the tailings vary widely. So the physical nature of tailings cannot be considered separately in isolation without considering the mill effluent from the beneficiation plant. In the present investigation, studies were carried out in three iron ore mines in India to evaluate the characteristics of tailings from the beneficiation plant. The tailing samples were collected from different locations during a time period of one year and analyzed. The results will be discussed and the characteristics of the tailings have been evaluated. The disposal system of tailings at the three mines was evaluated. The research on the settling characteristics of tailings indicates that the area provided in the existing tailing ponds in India is much more than that of the experimental results. This paper concludes that sedimentation tests are essential to design appropriate tailings pond with respect to size. In the proposed scheme it is suggested that an additional unit process, i.e. sedimentation tank be added to the process before the disposal of tailings to the tailing pond.

Jan 1, 2003

By Michael Young, Rodrigo Araya, Cesar Gomez, Le Huynh, James Finch

The performance of flotation machines is closely related to the gas dispersion characteristics achieved in their operation. It is known that frothers play a major role not only in the size distribution of the generated bubble population, but also on the swarm rise velocity. A program to characterize frother effects on two-phase hydrodynamics (air-water) of Jameson cells was undertaken in a fully instrumented pilot unit installed at McGill University. Gas dispersion parameters (gas velocity, gas holdup and bubble size distributions) of the gas-liquid mixture generated in the downcomer were measured after its release into the separation chamber. The unit was operated under different conditions including gas velocity, liquid jet velocity and frother concentration (MIBC). Bubble size distributions for concentrations above the critical coalescence concentration (CCC) were narrow with Sauter mean diameter between 0.6 to 1.3 mm. The bubble size followed the general trend as encountered in other types of flotation machines, i.e., rapid decrease with increasing frother concentrations below 10 µL/L, and almost no change in size for concentrations above 30 µL/L. The CCC was between 10 and 15 µL/L. Strong correlations were found between bubble size and air-to-liquid flow ratio (APR), and between the bulk density of the aerated solution in the separation tank and the vacuum inside the downcomer. An equation is proposed to estimate gas holdup in the separation tank from vacuum and bulk density readings.

Jan 1, 2014

By Bern Klein, Givemore Sakuhuni

Flotation is the dominant base metal extraction technique and it requires ore to be ground below the liberation size for effective recovery due to its dependence on mineral surface properties. For a poly-metallic ore, differential hardness of constituent minerals leads to over grinding of softer minerals before the harder minerals are fully liberated, which makes selection of a target grind size complex, and the complexity increases when treating leaner ores. Flotation recovery deteriorates above 100 and below 10 microns causing significant loses of valuable minerals in most flotation circuits. Although enhanced gravity concentration is known to be effective at recovering particles above 100 microns, which are difficult to hydro-dynamically transport, exploiting differential specific gravity for base metal separation was not an option in the past. This was partly due to the absence of high mass yield gravity techniques. The advent of continuous gravity concentrators with flexible mass yields, in the past two decades, has necessitated investigation of possible base metal applications. The Knelson CVD, one of the continuous discharge centrifugal concentrators, has an established application niche for sulphides associated with gold. However, its base metal applications have never been rigorously tested. Considering the dependence of flotation recovery on particle size, continuous discharge centrifugal concentration can be integrated in flotation circuits to capture flotation non-recoverable coarse particles. This paper explores how particle size, liberation and mineral hardness influence application of continuous centrifugal technology to regulate galena deportment in a polymetallic base metal plant. The results show that particle size affects galena recovery more than liberation. In view of these results, the CVD can be integrated in the primary grinding circuit to pre-concentrate the heavy sulphides at a coarser grind size with potential grinding energy savings.

Jan 1, 2014

By K. D. Mehta

Copper concentrate containing complex sulphides of Cu-Ni with 14.13% Cu and 10.04% Ni besides impurities such as 0.33% Co and 0.71% Mo, is not acceptable to the copper smelters due to the smelting problems. Selective bioleaching of nickel from the concentrate can be considered an an option to make it suitable for copper smelter. Effect of elemental sulphur and pyrite on the bioleaching of Ni and Co has been investigated using a mixed culture of Acidithiobacillus ferrooxidans (A.f.) and Acidithiobacillus thiooxidans (A.t.). Parameters such as pH, temperature, particles, size etc. were optimised to understand the role of additives in the bioleaching. With 2% S at 308K, 2.0 pH, 5% pulp density and 200µm size particles while shaking at 100 rpm, 79% Ni, 41% Co and 5% Cu were bioleached as compared to the leaching of 81% Ni and 45% Co in presence of 10 wt% pyrite in 60 days; 27% Ni, 24% Co and 4% Cu were recovered in sterile experiment. The selective dissolution of nickel and cobalt as compared to copper may be explained by galvanic interaction between chalcopyrite: pentlandite / pyrrhotite. Bioleaching was facilitated by strong oxidizing conditions due to higher Fe3+/Fe2+ ratio/redox potential (ESCE) as a result of improved interaction of minerals through biooxidation of sulphur and accelerated galvanic interaction by pyrite. The mechanism of selective leaching of Ni and Co was established by XRD phase identification.

Sep 1, 2012

By M. McClintock

"High Pressure Grinding Rolls (HPGR) is increasingly being chosen for comminution circuits for its energy-efficiency for hard rock mining operations. Despite HPGR’s growing popularity, geo-metallurgy programs are still prohibitively expensive because of large sample requirements. Traditionally, sizing HPGR application requires a significant amount of sample, in the order of 1 to 10 tonnes, depending on the level of study and complexity of the resource. Therefore, engineering and design of HPGR circuits typically are reliant on few composites that are required to be representative of the resource and mine plan. Three bench scale HPGR test methodologies were developed in 2015 at NBK Institute of Mining Engineering at the University of British Columbia (UBC) capable of defining the energy size breakage relationship of pilot HPGR. These methodologies were developed to allow bench scale testing for different levels of study from PEA (Preliminary Economic Assessment) to detailed design and operations. The Piston Press test defines the energy size reduction relationship by measuring specific energy consumption and corresponding feed and product size distributions. The advantage of Piston Press testing is the small sample requirement of 5 to 10 kg. Small sample requirements allow multiple samples to be tested across a deposit. Piston Press testing allows extensive geo-metallurgy programs that can determine ore variability and inform in optimization and production planning. To date, the Piston Press methodologies have been validated on a pilot-scale HPGR unit located at UBC. Recently, testing was conducted to calibrate and validate the Piston Press test against full-scale HPGR on a gold sulphide ore deposit. The program included the collection of operational data with the respective HPGR feed and recycle sample. Sample collected was used both for Piston Press testing and HPGR piloting. This paper discusses the test-work and results as carried out during this program, using the aforementioned methodologies."

Jan 1, 2018

By K. Price, G. Wang, L. R. Moore, Y. R. Xiong, J. Gu

"Sand is a naturally occurring granular material that is a result of the physical or chemical breakdown of rocks. Raw/natural sand contains various minerals, and it is this combination of minerals that brings the value to the construction aggregate industry. However, it is the quartz component of sand that is valued in its application as industrial silica in glass, high purity quartz, foundry sand, and oil and gas fracturing. Depending on the industry, the grade requirements can be the concentration of SiO2, particle size or shape, as well as contaminants such as the heavy refractory minerals, clays, fines, and alkali metals. With the ever declining mineral feed grade, processors continue to face challenges making the industry specific concentrate grade. As such, it becomes more important to tailor make the processing reagents to fit the mineral, processing conditions, and desired quality specification. This paper focuses on the application of novel collectors towards the beneficiation of sand for each of the major industrial silica applications. INTRODUCTIONQuartz is a mineral composed of a framework of silicon and oxygen atoms. It is the second most abundant mineral in the Earth’s continental crust. There are many species of quartz, which are defined by the particle size, shape, and color. Silica sand is one of those species, which is the byproduct of the physical or chemical breakdown of quartz containing rocks or crystals. Geologist define silica sand as a silica mineral with a particle size between 0.05 and 2 mm in diameter (Zdunczyk, 2001). Silica sand has significant value to society in recreation (beach sand and toy boxes), construction aggregate (concretes and asphalt), industry (iron & steel casting, fracturing proppants, and glass), as well as specialty products such as semiconductors or silicon wafers used in electronics.The aggregate business is the dominant consumer of silica sand products. Conversely, glass making, foundries (casting material), abrasives, and fracturing proppants are the dominant applications for industrial silica sand. Not much beneficiation is required for sand products going into the construction industry, other than sizing and scrubbing to remove clay. However, the industrial silica sand applications are dependent on particle size, shape, and degree of contaminants, which sets the market value (table 1). Therefore, as the more specialty products are produced, more sophisticated techniques and chemistry need to be applied to ensure the appropriate quality is met."

Jan 1, 2018

By O. Guven, I. Karaagaçlioglu, O. Ozdemir, M. Çelik

In this study, the effect of morphological properties such as shape factors and roughness of the glass beads on their flotation efficiency were determined using image analysis and profilometer methods. There is a considerable correlation between the shape factors, surface roughness and flotation behavior of the glass sample. These results indicate that as the particle shape is away from roundness and smoothness, it acquires more hydrophobicity.

Jan 1, 2014

By Ivan Bonchev, Elisaveta Djurova, Evgeni Visariev, Viktor Stoilov

"Two types of slag cooling is used world wide - cooling in ladles and cooling in pits. The flotation comparison of the type cooling regimes was done at Aurubis Bulgaria. In Copper smelter Aurubis Bulgaria (Pirdop), the residual copper in the slags from Flash furnace and Converters is recovered by flotation. The slag flotation plant processes blend of three parts FSF and one part converter slag. The slag processing comprises: pit cooling, crushing, two stages grinding (autogenous and ball) and flotation. The copper content in flotation concentrate is over 20 % and in flotation tail is under 0.5 %. In order to improve the copper recovery, different cooling rate of the liquid slags should be achieved. Three cooling experiments with flash smelting furnace (FSF) slag have been done: two of them in pits with different volume and one cooling experiment in metallic ladle. The ambient air temperature during the cooling tests was from -4 to +4 OC. Laboratory flotation tests in open cycle were done with an aim to compare the slags behaviors cooled at different conditions. In this paper the effect of different cooling on the flotation process will be presented. The copper recovery and distribution by sizes in the tails, copper content in the rougher concentrate, and the concentrate yield are discussed."

Jan 1, 2014

By Massimiliano Zanin, William Skinner, Keith Quast

With the scarcity of fresh potable water on mining sites, particularly those in arid and remote regions, the use of sea water is being contemplated and actually used at a number of mineral concentrators worldwide. The use of sea water, with its associated high dissolved solids content and ionic strength will affect all aspects of processing (e.g. frothing and froth stability, buffering capacity for pH control, solubility of minerals and corrosion of media and equipment). Selective depression of pyrite is more difficult in sea water due to the increased lime addition to achieve a high pH. This results in a reduction in molybdenite recovery. The use of grinding media with various chromium contents and its effect on subsequent Cu and Mo flotation has been investigated in the current study, with both grinding and flotation testing conducted in artificial sea water. Pulp chemistry (pH, Eh and DO) have been logged during grinding, and flotation testing was conducted on subsequent slurries. The flotation results show that the system is remarkably robust and able to cope with wide ranges of pulp chemistry parameters. The effects of artificial sea water on the surfaces of minerals during grinding have been examined using X-Ray Photoelectron Spectroscopy

Jan 1, 2014

By M. Wierink

The lifecycle of a resource, specifically a metal, is often depicted as circle from primary production from mining to minerals processing and via manufacturing and use to disposal and recycling. When a material is recycled, it is considered as saved energy compared to producing the material from primary source. This is, however, a simplified representation. At each step of the life cycle, there are changes in the quality of the material due to different materials chemical and physical interaction, phase transformations due to diffusion, radiation etc. Due to these physical and chemical phenomena there are losses of material, which as determined by thermodynamics are not recoverable within the life cycle. This can be considered as an increase embodied energy per mass unit of substance and loss exergy from the material cycle. With low cyclic loss, the material will remain in the recycling loop for a long time, while with high cyclic loss it will disappear fast. Where you put the material determines how fast it disappears from the resource cycle. In this research paper the phenomena of quality degradation of material and idea of increase of embodied energy and exergy loss within a resource cycles is portrayed. The cyclic nature of material states and transitions and energy consequences of material quality degradation as well as the parameters affecting the length of the time scales of the material cycles are described. Currently the cyclic embedded energy and the history of material is not taken under consideration in generic life cycle assessment. This easily results into a conclusion that lifetime of product is optimized only in terms of how much energy an appliance uses during its lifetime, while ignoring the loss of embedded energy due to the effects caused by the environs the material is subjected to during the cycle. Better description of the time properties of material give us a possibility to design toward optimum lifetime of products in terms of resource efficiency. Keywords: product lifetime, material lifetime, recycling, cyclic resource efficiency, material state, material quality, exergy, embodied energy

Sep 1, 2012

By Gaurav Soni, Roe-Hoan Yoon

A first principle flotation model has been derived from the basic mechanisms involved in flotation. The model consists of a series of analytical equations for bubble generation, bubble-particle collision, attachment, detachment, and froth phase recovery. The model can predict flotation from both hydrodynamic and surface chemistry parameters such as bubble size, particle size, energy dissipation rate, ?-potentials, contact angles (?), etc. The model predictions have been compared with the results of a series of laboratory flotation tests conducted using mono-sized glass beads of different particle sizes and contact angles under different operating conditions. Some of the experiments were conducted using a 1:9 mixtures if hydrophobic silica and hydrophilic magnetite, so that the product graded could be readily determined by simple magnetic separation. The model predictions are in good agreement with the experimental results both in recovery and grade. A computer simulator developed on the basis of the model is capable of designing a preliminary flotation circuit on the basis of feed characteristics, e.g., size-by-liberation data, and the hydrodynamic information on the flotation cells to be incorporated in the circuit.

Jan 1, 2014

By S. Banisi, Seyyed Omid Mousavi, Alireza Ghasemi

Discrete Element Method (DEM) is extensively used to simulate the behaviour of particles in various processing units. This method is based on modelling the forces acting between particles in any contact and consequently calculating the new position of particles. The high number of elements and numerous equations needed have made the required computation time very long even using computers with very fast processors. The required computation time mainly depends on the time step. If the time step is chosen very short, the computation time will significantly increase. On the other hand, if the time step is chosen very long, the simulation will be in error due to not fully observing the contacts. In DEM calculations, the time step is chosen as a fraction of the collision time. In this research, a relationship was proposed between the contact time fraction and the error of simulation. In other words, to select the time step in addition to physical parameters, the accuracy of the simulation was also accommodated in the frequently-used relationships. The simulation error was defined as the relative difference between the initial restitution coefficient and the restitution coefficient obtained after simulation. The required time steps to arrive at the simulation error of 5% were calculated for two common contact force models namely, Hertz-Mindlin and linear spring-dashpot. The simulation was performed for two particles with radius of 3 cm, elasticity modulus of 210 GPa, Poisson&apos;s ratio of 0.3 and relative velocity of 0.5 m/s. The required time steps were found to be 2.3 and 0.19 µs for the Hertz-Mindlin and the linear spring-dashpot contact force models, respectively. Results showed that with a scaling-down the modulus from 210 GPa to 2.1 MPa, the required time steps for the Hertz-Mindlin and the linear spring-dashpot contact force models, with an equal simulation error, increased by 100 and 316 times, respectively.

Jan 1, 2014