Search Documents

Building block cave operations involves significant commitments to capital and time before arriving at completion or full scale production. The focus of the mine planners is habitually to minimize the upfront capital and achieve target development rates and construction dates. However, management of the construction of block cave projects requires understanding of the strategic/tactical decisions that the project planners can take to increase the value of the project and cope with schedule risks. As block cave projects involve development of vertical and lateral access to the underground footprint and ore body, more attention to uncertainties is required as most of those projects need to confront delays leading to a decrease in value. The standard process of development and construction scheduling is often not planned with provisions to realise uncertainties in the construction phase. Uncertainties and risks such as poor ground conditions and equipment breakdowns are the primary reasons for schedule slippage. As a result most of these projects experience delays. In order to mitigate risk to the project schedule, we centre our attention on increasing the value of the project by spending more money upfront and improving the schedule rather than focus on minimizing the upfront capital. Scenario-based simulations along with project appraisal methods hold great promise to address those challenges and enable mine planners to make better decision and enhancing project value. This paper will present a hybrid simulation based approach to model the process of development and construction based on discrete event (DES) and continuous based simulations (CBS). Such simulations will be conducted for scheduling multiple development headings (Apex, undercut, extraction, and haulage levels). Two scenarios are utilized to investigate the impact of risks on schedule and how flexibility in construction can be incorporated to minimize delays: (1) The bench mark scenario is created based on actual data from an existing case study. (2) The flexible scenario is an alternative way to use primary production infrastructure (Ore Handling System) including conveyer drifts, conveyer belt, and a gyratory crusher in the construction process and that requires the system to be in place before starting lateral development in the footprint. As a result, a truck-load system versus a conveying system will be studied to observe the gain in construction speed and to estimate the premium cost for mucking rock from headings and drawbells and transporting it to surface. Such strategic flexibility would be part of the original budget to make the tactical flexibility available during construction. Continuous cash flow models are used to evaluate and compare the scenarios. This study is an example of how new techniques in planning can be used to support decision-making processes in block cave design.

Aug 1, 2013

By Lei Wang, Xiancong Zhao, Hongxu Li, Lifeng Zhang

Recovery of Fe from red mud has been studied using statistically designed experiments. The effects of three factors, namely: reduction temperature, reduction time and proportion of additive on recovery of Fe have been investigated. Experiments have been carried out using orthogonal central composite design and factorial design methods. A model has been obtained through variance analysis at 92.5% confidence level.

Jan 1, 2015

By Autumn Fjeld

The filling and solidification of an industrial scale composite casting has been simulated to investigate the re-melting of the outer shell material during casting. Production of a composite requires a two-step casting process which includes centrifugal casting of a hard outer shell material and subsequent conventional ingot casting of a tough core material into the shell-mold assembly. During casting of the core, a thin layer of the shell material re-melts and mixes with the core material producing a bonding layer of intermediate composition. The present numerical model employs the volume of fluid method and an enthalpy-porosity technique to couple the filling of the core material and re-melting of the shell material. The interface between the two materials is tracked and can be used as a guide to examine the extent of melting and, to some degree, mixing of the two fluids and species diffusion.

Jan 1, 2008

By V. Buddhiraju

Liberation of valuable minerals from low-grade ores is an important challenge faced by the minerals industry. It is necessary to grind the ore to extremely fine size in order to liberate the mineral from such ores. As it is generally difficult to produce ultra-fine particles of diameter less than 1-10 mm using traditional grinding equipment such as the ball mill and the rod mill, several new types of grinding equipment are being explored to produce ultra-fine particles. The planetary ball mill is promising in that it makes grinding to submicron sizes possible by imparting high energy to the ground powder. It consists of four pots mounted on a disc. The pots and the disc rotate in opposite directions. The high-speed rotation of the pots and the disc makes the grinding media move at high centrifugal velocities leading to high impact energy and thus improving the grinding performance. The planetary mill is currently being explored for its capability to produce ultra-fine particles on a large scale. In this context, there is a need to understand the dynamics of the ultra-fine grinding process within the planetary ball mill. Here, we address this need by studying the movement of media inside the pots and the grinding kinetics, employing both experimental and modeling techniques. The experiments were conducted on a laboratory scale planetary ball mill. The particle size distribution of the powder was monitored to allow an estimate of the grinding kinetics to be made while video images of the media profiles inside the pots were taken to allow quantification of the media movement as a function of process variables such as the media size, percent loading, etc. The population balance modeling (PBM) framework was used to simulate the kinetics of grinding while the discrete element modeling (DEM) approach was utilized to simulate the media motion. The parameters in the population balance model, namely, the specific breakage rate and the breakage distribution function were determined by fitting the batch grinding equation to the experimental data. The predicted powder size distributions were found to be in good agreement with the experimental data. The experimental results have shown that the mean size of product powder decreases as the size of the grinding media decreases. Keywords: population balance modeling, grinding, planetary mill, grinding media, grinding equation

