Search Documents

By Janusz S. Laskowski

It is customary to assume that wettability plays an important role in coal flotation, and it is well established that while hydrophobic metallurgical coals float easily, low-rank less hydrophobic coals are difficult to float. However, coal surface wettability plays a dominant role in many coal preparation unit operations. As it will be discussed in this paper, it affects coal-dewatering, selection of flocculants, fine coal pelletization and briquetting, and even coal handleability.

Sep 1, 2012

By L. Stamenov, H. Bouzahzah, S. Gaydardzhiev, D. Naumov

"The Chelopech epithermal arsenic bearing high-sulphidation, copper-gold deposit located in Bulgaria is often referred as the richest European gold deposit currently under exploitation. Key points from the flotation circuit of the processing plant have been sampled and selected physicochemical parameters measured. Plant survey has followed with selected samples being submitted to chemical and mineral liberation analysis (MLA). Quantitative analysis has enabled to link the flotation response of the main copper bearing minerals with the physiochemical parameters (Eh, DO, EDTA extractable Cu and Fe) being measured. Data from mineral liberation analysis have highlighted an important process-relevant information explaining the variable floatability in the various granulometric classes.The low chalcopyrite recovery could be linked to the non-sufficient amount of dissolved oxygen measured in the copper selective flotation circuit. Given the finely-disseminated copper bearing minerals, the observed copper losses inside the flotation could be also explained by the insufficient degree of liberation of the mineral grains. Suggestions for improvement in metallurgical recovery are therefore proposed.INTRODUCTIONThe Chelopech deposit is a high sulphidation (HS) volcanic-hosted epithermal gold-copper one (Bonev et al., 2002; Kouzmanov et al, 2009; Strashimirov, et al. 2002). It is considered as one of the largest and economically most important ore deposits in Bulgaria, and often is referred to Europe’s largest currently mined gold resource. It is located in the Srednogorie zone which falls within the central sector of the extensive Tethyan–Eurasian copper belt formed during Mesozoic–Cenozoic time along the southern margin of Eurasia. The Srednogorie zone has the features of a subduction-related, Late Cretaceous island arc, with development of a calk-alkaline to sub-alkaline volcano-plutonic magmatism and associated hydrothermal ore deposits (Bonev et al., 2002).Ore mineralisation occurs as veins, mineralised breccia, and massive sulphide replacements and disseminations. Three mineralisation stages have been identified at Chelopech deposit. An early pyrite–marcasite mineralisation (stage I), characterised by disseminated to massive pyrite followed and partly brecciated by enargite– luzonite and tennantite–chalcopyrite–bornite (stage II) mineralisation. The gold is met in ores belonging to the stage II. Sphalerite and galena veins characterise stage III, best developed below and adjacent to the pyrite–enargite– chalcopyrite ore bodies. Advanced argillic altered volcanoclastic rocks commonly host this mineralisation (Chambefort and R. Moritz, 2014)."

Jan 1, 2018

By T. Kalliomäki, B. P. Wilson, M. Lundström, J. Aromaa

"The diffusion rate of cupric ions and limiting current density are important factors that determine the rate of copper mass transport in electrolytes proximal to the electrodes used in electrorefining. As a result, diffusion can critically affect both the rate and morphology of cathodic copper deposition thus it is crucial that it be optimized to improve production. The diffusion coefficient of cupric ion (DCu(II)) in copper electrolytes has been previously researched but the values detailed in the literature show significant variations and moreover arsenic - a common impurity - is generally excluded. Furthermore, in previous studies the models have only been validated with model solutions rather than real industrial electrolytes. In this study, the limiting current density and diffusion coefficient of cupric ion were defined for synthetic electrolytes containing copper (40–60 g/dm3), nickel (0–20 g/dm3), sulfuric acid (130–220 g/dm3) and arsenic (0–45 g/dm3). The cathodic limiting current densities were measured at three temperatures (50, 60, 70 °C) with linear sweep voltammetry using the rotating disc electrode (RDE) as the cathode. By utilizing the limiting current density, models for the DCu(II) were constructed using Levich equation. The validity of the DCu(II) model for technical use was confirmed through their application to industrial electrorefining electrolytes with known element concentrations. Consequently, this study provides improved models for limiting current density and diffusion coefficient of cupric ion that can be utilized for industrial process optimization.INTRODUCTION For an efficient copper electrolysis process, the mass transport of copper ions has to be effective. The main mass transport types of copper ions in copper electrorefining are diffusion and convection while the effect of migration is negligible. In order to maintain good mass transport, it is beneficial to keep both diffusion and convection high to maximize the limiting current density (jlim) (Bard and Faulkner, 2001). jlim is the highest current density that can flow in the cell, thus by maximizing jlim the deposition time can be minimized and the good quality of cathodic copper can be ensured. Thus, it is important to research possibilities that increase and optimize diffusion since slow diffusion in the electrolyte limits the rate of the process (Levich, 1947; Moats et al., 2000; Subbaiah and Das, 1989)."

