Search Documents

By A. Della Rocca, P. Jonker, B. Henning

Successful prediction of the collection efficiency of cyclone separators is a necessary step in cyclone design and verification. Aside from standard empirical relations, CFD modelling is widely adopted to supplement the design of cyclone units in relation to the desired separation performance. Nevertheless, effective adoption of numerical modelling is sometimes prevented by large modelling errors, especially for small size particles which are more subject to turbulence effects than larger particle classes. In the past, this motivated the adoption of more refined turbulence modelling approaches, such as large-eddy simulations. These approaches provide greater insight into time dependent behaviour, while leading to larger computational cost with respect to more standard approaches. In this paper, an effective Eulerian-Lagrangian Reynolds-averaged Navier-Stokes (RANS) approach is investigated, which is able to predict the complete separation efficiency curve with satisfactory accuracy, without requiring prohibitively long computational times. An assessment of the proposed approach is carried out against experimental data from large laboratory scale cyclones, in the form of mean and fluctuating velocity profiles from LDA measurements (Hoekstra, 2000), as well as with respect to the experimental separation curve (Zhao, 2005). The final application to a full-scale industrial cyclone unit is then carried out with the aim of tuning the cyclone design to the separation efficiency curve as required by downstream processes. Keywords: Cyclones, separation, CFD, RANS, RSM, Lagrangian

Jan 1, 2020

By S. Purushotham, Y. R. Murthy, M. Narasimha, S. K. Tripathy, M. Kumar

Spiral concentrators are gravity-based separators widely used to beneficiate coal, heavy minerals and beach sand. It consists of a helicoid trough that spirals downward in the helical configuration around an axial column. The particle segregation takes place based on size and density under the influence of gravity as the particles spiral down within the thin water layer. In the last century, most design modifications of this equipment were based on trial and error based lab-scale experimentations and empirical approach. The use of Computational Fluid Dynamics (CFD) has gained interest in the past two decades to understand the underlying flow physics in the spiral concentrators. In this paper, we used CFD to model the fluid flow and investigate the particle segregation inside a conventional LD9 spiral concentrator. Firstly, the Volume of Fluid model was used to resolve the liquid-air interface. The liquid depth and flow field predictions for varying flowrates were validated against the experimental literature data. One-way coupled VOF-LPT model was then used to describe the particle dynamics in the spiral flow field assuming dilute particulate flow. Moderate dense particulate flow (20 vol.%) was then modelled using the Algebraic Slip Mixture (ASM) model. Four particle-phases representing particles of different sizes of dual-density were simulated, and the component particle segregation was analysed. Keywords: Spiral concentrator, two-phase flow, multiphase, discrete phase model, ASM model

Jan 1, 2020

By R. Kuyvenhoven, A. Parbhakar-Fox, X. Prieto, J. Quevedo, L. V. Danyushevsky

This paper presents the results of a study focussed on determining the distribution of cobalt in different types of pyrite of the La Estrella IOCG deposit in Chile. In the IOCG mineral deposits, where copper is the main product, the quantification and distribution of cobalt is required to correctly define the technical and economic viability of recovering cobalt as a by-product. In the La Estrella mineral deposit, seven alteration and mineralisation events were recognised through both optical and automated electronic mineralogical analysis. Of these seven events, five different types of pyrites were identified, each with unique characteristics of the pyrite crystals. Through advanced elemental analysis of these five types of pyrites, it was determined that the concentration and distribution of cobalt within the pyrite crystal is dependent on the mineralisation event. The results indicate high cobalt concentrations in the early pyrite mineralisation events, associated with relatively fine pyrite grains, and cobalt generally occurring in a solid form in the pyrite core. The late pyrite mineralisation event exhibits the lowest concentrations of cobalt, with cobalt predominantly present along the outer periphery of the pyrite. These results are relevant from a metallurgical point of view because the preferential location of cobalt within the pyrite crystal will influence the maximum metallurgical recovery of cobalt. For this, a cobaltiferous pyrite concentrate must be formed which, according to the characteristics of pyrite, must be processed with the optimal solubilisation method. The overall objective is to ensure both good cobalt recovery and a high concentration of cobalt in the leach solution resulting from pyrite leaching. Keywords: Cobalt, pyrite, IOCG, copper, metallurgical recovery, solubilisation, Chile

