Search Documents

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

Cyanidation tailings which are disposed in a surface impoundment experience a loss of cyanide due to natural degradation, frequently reducing the cyanide concentration to very low levels. Quantifying cyanide losses in terms of pond geometry, local weather conditions and feed solution chemistry has been largely empirical though in many cases mining operations rely on surface impoundments to reduce cyanide to below an internally regulated concentration or an effluent limitation. To permit a more quantitative evaluation of cyanide losses in an impoundment, a computer simulation was developed to estimate the losses of free, WAD and total cyanide due to dissociation, photolysis and volatilization. Results of the model are compared against data collected for a North American tailings impoundment during 1998.

Jan 1, 1999

By M. M. Botz

Cyanidation tailings which are disposed in a surface impoundment experience a loss of cyanide due to natural degradation, frequently reducing the cyanide concentration to very low levels. Quantifying cyanide losses in terms of pond geometry, local weather conditions and feed solution chemistry has been largely empirical though in many cases mining operations rely on surface impoundments to reduce cyanide to below an internally regulated concentration or an effluent limitation. To permit a more quantitative evaluation of cyanide losses in an impoundment, a computer simulation was developed to estimate the losses of free, WAD and total cyanide due to dissociation, photolysis and volatilization. Results of the model are compared against data collected for a North American tailings impoundment during 1998.

Jan 1, 1999

By Troy Williams, Chris Preston, Ian Lipchak

"Underground blasting operations are challenging from the standpoint of the distribution of explosivesenergy representative of ring blasting. Energy from both shock and pressure regimes of commercialexplosives may appear concentrated in the collar region of a typical ring blast and diluted at the toes of holes due to the oblique geometries of blastholes. The non-homogenous nature of ore in whichexplosives are distributed via drillholes, adds to the complexities of generating particulate profiles from fragmented material with consistencies that are predictable from blast pattern to blast pattern - well suited for specific underground handling equipment and mill processing. In an ideal world, it would be the blasting operations themselves that represent the primary crushing mechanism, or at least mitigate mechanical crushing that can comprise a large component of the cost in generating suitable muck. This paper presents a dynamic break view in 3D that allows a planner to visualize the potential break zone around a blasthole generated by an explosive load using a Kleine field. Simple as it sounds, this methodology provides information that can be used in conjunction with cavity monitoring surveys (CMS) to potentially judge dilution due to overbreak as well as recovery for a typical blast. As examples, there are two break geometries that are examined regarding circular breaks and elliptical breaks around blastholes. Using a Kleine field to define break, a planner generated isosurface can be generated and compared to CMS data for calibration and prediction, using AEGIS 3D ring design software."

Jan 1, 2016

By M Reza Afshar

In this work a mathematical model was developed to study the effect of electromagnetic force on the magnesium melt containing nonmetallic inclusion in an electrical isolated channel. The electromagnetic force exerted on the inclusion induced by DC current field was calculated. A three dimensional turbulent flow equations were solved numerically to compute the velocity field within the channel. On the basis of computed velocity and magnetic fields, the trajectories of inclusions were calculated using the equations of motion for inclusions. Parametric studies were carried out to evaluate the effect of various parameters such as current density, channel diameter and particle density on the removal efficiency of inclusions from magnesium melts.

Jan 1, 2004

By Nagy El-Kaddah, T. T. Natarajan

This paper describes a novel approach for modeling three-dimensional electromagnetically driven flow in inductively stirred systems, and its application to modeling electromagnetic and flow phenomena in sub-mold rotary stirring of billets. It is based on a new T-\jf potential formulation of electromagnetic field, which does not require the calculation of the field in free space. A finite element algorithm has been developed for solving the electromagnetic field and Navier-Stokes equations in the fluid flow domain. The computed results have shown that the electromagnetic force field generates a strong rotational flow within the vertical section covered by the stirrer, and a relatively strong secondary flow beyond the stirrer. It has also been shown that the rotational and secondary flows were driven primarily by the vorticity of the force field at the billet corners. The induced flow in the molten pool was found to be turbulent and the effective mixing region in the molten pool was about three times the length of the stirrer. The mathematical approach presented in this paper provides the theoretical framework for quantifying these effects for any given stirrer configuration and operating parameters.

