Search Documents

By A. V. Rose

A redevelopment in the City of London has been constructed with a piled foundation system that is believed to be the first of its kind. The purpose of the research is to explore this novel foundation method, which involves the installation of closely spaced piles in a perimeter arrangement with a body of soil enclosed in the centre. Both centrifuge modelling and numerical modelling techniques are used.

Jan 1, 2011

A model has been developed as part of the AMIRA P420 Gold Processing Technology Project to simulate the recovery of gold by centrifugal gravity devices. The model (which is available online) uses the simulation methodology developed by Laplante et al (1995) to predict gravity recoverable gold (GRG) recovery under varying operating conditions and circuit configurations. In addition to predicting GRG recovery, the model also predicts the GRG content and size distribution of the cyclone overflow and underflow. The predicted overflow GRG content and size distribution can be used as inputs for a leaching model in which size dependent leaching rates are applied for each size class of GRG and non-gravity recoverable gold (non-GRG) in the cyclone overflow, allowing size-by-size leach recoveries to be determined for GRG and non-GRG. The application of the gravity model and leaching model are demonstrated using a case study. The simulated gravity recovery and leaching performance showed excellent agreement with the measured values. The developed models can be used to assess the sensitivity of gravity gold and leach recoveries to varying feed GRG content, GRG size distribution, gravity circuit feed rate, classifier performance and leach residence time. The described methods can be used to optimise new and existing gravity and leaching circuits.

Aug 1, 2010

By V. G. Papangelakis, M. Cheung, M. Dry

The high pressure acid leaching of laterite was modelled and simulated using a combination of two commercially available software packages, namely Aspen Plus 11.1 and the OLI Stream Analyzer. New kinetic parameters for the acid-controlled metal dissolution and precipitation reactions at temperature were first determined from previous experimental data on the leaching of limonite. A grain model with cylindrical pellets was used to model the solid-liquid dissolution reactions. Precipitation reactions were modeled via aqueous equilibrium chemistry. A dedicated thermodynamic database previously developed at the University of Toronto was used in conjunction with the OLI software to model the solution chemistry and calculate the concentration of hydrogen ions in solution at the leaching temperature. Mass balances were handled by Aspen Plus. The simulation gave good predictions of nickel dissolution under industrial and laboratory conditions, for the dissolution of limonite and limonite/saprolite blends. Computational

Jan 1, 2005

By P Holtham, R Dunne, D Stoll

Flash flotation is a process designed to remove coarse, liberated, fast floating minerals from hydrocyclone underflow. Irrespective of liberation, coarse particles would normally remain in a grinding circuit until they are reground fine enough to report to the hydrocyclone overflow. Additionally, minerals of high density also report to hydrocyclone underflow despite being fine in size and/or liberated because hydrocyclones separate on the characteristics of both particle size and density. Dense particles continually reporting to the hydrocyclone underflow can lead to over-grinding and the production of ultra fine particles. These particles exhibit slow flotation kinetics and thus are difficult to float in conventional flotation circuits. Therefore, removing coarse, liberated particles from the hydrocylcone underflow by flash flotation has the potential to improve recovery and can improve dewatering unit operations and reduce grinding costs.The impact on the overall performance of the flotation circuit by removing coarse, fast floating minerals prior to conventional flotation is not well understood. Ultimately, predicting the overall recovery improvement provides the justification for the installation of a flash flotation cell within the milling circuit. This paper will describe attempts to estimate flash flotation recovery improvements to the overall circuit using kinetic data from simplified batch cell flotation tests combined with limn-based modelling techniques.CITATION:Stoll, D, Dunne, R and Holtham, P, 2012. Modelling the impact of flash flotation from batch cell flotation tests using limn, in Proceedings 11th AusIMM Mill OperatorsÆ Conference, pp 179-184 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Oct 29, 2012

By S S. Kanchibotla, I Brunton

