Search Documents

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

Jan 1, 1993

By Russell Barr, Robert E. Cook

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

Jan 1, 2007

By L. G. Austin

The procedure for scaling breakage parameters determined in a laboratory mill to values for a full-scale mill is briefly presented. A simulation model of a closed circuit mill also requires a model of residence time distribution, a mass transfer law to predict over-filling as a function of flow rate through the mill, and a model to describe classification action. For fine dry grinding, or viscous slurry grinding, it is necessary to allow for slowing-down of breakage rates for a fine mill charge. Four examples of comparing simulation models to actual plant data are presented. The first is wet grinding of copper ore in normal closed circuit at 50 tph (metric) production, under conditions where the mill does not overfill. The second is the wet grinding of copper ore in reversed closed circuit at production of 87, 102 and 108 tph, with flow through the mill of approximately 150, 200 and 300 tph. The third is dry grinding of coal with a hot air sweeping stream to remove all of the coal as air-borne material, in open and normal closed circuit. The fourth is grinding of cement clinker in a two-compartment, air-swept tube mill in normal closed circuit with two different classifiers, producing 120 tph of 3250 Blaine cement. In all cases the simulated capacities and product size distributions were in reasonable agreement with actual plant performance.

Jan 1, 1984

By Paul K. Krueger

SUMMARY Since 1970 the Federal Preparedness Agency has been using computer models to help establish the amounts of strategic and critical materials to be held in the National Stockpile. Fifty-three of these materials are metals or minerals. The following diagram shows the linkage of models used to estimate material requirements. [ ] The starting point for the models are estimates of the major components of GNP for some future period. These may be estimates purchased from a forecasting service or may be produced in-house. The projections are translated into industrial final demands using DITT (Demand Impact Transformation Tables), a system developed at FPA. Then using a version of the current Bureau of Economic Analysis input-output model, estimates are made of the gross output of each industry necessary to support the level and composition of GNP as specified by the macro-economic model.

Jan 1, 1977

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

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

Jan 1, 2000

By Donald Thompson, Jessie Bonner, Ileana Tibuleac

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

Jan 1, 2005

By N. El-Kaddah, V. L. Acoff, M. Arenas

"Control of the weld penetration, weld pool shape and temperature field in arc welding is critical to achieve the desired microstructure and mechanical properties in the resultant welds. This is particularly important in welding gamma titanium aluminides (y-TiAl) where the weld pool is relatively shallow and the temperature gradient in the heat-affected zone is quite steep. This paper investigates theoretically and experimentally the evolution of the shape of the weld pool for stationary gas tungsten arc welds in y-TiAl alloys. A mathematical model has been developed for predicting the electric current flow, the resultant electromagnetic forces in the metal, and the transient development of the temperature and velocity fields in the molten pool generated by a spatially variable heat flux at the arc-striking zone. The governing equations were solved numerically using the control volume technique. Welding experiments were carried out to determine the shape of the weld pool and the extent of the heat affected zone at different welding times for model validation. The numerical predictions were found to be in reasonable agreement with measurements. It is suggested that the model presented in this paper is likely to provide the information needed for predicting the microstructure and mechanical properties of y-TiAl welds. IntroductionIn recent years, there has been considerable interest in the welding of gamma titanium aluminides as a means of manufacturing high-temperature, high-performance components for turbine and automobile engines. These alloys are highly reactive, and must be welded in vacuum or inert gas atmosphere, rendering electron beam, and inert gas tungsten arc welding (GTA) to be the only suitable welding techniques for joining these materials [1-3]. The GTA welding is especially attractive because of its flexibility and relative cost-effectiveness. However, previous studies on the weldability and weld properties of these alloys using this process revealed a decrease in the mechanical properties in the fusion zone [4]. This is generally attributed to the change of the original duplex structure to a dendritic structure amplified by the shallowness of weld spot. Quantifying the effect of arc operating parameters on the geometry and temperature profile of the welding pool prior cooling is the first step for controlling the fusion zone microstructure and mechanical properties."

Jan 1, 1999

By Peter Klaver

