Search Documents

By Ravi S. Bellur-Ramaswamy, Daniel A. Tortorelli, Robert Haber, Nahil A. Sobh

We optimize continuous quench process parameters to produce functionally graded aluminum alloy extrudates. To perform this task, an optimization problem is defined and solved using a standard nonlinear programming algorithm. Ingredients of this algorithm include 1) the process parameters to be optimized, 2) a cost function: the weighted average of the precipitate number density distribution, 3) constraint functions to limit the temperature gradient (and hence distortion and residual stress) and exit temperature, and 4) their sensitivities with respect to the process parameters. The cost and constraint functions are dependent on the temperature and precipitate size which are obtained by balancing energy to determine the temperature distribution and by using a reaction-rate theory to determine the precipitate particle sizes and their distributions. Both the temperature and the precipitate models are solved via the discontinuous Galerkin finite element method. The energy balance incorporates nonlinear boundary conditions and material properties. The temperature field is then used in the reaction rate model which has as many as 105 degrees of- freedom per finite element node. After computing the temperature and precipitate size distributions we must compute their sensitivities. This seemingly intractable computational task is resolved thanks to the discontinuous Galerkin finite element formulation and the direct differentiation sensitivity method. A three-dimension example is provided to demonstrate the algorithm.

Jan 1, 2001

By S. Zhang, Z. Hua, D. Feng, X. Tong, X. Xie

The separation and beneficiation of low grade collophanite has long been a technical challenge in phosphorous production. It requires the joint efforts of process design, reagents optimisation and operational control. In this study, the statistical response surface method was applied for determining the optimum reagent dosing in the flocculate flotation of a low grade (18.59%) collophanite ore. The dosage of three reagents, including a pH regulator (Na2CO3), depressant (Na2SiO3) and flocculant (modified starch), were input parameters of the modelling design. Two quadratic mathematical models were derived for prediction of grade and recovery based on the statistical design. Statistical analysis of the experimental results validates the reliability of the models. Through further calculation, the optimum dosage of the three reagents were found to be 3,020g/t Na2CO3, 1,990g/t Na2SiO3 and 880g/t modified starch, with which a maximum recovery of 75.81% and P2O5 grade of 30.08% was obtained from flocculate flotation test. Keywords: Process optimisation, modelling, froth flotation, collophanite ore, response surface methodology

Jan 1, 2020

By Ioannis K. Kapageridis

The lignite deposits of the Kozani-Ptolemais-Amyntaio basin in western Macedonia (North-West Greece) present a difficult modelling problem. Each deposit consists of several thin lignite layers ranging from a few centimetres to a few tens of centimeters in thickness, and thin intercalated sterile layers of marly limestones, marls, clays, and sands. These layers are practically impossible to model individually due to fast lateral transition of lignite layers to humus clay and vice versa and the lack of reliable and detailed stratigraphic correlation, leading to compositing methods being applied prior to any interpolation and modelling. The complexity of the deposit is further increased by the presence of a number of faults. This paper presents the modelling and estimation procedures applied during the most recent study of the South Western Field ? a new deposit that is considered for development and lignite production over the next years and will support the operation of a new power station in the area. Special focus is given to drillhole database structuring, fault modelling, stratigraphic correlation, mineable lignite intervals compositing, and resource estimation. The use of specialised mine planning software in all stages of modelling and estimation is thoroughly explained.

Oct 1, 2009

By S. Bellec

A leaching model for gold ore that takes into account the impact of cyanide and oxygen concentrations as well as gold surface exposure on leaching performance is developed. Cyanide consumption is estimated by an empirical model that is a function of ore particles size distribution as well as cyanide, copper and sulphur concentrations. The proposed model has been validated by data from laboratory leaching tests on the Sleeping Giant gold mine ore. Keywords: gold leaching, modelling, exposure surface.

Jan 1, 2009