Jan 1, 2018

By Ghislain Bournival, Seher Ata, Yesenia Saavedra Moreno

"In froth flotation, frothers are commonly added to the pulp to modify the interfacial properties of bubbles. These properties may be modelled using commercial software packages. The volume of fluid (VOF) method was used to track the interface motion during the post-rupture oscillation of the newly formed bubble. The aim of this research was to provide an insight into the modelling of viscoelastic phenomena through the coalescence dynamics of two bubbles in solutions of various surface tensions using a three-dimensional computational study. Different surface tensions were used to simulate changes in methyl isobutyl carbinol (MIBC) concentration. The computational results were evaluated by comparing the oscillatory motion of the newly formed bubble for the different MIBC concentrations. For known MIBC concentrations, the deviations of the computational damping constants from the experiments were quantified, revealing possible deviations in the VOF surface tension model. In addition, the analysis provided insight into the flow patterns of the fluids and the effect of frother concentration on the velocity magnitude in the liquid and air phases. This analysis provided a further understanding of coalescence dynamic and post-rupture oscillation in froth flotation. Such a process is believed to affect the stability of particles at the interface due to the release of surface energy when the bubble surface is reduced.INTRODUCTIONIn froth flotation, bubble coalescence is an important phenomenon that is highly dependent on the interfacial properties of the liquid film between bubbles. It is known that the frother adsorption at the liquid-air interface changes the interfacial properties, resulting in an increase of surface elasticity and reduction of surface energy (Pugh, 2016). However, in a two-phase system, it is the frother type, concentration, frother molecule conformation at the interface and viscoelasticity properties that dictate film stability (Pugh, 2016). The viscoelasticity phenomena can be explained in terms of surface viscosity, which quantifies the energy dissipation at the interface and surface elasticity that represents the energy absorbed at the interface as result of bubble oscillation (Pugh, 2016). In the case of systems stabilized by nonionic frothers, it has been shown that surface forces alone do not account for film stability, but the viscoelasticity represents the ability to damp external disturbances in the thin film separating the bubbles and prevent coalescence (Santini et al., 2007)."

Jan 1, 2018

In this study, the effect of different regrinding conditions on the depression of fine pyrite at the cleaner stage was investigated. The regrinding condition was manipulated by different grinding media (stainless steel and mild steel). It was found that, after regrinding with mild steel grinding media, the flotation recovery of pyrite is much lower than with stainless steel grinding media, and it is more difficult to restore by adding more collector. This means, regrinding with mild steel grinding media may depress the pyrite flotation more strongly. The copper activation of pyrite during regrinding was also studied. It was found, that the mild steel grinding media is more beneficial to copper activation than stainless steel media. Additionally, the analysis of pyrite recovery based on the particle size indicated, that the fine particles are more susceptible to different grinding environments than intermediate and coarse particles and therefore are more easily depressed with a much lower flotation recovery.

Sep 1, 2012

By Shankar Narasimhan

Data Reconciliation is a technique that has been developed more than fifty years ago and has been refined over the past few decades and is now being deployed routinely in several process industries such as refining, petrochemicals, mineral processing and hydro-metallurgy. The main aim of data reconciliation is to obtain accurate estimates of all process variables and parameters using plant measurements that are invariably corrupted by random errors and may also contain biases and other gross errors. These reconciled estimates are subsequently used either in offline exercises such as simulation, retrofitting, and de-bottlenecking or in online optimization and control applications. In this presentation, a brief history of data reconciliation, its assumptions, and uses are reviewed. New developments in this area are highlighted such as robust data reconciliation to handle model uncertainties, and reconciliation using identified models derived purely from operating data. Several case studies in the mineral processing industries reported in the literature are also highlighted to indicate the potential benefits of this technique.