Most blast fragmentation models assume the rock mass properties, explosive properties and blast design variables to be constants and uniformly distributed within a blast. However, in reality all these input variables vary within a blast resulting in variation in the resulting fragmentation size distribution. A stochastic modelling approach is introduced in this paper to quantify this variation. This technique takes the input variables as statistical distributions rather than constants and through several thousand iterations, generates a statistical representation of the expected fragmentation resulting from a production blast. A case study of three production blasts from a large open pit mine are presented and the modelled fragmentation æenvelopeÆ shows good agreement with the fragmentation æenvelopeÆ estimated from Split image analysis. The various blast-related parameters influence different parts of the fragmentation distribution, eg rock strength and explosive velocity of detonation have most impact on the fines. The technique is used to identify the parameters that have the greatest influence on various size fractions. Such an analysis will be useful to direct resources to efficiently minimise the variation.

Jan 1, 2001

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 R Wolff

One critical factor in open pit mining projects is the estimation of the recoverable reserves. The reason for this is that at the valuation stage there is a lack of information about future metal prices and production costs. Consequently, the estimation of the recoverable reserves needs to be done based on a fixed break-even cut-off grade (BECG), in which fixed expected metal prices and production costs are assumed throughout the operating life of the mine (OLM). In this paper, an alternative technique for estimating the recoverable reserves of an open pit mining project is presented and explained in detail. The main characteristic of this technique is that it estimates the recoverable reserves of the project using mining costs as reference: this is done due to the fact that mining costs are easier to model than metal prices. Another important characteristic of this technique is that it uses the uncertainty of the geology of the orebody, given by simulations of the deposit, to make final strategic decisions, such as the selection of the operating life of the mine and the optimal production rate, which minimise the risk of not achieving future production targets while maximising the upside potential of future rewards. As it will be shown, the final outcome of the proposed technique is the generation of a robust open pit design, that include cutbacks and ultimate pit limits, called the marginal open pit mine design (MOPMD). One characteristic of the MOPMD is that, at each production period, it will be characterised by their respective geological uncertainty, given as probability distributions of ore, waste and metal quantities, and by marginal economical project indicators, such as cut-off grades and metal prices, among others. The details follow in the paper.

Jan 1, 2006

By Y. Houbaert, L. Kestens, R. Decocker

"A magneto-elastic model is presented to calculate the orientation dependence of the magnetostrictive strain, observed at saturated magnetisation in ferromagnetic materials with a cubic crystal structure and an arbitrary crystallographic texture. The formula of Becker and Doring is used to express the anisotropy of magnetostriction for a single crystal. In order to simulate the macroscopic average magnetostriction of a polycrystalline aggregate (with an arbitrary texture) the Reuss assumption of the elasticity theory was applied. According to this assumption the various orientations of the polycrystal can defonn without constraints, producing local strain incompatibilities in the microstructure but observing total stress equilibrium. The macroscopic strain is calculated as the weighted average of the individual strains of all orientations composing the polycrystal. The weight factors are determined by the volume fractions of the corresponding orientations of the given texture, which can be measured by standard X-ray diffraction techniques. The model is applied to simulate the variation of magnetostriction (at saturation) with respect to the rolling direction for a standard grade of Goss-oriented electrical steel. A brief comparison with experimental data allows validating the basic model assumptions. IntroductionBecause of tougher environmental constraints, the noise reduction in transformer cores has become increasingly important. It is generally accepted that the main source of noise is magnetostriction, which is known to be caused during the magnetisation process by the movement of 90° magnetic domain walls and the rotation of the magnetisation vector out of the easy directions. Therefore, it is obvious that when a ferromagnetic body is magnetised to saturation, the resulting magnetostriction will be dependent on the distribution of magnetic domains in the initial state.In the literature two coefficients are used to describe the resulting magnetostriction at a saturating magnetic field. The total strain, measured from an arbitrary demagnetised state to technical saturation, is called the saturation magnetostriction, and is denoted by (~)s. The tenn magnetostriction coefficient will be employed to designate the strain measured from the ideal demagnetised state (i.e. one in which no preferred orientation of domains exists) to technical saturation. Failure to recognise the difference between these two quantities has led to considerable confusion in the literature.The magnetostriction coefficient for a single crystal is a well-known property for ferromagnetic cubic materials. In 1939 Becker and Doring [1] have derived the following relation, describing the magnetostriction coefficient for single crystals with a cubic symmetry:"

