Search Documents

By R. K. Rajamani, A. Datta

An appropriate scale-up methodology is very useful for ball mill design and optimization of the grinding circuit. To date, the scale-up is based on the selection and breakage function model using lab mill data. However, a major drawback of that approach is that the charge motion in a lab mill is vastly different from that in a large industrial mill. Therefore, in the methodology proposed here, first the energy of individual collision in the mill is computed employing the discrete element method (DEM). Then the breakage of particles as a function of collision energy is studied in a separate experiment. Together, this information constitutes what is known as collision energy based population balance model. The appropriateness of this model is shown for different batch grinding mills.

Jan 1, 1999

By A. Datta

An appropriate scale-up methodology is very useful for ball mill design and optimization of the grinding circuit. To date, the scale-up is based on the selection and breakage function model using lab mill data. However, a major drawback of that approach is that the charge motion in a lab mill is vastly different from that in a large industrial mill. Therefore, in the methodology proposed here, first the energy of individual collision in the mill is computed employing the discrete element method (DEM). Then the breakage of particles as a function of collision energy is studied in a separate experiment. Together, this information constitutes what is known as collision energy based population balance model. The appropriateness of this model is shown for different batch grinding mills.

Jan 1, 1999

By N. R. Green, N. J. Humphreys, D. McBride, T. N. Croft, M. Cross, D. M. Shevchenko

"Centrifugal casting offers a route to high quality products in difficult to cast high temperature low superheat alloys and thin section molds. Under centrifugal forces metal is forced into thin sections and can fill thicknesses of less than a millimeter. However, due to the high liquid metal velocities there is a high risk of surface turbulent flow and air entrainment within the liquid metal. The combination of interacting flow-thermal-solidification phenomena and associated defects is a challenging modeling task which the authors have previously described and validated. Capturing the metal-air interface, on what are inevitably complex three dimensional geometries, results in highly computationally expensive simulations and simulating a single cast can take weeks on a single processor. This contribution reports on modeling a complex centrifugal cast, gas entrainment, bubble transport and solidification, employing meshes of up to a million elements and investigates the scalability of the model on high performance clusters.Introduction The ultimate goal of any casting procedure is to produce a defect-free product. For difficult to cast alloys, such as Titanium, centrifugal casting has a number of advantage over conventional gravity methods [1]. In the casting of thin sections, such as turbine blades, the metal is forced into the thin sections under centrifugal forces. The centrifugal force can overcome the backpressure due to surface tension and the metal can fill mould sections with thicknesses substantially less than a millimeter. However, the materials used in these castings [2] tend to be highly reactive and requires low superheat. High liquid metal velocity gives rise to a high risk of free surface turbulent flow and of associated entrainment through free surface turbulence [3] and trapping of gas present in the flow [1]."

Jan 1, 2012

By A. Sand

Within stamp-charged coke making a large volume of coal is compacted to one single coal cake before entering the coke oven chamber. This is done by means of several falling stampers in a stamping machine having a mold nearly of the oven?s dimensions. Producing a high quality coke also from inferior coals requires a cake density of approximately 80% of the coal density. Besides that the industrial stamp-charging process demands a minimum mechanical strength of the coal cake to move it from the stamping mold into the oven chamber without failure. Densification and the build-up of cake strength were investigated earlier in stamping tests using a micro-stamping device. The quantities derived from these tests (as e.g. cake density or porosity) represent average values for the entire cake. Statements describing the local compaction state at different heights are difficult to achieve without destructing the cake. In order to gain better insight into the densification process and the inner structure of the coal cake a computational model based on the discrete element method (DEM) has been set up using 2 and also 3-dimensional simulation software. The stamper?s position and velocity as well as the force acting on the stamper were monitored and the model?s response was compared against measurement data from laboratory stamping tests. It was possible to reproduce the force and displacement pattern of the stamper in response to the viscoelastic properties of the cake using standard DEM bonding and contact models. Furthermore, the rearrangement of particles in response to the compaction by the stamper was tracked by calculating their displacement at the point the stamper hit the coal surface. The latter can also be used as indicators of particle deformation or breakage. By defining control points at different heights the particle displacement, stress and strain rates, porosity could be studied at different heights. Keywords: coke making, coal stamping, compaction, agglomeration, DEM, simulation