Sep 1, 2012

By A. Oleynik

The paper presents the modeling and simulation experience in field of research into the mineral row material separation. Features of the subject domain generate a need of different approaches using for modeling and simulation. The general scheme of separation modeling process and some examples of developed simulators are presented. Due to different features and standards of knowledge, different formal models should be used to describe and analyze separation processes in the frame of general scheme. The main types of simulators are modeling tools based on statistical, analytical or expert (logical) methods. Advantages and disadvantages of each method above are mentioned in the paper. Based on the studies in the field considered, an aggregated mathematical model of a separation process has been generated. The aggregated model is presented by two combined equations. The actual apparatus model is formed by substitution of the particular functional expressions of the forces acting in the apparatus into the aggregated model. The developed simulation tools are based on integration of limited experimental data and on ex-pert knowledge available on separation processes. Simulation does not allow producing an optimal, in terms of mathematics, solution. But it gives an opportunity to produce and compare quite adequate solution variants in real time. The results obtained are the theoretical and methodological base for integrated information technologies to support studies in the field of mineral processing.

Jan 1, 2014

By D. Minson, G. Lindner, S. Miskovic

Thickening is an essential process in the mining industry and currently one of the most costeffective technologies for recovering water from tailings streams. Improved performance of thickeners can lead to significant energy savings, and increased recovery of water close to the mill. Its performance is directly affected by the rheological properties of slurry at the raking zone. To further understand the dynamics of slurry inside thickeners, a multiphase particle-in-cell method (MP-PIC), is used to investigate its behaviour on a sliding channel. Parameters like particle-size distribution (PSD), rheological properties, and domain discretisation were benchmarked to study its influence on highly dense flows. The slurry concentration by weight of the samples ranged between 60 to 70% of solids by mass, under which conditions it is known to behave like a non-newtonian fluid. An experimental flume, 0.2 m wide and 2 m long, was used for validation under different angles and slurry composition. Visual analysis of the slurry height and shape, as well as velocity measurements, were compared with numerical data. The preliminary results show the ability of MP-PIC method to represent different types of highly dense flows, while still being computationally affordable and practical for industrial use. Keywords: High density thickener, CFD, MP-PIC, multiphase-flow

Jan 1, 2020

By Yinlun Huang, Cristina Piluso

"Industrial sustainability looks to improve material and energy efficiencies, product quality and variety, and productivity within an industrial zone, thus pursuing the long-term sustainable development of the region. This paper introduces a systems-based modeling methodology, which provides industrial leaders with the ability to evaluate future production options to determine the path with the highest positive impact on both their own sustainability and the sustainability of the overall industrial zone. A comprehensive study on sustainable development of an auto-manufacturing focused industrial zone under uncertainty is illustrated. Introduction Sustainable development refers to a continuous process of improvements that must be followed in order to achieve a state of sustainability. In its most general form, a widely accepted definition of sustainable development was identified in the Brundtland Report as “development that meets the needs of the present without compromising the ability of future generations to meet their own needs” [1]. Sustainability will occur when the material and social conditions can be maintained or improved for human health and the environment over time without exceeding the ecological capabilities that support them [2]. An industrial zone is defined as a geographic area comprised of a network of industrial sectors, each composed of a number of entities. Such a zone is highly integrated and complex, and contains numerous uncertainties due to each industry’s extreme dependencies on its suppliers and customers throughout the product supply chain. In the context of an industrial zone, sustainable development refers to the continuous process of advancements in the zone, and the entities that comprise it, creating more value, wealth, and profits (economically viable dimension), providing environmentally benign products with less raw material consumption and waste generation (environmentally compatible dimension), and collaboration and a direct relationship among all involved in the supply chain, i.e. the corporations, employees, suppliers, investors, communities, customers, etc. (socially responsible dimension). These are the triple bottom lines of industrial sustainability. In order for economic, environmental, and social sustainability improvements to be successfully implemented within an industrial zone, the industrial decision makers must possess system-wide analysis abilities. The analysis of the sustainable development of an industrial zone involves decisions at the local, regional, or national levels, which are complex and involves handling of ill-defined parameters with a high degree of uncertainty due to the imperfect understanding of the underlying issues [3]. This work focuses strictly on the uncertainties in future zone planning and a method and a scenario based approach that can be used to evaluate the sustainable development of an industrial zone among these uncertainties."