Deep tunnels have been used for decades for the temporary storage of large volumes of urban runoff, both storm water and combined sewer overflow (CSO). Tunnels are attractive because they reduce surface disruption and do double duty as storage and conveyance systems, collecting overflows from multiple locations and transporting them for treatment to various facilities. Like any hydraulic system, however, there are dynamic effects that must be fully understood during the design phase so that the final constructed product does not have unintended negative consequences. Preliminary sizing of storage tunnels is typically based on the total volume required. For example, if studies have determined that a volume of 56.8 ML (15 million gal) will provide the intended level of control, then that figure can be used, along with a preliminary idea of the desired tunnel length, to compute the required diameter. Of course, filling a long tube with water from various points along the way is not at all the same as pouring water into a bucket, and these dynamic effects require additional attention from design engineers. The rapid filling of storage tunnels can lead to the formation of bores, which are essentially moving walls of water. An open channel bore can build in height until it reaches the crown of the tunnel, at which point it becomes a pipe-filling bore. Figure 1 depicts both kinds of bores moving upwards in a rapidly filling tunnel. Bores can strike the ends of tun-nels with tremendous force, driving water and entrained air upwards into dropshafts and creating strong pressure surges in portions of the tunnel that are already filled. If the hydraulic surges are high enough to reach grade, they can cause overflows and/or damage structures at the surface. Trapped air pockets can occur through several phenomena, including the reflection of bores off dropshaft structures and the occurrence of ?premature pressurization.? The latter refers to the situation where the inflow to a tunnel is increasing more rapidly than the tunnel can convey it, such that the pipe becomes surcharged at an upstream point while points downstream are still flowing as an open channel. These air pockets are a great nuisance as they migrate toward points in the tunnel where they can rise to the surface, which are often shafts that are already filled with water (Fig. 2). The ?geyersing? that results can be hazardous and costly.

Jan 3, 2010

By Angelus S. C. Owili-Eger

Although the foundation of a unified treatment of fluid flow in variably saturated porous media was established more than seventy years ago, no serious attempt was made to develop it until the early 1970's. Ever since, a number of unsaturated-saturated numerical models have been reported in the literature. However, none of these models have treated multidimensional, multiphase flow in deforming media, parts of which may undergo extensive secondary fracturing. A situation such as this is not uncommon in coal mining operations in which the longwall method of extraction is employed. The present model has been structured to handle such physical systems. This paper presents the mathematical formulation, the numerical algorithm, and a further validation and evaluation of the modified model. The model treats a multiphase flow in a three- dimensional regime and handles mining excavations of complex and changing configurations.

Jan 1, 1980

By A. S. C. Owili-eger

The computer-aided simulation model re¬ported here is an extension of the one developed by the author while at The Pennsylvania State University. The present model, however, includes consideration of the deformability and secondary fracturing of the porous system. The modifications enable simulation of multiphase shallow unsaturated-saturated and/or artesian ground water systems in dynamic mining environments. The model has been validated, verified, and evaluated on a number of field case studies and the results have been excellent.

Jan 1, 1984

By Tsvetan Damyanov, Nikolay Hristov, Evgeni Kostadinov, Lubomir Kuzev

"The process of grinding in drum mills is directly dependent on the dynamic forces occuring in the free movement of the grinding bodies in the volume. The subject of discussion are the potential impacts between grinding bodies having different shape and equal weight in a drum mill that has a diameter of 1200 mm. The study utilised ball grinding bodies having a diameter of 120 mm and tetraballpebs having a conditional diameter of 120 mm. The study was carried out using the Solid Works Simulation software. The following cases were examined: • Impact between two identical grinding bodies; • Impact between a grinding body and the shell surface. The study provided numerical values for impact force, kinetic energy of impact, linear speed of impact, linear displacment of the center of gravity and duration of impact deformation. The impact force of the balls is approximately 179000 N. The impact force of the tetraballpebs with typical vertices, edges and spherical surfaces depending on the area of impact is 177000 to 185000 N. Upon impact between two spherical surfaces of tetraballpebs, the contact area is 263.55 mm2, which is more than twice the contact area upon impact between two balls, which is 114.39 mm2. This determines better conditions for the transmission of kinetic energy by reducing opportunities for overgrinding."

Jan 1, 2014

By C. Pomeroy, J. Weyer, J. Brune

"Mining began in the Erzgebirge (“ore mountains”) region of Saxony, Germany more than 800 years ago, and extensive underground workings cut through the entire area. Some of these old workings are still used today: not for production, but for offering an underground venue for research, tourism and university teaching. One example is the Reiche Zeche mine at the Technische Universität Bergakadamie Freiberg which hosts one of the oldest mining engineering programs in the world. It is imperative to keep all active mine workings well ventilated to make the mine safe for students, visitors and researchers. Since the mine workings at Reiche Zeche are so vast, extensive and, to a large extent unexplored and unmapped, it has been difficult to provide a complete, functioning and adequate mine ventilation model for this mine that accounts for all known and unknown leakage air flows. Researchers at the Colorado School of Mines have cooperation with their colleagues at Freiberg to establish a mine ventilation network model based on a limited number of air readings in known and accessible drifts. In order to balance the measured air quantities, leakage flows were quantified and mapped by comparing input and output flows of mapped drifts, and modelled using Ventsim Visual, a powerful visualization software package. The resulting mine ventilation model is used to locate and quantify the leakage air flows and to examine how much air is lost in each area. With this model, ventilation changes and improvements underground can be documented reflected in the model. This enables the Bergakadamie Freiberg to project air changes required for scientific experiments and to set up mine ventilation in specific areas and for specific needs. The model also enables better ventilation control and monitoring than before. Finally, the model can be used to predict ventilation changes and escape way contamination in case of a fire, ground fall or inundation in the mine. The latter becomes important as the lower mine workings are often flooded during spring runoff, which has created uncertainty about mine ventilation in the past. HISTORICAL BACKGROUND The discovery of silver, lead, zinc, tin and other valuable metal ores in the Erzgebirge region of Saxony, Germany began in the12th century and started a mining boom than lasted for centuries. Surface mines quickly developed underground to access the higher grade veins. Explored by the early geologists and driven by the ingenuity of the miners, the underground works expanded vastly throughout Freiberg and the surrounding region. Development followed the highest grade veins, and engineers developed solutions to drain water and provide limited ventilation to the underground workings. Still, engineered and planned ventilation was usually not a priority in the development of a new drift, which today makes ventilation modeling a significant challenge."