Modelling and Simulation to Design a Flotation Circuit for Gold Concentration

Sep 13, 2010

By S Hancock, B H. Brady, H Muhlhaus

The aims of this paper are to present the governing science, verify the performance and demonstrate the potential engineering applications of a computational simulator for in situ leach mining of non-ferrous minerals. Assuming that the target mineral is deposited either uniformly within rock matrixes (matrix-hosted) or on fracture walls (fracture-hosted), in situ leaching can be treated as reactive flow and transport in geological media (porous media for matrix-hosted ore deposits and fractured media for fracture-hosted ore deposits). Under this treatment, the processes of in situ leach mining for non-ferrous minerals can be quantified through several sets of equations and these coupled equations have been numerically solved. The computational simulator has been applied to the simulations of one-, two-, and three-dimensional in situ leaching cases. The computational simulator, together with laboratory and in situ leaching tests, can be used to:optimise the leaching chemistry; optimise the in situ scale mining operations; scale up the in situ scale test operation to a full-scale mining operation; and evaluate the post-mining environmental restoration plan.

Jan 1, 2000

By O Handayani, R Barr, E Isaaks

The mining industry has focused a great deal of attention recently on the impact of blast movement on grade control. Although the problem of blast movement and its contribution to dilution and ore loss is well known, solutions to the problem are less than satisfactory. Perhaps the most recent innovation is the application of blast movement monitoring (BMM) devices to quantify blast movement. These devices emit an electronic signal and are planted at surveyed locations within the preblast material. Following the blast, their new locations are found using an electronic scanner. The magnitude and direction of the blast movement is given by the vectors between the pre- and post-blast BMM locations. However, it’s not clear how the BMM vector data can be best used to minimise dilution and ore loss. One common solution is to displace preblast dig line polygons by distances indicated by the BMM vectors. However, studies of BMM vectors show that blast movement not only displaces material laterally but also mixes material internally within the bench (Thornton, 2009a).This paper presents a method for modelling the post-blast muck pile that accounts for displacement and internal dilution using simulated annealing. The preblast ore control block model is re-blocked to 1?×?1?×?1?m sub-blocks. Each sub-block is displaced by simulated BMM vectors conditional to post-blast surface topography and blast initiation sequence data in addition to BMM vector data. The displaced sub-blocks are aggregated into new ore control model blocks whose grades are calculated from the contained sub-blocks. New dig lines are then designed on the post-blast muck pile model using the new ore control block model grades. A case study based on actual data is provided.CITATION:Isaaks, E, Barr, R and Handayani, O, 2014. Modelling blast movement for grade control, in Proceedings Ninth International Mining Geology Conference 2014 , pp 433–440 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 18, 2014

By E J. Sellers, A Cocker

