Search Documents

By Evgueni Jak

Recent advances in analytical, experimental techniques, and computer-based theoretical modelling of fundamental properties and elemental processes, provide new opportunities to develop the next level of whole-of-reactor pyrometallurgical furnace models. These models have the potential to significantly improve the prediction of, and adding value to, industrial operations. In non-ferrous smelting, the starting point of these models is the development of multicomponent thermodynamic databases for gas-slag-matte-speiss-metal-solids phases supported by systematic experimental research. The whole-of-reactor-models additionally should take into account kinetic processes taking place at micro- and macro- scales, and other key factors. Examples of applications of the latest research tools and modelling approaches to analysis of industrial flash and top submerged lance (TSL) sulphide smelting processes are presented. Different levels of industrial modelling are discussed from elemental local reactions, through general and more detailed whole-of-reactor-models, to plant sections and further to whole plant operation models. Some principles for development of pyrometallurgical reactor models are discussed.

By R Steinhaus, D M. Frantois-Bongarton, P W. Cleary, G K. Robinson, M D. Sinnott

Multotec Process Equipment and Dominique Frantois-Bongarton have a patent for a mineral sampling device, which they call a æTrue-Belt SamplerÆ. This uses a plough, which is rotated about an axis parallel to a conveyor belt to throw material off the belt. A fixed sample cutter takes a sample from the material thrown off the belt.Discrete element modelling has been used to consider a range of plough designs, plough speeds, positions for the cutter relative to the plough, cutter apertures, cutter designs and frictional properties, prior to the building of a physical prototype. This modelling work will be discussed, both as an illustration of the role of modelling in the development of new materials-handling devices and in order to explain the likely strengths and weaknesses of this new type of sampling device.CITATION:Robinson, G K, Sinnott, M D, Cleary, P W, Steinhaus, R and Frantois-Bongarton, D M, 2012. Modelling of a ætrue-belt samplerÆ - A new device for sampling from a conveyor belt, in Proceedings Sampling 2012 , pp 33-40 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 21, 2012

By P Edmondson, S Narendranathan, S Setchell

The complexities of mining in current economic conditions, where ore extraction schedules are varied on a weekly basis so as to ensure constant supply to the process plant, can result in a number of scenarios that pose geotechnical challenges. This paper presents a case study where alterations in the mining schedule of a particular mine resulted in the necessity to mine under a large uncompacted embankment, ie integrated tailings and mine waste landform. To effectively ensure that ore recovery was maximised, a well-designed slope monitoring system to evaluate design performance and to provide timely warning of impending failures was imperative. In order for an appropriate slope monitoring strategy to be developed, a thorough understanding of the potential slope failure mechanisms, capabilities and limitations of available monitoring techniques is required. Communication of the monitoring results and response strategy to alarms generated by monitoring systems must also be carefully considered to maximise effectiveness. The results of preliminary stability modelling can be used to guide the design of a monitoring system, which in turn can then generate data which can be used to iteratively refine the geotechnical model. This paper illustrates the development and implementation of a slope monitoring system and associated communication system, using a case study developed at a large iron ore open pit mine in the Pilbara region of Western Australia.CITATION:Narendranathan, S, Lee, E C, Setchell, S and Edmondson, P, 2015. Modelling of an integrated waste storage facility, in Proceedings Tailings and Mine Waste Management for the 21st Century , pp 183–190 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 27, 2015

By J. Sjöberg, T. Savilahti, S. B. Umar

