Search Documents

By Malsuzaki. Kenji

A mathematical simulation model of the Imperial Smelting Furnace operation has been developed for purpose of improving the furnace operation. The model consists of top, shaft and bottom regions. Analysis in the shaft region is based on differential heat and materials balances taking account of reaction kinetics and heat transfer rates, and those in the top and the bottom regions are based on heat and materials balances under a steady-state condition. Longitudinal distributions of process variables such as temperatures of gas, sinter and coke, gas composition and fractional conversion under a practical operating condition were calculated. With this model, the furnace performance under various operating conditions also can be shown reasonably predicted.

Jan 1, 2003

The paper examines the design issues associated with the integration of mineral processing within underground mining systems. It demonstrates conceptual integrated system design options appropriate to a range of geological scenarios. Particular reference is then given to a case study of stoping in narrow base metal veins. The performance of such integrated systems in terms of productivity, economic and environmental criteria is demonstrated by means of simulation. Particular attention is given to the potential benefits to mine waste management. Processing technologies appear to offer genuine potential benefit, whilst underground mining methods may need to evolve further to facilitate full integration. This paper is based on industry supported research aimed at developing selective and clean underground mining systems.

Jan 1, 2003

By J. R. Bredell, J. Branehog, C. Coetzee, J. J. van Tonder

A nutating grinding mill, exemplified by the commercial HICOM mill, achieves efficient particle size reduction through the accelerated motion of the material charge within the nutating chamber. This design, characterised by high particle velocities and rapid charge circulation, offers a compact solution compared to tumbling mills. There have been concerns regarding the mechanical reliability of kinetic joints for certain specialised mining applications. Recognising the need to model internal forces in such systems, this study focuses on the uncoupled simulation using multibody dynamics and the discrete element method. The investigation employed a small-scale experimental test mill, referred to as the NuMILL, featuring a cylindrical chamber (140 mm internal diameter, 245 mm height, 6° nutation angle). Instrumented with an array of sensors, the NuMILL provided internal load data for validating the simulations. Experiments were conducted with an empty chamber or filled with charge material. Charge materials included coarse sand or steel pellets with a charge-to-chamber mass ratio of 0.25. The mill speed ranged between 100 to 700 RPM. Multibody dynamics modelled resultant inertial and frictional loads within 5%, while the discrete element method simulated the internal material contact loads with accuracies ranging from 1% to 19%. The combined simulation approach predicted the resultant forces in the critical kinetic joint within a 20% margin. Representing the charge using rigid lumped mass elements was also studied. The presented simulation techniques, validated against experimental data, offer a robust foundation for analysing larger systems, including the industrial-scale HICOM.

Mar 5, 2025

By Guangyuan Li

A 3-D block-spring model (BSM) has been developed for simulating the behaviour of jointed rocks. The model uses a discontinuous formulation and an iterative procedure to solve the large displacements in jointed and faulted rock masses. The 3-D BSM has been used to simulate mining excavations in jointed rocks around an open pit and underground settings. The applications have demonstrated that it is a practical tool to solve complicated problems of rock mechanics in civil and mining engineering. This model has recently been modified to incorporate complex geological conditions, such as underground water flow (seepage), unstable blocks and the time-dependent behavior of backfill materials. To demonstrate the capability of the 3-D model this paper presents the fundamentals of the 3- D BSM and its application for open-pit mining. A complete computer package BSM3D has been developed based on the 3-D block-spring model. The package consists of a solving engine for numerical analysis and a Windows based processor for mesh generation and results display. The multi-functioned pre-processor provides tools and menus to help users setting up models easily and conveniently. The fully featured post-processor allows users to view the results in different windows and orientations. The graphical results include contours, vectors, animation of rock movement, and the failure state of the rock structures.

May 1, 2001

By A. De1

Because of the extremely low utilization of strain energy in the breakage of solid particles, the energy efficiency prevalent in comminution technology is extremely poor. Schoenert considered that single-particle comminution under slow compression is the most energy-efficient mode of comminution and that energy efficiency should be compared to that. The usual method for grinding minerals and cement has been the ball mill, in which the hit-and-miss aspect of loose-bed comminution severely reduces energy efficiency. The fine grinding of minerals is a particularly energy-intensive process. Schoenert showed that the energy consumed in the comminution of particles in confined beds is second only to that of single-particle comminution under slow compression and developed the method of conducting particle-bed comminution continuously in a high-pressure roll mill. Since its first industrial installation in 1985, the high-pressure roll mill has established itself as the comminution device of choice for significantly improving the productivity and energy efficiency of size reduction from coarse size feeds down to products finer than 20 microns. The first of several hundred installations were in cement plants. The preliminary design and optimization of the industrial operation of high-pressure grinding roll (HPGR) systems, owing to the continuous nature of the device, cannot be efficiently done through laboratory experiments and needs to be achieved through proper system modeling. While the energetics of size reduction in high-pressure grinding rolls has been modeled very well, no model exists for the macroscopic variables in HPGR comminution. In high-pressure roll mill grinding, size reduction results from confined particle-bed comminution under compressive stresses such that the daughter particles produced are highly fractured, stressed, or otherwise weakened. This paper attempts to relate the comminution / compaction behavior for confined particle-bed under compression in a piston-die press with the macroscopic variables in HPGR comminution. A model describing the pressure-densification relationship is presented and used to accurately simulate the macroscopic variables in high-pressure grinding rolls, utilizing experimental data for particle-bed comminution with the piston-die press. The accuracy of the predictive simulations strongly suggests that, starting with a few batch-mode laboratory-scale experiments, it should be possible not only to identify the optimal operating conditions for existing operations but also to successfully design new industrial high- pressure grinding roll systems.

