Search Documents

In recent years there have been significant advances in computing power and programming technology, enabling more complicated processes to be modelled in greater detail. This has led to increased understanding of these complex processes, including mineral processing circuits. JKSimMet was developed over the past 35 years, incorporating outcomes from various research projects conducted at the Julius Kruttschnitt Mineral Research Centre (JKMRC), part of the Sustainable Minerals Institute at The University of Queensland. It is a steady-state software program that allows the user to mass balance, model fit and simulate crushing, grinding and classification circuits. JKSimFloat has been developed over the past 20 years, incorporating outcomes from the JKMRC and other institutions around the world through the AMIRA International P9 project ? The Optimisation and Simulation of Mineral Processing Circuits. It is also a steady-state software program, similar to JKSimMet, but is focused on analysing flotation circuits. Both programs have been used in a wide variety of optimisation and design projects over many years. This paper provides an overview of both JKSimMet and JKSimFloat and examples of how both programs have been used to better understand, develop and optimise mineral processing circuits around the world.

Feb 27, 2013

By Luis Bernal, Pertti Lamberg

"SkimAir®, Outotec’s Flash Flotation unit, is placed in the grinding circuit to float the valuable liberated mineral particles from the circulating load of the grinding process while still relatively coarse. Esperanza copper-gold deposit located near Sierra Gorda, Chile, hosts 535 Mtn ore grading 0.57 % copper, 0.27 g/t gold. A series of laboratory and pilot flotation tests were conducted to find out the flotation performance of different ore types both with and without the Flash Flotation. This data was used in modeling and simulating to find the optimal flow sheet for the highest Au and Cu recoveries with required Cu concentrate grade. Simulations with the HSC Chemistry 6.1 software showed that with the SkimAir® and Flash Cleaning it is possible to achieve the required copper grade of 27% Cu with 5% higher recovery for gold and 2% higher recovery for copper compared to the circuit without the Flash Flotation. The pilot tests with simulations were the back up for the installation of four Outotec SkimAir® SK 2400 Flash Cells in the Esperanza circuit for the processing 95.000 tpd of copper-gold ore.INTRODUCTIONEsperanza copper-gold deposit is located near the town of Sierra Gorda, 2nd Administrative Region in Chile. Deposit hosts 535 Mtn ore grading 0.57 % copper, 0.27 g/t gold. The main copper minerals are chalcopyrite and bornite. Gold is carried by native gold and solid solution of copper.Antofagasta Minerals S.A. conducted a series of flotation studies at SGS Minerals Services. The studies included both laboratory and pilot scale work. In the pilot scale study the applicability of the Flash Flotation technology was studied with High-grade Andesite and High-grade Porphyry samples (SGS Mineral Services). The Flash Flotation technology gave significantly higher recovery for gold than the conventional flotation circuit. The recoveries for copper were slightly higher using the conventional flotation."

Jan 1, 2009

By V. Buddhiraju

Liberation of valuable minerals from low-grade ores is an important challenge faced by the minerals industry. It is necessary to grind the ore to extremely fine size in order to liberate the mineral from such ores. As it is generally difficult to produce ultra-fine particles of diameter less than 1-10 mm using traditional grinding equipment such as the ball mill and the rod mill, several new types of grinding equipment are being explored to produce ultra-fine particles. The planetary ball mill is promising in that it makes grinding to submicron sizes possible by imparting high energy to the ground powder. It consists of four pots mounted on a disc. The pots and the disc rotate in opposite directions. The high-speed rotation of the pots and the disc makes the grinding media move at high centrifugal velocities leading to high impact energy and thus improving the grinding performance. The planetary mill is currently being explored for its capability to produce ultra-fine particles on a large scale. In this context, there is a need to understand the dynamics of the ultra-fine grinding process within the planetary ball mill. Here, we address this need by studying the movement of media inside the pots and the grinding kinetics, employing both experimental and modeling techniques. The experiments were conducted on a laboratory scale planetary ball mill. The particle size distribution of the powder was monitored to allow an estimate of the grinding kinetics to be made while video images of the media profiles inside the pots were taken to allow quantification of the media movement as a function of process variables such as the media size, percent loading, etc. The population balance modeling (PBM) framework was used to simulate the kinetics of grinding while the discrete element modeling (DEM) approach was utilized to simulate the media motion. The parameters in the population balance model, namely, the specific breakage rate and the breakage distribution function were determined by fitting the batch grinding equation to the experimental data. The predicted powder size distributions were found to be in good agreement with the experimental data. The experimental results have shown that the mean size of product powder decreases as the size of the grinding media decreases. Keywords: population balance modeling, grinding, planetary mill, grinding media, grinding equation