"Sublevel caving in Malmberget mine results in mining-induced surface deformation. One of the currently mined orebodies is the Printzsköld orebody. As mining deepens there is need to assess the behaviour of the cave formed in the subsurface above this orebody. Numerical analysis was used to assess the effects of extraction to deeper levels and perform strength parametric studies. Stress redistribution was studied, and the results showed high stress buildups in the hangingwall and the crown pillar. Two failure mechanisms have been identified – shear and tensile. Reducing cohesion by 50% increased the area of the yielded zone by more than 100% in the hangingwall.IntroductionThe Malmberget mine is owned and operated by the Luossavaara-Kiirunavaara AB (LKAB) mining company. Iron ore is extracted by large-scale sublevel caving. The mine comprises a total of 20 orebodies, about 10 of which are in production today (Figure 1). The mine is located in the municipality of Gallivare, about 100 km north of the Arctic Circle (Figure 2).Successful sublevel caving requires progressive caving of the orebody after initial blasting. This results in caving of the hangingwall, sometimes all the way to surface, resulting in surface deformation. As mining progresses to deeper levels a larger area is affected by deformation. Over the years, extraction from the Malmberget orebodies has led to deformation of the ground surface, affecting residential areas and existing infrastructure.Forecasting the caving process is not straightforward as it can be rapid or slow, depending on the rock mass conditions, stress conditions, the failure mechanisms at work, etc. Moreover, several of the orebodies in Malmberget are non-daylighting, which makes reliable caving analysis even more difficult, primarily since there is less available experience from such cases. The present study was designed to improve the understanding of the caving and deformation mechanisms at Malmberget, and in particularly for the Printzskold orebody"

Jan 1, 2017

By J. Sjöberg, T. Savilahti, S. B. Umar

"Sublevel caving in Malmberget mine results in mining-induced surface deformation. One of the currently mined orebodies is the Printzsköld orebody. As mining deepens there is need to assess the behaviour of the cave formed in the subsurface above this orebody. Numerical analysis was used to assess the effects of extraction to deeper levels and perform strength parametric studies. Stress redistribution was studied, and the results showed high stress buildups in the hangingwall and the crown pillar. Two failure mechanisms have been identified – shear and tensile. Reducing cohesion by 50% increased the area of the yielded zone by more than 100% in the hangingwall.IntroductionThe Malmberget mine is owned and operated by the Luossavaara-Kiirunavaara AB (LKAB) mining company. Iron ore is extracted by large-scale sublevel caving. The mine comprises a total of 20 orebodies, about 10 of which are in production today (Figure 1). The mine is located in the municipality of Gallivare, about 100 km north of the Arctic Circle (Figure 2).Successful sublevel caving requires progressive caving of the orebody after initial blasting. This results in caving of the hangingwall, sometimes all the way to surface, resulting in surface deformation. As mining progresses to deeper levels a larger area is affected by deformation. Over the years, extraction from the Malmberget orebodies has led to deformation of the ground surface, affecting residential areas and existing infrastructure.Forecasting the caving process is not straightforward as it can be rapid or slow, depending on the rock mass conditions, stress conditions, the failure mechanisms at work, etc. Moreover, several of the orebodies in Malmberget are non-daylighting, which makes reliable caving analysis even more difficult, primarily since there is less available experience from such cases. The present study was designed to improve the understanding of the caving and deformation mechanisms at Malmberget, and in particularly for the Printzskold orebody"

Jan 1, 2017

By Amar Prakash, V. M. S. R. Murthy, Ajoy Kumar Singh, Avinash Paul, John Loui Porathur

Weak, layered, and fragile rock mass formation, if not supported properly, is subject to roof failure thereby affecting safety and productivity in underground coal mines. Though Central Mining Research Institute – Indian School of Mines Rock Mass Rating (CMRI-ISM RMR) based well-defined support design guidelines are established, still occurrence of roof failures in underground coal mines is a real matter of concern for mining engineers and researchers. Numerical modeling techniques are successfully used by several researchers by simulating rock mass condition for stability assessment of mine openings. The analysis becomes crucial in case of weak and fragile rock formation. The present research envelops the determination of 31 cases of rock load by CMRI-ISM RMR under different geo-mining conditions followed by the development of modified Rock Mass Rating (RMR), i.e., rock mass rating dynamic ( RMRdyn) by incorporating P-wave velocity as a new parameter. The rock load determined using RMRdyn and numerical models was correlated and found in close agreement. In addition, the deviation in rock load determined by all the three approaches, i.e., CMRI-ISM RMR, numerical modeling, and RMRdyn, was compared with the actual field data. The percentage deviation obtained in RMRdyn and numerical modeling is less compared to CMRI-ISM RMR.

