Search Documents

By J. P. Gu, C. Beckermann

"Melt convection and macrosegregation in casting of a large steel ingot is numerically simulated. The simulation is based on a previously developed model for multicomponent steel solidification with melt convection, and involves the solution of fully coupled conservation equations for the transport phenomena in the liquid, mush, and solid. Heat transfer in the mold and insulation materials, as well as the formation of a shrinkage cavity at the top are taken into account. The numerical results show the evolution of the temperature, melt velocity, and species concentration fields during solidification. The predicted variation of the macrosegregation of carbon and sulfur along the vertical centerline is compared with measurements from an industrial steel ingot that was sectioned and analyzed. Although generally good agreement is obtained, the neglect of sedimentation of free equiaxed grains prevents the prediction of the zone of negative macrosegregation observed in the lower part of the ingot. It is also shown that the inclusion of the shrinkage cavity at the top is important in obtaining good agreement between the predictions and measurements.IntroductionMacrosegregation in large steel ingots is one of the most well-known and now classical problems in the field of solidification and casting [1]. Commonly found macrosegregation patterns in steel ingots are a positively segregated zone at the top, negative segregation near the bottom in the equiaxed zone, inverse segregation near the ingot surface, V-segregates along the centerline, and A-segregates in the columnar zone [2, 3]. It is well established that the macroscopic transport of species leading to macrosegregation can occur by a variety of mechanisms, including thermally and solutally driven natural convection of the melt in the mushy zone, flow due to solidification contraction, and the sedimentation of free equiaxed crystals. Fredriksson and coworkers [ 4,5] provide recent overviews of the solidification transport phenomena and macrosegregation in ingot casting."

Jan 1, 1999

By H S. Mitri, A Sainoki

Destress blasting techniques have been widely employed in underground mines for the purpose of preconditioning highly stressed rock mass. Appropriate applications of the technique successfully contribute to increasing mine safety (eg mitigating the risks for rock bursts). This paper presents a simulation technique to estimate the extent of damage zones induced by destress blasting. A 3D numerical model with a single blasthole is constructed with FLAC3D code, whereby static analysis is carried out in which horizontal and vertical stresses are applied to the model in order to simulate stress conditions in a deep underground mine. Afterwards, dynamic analysis is performed whilst applying time-varying blast pressure along the blasthole. The blast pressure profile takes into consideration an instantaneous increase and gradual decrease in blast pressure as well as the propagation of detonation along the blasthole. In this paper, constitutive models allowing for shear and tensile strengths of rock mass that are dependent upon strain rate are newly implemented into FLAC3D code with C++ programming language, whereby blast-induced damage is assessed by the extent of tensile and shear failure zones. The results obtained from the dynamic analysis show significant decrease in the extent of yielding zones induced by destress blasting with increasing the magnitude of in situ stress. They indicate that considering in situ stress state is indispensable in order to propose proper blast design for each level in underground mines. Based on these results, the relationship between the extent of the damage zone and depth at which destress blasting is performed is established. Furthermore, stress dissipation induced by destress blasting is examined with static analyses, in which the modulus of elasticity is decreased according to a rock fragmentation factor, a, only for the area where tensile and shear failures takes place during the dynamic analyses. A model parametrical study with respect to a is undertaken. The relationship between the rock fragmentation factor and the stress dissipation in the yielding zones is developed from the parametrical study. The developed relationship can be useful when destress blasting is designed for massive rock masses.CITATION:Sainoki, A and Mitri, H S, 2015. Simulation of damage zones induced by destress blasting, in Proceedings 11th International Symposium on Rock Fragmentation by Blasting, pp 245–256 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 24, 2015

By Ji Hoe Kwon

Hindered-settling column process has been used for the separation of mineral particles for many years. In this study, a model based on DEM (Discrete Element Method) was used to simulate the movements of particles in a settling column separation process. DEM model is capable of describing the motion of particles and calculating the impact forces between particle and particle, and particle and wall contact under various operating conditions. Drag force, rotation induced Magnus lift force, contact force and the gravitational force acting on individual particles were considered for establishing the mathematic models. The overall motion of the particles in water is satisfactorily simulated by the DEM, and the result shows good agreements with the result of lap-scale experiment.

Jan 1, 2014

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