Sep 1, 2012

By A. Oleynik

The paper presents the modeling and simulation experience in field of research into the mineral row material separation. Features of the subject domain generate a need of different approaches using for modeling and simulation. The general scheme of separation modeling process and some examples of developed simulators are presented. Due to different features and standards of knowledge, different formal models should be used to describe and analyze separation processes in the frame of general scheme. The main types of simulators are modeling tools based on statistical, analytical or expert (logical) methods. Advantages and disadvantages of each method above are mentioned in the paper. Based on the studies in the field considered, an aggregated mathematical model of a separation process has been generated. The aggregated model is presented by two combined equations. The actual apparatus model is formed by substitution of the particular functional expressions of the forces acting in the apparatus into the aggregated model. The developed simulation tools are based on integration of limited experimental data and on ex-pert knowledge available on separation processes. Simulation does not allow producing an optimal, in terms of mathematics, solution. But it gives an opportunity to produce and compare quite adequate solution variants in real time. The results obtained are the theoretical and methodological base for integrated information technologies to support studies in the field of mineral processing.

Jan 1, 2014

By D. Minson, G. Lindner, S. Miskovic

Thickening is an essential process in the mining industry and currently one of the most costeffective technologies for recovering water from tailings streams. Improved performance of thickeners can lead to significant energy savings, and increased recovery of water close to the mill. Its performance is directly affected by the rheological properties of slurry at the raking zone. To further understand the dynamics of slurry inside thickeners, a multiphase particle-in-cell method (MP-PIC), is used to investigate its behaviour on a sliding channel. Parameters like particle-size distribution (PSD), rheological properties, and domain discretisation were benchmarked to study its influence on highly dense flows. The slurry concentration by weight of the samples ranged between 60 to 70% of solids by mass, under which conditions it is known to behave like a non-newtonian fluid. An experimental flume, 0.2 m wide and 2 m long, was used for validation under different angles and slurry composition. Visual analysis of the slurry height and shape, as well as velocity measurements, were compared with numerical data. The preliminary results show the ability of MP-PIC method to represent different types of highly dense flows, while still being computationally affordable and practical for industrial use. Keywords: High density thickener, CFD, MP-PIC, multiphase-flow

Jan 1, 2020

By Neale R. Neelameggham, Ganesan Subramanian

Anthropogenic Heat Flux (AHF) is defined as heat generated from use of non-renewable energy by humans in residential, commercial and industrial activity. State of the art Global Circulation Models (GCM) used to predict global impacts of climate change to help set related policy do not include this heat flux. AHF is considered to be a much smaller an insignificant forcing, solely based on the global mean values. Various Climate models have discussed the “Heat Island Effect ” in high AHF areas. Through the use of simple climate models, it is shown that significant pockets of AHF can impact climate locally in ways significant enough to impact global climate. Inclusion of spatial AHF distribution, rather than a global mean value is important to accurate modeling either on a regional or global scale.

Mar 1, 2017

By T. Doshi, S. Ennis

The interactions of mineralogy, chemistry and mechanics in an ore concentrator are highly complex, so it can be particularly hard to measure the impact of a specific plant change. This paper describes a method for modeling a baseline metal recovery so that the value of a major process improvement can be measured. This method only considers variables that cannot be optimized, such as mineralogy or equipment availability. Linear regression, random forests, extreme gradient boosting (XGB), and simple two-layer neural network models of process value were tested as candidate baseline model types to assess the impact on throughput and recovery in a copper concentrator. XGB performed the best in terms of accuracy for both throughput and recovery models. Using an XGB model of nonoptimizable process variables provides an accurate, dynamic baseline that can be compared with continuous production data to determine the value added by major concentrator improvements in real time.

Mar 2, 2022

By Cristian Castro, Patricio Gomez, Jason Furtney, Augusto Lucarelli

Understanding the effect of blast-induced damage on wall stability is critical to safe and economical operation of a large open pit mine. In this paper, a numerical modeling methodology is used to investigate damage to rock at distances of greater than 10 m (30 ft) from a blast. Two applications of the model are demonstrated. First, the method is applied at a zinc mine to investigate the effects of pit blasting on underground stopes. The model predictions are compared with the scaled distance formula. Next, the model is used to investigate a common problem in which unfavorable jointing creates a rock wedge on the bench above a blast. The numerical predictions of wedge stability are compared to the wedge energy model described by Wong & Pang (1992).