Sep 1, 2023

By B J. Monaghan, P Zulli, X F. Dong

In the lower zone of blast furnace (BF), slag trickles down through a packed bed of carbonaceous particles in the form of films, rivulets or droplets. During its downward flow, there are significant interactions occurring between slag and other phases, ie gas and packing particles. In particular, the interaction between gas and slag is closely associated with the key flow phenomena such as loading, flooding or channelling in the lower zone of the BF. With gas introduced laterally through tuyeres and slag flowing downward from the cohesive zone to the hearth, this interaction varies throughout the lower zone and in the cohesive zone:  strong cross-flow of gas relative to the slag occurs in front of the raceway  counter-current flow of gas and slag occurs in the upper parts of the lower zone  both strong cross- and counter-current flow of gas and slag around and within the cohesive zone. These interactions can be linked to furnace irregularities, affecting smooth operation and limiting production. Therefore, understanding the gas-slag interaction and its influencing factors remains critical for BF process control. In the current review, the physical behaviour and key influencing factors of gas and slag flow are initially examined. Numerical models describing gas and slag flows in a packed bed at micro-, mesoand macro-scopic levels, are reviewed. Each simulation approach has its advantages and limitations, as well as unique assumptions when used in BF process modelling. Applications of these approaches to gas-slag flow in key identified functional zones will highlight the differences between approaches and modelling outcomes. As no single approach can fully address all aspects of gasslag interactions, an integration of different approaches is considered an effective and beneficial method in providing a more complete understanding of gas-slag flow in the BF process.

Jun 19, 2024

By N. Syna

The selection and preparation of activated carbon suitable for specific application can often involve exhaustive adsorbent preparation and testing. This work examined an emerging technique based on modelling the correlation of observed adsorption capacity to the carbon characteristics as a means of decoupling the roles of carbon textural and surface chemical properties. The insights to these specific roles are then used in tailor-making adsorbents and in establishing the optimal carbon properties for specific application. In this study, a series of activated carbons prepared from bagasse by physical and chemical techniques and commercially available carbons were used in developing the models. These carbons were tested for the adsorption of thiourea-gold (I) chloride complexes. Thiourea is an important alternative gold lixiviant with the potential to address current ecological and safety concerns in the extraction of gold. A four-component model, which correlated the adsorption gold capacity to carbon textural and surface chemical properties, was proposed. The model fit was assessed using multi-linear regression and a partial F-test. This was followed by goodness of fit test achieved through r2, adjusted?r2 and mean square error (MSE) values. The results suggest the adsorption of the gold-thiourea complex is strongly influenced by the internal surface area and pH of the carbon. Tailor-made activated carbon from bagasse showed adsorption gold capacities of up to 198 mgAu/g adsorbent, demonstrating the potential of bagasse as a precursor for gold adsorbents and of this empirical modelling approach in developing adsorbents for metal recovery.

Jan 1, 2014

By Pedro Lobarinhas

In the present work a two dimensional computational model is derived for the description of the time dependent heap leaching process. The flow inside the heap is calculated using the diffusion equations for the flow in porous media. The soil conductivity is also estimated assuming a close packed structure for the grain packing inside the heap. These equations, with appropriate boundary conditions, are solved using a time step splitting scheme. From this calculation the amount of leaching reagent in contact with the ore at any position inside the heap can be calculated and, once coupled with appropriate reaction equations, the amount of metals extracted can be computed and leaching agent concentration updated. The model is tested with available experimental data and the results have shown to be in a very good agreement. The model predicts all the transient features of the data.

Jan 1, 1999

By D Adhikary, M Khanal, B Poulsen, H Guo

