Search Documents

By Norman Paley, Zachary Pickett

The Red Dog Mine is a high-grade zinc-lead mine located in northwest Alaska which began operations in November 1989. In June 2014 the upper half of a charge in a trim shot in the Aqqaluk pit deflagrated. No other holes in the area were observed to react and the remainder of the pattern was shot without incident. This was the only known case of a deflagration in the 25 years of operations. As part of the investigation into the cause of this incident rock samples were obtained from the area of this hole and from other areas in this pit. The samples were assayed and tested for reactivity using the Reactive Ground Catalyzed Adiabatic test. Approximately half of the samples were found to be reactive with standard blasting agents. Unexpectedly, some of these were found to be reactive with blasting agents doped with urea to act as an inhibiter. Relationships between element grades and reactivity and between Specific Heating Capacity and reactivity were developed; however, the small sample size led to these models having an unrealistically large number of geological model blocks tagged as reactive. A program to better understand the behavior by testing more samples was initiated. Drill core was selected based on the predicted level of reactivity using the aforementioned relationships. The samples covered a greater range of grades and lithologies and included rock from the soon to be developed Qanaiyaq pit. The samples were dry pulverized then split for elemental assay, Reactive Ground Catalyzed Adiabatic testing, Specific Heating Capacity (SHC) determination, and X-ray Diffraction (XRD) mineralogical analysis. From this data, a Red Dog Mine specific assay based classification for separating reactive from non-reactive ground was developed. Furthermore, it was possible to divide the reactive ground into two subgroups: i) rock which reacts with ANFO-emulsion blends even when manufactured specifically with a reactivity suppression agent like urea and ii) rock which reacts with ordinary ANFO-emulsion blends only when a reactivity suppression agent like urea is absent.

By Norman Paley, Zachary Pickett

The Red Dog Mine is a high-grade zinc-lead mine located in northwest Alaska which began operations in November 1989. In June 2014 the upper half of a charge in a trim shot in the Aqqaluk pit deflagrated. No other holes in the area were observed to react and the remainder of the pattern was shot without incident. This was the only known case of a deflagration in the 25 years of operations. As part of the investigation into the cause of this incident rock samples were obtained from the area of this hole and from other areas in this pit.

Feb 1, 2020

By J S. Kuszmaul

"A computational model of rock blasting is being developed to examine the blasting problems associated with in situ oil shale processing. This model, however, will also be useful as a design tool for the traditional problems in rock blasting. The model includes fundamental treatment of both shock-wave propagation and the accumulation of brittle fracture in the rock. As a result, the model accurately predicts the degree and extent of fracturing as functions of design parameters. The model has proven useful for making parametric studies and for evaluation of alternate blast designs. This paper demonstrates the use of the numerical model to simulate the fracturing induced by the detonation of a vertical explosive column near a bench. The fracturing induced by three differentexplosives indicate that (in the chosen geometry) the most efficient breakage is done by a column of ammonium nitrate and fuel oil mixture (ANFO) used with a toe charge of aluminized ANFO. There was too much unfractured rock left when ANFO was used alone; aluminized ANFO used for the entire explosive column caused excessive fracturing. A final case involves ANFO used alone to fracture a different rock type. This case points out that in a different rock type, the ANFO will not leave excessive unfractured rock."

Jan 1, 1987

By Danqing Gao, Brijes Mishra

Spatial variance of rock properties is an important characteristic of rock which addresses the heterogeneity of mechanical and physical properties and their effect on underground mining. Assessing spatial variance can be useful for both locating potential difficult grounds and performing reliability analysis in the mining support system. This study shows a microscale numerical model for rock material based on continuum mechanics. To consider the spatial correlation factor of mechanical properties in realistic rock mass, this study introduces a spatial correlation concept into the traditional type III extreme value distribution model. Based on the improved microscale random model, we conducted a series of numerical uniaxial compressive tests and one longwall mine model to investigate the progressive failure of rock. The results show that spatial correlation factor can significantly affect the load-deformation curves and fracture patterns. The spatial variance model was applied to a longwall mine that included several cutting sequences to observe the influence of the spatial variance on roof behavior. The cohesion and friction random fields induced failure that propagated randomly through the immediate roof in different entries. The results showed that the spatial variance provides an accurate prediction of erratic roof falls in coal mines.