Jan 1, 2020

By S. Shi, X. Xia, F. Xuesai, Z. Shen, D. Han

The hydrodynamic characteristics of a mechanical flotation cell’s transportation zone have significant impact on its metallurgical performance, especially on coarse particles recovery efficiency. However, how to define and adjust the transportation zone was seldom investigated by researchers. In this paper a method will be introduced to define the boundary of the transportation zone by PIV measurement. It is found that the transportation zone can be extended upwards by enlarging the entire rotor-stator assembly from the bottom of tank and increasing the jet angle of the impeller. Two control methods have been validated in industry tests on KYF-160m3 and KYF-320m3 air-forced mechanical flotation cells. The result shows coarse particles recovery are both improved distinctly. It is suggested that enlarging the top transport boundary is able to helpful to reduce the floating distance of bubbleparticle aggregates, which contribute to improve coarse particle recovery in large flotation cells. Key Words: Flotation cell; Transportation zone; Coarse particles flotation

Jan 1, 2020

By X. Yang, G. D. Zhai

The prop is the most important part to ensure the proper functioning of hydraulic support. When the hydraulic support is impacted by the roof, the prop is prone to extrusion deformation, expansion, and even bursting. In order to study the stress on a doubly-telescopic hydraulic support prop under impact load, an impact simulation was carried out based on the drop weight method. First, the impact model of the hydraulic system of the prop was established in AMESim software and the dynamic response curves of the bottom and middle cylinder of the prop obtained. Then, according to the conservation of energy, the pressure formulaes in the hydraulic cylinder of the prop under impact load were derived and verified by the AMESim simulation results, and the maximum pressure on the inner surface of the hydraulic cylinder of the prop obtained. Lastly, the transient dynamic simulation of a hydraulic support prop was carried out in ANSYS Workbench software, and cloud diagrams of stress and deformation of the prop obtained using Workbench simulation. The AMESim simulation process describes the change in fluid pressure in the prop, and the pressure formulae can be used to estimate the internal pressure of the prop under impact load. The finite element analysis results show that the stress of the middle hydraulic cylinder is much greater than that on the bottom hydraulic cylinder under impact load, which can provide a reference for prop design.

Jul 2, 2022

By E. Mokgwamme, T. Goff, M. Mvalelwa, P. Spratt, P. Bezuidenhout, P. Nkosi, H. Joubert, G. de Villiers

Sibanye-Stillwater recently completed a partial sidewall rebuild on Furnace 1 at its Marikana smelter complex. The project involved matte taphole maintenance, which required the removal and subsequent replacement of 18 composite copper-graphite coolers and a section of the lower refractory sidewall. These composite copper-graphite coolers are a newly developed cooling design by Tenova Pyromet, and this rebuild project provided valuable insight into the performance of the composite coolers. The novel composite cooler design was implemented to serve as a safeguard against localised wear and corrosion in the furnace sidewall, particularly at the slag-concentrate interface, where chloride-accelerated sulphidation has been identified as a dominant degradation mechanism affecting the copper cooling elements. The composite cooler design is unique in that all the exposed sides of the water-cooled copper elements are fully encapsulated in graphite, which protects the copper from corrosive sulphur- and chlorinebearing gases and condensates. This sidewall cooling system was designed to improve corrosion resistance, stabilise freeze lining formation, and extend the campaign life of the furnace sidewall. The partial sidewall rebuild was executed during a planned maintenance shutdown with the main objective of replacing the matte taphole lintel bricks. This paper serves as a case study that compares and evaluates the operational performance of the composite coolers before and after the rebuild. The installation method is discussed to highlight the approach taken to ensure a sound installation while minimising downtime, as well as to present some lessons learned in the process.

Mar 6, 2026

By S. Annavarapu