Sep 1, 2012

By Aysan Molaei, Kristian Waters

Solvent extraction is the one of the most common techniques for recovering metal ions from aqueous systems. However, solvent extraction is generally limited to relatively high concentrations and requires a high organic/aqueous ratio. Predispersed solvent extraction (PDSE) is a modified solvent extraction technique in which organic phase pre-comminute into colloidal liquid aphrons (CLAs). In this study, the extraction of copper by the PDSE technique was investigated. The main focus of this study was on the characterization and application of CLAs under different operational conditions, and to investigate the effects of different parameters on the extraction of copper from dilute streams. The effect of phase volume ratio (PVR), initial extractant concentration and initial copper concentration was studied and optimum conditions for copper removal using CLAs generated by LIX 984N as extractant, oil soluble and water soluble surfactants were proposed. Experimental observations indicated that PVR and initial LIX 984N concentrations had a significant effect on copper extraction.

Jan 1, 2014

By Z. K. Liu

Production of the amines collectors and cation surfactant gives out the liquid waste containing high COD (Chemical Oxygen Demand) and ammonium & nitrogen, which is difficult to be treated with conventional biological method to reach waste water discharge standard. In this paper an engineering practice of treatment on the liquid waste by the technology combined high effective microbe with A/O (Anaerobic& aerobic) process is presented. The quality of treated water comes up to class A national standard for discharging.

Jan 1, 2014

By Khoi Nguyen, Keeney. Luke

An outcome of geometallurgical mapping and modeling is the generation of continuous down hole profiles of quantitative geometallurgical attributes at assay scale. These attributes vary from traditional routine data collection methods such as assays and geotechnical logging (e.g. RQD, Fracture Frequency), to new measurements including petrophysical attributes, EQUOtip hardness data and modelled estimates of metallurgical performance indices (e.g. A*b, BMWi). Highly variable continuous geometallurgical data with high spatial resolution can make it difficult to identify spatially continuous domains. An automated method based on the well known time series analysis technique of cumulative summation (CuSum) has been developed which uses statistical analysis to identify domain boundaries in down hole profiles. A bootstrap analysis automatically identifies potential domain boundary locations with associated confidence levels and a series of t-tests determine the significance of defined populations. A hierarchical clustering algorithm is applied to results to enable data sensitivity to be taken into account. The method produces statistically significant populations and is suitable for application on either single or multiple drill hole datasets. Domains defined preserve the inherent variability of the continuous down hole data and provide a method for scaling up variability data (i.e. assay interval scale) to mine scale (i.e. bench scale, block model). A major advantage of the method is its ability to provide suitable input parameters for use in spatial modeling of non-additive and non-linear geometallurgical attributes through conditional simulation techniques.

Jan 1, 2014

By P. Somasundaran, Tarun Bhambhani, D. R. Nagaraj, P. Patra

"Certain non-sulfide gangue minerals are known to have a detrimental impact on flotation performance. When present in large amounts, they can negatively affect recoveries, grades and throughput at current operations, and dictate the viability and economics of a greenfield project. While the problem has been investigated for a long period of time, a scientific understanding of the mechanisms by which these minerals affect performance has not been obtained, and no acceptable solutions have been developed. This study focused on the role of physical aspects (e.g. shape, size and morphology) of these minerals as the underlying factors for poor performance. Flotation tests were conducted using model mineral systems to understand the impact of one phyllosilicate mineral (muscovite) on flotation kinetics, recoveries and grades. The study reveals that there is minimal impact on flotation kinetics, no impact on recoveries, but a very strong detrimental impact on concentrate grade. The traditional mechanisms of transport such as true flotation and entrainment are not sufficient in describing the extent of transport. The transport is attributed to the unique platy morphology of muscovite, which allows for action of buoyancy forces. Data also shows that the muscovite transport to the froth phase depends on particle size; coarser sized muscovite (100-150@) report to the froth phase in much larger amounts than medium sized muscovite (50-100@). At finer size fractions (<50@) the transport is also increased, and this is likely due to entrainment."

Jan 1, 2014

By Bruno Miotto, Rômulo Faria, Juliano Tessinari