Jul 20, 2020

By Sidney L. Hansen

This report was written because of concerns of the possible error involved when stopping leakage is incorrectly assumed to be turbulent. It is common practice to use turbulent flow to characterize each airway branch for computer models of mine ventilation circuits. This report addresses computer methods to accurately model air-leakage paths. A graphical method is described for deriving the equations of leakage through stoppings. The significance of modeling leakage according to the proper flow regime (turbulent or laminar) is demonstrated by the following method. A hypothetical mine was developed to represent present-day conditions. The stoppings of the mine were chosen to leak according to laminar flow properties. The mine was represented by constructing two base models. One base model used turbulent leakage equations for each stopping, while allowing the corresponding resistances to vary. The other base model used the correct laminar leakage equations and held corresponding resistances constant. The mine circuit was then expanded to demonstrate future development. The expanded portion was designed with stoppings that were identical to those in the base mine. The mine in the future state was represented by two projection models. Each projection model was an extended version of its associated base model; one used turbulent leakage equations and one used laminar leakage equations for stoppings. The laminar projection was theoretically correct and is used as a basis for comparison of differences in airflow quantities. Rather than modeling every leakage path for stoppings, a more practical approach is also discussed. A good level of accuracy is maintained by combining leakage paths if the proper flow equations are used. This method was used on the hypothetical mine projection and the flows were very similar to the theoretical values.

Jan 1, 1989

By Shoma Sinha, Nirmal K. Sinha

Deformation in a material at any temperature is known to consist of two major components - a reversible (elastic) part and an irreversible part (often called plastic). What is not universally known is the fact that at high temperatures greater than about 0.4Tm, where Tm, is the melting point in Kelvin, the reversible part includes not only elastic but also a delayed elastic (recoverable but time dependent) component. A very simple computational model is presented here for a constant-strain stress relaxation test (SRT) on the basis of a three-component constant-stress rheological equation derived from strain relaxation and recovery test (SRRT) data. It is shown that the rapid stress reduction during the early periods of relaxation when the stresses are high is governed primarily by delayed elasticity irrespective of the imposed strain level. This is in contradiction to the popular concept prevailing today. The asymptotic approach to a quasi-stable stress rate at longer times and lower stresses is, however, controlled by viscous flow (power-law creep). The validity of the theory is examined by comparing the predictions with experimental results obtained for a titanium-base alloy Ti-6246 at 0.45Tm.

Jan 1, 2004

By André C. Zingano

The domino effect in a violent pillar collapse occurs due to the redistribution of stress in the remaining pillars near the collapse, and this stress transfer occurs in the immediate and main roofs. The remaining pillars do not have enough load capacity to support this overstress and also fail, and the domino effect propagates further to the entire panel. The domino effect does not occur only due to low load capacity of the pillars, but also the roof rock mass quality, because the stress is transferred by the roof. In high competent roof without discontinuity (fault), like massive sandstone seam, the deformation of the roof is the cause of pillar collapse extension and propagation. Depending on how close it is approaching the collapse limit; the immediate and main roofs behave differently. When it is close to the collapse limit roof convergence and pillars compression occur. Conversely when it is far from the collapse limit divergence and pillars decompression occur. This is called the arc-effect. This paper studies this behavior of the roof (immediate and main roofs) using numerical modeling to simulate the roof deformation and stress distribution in the roof and pillars. This simulation represents the progression of pillar collapse and shows the zone of convergence and divergence of the roof, and also the effect of this behavior on pillar confinement. Two-dimensional plain strain model was applied to simulate the arc-effect. The simulation results were compared to those of the instrumentation that was installed before the collapse. Keywords: pillar collapse, room and pillar, numerical modeling

Jan 1, 2005

By S. Sundarraj

