Search Documents

By C. Hennesmann, D. Maijer, S. L. Cockcroft, P. Yo

"Die cast aluminum wheels are one of the most difficult automotive castings to produce because of stringent cast surface and internal quality requirements. As part of a collaborative research agreement between researchers at the University of British Columbia and Canadian Auto Parts Toyota Incorporated, work has been underway to predict heat transport and porosity formation in die cast A356 wheels. The basic heat transfer equations have been solved using the commercial finite element package ABAQUS to which specialized routines have been added to predict porosity. Preliminary work on predicting porosity formation focused on using the Niyama parameter as a measure of the probability of porosity. The latest version of the model incorporates a new fundamentally based porosity criterion, which takes into account the effect of hydrogen concentration. The development of this model together with its application to a series of cylindrical test castings as well as a production cast wheel are presented.I. IntroductionDie cast aluminum wheels are one of the most difficult automotive castings to produce because of stringent cast surface and internal quality requirements. Microporosity, in particular, is a common defect that can adversely affect wheel quality. The mechanism of microporosity formation in aluminum alloys is complex owing to its dependence on the interaction of a number of phenomena occurring during solidification including the decrease in hydrogen solubility, the decrease in volume, the development of the solidification structure, and the increase in difficulty in interdendritic feeding [1]. As a result, the amount, size, and distribution of microporosity is affected by a large number of casting parameters including thermal variables, melt composition, grain structure, modification, and inclusion content, making optimization by trial-and-error methodologies difficult.The use of criteria functions, based on thermal conditions during casting, offers a relatively simple approach to predicting microporosity. In practical terms, they are useful because they provide a straightforward method of assessing the impact of varying casting parameters on porosity through various thermal parametersl21 such as thermal gradient, cooling rate, solidification time, and solidus velocity. Drawbacks of criteria functions include insensitivity to varying hydrogen concentration and the limitation to qualitative predictions. Consequently, in aluminum alloys, attempts have been made to develop models incorporating the effects of hydrogen gas and microstructure evolution. However, these types of models represent a significant increase in complexity and sacrifice the simplicity and ease of application of the thermally based criteria functions."

Jan 1, 2001

By A. McLean, L. Gao, Z. Shi, Y. D. Yang, K. Chattopadhyay, G. F. Zhang

Levitation melting of metals is one of the emerging applications of magnetic and electric fields in liquid metal refining. The varying magnetic field generates induced current inside the metal droplet and leads to the Joule heating effect and Lorentz force against gravity. This contactless process can avoid contamination from the crucible and can be used to produce high-purity materials. Visualization of the levitation melting process is difficult, and thus generally investigated by mathematical modeling. In the present study, a three-dimensional mathematical model of metal levitation was developed using finite element method for proper understanding of the levitation process. The effects of harmonic magnetic field frequency, coils power input, and sample size on levitation nickel droplet were investigated. Temporal evolution of temperature fields were predicted directly though the model. Levitation and balance condition of the droplet were investigated as well.

Jan 1, 2015

By J. R. Bhamidipati, Nagy EI-Kaddah

"Over the years, mathematical modeling has been established as a viable tool for design and analysis of melting and solidification processes. Modeling of electromagnetic casting and melting systems involves three coupled problems -- determination of the containment field, the shape of the meniscus, and the temperature field -- which have to be solved simultaneously. One of the difficulties in numerical simulation of these systems is in defining and gridding the solution domain as the meniscus shape of the liquid region is not known a priori. This paper describes a methodology for solving these coupled problems in a two-dimensional electromagnetic containment field. This method is based on the solution of the electromagnetic, free surface dynamic, and heat transfer equations in a fixed domain in computational space by mapping the physical space into an invariant computational space during the search for the equilibrium shape using a variable non-orthogonal coordinate transformation. This approach, also, eliminates the need for regridding the liquid and solid regions during phase change. This methodology is demonstrated for two electromagnetic confinement systems -- the Magnetic Suspension Melting process and the electromagnetic casting of aluminum. The model predictions are in good agreement with available data for these systems. It is suggested that the methodology presented here is likely to provide a useful tool in the design, analysis, and optimization of electromagnetic confinement systems. IntroductionDuring the past decade, mathematical modeling has emerged as a viable tool for the design, analysis, and optimization of casting and melting processes /1-2/. It is generally recognized that the usefulness of modeling to study the interplay between the various process parameters without resorting to extensive experimentation relies heavily on accurate representation of the physical phenomena in these systems. In modeling electromagnetic confinement casting and melting systems, it is important that the model address all aspects of electromagnetic field interaction with the metal, if accurate predictions are to be obtained /3/. The electromagnetic field affects the heat transfer phenomena primarily through joule heating. In addition, the electromagnetic forces play a key role in heat redistribution in the molten pool via melt stirring, and, in confinement systems, support the molten metal pool."

