Search Documents

By Alejandro Gutiérrez, David Ahues, Paola Merino, Felipe González

A Discrete Elements Method (DEM) study of load movement within a SAG mill 12,192 × 7620 (40′ × 25′) in normal operating conditions is performed. The purpose is to determine how the geometry of the lifters affects the kinetics of the load inside the mill. Specifically, the changes produced by the incorporation of inclined lifters and spiral or helicoidal pulp lifters to the kinetics are evaluated. Four study cases are identified: first, the base case corresponds to a conventional geometry; the second case corresponds to inclined lifters incorporation; third case corresponds to a mill with helicoidal discharge; and the fourth case corresponds to a mill with inclined lifters and helicoidal pulp lifters. To save computational resources, the water that enters the mill has been modeled by varying the coefficient of restitution. The results of a series of experiments are presented to determine how this coefficient varies when having impacts on surfaces with different liquid levels. A comparative analysis is made of the following variables: mean transport speed and material flow, for normal operating conditions. Taking as reference the mill with straight lifters and radial pulp lifters (case 1), the case with the inclined lifters and radial pulp lifters increases the material discharge flow less to 76 mm of size by 5.1%, and this value increases even more when the helicoidal pulp lifters are incorporated reaching a 58.6%. The most significant change is seen in the flow through the grates and the pulp lifters discharge, which can increase substantially improving the operation.

By J. C. Kurnia

Methane is recognized as one of the prime causes of underground mine disasters around the globe. A cost-effective mine ventilation system and methane control is required to ensure a safe and productive environment in an underground mine whilst keeping energy usage and operating cost at minimum. Typically methane enters a mine environment or tunnel from localized distributed sources at high initial concentration. Hence, an understanding of its dispersion from discrete sources is important to design a methane management system via effective ventilation. This study employs computational fluid dynamics to model methane dispersion in an underground tunnel which has a number of discrete sources of methane and effective ways to dilute and or disperse it. Air flow distribution within the tunnel as well as methane dispersion is simulated by solving the governing transport equations subject to appropriate boundary conditions. First, the effectiveness of ventilation duct placement is examined. Several possibilities are explored and investigated. It is found that the common configuration utilized in underground mine performs best as compared to other configurations in maintaining the concentration of methane below permissible levels. Next, a novel intermittent airflow ventilation system is proposed and evaluated via simulation with the goal of reducing the energy cost whilst maintaining methane level in the mining face below the allowable level. Parametric studies are conducted to investigate effects of various factors influencing the effectiveness and performance of this novel intermittent ventilation system. This study will provide some new ideas for designing an efficient and effective underground mine ventilation system by lowering the net average amount of air that needs to be circulated in the tunnel for a safe operation. KEYWORDS

Aug 1, 2013

By J. D. Zhu

A numerical model for predicting microporosity formation in aluminum castings has been developed, which describes the redistribution of hydrogen between solid and liquid phases, the transport of hydrogen in liquid by diffusion and Darcy flow in the mushy zone, One of the key features of the model is that it uses a two-stage approach for porosity prediction, In the first stage, the volume fraction of porosity is calculated based on the reduced pressure, whereas, in the second stage, at fractions solid greater than the liquid encapsulation point, the fraction porosity is calculated based on the volume of liquid trapped within the continuous solid network, which is estimated using a correlation based on the Niyama parameter, The porosity model uses thermal data calculated using the commercial finite element package ABAQUSTM. The model has been applied to examine the evolution of porosity in a series of experimental samples cast using unmodified A356 in which the initial hydrogen content was varied from 0,048 to 0.137 [cc/100g], A comparison between the model predictions and the experimental measurements indicates good agreement in terms of the variation in porosity with distance from the chill and the variation resulting from initial hydrogen content.

Jan 1, 2005

By Jie Ma, Shun-Li Zhao, Jian Zhao, Jie-Yu Zhang

