Search Documents

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 Roland B. Drugge

The adaptation of pelletizing processes involves the coordination of various factors such as temperatures, gas flows, productions, etc., in order to obtain optimum performance. The consequent problems are considerable due to the large number involved as well as the lack of knowledge of their correlation. It is therefore essential to establish a fundamental theory which takes into account the heat transfer phenomena and the secondary factors which are relevant to the properties of the product. Straight grate and grate kiln pelletizing processes involve heating of a stationary bed on a travelling grate. The heat transfer conditions are therefore of primary importance. This is why this article pays particular attention to conditions which obtain when a stationary bed is heated, in which certain special reactions occur. The conclusions drawn from this penetration of the heat transfer phenomena can be expressed in postulates acting as guide-lines for the establishment of proto-typal processes. Examples of elaborating such processes are also given in the article. Some operational experience at LKAB are also given.

Jan 1, 1975

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 H. Lind

The consumption rate of Estonian oil shale increases due to the growing rate of exported oil-shale oil. Estonian oil shale active resources in working mines are predicted until the year of 2025 when the consumption rate stays at the level of 12 tons per year. If the consumption rate stays at that level, new mines should be opened after 20 years, if the rate increases new mines should be opened even earlier. The predicted active resources in Estonian oil shale deposits are calculated for 60 years, based on the parameters of the power station technicoeconomics. As a result, there is a need to open new mines. One of the potential areas for surface mine is the Tammiku-Kose prospect area (Rheinbaun Engineering, 1996).

Jan 1, 2006

By M. Morales, J. Yianatos

"In this work, an alternative approach to evaluate the flotation performance along a flotation rougher row is presented. The collection zone as well as the froth zone performances were separately evaluated. For this purpose, sampling for metallurgical balances along with sampling of top-of-froth (TOF) grades in each cell of a flotation row were employed. The modelling of the rougher circuit consisted of using TOF grades for evaluating the collection zone recovery, because they are a good representation of the bubble loading, which identifies the true flotation process, as shown in previous studies for shallow froths. In addition, the feed characterization, size-by-liberation, allowed the identification of the mineral distribution and its evolution along the circuit. The latter result was obtained by fitting the recovery profiles as well as the TOF grades down the flotation bank. The flotation rates and the maximum recoveries in terms of size-by-liberation (size classes: -45 µm, +45-150 µm, +150 µm, and liberation classes: -20%, +20%-80%, +80%) were then estimated. The results were in good agreement with those reported in literature by other authors and showed the new approach allows flotation kinetics to be characterized in industrial circuits.INTRODUCTIONMineral recovery in industrial flotation circuits is a function of the mineral characteristics, such as particle size distribution and liberation (Trahar, 1981, Jameson, 2012), and operating conditions such as air flowrate, pulp flowrate, solids percent and reagent conditioning (Jameson 1977; Klimpel, 1980; Mehrotra and Kapur, 1974). The cell type and size (self-aerated or forced air) and the circuit arrangement (single cells or banks of cells in series) are also relevant factors. Sampling of circuits (e.g. rougher, cleaner or scavenger) are commonly carried out using standard automatic samplers for metallurgical balance assessment. However, the flotation plants are generally not provided with proper accesses for local streams, which hinders the recovery estimations and the flotation rate characterization."

Jan 1, 2018

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 Y Chung, B Choi, H Lee

Numerous studies have explored refractory corrosion in the presence of slag and/or metal, utilising the Finger Rotating Test (FRT). Recent experimental results revealed that corrosion initially increased with rising rotation velocity but was then suppressed above a critical rotating speed, under conditions involving a two-liquid solution of slag and molten steel. This phenomenon indicated that the liquid/liquid interface movement occurred during spindle rotor. However, understanding of the phenomenon, known as the rod climbing effect, remains inadequate in the case involving slag, molten steel and refractories. To gain deeper insights, a simulation experiment was conducted on the phenomenon using silicon oil/water and the phenomenon was analysed using computerised simulation software (Ansys Fluent). Experimental modelling demonstrated that interfacial movement occurs even at low rotation velocity, as low as 200 rev/min, when the viscosity of silicone oil is at 100 cP. Computational analysis utilising the Volume of Fluid (VOF) method and a realisable k-ε turbulence model elucidated that interface elevation during rotor rotation arises due to viscosity disparities between the two fluids. Furthermore, this movement is more pronounced when the viscosity difference is lower. Additionally, the study findings suggest that higher density differences between two liquids suppress the phenomenon of interface movement.