Jan 1, 2008

By T. Doshi, S. Ennis

The interactions of mineralogy, chemistry and mechanics in an ore concentrator are highly complex, so it can be particularly hard to measure the impact of a specific plant change. This paper describes a method for modeling a baseline metal recovery so that the value of a major process improvement can be measured. This method only considers variables that cannot be optimized, such as mineralogy or equipment availability. Linear regression, random forests, extreme gradient boosting (XGB), and simple two-layer neural network models of process value were tested as candidate baseline model types to assess the impact on throughput and recovery in a copper concentrator. XGB performed the best in terms of accuracy for both throughput and recovery models. Using an XGB model of nonoptimizable process variables provides an accurate, dynamic baseline that can be compared with continuous production data to determine the value added by major concentrator improvements in real time.

Mar 2, 2022

By Cristian Castro, Patricio Gomez, Jason Furtney, Augusto Lucarelli

Understanding the effect of blast-induced damage on wall stability is critical to safe and economical operation of a large open pit mine. In this paper, a numerical modeling methodology is used to investigate damage to rock at distances of greater than 10 m (30 ft) from a blast. Two applications of the model are demonstrated. First, the method is applied at a zinc mine to investigate the effects of pit blasting on underground stopes. The model predictions are compared with the scaled distance formula. Next, the model is used to investigate a common problem in which unfavorable jointing creates a rock wedge on the bench above a blast. The numerical predictions of wedge stability are compared to the wedge energy model described by Wong & Pang (1992).

By Cristian Castro, Patricio Gomez, Jason Furtney, Augusto Lucarelli

Understanding the effect of blast-induced damage on wall stability is critical to safe and economical operation of a large open pit mine. In this paper, a numerical modeling methodology is used to investigate damage to rock at distances of greater than 10 m (30 ft) from a blast. Two applications of the model are demonstrated. First, the method is applied at a zinc mine to investigate the effects of pit blasting on underground stopes. The model predictions are compared with the scaled distance formula. Next, the model is used to investigate a common problem in which unfavorable jointing creates a rock wedge on the bench above a blast. The numerical predictions of wedge stability are compared to the wedge energy model described by Wong & Pang (1992).

Feb 1, 2020

By A. Jolly, A. R. Kumar, T. Novak

"Coal dust in an underground mining environment has always been a health and safety issue, especially with longwall faces. Efforts to integrate a flooded bed scrubber within a longwall shearer are underway at the University of Kentucky (UK) in order to reduce worker exposure. Small scale physical modeling, along with numerical modeling and full-scale experimentation, are engineering tools being used by the authors for simulating operating conditions underground. Recent advances in affordable 3D printers offer a means for fabricating accurate, scale models for obtaining timely insights prior to building large, expensive test setups. Comparing scale-model results with those of numerical modeling, such as computational fluid dynamics (CFD), lends guidance and confidence in a proposed design. This paper discusses the small scale physical modeling of a longwall shearer with an incorporated flooded-bed dust scrubber in its operating environment. A tracer gas was used in the model in lieu of dust to measure dust-capture efficiency. Results obtained from the physical model are presented and compared with those of CFD model. The small scale model and the CFD results demonstrate good agreement. Dust capture results show the usefulness of the approach in advance of actual design implementation. INTRODUCTION Dust generation is an undesirable, but innate, byproduct of all mining processes. It is particularly problematic in underground coal mining, where coal dust poses potential hazards from two fronts – health and safety. Coal workers’ pneumoconiosis (CWP), commonly referred to as black lung, is a chronic lung disease that results from the long-term overexposure to respirable coal dust. Although the prevalence of CWP has steadily and significantly declined over the three decades following the Federal Coal Mine Health and Safety Act of 1969, a study by the National Institute for Occupational Safety and Health (NIOSH, 2008) indicates that significant health hazards associated with respirable dust still exist within the coal mining industry. In addition to being a health hazard, coal dust can also create a potentially catastrophic safety hazard. Dust that is not captured or removed at its source is suspended and transported throughout a mine by the ventilation airstream. This float dust is subsequently deposited on all mine surfaces along its travel path. Without adequate application of rock dust to these surfaces, forces from a localized methane ignition can reintroduce settled coal dust into the airstream creating an explosive mixture. Therefore, capturing and removing coal dust at its source mitigates the probability of two major hazards that continue to plague the mining industry."