Movement of the fragmented rock mass is an inevitable consequence of open cut blasting. Uncontrolled, unaccounted movement from blasting commonly causes ore loss and dilution from waste mixing. Therefore, the prediction of the post-blast location of both ore and waste is vital in grade control and of direct interest to mining geologists. Direct measurement techniques have matured significantly over the past decade, but their use remains expensive and labour intensive. Access and safety issues are further disincentives to regular movement monitoring. What is needed is a combination of direct measurement to establish baselines and reliable predictive modelling.The W H Bryan Mining and Geology Centre of the University of Queensland has pioneered 3D blast movement modelling using ‘physics engines’ to handle in-flight collisions and predict the final resting place of blasted material. Open cut metalliferous mines are a demanding environment for this style of modelling as they can involve complex geometries and variable buffering, which leads to inconsistent confinement. A version of the model designed to deal with these complexities is under development. A well-controlled field trial was needed to validate this new version.A validation trial was completed at a gold mine in South-East Asia. The aim was to increase throughput by increasing the blast powder factor without losing grade through dilution and mixing with waste. Four blocks were blasted to identify the effect of bench height and face confinement on ore movement, fragmentation and process throughput. The powder factor and blast timing were kept constant. Movement of the broken rock was measured by placing blast movement monitors at known positions in the block and then locating them after the blast. The majority of sensors were recovered, providing a unique insight into blast movement dynamics and generating a substantial data set for model validation.The model (called ‘MetaMove ‘) was calibrated for one blast and applied to investigate the extent of grade loss and possible mitigation strategies. Conventional mining techniques that involved firing across the strike of the orebody, and no post-blast boundary adjustment of the ore/waste boundaries prior to digging, were found to cause up to 44 per cent loss of high-grade ore, depending on the grade distributions. Post-blast adjustment of the ore/waste boundaries can negate this loss. The model indicates that altering timing and ensuring that the block boundaries are well confined can further constrain the movement of ore/waste boundaries. Thus, the model can be successfully applied to minimise grade dilution and limit misclassification of ore as waste whilst using higher powder factors to increase mill throughput. Application of a calibrated model to move ore/waste boundaries prior to digging will enable the production geology team to add significant value to the mining operation.CITATION:Cocker, A and Sellers, E J, 2014. Modelling blast movement for grade control at an open cut gold mine, in Proceedings Ninth International Mining Geology Conference 2014 , pp 377–386 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 18, 2014

Tonnage and grade are the principal reported components of mineral resource estimates but tonnage is a derivative of two more fundamental estimates; volume and bulk density. Typically, most of the effort in a resource estimation study is directed towards grade and volume, yet poor application of bulk density data to resource models can result in significant biases. The impact on tonnage estimates is easily recognised but poor density estimates may also introduce subtle biases into the grade estimates. Methods for measuring bulk density have been discussed previously. This paper focuses on the manipulation of bulk density data and the extension of bulk density from essentially point measurements to block estimates within resource models. Bulk density is a continuous, geologically controlled, spatial variable and should be treated as such. Issues such as combining data from multiple sources, sample support and the interdependence of bulk density and grade estimates must be resolved. The statistical and geostatistical techniques used for grade estimation are, under the right circumstances, equally applicable to bulk density. The importance of bulk density needs to be recognised. With robust raw data and correct estimation procedures, better estimates of bulk density can be obtained and both local and global tonnage and grade estimates can be improved.

Jan 1, 2000

By N. Shi, T. G. Joseph

ABSTRACTA series of algorithms were developed to perform geometric and kinematic analyses of cable shovel dipper-ground interactions for a given set of mining face conditions. The routines examine and analyze (a) the position of the dipper in relation to the shovel; (b) the dipper in relation to the face conditions; (c) the dipper position in relation to possible interference or collision with crawler tracks during digging cycles; and (d) kinematic and geometric considerations with respect to any point on the dip-per as it conducts a duty cycle. In all cases, field measurement verifcation provided good agreement(within 4% variance).

Jan 1, 2011

Modelling Channel Scale Flow in Heap Leaching Using Smooth Particle Hydrodynamics (SPH)

Sep 13, 2010

By K Williams, T Donohue, W Chen, T Roessler, A Katterfeld

