Search Documents

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

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 L. Katgerman, E. N. Straatsma, Suyitno

"Single-roll strip casting process is studied using computational fluid dynamics which takes into account fluid flow, heat transfer and solidification. The computational model is based on the conservation Jaws of momentum, energy and continuity. Correlation of strip thickness, solidification zone length, and puddle shape on single-roll strip casting of Al-1 %Mn is explored. Feeding pressure and roll speed are investigated. The results show that the··strip thickness slightly increases with pressure, similarly the length of solidification zone. The variations of strip thickness and solidification zone length with roll speed show a sharp decline for strip thickness and slight decline for the length of solidification zone with increasing roll speed. The constraint and free-jet flow have different feature in term of flow field. Severe turbulence is observed on the free-jet flow.IntroductionThe near-net shape casting of strip metals is a demand nowadays. The demand is particularly stimulated by the challenge to become more competitive through enhanced performance. The improvement in the existing processes and new innovative technologies are introduced by the metal industry to gain better product quality, more efficiency in energy consumption, low throughput, yield and operating costs [1-3].There are many advantages that can be obtained by the direct casting of strip. The most prominence is the elimination of hot rolling. Beside that the high rate of heat extraction associated with strip casting offers the possibility of a finer, Jess segregated cast structure. The microstructure will result in the improvement of mechanical properties compared to conventional material, or alternatively the opportunity of using grades of metal outside the normal range for continuous casting.The concept of single-roll strip casting has been difficult to commercialize, because of the complex interactions between fluid flow, solidification, shrinkage, and stress, which lead to serious quality problems. Besides that, the quality and the thickness of the strip are affected by the thermal properties of the strip, roll cooling conditions, roll geometry, and other process parameters. To understand how these parameters affect the strip formation, some researchers have modeled the process [4]."

Jan 1, 2004

By Kwan Ho Kim

In this study, free fall tests are conducted with the objective of understanding the breakage properties of waste concrete for the guide of the direction of the better comminution process. Waste concrete of 75mm×53mm is tested at various drop heights with or without heat pretreatment. Test results indicate that a block of waste concrete does not break by free fall, but eventually breaks by repeated free falls. The number of free falls required for complete separation of large aggregates from the cement mortar considerably reduced when the samples were heated at 400?. Therefore, a great energy saving can be obtained when heating the waste concrete prior to crushing. The energy saving was found to be dependent on the height of free fall but in all cases, was still large enough to offset the additional energy required for heating. Also, it was found that the heat pre-treatment gives an additional advantage in that it produced cleaner recycled aggregates. This may result from easier detachment of cement mortar from the aggregates due to weakening of bonding forces by heating. By analyzing the size distribution of the crushed product, a free-fall breakage model is developed using the particle breakage ratio and the breakage pattern. The developed model can be used to predict and to analyze waste concrete processing circuit.

Jan 1, 2014

By J. L. Hill, M. A. R. Sharif, J. Lin

"A numerical algorithm has been developed to simulate the movement, breakup, and entrapment of the oxide film inclusions encountered during the mold-filling processes of aluminum castings. The flowfield is solved using the well-known Marker and Cell (MAC) method by a time marching process. The Volume of Fluid (VOF) method is used to track the free surface boundaries. A kinematic approach is employed to track the movement and breakup of the oxide films on the free surface or in the bulk liquid metal. The computer program based on the proposed algorithm is able to model the flow behavior and oxide film movement, breakup, and entrapment in mold cavities for aluminum castings.IntroductionMetallographic and fractographic studies of the aluminum castings show tangled networks of oxide films in the castings. These oxide films can form in the furnaces or during the mold filling processes. Green and Campbell [I] concluded that the oxide film defects dominate the strength of aluminum castings. The main cause of oxide film entrapment into the bulk molten metal during the casting processes is the surface turbulence. In film-forming alloys, if the surface breaks, due to a drop, a fountain, or a breaking wave, then the surface film is folded over and incorporated into the bulk liquid. The incorporation of surface films leads to cracks in the casting and other defects that are nucleated by the presence of the folded films, such as porosity and hot tears. Campbell [2] estimated that up to 80% of scrap in light-alloy casting is the result of surface turbulence."

Jan 1, 1999

By Nickolas J. Themelis

Nearly thirty two million tons of municipal solid wastes (MSW) are combusted annually in over one hundred U.S. Waste-to-Energy (WTE) power plants, thus reducing the use of fuel oil by about 1.6 billion gallons of fuel oil One of the advanced WTE processes is the 0.9-million ton/y SEMASS facility at Rochester, Mass. Most of the combustion in the three giant combustion chambers occurs while the injected shredded wastes are in suspension, i.e., in flash combustion mode. The velocity, temperature, and concentration profiles in the SEMASS combustion chamber were simulated by representing the combustible fraction of MSW by the simplified formula C6HIO04 and using the CFD program Fluent to solve the turbulent energy and mass transport equations. After validation against plant data, this model will be used to examine options for increasing the productivity of WTE combustion chambers, such as optimum distribution of primary and secondary air and oxygen enrichment.

Jan 1, 2003

By Nickolas J. Themelis

Nearly thirty two million tons of municipal solid wastes (MSW) are combusted annually in over one hundred U.S. Waste-to-Energy (WTE) power plants, thus reducing the use of fuel oil by about 1.6 billion gallons of fuel oil. One of the advanced WTE processes is the 0.9-million ton/y SEMASS facility at Rochester. Mass. Most of the combustion in the three giant combustion chambers occurs while the injected shredded wastes are in suspension, i.e .? in flash combustion mode. The velocity. temperature, and concentration profiles in the SEMASS combustion chamber were simulated by representing the combustible fraction of MSW by the simplified formula C6HIO04 and using the CFD program Fluent to solve the turbulent energy and mass transport equations. After validation against plant data. this model will be used to examine options for increasing the productivity of WTE combustion chambers, such as optimum distribution of primary and secondary air and oxygen enrichment.

Jan 1, 2003