Jan 1, 1994

By Olusegun J. Ilegbusi

"A two-fluid model is used to calculate flow and heat transfer characteristics of a plasma jet. This model accounts for both gradient-diffusion mixing and unmixing that may result from pressure-gradient-body-force imbalance in the system. The unmixedness is considered to be essentially a Rayleigh-Taylor kind instability. Separate transport equations are solved for the plasma and ambient gas velocities, temperatures and volume fractions. The two fluids allowed to exchange mass, momentum and energy through empirical interaction parameters. The empirical coefficients are first established by comparison of predictions with available experimental data for shear flows. The model is then applied to an Argon plasma jet ejecting into stagnant air. The predicted results show the significant build-up of unmixed air within the plasma gas, even relatively far downstream of the plasma torch. By adjusting the inlet condition, the model adequately reproduces the available experimental data on the temperature profile.IntroductionPlasma jets have important applications in materials processing, some of which include spray deposition, melting and refining, heat treatment and materials synthesis. In a typical reactor, the hot plasma stream entrains a surrounding gas and the resulting heat and mass transfer and the condition of operation of the torch, determines to a large extent, the performance of the unit. The entrainment may prevent uniform mixing and produce regions of unmixed hot/cold gases and non-uniform deposits as a result of insufficient melting of deposit in cold regions. This phenomenon may also affect chemical reaction rate in plasma systems for NOx reduction in exhaust gases.The study of plasma phenomena is often conveniently divided into three parts namely, the plasma torch, the plasma jet and analysis of the particles carried in the jet. The present work concerns processes occuring in the near-field of the plasma jet i.e. just downstream of the torch."

Jan 1, 1994

By Amjad Asad, Boyd Davis, Kinnor Chattopadhyay, Rüdiger Schwarze, Christoph Kratzsch, Donghui Li, Lei GAO

Sodium (Na) and Lithium (Li) are produced using molten salt electrolysis. The electrochemistry of the electrolyte is well-researched; however, there are benefits to understanding the melt flow and implications on it for cell design modifications. The basic configuration of alkali metal cells is the Downs cell. This consists of a central anode surrounded by a cathode, and this geometry was the basis for this mathematical modeling study. The behavior of gas bubbles in molten electrolyte was studied in both Na and Li cells through the use of computational fluid dynamics (CFD) techniques. The distance between the anode and the cathode was varied in the CFD model to ascertain whether strong circulatory flows would change significantly in the cell. The standard k-ε turbulence model was used to mimic turbulent flow, and a two-way coupled Discrete Phase Model (DPM) was adopted to simulate flotation behavior of chlorine bubbles and liquid metal droplets. The liquid metal distribution on the free surface was predicted using the Volume of Fluid (VOF) multi-phase model.

Mar 1, 2017

By Han-Tao Hu

Keywords: RH refining process, flow of molten steel, two-phase flow, gas holdup, circulation rate, mathematical modeling A three-dimensional mathematical model for the flow of the molten steel in the whole degasser during the RH refining process has been proposed and developed with considering the physical characteristics of the process, particularly the behaviors of gas-liquid two-phase flow in the up-snorkel and the momentum exchange between the two phases. The fluid flow fields and the gas holdups of liquid phases and others respectively in a 90 t RH degasser and its water model unit with a 1/5 linear scale have been computed using this model. The results showed that the flow pattern of molten steel in the whole RH degasser could be well modeled by the model. The liquid can be fully mixed during the refining process except the area close to the free surface of liquid and zone between the two snorkels in the ladle, but there is a boundary layer between the descending liquid stream from the down-snorkel and its surrounding liquid, which is a typical liquid-liquid two phase flow, and the molten steel in the ladle is not in a perfect mixing state. The lifting gas blown is rising mostly near the up-snorkel wall, which is more obvious under the conditions of a practical RH degasser, and the flow pattern of the bubbles and liquid in the up-snorkel is closer to an annular flow. The calculated circulation rates for the model unit at different lifting gas rates are in good agreement with the determined values.