Water has become a scarce resource in quality. Companies in the ornamental stones sector use water intensively in their processes. Given the current situation, it is imperative to rationalize this feature that at first is related to the reduction of consumption, through no waste. Considering the importance of environmental preservation of this resource and also to comply with the legislation that governs the procedure for processing of ornamental stones, new alternatives should be available for the maintenance of current activities and promoting the development / growth of the sector on the principles of environmental sustainability and economic in future scenarios. The ornamental stones processing industries, especially for micro and small enterprises, require the implementation of a management program that optimizes the use of this feature, since the techniques individually processed operate with three basic processes dependent on water use: split blocks, polishing and finishing slabs. This technology requires a large volume of water resources and control measures involving operating recirculation water can minimize power consumption, making the company more profitable and sustainable. This is a survey of literature, associated with technical reports provided by partner companies and government agencies. This study aimed to the organization and provision of theoretical knowledge about the subject matter coupled with management practices, which may be reproduced or re-purposed in ornamental stones processing companies. It should be noted that the industry in question is lacking in technical expertise for the implementation of sustainable attitudes relevant to water reuse.

Jan 1, 2014

By S. S. Indi, D. Santhiya, S. Subramanian, K. A. Natarajan

The interaction of sphalerite- and galena-adapted Bacillus polymyxa (B. polymyxa) cells with sphalerite and galena has been investigated through adsorption, electrokinetic and flotation studies. The amount of sphalerite-adapted cells adsorbed onto sphalerite and galena is found to decrease with increase in pH. On the other hand, the adsorption density of galena-adapted cells onto galena shows a feeble maximum around pH 8, and the amount adsorbed is marginally higher in the pH range of 8 . 11.5. In the case of sphalerite the amount adsorbed is higher in the pH range of 5 to 7. The adsorption density of both the adapted cells is higher for galena vis-à-vis sphalerite. Electrokinetic studies indicate significant surface chemical changes on the chosen sulphide minerals and the bacterial cells, consequent to microbe-mineral interactions. Selective flotation tests on a synthetic mixture of galena and sphalerite confirm that sphalerite could be preferentially floated from galena, using galena-adapted cells, but not with the sphalerite-adapted cells. The cells adapted to galena assay a higher polysaccharide content and are more hydrophilic, while the sphalerite-adapted cells have a higher protein content and possess greater hydrophobicity. Transmission electron microscopic studies have provided evidence in support of the differences between the cell surface characteristics of sphaleriteand galena-adapted cells. Possible mechanisms of interaction between the minerals and the bacterial cells are discussed.

Jan 1, 2003

The multiphase reactive flow in process equipments, e.g., fluidized beds, involves complex coupling among flow, heat/mass transfer and reaction kinetics over scales ranging from molecules, sub-particle eddies up to equipments. Conventional two-fluid model (TFM) ignores the effects of meso-scale structures, and hence is not suitable for describing the intrinsically multi-scale heterogeneity within multiphase systems. To be consistent with the meso-scale structures, we proposed a structure-dependent multi-fluid model (SFM), and closed its drag force by invoking the stability condition defined in the energy-minimization multi-scale (EMMS) model. The SFM reduces to the conventional TFM if homogeneous distribution is assumed within each grid, and reverts to the original EMMS model if it is used to describe steady state, global behavior of a fluidized bed. The multi-scale computational fluid dynamics (CFD), which is characterized by coupling of TFM and EMMS drag coefficient, can then be viewed as a simplified solution of the SFM. We demonstrate that the SFM significantly improves the accuracy and at the same time enables higher efficiency by allowing for much coarser grid used than in TFM. Recent application of it in 3D, full-loop, transient simulation of industrial processes, e.g., fluid catalytic cracking, coal combustion, gasification, kaolinite clay calcinations, shows its advantages in helping design and troubleshoot multiphase process equipments. Keywords: meso-scale, SFM, multiscale, CFD, simulation, structure, EMMS, multiphase

Sep 1, 2012

By Rod Hocking, Nenad Djordjevic, Randolph A. Pax