During longwall mining, the ability to reliably estimate the height of longwall mining-induced fractures is critically important in order to effectively manage mining-induced water inflows and aquifer interference. The extent of mining-induced subsurface deformation, fracture, surface subsidence and aquifer interference is largely controlled by local lithological conditions, mining methods and mine layouts. Therefore for the purpose of quantifying mining impacts it is imperative to incorporate the local geology and mine plan parameters in a mine scale model. Two different numerical models based on 3D continuum (Finite Element Method) and 2D discrete (Discrete element method) domains, respectively, are developed to incorporate the mine site parameters of a study site in Australia. Both the models are calibrated and validated with measured mine data. The continuum numerical modelling method is used to investigate the change in permeability of the overburden strata due to underground extraction considering the full 3D nature of geology and mine layout. The discrete element modelling method is used to simulate the initiation and propagation of fractures, allowing the caving of overburden strata to be represented in a realistic manner, and the permeability enhancement of the strata overlying the goaf to be calculated from first principles. Permeability enhancement in the discrete element model is a direct function of the fractures induced in the deformation/caving process; fracture location, extent, aperture and connectedness. From the numerical modelling under the simulated condition for the mine studied here, it is predicted that the permeability of the caving zone from the goaf to five times the mining height increases by almost six orders of magnitude. In the fracture zone to approximately ten times the mining height, the change in permeability is predicted to be around two to three orders of magnitude. At the height of more than ten times the mining height the magnitude of permeability increased is predicted to be one to two orders of magnitude. CITATION:Khanal, M, Poulsen, B, Adhikary, D and Guo, H, 2018. Modelling of height of connective fractures during longwall mining, in Proceedings The Fourth Australasian Ground Control in Mining Conference (AusRock), pp 254–265 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 30, 2018

By B. Csoke

To improve the efficiency of the crushing/screening plant of Bazalt-Középko Kobányák Ltd. (Hungary), the implementation of computer aided quality control is planned. In the current phase, the R & D work aims at the solution of high-priority tasks: 1) Developing the process models, validation of the models by means of industrial measurements. 2) Developing the algorithm for the process control of the crushing/sizing plant. 3) Developing the suitable computer program. To develop the mathematical models of crushing (breakage and probability functions) and to determine the model parameters, the experimental industrial trials were carried out on hydrocone crushers of types Svedala H 4800 EF and H 4800 F, changing the gap size and crushing relatively narrow grain size fractions. Additionally, the parameters of the separation function for the screens operating in the system were determined. It was established that the crushing probability depends on the grain size of the feed (Y) and the gap size (s), so the breakage probability S, for each cone crusher, can be well characterised by the function formed using the ratio of these two parameters: SfYsf==. The distribution of the broken grains can be described reliably using the relative grain size distribution function B(? = x/Y) (usually Schuhmann-Gaudin- or Rosin-Rammler-shaped function). For the screens, the probability function ()2()111/NTx = xd ???---???? is suitable for cha-racterising the separation (where x? grain size, ?-open-area ratio, d ? rectangular screen aperture size, N ? number of appearance of the particle on a screen surface). Using these functions, the comminution circuits were described by means of vector/matrix equations. The calculations and measurements verified the reliability of the developed models and algorithms.

Jan 1, 2014

By T Krull, T J. Donohue, D Ilic

"Problematic issues in materials handling applications are often amplified in the iron ore industry due to the relatively high density of the ore as well as the abrasive nature of the product. One of the main issues arising from the handling of iron ore is excessive wear, with two applications typically suffering accelerated wear issues being transfer chutes and conveyor belts. These wear issues are problematic due to additional capital costs, but mainly due to the associated downtimes resulting in the accelerated wear. This study will focus on the modelling of iron ore in practical applications where design work was undertaken in order to reduce the wear rates of different materials handling equipment.CITATION:Donohue, T J, Ilic, D and Krull, T, 2013. Modelling of iron ore in materials handling applications with a focus on wear, in Proceedings Iron Ore 2013, pp 355-362 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Aug 12, 2013