By Cristian Castro, Patricio Gomez, Jason Furtney, Augusto Lucarelli

Understanding the effect of blast-induced damage on wall stability is critical to safe and economical operation of a large open pit mine. In this paper, a numerical modeling methodology is used to investigate damage to rock at distances of greater than 10 m (30 ft) from a blast. Two applications of the model are demonstrated. First, the method is applied at a zinc mine to investigate the effects of pit blasting on underground stopes. The model predictions are compared with the scaled distance formula. Next, the model is used to investigate a common problem in which unfavorable jointing creates a rock wedge on the bench above a blast. The numerical predictions of wedge stability are compared to the wedge energy model described by Wong & Pang (1992).

Feb 1, 2020

By S. Dixon

The Phelps Dodge - PD/PERI Heap2D Model is a phenomenological stockpile leaching model that has been used to investigate the best operating practices for leaching chalcocite. The model was used to investigate the interrelationships of raffinate application rate, aeration rate, raffinate iron grade, acid cures, lift height, crush size, bulk density, air injection spacing and raffinate drip emitter spacing.

Jan 1, 2006

By Nader Golchinfar

The V-shaped breakout zones around excavation boundary have been captured by a variety of numerical methods such as FEM and DEM methods and failure criteria such as linear Mohr-Coulomb and non-linear Hoek-Brown failure criteria. Diederichs (2007) used Hoek-Brown parameters in his spalling parameter approach; however, it is convenient to use the Mohr-Coulomb failure criterion in numerical modeling because the material parameters can easily be determined from triaxial tests. It is shown in this study that when peak and residual strength parameters as well as dilation angle are properly defined, the brittle Mohr-Coulomb failure criterion can be used to model depth of failure near underground excavations satisfactorily.

Jun 1, 2012

By B. R. Rakishev

When transferring from the railway transport to the combined automobile-conveyor transport by stripping works at coal mines, it was proposed to form at the working zone a concentration horizon for reloading part of the stripping volumes at the railway transport. The procedure is adduced for substantiating parameters of bloc-panels, depending on the step of transferring of reloading points in interrelation with the intensity of the working zone?s development with observing the balanced movement of parts of the stripping zone relative to the temporary concentration level.

Aug 1, 2013

By Xiaozhen Song, Ashish Ranjan Kumar, Shimin Liu

Diesel particulate matter (DPM) presents considerable occupational health hazards in large-opening underground mines (LOUM). Previous studies have primarily focused on single-factor ventilation optimization, while this study innovatively integrates booster fan configurations and loading truck placements utilizing a coupled CFD turbulence– particle tracing model. Three distinct ventilation scenarios were simulated to evaluate their effects on DPM dispersion and dilution efficiency. Results indicate that dual booster fans substantially enhance airflow distribution, increasing air exchange efficiency; The migration pathways of DPM closely align with airflow trajectories, and air turbulence is the predominant factor influencing particle behavior within LOUM; Loading truck placement critically influences local DPM accumulation; Airflow effectiveness diminishes markedly with distance from booster fans, causing stagnation at distant working faces. The coupled CFD approach accurately captured airflow–particle interactions, offering valuable insights for optimizing mine ventilation strategies and effectively mitigating miners’ exposure to DPM.

Feb 22, 2026

By R. Panchadhara, P. A. Gordon

"A non-local formulation of the classical continuum mechanics theory called peridynamics is used to study driven dynamic fractures in the near wellbore region. The controlled application of pressure on the inner walls of the borehole can create multiple fractures that initiate from the borehole and propagate radially outwards. The number of fractures and their final length is primarily controlled by the loading conditions and the in-situ stress state. It is the objective of this study to use a novel numerical technique called peridynamics to investigate: 1) the influence of loading conditions on the fracture pattern, and 2) the effect of in-situ stress anisotropy on fracture propagation. In peridynamics, the momentum equation from the classical theory of solid mechanics is replaced with the integro-differential momentum equation of peridynamics. As a result, discontinuities in the body, such as fractures, occur spontaneously. The material is represented by a continuum body with material points that interact non-locally with all other points within a specified distance. Fractures evolve using a damage evolution law, wherein, damage at any point is represented as the ratio of the number of broken bonds in the current configuration to the total number of bonds in the undeformed configuration.Simulations are performed on simplified test cases and the results are compared to relevant experimental and numerical studies found in the literature. It is shown that the loading rate significantly influences the number of fractures initiated from the borehole as well as their extent. Low loading rate produces fewer fractures, whereas, high loading rates produce longer and a greater number of fractures from the borehole. The numerical method successfully predicts the fracture pattern obtained from a wide range of loading conditions. Predictions also indicate that the resultant fractures are shorter when in-situ stress is applied. Numerical results show that the fracture extension is largest in the direction perpendicular to the minimum compressive stress for the range of loading rates used in this study."