Jan 1, 1999

By Shawn A. Cefalu, Matthew J. M. Krane, Kent J. VanEvery

"Segregation during electroslag remelting (ESR) of a nickel-chromium-molybdenum alloy (Ni-8 wt% Cr-25 wt% Mo) has been modeled. Features of the model include species transport and AC electromagnetic effects coupled with fluid flow, heat transfer, and solidification. The effect of varying ingot diameter and melt rate, on macrosegregation patterns, local solidification time, and liquid pool profiles obtained from this model are demonstrated. The results show that the macrosegregation in the ingots correlates to the filling velocity.IntroductionElectroslag remelting (ESR) is a secondary melting process used to refine the chemical composition and the microstructure of metals. Macrosegregation of solute elements causes defects in ESR ingots, including freckles and the maldistribution of phases. These defects are not correctable through heat treatment due to the large length scales over which they occur. However, the frequency and distribution of these defects may be minimized, or even prevented, through process control based on an understanding of the conditions that lead to their formation.In the past, studies have been conducted on the effect of the process parameters on the shape of the melt pool and the thickness of the mushy zone [1-7]. In these studies, the shape of the liquid pool is determined by consideration of the heat and fluid flow induced by electromagnetic (Lorentz) force and thermal buoyancy. The transport of solute and the resultant solutal buoyancy effects were not considered in these macroscopic models. Defect formation is tied to formation of secondary phases (e.g., freckling and white spots in nickel superalloys) [5,8]. In cases where secondary phase formation is calculated (e.g., [6]), phase fractions are determined using a micro-scale model given the local solidification time (LST) and local thermal history from the macroscale model.Using the local thermal history from the macroscopic transport model to predict microstructure increases the importance of the accuracy of those results. This investigation models the macrosegregation of molybdenum during the ESR of a Ni-Cr-Mo alloy. The effect of macrosegregation on pool shape and depth is demonstrated, along with its effect on the local solidification time. By comparing the relative effects of the melt rate and ingot radius, a better understanding of the factors that affect macrosegregation, and ultimately secondary phase formation may be obtained."

Jan 1, 2004

By S. Esfahani, M. Barati

Heat recovery from blast furnace slags is often contradicted by another requirement, to generate amorphous slag for its use in cement production. As both the rate and extent of heat recovery and slag structure are determined by its cooling rate, a relation between the crystallization kinetics and the cooling conditions is highly desired. In this study, slags with different basicities were studied during non-isothermal cooling. Their continuous cooling transformation (CCT) diagram were plotted and compared with each other. Then, a mathematical model was developed that predicts the amount of crystallization during cooling. The results were validated by experiments performed by Single Hot Thermocouple Technique and image analysis. A second mathematical model that calculates temperature history of slag during its cooling was coupled with the above model, to allow optimizing the slag/air ratio for both heat recovery and slag glass content.

Jan 1, 2015

By Leslie C. Taylor, Anthony L. Kee

Computational modeling of explosively loaded saturated sand on a suspended plate is an inherently complex problem. In this study a computational method is used to predict the dynamic pressure load/impulse from a buried charge. A fully coupled code modeling the explosive charge, sand and air with a Eulerian mesh and the plate structure with Lagrange elements, such as Gemini-Dyna, is an appropriate computational tool. In water-saturated sand, the loadings are soil and water-transmitted shock and bubble or blast effects depending on the depth and the charge size. The computations provide the capability to determine the total impulse load placed on a horizontal plate on or above the soil due to an explosion in the saturated sand for different burial depths in saturated sand and different plate stand off distances. The following test case will be presented: (a) 4.5 Kg (10 lbs) TNT, 10.16 cm (4 in) depth of burial to top of charge, and 20.32 cm (8 inch) stand off distance (Test 4).