Jan 1, 2005

By Mike Trovant, Stavros A. Argyropoulos

"A fixed grid numerical model for phase change problems is proposed. The model solves for heat transfer and fluid flow, via the coupled solution of the energy, continuity and momentum equations. In addition, the numerical scheme accounts for temperature varying thermophysical properties and determines the shrinkage profile resulting from phase and density change. Simulations are run to determine the effect of varying model parameters and wall cooling rates on the location and shape of the shrinkage void formation in cylindrical casts. Results are shown for both aluminum and magnesium.IntroductionThe mathematical modeling of phase change processes (be they solidification or melting problems) entails solving for many interrelated phenomena. A comprehensive numerical simulation of the casting of a metal, for example, would have to account for: (a) heat transfer in both the solid and liquid portions of the metal and in the mold, (b) the release of latent heat during phase change, (c) fluid flow behavior in the liquid portion of the metal, (d) temperature induced shrinkage caused by liquid contraction and solidification contraction, and (e) the development of thermal stresses with accompanying shrinkage in the solid portion of the cast and mold. A schematic of the relationships that exist between the coupled equations which govern solidification is shown in Figure 1."

Jan 1, 1994

By F. Chabchoub

A mathematical model has been developed to study the influence of process variables on the horizontal Ohno continuous casting process. The model simulates the heat transfer and fluid flow in a wide plate ingot of pure tin. A two dimensional numerical analysis based on a control-volume finite difference approach, the Semi-Implicit Method for Pressure-Linked Equations algorithm (SIMPLER), and an Enthalpy Transforming Model was used. The shape and position of the solid-liquid interface were calculated at different casting speeds, ingot thicknesses, and cooling device characteristics, as well as various mould temperatures. Results of the mathematical model were compared with

Jan 1, 1992

By B. Abedian, L. M. Racz, S. R. Berry, R. W. Hyers

"The presented work represents an effort to improve understanding and prediction of turbulent flow inside electromagnetically-levitated droplets. We show that the flow field in a test case, a nickel droplet levitated under microgravity conditions, is a low Reynolds number turbulent flow, i.e. in the transitional regime between laminar and turbulent. Past research efforts have used laminar, enhanced viscosity, and k - E turbulence models to describe these flows. In our study, we compare the RNG algorithm to laminar and k - E turbulence model results. We show that accurate description of the turbulent eddy viscosity µr is critical in order to obtain realistic velocity fields, and that µr cannot be uniform in levitated droplets. In the RNG method there are no characteristic length or time scales associated with the flow, thus allowing such anisotropic features to be captured. We perform calculations and analyses for two cases (1) a small, non-deforming, spherical nickel droplet and (2) a deforming, oscillating droplet. In addition to flow field calculations, we compare numerically-obtained surface tension and viscosity values to those obtained experimentally by the oscillating drop technique.IntroductionElectromagnetic levitation (EML) is a containerless processing technique for conductive materials. Its advantages are the avoidance of contamination or chemical reaction at container walls and the ability to perform long-term experiments in a metastable state. Because of its current limitation to small sample sizes, its primary use is in the laboratory. One application of interest is the measurement of thermophysical properties of undercooled liquid metals [1, 2, 3, 4]. In EML, a sample weighing approximately 1 gram is positioned inside an electromagnetic field generated by alternating currents flowing through an axisymmetric water-cooled copper coil. The coil must be designed appropriately so that thesample levitates in a stable manner.Traditionally, mathematical modeling has been performed in tandem with both earth-and spacebound electromagnetic levitation experiments, in order to determine the effect of the magnetic field on material and flow characteristics [5, 6, 7, 8]. These calculations are important because they allow design, planning, and correct interpretation of experiments. In addition, the calculations are of fundamental scientific interest, as they can be applied to other electromagnetic materials processing applications [9]."