Sep 1, 2012

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

By J. Rajek

The novel multi-component, multi-particle, multi-phase precipitation model developed recently is used to predict the precipitation kinetics in complex 9-12% Cr steels investigated within the European COST 522 project. These steels are used for tubes, pipes, casings and rotors in USC (ultra super critical) steam power plants for the 21st century. In the computer simulations, the evolution of the precipitate microstructure is monitored over the entire production heat treatment starting with the casting process, including austenizing, several annealing treatments and in-service conditions. The main interest lies on the concurrent growth and coarsening interaction of different types of precipitates. Keywords: Precipitation kinetics, heat treatment, power plant steels, simulation

Jan 1, 2003

By M. Mathey

The application of the Monte Carlo technique to production planning and everyday economic decisionmaking in mine production management is demonstrated. The logic is detailed using an example of underground production with continuous miners (CMs) and truck haulage. It is argued that availability of equipment and personnel are the predominant variables influencing mine output and productivity and that those availabilities may be well represented by binomial probability distributions. The probabilistic model is implemented in a standard Excel® spreadsheet with Palisade’s @Risk add-on to facilitate simulations. Starting from model calibration against data obtained from a mine’s annual reports, some general interdependencies of availability, utilization, productivity, and costs of production processes are outlined. Finally, several possible options and their consequences as regards production improvements are explored.

Dec 13, 2022

By Paul W. Cleary

Keywords: Reactive pellets, CFD, Lagrangian SPH, Pyrometallurgy In several pyrometallurgical processes, reactive pellets are added to a melt. Heat is transferred to the pellets and reaction products such as gases, metal and slag are formed. For many pyrometallurgical processes traditional grid based CFD methods, such as using the Finite Element Method (FEM) and Finite Volume (or Control Volume methods), produce suitable predictions. However, the inclusion of solids in the bath represents significant challenges for these methods. Smoothed Particle Hydrodynamics (SPH), a Lagrangian simulation method, is able to simulate both the fluid component of the bath and any immersed solid materials, tracking the motion of these solids, and their interaction with the fluid, including gas generation. In this paper, we will summarize the key aspects of the SPH method and will show a series of structured examples illustrating the behavior of the different physics sub-systems that are used in the model of the overall system.

Jan 1, 2005

By S. Kelebek, A. Khan, J. Caceres, Beckett. D.

"Pyrrhotite in nickel-copper sulphide ores is nickeliferrous and intimately associated with pentlandite. Due to its low nickel content, pyrrhotite rejection is an important step in processing of these ores. In our laboratory, various reagents were investigated for their effect on pyrrhotite depression. Among these, a combined use of triethylenetetramine and sodium metabisulphite was found highly promising.The work described in this paper has two parts. The firs part is a summary of results obtained on process samples from the pyrrhotite rejection circuit at Falconbridge’s Strathcona Mill. The second part involves simulation of the circuit performance based on the reagent conditions used in the experimental work and the concept of split factors.Two sets of reagent conditions were evaluated. The first one involved baseline conditions, simulating the current reagent performance. The second one simulated the impact these new reagents can have on the circuit performance.The simulation results obtained under baseline conditions indicated a good agreement with the current circuit performance. The simulation results suggest that at practically the same pentlandite recovery and tailing grades, the circuit concentrate nickel grade can be doubled through effective pyrrhotite rejection by these new reagents."

Jan 1, 2003

By Johannes Labenski, Christian Moormann, Johannes Aschrafi