"Microstructural morphology evolution of Ni-45wt.%Mo droplets during rapid solidification process was investigated using an advanced modified 3D cellular automaton-finite element(CAFE) model by CALCOSOFT software. The model was modified by User-defined functions (UDFs) and validated with experimental solidification microstructural morphology of Ni-45wt.%Mo droplets available in the literature. Then, the influences of convection heat transfer coefficient, superheat of the alloy droplets and droplets size were investigated during the rapid solidification process. The predictions show that microstructural morphology is obviously sensitive to convection heat transfer coefficient and droplet size. With the decrease of convection heat transfer coefficient, the grain density decreases. In addition, the decrease of the droplet size could obtain refined microstructure. However, decrease the superheat ofNi-45wt.%Mo droplets has little effect on microstructural morphology.IntroductionFor several decades, much attention has been paid to undercooled alloy droplets of rapid solidification process for which can obtain refined microstructure, decreased microsegregation and extended solubility of alloying elements [1-3]. Extensive prior work shows that solidification of droplets plays an essential role in advanced material properties [ 4]. Consequently, it is very important to investigate the microstructural morphology evolution, especially for dendrite growth during non-equilibrium solidification process, for dendrite growth reflected the kinetics of the crystal growth process."

Jan 1, 2013

By Hebi Yin

A coupled cellular automaton (CA) - finite element (FE) model was developed to calculate the growth of dendrites during the solidification of cubic and hexagonal metals. The model solves the conservation equations of mass, energy and solutes in order to calculate the temperature field, solute concentration and the growth morphology of dendrites, including the grain structure and the dendritic microstructure. Validation of the model was performed by comparing the simulation results with calculation results from previously published works, showing qualitatively good agreement in the dendritic morphology. Application to magnesium alloy AZ91 illustrates the difficulty of modeling dendrite growth in hexagonal systems, observed as deviations in growth direction caused by mesh-induced anisotropy. The model was applied to the simulation of small specimens with single- and multiple- equiaxed grain growths and columnar grain growth in directional solidification. The influence of cooling rate and some kinetics parameters on the grain morphology are also discussed.

Jan 1, 2009

By Trutwin W, Tracz J

The outbreak of open fires underground cause disturbances to the conditions of the mine ventilation process. Thermal and chemical effects due to fire create unsteady, transient states in airflow and the spread of combustion products in the mine ventilation network. Changes in airflow have an influence on the fire, constituting a sort of feed-back between the fire and the ventilation process. This paper deals with a numerical simulation method, which enables the determination of transients caused by fires. The mathematical model, of the simulation method considered, has been based upon a: simple model of combustion of fuel in the airway, unsteady exchange of heat between the fire, airflow and surrounding rocks, airflow in mine ventilation networks and the theory of heatflow in the rocks. Initial and boundary conditions allow the determination of a set of parameters describing the airflow, the propagation of combustion products, the tempera- ture of the fire and flow etc. due to the fire. The simulation method, facilitates a graphical presentation of the spread of combustion products throughout the network. The simulation method has been prepared for an IBM PC/AT Computer.

Jan 1, 1988

By Husheng Li, Zhongli Wu

This paper discusses the fundamental strategies adopted in the simulation package FIRES developed for ventilation simulation under mine fires. These fundamental strategies include (1) fire characteristic curves specifying the development of a mine fire with time, (2) simulation of fume spread in a ventilation system, (3) determination of the aerodynamic effects of mine fires on flows; and (4) dynamic simulation of flow state under mine fires. The major functions of the simulation package are also addressed. The package has been used to simulate an experiment of mine fire. A comparison is made between the simulation results and the measured data.

Jan 1, 1993

By P. S. Heyns, M. Heyns

A computer model has been developed of mine hoisting systems in deep South African gold mines. Particular attention was paid to the contribution of different modelling assumptions for the guides and guide-roller assemblies. Objective functions indicative of guide-roller tyre endurance, slipper-plate forces and perceived vibration levels were introduced to quantify dynamic performance. These objective functions were used to compare models. It was found that the choice of guide and guide-roller models did not influence the forces on the guide rollers significantly for typical shaft and conveyance characteristics at typical levels of viscous damping coefficients. The rigid and finite-element guide models gave the same results as those of the effective-stiffness guide model, and while the latter is suitable for simulation for design purposes, the finite-element model should be used for unusual guide geometries

Jun 19, 1905

By J. A. C. Diering

The simulation of mining in faulted ground was recently described by Crouch. His method is extended here to enable a wide range of non-homogeneous mining problems to be modelled. Essentially, the boundary conditions at geological interfaces are incorporated into the overall iterative solution, which is necessary if faults and joints are to be modelled realistically. Two examples are given to demonstrate the method: the mining of coal a seam below a thick dolerite sill, and the mining of a flat reef in close proximity to a vertical fault where hangingwall and footwall rocks have different properties.

