Modeling the dynamics of Magnetic Semilevitation Melting

In semilevitation melting, a cylindrical metal ingot is melted by a coaxial a.c. induction coil. A watercooled solid base supports the ingot, while the top and side free surface is confined by the magnetic forces as the melting front progresses. The dynamic interplay between gravity, hydrodynamic stress, and the Lorentz force in the fluid determines the instantaneous free surface shape. The coupled nonstationary equations for turbulent flow, heat with phase change, and high-frequency electromagnetic field are solved numerically for the axisymmetric time-dependent domain by a continuous mesh transformation, using a pseudospectral method. Results are obtained for the two actually existing coil configurations and several validation cases.

Bulk synchronous parallelization of industrial electromagentic software

The parallelization of existing/industrial electromagnetic software using the bulk synchronous parallel (BSP) computation model is presented. The software employs the finite element method with a preconditioned conjugate gradient-type solution for the resulting linear systems of equations. A geometric mesh-partitioning approach is applied within the BSP framework for the assembly and solution phases of the finite element computation. This is combined with a nongeometric, data-driven parallel quadrature procedure for the evaluation of right-hand-side terms in applications involving coil fields. A similar parallel decomposition is applied to the parallel calculation of electron beam trajectories required for the design of tube devices. The BSP parallelization approach adopted is fully portable, conceptually simple, and cost-effective, and it can be applied to a wide range of finite element applications not necessarily related to electromagnetics.

Magnetic levitation fluid dynamics

This work is concerned with the accurate computation of flow in a rapidly deforming liquid metal droplet, suspended in an AC magnetic field. Intense flow motion due to the induced electromagnetic force distorts dynamically the droplet envelope, which is initially spherical. The relative positional change between the liquid metal surface and the surrounding coil means that fluid flow and magnetic field computations need to be closely coupled. A spectral technique is used to solve this problem, which is assumed axisymmetric. The computed results are compared against a physical experiment and "ideal sphere" analytic solutions. A comparison between the "magnetic pressure" approximation and the full electromagnetic force solutions, shows fundamental differences; the full electromagnetic force solution is necessary for accurate results in most practical applications of this technique. The physical reason for the fundamental discrepancy is the difference in the electromagnetic force representation: only the gradient part of the full force is accounted for in the "magnetic pressure" approximation. Figs 9, Refs 13.

Dynamics of magnetically suspended fluid

Magnetic suspension is a technique for processing pure or reactive materials without contact to walls. This work is concerned with the flow in the rapidly deforming liquid volume, suspended in an AC magnetic field. Intense flow motion due to the induced electromagnetic force distorts dynamically the droplet envelope. The relative positional change between the liquid surface and the surrounding coil means that fluid flow and magnetic field computations need to be closely coupled. The computed results are compared against a physical experiment and nearly spherical analytic solutions. A comparison between the "magetic pressure" approximation and the full electromagnetic force solutions shows fundamental differences; the full electromagnetic force is necessary for accurate results in most practical applications of this technique. The physical reason for the fundamental discrepancy is the difference in the electromagnetic force representation: only the gradient part of the full force is accounted for in the "magnetic pressure" approximation.

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.

Experimental and numerical study of the cold crucible melting process

The cold crucible, or induction skull melting process as is otherwise known, has the potential to produce high purity melts of a range of difficult to melt materials, including Ti–Al and Ti6Al4V alloys for Aerospace, Ti–Ta and other biocompatible materials for surgical implants, silicon for photovoltaic and electronic applications, etc. A water cooled AC coil surrounds the crucible causing induction currents to melt the alloy and partially suspend it against gravity away from water-cooled surfaces. Strong stirring takes place in the melt due to the induced electromagnetic Lorentz forces and very high temperatures are attainable under the right conditions (i.e., provided contact with water cooled walls is minimised). In a joint numerical and experimental research programme, various aspects of the design and operation of this process are investigated to increase our understanding of the physical mechanisms involved and to maximise process efficiency. A combination of FV and Spectral CFD techniques are used at Greenwich to tackle this problem numerically, with the experimental work taking place at Birmingham University. Results of this study, presented here, highlight the influence of turbulence and free surface behaviour on attained superheat and also discuss coil design variations and dual frequency options that may lead to winning crucible designs.

Maximising heat transfer efficiency in the cold crucible induction melting process