Jan 1, 2005

By Kenneth P. Michalek, James E. Kelly, Jonathan A. Dantzig

"Many investigations have been reported regarding the simulation of continuous casting with emphasis on in-mold heat transfer. These analyses have been limited, however, because they require either estimation or measurement of the characteristic heat transfer between the ingot and the mold. In this paper, a self-consistent model 'is presented for calculating the heat transfer between a solidifying continuous casting and the mold, in the context of the continuous casting of steel. Using both thermal and stress analyses, the air gap formed during continuous casting is calculated and used to compute the heat transfer coefficients. Iteration between the thermal and stress analyses is used to determine a self-consistent solution. IntroductionThe continuous casting process for steel is illustrated schematically in Fig. 1. Molten steel is tapped into a tundish from which it flows at a controlled rate into the continuous casting mold. The initial shell of the casting forms in the mold by transfer of heat from the metal through the mold wall to the cooling water. The remainder of the casting is solidified by application of water sprays, and/or radiation cooling below the mold.Heat transfer in the mold portion of the process is important because it controls the formation of various defects which are deleterious to subsequent processing. Examples of defects formed in this region are surface irregularities, such as cold folding and oscillation marks, and more serious defects, such as longitudinal or transverse cracks. A more catastrophic problem which may occur is a breakout, in which the shell ruptures and molten steel is discharged in an uncontrolled way."

Jan 1, 1986

By Gordon Irons, Kevin Graham, Krishnakumar Krishnapisharody

Ladle or secondary metallurgy in steelmaking is an essential process step for the production of high-quality steel. Professor David Robertson was a pioneer in the modeling of these processes, in particular gas-metal fluid dynamic interactions and the rates of slag-metal reactions. This is now a highly active area of research, and some of the recent developments in this area will be reviewed.

Jan 1, 2014

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 K. P. Nishad

The origin of many casting defects can be traced at the mold filling stage. Entrapped air/gas pockets and oxide film formation are notable examples. In the present work, the mold-filling process is simulated using a mathematical model, built on a multi-physics based commercial package PHYSICA. The simulation code has the capability to solve multi-phase problems. The numerical model addresses a complicated dynamics of filling operation including free surface flow and subsequent flow behavior and is based on conservation of mass and momentum. The Free surface is modeled using the scalar equation method (SEM). The model is validated with benchmark experimental results from the literature. The role of turbulence and the pouring rate on variables affecting sites prone to oxide film formation are investigated using a simple but novel technique. Key findings of the work are presented.

Jan 1, 2006

By M. M. Botz

Cyanidation tailings disposed of in a surface impoundment experience a loss of cyanide due to natural attenuation, which frequently reduces the cyanide concentration to very low levels. Quantifying cyanide losses in terms of impoundment geometry, local weather conditions and feed-solution chemistry has been largely empirical in spite of the fact that, in many cases, mining operations rely on surface impoundments to reduce cyanide to below an internally regulated concentration or below an effluent limitation. To permit a quantitative evaluation of cyanide losses in an impoundment, a computer simulation was developed to estimate the losses of free, weak acid dissociable (WAD) and total cyanide due to dissociation, photolysis and volatilization. Results of the model are compared with data collected for a North American tailings impoundment in 1998.

Jan 1, 2000

By B. Farouk, D. Apelian