The semiautogenous grinding (SAG) mill is one of the most important pieces of equipment in the comminution of materials. Its objective is to reduce the size of the ore until the ore can be released and continue on its transformation process. Although it is a simple process that involves impact, abrasion and attrition, its energy consumption is extremely high and inefficient, with no more of 17 percent of the energy supplied going toward size reduction. This gives rise to an important challenge because any improvement generates great economic benefits.

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 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 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 Jie Ma, Bo Wang, Jieyu Zhang, Jian Zhao, Zheng Chen

"In this paper, solidification process of steel ingot was simulated by ProCAST software, based on cellular automaton-finite element (CAFE) model. A 22 t steel ingot was divided into continuous small zones (100x150x180 mm) in radial direction, then the temperature field, shrinkage cavity and porosity, particularly, the microstructure were simulated. The results show that: there was a columnar-to-equiaxed transition in microstructure on radial direction. Zones, closing to the center of the ingot, had a more developed and fined equiaxed zone, while zones closing to the periphery of the ingot had a more developed columnar zone. It indicated that temperature gradient has a great influence on microstructural morphology. Thus, the simulation results were highly agree with the experimental results, verifying the feasibility of this simulation method. It significantly reduces the time of simulation, and has important significance for predicting the microstructure of large steel ingot.IntroductionSolidification microstructure of steel ingot have a significant effect on the quality of the final steel products. Although many experimental works have been carried out in order to obtain the optimum condition for solidification process of continuous casting, it is also very desirable to establish mathematical models to simulate the evolution of solidification structure. The simulation methods can be divided into deterministic models and stochastic models, which have been developed rapidly over the last couple of decades. Cellular automaton method is one of the stochastic models. Rappaz and Gandin [1, 2], Nasatac and Stefanescu [3, 4] are the earliest researchers engaged in the this field. Moreover, Rappaz and Gandin simulated the solidification structure and CET (columnar grain to equiaxed grain transition) of metals using the CAFE method (cellular automaton-finite element method) for the first time [1, 2]. Recently, the CA (cellular automaton) model also considered the effects of fluid flow on the solidification process. Zhu, et al. [5-11] improved method of CA , they used the KGT model to describe the undercooling needed during the growth of the dendrite tip."

Jan 1, 2014

By G. Cailletaud, I. Guillot, D. Larouche, R. Martinez, D. Massinon

The precipitation of Al2Cu particles in a 319 T7 aluminum alloy has been modeled. A theoretical approach enables the concomitant computation of nucleation, growth and coarsening, based on an implicit scheme using the finite differences. The equation of continuity is discretized in time and space in order to obtain a matricial form. The inversion of a tridiagonal matrix allows us to determine the evolution of the size distribution of Al2Cu particles at t + dt. The fluxes of solute in-between the boundaries, as well as the ones at the boundaries, are computed in order to respect the conservation of the mass of the system. The main results of the model are compared to TEM measurements. Simulations provide quantitative features on the impact of the cooling rate on the size distribution of Al2Cu particles in good agreement with the TEM measurements. This metallurgical model is integrated into FEA software in order to perform mechanical computations. This kind of multiscale approach allows new perspectives to be investigated in the process of designing cylinder heads. It enables a more precise prediction of the microstructure, its evolution as a function of continuous cooling rates and its impact on the integrity of the structure.

Jan 1, 2015

By A M. Wala, H W. Wu

A research project incorporating a number of mine site exercises to reduce the effects of fire incidents and possible consequent health and safety hazards has been undertaken focused on the application of mine fire simulation software packages for contaminate tracing and fire modelling in coal and metalliferous mines. This paper examines aspects of introduction of inert gases to underground workings to aid recovery of a mine following a fire. Broad conclusions from work undertaken at individual Australian coal mines are discussed as examples. The effort is built around the introduction of the fire simulation computer program æVentgraphÆ to the Australian mining industry and the consequent modelling of fire scenarios in selected different mine layouts. Case studies have been developed to examine usage of inertisation units and particularly application of the GAG jet engine unit. One example has focused on selection of the best surface portal location for placement of the GAG for most efficient suppression of a fire. A second has examined a situation with significant seam gas being emitted on the face. This has shown that under certain face dip angles stopping the mine surface fan to reduce dilution of GAG exhaust gases will cause reversal of face air and consequent mine explosion as gas laden air is drawn across a fire. A third examines inertisation and dilution issues in mains headings. mains headings present a complex ventilation network with often numerous parallel headings, hundreds of cut-throughs and a variety of ventilation control devices. In such a complex system (with additional interference from a fire), maintaining control of the movement of inert gas is more difficult than elsewhere in the mine. Some illustrations of this issue are given.

