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

By M Phil Schwarz

"A CFD model for the gas-side of an impinging jet iron bath reactor has been developed to assist in the interpretation and design of oxidizing gas-melt reaction rate experiments. Heat transfer and reaction at the surface have been included. A free surface model, which solves for the liquid motion within the bath as well as the gas motion above the bath, has also been developed. This model also predicts extent and shape of the dimple formed on the liquid surface by the impinging jet. Both models have been run for conditions typical of jet impingement experiments using CO2 , H2O and O2• The main conclusions are• Results from the model for gas-side mass transport are in good agreement with the empirical correlation of Belton and Belton (1980)• Experimentally determined rates of decarburization by CO2 should be corrected for mass transfer because the experiments are often carried out under mixed rate control, eg for gas flow rate 10 L/min when the sulphur concentration is lower than 0.045• If the rate constant for oxygen reaction is greater than about 5 mol/m2.s.atm, it cannot be determined using a gas flow rate of 10 L/min or lower because of mass transfer control.• A cavity of depth 4.5 mm is predicted for a gas flow rate of 10 L/min and this substantially changes the flow field (and hence reaction rate) compared with a simulation in which the surface is assumed to be flat."

Jan 1, 1999

By John ?gren, Martin Schwind, Henrik Larsson

"A new model for internal oxidation, as well as similar processes, has been implemented on computer. It is based upon a combination of random walk and cellular automata and is capable of treating 3-dimensional cases. To model diffusion, it utilises quantised units of matter. All steps of the process are included: nucleation, growth and coarsening. The effect of grain boundaries on diffusion as well as nucleation is taken into account. Further possible development is briefly discussed; especially the treatment of nucleation needs further consideration. IntroductionComplex diffusion-controlled processes may now be simulated on a routine basis when concentration varies along one co-ordinate only. One-dimensional simulations capture much of the information of interest but not all. For the case at hand, especially the spatial and the size distribution of precipitates, as well as the impact of grain boundaries, would be interesting information not accessible from 1-D simulationsSolving the full Fick-Onsager law in three dimensions does not seem feasible with the computers of today. In the model described in this article, random walk by quantised units of matter are considered. The model should be useful for a number of cases involving dispersed systems. With further development, it would indeed seem possible to simulate most types of phase transformations.Previous workWe have previously published two articles [1,2] on the new approach discussed in this work. In those articles we described the basic idea and showed, among other things, that very good agreement with an experimental concentration profile can be obtained, in a single-phase system, using independently assessed thermodynamic and kinetic data. Also, Schwind used the same random walk based concept to model coarsening [3]."

Jan 1, 2001

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

Under normal circumstances, methane released along a longwall face represents only 10%-20% of the total me- thane emitted from the entire longwall panel. However, there might be a sudden surge of methane at the face which may cause a potential problem because the friction from shearer cutting, support movement, and scraper conveyor dragging can be ready ignition sources. Therefore, to predict the distribution of methane will be critical to prevent possible explosions. In this paper, methane dispersion patterns along a lognwall face are simulated under normal and abnormal conditions. With the simulation of methane released from the face, around shearer drum, and from the broken coal on the conveyor, places with possible high methane content can be identified by different phases of working process and shearer positions. Local ventilation measures are evaluated so precautionary measures can be planned ahead.

Jan 1, 2012

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 Baicheng Liu, Houfa Shen, Wensheng Li

"In order to understand the formation of macrosegregation in steel ingots, a mixture model has been developed which accounts for coupled mass, momentum, species, and energy transport to describe the sedimentation of free equiaxed grains during solidification. In this model, the mushy zone is divided into the slurry and porous regimes according to a solid packing fraction. Settling velocity of equiaxed grains in the slurry regime is given by an explicit expression. The solidification of a binary steel ingot (Fe-0.36wt.%C) with a two-dimensional geometry is numerically simulated. It is demonstrated that the predicted macrosegregation patterns agree to some extent with classical experimental results. Further refinement of this model is still necessary for realistically simulating steel ingot solidification.IntroductionModeling macrosegregation in steel ingots is still a challenging study, mainly due to the length scale of the products and the complexity of global multiphase transport phenomena coupled with microscopic solidification kinetics [l]. The thermal-solutal convection and the motion of equiaxed grains are the two main mechanisms responsible for macrosegregation. But there are few models that deal with the motion of equiaxed grains, which is confirmed as the main phenomenon that explains the negative macrosegregation along the centerline of large steel ingots [2]. A model which takes into account melt convection and grain sedimentation was first proposed by Beckermann and co-workers ([3], and references therein). More recently, similar multiphase solidification models with advanced nature have been presented by Wu and Ludwig [1,4] and Combeau and co-workers [5,6]. It is known that the application of these sophisticated models to steel ingot casting is at its beginning.In this paper, a solidification model based on continuum mixture approach, using much less sophisticated relations to obtain the settling velocity of free equiaxed grains, has been proposed in order to predict the macrosegregation in steel ingots. The model is applied to simulate the solidification of a binary steel ingot with a two-dimensional geometry, and preliminary modeling results are presented."

Jan 1, 2011

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

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