Jan 1, 2016

By Nader Golchinfar

The V-shaped breakout zones around excavation boundary have been captured by a variety of numerical methods such as FEM and DEM methods and failure criteria such as linear Mohr-Coulomb and non-linear Hoek-Brown failure criteria. Diederichs (2007) used Hoek-Brown parameters in his spalling parameter approach; however, it is convenient to use the Mohr-Coulomb failure criterion in numerical modeling because the material parameters can easily be determined from triaxial tests. It is shown in this study that when peak and residual strength parameters as well as dilation angle are properly defined, the brittle Mohr-Coulomb failure criterion can be used to model depth of failure near underground excavations satisfactorily.

Jun 1, 2012

By B. R. Rakishev

When transferring from the railway transport to the combined automobile-conveyor transport by stripping works at coal mines, it was proposed to form at the working zone a concentration horizon for reloading part of the stripping volumes at the railway transport. The procedure is adduced for substantiating parameters of bloc-panels, depending on the step of transferring of reloading points in interrelation with the intensity of the working zone?s development with observing the balanced movement of parts of the stripping zone relative to the temporary concentration level.

Aug 1, 2013

By Tracie Zoeller, G. Paul Neitzel, Jr. Sanders

"As-cast aluminum alloy microstructures are highly segregated. Coring occurs when solute is rejected during solidification resulting in a non-uniform solute distribution across the dendrite arms. Furthermore, coarse constituent phases are formed in the interdendritic channels as a result of coring. These constituents are either soluble or insoluble and can form from the major solute elements, or the impurity elements iron and silicon. The dissolution of these constituents, or precipitates, is an important consideration in determining the microstructural evolution of an alloy. The properties of an alloy are significantly affected by the morphology of the precipitates. Precipitates can be either beneficial or detrimental to the properties of the alloy. Regardless of the nature of the precipitates present, it is useful to understand the dissolution process. A diffusion-limited dissolution model for planar precipitates is presented and compared to an analytical solution and to models in the literature. Experimental data for the dissolution of AhCu precipitates in a dilute alloy are compared to the proposed model, and found to have good agreement.IntroductionDuring the processing of precipitation hardening alloys, beginning with the solidification step and ending with the final solution heat treatment step, there are a variety of particles that form. Some of these particles must be dissolved and solute redistributed and re-precipitated in a more homogeneous distribution than as it first appears in the as-cast microstructure. The particles can be classified by size as microstructure develops during processing. Coarse particles that form during solidification can either be soluble, partially soluble or insoluble. These particles can be the result of coring of the major solute elements as well as the limited solubility of elements such as iron and silicon. Iron and silicon can combine with the major solute elements such as copper in 2xxx alloys forming coarse phases often referred to as constituents. Figure l(a) shows the as cast microstructure with the constituent phases, Al2Cu and the insoluble phases containing Al-Fe-Mn-Cu-Si [7]. During preheating the soluble phases such as AhCu that were formed by eutectic decomposition must be taken back into solution, whereas the insoluble phases remain essentially unchanged by most thermal operations. Figure l(b) and (c) show some of these constituents after subsequent preheating steps. In Figure l(b), AhCu platelets are shown and in Figure l(c) the remnants of the AhCu phase are noticeable, along with Mn-bearing dispersoid phases. Reduction in the volume fraction and size of these phases improves the ductility and fracture toughness as well as reducing the inhomogeneities that occur during coring. Also adjacent to the dendrite arm boundaries are smaller AhCu phases that precipitate by a solid-state transformation on cooling to room temperature."

Jan 1, 2004

By Xiaozhen Song, Ashish Ranjan Kumar, Shimin Liu