Wet sticky iron ore poses particular difficulties in every aspect of the material handling process. Material hang-ups in bins, pipes and chutes; carry-back and spillage due to insufficient belt cleaning; clogging of screens; malfunctioning of crushers as well as rat holing of gravity reclaimed stockpiles are all predominantly due to the inherent cohesion and adhesion characteristics of the wet, sticky ore (Roberts, 1993). Limited understanding and quantitative measures on the stickiness of such a material significantly lower the mining operation efficiency.This research aims to investigate the cohesion and adhesion characteristics of the wet, sticky ore on physical modelling, experimental techniques and numerical implementation. Initially, the suitable bridging mechanism that contributes to the cohesive and adhesive forces within the iron ore material was identified (Hering, Martin and Stohrer, 1989; Rumpf, 1958). A cohesion/adhesion energy density factor based contact model was then introduced to describe the stickiness of the material. Secondly, an experimental system was proposed to directly test the stickiness of the material. Such a system is able to obtain both the cohesive/adhesive shear force and tensile force of the material at a specific moisture content. Lastly, the implementation of the proposed cohesion/adhesion model into the discrete element modelling (DEM) software was also discussed. A cohesion/adhesion calibration for modelling wet sticky iron ore in DEM was then conducted to ensure the material modelled in the simulation is reflective of the reality.From the study, understanding of the cohesion/adhesion characteristics of the wet iron ore material is significantly enhanced. Experimental and numerical outcomes can be directly applied for better designing or optimising the material handling processes that involve wet sticky iron ore. Consequently, the mining operating costs can be greatly reduced.CITATION:Chen, W, Donohue, T, Williams, K, Katterfeld, A and Roessler, T, 2015. Modelling cohesion and adhesion of wet, sticky iron ores in discrete element modelling for material handling processes, in Proceedings Iron Ore 2015, pp 387–392 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 13, 2015

Traditionally within the mining industry, single models for both grade and geology of orebodies are created upon which all mine development decisions are based. These models provide a single interpretation of the extent and continuity of the mineralization envelope based on solids and sections interpreted from relatively wide spaced drilling i.e. the models infer behaviour of grade and geology at smaller/local scales based on data available at wider scales. The variable behaviour of grade and geology cannot be derived from single kriged resource models. In recognition of this, Newmont Mining Corporation typically uses Sequential Gaussian Simulation (SGS) to quantify grade variability and to determine its effect on some of these mine development decisions. In most of these studies the mineralized envelope remains fixed and grade simulations are generated within a fixed volume. To account for uncertainty in the geology/ mineralization envelope Newmont Mining Corporation has been experimenting with Multiple Point Statistics (MPS). Multiple-Point Statistics is an emerging technique, which can be used to generate multiple, geologically realistic, realizations of data that displays complex non-linear features. MPS uses a conceptual model of the geology, called a Training Image (TI), to infer these high-order spatial relationships. The multiple realizations allow for an assessment of the geologic uncertainty, which can be further combined with grade simulations to generate a more complete picture of the true uncertainty of a mineral deposit.

Jan 1, 2011

By S. T. Johansen, C. Colomb, B. Ravary

Combustion and thermal NOx formation have been modelled using the Computational Fluid Dynamics (CFD) commercial code Fluent 6.1. The aim was to predict and give a better understanding of the combustion in electric arc furnaces for the production of ferro-silicon and silicon metal. Another important objective was to assess the impact of a number of operating parameters on the NOx formation. The CFD tools used are de-scribed with a special focus on models for species and turbulence transport. Sensitivity studies have been carried out to study the influence of several parameters: sub-models, process gas composition and distribution at the charge surface and influence of the presence of SiO. The main conclusions are: ?Using the combination of a realistic design, boundary conditions and models, the simulated NOx level is in quantitative agreement with measurements. The presence of water vapor is particularly important because of its influence on the specific heat of the mixture and resulting temperature. ?Release of significant amounts of SiO from the charge surface results in an important increase of NOx, even if we applied a simplified combustion model for SiO gas. However, because of lack of thermal and compositional data for blowing situations (when SiO flows directly from the arc cavity through open channels to the charge surface), no quantitative validation is possible. The way SiO is distributed (uniformly or as a jet) has a limited impact on the result in our computations.

Jan 1, 2007