The morphological transition between columnar and equiaxed structures is an important solidification phenomenon frequently observed in castings, especially complex parts with varying cross-section, The mechanical properties and the occurrence of defects such as hot tears are strongly dependent on the structure formed in the casting. Hence, it is important to understand the conditions leading to columnar or equiaxed structures. Although extensive work has been done in modeling solidification systems, almost all of the existing models are based on the assumption that the structure formed in the casting is either purely columnar or purely equiaxed. Prediction of the columnar-equiaxed transition (CET) requires a solidification model which accounts for the formation of both columnar and equiaxed structures, along with related microscopic and macroscopic solidification phenomena. In this paper a CET model is presented which: (i) takes account of both columnar and equiaxed structures, (ii) accounts for macroscopic heat transfer coupled with microscopic phenomena such as nucleation, growth and microsegregation, (iii) satisfies solute conservation, and (iv) is easy to implement in commercial finite element and finite difference software. The proposed model predictions of CET have been validated with available experimental results for a uni-directional solidification of an aluminum-copper alloy. A sensitivity analysis has also been carried out to determine the effects of various parameters on the CET model predictions. The model is being extended to other alloy systems and more complex geometries. This work will be published in due course.

Jan 1, 1996

By Mohsen Mohammadi, Edward Cyr

"To date, several studies have been performed on the static mechanical properties of SLM-AlSi10Mg alloy. There have also been more recent studies on the relation between scan strategy, texture, and mechanical properties of SLM metals and composites. However, there is a huge knowledge gap in the available literature regarding texture and mechanical behavior of SLM-AlSi10Mg alloy under high strain rates. Therefore, in this study, dynamic mechanical behavior of AlSi10Mg_200C alloy manufactured by DMLS technique was investigated at high strain rates ranging from 900/s to 1700/s using Split Hopkinson Pressure Bar. Texture using X-ray diffraction technique of the initial samples were measured. The yield strength, peak flow stress, and ductility increased in all samples with increase of strain rate. TEM studies have also revealed evidence of continuous dynamic recrystallization in grains after 900/s–1700/s strain-rate tests. A Taylor polycrystal crystal plasticity constitutive model was then developed to capture and predict the texture evolution and mechanical response under dynamic loading, using plastic work as a critical parameter for the onset of recrystallization.INTRODUCTION Additive manufacturing (AM) powder techniques are among the fastest advancing breakthroughs in manufacturing in recent years, defined by a bottom-up layer by layer process of joining with energy input from lasers or electron beam (Herzog et al., 2016; Gibson et al., 2015). Among alloys fabricated by SLM, AlSi10Mg is of high interest to industry and academia due to its low density, high specific strength and excellent corrosion resistance, making this alloy a viable contender for aerospace, automotive, and marine applications. Circumstances such as car crash or bird aircraft strike hazards subject potential additively manufactured AlSi10Mg parts to impact loadings. This drives an interest in understanding and predicting dynamic behaviour of these alloys. AlSi10Mg is the AM counterpart of A360 die-cast aluminum, which possess hypoeutectic microstructure (Li et al., 2015). Up to the present time, many studies have presented the mechanical properties of SLM-AlSi10Mg alloy. The mechanical properties and microstructure of AlSi10Mg produced by selective laser melting (SLM) or DMLS techniques can be controlled using different parameters, such as feedstock powder (Asgari et al., 2017), and process parameters (Olakanmi et al., 2015). The microstructure and mechanical properties of AlSi10Mg can also be controlled using post heat treatment after SLM process Li et al., 2016). The mechanical properties of SLM-AlSi10Mg are sensitive to strain rate and can both change under quasi-static (Rosenthal et al., 2017) and high strain rates (Asgari et al., 2018). However, a significant knowledge gap in available literature regarding modeling these behaviours, and especially micromechanical modeling of dynamic behaviour."

Jan 1, 2018

By Liming Yuan