Jan 1, 1999

By Elliot Schwartz, Julian Szekely

"In this paper we describe an experiment that was recently performed on the IML-2 Space Shuttle mission in which the thermophysical properties of molten metals were measured using electromagnetic levitation. In the planning of these experiments calculations had to be performed to predict the levitation and stirring forces, heating rates, sample deformation, and internal flows within the sample. These predictive calculations were essential for both the planning of the experiments and for the re-planning of experimental runs during the actual Space Shuttle mission.Introduction""Materials Processing in the Computer Age"" involves the use of mathematical modeling to obtain a quantitative understanding of the phenomena that govern materials processing operations. Microgravity, or space shuttle experimentation represents a specific case in point here; the very high cost of the space shuttle flights, the very limited flight opportunities that are available, and the tight resources aboard a shuttle clearly mandate that these experiments be very precisely planned and that all possible opportunities for mathematical modeling be fully exploited.The purpose of this paper is to present the results of calculations performed for the planning of a space shuttle experiment involving electromagnetic levitation and wherever possible present a comparison of the computational predictions with both experimental measurements and asymptotic analyses. Finally we shall seek to discuss the practical implications of both the computational methodologies and the potential experimental findings."

Jan 1, 1994

By Y. Y. Sheng

Plume flows created by gas injection into liquid metals have been studied extensively, however the dynamics of bubbles within the plume have not received adequate attention, particularly for large, irregular bubbles, characteristic of the processing of liquid metals. The present model uses a combined Eulerian-Lagrangian approach in which an estimated void fraction distribution is used to calculate the liquid flow field with a conventional Eulerian scheme. The trajectories of bubbles are then computed in a Lagrangian manner using the estimated flow field, and experimentally measured information on bubble drag coefficients', lateral migration due to lateral lift forces and interaction with turbulent eddies, and variation in bubble size due to break-up. Turbulence in the two-phase zone is modeled with a modified k-e model with extra source terms to account for the second phase which have been recently determined by the present authors. The computed void fraction and turbulent liquid flow field distributions are in good agreement with experimental measurements.

Jan 1, 1993

By J. Zou

With the increasing applications of the Chinese metal magnetic liners, the requirement for their technical specifications has become stricter and stricter. It is necessary to deeply analyse the metal magnetic liners in order to improve their reliability and economics. The key issue of metal magnetic liners is magnetic forces: their attractions to mill shell and protection layer. Mathematical models were developed using finite-element analysis software and systematic simulation analyses of the two most critical indexes of magnetic liners were performed to find out the most important parameters and their effect on the liner?s characteristics. The relationships between magnetic forces and installation gaps, thickness of permanent magnets and thickness of protection layer may be used to guide the design of metal magnetic liners for ball mills in order to achieve good economic indexes.

Jan 1, 2014

By M. I. Kalashnikova, Y. M. Shneerson, E. E. Zhmarin, Y. V. Degtiarev, E. M. Vigdorchik

Separation of a magnetic fraction from nickel concentrate following matte flotation, and its treatment in a separate circuit enable reduces precious metals losses. The magnetic fraction is treated in two stages of oxidative leaching: atmospheric and pressure leaching. A study of the atmospheric stage, using mathematical modelling, is described in this paper. The mathematical model and experimental data were used to carry out industrial process calculations. The model enabled extrapolation of the results of batch laboratory test to industrial-scale continuous processing and the selection of optimal process conditions and large-scale equipment. In December 2004, commercial oxidative leaching tests were performed on the magnetic fraction at the Severonickel Plant (Monchegorsk, Russia). A comparative analysis of the mathematical modelling results and the actual process parameters is presented.

Jan 1, 2005

By Vladimir Panjkovic, Ostrovski. Oieg, Campbell Cripps Clark, John Truelove