Diesel particulate matter (DPM) presents considerable occupational health hazards in large-opening underground mines (LOUM). Previous studies have primarily focused on single-factor ventilation optimization, while this study innovatively integrates booster fan configurations and loading truck placements utilizing a coupled CFD turbulence– particle tracing model. Three distinct ventilation scenarios were simulated to evaluate their effects on DPM dispersion and dilution efficiency. Results indicate that dual booster fans substantially enhance airflow distribution, increasing air exchange efficiency; The migration pathways of DPM closely align with airflow trajectories, and air turbulence is the predominant factor influencing particle behavior within LOUM; Loading truck placement critically influences local DPM accumulation; Airflow effectiveness diminishes markedly with distance from booster fans, causing stagnation at distant working faces. The coupled CFD approach accurately captured airflow–particle interactions, offering valuable insights for optimizing mine ventilation strategies and effectively mitigating miners’ exposure to DPM.

Feb 22, 2026

By R. Panchadhara, P. A. Gordon

"A non-local formulation of the classical continuum mechanics theory called peridynamics is used to study driven dynamic fractures in the near wellbore region. The controlled application of pressure on the inner walls of the borehole can create multiple fractures that initiate from the borehole and propagate radially outwards. The number of fractures and their final length is primarily controlled by the loading conditions and the in-situ stress state. It is the objective of this study to use a novel numerical technique called peridynamics to investigate: 1) the influence of loading conditions on the fracture pattern, and 2) the effect of in-situ stress anisotropy on fracture propagation. In peridynamics, the momentum equation from the classical theory of solid mechanics is replaced with the integro-differential momentum equation of peridynamics. As a result, discontinuities in the body, such as fractures, occur spontaneously. The material is represented by a continuum body with material points that interact non-locally with all other points within a specified distance. Fractures evolve using a damage evolution law, wherein, damage at any point is represented as the ratio of the number of broken bonds in the current configuration to the total number of bonds in the undeformed configuration.Simulations are performed on simplified test cases and the results are compared to relevant experimental and numerical studies found in the literature. It is shown that the loading rate significantly influences the number of fractures initiated from the borehole as well as their extent. Low loading rate produces fewer fractures, whereas, high loading rates produce longer and a greater number of fractures from the borehole. The numerical method successfully predicts the fracture pattern obtained from a wide range of loading conditions. Predictions also indicate that the resultant fractures are shorter when in-situ stress is applied. Numerical results show that the fracture extension is largest in the direction perpendicular to the minimum compressive stress for the range of loading rates used in this study."

Jan 1, 2015

By Trevor Saldanha, Stephen Grano

A key uncertainty in flotation modelling is determining, a priori, the proportion of each mineral which will develop a certain rate constant for different feeds of variable liberation and hydrophobicity. First and foremost is the ability to predict the proportion of minerals that will not be recovered, even at very long flotation times, and which therefore report to the tailing stream as the non-floatable fraction (i.e., having a zero rate constant). Recent empirical evidence has shown there are at least two factors which contribute to heterogeneity in flotation rates and, in particular, the appearance of a non-floatable fraction - heterogeneity in the distribution of adsorbed collector amount on a particle-by-particle basis, even for the same liberated mineral in the same size fraction, and the existence of a critical contact angle below which the bubble-particle union is unstable. In this paper, inferred measures of the critical contact angle for liberated particles are compared with stability model calculations. This has demonstrated that the dependency of the critical contact angle with particle size and mineral specific gravity are explicable by bubble-particle stability considerations. This argument is further extended to the flotation behaviour of composite particles, as a function of liberation class, in which the wetting perimeter of the attached bubble on the composite particle surface is constrained by the presence of a surrounding hydrophilic domain for the case of simple binary composites. Reasons for the different flotation behaviour of composite particles having different textures, as noted in the literature, are also discussed in terms of the stability of the bubble-particle union.

Jan 1, 2014

By L. Barry Phelps

Selective placement of toxic mine spoil has been advocated as a means of reducing environmental pollution at the source. However, methods of burial to control environmental damage are not well understood. Using the finite element method, this paper investigates the burial of toxic material within low permeability packages below the water table. The effects of differing package and spoil permeabilities, anisotropy, and infiltration are modeled to provide information about their effect on the flow system. Projection can then be made which will assist in the design of the toxic spoil packages to limit the amount of groundwater flowing through them.

Jan 1, 1983

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