By M. W. Bates, A. L. Mular

This paper, based on work carried out at Queen's University, describes the development and use of mathematical models of the parallel cyclone classifiers being used in the Strathcona mill of Falconbridge Nickel Min.es Ltd.

Jan 1, 1971

By A D. Campbell, G R. Power

The increasing availability of marker data in sublevel caving operations is changing the way draw can be managed to achieve positive recovery outcomes in these mines. This paper describes how marker data has been incorporated into flow simulations to help calibrate full mine ore recovery simulations for the forecasting and optimisation of recovery.Some flow modelling techniques do not use marker data, relying instead on the outcomes of physical and numerical experiments to guide assumptions about an idealised recovery flow shape. Other mines have access to marker trial data, but find it difficult to specifically link this marker trial information to their recovery models. This can make it difficult to create a link between the marker trial information and recovery forecasting and limits the potential value of the marker trials at the mine.In the simulations shown here, the marker rings are fully incorporated into the flow model, and can be tracked in the model. The model is compared against the real underground marker recovery information to give the user an indication of how accurately the flow settings in the model are reflecting the conditions in the experimental marker rings. An example of use of this system at an operational mine is included.CITATION:Power, G R and Campbell, A D, 2016. Modelling of real-time marker data to improve operational recovery in sublevel caving mines, in Proceedings Seventh International Conference and Exhibition on Mass Mining (MassMin 2016), pp 105–110 (The Australasian Institute of Mining and Metallurgy: Melbourne).

May 9, 2016

By J. Yianatos

The regrinding stage is necessary to achieve the particle size (liberation) for the final upgrade in the cleaning flotation process. At present, vertical mills (Vertimill) are employed at industrial scale, using small steel balls as grinding media. The regrinding process characterization requires mixing regime knowledge along with the mean residence time of the liquid and solid phases. The residence time distribution (RTD) and the actual mean residence time estimations from theoretical relationships are complex due to the presence of three components (solid mineral, balls and water), mass transport characteristics and segregation affecting solid holdup. Therefore, the experimental RTD determination is suitable for industrial mills. In this work, RTD measurements were performed in two different industrial flotation plants. Both plants operate with Vertimills of 1000-1500 HP in the regrinding process with different circuit arrangements according to the plant treatment capacity. The RTD was measured using the radioactive tracer technique. This technique allows for non-invasive tracer detection where the procedure consists of selecting a liquid (Br82 in aqueous solution) or solid (irradiated solid) tracer for RTD data acquisition in real time. The RTD was modelled by the Large and Small Tanks in Series (LSTS) model, which allowed the characterization of the different mixing regimes for both liquid and solids in the regrinding mills. Effective mean residence times in the range 2-11 min were obtained with RTD shapes similar to those obtained in tumbling ball mills. It was found that residence times of liquid and fine particles were similar, which is consistent with results observed in other mixing systems, e.g. mechanical flotation cells. Also, the mean residence times estimated from the RTD allowed the mass balance adjustment to be validated from metallurgical sampling data. RTD estimations are useful for better understanding the Vertimill behaviour and to identify the actual breakage rates. Keywords: vertical mill, residence time distribution, particle size, modeling

Sep 1, 2012

By R. Nilsson

In the present work, an assessment of the sulphide capacities of a number of binary and ternary slag systems have been carried out in the light of a sulphide capacity model. The model is based on the expression: {CS- exp (-A G °/R T) ? ao2 / fs2-} where AGO is the standard Gibbs energy change for the above reaction, ao2 is the activity of oxygen ions and f, is the activity coefficient of the sulphide ion in the slag. In the present model, the ratio (ao2/fs2) is optimized for the experimental data available in literature. The systems assessed are the binaries MnO-SiO2, CaO-SiO2, Al2O3-CaO, as well as the ternaries CaO-Al2O3-SiO2 and CaO-MnO-SiO2. In the case of all these systems, the experimental data available in literature were complemented by measurements carried out at the Division of Theoretical Metallurgy. The sulphide capacities of some FeO- containing binary, ternary and quaternary slags are also examined. The potentiality of the model to estimate the sulphide capacities of complex slag systems from binary and ternary data is demonstrated.