"Semi Autogenous (SAG) mills exhibit complex interactions between the competing processes occurring inside them. Traditional practice is to measure the operation of the SAG mill at an operating point in terms of its load, power draw, steel load, ore throughput rate, water addition rate and more recently the feed size distribution. However significant departures from standard operating conditions can occur due to changes, for example, in ore hardness, the quest for higher throughputs and mill internal geometry.Fundamental to SAG mill operation are the internal processes that ensure that a given ore throughput and discharge grind size are obtained and maintained. These processes can be measured using non-contact acoustic emissions from the SAG mills which utilise the information rich noise that is generated by the grinding, flow and charge movement inside a mill. Non-contact acoustic in-mill measurements can be used for the validation of discrete element modelling (DEM). The approach being reported in this paper is to model the individual physical processes and then to simulate the assemblage of a physical mill.By modelling the assemblage of steel balls in a 1 m pilot mill and duplicating the experimental arrangement exactly with the PFC3D model, it was possible to replicate the pattern of falling balls in the experimental mill. Furthermore, the key features of the acoustic emission spectrum that was calculated compared favourably to the measured spectrum at low frequencies. The experimental spectrum is shown in Figure 1."

Jan 1, 2003

By Hermann Wotruba, Günther Petzold, Fabian Riedel

"Increased mechanisation in underground stoping operations results in higher waste percentages in the run of mine. A new method of reducing waste is optical sorting. Optical sorters have low maintenance costs and can be installed underground. Rejected waste can be backfilled directly without having to be hoisted, therefore increasing the shaft capacity for a more valuable product.Pre-Concentration by Optical SortingA series of tests with an optical sorter at the laboratories of the Mineral Processing Department of the RWTH Aachen have shown that it is technically possible to sort almost any kind of ROM which has even slight differences in colour (Platinum Reef, Gold Reef, Kimberlites, Industrial Minerals, Aggregates, Coal, etc.).By employing optical sorters prior to conventional sorting plants, it is possible to reduce the amount of waste to be processed. Figure one shows test results of ore which is quite similar in colour to the waste (grey). With increasing waste in the feed, the recovery decreases according to the graph. In figure two, the colour ratio parameter of the sorter was changed with a similar waste percentage in the feed. With increasing loss of ore, a cleaner product can be achieved. Before designing a sorting plant, throughput tests are carried out. Throughputs vary depending on the size range. Sorters have a limit of about 250mm top size. In this upper size range, throughputs of 150t/h can be achieved."

Jan 1, 2003

By Braulio Fernández, Gonzalo Montes-Atenas

The mechanisms taking place during the leaching of chalcopyrite when it is in contact with other sulphide minerals, particularly pyrite, is still a matter not fully solved. The GalvanoxTM process has been suggested as a technical solution to extract copper from chalcopyrite at commercial level inhibiting passivation phenomena widely reported in the field of hydrometallurgy. Although there are some rules of thumb to implement the process in fully flooded stirred conditions such as the optimum pyrite/chalcopyrite mass ratio should range between 2 and 4, among others; not many research studies have been focused on studying the effect of the surface contact area between chalcopyrite and pyrite to reach best copper extraction conditions. The aim of this research work is to determine the effect of the mass and surface area ratio between the minerals and the surface area of each of the mineral species effectively exposed to the leaching solution. Different particle size distributions of pure pyrite and chalcopyrite were contacted in a small column and subjected to leach tests using sulphuric acid in concentrations close those used in industrial leaching processes. The mass and surface area ratios which optimize the extraction of copper were obtained. It was found the best leaching conditions were attained when reducing the particle size distribution of pyrite regardless the particle size distribution of chalcopyrite. Further analysis concerning the impact of this conclusions on large scale operations is presented and discussed.

Jan 1, 2014

By Arkady Senchenko, Sergey Rybkin, Svetlana Lagutina, Alexander Aksenov

