Obtaining data to validate a model of an induction skull melting furnace

A casting route is often the most cost-effective means of producing engineering components. However, certain materials, particularly those based on Ti, TiAl and Zr alloy systems, are very reactive in the molten condition and must be melted in special furnaces. Induction Skull Melting (ISM) is the most widely-used process for melting these alloys prior to casting components such as turbine blades, engine valves, turbocharger rotors and medical prostheses. A major research project is underway with the specific target of developing robust techniques for casting TiAl components. The aims include increasing the superheat in the molten metal to allow thin section components to be cast, improving the quality of the cast components and increasing the energy efficiency of the process. As part of this, the University of Greenwich (UK) is developing a computer model of the ISM process in close collaboration with the University of Birmingham (UK) where extensive melting trials are being undertaken. This paper describes the experimental measurements to obtain data to feed into and to validate the model. These include measurements of the true RMS current applied to the induction coil, the heat transfer from the molten metal to the crucible cooling water, and the shape of the column of semi-levitated molten metal. Data are presented for Al, Ni and TiAl.

Numerical simulation of vacuum dezincing of lead alloy

Removing zinc by distillation can leave the lead bullion virtually free of zinc and also produces pure zinc crystals. Batch distillation is considered in a hemispherical kettle with water-cooled lid, under high vacuum (50 Pa or less). Sufficient zinc concentration at the evaporating surface is achieved by means of a mechanical stirrer. The numerical model is based on the multiphysics simulation package PHYSICA. The fluid flow module of the code is used to simulate the action of the stirring impeller and to determine the temperature and concentration fields throughout the liquid volume including the evaporating surface. The rate of zinc evaporation and condensation is then modelled using Langmuir’s equations. Diffusion of the zinc vapour through the residual air in the vacuum gap is also taken into account. Computed results show that the mixing is sufficient and the rate-limiting step of the process is the surface evaporation driven by the difference of the equilibrium vapour pressure and the actual partial pressure of zinc vapour. However, at higher zinc concentrations, the heat transfer through the growing zinc crystal crust towards the cold steel lid may become the limiting factor because the crystallization front may reach the melting point. The computational model can be very useful in optimising the process within its safe limits.

Macro-segregation of multi-component alloys in shape casting simulation

The objective of this work is to present a new scheme for temperature-solute coupling in a solidification model, where the temperature and concentration fields simultaneously satisfy the macro-scale transport equations and, in the mushy region, meet the constraints imposed by the thermodynamics and the local scale processes. A step-by-step explanation of the macrosegregation algorithm, implemented in the finite volume unstructured mesh multi-physics modelling code PHYSICA, is initially presented and then the proposed scheme is validated against experimental results obtained by Krane for binary and a ternary alloys.

A comprehensive copper stockpile leach model: background and model formulation

The first phase in the sign, development and implementation of a comprehensive computational model of a copper stockpile leach process is presented. The model accounts for transport phenomena through the stockpile, reaction kinetics for the important mineral species, oxgen and bacterial effects on the leach reactions, plus heat, energy and acid balances for the overall leach process. The paper describes the formulation of the leach process model and its implementation in PHYSICA+, a computational fluid dynamic (CFD) software environment. The model draws on a number of phenomena to represent the competing physical and chemical features active in the process model. The phenomena are essentially represented by a three-phased (solid liquid gas) multi-component transport system; novel algorithms and procedures are required to solve the model equations, including a methodology for dealing with multiple chemical species with different reaction rates in ore represented by multiple particle size fractions. Some initial validation results and application simulations are shown to illustrate the potential of the model.

The multiphysics modeling of solidification and melting processes

Solidification and melting processes involve a range of physical phenomena and their interactions (i.e., multiphysics). Computational modeling of such processes presents a significant challenge, both in representing the physics involved and in handling the resulting coupled behavior. Two methods for the computational modeling of multiphysics processes in complex geometries are highlighted in the context of four challenging applications

Computational modelling of freeze layers in smelting processes