Jan 1, 2005

By Gary E. McIntosh

This paper will discuss a U.S. Bureau of Mines project designed to assess and predict the impacts which surface coal mining will have on the regional shallow groundwater systems in the Powder River Basin. The primary objective will be met through the accomplishment of the following: (1) The design and establishment of a regional network for monitoring of shallow groundwater systems in the Basin. (2) Monitoring of the direction and rate of flow, water quality, and aquifer characteristics of the shallow groundwater systems. (3) Development of a regional groundwater simulation model with appurtenant subroutines to determine local effects.

Jan 1, 1978

By A. V. Reddy, C. Beckermann

"The DC continuous casting of a grain refined AI-4.5%Cu round ingot with realistic boundary conditions is numerically simulated. Of particular interest are the effects of the transport of free solid in the form of equiaxed grains on the macro segregation patterns and the final grain size distribution. The study uses a previously developed two-phase model that incorporates descriptions of heat transfer, solute redistribution, melt convection, and solid transport on the system scale with microscopic models of solutal undercooling, microsegregation and grain impingement. The numerical results indicate that grain transport inhibits the formation of inverse segregation in the ingot shell and promotes negative centerline segregation. This is explained by the reduction in the solid coherency next to the mold and the advection of solute-lean grains towards the centerline. Comparisons with the experiments of Finn et al. [5] illustrate the need for considerable improvements in the model. IntroductionMacrosegregation continues to be one of the concerns in DC continuous casting of aluminum ingots. The variation in composition over the ingot cross section results in the need to adjust alloy composition in order to meet mechanical property specifications in finished products. The extent of macrosegregation is greatly affected by the casting conditions, such as casting speed, superheat, cooling rate, mold size, alloy composition, and grain refining practices. Thus, it is of great interest to develop a comprehensive mathematical model to simulate the casting process, enabling the selection of conditions that minimize macrosegregation while maintaining a high productivity."

Jan 1, 1994

By J. Samir

The objective of this study is to simulate the isothermal filling up of moulds by a RTM process when adopting CV/FEM and VOF method and when also taking into account the presence of obstacles coupled to the effect of the thickness and reinforcement layers variation. The codes developed intended to ensure an optimal positioning of the injection points and of the vents and also to determine optimal injection pressures, in order to minimize the cycle times. These codes use a mesh generator proper to the team, allowing a discretization of the domain in triangular unstructured elements, with the possibility of a local refinement and inserts inclusion.

Jan 1, 2009

By Tie-jun Tao, Ming-sheng Zhao, En-an Chi

"Karst landform is a special kind of geomorphic unit which is often accompanied by a large number of geological defects: caves. In blasting engineering, if there are karst caves around the blast-holes, it may result in higher boulder yield and the generation of flying rocks due to the change of burden which is very hazardous for the safety of surroundings. The models of rocks with karst caves from blasting holewere built by the software of ANSYS/LS-DYNA. The results of computation simulation show that thecave can increase more free surface for blasting and leads to a broken zone around itself. Thus rockmass around the cave will fly toward inside and it will result in undesirable result of whole rock mass to be blasted. In addition, the models of rocks with karst caves in different shapes were simulated and thecomputation results show that the caves in round shape have higher load resistance capability than that the caves in rectangle shape. All the simulation results can offer effective guidance for blasting design. The models of rocks with karst caves were built in different distances and shapes from blasting hole by the software of ANSYS/LS-DYNA. The effect of different karst caves on blasting stress and blastingresult were studied."

Jan 1, 2017

By B. DANIEL MARJAVAARA, I. A. SOFIA LARSSON, ANNA-LENA LJUNG

One way to upgrade iron ore is to process it into pellets. The pelletizing process is multifaceted, involving complex fluid dynamics with highly turbulent, high-temperature flow in intricate geometries. This makes it hard to perform measurements during operation, so a useful tool for better process understanding is computational fluid dynamics (CFD), where the challenge is to perform efficient and trustful simulations of complicated problems. A validated simulation model can be valuable for increased process control and development. In this paper, the flow field and coal combustion process in a pilot-scale iron-ore pelletizing kiln is simulated using a CFD model. The objective of this work is to apply a coal combustion model to enable an investigation of thermal effects from the flame on the surrounding flow field.