This paper presents a case study where Cenozoic age detritals and an ancient buried landslide have been modelled in the final walls of an iron ore pit. Orebody 25 is located in Western Australia’s Pilbara region and is an operating pit where a redesign of the final wall cutback was required. The existing slope design utilised a simplistic geological model of mixed conglomeratic, calcretic and clay –dominated detritals overlying deformed bedrock sequences. This simplicity resulted in a conservative design that was not reflective of the slopes performance to that point. Further field investigation was carried out with the results assessed in conjunction with historic drilling, mapping and laboratory testing. A detailed geotechnical model was developed including three conglomerate and two clay dominated units and laterstage calcrete overprinting. The clay units were of significant importance to the stability assessment and resulting slope design due to the thickness and poor strength conditions indicated by the laboratory and field testing. However, the drilling investigation shed only limited understanding as to the origin of the clays and thereby how they should be treated as a component in the stability analysis. A comparison of the ratio of various geochemical elements measured by assay was made between the clays and 1) shale bedrock and 2) clays of known detrital origin from other nearby sites. This indicated the clays had a strong affinity with weathered bedrock. This and other textural features resulted in the interpretation that the clay units were relic debris preserved following ancient landslides in bedrock shale. Using laboratory testing data and calibration of results for a landslide origin, stability analyses were run, resulting in a steeper slope design being implemented. Excavation of the final cutback is now in progress with the model confirmed as appropriate.

Jan 1, 2015

By T Ma, F T. Suorineni, Q He, J Oh

Directional hydraulic fracturing was invented by Polish experts and was first applied in coal mining to prevent hard roof rock bursts. Based on its successful application in coal mining, directional hydraulic fracturing has the potential to reorientate the hydraulic fracture trajectory. Recently, this technique has been introduced into cave mining for the purpose of lowering breakdown pressure. The hydraulic fracture is much easier to initiate with the assistance of an initial notch. No systematic theoretical or experimental studies have been carried out to investigate the influencing factors in the hydraulic fracture reorientation process. In this paper, the numerical modelling software Realistic Failure Process Analysis3D (RFPA3D) – Flow Version is used to study the influence of the in situ stress condition and rock mass heterogeneity on directional hydraulic fracturing in brittle heterogeneous rock masses. Firstly, the reliability of RFPA3D in simulating hydraulic fracturing is validated by comparing its results with experimental results and the simulation results from the extended finite element method. The hydraulic fracture trajectory is then observed in a three-dimensional stress environment. The results indicate that for the in situ stress condition, when the borehole lies in the ?2–?3 plane, only the difference between ?2 and ?3 has an important influence on hydraulic fracture reorientation. Furthermore, rock mass heterogeneity is a key factor that must be accounted for when fracture reorientation is required.CITATION:He, Q, Suorineni, F T, Ma, T and Oh, J, 2016. Modelling directional hydraulic fractures in heterogeneous rock masses, in Proceedings Seventh International Conference and Exhibition on Mass Mining (MassMin 2016), pp 369–384 (The Australasian Institute of Mining and Metallurgy: Melbourne).

May 9, 2016

By S Bhandari

A model for predicting dust plume dispersion from blasting operations using the Eulerian mathematical algorithm has been developed. This model considers atmospheric stability, and wind velocity and direction for computing dust concentrations at different distances from the blast. Software has been developed which simulates the movement of a dust plume at a given horizontal distance from the blast location and time elapsed after the blast. The concentration of dust is also given in vertical planes. It takes into account ground profile and atmospheric conditions such as temperature, pressure, humidity, wind velocity, wind direction and shear. Experiments were carried out in two limestone mines to calibrate the accuracy of the predictions made by such model. High volume dust samplers were used to collect dust samples. Difficulties encountered in satisfactory dust sampling and quantification of results have been described.

Nov 8, 2011

By Y Wei, H-B Mnhlhaus

Particles are by far the easiest objects to simulate, because their dynamics is simple. Determining the motion of particles and their interaction, a wide variety of deformable structures can be constructed by connecting particles that have mass, position, and velocity, and that respond to forces. Starting from this simple idea we have derived what we call the particle-lattice method for the solutions of boundary value problems of deformation and fracture of solid mass. The problem of rock breaking by blasting operations as well as fracture and flow of rocks fragments involve momentum transfer between solid and fluid; gas flow and reactive mass transfer. In this paper we outline the method and describe briefly the salient features of the algorithm. We derive relationships between the contact properties for the discrete model and the properties of the underlying continuum and fluid flow coupling procedure. Finally we illustrate the potential of the method by means of example solutions related to the dynamic fracture and flow in jointed rock by blasting.