"Open steel pipe piles are used for various applications in costal engineering, port structures and increasingly important for offshore structures. During the installation of open steel pipe piles the formation of a plug has an influence on the installation process of the pile and also on the final bearing capacity. The formation of the plug depends on different factors, e.g. the pile diameter and the installation method. This paper starts with a structured overview about analytical methods (API, EAP-EAU, UWA-05, HKU-12, Modified API) to calculate the axial bearing capacity of open steel pipe piles. The numerical simulation of the installation process of open steel pipe piles has to fulfill high demands, i.e. to accurately represent the penetration of the pile inside the soil, the formation of the plug and also the stresses and strains during the installation. A numerical back analysis of a comprehensive field test, using the Coupled Eulerian-Lagrangian (CEL) method and a hypoplastic constitutive model for sand, reveals the stresses and displacements of the soil and of the piles. The analysis provides a comprehensive understanding of soilpile interaction.INTRODRUCTIONThe particular difficulty in the numerical simulation of plug formation, is that the effect of plugging inside open ended steel pipe piles has not been fully investigated. In scientific papers different approaches are documented to explain this phenomenon. Some are based on in-situ tests and others are based on numerical simulations. One of the challenges is to determine the end bearing capacity and shaft friction of the steel pipe pile. Different approaches have been developed over the years to determine them e.g. from Randolph et. al. (1991, 1994) and Paik et. al. (2003). Recent approaches are the ICP-05 Method (Jardine 2005), the UWA-05 Method (Lehane 2005) and the HKU-12 Method (Yu & Yang 2012a). A comparison of some methods can be found in Schneider et al. (2010). A plug inside an open steel pipe pile can be visualized as a spatial bracing of soil between the inner surfaces of the pile, cf. Lüking (2010).. However according to Lüking & Kempfert (2012) this bracing only exists for the height in the end of the pile equal to twice its diameter.To identify a plug inside a pile, there are different measurable parameters. One parameter is the development of stresses around the pile during and after its installation. Henke (2013) measured the inand outside horizontal stresses at the pile tip during the installation. The results showed higher stresses on the inside compared to the outside. To validate his measurements he performed a numerical back-analysis of the measured installation process and got matching results."

Jan 1, 2017

By SenYang Qian, ShuangYu Du, Bo Wang, Jieyu Zhang, Jian Zhao

During the solidification process, the flow pattern of the molten steel has significant influence on the formation of solidification microstructure. In present study, the solidification microstructure of 30Cr2Ni4MoV steel ingot under different intensities of mechanical oscillation(forced convection) was simulated by using ProCAST®software. The results showed that the increase of the intensity of oscillation could restrain the development of the columnar grains and refine the grains. As the intensity of oscillation increased, the proportion of columnar grains decreased and the proportion of equiaxed grains increased. When the intensity of oscillation was at a low level, the effect of refining grains could be improved obviously by increasing the intensity. However, as the intensity increased to over 0.833, the effect of refining grains would not be improved.

Jan 1, 2016

An understanding of the way that solutes are transported in mining or natural environments is essential to the development of sound strategies for managing waste disposal. Research on miscible pollutant movement in unsaturated porous media is carried out both experimentally and numerically. Results of numerical modelling are obtained by solving the groundwater flow equation and the solute transport equation. The finite element method is employed and comparisons between simulation results and experimental data are found in remarkable agreement.

Jan 1, 1995

By A. C. Smith

Although only used in a few U.S. longwall mines, bleeder-less ventilation systems can be an effective spontaneous combustion control method in mines having a demonstrated history of spontaneous combustion. To provide insights for the optimization of bleederless ventilation systems for U.S. underground coal mines, a computational fluid dynamics (CFD) study was conducted to model the spontaneous heating in longwall gob areas using a bleederless ventilation system. A single longwall panel with a bleederless ventilation system was simulated, and typical longwall mine ventilation data were used in the simulations. The permeability and porosity profiles for the longwall gob were estimated using a geotechnical model and were used as inputs for the CFD modeling. The effects of gob permeability and resistance of the collapsed entries on the spontaneous heating were studied. The effectiveness of using nitrogen injection to prevent spontaneous heating in the gob was also examined.1

Jan 1, 2008

By V. K. Karra

"A simplistic algorithm jar simulating the rod or ball mill Bond grindability test applied to ores and minerals is presented. A performance comparison between the present algorithm and a previously published algorithm by Kapurot, modified in this work jar a fair evaluation, is also given. A worked example is furnished jar leading the reader through the various calculations.IntroductionThe Bond grindability test is widely used for the prediction of the power requirement of an industrial grinding mill. The procedures for conducting this test in either a ball or a rod mill are available(1). In the case of a ball mill, batch grind abilities have been correlated to Bond grindabilities to help monitor day-today variations in ore grindability for the purpose of plant control( 2,3). Of special interest is the algorithm presented by Kapur(3) for the simulation of the Bond test. In practice, the test lasts from 6 to 14 cycles (4 to 10 hours). With this algorithm, one can predict the Bond grind ability using information available from only the first two cycles.Sometimes the prepared feed for the grindability test contains too much fines compared to the total product desired at the mesh of grind under consideration. In such a case, the operator separates the fines from the test feed and adds to the remaining coarse fraction an equal amount of fresh feed. This combined feed then forms the feed to the first cycle of the test. In the present work, Kapur's algorithm is modified to include this special case. A new and relatively simpler algorithm is also presented. Statistical equations relating work index to pertinent variables are also given."

Jan 1, 1981