"Effects of properties of raw materials and blast on the performance of an iron bath reactor have been examined by a mathematical model. The model consists of three modules: 1) chemical reactions in the gas phase and heat transfer from gas to droplets ejected from the bath; 2) heat losses from the vessel and radiative heat transfer from the freeboard to the bath; and 3) heat and mass balances for the bath. The following parameters of the process are calculated: coal, ore, flux and oxygen consumption, post-combustion ratio and heat-transfer efficiency, temperature and composition of the off-gas, and slag and metal chemistry.IntroductionOne of the anticipated advantages of smelting reduction over conventional ironmaking is the use of lower-grade raw materials. In this work, the effect of-properties of raw materials on the performance of an iron-bath reactor is examined with a numerical model. The model is also used to assess the effect of properties of the blast (flowrate, temperature and oxygen enrichment) and ejected droplets (generation rate, size and ratio of total weights of slag and iron droplets) on the key parameters of the process. In the process coal, oxygen, ore (which may be partly prereduced) and fluxes for specified slag composition are charged to the bath. Slag and iron droplets are ejected from the bath into the gas phase, where they are heated before returning to the bath. Iron droplets undergo decarburisation and reoxidation with CO2 and H2O. Dust is also generated from the bath, and oxygen or preheated air for post combustion is injected from the top of the vessel.The model consists of three modules: 1) reactions in the gas phase and heat transfer from gas to droplets ejected from the bath; 2) heat losses from the vessel and radiative heat transfer from the freeboard to the bath; and 3) heat and mass balances for the bath. The structure of the model and the data exchange between the modules are shown in Fig. 1. The following parameters of the process are calculated: coal, ore, flux and oxygen consumption, postcombustion ratio (PCR), heat transfer efficiency (HTE), off-gas temperature and composition, slag and metal chemistry, and rates of decarburisation and reoxidation of the iron droplets."

Jan 1, 1996

By Shoji Taniguchi, Lifeng Zhang, Kaike Cai

"In the present paper, based on the two-phase (Eulerian-Eulerian) model, the 3D fluid flows in gas-stirred systems, i.e., air-stirred water vessel and argon-stirred liquid steel ladle, are simulated. In the Eulerian-Eulerian two-phase model, gas and liquid are considered to be two different continuous fields. The phases are assumed to share space in proportion to their volume fractions. The exchange between the phases is represented by source terms in conversation equations. Turbulence simulated by the k- e model is assumed to be the property of the liquid phase. The following effects of mathematical treatments and operational factors on the fluid flow are discussed, such as interphase drag force, turbulent model, the size of bubble, and gas injection mode.Some interesting results are derived. Except the interphase drag force, the interphase lift force should be taken into accounted in order to exactly simulate the fluid flow in gas-stirred systems. Suitable turbulent model and drag coefficient have effect on the mathematical simulation. Injecting small bubbles can realize a well mixing flow condition. The distance between the two gas-injection nozzles has effect on the fluid flow."

Jan 1, 2000

By Saul Garcia Hernandez

The research work objective is to study a new ladle shroud internal design to reduce the slag emulsification in the tundish during ladle transient periods. A 3D mathematical model of the tundish was created to solve the fundamental equations, a turbulence model, and a multiphase model. The results show that using a conventional ladle shroud, the fluid is delivered with an excessive amount of kinetic energy being dissipated inside the tundish bath, generating strong mixing flow patterns and entrapping a massive amount of slag. In the other hand, the proposed ladle shroud dissipates the kinetic energy before the fluid enters the tundish promoting less intense mixing patterns; thus, the amount of slag emulsification is reduced significantly. Consequently, if the internal design of the ladle shroud is able to dissipate the kinetic energy before the fluid enters the tundish, it shall be possible to reduce considerable the slag emulsification and the slag opening area.

Aug 17, 2017

By A S. Kurilko

For mines located in the cryolitic zone, thermal conditions are an important factor for all kinds of mining works, especially for ventilation of workings. In the present paper, some models of mine air heat transfer considering air distribution in a net of mine workings are developed. A combined solution is suggested for heat transfer and air distribution problems in mine workings and other underground structures located in the cryolitic zone. A software package is also presented for calculating temperature conditions in an underground refrigeration plant and for optimisation of ventilation conditions. This software uses a special 3D mathematical model to describe heat transfer in a freezing chamber between surrounding rock massif and stored product as well as between atmospheric air and the daylight surface of rock massif. A problem requiring strengthening of a backfill material and its thermal interaction with surrounding rocks during underground deposit mining by horizontal layers in conditions of the cryolitic zone is also discussed. The developed program takes into account cement consumption, geometrical dimensions of workings filled and a sequence of deposit mining.