Jan 1, 2015

By R Yang, M Lownds

This paper reports studies showing the effects of delay scatter on the peak particle velocity (PPV) of blast vibration using a multiple seed wave (MSW) vibration model developed in recent years. The case studies also show the advantage of accurate delay timing on frequency control of blast vibration. The site parameters used in one of the studies are from a case in which the MSW model has been tested and confirmed to yield accurate predictions. This study uses an open pit blast design and a tunnel design in which blast vibrations using electronic detonators are compared with those from pyrotechnic detonators. Mean peak particle velocity (PPV) and the 97.5% upper bound PPV are compared. It is found that both these values increase significantly with the increase of the coefficient of variation of the site parameter. Using an open pit blast design, the modeling shows that electronic initiation shifts vibration energy to the expected frequency for a selected point of interest. However, for the same design of nominal initiation timing pyrotechnic detonators fail to shift vibration to the expected frequency due to its large scatter of timing.

Jan 1, 2011

By A. Saigal

A finite element analysis based investigation using DEFORM software was undertaken to model the extrusion of 15% alumina and 15% silicon carbide reinforced aluminum metal matrix composites. The study dealt with predicting the extrusion load as a function of ram displacement and initial billet temperature. The finite element simulation was also used to determine the maximum extrusion load, temperature distribution, maximum temperature, strain rate, velocity and stress distributions in the billet during the hot extrusion process. The constitutive equation used was of the type: Z = E exp(Q RT) = A{sinh(ao-)}" where Z is the Zener-Hollman parameter (temperature compensated strain rate), s is the strain rate, Q is the activation energy, R is the gas constant, T is the absolute temperature, a is the flow stress and A, a and n are constants. Based on the work of Davies et al and McQueen and Konopleva, the values of Q (KJ/mol), a (MPa-1), n and A for aluminum 6061 composites reinforced with 15% alumina and 15% silicon carbide used are 181, 0.024, 4.15 and 1.33x10'2 and 233, 0.052, 2.60 and 1.26xl0'4, respectively. Even though the individual parameters are quite different, it was found that the predicted maximum extrusion load, temperature distribution, maximum temperature, strain rate, velocity and stress distributions are quite similar and differ by no more than 1-2 percent between the two composites indicating that the deformation behavior is primarily controlled by the matrix material and is not a function of the type of reinforcing ceramic particulates.

Jan 1, 2004

MODELLING FOR DESIGN AND CONTROL AT KALGOORLIE NICKEL SMELTER Modelling at KNS has been driven by a need to have better tools to answer questions about control, design and optimisation of the plant. This paper describes a number of modelling tools developed for use at KNS and how the development of modelling tools fits with a design and control strategy. In the area of control, a thermodynamic model has been developed to allow process parameters to be calculated. The model allows for the automatic selection of data from a laboratory data base, desired operating set points are input and the operating parameters are calculated. The model has the ability to adjust the burner efficiency and the furnace heat loss depending upon conditions. To assist with the testing of proposed furnace design changes, various mathematical models of the burner and shaft have been created using PHOENICS. These models are particularly well suited to combustion heat transfer and flow aspects of both the burners and shaft. Validation of these models have been based on novel sampling and measuring techniques in the shaft and also by using scale and cold models. Overall plant energy and cost structures can be analysed using a flowsheet simulation, developed to evaluate potential alterations to the plant flowsheet. This model system allows the global evaluation of process changes.

Jan 1, 1993

By Russell Barr, Robert E. Cook