Jun 19, 2024

By A. A. Hejja

"Liquidus temperature and electrical conductivity data measured on synthetic slags were modelled by the use of neural nets (NN). In this work the applied multilayer feedforward NNs were trained by a conjugate-gradient optimisation for which a three layer formulation was used. Very good fits were obtained for both the liquidus temperature and the conductivity data. The synthetic slags were prepared from pure oxides to represent a wide range of compositions likely to be encountered in ferromanganese and silicomanganese smelting. The slag constituents were in the following range: Mn0;5-30%, Ca0;20-35%, Mg0;5-15%, Si02;27-58%, Al20 3;5%. Liquidus temperatures varied from 1300° to 1380° and increased with increasing basicity ratio. The electrical resistivity of slags decreased with the increase of basicity ratio from 0.55 to 1.1 but above 1.1 basicity ratio the resistivity tended to increase depending upon the MnO content.Neural NetsVarious neural net (NN) topologies such as Hopfied nets, Hamming nets, Carpenter/Grossberg classifiers have been discussed in some detail by Lipmann (1). In this work, the applied multilayer feedforward NNs are not trained by backpropagation (BP) using the generalized delta rule (2) (GDR), but by a conjugate-gradient optimisation [3,4] - a procedure similar to that proposed by Barnard and Cole [5]. A neural net is capable ·of learning the non-linear functional relationship existing between a set of inputs and outputs [ 6]. In this capacity a neural net is used to determine the functional relationship that exists . between N dependent and independent variables, which characterize the system under consideration For a three layer neural net [1] the following function may be formulated, which relates X; inputs to the output Y 0 from node o:"

Jan 1, 1997

By Bo Zhou

During secondary aluminum recovery, complex aluminum scraps are melted and refined often in a rotary melting furnace and the aluminum metal has to pass through a molten salt layer and melt down in a bottom aluminum bath. Previous experimental study and industrial observations have indicated that salt shell formation and re-melting on the metal solids in the early stage plays a critical role in the overall melting process of the scrap. In the present paper, mathematical models are developed to simulate the process of salt shell formation and its re-melting, as well as aluminum melting in the molten melt. Influence of the salt layer properties, the residence time of the metal particles in the salt layer, the particle size, the temperature of bulk melt were investigated. Reasonable agreement with the measurement data was reached. The developed melting model will be used to construct a population balance model of the scrap melting behavior for a rotary melting furnace, taking into account the distributed characteristics of the scrap. Finally, an overall process model based on Computational Fluid Dynamics (CFD) for the rotary furnace will be constructed, coupling with the population balance model for the scrap melting.

Jan 1, 2003

By Bo Zhou

During secondary aluminum recovery, complex aluminum scraps are melted and refined often in a rotary melting furnace and the aluminum metal has to pass through a molten salt layer and melt down in a bottom aluminum bath. Previous experimental study and industrial observations have indicated that salt shell formation and re-melting on the metal solids in the early stage plays a critical role in the overall melting process of the scrap. In the present paper, mathematical models are developed to simulate the process of salt shell formation and its re melting, as well as aluminum melting in the molten melt. Influence of the salt layer properties, the residence time of the metal particles in the salt layer, the particle size, the temperature of bulk melt were investigated. Reasonable agreement with the measurement data was reached. The developed melting model will be used to construct a population balance model of the scrap melting behavior for a rotary melting furnace, taking into account the distributed characteristics of the scrap. Finally, an overall process model based on Computational Fluid Dynamics (CFD) for the rotary furnace will be constructed, coupling with the population balance model for the scrap melting.

Jan 1, 2003

By E Werner, J J. Donovan, A J. Morris

Coal mining in the Pittsburgh coal seam, in the eastern United States, began in the late-1700s. Since that time numerous mining technologies have been used throughout the basin. Environmental concerns, particularly mine discharge quality, have only been considered since the early-1970s. During the last 20 years, many of the larger underground mines have closed and are now either flooding or already flooded. In order to understand the consequences and potential future effects of this flooding, mine mapping of the basin has been compiled utilising a geographic information system (GIS) with an emphasis on locating barrier pillars between mines. Mine flooding has been monitored at numerous sites utilising pressure transducers and data loggers. These data have been combined with historical data from mining companies and structure contour maps to yield chronological flooding maps of the basin from 1980 to present. Utilising these maps and the indicated flooding rates, future discharge locations have been identified along with estimates of when these discharges are likely to occur. Related work has identified barrier pillar hydraulic conductivities and vertical infiltration rates as a function of overburden thickness. All of this information has been integrated into a groundwater flow model that simulates mine flooding. This model is useful in defining mine to mine interactions as the various mines fill with water, and the testing of various flooding scenarios. The model will be used in evaluating mine discharge rates and evaluating the options for regional mine drainage treatment.