Jan 1, 2001

By M. Patel, A. W. Piper, R. A. Forsdick

"Lead bullion is refined in kettles to remove the contained impurities such as copper, antimony, silver and bismuth. These impurities are removed sequentially by the addition of one, or more, reagents to the bath. Impurity elements react selectively with the reagents to form a dross layer on the surface that is subsequently removed for further processing. Mixing the lead bath using specific impeller designs, which form a vortex to ensure that good contact takes place, carries out the reactions. It is important that the mixing is controlled to achieve the optimum bath conditions and hence good dross grades.This paper presents initial modelling work on the kettle refining process using both physical modelling (Water Models) and mathematical modelling (Computational Fluid Dynamics - CFD). Comparisons between measured plant data and simulation results are given.IntroductionThe refining of lead bullion takes place in hemispherical kettles of varying sizes. It is normal practice to remove the impurity elements sequentially, by the addition of reagents, which selectively react with the impurity element. The impurity element forms a dross layer on the surface, which can be removed for subsequent processing at the end of the reaction period.All the refining stages in the production of lead and lead alloys involve mixing the lead bath, although different conditions are required for each stage. At Britannia Refined Metals mixing is achieved using either 3 or 4 angled or flat blade impellers, which generate a vortex in the central region of the kettle. However, throughout the industry a wide range of mixer types are used which are all trying to achieve the same result."

Jan 1, 1999

By Y. Cao, C. Li, A. Song

Bubble size plays a key role in determining flotation efficiency. Bubbles can be formed by passing air through pores. Bubble detachment from a pore is driven by lifting forces including buoyancy force and momentum force. Enhancing the momentum force is the key to accelerate bubble detachment and to minimise bubble size. Recent studies show that replacing steady air with oscillatory air supply could significantly decrease bubble size by enhancing the momentum force. To date, the underlying mechanism for the improved momentum force with oscillatory air supply is not yet clear. This work, for the first time, developed a qualitative model for the momentum force as a function of frequency, on/off time ratio and gas flow rate. It was found that the momentum force was positively correlated with the frequency as well as the air flow rate, and was negatively related to the on/off time ratio. The relation between bubble size and these operating parameters was also investigated. Keywords: Flotation, bubble size, momentum force, modelling, oscillatory air supply

Jan 1, 2020

By R. K. Brummer

A simplified two-dimensional plane strain numerical technique for the non- linear modelling of a tabular mining excavation in brittle host rock under high stress conditions is described. The model has been developed using relatively straightforward assumptions based on extensive underground investigation of the fracturing and deformation of rock at the stope face in deep gold mines on the Witwatersrand, together with observed rock behaviour under laboratory conditions. The reef rock is modelled as a linear string of yielding elements which interact horizontally, while the elastic rockmass remote from the reef is modelled using boundary element techniques. This elementary technique is shown to produce results which agree closely with the observed behaviour of the fractured rock during the mining of a remnant pillar

Jan 1, 1987

By A. F. Raschenko, Y. A. Golodcov, V. V. Elshin

"The autoclave noble metal desorption in apparatuses with a dense stationary sorbent layer is perspective for the technology of sorption gold and silver extraction from solutions and pulps by carbon sorbents. The study of the mechanism and optimization of the autoclave desorption process involve some difficulties. To overcome them it is necessary to use mathematical and quantum-chemical modelling. The given paper deals with a mathematical model for the dynamics of gold and silver desorption from active carbons by alkali solutions based on the physical-chemical notions of the process. The kinetics and the mechanism of the process were investigated using the mathematical model for gold and silver desorption.The calculation of the efficient mass-transfer coefficient (ße) for various initial and intermediate noble metal concentrations in carbon (Cy) made it possible to find that the dependence of ße on Cy in case of gold desoprtion shows a clear extreme character. The value of the initial metal concentration in carbon therewith does not influence the character of the curve the inflection point of which is changed essentially at elevated process temperatures only. The revealed characteristic property, in the author's opinion can only be indicative of the change in the mechanism of gold desorption under autoclave desorption conditions. The assumed [Au(CN)(OH)]- type intermediate that has a higher diffusion coefficient in comparison with the adsorbed complex [Au(CN)2]- can be considered as an intermediate compound in the process of gold desorption."