"Fragment size distributions play an important part in the design and planning of block cave mining operations. Though several methods have been proposed for estimating fragment size distributions in advance of mining, the calibration of these estimating procedures to field measurements has been a challenge. The Block Size Estimator (BSE) program was developed for providing an assessment of primary fragmentation expected in a block cave mine using drill-core piece length data and joint characteristics from the exploration and geotechnical evaluation programmes. In order to validate the results from the BSE program, the drill-core piece length data and the joint characteristics in different rock types from the DOZ block cave mine in Indonesia were used to help generate fragment size distributions in different areas of the mine. The predicted fragment size distributions were compared with the fragmentation observed at the drawpoints during operation of the mine.This paper presents the details of the validation of the fragment size estimates by the BSE program using fragment size distribution data from the drawpoints at the extraction level at the DOZ block cave mine, and the challenges encountered in developing reasonable correlations between the estimated and measured fragment size distributions. IntroductionThe assessment of fragmentation is an important aspect of the design and planning of the block cave mining method, and forms an integral part of the design of the excavations at the extraction level and the selection of material handling systems for transporting the ore to the processing stations. Secondary blasting requirements can also be estimated based on the fragment size distributions developed for the block cave (Laubscher, 1994).Several methods of estimating fragment size distributions in block cave mines have been developed based on joint set parameters estimated from structural data from oriented core drilling campaigns and mapping of available excavations and outcrops. The BCF program (Esterhuizen, 1999) estimates the fragmentation using a combination of empirical, analytical, and rational methods to model the behaviour of materials during the primary and secondary fragmentation processes. The program JKFrag uses advanced tessellation of joint traces to create rock blocks and then generate fragment size distributions (Brown, 2002). Nickson, Coulson and Hussey (2000) presented the details of processing of the geotechnical data from the Mont Porphyre project, including an assessment of block sizes from the data collected on the core logging sheet. The distance between natural fractures was used to create distributions of block lengths, and a distribution of potential block sizes was estimated assuming that the blocks had equal dimensions in all directions. The results agreed closely with the results of the BCF program. Hadjigeorgiou, Grennon, and Nickson (2002) demonstrated the use of oriented borehole data to provide characteristic block size distributions for a mining project using a software package called Stereoblock, which generates a three-dimensional joint network for a given volume and calculates the volumes of blocks created by the intersections of these joints, simulated as circular planes. The distributions of joint orientation and joint spacing generated from statistical analysis of the oriented core data are used for the creation of the Stereoblock model. The diameter of the circular planes representing the joints is determined based on the mean and standard deviation of the joint trace length obtained from scan-line mapping."

May 1, 2019

By D. Orynbassar, N. Madani

This work addresses the problem of geostatistical simulation of cross-correlated variables by factorization approaches in the case when the sampling pattern is unequal. A solution is presented, based on a Co-Gibbs sampler algorithm, by which the missing values can be imputed. In this algorithm, a heterotopic simple cokriging approach is introduced to take into account the cross-dependency of the undersampled variable with the secondary variable that is more available over the entire region. A real gold deposit is employed to test the algorithm. The imputation results are compared with other Gibbs sampler techniques for which simple cokriging and simple kriging are used. The results show that heterotopic simple cokriging outperforms the other two techniques. The imputed values are then employed for the purpose of resource estimation by using principal component analysis (PCA) as a factorization technique, and the output compared with traditional factorization approaches where the heterotopic part of the data is removed. Comparison of the results of these two techniques shows that the latter leads to substantial losses of important information in the case of an unequal sampling pattern, while the former is capable of reproducing better recovery functions.

Aug 1, 2021

By A. J. S. Spearing, M. Sharifzadeh, K. V. Jessu