Barometric pressure changes affect air density, leading to change in the mass of the gas in the gob. When the barometric pressure decreases, the volume of gas in the gob expands, while the volume of gas contracts when the barometric pressure increases, causing the gob to breathe out and in. Although the concept of gob ?breathing? is simple, its effect on spontaneous heating of coal in the gob area is not clear. In this study, computational fluid dynamics (CFD) simulations were conducted to model the spontaneous heating of coal in longwall gob area under measured barometric pressure changes. A single longwall panel using a bleederless ventilation system was simulated. If there is no barometric pressure change, the intake airflow rate is equal to the return airflow rate. When the barometric pressure changes, these two airflow rates are no longer equal, and the difference between the two airflow rates represents the airflow rate the gob breathes in and out. The effect of inflow and outflow of gas on the potential spontaneous heating in the gob was investigated using the CFD model developed in our previous study. The effect of barometric pressure changes on the spontaneous heating was found to be dependent on the gob permeability and the coal oxidation rate. The effect of barometric pressure changes on oxygen concentrations in the gob was also examined.

Jan 1, 2010

By Liming Yuan, A. C. Smith

Barometric pressure changes affect air density, leading to change in the mass of the gas in the gob. When the barometric pressure decreases, the volume of gas in the gob expands, while the volume of gas contracts when the barometric pressure increases, a process which causes the gob to “breathe” out and in. Although the concept of gob “breathing” is simple, its effect on spontaneous heating of coal in the gob area is not clear. In this study, computational fluid dynamic (CFD) simulations were conducted to model the spontaneous heating of coal in the longwall gob area under measured barometric pressure changes. A single longwall panel using a bleederless ventilation system was simulated. If there is no barometric pressure change, the intake airflow rate is equal to the return airflow rate. When the barometric pressure changes, these two airflow rates are no longer equal, and the difference between the two airflow rates represents the airflow rate the gob breathes in and out. The effect of inflow and outflow of gas on the potential spontaneous heating in the gob was investigated using the CFD model developed in our previous study. The effect of barometric pressure changes on the spontaneous heating was found to be dependent on the gob permeability and the coal oxidation rate. The effect of barometric pressure changes on oxygen concentrations in the gob was also examined.

Jan 1, 2010

By Gabriel Walton, Sankhaneel Sinha

"Rib failures are a major hazard in underground mining. When pillars are subjected to loading, extensile cracks are generated parallel to the excavation boundary, which can manifest as rib falls. A common mitigation measure is to install supports, especially bolts, normal to the rib surface. The type and spatial distribution of these supports along the pillar is often based on experience rather than comprehensive analysis. Since conducting field experiments is often impractical and expensive, a convenient alternative is to perform 3D numerical modeling. In this study, FLAC3D has been used to capture the interaction of rib bolts and the failure processes governing the stability of pillars. An associated goal is to determine whether, at all, a continuum software is capable of modeling the effect of supports. The modeling involved two important steps: selection of a proper constitutive model, and, explicit simulation of rock bolts. Based on precursory work in the last decade, an S-shaped strength criterion was chosen that can describe the tensile fracture mechanism of rock at low confinement and shear failure mechanism of rock at high confinement. The built-in pile structural element in FLAC3D is a simplified representation of the complete bolt system. Hence, the different components of a rock bolt (i.e., steel, grout, steel-grout interface, and grout-rock interface) have been explicitly modeled within the rock mass and then tested for its ability to replicate the ground-support interaction behavior. INTRODUCTION Rib failure-related hazards and injuries have plagued the mining industry for decades (Mohamed et al., 2016; Iannacchione and Prosser, 1998). In highly stressed environments, excavation parallel cracks are the primary initiators of skin instability. In absence of any support, this failure can indirectly destabilize the excavation by increasing the exposed span of intersections and roadways. Over the last two decades, there have been numerous injuries (fatal and non-fatal) due to rib failures in underground coal/non-coal mines."

Jan 1, 2017

By Mick Lownds, Ruilin Yang

The potential advantages of electronic detonators with accurate timing have been increasingly recognized by the blasting community. However, in comparative blasts the results with pyrotechnic timing are sometimes comparable to those from electronic timing. Statistical variation in blast results is a plausible cause of these observations.

Jan 1, 2011

By Marianella Perez, Williams Calderón, Gonzalo Montes-Atenas