Induction heating is an efficient method used to melt electrically conductive materials, particularly if melting takes place in a ceramic crucible. This form of melting is particularly good for alloys, as electromagnetic forces set up by the induction coil lead to vigorous stirring of the melt ensuring homogeneity and uniformity in temperature. However, for certain reactive alloys, or where high purity is required, ceramic crucibles cannot be used, but a water-cooled segmented copper crucible is employed instead. Water cooling prevents meltdown or distortion of the metal wall, but much of the energy goes into the coolant. To reduce this loss, the electromagnetic force generated by the coil is used to push the melt away from the walls and so minimise contact with water-cooled surfaces. Even then, heat is lost through the crucible base where contact is inevitable. In a collaborative programme between Greenwich and Birmingham Universities, computer modelling has been used in conjunction with experiments to improve the superheat attainable in the melt for a,number of alloys, especially for y-TiAl intermetallics to cast aeroengine turbine blades. The model solves the discretised form of the turbulent Navier-Stokes, thermal energy conservation and Maxwell equations using a Spectral Collocation technique. The time-varying melt envelope is followed explicitly during the computation using an adaptive mesh. This paper briefly describes the mathematical model used to represent the interaction between the magnetic field, fluid flow, heat transfer and change of phase in the crucible and identifies the proportions of energy used in the melt, lost in the crucible base and in the crucible walls. The role of turbulence is highlighted as important in controlling heat losses and turbulence damping is introduced as a means of improving superheat. Model validation is against experimental results and shows good agreement with measured temperatures and energy losses in the cooling fluid throughout the melting cycle.

Experimental and numerical study of the cold crucible melting process

The cold crucible, or induction skull melting process as is otherwise known, has the potential to produce high purity melts of a range of difficult to melt materials, including Ti–Al and Ti6Al4V alloys for Aerospace, Ti–Ta and other biocompatible materials for surgical implants, silicon for photovoltaic and electronic applications, etc. A water cooled AC coil surrounds the crucible causing induction currents to melt the alloy and partially suspend it against gravity away from water-cooled surfaces. Strong stirring takes place in the melt due to the induced electromagnetic Lorentz forces and very high temperatures are attainable under the right conditions (i.e., provided contact with water cooled walls is minimised). In a joint numerical and experimental research programme, various aspects of the design and operation of this process are investigated to increase our understanding of the physical mechanisms involved and to maximise process efficiency. A combination of FV and Spectral CFD techniques are used at Greenwich to tackle this problem numerically, with the experimental work taking place at Birmingham University. Results of this study, presented here, highlight the influence of turbulence and free surface behaviour on attained superheat and also discuss coil design variations and dual frequency options that may lead to winning crucible designs.

The study of flow and temperature fields in conducting droplets suspended in a DC/AC combination field

Electromagnetic Levitation (EML) is a valuable method for measuring the thermo-physical properties of metals - surface tensions, viscosity, thermal/electrical conductivity, specific heat, hemispherical emissivity, etc. – beyond their melting temperature. In EML, a small amount of the test specimen is melted by Joule heating in a suspended AC coil. Once in liquid state, a small perturbation causes the liquid envelope to oscillate and the frequency of oscillation is then used to compute its surface tension by the well know Rayleigh formula. Similarly, the rate at which the oscillation is dampened relates to the viscosity. To measure thermal conductivity, a sinusoidally varying laser source may be used to heat the polar axis of the droplet and the temperature response measured at the polar opposite – the resulting phase shift yields thermal conductivity. All these theoretical methods assume that convective effects due to flow within the droplet are negligible compared to conduction, and similarly that the flow conditions are laminar; a situation that can only be realised under microgravity conditions. Hence the EML experiment is the method favoured for Spacelab experiments (viz. TEMPUS). Under terrestrial conditions, the full gravity force has to be countered by a much larger induced magnetic field. The magnetic field generates strong flow within the droplet, which for droplets of practical size becomes irrotational and turbulent. At the same time the droplet oscillation envelope is no longer ellipsoidal. Both these conditions invalidate simple theoretical models and prevent widespread EML use in terrestrial laboratories. The authors have shown in earlier publications that it is possible to suppress most of the turbulent convection generated in the droplet skin layer, through use of a static magnetic field. Using a pseudo-spectral discretisation method it is possible compute very accurately the dynamic variation in the suspended fluid envelope and simultaneously compute the time-varying electromagnetic, flow and thermal fields. The use of a DC field as a dampening agent was also demonstrated in cold crucible melting, where suppression of turbulence was achieved in a much larger liquid metal volume and led to increased superheat in the melt and reduction of heat losses to the water-cooled walls. In this paper, the authors describe the pseudo-spectral technique as applied to EML to compute the combined effects of AC and DC fields, accounting for all the flow-induced forces acting on the liquid volume (Lorentz, Maragoni, surface tension, gravity) and show example simulations.