Jan 1, 2003

By P B. Hills

This paper discusses the computerisation of the drill hole database, some historical records and the creation of a three-dimensional model of the gold bearing structures for resource estimation and mine design at the Beaconsfield Gold Mine in northern Tasmania. The phase of the work reported here commenced in November 1997 and was completed in March 1998. The historical mining on the Tasmania Reef had reached a depth of 450 m below surface when work ceased in 1914. Production up to this time was 854 000 ounces gold. Recent activity has concentrated on testing depth extensions of the previously mined ore structure, to a depth of 1000 m below surface, and the establishment of an underground mining operation to exploit this resource. The lode structure at Beaconsfield is a typical mesothermal vein style deposit. Gold is contained within narrow quartz, carbonate and minor sulphide reefs or veins. Reef width is variable from a few centimetres to several metres with æsplitsÆ or bifurcations being common. The reef averages 2.7 m true width in the economic zones. The key feature of this style of mineralisation is that reef continuity can be confirmed by drilling to a high level of confidence, but the gold grades are somewhat more erratic in nature. The resource estimates were based entirely on drill data. The drill intercepts on the reef were coded for each split directly into a database and these allocations validated graphically. The distribution of grades and widths show an increasing gold grade with reef width for both the drill hole and historical data. The same histogram patterns, mixed populations and overall arithmetic average grade for the historical mine (at 21.52 g/t > 3 g/t Au) and new drill data (at 21.55 g/t > 3 g/t Au) would indicate that there has been no major change in the overall style or grade of the ore from the historically mined areas to the new resource. The reef, host geology, area of historical workings and resource zones at the Beaconsfield Mine were defined as a series of closed wireframes. The Tasmania Reef was modelled as five mutually exclusive wireframes, one each for the Main Reef, Footwall split, Second footwall split, Third footwall split and Hangingwall split. The data for this exercise was extracted directly from the drill hole database, stored in Microsoft Access format. The formation of the solids was largely automated with manual validation and refinements allowing for multiple runs using different interpolation methods and combinations of drill holes. A rotated 3D-block model in the plane of the Tasmania Reef was created to contain the grade model and allow for tonnage, grade and reef width estimates. An inverse distance squared (ID2) interpolation method was used with ranges of 25 m along strike and 125 m down plunge with blank areas infilled with 50 m by 250 m. Each reef split was filled separately using only drill intercepts from that split. The method was independently validated used MIK methods. Longitudinal sections have been prepared showing gold grade, reef width and width +grade. Estimates were also made for arsenic, sulphur and other minor elements. The long sections show three distinct mining zones. A method is presented for the conversion of the resources to reserves by the addition of mining parameters for minimum mining width, minimum pillar, cut-off grade and internal and external dilution. The selection is done on columns within the block model, and allows for all splits to be considered at once. Multiple runs can be made for varying combinations of parameters to test different mining methods and their impact on reserves, mine life and NPV.

Jan 1, 1999

By M. Kozlowski, I. Candy, S. Kashani-Nejad

The production of ferronickel and matte from nickel laterite ores is dominated by the Rotary Kiln, Electric Furnace (RKEF) process. Over time, major improvements have been made to kiln operation. In particular, kilns are now required to perform substantial amounts of reduction of the metallic oxides, thereby reducing the load on the electric furnace and substituting less expensive fossil fuel for electrical energy. To the basic refractory lined cylinder have been added: internal features such as dams and lifters, aimed at changing residence times in different sections of the kiln and enhancing gas/solids interaction, on-board fans to introduce tertiary air and on-board scoop mechanisms to allow reductant coal additions to be made at advantageous locations along the kiln. Many of these advances in design have been made on a trial and error basis and by 'borrowing' technology from related kiln processes, supplemented by proprietary modeling to various levels of sophistication. This paper presents a conceptual review of modeling components for the development of a computational model for granular material flow, combustion and heat and mass transfer in a nickel laterite kiln. These components include: granular material flow and segregation models for the simulation of bed motion, heat transfer models that simulate heat transfer by turbulent diffusion and radiation and the distribution of this thermal energy within the bed, and fluid flow models for the prediction of combustion aerodynamics.

Jan 1, 2005

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 G. P. Demopoulos