The lack of fresh water sources, particularly in the Northern part of Chile, has driven the mine sites to search for new and sustainable paths to obtain water resources to run minerals processing operations. Among them, the use of seawater in concentrators appears as an interesting solution to this problem. However, the presence of seawater will have an impact onto the operation, not only because it will enhance corrosion reactions in the equipment installed in a plant, but also because it may affect how particles will effectively be processed across the numerous unit operations existing in the mine site. Among such unit operations, the grinding stage is a highly regarded as it defines the success of the separation of valuable minerals from gangue material in the froth flotation operation. Then, one important aspect to be addressed is to understand how the presence of seawater will affect the different mechanisms acting during the operation of a ball mill. The aim of this research work is to assess the effect of changing the aqueous matrix from process water to seawater in ball mills onto the mill performance. The methodology considered the use of Fluent® software and its capabilities to run numerical simulations with heterogeneous media. On one hand, from the data analysis, it is possible to infer that breakability mechanisms should not be significantly affected. On the other hand, the high ionic strength of the matrix should impact importantly particle transport mechanisms within the mill. Particularly, when transporting finer particles. Rheology results were in agreement with previous published research studies suggesting the suspensions to behave as Newtonian, Dilatant o Pseudo-plastic fluid depending on the amount of fine particles present in the grinding media. The latter will result in changes in the residence time of particles in the mill producing overgrinding or an excess in the circulating loads.

Jan 1, 2014

By L. M. Tavares, H. A. Petit, P. P. Cavalcanti

Iron ore pellets experience degradation during handling due low energy impacts and stresses. These impacts lead to volume or surface breakage of the pellets, creating fine particles that accumulate along the process. Some of the particles are in contact with the air and are dragged as dust into the workplace or into the atmosphere. The aerodynamic size of the particles allows them to be potentially inhaled by the workers or to be carried to nearby communities. Dust emission becomes critical at specific locations, such as transfer chutes and stockpiles. This work models the dust generation that occurs in an iron ore pellet stockyard using a methodology that combines experiments to quantify the degradation of the pellets to mathematical modelling of the dust emission mechanisms in the stockyard. First, laboratory tests were used to characterise the amount of mass and the particle size distribution of the dust generated by degradation under controlled conditions. Second, the flow over a stockpile of an iron ore pellet company was simulated with the aid of computational fluid dynamics (CFD) simulations. Third, semiempirical equations linked the data generated by simulations and laboratory tests to estimate the dust emission in the stockyard. Finally, results were validated with measured field data, obtained with a novel online measurement tool. Results suggest that the approach is able to predict dust emission with sufficient accuracy and can be used for further implementation in the management of operations and decision making in order to reduce the environmental impact of the process. Keywords: Iron ore pellets, abrasion, dust emission, mathematical modelling

Jan 1, 2020

Gold cyanide equilibrium adsorption isotherms in the loading range typical of operating carbon-in-pulp/carbon-in-leach plants have been established and the effect of various cations and the presence of copper cyanide species quantified. It was found that the equilibrium loading of gold was enhanced by the presence of sodium or calcium, with calcium having an order of magnitude stronger effect. The presence of copper has an adverse effect on the equilibrium loading of gold. A unified model was developed to predict the effect of sodium and calcium concentrations in solution and the competitive carbon adsorption of copper on the equilibrium loading of gold on carbon. A satisfactory agreement between the measured gold equilibrium loadings and the model results was obtained over a wide range of solution conditions.

Jan 1, 2010

By W. -S. Hwang, R. A. Stoehr, P. Ingerslev