Jan 1, 2003

By C. O. Aksoy, P. M. Ozfirat, M. K. Ozfirat, C. Tatar, C. Pamukcu

"Temporal and other factors contributing to the occurrence of occupational accidents at an underground coal mine in Turkey between 1992 and 2000 have been analyzed. The time period of the study spanned the transition from hand-worked mining methods (1992 to 1996) to hand-worked mining methods and mechanized mining methods (1997 to 2000).The mine selected for the study employed an average of 5,835 men per year on three separate eight-hour production shifts for 302 d/y (Sundays are not worked, and there are 11 holiday days per year). Production on an annual basis averaged 190,000 tons per year between 1992 and 1996, and 300,000 tons per year between 1997 and 2000.The purpose of the study was to determine which, if any, of the factors studied contributed to the occurrence of accidents so that the mine could take the necessary steps to reduce the risk to the workforce. Regression analysis was used to determine the strength of the relationship between each factor studied and the number of occupational accidents.The strongest relationships observed were between the age of the employee, the day of the week, and the total yearly labour."

Jan 1, 2006

By Mark Jolly

"Over the last ten years or so many man hours have been dedicated to developing methods for simulating shape casting processes. The majority of the methods developed have been devised by a combination of computer, mathematics and materials specialists often with limited knowledge of foundries and foundrymen. As a result of this it would appear that although there are many software packages for the foundryman to use there still exist fundamental misunderstandings as to their usefulness and suitability within the foundry.The paper aims to identify how well current software packages perform. The relationships between physical phenomena, practical defects and software capabilities are presented. Some discussion of the use of criteria functions is also presented. Finally issues arising from the postprocessing of results are discussed as the presentation of the results to the lay-person in simulation techniques is possibly one of the most important aspects in influencing the adoption of this type of software in the foundry community."

Jan 1, 1999

By J. Petersen, S. W. Robertson

Heap leaching is used extensively for the processing of low-grade copper, gold, and uranium ores. The hydrodynamics is often the most important factor that limits the rate of metal recovery. Over the past decades the author has performed physical and hydraulic testing of a large number of ore samples considered for heap leaching. The data were reviewed to investigate the effect of physical properties such as particle size distribution (PSD) and bulk density on the hydrology. The unsaturated flow through the bed is traditionally modelled with soil science models such as Brooks-Corey (BC) and Van Genuchten-Mualem (VGM), which relate the hydraulic conductivity to the degree of saturation, also known as the hydraulic conductivity function (HCF). However, most authors found a disconnect between the saturated hydraulic conductivity (KS) and the rest of the HCF, which is explained as resulting from a transition from capillary-to gravity-controlled flow. However, the author found no such transition, and found the flow to be governed by capillary flow throughout. Traditional BC or VGM fits were found to be appropriate, however it was necessary to modify the model to include a point of inflection, above which the hydraulic conductivity remains constant while the moisture content keeps on increasing. It was proposed that this is the result of unconnected or poorly connected pores, which fill with solution but do not create additional flow channels. The point of discontinuity was also found to correspond with the air entry point, measured independently from the air conductivity curves. On average, the air entry point corresponds to the experiential rule that the void saturation should remain below 65% for aerated heaps.

Dec 11, 2024

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