Empirical pillar design methods are commonly used in the mining industry. The parameters within which these methods are valid are frequently unknown to the user or ignored. In addition, empirical design may not consider essential parameters such as blasting effects, orebody dip and the presence of geological structures, which all adversely affect the stability of the pillars. This can result in potentially serious pillar design strength over-estimates. Although the commonly based tributary area method is generally conservative, as the spans are seldom that large relative to the depth, failing to consider other relevant parameters can result in errors. Problems associated with an under-designed pillar can range from a local pillar collapse to a catastrophic chain reaction collapse (or run). Over-designed pillars are generally safe but reduce the extraction of the orebody, thus adversely affecting the profitability of the mining operation. We used laboratory tests and numerical modelling to understand the effects of pillar orientation, blasting and the presence of discontinuities on pillar strength. Reduction factors were developed with these models to be implemented in conjunction with the existing empirical pillar design methods. For any pillar or mine design, once it is implemented, the actual performance of the system must be checked regularly by observation and monitoring and adjusted if needed. The pillar design approach outlined in this paper can better optimize the pillar mining method by considering other generally ignored but important parameters, thus improving safety, productivity, and economic aspects.

Oct 1, 2022

By T. Nghipulile, M. Moys, M. Bwalya, H. Simonsen

The quest for efficiency in comminution is an ongoing concern as comminution usually constitutes a major cost component in the metal production industry. Such improvements can be made by circuit optimization or development of more efficient equipment. A novel crusher, known as the rotary offset crusher (ROC), promises to deliver in this space. The ROC was invented in 2002 by Michael Hunt, Henry Simonsen, and Ian Sinclair, but failed to garner enough support to progress to production. The original design concept was recently rekindled, and a laboratory version of the crusher has been built and commissioned at the University of the Witwatersrand. The crusher is simple in design, with two cylindrical discs that are parallel to each other, and, as the name implies, there is an offset between the vertical axes of the discs. The top disc has a conical section on its lower face, and this creates a crushing zone between the opposite faces of the two spinning discs. Centrifugal motion transports particles through the crushing zone. Batch experiments have been conducted with quartz at various crusher settings (discs offset, rotational speed, and vertical exit gap) for various feed size distributions. The indications so far suggest that the disc speed is a key factor affecting the performance. Size reduction ratios as high as 11 were recorded from experiments with quartz at a speed of 830 r/min.

Jan 1, 2021

By D. F. Malan, T. E. Oates

A recent experimental pillar extraction project at a UG2 bord-and-pillar mine presented a unique opportunity to compile a new pillar database. Currently, the South African hard rock bordand- pillar mines are designed using the Hedley and Grant formula with a modified K-value. This empirically derived formula was developed for uranium mines in the Elliot Lake district of Canada. The use of this formula for the design of pillars in South Africa is questionable. Very few pillar failures have nevertheless been observed and its current calibrations for the various reef types are possibly too conservative. A new UG2 pillar database of 66 pillars, of which seven are classified as failed, was compiled by the authors. This enabled a revised ‘first-order’ calibration of the K-value for the Hedley and Grant formula. The new estimated value for the UG2 is K = 75 MPa. This gives a pillar strength that is more conservative than the PlatMine formula. This work should nevertheless be considered as only a preliminary calibration as the database was small. Further work is also required to determine whether the exponents in the formula for the width and height parameters are appropriate for UG2 pillars.

May 3, 2023

By J. W. van der Merwe, K. G. Moloko

Furnace wall tubes from 600 MW subcritical boilers at three coal-fired power stations were sampled and the fireside deposits examined to determine the mechanism of fireside corrosion. This involved an indepth investigation into the morphology and composition of the fireside deposits and the conditions of the furnace that enable this type of attack. SEM-EDS analysis revealed high concentrations of oxygen, iron, and sulphur, QEMSCAN and XRD analyses identified the presence of Fe3O4, Fe2O3, FeS, and FeS2. Differential thermal analysis showed thermal activities at temperatures of 500–600°C, 900–1100°C, and 1100–1250°C, which are associated, respectively, with FeS2 oxidation to FeS and Fe2O3, at 475–525°C, formation of aluminosilicates at 925–1100°C, and melting of FeS around 1190°C. The absence of sodium and potassium eliminates the contribution of molten alkali sulphates to the corrosion. The consistent coexistence of iron sulphide and iron oxide is indicative of the substoichiometric conditions in the furnace, while the detection of pyrite suggests that the coal is not completely combusted, which points to a poor combustion process. These observations were affirmed by gas analysis at one of the stations, where very high levels of carbon monoxide were measured at the furnace wall (≥ 14 000 ppm) and furnace exit (≥ 3500 ppm). The high CO concentrations are indicative of limited combustion caused by limited O2. These reducing conditions promote the formation of FeS-rich deposit, which is the corrosive species responsible for degradation.