Sep 1, 2012

By Yaming Chen, Gerry Papini

This paper presents a case study of three-dimensional groundwater modeling of a proposed open pit at a gold mine project with historic flooded underground (u/g) mine workings in British Columbia, Canada. The modelling approach for dewatering the open pit included passive dewatering by draining water under natural hydraulic gradients and active dewatering by using active dewater wells, together with residual water drains. Target water level drawdowns behind the pit wall were used to determine the number and depths of dewatering wells to meet the slope stability requirements and in estimation of pit inflow rates. Conceptualizing the u/g mine workings and assigning appropriate flow boundary conditions throughout the life of the mine posed some unique challenges.

Sep 30, 2022

By James P. Bennett

EAF slag chemistry and its control have been recognized as important in producing quality steel at a low price for a number of decades. EAF steel producers make lime and dolomite additions to furnaces, controlling the CIS ratios and sulfur pickup from the steel. Existing phase diagram data at 1600°C for CaO, MgO, FeO, SiO2, and Al2O3 was used to write a computer model optimizing slag chemistry for slag foaming and reduced refractory wear. The model, its development, and its use in an industrial steelmaking environment will be discussed along with how the model can improve energy efficiency. The importance of EAF mass balance to account for critical materials like Si in scrap and refractory wear using the model will also be emphasized.

Jan 1, 2003

By James P. Bennett

EAF slag chemistry and its control have been recognized as important in producing quality steel at a low price for a number of decades_ EAF steel producers make lime and dolomite additions to furnaces, controlling the CIS ratios and sulfur pickup from the steel. Existing phase diagram data at 1600°C for CaO, MgO, FeO, SiO2, and Al2O3 was used to write a computer model optimizing slag chemistry for slag foaming and reduced refractory wear. The model, its development, and its use in an industrial steelmaking environment will be discussed along with how the model can improve energy efficiency. The importance of EAF mass balance to account for critical materials like Si in scrap and refractory wear using the model will also be emphasized_

Jan 1, 2003

By Vaughan R. Voller

In many high temperature processes involving liquid melts the common practice is to solidify a fraction of the melt adjacent to the containment wall; a practice that can provide protection of the structural components of the process. The position, shape and transient behavior of these so called "skulls" and/or "ledges" are directly dependent on the thermal process conditions. This paper presents a model of a generic ledge process that can predict the steady shape of a ledge. Two numerical implementations are presented and compared: (1) based on a space grid of continuously deforming finite elements (2) based on a conjugate grid of fixed control volumes and deforming boundary elements.

Jan 1, 1996

By B. Q. Li, S. P. Song

"Numerical models are developed to represent complex electromagnetic, free surface deformation and fluid flow phenomena in floating zone refining processes with radio frequency inductive heating. The computational methodology is based on the coupled finite/boundary element method for the electromagnetic field distribution and the Galerkin finite element method for magnetically driven fluid flow and free surface profiles of the molten zone. With the models, the complex transport and free surface phenomena in a floating zone system can be studied as a function of various operating conditions including applied current, frequency, inductor position and shape, surface tension, floating zone diameter and height. Results show that the equilibrium shape is sensitive to the height of the molten zone and that the free surface shapes can have a significant effect on the melt flow in the molten zone. Our extensive simulations also suggest that to ensure accurate results, the same mesh should be used, whenever possible, for both the electromagnetic forces and fluid flow calculations.1. IntroductionThe floating zone processes have been widely used in the materials industry to grow ultra-pure rare earth metal or semiconductor single crystals for sensor or electronic applications that have a tough demand on purity. Developed based on the principle of magnetic melting/levitation, the processes possess a main advantage that the molten zone is supported in air by the electromagnetic forces generated by the coil (see Figure 1), and therefore is free from contamination from a container. Consequently, the final products made from the floating zone processes have higher inherent purity than those by other techniques [1-4].Complex physical phenomena occur in the floating zone processing of materials. These include free surface deformation, magnetically-driven melt flow, which may be turbulent or mildly turbulent [5] depending on the applied conditions, and solidification. Thermal stresses also develop in the single crystals, which are affected by the flow and temperature field and which can have direct implications to qualities of single crystals [6]. An understanding of the fundamentals governing these physical phenomena can be of critical value in elucidating the physics that controls the properties of the metal or semiconductor crystals being produced, thereby providing a rational basis for process improvement and quality control."