Jan 1, 2005

By H. G. Wang

The production process of a coal mine can often be interrupted by all kinds of events scheduled or unscheduled and thus often needs gas drainage for a temporary or long-term sealed area. Gas drainage has a relatively high risk, therefore attention must be paid to its design to keep the process safe. Based on the conservation laws of mass and constituents, and under the conditions of the allowable gas concentration limit in the return air route, the relationship between gas concentration and gas drainage flow rate as a function of time during the drainage process is discussed in detail to build a mathematical model of gas drainage from a sealed area. With respect to the gas drainage process, two different conditions need to be considered: one with a methane free fresh air supply, the other with a low gas concentration. This paper presents the derived mathematical models for these situations and then shows that they can be incorporated into one model with some modification. This combined model can be employed to determine the optimum gas drainage flow rate with the shortest drainage time. The application is coded using the C# programming language within the .NET framework. Finally, the application program developed by the author is used to validate the mathematical model through an example of gas drainage with varying rate.

Aug 1, 2013

By A. E. Hilder

Introductory The following pages present the results of field investigations, laboratory and commercial tests, drilling, surveying, and preliminary development of the refractory clays which have been discovered recently in the township of Kipling, at the foot of Long rapids on the west side of the Mattagami river about 57 miles north of Kapuskasing and approximately 7 miles north of Smoky Falls, the terminus of the privately owned railway operated by the Spruce Falls Power and Paper Company, Kapuskasing (see Figure 1). Investigations on the clays of northern Ontario have been carried on over a period of seven years by officers of the Ontario Department of Mines, and I am indebted to the information published by that Department in reports by W. S. Dyer and A. R. Crozier, as well as?by Professor R. J. Montgomery, of the Ceramic Department, University of Toronto. It was, in large measure, the data given in the 37th, 38th, and 42nd Annual Reports of the Ontario Department of Mines that led me to make a personal investigation of the Mattagami River clays, with the result that we now have a commercial undertaking that will develop eventually into an industry of considerable importance in supplying the central manufacturing districts of Canada with refractories that heretofore have been largely imported from the United States. An extract from the 42nd Annual Report of the Department, Part III, 1933, page 75, reads as follows: "As for the refractory clay, the deposits are undoubtedly very large, but what is urgently needed is a deposit near the railway under sufficiently shallow overburden for recovery by stripping'. Having previously made an investigation of the clays adjacent to the Temiskaming and Northern Ontario railway, and these clays proving unsuitable, I turned my attention to the outcrops of refractory clays on the Mattagami river, first reported by Keele in 1919. This area offered possibilities within the economically permissible distance from rail transportation.

Jan 1, 1935

Recovery of tin and copper from concen trates containing cassiterite, stannite and chalcopyrite was studied in laboratory-scale tests designed to simulate continuous operation of the matte fuming process in which tin is recovered as SnS fume while copper accumulates in the matte. Charges made up of 300g of mixtures of the concentrate and siliceous flux were added at a uniform rate in pelletized form to crucibles containing 300g of matte and 50g of slag (both of which were recovered from previous tests), while blowing through the matte with air or oxygen-nitrogen mixtures at a rate of 3 to 5 1/min. Accurate material balances enabled the effects of process operating conditions on tin recovery and sulfur utilization to be deter- mined. When operated correctly, the matte fuming process gave tin recoveries exceeding 90 percent in spite of wide variations in the copper content of the matte and in the operat- ing conditions. Deterioration of tin recovery was associated with undersaturation of the slag with respect to silica, and with high partial pressure of sulfur dioxide (greater than approx- imately 0.2) brought about by oxygen enrichment of the blowing air. The utilization of sulfur varied widely in low grade mattes. Copper in excess of approx- imately 10% promoted formation of SO2 by the commonly accepted reaction: FeS + 1-120 2= FeO + S02 while in copper-free and lower grade mattes a substantial fraction of the oxygen was consumed in the alternative reaction FeS + i 02 = FeO + z S2.

Jan 1, 1980