Dynamic melting model for small samples in cold crucible

Purpose – A small size cold crucible offers possibilities for melting various electrically conducting materials with a minimal wall contact. Such small samples can be used for express contamination analysis, preparing limited amounts of reactive alloys or experimental material analyses. Aims to present a model to follow the melting process. Design/methodology/approach – The presents a numerical model in which different types of axisymmetric coil configurations are analysed. Findings – The presented numerical model permits dynamically to follow the melting process, the high-frequency magnetic field distribution change, the free surface and the melting front evolution, and the associated turbulent fluid dynamics. The partially solidified skin on the contact to the cold crucible walls and bottom is dynamically predicted. The segmented crucible shape is either cylindrical, hemispherical or arbitrary shaped. Originality/value – The model presented within the paper permits the analysis of melting times, melt shapes, electrical efficiency and particle tracks.

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.

Manufacturing of uniformly-sized silicon particles for solar cell from molten metal jet by electromagnetic pinch force

Spherical silicon solar cells are expected to serve as a technology to reduce silicon usage of photovoltaic (PV) power systems[1, 2, 3]. In order to establish the spherical silicon solar cell, a manufacturing method of uniformly sized silicon particles of 1mm in diameter is required. However, it is difficult to mass-produce the mono-sized silicon particles at low cost by existent processes now. We proposed a new method to generate liquid metal droplets uniformly by applying electromagnetic pinch force to a liquid metal jet[4]. The electromagnetic force was intermittently applied to the liquid metal jet issued from a nozzle in order to fluctuate the surface of the jet. As the fluctuation grew, the liquid jet was broken up into small droplets according to a frequency of the intermittent electromagnetic force. Firstly, a preliminary experiment was carried out. A single pulse current was applied instantaneously to a single turn coil around a molten gallium jet. It was confirmed that the jet could be split up by pinch force generated by the current. And then, electromagnetic pinch force was applied intermittently to the jet. It was found that the jet was broken up into mono-sized droplets in the case of a force frequency was equal to a critical frequency[5], which corresponds to a natural disturbance wave length of the jet. Numerical simulations of the droplet generation from the liquid jet were then carried out, which consisted of an electromagnetic analysis and a fluid flow calculation with a free surface of the jet. The simulation results were compared with the experiments and the agreement between the two was quite good.

Structure and denaturation of 4-chlorobenzoyl coenzyme A dehalogenase from Arthrobacter sp. strain TM-1

The secondary structure of the trimeric protein 4-chlorobenzoyl coenzyme A dehalogenase from Arthrobacter sp. strain TM-1, the second of three enzymes involved in the dechlorination of 4-chlorobenzoate to form 4-hydroxybenzoate, has been examined. E(mM) for the enzyme was 12.59. Analysis by circular dichroism spectrometry in the far uv indicated that 4-chlorobenzoyl coenzyme A dehalogenase was composed mostly of alpha-helix (56%) with lesser amounts of random coil (21%), beta-turn (13%) and beta-sheet (9%). These data are in close agreement with a computational prediction of secondary structure from the primary amino acid sequence, which indicated 55.8% alpha-helix, 33.7% random coil and 10.5% beta-sheet; the enzyme is, therefore, similar to the 4-chlorobenzoyl coenzyme A dehalogenase from Pseudomonas sp. CBS-3. The three-dimensional structure, including that of the presumed active site, predicted by computational analysis, is also closely similar to that of the Pseudomonas dehalogenase. Study of the stability and physicochemical properties revealed that at room temperature, the enzyme was stable for 24 h but was completely inactivated by heating to 60 degrees C for 5 min; thereafter by cooling at 1 degrees C min(-1) to 45 degrees C, 20.6% of the activity could be recovered. Mildly acidic (pH 5.2) or alkaline (pH 10.1) conditions caused complete inactivation, but activity was fully recovered on returning the enzyme to pH 7.4. Circular dichroism studies also indicated that secondary structure was little altered by heating to 60 degrees C, or by changing the pH from 7.4 to 6.0 or 9.2. Complete, irreversible destruction of, and maximal decrease in the fluorescence yield of the protein at 330-350 nm were brought about by 4.5 M urea or 1.1 M guanidinium chloride. Evidence was obtained to support the hypothetical three-dimensional model, that residues W140 and W167 are buried in a non-polar environment, whereas W182 appears at or close to the surface of the protein. At least one of the enzymes of the dehalogenase system (the combined 4-chlorobenzoate:CoA ligase, the dehalogenase and 4-hydroxybenzoyl coenzyme A thioesterase) appears to be capable of association with the cell membrane.