"A numerical model has been developed to predict the temperature history of particles injected in a low pressure d.c. plasma jet. The temperature and velocity fields of the plasma jet are predicted as a free jet by solving the parabolized compressible Navier-Stokes equations using a spatial marching scheme. Particle trojectories and heat transfer characteristics are calculated using the predicted plasma jet temperature and velocity fields. Correction factors have been introduced to take into account non-continuum effects encountered in the low pressure environment. The plasma jet profiles as well as the particle/plasma interactions under different jet pressure ratios (from underexpanded to overexpanded) have been investigated.I. IntroductionIn the plasma spraying process a d.c. plasma jet is used as a heat source to melt and accelerate powder particles which subsequently impinge and solidify on a given substrate. During low pressure plasma deposition (LPPD) prooess, the plasma jet is expanded into a low pressure environment (40-80 Corr) which results in high plasma jet velocities with Mach numbers ranging from 2 to 3. Previous efforts to predict heat transfer and fluid flow in plasmas have made use of incompressible flow formulation employing turbulence models [1-3] for atmospheric plasma spraying. However, these theoretical calculations have been limited to the subsonic flow regime; hence, the compressibility effects and viscous heat dissipation were neglected. In the LPPD process, the plasma jet is operated at the supersonic flow regime, therefore, the compressibility effects and viscous heat dissipation are no longer negligible. This is the main motivation of this work.In addition, production of high quality dense deposits requires that a large fraction of the injected particles be in a molten state when they impact the substrate. An additional motivation for this work is that it is important to understand the plasma/particle heat and momentum transfer since these phenomena dictate both the particle thermal history and particle trajectory.Early computational works in supersonic nozzle field flow [4,5,6] consisted of patching methods that divided the flow field into an invisoid free stream and a viscous boundary and. mixing layers. Each was analyzed independently and coupled through appropriate boundary conditions. Later schemes [7,8] involved solving time-dependent, compressible Navier-Stokes equations over the entire computational domain. These schemes show a great promise in predicting flows with complex structure, and where the viscous effects became prevalent."

Jan 1, 1986

By Murali Thirukkonda, Rahul Shah, Srinivas Tangrila, Shankar Rachakonda, Jan Kajuch, Arthur Gurson

"High-temperature superconductors have potential for a variety of applications such as power generators, superconducting magnets for mine sweepers or ship propulsion motors, and magnetic levitation transportation systems. Interest in the use of superconducting wires and tapes increased dramatically after the discovery of high-critical-temperature (high-Tc) oxide superconducting compounds, which offer operating temperatures above that of liquid nitrogen. Several techniques have been proposed to produce commercially viable high-Tc oxide based superconducting wires. One of the promising methods involves multi-pass powder-in-tube (PIT) wire drawing followed by rolling and heat treatment. This work focuses on the development of a numerical model to simulate the PIT drawing process for fabrication of long lengths of silver sheathed Bi-2212 superconducting wires. The finite element based model was used to predict densification of the oxide powder, wire drawing forces, and silver sheath thickness during drawing. A cap-type pressure dependent constitutive equation was implemented in the model to simulate the powder behavior. The model incorporated experimentally obtained material data for the silver and powder. Results indicate that powder densification is increased by decreasing the initial silver thickness and increasing the silver strength. Data from wire drawing experiments were used to verify model predictions. IntroductionSuperconducting wires and tapes can be used to carry large dc currents with no measurable loss due to resistivity. When used in superconducting magnets for power generators, motors, levitation devices, magnetic detectors, or energy storage devices, they could potentially result in significantly reduced energy losses. In the past, most prototypes were tested in liquid helium at 4.2K [1]. However, the necessity to provide liquid-helium cooling systems proved to be a serious limitation in commercialization. In 1986 high-critical-temperature (high-Tc) oxide superconductors were discovered [2], which could allow for higher operating temperatures above liquid nitrogen (77 K). Since then, both academia and industry have concentrated their efforts on the development of high-Tc oxide superconducting wires, tapes and coils. As yet, large scale fabrication of long lengths of these wires with the required current carrying capacity and mechanical properties remains an enormous challenge."

Jan 1, 1994

By E. Kozeschnik

A new model is presented for the simulation of the evolution of precipitates in a multi-component multi-particle multi-phase environment The mathematical formulation as well as the physical background of the model are briefly introduced and its capabilities are discussed on the example of a precipitation sequence in a complex steel. Keywords: Modeling, precipitation kinetics, nucleation, growth, coarsening

Jan 1, 2003

By Michael Klichowicz, Holger Lieberwirth