Jan 1, 1980

By I. Ohnaka

Although many works have been reported on the semisolid casting, the mold filling behavior of semi-solid alloys has not been so well understood. This paper compares the physical and numerical simulation results on the mold filling behavior. Physical simulation was performed by using solid-water mixture. The numerical simulation was made under the assumption of one phase flow with an effective viscosity. It was found that the solid particles stabilize the flow greatly if the kaction solid exceeds about 0.4, because viscosity increases greatly. The overall flow behavior can be estimated by the numerical simulation. Interesting flow behavior was also found which can not be observed in the case of one phase flow and which is difficult to simulate by the one phase flow assumption.

Jan 1, 1997

By G. R. Marchant, B. K. Loveday

A mathematical model has been developed for flotation circuit simulation. This model was designed to investigate the effect of varying the circulating loads and to determine whether an improvement in overall efficiency could be gained by a size separation prior to flotation. Model parameters were obtained by regression to plant data. Highly correlated parameters in the model were combined to reduce the number of parameters so that the need for excessive data gathering was eliminated. Mass balance smoothing was used to obtain a consistent data set for regression. The results of the simulation indicate that for the same grade of concentrate, an increase in copper recovery of 0.6 per cent may be possible by increasing the circulating loads. The results of the simulation also indicated that recleaning of the plus 200-mesh material may be beneficial.

By G. R. Marchant, B. K. Loveday

"A mathematical model has been developed for flotation circuit simulation. This model was designed to investigate the effect of varying the circulating loads and to determine whether an improvement in overall efficiency could be gained by a size separation prior to flotation.Model parameters were obtained by regression to plant data. Highly correlated parameters in the model were combined to reduce the number of parameters so that the need for excessive data gathering was eliminated. Mass balance smoothing was used to obtain a consistent data set for regression.The results of the simulation indicate that for the same grade of concentrate, an increase in copper recovery of 0.6 per cent may be Possible by increasing the circulating loads. The results of the simulation also indicated that recleaning of the plus 200-mesh material may be beneficiaL."

Jan 1, 2014

By Ivan V. Klumpar

If components of a gas-liquid stream react with each other and with stationary solids, involving parallel and sequential, reversible and irreversible, and slow and instantaneous reactions, the process can be defined by a set of simultaneous differential rate equations and simulated as multistage semibatch re¬actors. Each stage is divided into finite time increments. Within each increment, even the simultaneous reactions are treated as sequential. The proposed method follows process evolution from non-steady state startup to normal operating conditions. It allows computing byproduct buildup in recycles, tracer response curves and, in case of underground processes, expansion to adjacent parts of the deposit.

Jan 1, 1979

By Shreyas Giridharan, Fursan Hamad, Sujith Gowda, Christian Moormann

"Pile installation in soil is a dynamic process which exhibits large deformation in the soil region. Experimental investigation of pile installation is both time and resource intensive. Numerical investigation can be a helpful tool for the analysis of the installation process. Furthermore, numerical analysis allows simulating various soils and loading conditions with ease. Finite Element Method (FEM) – a widely used framework for numerical simulations - fails to perform well for soil mechanics where large deformation and heavy material movement exists. In order to alleviate this shortcoming, Material Point Method (MPM), advancement to FEM has been used. In MPM, material points (particles) are used to represent the domain which carry the information and move over a fixed computational mesh where the governing equations are solved.In the present work, a variation of MPM, called CPDI (Convected Particle Domain Interpolation) is utilised which handles large deformation problems better than the classical MPM. Here, the material points are represented by a spatial domain which deforms according to the particle deformation. This eliminates stress oscillations while particles cross the grids. Implementation is validated and applied to pile installation in dry and saturated soil conditions. A 2D axisymmetric model is chosen, with penalty function method being used for evaluating contact between soil and pile. Sand is modelled via the hypoplastic constitutive description that takes into account large deformation, as well as small strain stiffness. From these simulations, various parameters like stress, void ratio and friction angle are studied.INTRODUCTIONOffshore piling is challenging in terms of testing and execution. Before the execution, predicting an efficient installation method requires study of various parameters, like soil conditions, type of loading, force requirement etc., and are the minimum parameters to be known beforehand. To obtain these parameters, testing a real scale prototype under different conditions is not feasible. Scaled model tests in lab can support in making certain decisions, but fail to replicate real behavior fully. Therefore, the need for a numerical simulation technique emerges, which can provide early predictions easily and accurately. This later can be extended to real scale tests, reducing the number of experimental investigation required.Simulation of displacement pile installation in sand involves large deformation due to penetration process, thereby making the application of FEM impossible, where large deformation causes heavy mesh distortions. Due to this, meshless methods are developed to treat large deformation in a better way. Among various meshless methods, Material Point Method (MPM) shows more applicability to geomechanical problems. This is mainly due to the ease of handling soil, structure and water in a unified manner."