"The flow of molten metal entering complex two-dimensional foundry molds has been modeled using the Simplified Marker and Cell (SMAC) and Solution Algorithm (SOLA) techniques. The molds were for a series of four and three-spoke wheel designs. The computer calculations were compared to physical experiments in which molten cast iron was poured into C02-bonded sand molds. The metal was photographed with a high speed motion picture camera focused through a Pyrex glass window on one face of the mold. To get good agreement between the calculations and the experiments, it was necessary to accurately represent the velocity through the ingate. This was calculated by application of the Bernoulli equation to flow through the feed system. Some difficulties were experienced with mathematical instability, expecially in the three-spoke design with a narrow ingate. This was overcome by using 100% upwind differencing, and over relaxation parameter near 1.0 (i.e. no over relaxation), and variable mesh spacing.IntroductionModeling the flow of molten metal entering molds has the potential for making significant contributions to the design of feed systems. It will supply information to the computer-aided-design systems of the future that will enable the designer to optimize the placement of gates, dimensioning of passages, and designation of casting speeds. Such design criteria have generally been based on application of energy balance (Bernoulli equation) principles to flow in the sprues, runners, and gates. Indeed, this is the easiest way to determine the proper dimensions to give a suitable filling time for a casting of a given size. A number of computer programs, including one by some of the present authors, have integrated such calculations into CAD systems. l Modeling the flow of metal inside the mold cavity requires other techniques that will handle free-surface transient flow. Stoehr and Hwang reported on the application of the Marker-and-Cell (MAC) technique to mold filling problems at the previous conference in this series at Hennicker, N.H. in 1983. 2 Although considerable progress has been made since that time in improving the accuracy, stability, and convenience of these techniques, the shapes to which they were applied remained limited. Also, physical verification of the results of the calculations was limited to water models.The work to be reported in this paper was undertaken with two objectives: (1) to apply these techniques to more complex and realistic shapes, and (2) to verify the results by comparison with flow patterns observed in molten metal poured into actual sand molds."

Jan 1, 1986

By M. G. Lipsett

One challenge in applying information technologies in Mining is building a business case. In Mine Maintenance, there is a complex flow of information that contributes to making informed business decisions. By understanding those information relationships and the cost of wrong decisions, analysts can estimate the value of information flow. Developers can then ensure that systems deliver key information in an accurate and timely way. In this paper, general maintenance decision-making is discussed, with examples of how model-based analysis can improve the information delivered to decision makers.

May 1, 2001

In low carbon TRIP steels, the formation of carbide-free bainitic ferrite is crucial for the stabilization of austenite. Based on experimental investigations modeling the bainite reaction is examined adopting diffusional and displacive mechanisms, as well as the semi-empirical JMAK approach. Transformation tests were carried out to study the isothermal austenite decomposition into bainitic ferrite for a model steel containing 0.6wt%C, 1.5wt%Mm and 1.5wt%Si which replicates the composition of remaining austenite at the start of the bainite reaction in a classical C-Mn-Si TRIP steel. In the diffusional approach, classical nucleation theory is employed and subsequent growth of bainite plates is modeled using the Zener-Hillert theory. In the displacive approach, the autocatalytic nucleation of bairtite plates is assumed to be rate controlling. The predictive capability of all models is evaluated for temperature ranges where different morphologies of bainitic ferrite have been detected, The shortcomings and inherent challenges associated with either of these modeling approaches will be delineated.

Jan 1, 2003

By Daniela Campos, Aldo Casali, Gonzalo Montes-Atenas

The understanding of the effect of the major variables influencing the froth flotation performance may allow both to predict and to identify best operating conditions towards optimizing the process efficiency. The aim of this research work is to model the effect of dissolved oxygen onto froth flotation performance. Experimental pilot plant results were obtained from a continuous circuit consisting of a ball mill and a rougher bank of four flotation cells. A series of flotation experiments were run with different conditions of dissolved oxygen keeping constant most of other operating conditions. Variables difficult to control were measured or estimated. A set of mathematical models suggested in the published literature were tested. It was found that the copper recovery results were highly sensitive to the concentration of dissolved oxygen. The latter confirmed the reactions mechanisms associated to the collecting capabilities of current flotation reagents used in the flotation of copper sulfide ores. On the contrary, the results of molybdenum recovery varied non-linearly with changes in dissolved oxygen which could not be explained directly in terms of the computed specific rate constants. The error analysis confirmed that Klimpel continuous model exhibited the best specific rate constant correlation with molybdenum recovery data. Lastly, the classic particle force balance equations indicated that the adjustments of estimated mineral liberation with experimental data exhibited a turning point when passing from highly liberated to low liberated particles. The values of recovery for such turning point correlated well with the concentration of dissolved oxygen. The results indicated corrections due to changes in the dissolved oxygen must be incorporated in future froth flotation descriptions.

Jan 1, 2014