"The application of autoclave leaching technology for sulfide concentrates containing base and precious metals is spreading nowadays. In the process of concentrate leaching silver accumulates in the forms of jarosite compounds, transforms into passive chemical state and is mostly found in solid phase. Effective technology development for silver recovery from solid phase of concentrate autoclave leaching raises the necessity for comprehensive research of silver forms of occurrence and their properties in this product. The possible form of silver presence in solid phase of concentrate from autoclave leaching is argentojarosite. This article is about the study of physical and chemical properties and structure of synthetic argentojarosite - AgFe3(SO4)2(OH)6. A method of argentojarosite synthesis has been developed under conditions simulating the process of autoclave leaching of sulfide concentrate, and includes silver dissolution in the form of Ag2SO4, Fe3+ sulfate hydrolysis, AgFe3(SO4)2(OH)6 build-up. X-ray data of synthesized argentojarosite completely correspond to those of natural mineral. The structural formula of this complex compound has been studied and confirmed. The chemical properties of synthetic argentojarosite, its solubility in acids and alkalis have been investigated. Phase changes when heating argentojarosite with the build-up of iron oxide and silver sulfate stop at the temperature equal to 650 ºC. Based on the results received the technological solutions for solid phase of autoclave leaching of concentrate have been developed with silver recovery of 95-98%."

Jan 1, 2014

By l. Bokányi

The Tóth adsorption isotherm-model was derived and has been widely used for the inhomogeneous solid surfaces. The three-parameter isotherm incorporates the effect of the interaction between the adsorbed substances as well. All classic isotherms can be clearly derived from this one. Bio sorption is the sequestering of dissolved components, i.e. metals from aqueous environment by living or dead biomass. The mechanism of the bio sorption is rather complex and predominantly refers to the affinity of the chemical functional groups: acetoamido groups of chitin, structural polysaccharides, amino and phosphate groups in nucleic acids, amino, amido sulfhydryl and carboxyl groups in proteins, hydroxyls in polysaccharides, sulphates, etc. towards the soluted substances, i.e. heavy metal ions. The physisorption, as well as the micro precipitation also play a role in the bio sorption mechanism. The description of such a complex process can be done by a model reflecting this complexity. According to our experimental results obtained using algae and immobilised algae for the sorption of heavy metal cations, the most appropriate model for such a description is the Tóth-isotherm, the parameters of which can be interpreted. Flotation process is also a very complex one, the overall result of which depends on the partial results of many micro processes and parameters. The particle interaction which very often takes place in the flotation systems makes the phenomenon even more complex. The flotation kinetic model with particle interaction was already introduced by author in 1994 and further developed in 2000. In this paper based on the experimental results obtained with metallic particles in case of recycling flotation the applicability of the Tóth-equation as an alternative flotation kinetic model for the systems with particle interaction is proved. The applicability of Tóth-model in mineral processing demonstrates again its universal character and the great talent of the late József Tóth, great Hungarian chemist. Keywords : Toth isotherm, biosorption, Langmuir, Freundlich

Sep 1, 2012

By Prasad Kopparthi, Narasimha Mangadoddy, Balraju Vadlakonda, A. K. Mukhurjee

"Column flotation (CF) is widely using equipment in mineral processing since many years. ERT technique facilitates non-invasive and nonintrusive visualization of different parameters in a column flotation. The hydrodynamic investigation of two-phase and three-phase are studied using electrical resistance tomography (ERT) coupled with pressure transducers in laboratory column flotation. Bubble size distribution, bubble rise velocity and gas and solids hold-up hydrodynamic parameters are extensively used for column flotation design and performance evaluation. The effect of superficial gas velocity, feed flow rate and slurry height in the column on mean gas holdup and its radial distribution has been analyzed for both two-phase and three phase flows. Influence of solid particles on gas hold-up and its complete elucidation of dispersion in the column are assessed. Experimental results show ERT coupled with PT is suitable as an online monitoring tool to provide useful information on the hydrodynamic parameters of the column flotation. The tomography images, which were generated using a modified sensitivity back projection algorithm, were employed to explore the parameter study in two phase flow in column flotation.Furthermore, the results also indicate that ERT is a very powerful tool for diagnosing the `inside' flow behaviour of gas-liquid two-phase and gas–liquid–solid three-phase flow in columns operating for a wide range of flow regimes.INTRODUCTIONColumn flotation is being widely used in mineral processing industries in last few decades. There has been a lot of attention in understanding the various mechanisms that take place during operation. Due to many reasons column flotation has been found to yield better performance than conventional flotation cells, particularly with fine particles (Honaker and Mohanty, 1996). Due to complex flow behavior in nature, the performance of column flotation influenced by on flow characteristics of the process, (Jena et al., 2008). Previously basic research is carried out to evaluate the effect of the hydrodynamic parameters on the performance of column flotation (Finch and Dobby, 1991)."

Jan 1, 2018

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