Jan 1, 2017

By S. Kenzap, V. Kazakidis

Risk assessment is an essential component of mine feasibility studies. However, it is often challenging to identify the appropriate analysis and derive the required input data. The uniqueness of a mining project and the associated technical issues require customization of the risk assessment approach to meet both requirements and characteristics. The uncertainty driving the operating risk in mine feasibility studies and its assessment through Monte Carlo simulation as a component of a cash flow analysis are discussed in this paper. The paper introduces the use of the analytic hierarchy process for derivation of qualitative input and conducts quantification of the three main sources of operating performance delays, equipment, ground and schedule problems, in a cash flow simulation. A methodology for the assessment of operating risk in mine feasibility studies is presented.

Jan 1, 2011

By Luiz R. P. De Andrade Lima

This paper presents a mathematical model for the oxide copper ore heap leaching process. The data used in the model include physical-chemical, .geometrical and operational data, as: leachable metals ore contents, .flow rate and acid concentration, parameters of passivity, ore size distribution, the average residence time of the solution in the heap, height, irrigated area and weight of ore in the heap. In this model the heap is divided into horizontal and plane layers of constant area. The flow of the liquid in the heap is considered unidirectional and the solid-fluid reaction model is considered under diffusive control. The average residence time of the solution into the heap and the apparent effective diffusivity of the acid through the ore particles have .a significant influence on the rate of the copper extraction and its concentration in the pregnant solution; therefore these parameters are used in the model calibration. When this algorithm was applied to a laboratory test, the results showed that the model is adequate to predict the process performance.

Jan 1, 1999

By Fedotov Pavel

"Today the main objective of improving the process of ore dressing is to reduce energy consumption and performance of the basic rule of auxiliary enrichment processes. One of the modern means of ore destruction with a minimum of energy consumption is the interparticle fracturing, for example in the roller press. There still remain the unsolved problems of the energy consumption determination and the particle size distribution of the crushing product .Existing methods for determining these parameters are empirical and require considerable expenses and time. The improved finite element was developed to solve these problems using finite element contact interaction of the rock. Also it was developed to clarify the stress-strain and the rock fracture energy in a layer. These elements are used to study ore destruction in the layer and allows to create a model taking into account the interaction of crushing particles surfaces. The same algorithm was developed determining crushing product particle size distribution in the layer on the base of the Bond’s theory of crushing. The presented finite element model allows to determine the quantity of strain energy and particle size distribution of crushed product during the destruction processes of the material in the roller press with fair accuracy."

Jan 1, 2014

By J. R. Percival, L. Ip, M. D. Piggott, S. J. Neethling

Component wear due to particle-laden flow is a significant challenge in the minerals processing industry in the handling and transportation of abrasive mineral slurries. Studying wear through experimentation is time-consuming and expensive; therefore possessing the ability to accurately simulate, understand and predict wear in these complex environments is highly desirable. The ability to predict wear is important in both the design of equipment in order to maximise the performance over the wear life of the components, as well in minimising potentially costly and dangerous failures and unscheduled maintenance.. Particle-induced wear is a complex problem, requring an understanding of coupled fluid-solid motion, geometry deformation in response to fluid-solid behaviour, moving geometries and wear modelling. Although the literature demonstrates that rebounding particles and changes to surface profiles can have a significant effect on subsequent wear rates and patterns, the challenge of implementing a model to account for this effect in complex geometries has been addressed by only a few authors in the literature. This paper describes a new wear simulator, leveraging fluidity, an open-source, massively parallel, finite element based CFD simulator capable of unstructured adaptive meshing. The wear simulator includes the ability to model particle-laden flow using representative Lagrangian particles and dynamic remeshing of geometry in response to wear and component movement. These capabilities are demonstrated using wear and the associated dynamically deforming boundaries within a Coriolis tester arm, with the results comparing favourably with the experimentally observed wear patterns. Keywords: CFD, wear, two-phase flows, anisotropic mesh adaptivity, boundary deformation