Jan 1, 1997

By M. G. Basavaraj, S. Sarkar, G. S. Gupta, P. D. Patil

"Previous experimental study on raceway size hysteresis on two-dimensional cold model showed that the inter-particle and wall-particle friction had a very large effect on the raceway size. Existing literature correlations for raceway size ignore the frictional effects. It has also been shown in the present study that their applicability to the ironmaking blast furnace is questionable. To take into account the effect of friction on the raceway size a stress analysis has been done for the raceway region. Two-dimensional fluid flow equations under turbulent conditions have been developed and solved computationally in order to get the pressure forces. The frictional forces were obtained in terms of stresses. To predict the raceway size a force balance was done for the raceway zone considering the pressure force, the frictional force and the bed weight. The resulting equations from the force balance are able to describe the raceway hysteresis phenomena reasonably well along with the raceway size. A two dimensional experimental set up has been fabricated in order to validate the computer predictions. Predicted values agree well with the experimental values.IntroductionWhen a high velocity blast of air is introduced into a packed bed through a tuyere, a cavity is formed in front of the tuyere. This cavity, formed in the iron-making blast furnace, is called the raceway. Coke tumbles down and bums in this raceway, generating heat and reducing gases for the iron-making process. The raceway is responsible for producing most of the reducing gasses in the blast furnace. In addition, the size and shape of the raceway determines the movement of solids within the blast furnace, the distribution of gases and liquids (slag and metal), setting the location where the primary chemical reactions occur and the temperature profile in the lower part of the blast furnace.Due to the importance of the raceway, many investigators l - 6 have studied it and have predicted the raceway size using the empirical correlations based on experimental (hot & cold) and mathematical models. A number of correlations have been reported in the literature to predict the raceway size by various investigators3 - 6 which are different to each other. Each of the correlations shows a different dependence of raceway diameter (DR) on bed height, particle size, solid density, bed voidage and blast velocity. None of the correlations includes the effect of solid stresses on the raceway size and shape. The magnitude of these stresses will depend on the bulk solids properties (internal angle of friction, wall angle of friction) and the bed dimensions."

Jan 1, 2001

By S. A. Briggs

Bioremediation has been accepted as a treatment technique for groundwater contamination in subsurface soils and shows promise for contaminated fractured rock environments. Biological growth in fractured rock is expected to occur predominantly as biofilms attach to the fracture surfaces. Biofilms in rock fractures are subject to a complex system of forces and other phenomenon due to the dynamics of the bulk fluid in which they grow. In this paper, through the applications of computational fluid dynamics (CFD) to rock fractures, where the boundaries are rough and the flow is complex, a precise analysis was conducted of the interaction of a fluid flow and the rock fracture. Specifically, hydraulic parameters and velocity profiles of an actual rock fracture were calculated and compared to a fracture of equivalent aperture. From the analysis it is clear that it is important to use more complex models such as the Lattice Boltzmann Method used is this paper to describe fracture flow.

May 1, 2009

By Jaime Herrero, Santiago Gómez, Pablo Segarra, José Ángel Sanchidrián

This paper presents a mathematical model for the study of the 2D movement of blast bench faces based on the high-speed video analysis of three blasts at El Aljibe quarry (Spain) recorded at 2,000 frames per second. In these blasts, a set of bags of 0.75 × 0.75 m were placed in front of a selected blasthole defining a reference plane in which to study the kinematics of the movement. The burdens range from 3.86 m to 5.54 m, with a bench height of 12–13 m. An affine transformation system was calculated to convert coordinates from the video images into actual spatial coordinates. The problem is established from a multi-objective minimization approach in order to simultaneously optimize two defined objectives: trajectory fitting using Total Least Squares regression, and time–displacement regression using Ordinary Least Squares. Muti-objective minimization formulation produces Pareto plots of non-dominated solutions, i.e., optimal combinations of initial velocities and response times, (i.e. time between detonation and the onset of bench face motion). From the Pareto frontier, the most balanced solution is selected. Blast response time lasts only a few milliseconds and appears so brief at the frame rate used that it is difficult to detect visually. In this context, the bi-objective optimization-based approach proves to be a more robust tool than frame-by frame survey for quantifying such a rapid event while also maintaining a strong performance in the trajectory fitting task. Root Mean Square Error (RMSE) for the two minimizations average 0.22 m and 0.13 m, respectively. The initial velocities obtained range from 4 to 17 m/s, and the response times 0.7 ms to 21.7 ms. An analytical model is also presented for predicting fragments initial velocity and response time as functions of blast geometry and rock mass properties, explaining approximately 87% of the total variability.