As most mining companies are taking advantage of the high market prices for base metals such as steel, copper and zinc, numerous challenges still face them. While demand is strong from countries like China and the United States, demand is still mixed in Europe. This leads to the question of how long these high prices and overall high global demand will last. The tremendous market prices of most precious and base metals have enabled many companies to pay down a significant amount of debt accrued during acquisition and capital expenditure phases. Additionally, mining companies have been able to internally fund and engage in large projects to address future growth, (Less, 2006), thus maximizing the performance of the existing asset base, driving maintenance costs down, and improving and empowering the work-force (Mining Journal, 2006). Market prices and general efficiency gains have addressed these issues to some extent. However, critical issues still face operating managers in the mineral pro-cessing side of the industry. Energy prices, for example, continue to have a negative impact on cash costs for producing an ounce of gold or ton of a given base metal (Cummins, 2003).

Jan 1, 2007

By Raghavan Srinivasan, Y. V. R. K. Prasad, William G. Frazier, Steve C. Medeiros

"Aerospace gas turbine disks operate in an environment of relatively high stresses caused by centrifugal forces and elevated temperatures. These severe conditions necessitate the need for materials with good high temperature strength and low cycle fatigue resistance. One class of alloys used for this task is the nickel base superalloys, out of which, lN 718 is the most widely used in the aerospace industry [I]. The properties of IN 718 are attributed to the combined effects of the chemistry, heat treatment, and microstructure. The chemistry is tailored not only for solid solution strengthening but also for precipitation hardening developed during heat treatment, which combined with a fine grained microstructure lead to excellent mechanical properties such as low cycle fatigue resistance and elevated temperature strength. The properties of a gas turbine disk are sensitive to the microstructure, in particular the grain size, which is dependent on the processing history. The ability to precisely control the microstructural development during forging is dependent on controlling the process so that the workpiece is deformed within a ""safe"" region where no microstructural damage or flow instabilities occur. The microstructural mechanisms during deformation may themselves vary within the ""safe"" region and it is desirable to determine them within the range of parameters that are commonly used in industrial processing. The objective of this work is to establish a relationship between the grain size and the process control parameters i.e., temperature and strain rate, in the hot working of IN 178.BackgroundIN 718 belongs to a class of nickel-iron base superalloys with a face-centered cubic (FCC) austenitic matrix strengthened by precipitation of ordered intermetallic or carbide precipitates and is a nickel-rich alloy strengthened by ordered body centered tetragonal (BCT) precipitates, Y"". The alloying of nickel and iron results in the formation of the austenitic matrix, y, while the addition of chromium and molybdenum cause solid solution strengthening. Also, alloying with niobium causes precipitation of the metastable hardening constituent Y"" (Ni3(Nb,Ti)). Extended aging treatments in the temperature range 650 °C - 980 °C allow the more stable orthorhombic 8 phase (Ni3Nb) to form. Trace amounts of boron are also added to serve as grain boundary strengtheners for high temperature creep resistance. Hot deformation of IN 718 is performed between the Y"" and 8 solvus temperatures to facilitate recrystallization and to control grain growth. While deformation below the Y"" solvus (""""962°C) produces no change in grain size, deformation above the Y"" solvus but below the 8 solvus will lead to a dynamically recrystallized microstructure. The presence of the 8 precipitates will pin the grain boundaries after recrystallization and restrict grain growth so that a fine grain microstructure is obtained. Deformation above the 8 solvus ("""" 1038°C) will lead to a microstructure which has undergone grain growth during deformation."

Jan 1, 2000

By Donald Thompson, Jessie Bonner, Ileana Tibuleac

Predicting ground motion from complicated mining explosions is important for mines developing blasting programs in regions where vibrations must be kept below certain levels. Additionally, predicting ground motion from explosions in complex 3D media is important for source size estimation for nuclear treaty monitoring. Both problems have been addressed under the development of a new series of programs including: 1) Generalized Fourier Methods to generate synthetic seismograms (Green’s functions) in 3D media for a moment tensor source implementation and 2) MineSeis3D, a program that simulates seismograms for delay-fired mining explosions with a linear relationship between signals from small size individual shots. To test the programs, three-component seismic recordings of three production shots at an iron mine in northern Minnesota were compared to synthetic waveforms in 3D media. Comparison of the observed and synthetic waveforms shows promising results. The shape and arrival times of the normalized synthetic and observed waveforms are similar for most of the stations. The synthetic and observed waveform amplitude fit is best for the vertical components in the mean 3D model but not as promising for the transversal components. The observed effect of spall on the waveform spectra was weak in the case of fragmentation delay-fired mining explosions. Commercial applications of the code could provide data needed for designing explosions that do not exceed ground vibration requirements posed by the United States Department of Interior, Office of Surface Mining.

Jan 1, 2005