Jan 1, 1999

By P. Lundman, L. Rosengren, T. Brandshaug

Tunnel 99’, the Swedish National Road Administration´s (SNRA) general technical specification for road tunnels requires that the load bearing system must be designed to withstand specified design explosion loads. This paper presents the results of an SNRA study to provide a perspective of design analyses tailored to comply with the new requirement. The paper provides an example of a design analysis for assumed tunnel conditions using a numerical model. Dynamic predictions of potential damage to the rock mass and ground support system (rock bolts and shotcrete) from the specified explosion loads are presented including the nonlinear response of the shotcrete. “Tunnel 99”, qui fixe les règles de l’Administration nationale suédoise des Routes (Vägverket) concernant la conception et la construction des tunnels routiers, requiert entre autres que la structure porteuse soit capable de résister à certaines contraintes spécifiques résultant d’explosions. Cette étude présente les résultats d’une étude effectuée par Vägverket en vue d’obtenir une perspective d’analyses de conception spécialement élaborées pour se conformer aux nouvelles règles. Il contient un exemple d’analyse en fonction de conditions déterminées, basée sur un modèle numérique. Des prévisions de’études dynamiques de dommages potentiels subis par la masse rocheuse et les systèmes de soutènement (boulons d’ancrage et béton projeté) du fait des contraintes libérées par la force d’explosion spécifiée sont également présentées, y compris dans le cas d’une réaction non linéaire du revêtement en béton. “Tunnel 99” - die Vorschriften für Strassentunnel der Nationalen Schwedischen Behörde für Straßenwesen - (Vägverket) - fordern vom Tragsystem eine Auslegung in der Form, dass es spezifischen Explosionsbelastungen standhält. Dieser Artikel präsentiert die Ergebnisse einer Studie der Behörde, um einen Ausblick über konstruktive Auslegungen zu geben, die geeignet sind, die neuen Vorschriften zu erfüllen. Der Artikel zeigt ein Beispiel einer konstruktive Auslegung für eine angenommene Tunnelsituation mit Hilfe eines numerischen Modells. Voraussagen unter Beachtung der vollständigen dynamischen Berechnung zu potentiellen Schädigungen des Gebirges und des Ausbaus (Anker und Spritzbetonauskleidung) durch die spezifizierten Explosionsbelastungen, einschließlich des nichtlinearen Verhaltens der Spritzbetonauskleidung, werden gezeigt.

Jan 1, 2003

By R. P. Jensen

This is a sampling of the papers presented at the conference, the proceedings can be bought from the SAIMM Arbitrary Lagrangian Eulerian (ALE) computational techniques allow treatment of gases, liquids and solids in the same simulation. ALE methods include the ability to treat shock waves in gases, liquids and solids and the interaction of shock waves with each other and with media from one of the other categories. ALE codes can also treat explosive detonation and the expansion of the explosive gases and their interaction with air and solids. ALEGRA is a 3-D ALE code that has been developed at Sandia National Laboratories over the past few years. ALEGRA has been applied to a 2-D simulation of presplitting using decoupled explosives in rock blasting with very interesting results. The detonation of the explosive at the bottom of the hole sends a shock wave up the borehole driven by the explosive gas expanding into air. The explosive gas compresses the air against the stemming column where it rebounds and recompresses at the bottom of the borehole. This type of ringing takes several cycles to damp out. The explosively induced expansion of the borehole is also treated by ALEGRA as well as the shock wave imparted to the rock. The presentation of this paper will include several computer animations to aid in understanding this complex phenomenon.

Jan 1, 2000