Jan 1, 2015

By V. Bojarevics, K. Pericleous

"A number of different methods have been developed for noncontact electromagnetic treatment of liquid metal and to investigate the melting/solidification processes. Applying AC magnetic field in terrestrial conditions, along with the buoyancy and thermo-capillary forces, results in turbulent flow. The use of a homogenous DC magnetic field allows damping both of the turbulence and the large scale flow at different rates leading to surprising results at some stages. The dynamic interaction of the turbulent flow with the oscillating interface needs to be accounted if it is required to levitate liquid metal of up to few kilograms without the contact to container walls. At the high values of magnetic field some oscillation modes are damped quickly, while others are modified with a considerable shift of the oscillating droplet frequencies and the damping constants from the non-magnetic case. Numerical models are used to investigate behavior of liquid metal suspended in high DC magnetic fields providing microgravity-like conditions with laminar viscosity and heat transfer to the surroundings.IntroductionSince the early Electro Magnetic (EM) levitation experiments [1] in AC magnetic field it became apparent that the levitated liquid metal is prone to oscillation and instability. The EM and electrostatic levitation experiments with liquid metal droplets show difficulties related to confinement stability and a need for complex correction functions to establish a correlation between the measurements and the droplet material properties [2,3]. Intense internal fluid flow is visually observed, apparently being in the turbulent regime for earthbound conditions."

Jan 1, 2012

"In the conceptual stages of engineering, benchmarking and empirical ratios are commonly the only references for estimating the production rate of greenfield projects. For underground mining projects, the lack of similar cases often makes the comparison difficult due to each project’s unique geology, rock quality or other singular characteristic. This paper discusses an analytical alternative for determining the maximum achievable production rate for underground stoping mining by modeling production as a continuous process. The results show that long-term maximum capacity depends exclusively on the achievable development rate and the development intensity of the design. In operations constrained by backfill, the development rate could be limited by the backfilling rate. Extraction rate and stope size only have impacts on capacity estimates during the ramp-up period.These findings allow shaping of the maximum capacity profile of a project. They can be used as a tool for quick assessment of conceptual capacities, and sensitivities of mining rates. Benchmarking tests compared favorably with projects, mines in operation and scheduling exercises. IntroductionA key consideration in the value proposition of a mining project is its capacity or production rate (Tatman, 2001). At a conceptual level, it is common practice to refer to empirical references (Taylor, 1986; Long, 2009), rules of thumb (de la Vergne, 2014), benchmark ing studies or consultant experience to estimate or validate a target value. There is, however, a lack of analytical methodology to guide an early estimate of the maximum capacity.Araneda and Gaete (2004) presented an interesting mathematical approach to estimate the capacity of panel caving mines. Using a continuous model, they confirmed that the mine capacity is essentially determined by the rate of area incorporation, in terms of area per period, and ore column height. This simple relationship has been confirmed by multiple detailed scheduling works"

May 1, 2016

By F. M. Doyle

Metal removal kinetics and metal ion selectivity in ion flotation were studied from a thermodynamic perspective. Surface tension data could be used with the Gibbs adsorption equation to quantitatively predict the metal-ion removal kinetics in alkylsulfate-Cu systems. However, surface tensions predicted cation selectivity coefficients in the sodium dodecylsulfate (SDS)-Cu-Ca and SDS-Cu-Pb systems rather poorly, although the predicted order of the selectivity (Cu2+ < Ca2+ < Pb2+) was correct. Another thermodynamic selectivity model, which used the Grahame adsorption equation and a geometrical analysis, predicted cation selectivity coefficients that agreed very well with experimentally measured values.

Jan 1, 2002

By F. M. Doyle