The solubilities of calcium sulfate in HCl-CaC1z-H20 solutions up to 100°C were successfully modelled with the aid of OLI software after bringing proper modification to it. Modelling, more specifically, was carried out with the help of the Environmental Simulation Program, V-6.6.0.4 (www.olisystems.com). New model parameters were obtained by regressing the solubility data. A private databank was established by adding the newly defined model parameters and the calcium sulfate hemihydrate phase that is not included in the original OLI database. The McGill-modified OLI software was successfully applied not only to modelling solubilities but also the transformation of the various CaS04 phases, namely dihydrate, hemihydrate and anhydrite. In this paper, the modelling methodology will be discussed along with a comparison of the capabilities of the newly developed model vs the original OLI model.

Jan 1, 2005

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 G. Brooks, J. Naser, M. A. Sattar

CFD models have been developed and numerical simulations have been carried out to predict the formation of foam in oxygen steelmaking. Foam was considered as a separate phase comprising a mixture of gas and liquid. Bubble break up and coalescence models have also been incorporated in a CFD model to predict the number density of individual bubble classes. A population balance equation was used to track the number density of each bubble class. Decarburization with heat generation from chemical reactions was integrated in the process. User subroutines were written in FORTRAN to incorporate the foam formation, the bubble break up and coalescence rate and decarburisation in the main program. The model predicted the foam height, bubble number density, velocity of phases, decarburization, and turbulence. The result from the model has been compared with available data from literature and found to be in reasonable agreement with the experimental and plant data.

Jan 1, 2014

By J. L. Liow

Studies have been carried out on the causes and effects of slopping in an annular cylindrical geometry such as the top submerged combustion reactor. Properties of standing waves formed such as the wave profile, period, pressure distribution and velocity have been derived analytically. The wear profile due to wave erosion was calculated based on the assumption of linear motion and that the wear is a dissolution process. The use of the j, factor showed that the wear was proportional to the square root of the wave velocity at any point. Using the dimensions of the lsasmelt reactor given by Errington et al. (1987)[11, the wear profile due to wave erosion was calculated using the analytical solution. The wear profile calculated was shown to have a maximum just below the initial bath surface. The wear was found to be concentrated in the region where the first asymmetric standing wave was most active. The results showed that the refractory near the bath surface will be affected by wave erosion, more so than any other parts of the reactor. The contribution of wave erosion on wear must be accounted for in assessing the total refractory wear of a vessel in which slopping exists.

Jan 1, 1991

By A. Chagnes

This paper presents a modelling tool that can be used to predict the behaviour of solvent extraction processes. In particular, this model allows calculating metal extraction efficiency as a function of pH, solvent flowrate, stripping flowrate and aqueous feed flowrate in the mixers-settlers. This model is applied to the recovery of cobalt from acidic chloride media by using CYANEX® 272 as extractant diluted in kerosene. Keywords: Cobalt, mixer-settler, modelling, solvent extraction

Jan 1, 2020

By P. Myllykoski, J. Nylander

"A commercial program package called Crownoff for the simulation of temper rolling in four high mills has been verified against measured data. The agreement between measured and calculated transversal strain distributions of the strip was found to be good. Effects of the strip crown and work roll bending force on the strain distribution were simulated. The results showed that when a single ground camber of work rolls is used for all strip widths; the strip crown should be parabolic, excluding the edge drop. An increase or decrease of the strip crown by 30 % during a coil can be compensated by a 2 to 5 t change in the work roll bending force. The ability to control the strip crown by work roll bending was insensitive to strip width. IntroductionCold rolled steel strips are usually temper rolled after recrystallization annealing, which follows cold rolling. The purpose of temper rolling is to suppress the yield point, finish the strip surface and correct possible flatness errors. The presence of a yield point might otherwise cause problems during subsequent forming processes. The in homogenous yielding in non-temper rolled steels leads to the formation of Lüders bands in the strip surface.The rolling force is influenced by the rolling speed and changes in frictional conditions in the roll bite. Therefore it is very important to plan carefully the mill settings and control algorithms during rolling. The purpose of mill control algorithms is to ensure constant reduction and the best possible strip flatness during the entire coil.Correctly predicting the rolling force is important for determining the mill presettings. Torque, forward slip and power requirements are not usually critical, if the reductions are small. However, the large elastic deformations occurring in the roll bite make the physical description of the process difficult. Physical rolling models developed for cold rolling contain assumptions that lose their validity as the reductions decrease to the level of one per cent."

Jan 1, 1994