The use of computational modelling in examining process engineering issues is very powerful. It has been used in the development of the HIsmelt process from its concept. It is desirable to further water-cool the HIsmelt vessel to reduce downtime for replacing refractory. Water-cooled elements close to a metal bath run the risk of failure. This generally occurs when a process perturbation causes the freeze and refractory layers to come away from the water-cooled element, which is then exposed to liquid metal. The element fails as they are unable to remove all the heat. Modelling of the water-cooled element involves modelling the heat transfer, fluid flow, stress and solidification for a localised section of the reaction vessel. The complex interaction between the liquid slag and the refractory applied to the outside of thewater-cooled element is also being examined to model the wear of this layer. The model is being constructed in Physica, a CFD code developed at the University of Greenwich. Modelling of this system has commenced with modelling solidification test cases. These test cases have been used to validate the CFD code’s capability to model the solidification in this system. A model to track the penetration of slag into refractory has also been developed and tested.

Wetting and Reaction of Sn-2.8Ag-0.5Cu-1.0Bi Solder with Cu and Ni Substrates

The wettability of newly developed Sn-2.8Ag-0.5Cu-1.0Bi lead-free solder on Cu and Ni substrates was assessed through the wetting balance tests. The wettability assessment parameters such as contact angle (ϑc) and maximum wetting force (Fw) were documented for three solder bath temperatures with three commercial fluxes, namely, no-clean (NC), nonactivated (R), and water-soluble organic acid flux (WS). It was found that the lead-free Sn-2.8Ag-0.5Cu-1.0Bi solder exhibited less wetting force, i.e., poorer wettability, than the conventional Sn-37Pb solder for all flux types and solder bath temperatures. The wettability of Sn-2.8Ag-0.5Cu-1.0Bi lead-free solder on Cu substrate was much higher than that on Ni substrate. Nonwetting for Sn-2.8Ag-0.5Cu-1.0Bi and Sn-Pb solders on Ni substrate occurred when R-type flux was used. A model was built and simulations were performed for the wetting balance test. The simulation results were found very close to the experimental results. It was also observed that larger values of immersion depth resulted in a decrease of the wetting force and corresponding meniscus height, whereas the increase in substrate perimeter enhanced the wettability. The wetting reactions between the solder and Cu/Ni substrates were also investigated, and it was found that Cu atoms diffused into the solder through the intermetallic compounds (IMCs) much faster than did the Ni atoms. Rapid formation of IMCs inhibited the wettability of Sn-2.8Ag-0.5Cu-1.0Bi solder compared to the Sn-Pb solder.

Busbar sizing modeling tools: comparing an ANSYS® based 3D model with the versatile 1D model part of MHD-Valdis

The main goal of a cell stability MHD model like MHD-Valdis is to help locate the busbars around the cell in a way which leads to the generation of a magnetic field inside the cell that itself leads to a stable cell operation. Yet as far as the cell stability is concerned, the uniformity of the current density in the metal pad is also extremely important and can only be achieved with a correct busbar network sizing. This work compares the usage of a detailed ANSYS based 3D thermo-electric model with the one of the versatile 1D part of MHD-Valdis to help design a well balanced busbar network.

Modelling the tilt-casting process for the tranquil filling of titanium alloy turbine blades

The tilt-casting method is used to achieve tranquil filling of gamma-TiAl turbine blades. The reactive alloy is melted in a cold crucible using an induction coil and then the complete crucible-mould- running system assembly is rotated through 180degrees to transfer the metal into the mould. The induction current is ramped down gradually as the rotation starts and the mould is preheated to maintain superheat. The liquid metal then enters the mould and the gas within it (argon) escapes through the inlet aperture and through auxiliary vents. Solidification starts as soon the metal enters the mould and it is important to account for this effect to predict and prevent misruns. The rotation rate has to be controlled carefully to allow sufficient time for gas evacuation, but at the same time preserve superheat. This 3-phase system is modelled using the FV method, with a fast implicit numerical scheme used to capture the transient liquid free surface. The enthalpy method is used to model solidification and predict defects such as trapped bubbles, macro-porosity or surface connected porosity. Modeling is used to support an experimental program for the development of a production method for gamma-TiAl blades, with a target length of 40cm. The experiments provide validation for the model and the model in turn optimizes the tilt-casting process. The work is part of the EU project IMPRESS.