"For a detailed prognosis of mineral comminution processes, simulation tools like DEM still pose various challenges. Less realistic mineral structures can affect the results of simulation models such as bonded-particle models (BPM). Hence, in order to design mineral processes with these models as tools, representative mineral input data should be used. The article shows, how realistic 2D-structures can be modelled by taking advantage of the parameters identified by the Quantitative Mineralogical Analysis (QMA). Important quantitative parameters can be derived by analyzing three orthogonal thin or polished sections. For this reason point, line, and area analysis are applied. The mentioned parameters precisely describe the texture and structure of minerals. The newly developed software is able for the first time to semi-automatically analyze images of mineral structures using QMA. Moreover, statistically equivalent structures of minerals can be modelled using the QMA- parameters. Generating synthetic, but equivalent mineral structures opens the door for more reliable BPM based simulations. Within future research, the procedure will be adopted to generate 3D-structures. Therefore, QMA-based routines that already have been applied to generate 2D-models shall be utilized as basis. The aim of this procedure is to synthesize realistic microstructures based on certain statistical figures in order to improve comminution simulations.INTRODUCTION Due to its increasing capability, the use of computer technology is permanently intensifying in the design of mineral processes. The recent past has impressively shown that computer models can enhance the understanding of mineral processing. There is in fact a rising number of research activities, which are using simulation approaches to study various aspects of it. Weerasekara et al. (2013) give a comprehensive overview over the application of discrete element modelling (DEM) in the science of comminution.Nonetheless, using simulation tools such as DEM in general still poses a challenge. For a detailed simulation of a mineral comminution process with DEM, bonded particle models (BPM) appeared suitable in the past. The following issues, however, affect an appropriate prediction when modelling with BPM:• Simulations with high numbers of particles may need a lot of time for computation. • Internal modeling parameters like the spring and damper constants in DEM have to be adjusted by comparing its specific crushing behavior in the simulation to reality.• In order to use DEM for prediction of mineral processing, it is essential to validate the correctness of the models and assumptions, which is sometimes difficult. As Weerasekara et al. (2013) have stated it is crucial that the DEM method as well as the DEM code work properly.• It is challenging to implement the complex geometric structures of real minerals into the model."

Jan 1, 2016

By Anil Saigal, Jr. Fuller, Said Jahanmir

"All-ceramic crowns, including glass-infiltrated spinel composites, are coming into widespread use because of their superior aesthetics and chemical inertness. This study investigates the residual stresses that are developed in these composites as a result of cool down from the glass-infiltration temperature to room temperature due to slight mismatch in the coefficients of thermal expansion and its effect on the mechanical behavior of these composites. Two-dimensional finite element simulations were performed using an object oriented finite element program OOF. The OOF program is a combination of two programs. The first program, PPM200F, is designed to read an image file such as a micrograph. The individual pixels that constitute the micrograph may be collected into groups and their material properties assigned. PPM200F is then used to create the finite element model/mesh that OOF then reads. The average residual stresses are found to be tensile in the spinel matrix and compressive in the infiltrated-glass. There is large variation in residual stresses and strains from location to location with presence of locations at which the glass is under tensile stress. The crack initiation and initial propagation in glass-spinel composites is at the glass-spinel interface in both the glass and the spinel. The presence of residual stresses can lead to lower crack initiation stresses and degrade the mechanical properties of the composites.IntroductionGlass-infiltrated spinel ceramics are in clinical use as core materials in all-ceramic crowns and as inlays and onlays (1-3). Advantages of ceramic composites include near net-shape forming with high strength in the infiltrated state. Besides, spinel is quite translucent and may potentially be used without a veneer. However, these ceramic composites are subject to damage accumulation from repeat oral contact loading and often fail after a few years of use (4). This paper examines the residual stresses that are developed in these composites, after cooldown from the glass infiltration temperature to room temperature, due to the slight mismatch in the coefficients of the thermal expansion of the two phases."

Jan 1, 2001