Metal removal kinetics and metal ion selectivity in ion flotation were studied from a thermodynamic perspective. Surface tension data could be used with the Gibbs adsorption equation to quantitatively predict the metal ion removal kinetics in alkylsulfate-Cu systems. However, surface tensions predicted cation selectivity coefficients in the sodium dodecyl-sulfate (SDS)-Cu-Ca and SDS-Cu-Pb systems rather poorly, although the predicted order of the selectivity (Cu2+<Ca2+<Pb2+) was correct. Another thermodynamic selectivity model, which used the Grahame adsorption equation and a geometrical analysis, predicted cation selectivity coefficients that agreed very well with experimentally measured values.

Jan 1, 2000

By K. D. Luxbacher

Understanding methane emissions in underground coal mines is critical for a safe and productive mine. In addition to reasonable estimation of initial coalbed reservoir parameters, it is also crucial that changes in effective stress due to mining and pore pressure reduction are taken into account due to their effects on porosity and permeability. Primary parameters for estimation of emissions or modeling of the mining environment for this purpose are porosity and permeability which can change dramatically as a result of stress redistribution associated with mining and gas desorption from a large coal volume. These parameters affect the emission rates and ventilation requirements, as well as water inflow into the working environment. Stopping leakage, on the other hand, is a secondary stress dependent factor in estimation of emissions, as convergence of the roof and floor strata, compromising the integrity of the stopping, may result in leakage, making prediction of ventilation requirements difficult. This paper aims to examine the effects of porosity and permeability changes of the coal seam on methane emissions in an underground continuous miner section. The models were developed and executed in a dynamic fashion to simulate an advancing section. Through this process, the changes of effective stress in coal, particularly their change paths, on porosity and permeability were incorporated into the models and methane emissions, concentrations, air requirements, water inflow and possible leakage from the stoppings were investigated using a conventional coalbed methane reservoir model.

Jan 1, 2009

By Vaughan R. Voller

"Over the last ten years an increasing number of heat and mass transfer models of solidification processes have been presented in the literature. A key component in these models is the coupling between thermal and solute fields. This is achieved by developing relationships based on the behavior at the microscopic scale in the casting (typically a secondary dendrite arm space). An important class of models in this respect are models that describe the local scale redistribution of solute alloy components-collectively referred to as microsegregation models. This paper explores this class of models.IntroductionModeling takes a key role in processing materials for properties. This is particularly time in solidification processes where current computer techniques allow for a detailed treatment of the critical macroscopic phenomena. Typically these macro scale models consist ·of three main coupled components.1. A model of flow in the melt and solid + liquid mushy region.2. A model of heat transfer.3. A model of solute transport.The focus of this paper is the solute transport component, in particular the manner in which process scale transport is coupled to transport ·at the local scale of the solid-liquid interface. This step requires a modeling of the redistribution of solutes at the scale of the secondary arm spaces in the dendritic mushy region-a process referred to as microsegregation. Basic microsegregation models use the limit assumptions of no mass diffusion in the solid (Gulliver-Scheil [l]) or complete diffusion in the solid (equilibrium lever rule) in a fixed and space. Recently, however, more sophisticated microsegregation models that account for finite mass diffusion and coarsening of the arm space have been developed [2-19]. In this paper governing equations for a general microsegregation model for a multi-component alloy in a coarsening and spacing _are introduced. From these general equations various models for a binary system are developed and tested"

Jan 1, 2000

By Jan M. Mutmansky, Y. J. Wang

In the mathematical analysis of electrical networks, the electrical devices are often modeled using idealized network elements such as current sources, voltage sources and resistors. An analogous approach is outlined for the modeling of mine ventilation networks, which can be considered nonlinear resistive networks. The proposed modeling technique uses the following five basic (nonreducible) elements: H, Q, D, R and N. These elements are characterized by their pressure gain or loss, their air quantity and their pressure-vs.-quantity relationships, each of which may be known or unknown. In addition to the definition of the basic elements, the necessary conditions for which there exists a unique solution are discussed. This is performed using theoretical developments derived from electrical circuit theory and illustrated in example networks. The proposed procedures are useful in the solution and characterization of mine-ventilation networks.

Jan 1, 1998

By Y. J. Wang

In the mathematical analysis of electrical networks, the electrical devices are often modeled using idealized network elements such as current sources, voltage sources, and resistors. An analogous approach is outlined here for the modeling of mine ventilation networks, which can be considered nonlinear resistive networks. The modeling technique proposed uses five basic (non-reducible) elements: the H element, the Q element, the D element, the R element, and the N element. These elements are characterized by their (a) pressure gain or loss, (b) air quantity, and (c) pressure-versus-quantity relationship, each of which may be known or unknown. In addition to the definition of the basic elements, the necessary conditions for which there exists a unique solution are discussed. This is performed using theoretical developments derived from electrical circuit theory and illustrated in example networks. The proposed procedures are useful in the solution and characterization of mine ventilation networks.

Jan 1, 1996