Jan 26, 2026

By Danqing Gao, Brijes Mishra

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

Jul 20, 2020

By Shoma Sinha, Nirmal K. Sinha

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

Jan 1, 2004

By Gabriel Walton, Sankhaneel Sinha

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

Jan 1, 2017

By Alaina B. Hanlon

Keywords: Undercooling, Convection, Dendrites, Stainless Steel, Levitation Certain Fe-Cr-Ni stainless steel alloys solidify from an undercooled melt by a two-step process in which the metastable ferrite phase forms first followed by the stable austenite phase. Recent experiments have shown that the lifetime of the metastable phase is strongly influenced by flow within the molten sample. The current research will provide insight to why flow affects the metastable phase. This new understanding will lead to the use of convection to control microstructural evolution in spray-forming of stainless steels. Simulations using a commercial computational fluid dynamics package, FIDAP, yield a range of convective flow velocities that cause mechanical damage to the dendrites which is compared to experimental results. If the convective velocities are great enough that the dendrites bend, then low angle boundaries form at the points of collision. These result in high energy sites that could serve as nuclei for the stable phase.

Jan 1, 2005

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

Jan 1, 2010

By M. G. Lipsett

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

May 1, 2001

By H. Askari-Nasab, B. Osei, K. Awuah-Offei

The U.S. coal mining industry consumes approximately 142 billion kWh/a of energy. The U.S. Department of Energy estimates that the industry’s annual energy consumption could be reduced by 49%. This constitutes nearly $3.7 billion of potential savings on coal production costs at 5.3¢/kWh of energy. Additionally, with climate change regulation on the horizon, any benefits from energy savings in the near future are compounded by associated reductions in CO2 emissions. Research indicates that operating conditions significantly affect energy efficiency. Modeling and simulation have been shown to be a cheap and reliable means of evaluating the effects of process interactions and operating conditions. The goal of this work was to apply stochastic process simulation to model the energy efficiency of a typical truck-and-shovel mining system and use the model to evaluate production strategies to improve energy efficiency. The research team conducted energy audits of truck-and-shovel overburden removal and used the results to develop regression models describing truck and shovel fuel consumption. The team then built a stochastic simulation model of the truck-and-shovel overburden removal operation and used it to assess a variety of improvement measures by simulation experimentation. Valid fuel consumption models for shovel loading and truck haulage have been formulated, based on the energy audit results. Valid stochastic process models of truck-and-shovel operations have been formulated to study energy efficiency. The following strategies, in decreasing order of impact, provide the most energy savings for truck-and-shovel overburden removal at the mine: (1) shorten haul roads, (2) increase shovel capacity and (3) increase shovel utilization through optimal truck matching. Additional data will be required to adequately describe operator effects on energy efficiency.

Jan 1, 2011

By H. P. Kritzinger, T. C. Kingsle

"Rotary kilns are used for a variety of mineral processing operations. Hatch makes use of a kiln model developed from first principles to evaluate designs for its clients. This tool has been applied to a variety of pyrometallurgical applications, including ferrovanadium, nickel carbonate, nickel laterite, iron ore reduction, and spodumene (lithium) production.This paper illustrates the application of numerical optimization techniques in combination with the kiln model in the interrogation of a generic iron ore reduction process. The fundamental modelling concepts are explained, followed by a description of the optimization approach. The results show how the combination of the two methods, with modern computing power, can generate a large number of viable design and operating candidates."

Jan 1, 2015