Jan 1, 2020

By A. Henchiri, G. A. Gruber

INTRODUCTION Phosphatic resources may be quantified in terms of ore tons and phosphate analyses; however, to quantify phosphate reserves it is necessary to estimate the tons and quality of phosphate rock that can be extracted from the ore. The yield and grade of phosphate rock that can recovered from sedimentary phosphate deposits are usually difficult to estimate because of variations in stratigraphy and/or mineralogy. Prospecting phosphorite ores for purposes of estimating reserves and preparing mine plans requires considerably more effort than defining ore thickness and P205 content. In Florida, for example, it is standard practice to prospect by core drilling and to subject each strata (split) to beneficiation testing, and to perform detailed chemical analyses on the beneficiated products from each strata. For a completed drilling pattern the area of influence of each hole is 2.5 acres. The number of strata per drill hole typically varies from two to five, and thus complete prospecting of a 40acre tract requires drilling 16 core holes and beneficiation testing of about 48 ore samples. Florida practice is not appropriate for all phosphorite deposits. Functional drill rigs and trained drillers may not be readily available. Water for core drilling may be unavailable and the characteristics of either the overburden or ore may require careful attention and excessive cost to obtain representative core samples. Alternate ore sampling methods include auger drilling, rotary drilling to obtain cuttings, and the excavation of trenches, pits, and shafts. Drill cuttings samples are generally the least expensive to obtain, but they are usually unsuitable for beneficiation testing because of particle size reduction during drilling. Ore samples from surface trenches and pits are expensive and may not be representative of deeper ore. Ore samples from shafts are the most expensive to obtain, but if taken properly they are reliable for use in laboratory and pilot plant testing. Regardless of the prospecting technique used, reliable estimates of the yield and grade of phosphate rock that can be extracted from an ore deposit are required: • to justify capital expenditures for a new mine • for mine design and mine planning • for beneficiation plant design • for production forecasts, quality control, and marketing. This paper describes models that predict the beneficiation response of ores from three phosphate ore deposits. The purpose of these models is to obtain maximum utilization of available prospect data with the least cost. ORE DEPOSITS Simple beneficiation models have been developed for ores from three sedimentary phosphate deposits around the Mediterranean area. The first deposit, herein described as "A" ore, is located north of the Mediterranean Sea. The second deposit, designated as "B" ore, is located east of the Mediterranean Sea, while the third deposit, "C" ore, is located south of the Mediterranean Sea. Ore "C" is from the Sra Ouertane deposit in northwestern Tunisia.

Jan 1, 1993

By A. Paulin, J. Roth, S. Spaic

This experiment was carried out to see whether copper-iron ingots from Late Bronze Age hoards could have been produced by smelting ferruginous copper ores. Since ferruginous copper-ores are chalcopyritic, chalcopyritic concentrate from a present day copper mine was dead roasted in our experiment in lieu of a rich chalcopyritic ore. Charcoal was used as a reducing agent in the reduction of calcine (i.e., a product of dead-roasting), and quartz sand was used as a flux. Two separate metallic solutions, a copper-base and an iron-base alloy, were obtained as a result of two factors: the presence of carbon and the miscibility gap in the Cu-Fe-C system. The metallic melts appeared in a form of bigger or smaller droplets, enclosed by a molten slag which, during the cooling, solidified into metallic beads. The presence of both metals was noticed in the majority of the examined beads. It was therefore concluded that calcine consisted of an intimate mixture of oxides of both metals which were simultaneously reduced during the reduction process. Due to the miscibility gap in the Cu-Fe-C system, separate iron and copper regions were formed in the beads. In both regions, sulphide inclusions were found. The formed molten copper and molten iron solution were found to be saturated solutions. The metallographic analysis revealed the precipitation of iron out of copper and that of copper out of iron during the cooling and solidification of melts. The formed slag contained little iron and no copper oxide.

Jan 1, 2001