Experimental studies on viscosity of titanium bearing slags have been reviewed and the assessed data were used in extending the capability of the structurally related viscosity model for silicate melts, developed by the present authors. Experimental data on ternary and higher order systems containing CaO, MgO, MnO, FeO, Al2O3, TiO2, Ti2O3, SiO2were assessed and used. The newly extended viscosity model represented the viscosity behaviour over broad temperature and composition ranges. The structural role played by TiO2and Ti2O3in silicate network structure is discussed with respect to their effect on viscosity of melts. Keywords: viscosity model, TiOx containing silicate melts, silicate network structure.

Jan 1, 2004

By W. Wittke, J. Gattermann

ABSTRACT: One alternative to close a radioactive waste repository in rock salt is the construction of a “Cross Section Closure (CSC)“ sealing barrier. The proposed material for the CSC is highly compacted bentonite. To investigate the swelling behaviour of a highly compacted bentonite a large number of laboratory tests were performed. In addition, large scale model test were carried out to demonstrate the development of a nearly homogeneous and isotropic swelling pressure. The results of the large scale model tests were interpreted numerically based on the models for stress-strain behaviour and water uptake which are implemented in the finite element codes FEST03 and HYDOPO. The investigations show a good agreement of the results of the model tests and the corresponding analyses and show, that the theoretical models are capable of describing the behaviour of sealing structures based on highly compacted bentonite. ZUSAMMENFASSUNG: Eine Möglichkeit für den Verschluß eines Endlagers für radioaktive Abfälle im Steinsalz ist die Konstruktion einer Querschnittsabdichtung aus hochverdichtetem Bentonit. Zur Untersuchung des Quellverhaltens von hochverdichtetem Bentonit wurden eine große Anzahl von Laborversuchen durchgeführt. Um die Entwicklung eines nahezu homogenen und isotropen Quelldrucks nachzuweisen, wurden darüber hinaus auch Modellversuche in größerem Maßstab durchgeführt. Die Ergebnisse dieser Modellversuche wurden im Hinblick auf die Wasseraufnahme und das Spannungs-Verformungsverhalten mit Hilfe numerischer Berechnungen mit den finite Element Programmen FEST 03 und HYDOPO interpretiert. Die Untersuchungen zeigen gute Übereinstimmungen zwischen Versuchs- und Berechnungsergebnissen. Die numerischen Modelle sind somit geeignet, das Verhalten von Abdichtungsbauwerken aus hochverdichtetem Bentonit zutreffend zu beschreiben. RÉSUMÉ: Une possibilité pour fermer un dépôt final de déchets radioactifs dans sel gemme est la construction d'un étanchement de coupe transversale, fabriqué de bentonite fortement compacté. Pour étudier le comportement de gonflement de bentonite fortement compacté, de nombreux essais en laboratoire ont été exécutés. En outre, des essais sur modèle sur une plus grande échelle ont aussi été exécutes afin de prouver le développement d'une pression de gonflement à peu près homogène et isotrope. Les résultats de ces essais sur modèle étaient interprétés en vue de l'absorption d'eau et du comportement contraintes-déformations par le moyen de calculs numériques avec les programmes d'éléments finis FEST03 et HYDOPO. Les études montrent de bonnes conformités entre des résultats d'essai et de calcul. Les modèles numériques sont ainsi capables de décrire justement le comportement d'ouvrages d'étanchement fabriqués de bentonite fortement compacté.

Jan 1, 2003

By Pengfu Tan

The control of converter slag chemistry and slag viscosity has played an important role in optimizing converter operations. A thermodynamic model of copper P-S converter and a viscosity model of converter slag have been used to predict the chemistry and viscosity of converter slag. The effect of oxygen potential, SiO2/Fe in slag, detection of slag blow endpoint, and temperature on slag viscosity, liquidus temperature of converter slag and skimming operation has been discussed. Some applications of the industrial operations have been presented.

Jan 1, 2006

By Andrew L. Mular

Methodology to model, simulate and optimize full scale mineral pro- cessing circuits in production environments is reviewed and evaluated. Strategies for data acquisition and adjustment are considered along with techniques to fit models to adjusted data. Examples are considered for clarity.

Jan 1, 1989