Jan 1, 1996

By M. R. Khalesi, S. Bellec, C. Bazin, D. Hodouin

Batch grinding tests conducted on a gold ore indicated that gold is not randomly distributed within the size fractions of the ground ore. Since gold leaching performances were found to be dependent of the size of the gold carrying particles an integrated gold ore grinding-leaching simulator should incorporate a method to predict the gold size distribution of a ground ore. The paper describes the results obtained from batch grinding tests showing the non random distribution of gold within the size intervals of a ground ore and the initial developments of a grinding model to predict the gold size distribution in a batch grinding product.

Jan 1, 2011

By Q. H. Liao, B. G. Monaghan

The options for designing slopes in complex geological conditions are reviewed. It is shown that modelling, especially numerical modelling, allows a degree of analysis that cannot otherwise be achieved. Examples are presented of how numerical modelling can be used to investigate fundamental aspects of rock mass behaviour. It is demonstrated that an important consequence of change in scale is that the mechanisms of failure also change. The complex failure of the footwall slope at Aznalcollar, a large open-pit mine in southern Spain, is discussed. Numerical modelling with UDEC (universal distinct-element code) has allowed the realistic simulation of displacements measured in the field

Jun 18, 1905

By Markus A. Reuter, Yongxiang Yang, James R. Post

The packed coke bed in the lower part of an iron making blast furnace, the dead man, dissolves slowly in a pool of hot metal, which is tapped regularly. The coke particles of the dead man not only exhibit a size distribution, but also their size distribution is non-uniform throughout the dead man and as a result, the porosity of the dead man is non-uniform. This causes the flow of hot metal through the dead man to be irregular and non-uniform, bringing about channelling in the coke bed. Channelling is detrimental to the refractory lining life of the hearth of a blast furnace. Using population balance modelling and a porosity model it is possible to keep track of the local changes in size distribution and bed porosity during dissolution of coke particles. Numerical simulations show that carbon dissolution from a packed coke bed (with particles which exhibit a size distribution) increases the average porosity, decreases the mean particle size and changes the size distribution to a distribution with a smaller size range. Influence of the type of initial size distribution (mean, variance) is investigated. The calculated changes in porosity are used to dynamically adapt the flow and the mass transfer (kinetics) during calculations. This approach provides the possibility to implement dynamic porosity calculations for CFD applications, and thus simulate the results of various coke types and particle size distributions on dissolution and flow behaviour.

Jan 1, 2003

By J. Gebhardt, D. McBride, C. R. Bennett, Croft. N., M. Cross

A comprehensive model of solution flow through heap leach systems requires a numerical solution with the ability to predict variably saturated-unsaturated flow through porous media domains that contain materials with spatially variable properties. The types of flow associated with stockpile leach processes lead to flow problems, such as infiltration into dry soil, drainage, perched water tables and flow through heterogeneous materials. Computational methods for modelling this class of flow behaviour are conventionally based on the classical Richards' equation. However, the governing equations are highly non-linear, difficult to solve, and require iterative numerical solution methods. The pressure-based form of the Richards' equation suffers from poor mass balance, while the mixed form can possess convergence difficulties. An adaptive transformed mixed algorithm is described, which reduces the non-linearity of the problem, optimizes the time step size, and provides a fast, numerically robust scheme that significantly reduces computation (or CPU) time. This method is shown to give fast, accurate solutions on a number of complex flow problems. The utility of this algorithm is illustrated through its implementation within the PHYSICA computational modelling software, which incidentally, provides a framework for a comprehensive heap leach model.

Jan 1, 2005