Jun 1, 2021

By J. P. Gouvea, W. Geel

Resin bolting in South Africa has been developed and refined in underground coal mines (i.e., soft rock) since the late 1940s. The extension of resin bolting to hard rock mines started in 2003 and has required significant adaptations, challenging the status quo from both a supplier and applicator perspective. As far as resin bolting is concerned, the most affected aspect in coal mining and hard rock mining environments, is the rock drills or machinery used for installation. Hard rock mines mostly use hand-held airlegs with limited thrust, low torque, and rotation speed, relative to coal mine mechanised bolters. Resin bolt development for hard rock mines was, and still is today, predominantly focussed on compensating for the ‘shortfalls’ of using hand-held airlegs. At face value, hydro-powered rock drills, which are more powerful than traditional handheld equipment, provided hard rock mines with a means to attain more consistent and reliable resin bolt installations. This paper assesses the quality and performance of hydro-powered resin bolting at an intermediate to deep level platinum mine in South Africa and its potential contribution to rock-related instabilities when applied erroneously. The findings of this paper necessitate a revision of current quality control and quality assurance measures related to resin bolting for underground support.

Nov 24, 2025

By E. Córdova, A. Anani, W. Nyaaba

Optimization of panels, haulage fleet, and waiting area involves deterministic and low-fidelity methods and experiential knowledge. The process is challenging because coal recovery and operational capabilities must be considered in the solution. The approach in this manuscript comprises the development of an integrated stochastic simulation model of a coal room-and-pillar system that addresses these challenges. The decision variables evaluated are panel width, number of shuttle cars, and change-out time (COT). The results show that the mine should implement the shortest possible COT, decreasing the cycle time and thereby increasing productivity and continuous miner (CM) utilization. The highest productivity and CM utilization for a fleet size of three shuttle cars is found in the 15-entry panel width. For the evaluated fleet sizes, the 19-entry panel width is optimal for the four and five shuttle cars. Among the three variables studied, panel width and fleet size had the most significant effects (5% increase) on the CM productivity, cycle time, and utilization.

Apr 2, 2022

By M. Royle, A. Bloem, R. Uken

Rock mass damage induced by mining increases the fracture frequency and joint aperture width in the pit walls and rock mass behind open pit mine slopes. The depth and intensity of this mechanical damage relative to the pre-mining state is not well understood and yet is a major factor in the development of slope stability models and operational guidelines. The induced damage can be partitioned into domains or zones. The blast damage zone is essentially a ‘free flow’ rock mass volume, which transitions quite rapidly into the excavation damage zone, which extends beyond that. In both zones, transient high in situ water pressures are often linked to unstable slope conditions. Whereas rock damage can enhance passive depressurisation in the rock mass, high infiltration rates on the other hand, as one could expect during rainstorm events, will decrease slope stability within the dilated pit rock mass, especially in the blast damage zone-excavation damage zone, by increasing pore pressures for short periods of time. Ice jacking and freshet are additional factors to consider in (sub) arctic environments and can be exacerbated by unfavorable rock fabric orientation. On a local scale, brittle deformation damage zones, in the form of faults, can contribute to slope instability, or enhanced stability because of slope depressurisation, depending on their location and orientation relative to the pit shape and the water sources. In this study, the shapes of the damage zones and their anticipated impact on infiltration and/or passive depressurisation are considered for different pit morphologies, rock mass conditions, in situ stress, and hydrogeological settings. Some unusual hydrological situations, like mining towards rivers and lakes and the potential effects on pore pressures, are also considered. This investigation has resulted in the development of a method to adjust initial hydraulic conductivity values used to estimate transient pore pressure conditions found in rock masses exposed in different pit morphologies and stress regimes as well as discontinuities’ orientation relative to open faces within the open pit mine. The changes in K values and pore pressures can then be used to determine hydrogeological influence more accurately on future pit slope stability.

Sep 1, 2025