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.

Modelling induction skull melting design modifications

Induction Skull Melting (ISM) is used for heating, melting, mixing and, possibly, evaporating reactive liquid metals at high temperatures when a minimum contact at solid walls is required. The numerical model presented here involves the complete time dependent process analysis based on the coupled electromagnetic, temperature and turbulent velocity fields during the melting and liquid shape changes. The simulation is validated against measurements of liquid metal height, temperature and heat losses in a commercial size ISM furnace. The often observed limiting temperature plateau for ever increasing electrical power input is explained by the turbulent convective heat losses. Various methods to increase the superheat within the liquid melt, the process energy efficiency and stability are proposed.

Nonlinear MHD stability of aluminium reduction cells

An industrial electrolysis cell used to produce primary aluminium is sensitive to waves at the interface of liquid aluminium and electrolyte. The interface waves are similar to stratified sea layers [1], but the penetrating electric current and the associated magnetic field are intricately involved in the oscillation process, and the observed wave frequencies are shifted from the purely hydrodynamic ones [2]. The interface stability problem is of great practical importance because the electrolytic aluminium production is a major electrical energy consumer, and it is related to environmental pollution rate. The stability analysis was started in [3] and a short summary of the main developments is given in [2]. Important aspects of the multiple mode interaction have been introduced in [4], and a widely used linear friction law first applied in [5]. In [6] a systematic perturbation expansion is developed for the fluid dynamics and electric current problems permitting reduction of the three-dimensional problem to a two dimensional one. The procedure is more generally known as “shallow water approximation” which can be extended for the case of weakly non-linear and dispersive waves. The Boussinesq formulation permits to generalise the problem for non-unidirectionally propagating waves accounting for side walls and for a two fluid layer interface [1]. Attempts to extend the electrolytic cell wave modelling to the weakly nonlinear case have started in [7] where the basic equations are derived, including the nonlinearity and linear dispersion terms. An alternative approach for the nonlinear numerical simulation for an electrolysis cell wave evolution is attempted in [8 and references there], yet, omitting the dispersion terms and without a proper account for the dissipation, the model can predict unstable waves growth only. The present paper contains a generalisation of the previous non linear wave equations [7] by accounting for the turbulent horizontal circulation flows in the two fluid layers. The inclusion of the turbulence model is essential in order to explain the small amplitude self-sustained oscillations of the liquid metal surface observed in real cells, known as “MHD noise”. The fluid dynamic model is coupled to the extended electromagnetic simulation including not only the fluid layers, but the whole bus bar circuit and the ferromagnetic effects [9].

Cold crucible melting of reactive metals using combined DC and AC magnetic fields

Cold crucible furnace is widely used for melting reactive metals for high quality castings. Although the water cooled copper crucible avoids contamination, it produces a low superheat of the melt. Experimental and theoretical investigations of the process showed that the increase of the supplied power to the furnace leads to a saturation in the temperature rise of the melt, and no significant increase of the melt superheat can be obtained. The computer model of theprocess has been developed to simulate the time dependent turbulent flow, heat transfer with phase change, and AC and DC magnetohydrodynamics in a time varying liquid metal envelope. The model predicts that the supermimposition of a strong DC field on top of the normal AC field reduces the level of turbulience and stirring in the liquid metal, thereby reducing the heat loss through the base of the crucible and increasing the superheat. The direct measurements of the temperature in the commercial size cold crucbile has confirmed the computer redictions and showed that the addition of a DC field increased the superheat in molten TiAl from ~45C (AC field only) to ~81C (DC+AC fields). The present paper reports further predictions of the effect of a dDC field on top of the AC field and compares these with experimental data.

Droplet formation with electromagnetic pulse force

Methods for serial generation of droplets from a liquid jet are shortly reviewed. A method of liquid metal droplet generation based on AC high frequency magnetic field is considered in detail. Numerical model for direct simulation of the time dependent droplet generation process is presented. Computed examples demonstrate the liquid silicon droplet formation for the cases of 500-1500 μm diameter.

Multiscale modelling of dendrite growth during vacuum arc remelting

Vacuum Arc Remelting (VAR) is the accepted method for producing homogeneous, fine microstructures that are free of inclusions required for rotating grade applications. However, as ingot sizes are increasing INCONEL 718 becomes increasingly susceptible to defects such as freckles, tree rings, and white spots increases for large diameter billets. Therefore, predictive models of these defects are required to allow optimization of process parameters. In this paper, a multiscale and multi-physics model is presented to predict the development of microstructures in the VAR ingot during solidification. At the microscale, a combined stochastic nucleation approach and finite difference solution of the solute diffusion is applied in the semi-solid zone of the VAR ingot. The micromodel is coupled with a solution of the macroscale heat transfer, fluid flow and electromagnetism in the VAR process through the temperature, pressure and fluid flow fields. The main objective of this study is to achieve a better understanding of the formation of the defects in VAR by quantifying the influence of VAR processing parameters on grain nucleation and dendrite growth. In particular, the effect of different ingot growth velocities on the microstructure formation was investigated. It was found that reducing the velocity produces significantly more coarse grains.

The tilt casting of titanium aluminides

TiAl castings are prone to various defects including bubbles entrained during the turbulent filling of moulds. The present research has exploited the principles of the Durville tilt casting technique to develop a novel process in which the Induction Skull Melting (ISM) of TiAl alloys in a vacuum chamber has been combined with controlled tilt pouring to achieve the tranquil transfer of the metal into a hot ceramic shell mould. Practical casting equipment has been developed to evaluate the feasibility of this process in parallel with the development of novel software to simulate and optimize it. The PHYSICA CFD code was used to simulate the filling, heat transfer and solidification during tilt pouring using a number of free surface modelling techniques, including the novel Counter Diffusion Method (CDM). In view of the limited superheat, particular attention was paid to the mould design to minimize heat loss and gas entrainment caused by interaction between the counter-flowing metal and gas streams. The model has been validated against real-time X-ray movies of the tilt casting of aluminium and against TiAl blade castings. Modelling has contributed to designing a mould to promote progressive filling of the casting and has led to the use of a parabolic tilting cycle to balance the competing requirements for rapid filling to minimize the loss of superheat and slow filling minimize the turbulence-induced defects.

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.

Multiscale modeling of the onset of freckle formation during vacuum arc remelting

A multiscale model for the Vacuum Arc Remelting process (VAR) was developed to simulate dendritic microstructures during solidification and investigate the onset of freckle formation. On the macroscale, a 3D multi-physics model of VAR was used to study complex physical phenomena, including liquid metal flow with turbulence, heat transfer, and magnetohydrodynamics. The results showed that unsteady fluid flow in the liquid pool caused significant thermal perturbation at the solidification front. These results were coupled into a micromodel to simulate dendritic growth controlled by solute diffusion, including local remelting. The changes in Rayleigh number as the microstructure remelts was quantified to provide an indicator of when fluid flow channels (i.e. freckles) will initiate in the mushy zone. By examining the simulated microstructures, it was found that the Rayleigh number increased more than 300 times during remelting, which suggests that thermal perturbation could be responsible for the onset of freckle formation.

A unified 3D model of the VAR process

A 3D time-dependent model of the VAR process has been developed using CFD techniques. The model solves the coupled field equations for fluid flow, heat transfer (including phase change) and electromagnetic field, for both the electrode and the ingot. The motion of the electic arc 'preferred spot' can be specified based on observations. Correlations are sought between the local gap height, resulting from instantaneous liquid pool surface shape and electrode tip shape, and the arc motion. The detailed behaviour of the melting film on the electrode tip is studies using a spectral free surface technique, which allows investigation of the drops' detachment and drip shorts.

Quantized conductance of multiwalled carbon nanotubes

We report calculations of the transport properties of multiwalled carbon nanotubes based on a scattering-theoretic approach that takes into account scattering within each tube, between tubes, and at the metal contacts. We find that only the outer tube contributes to the conductance, as has been implied by experiments. Referring to experiments performed with liquid-metal contacts, we also explain why the measured conductance is close to an integer number of conductance quanta, when the tubes are immersed in the liquid metal for several hundreds of nanometers and is not an integer when they are immersed for only a few nanometers. Finally, we propose that the observed conductance of only one quantum

Effects of temperature and strain rate on the tensile behaviors of SIMP steel in static lead bismuth eutectic

In order to assess the susceptibility of candidate structural materials to liquid metal embrittlement, this work investigated the tensile behaviors of ferritic-martensitic steel in static lead bismuth eutectic (LBE). The tensile tests were carried out in static lead bismuth eutectic under different temperatures and strain rates. Pronounced liquid metal embrittlement phenomenon is observed between 200 °C and 450 °C. Total elongation is reduced greatly due to the liquid metal embrittlement in LBE environment. The range of ductility trough is larger under slow strain rate tensile (SSRT) test.

Experimental model of the interfacial instability in aluminium reduction cells

A solution has been found to the long-standing problem of experimental modelling of the interfacial instability in aluminium reduction cells. The idea is to replace the electrolyte overlaying molten aluminium with a mesh of thin rods supplying current down directly into the liquid metal layer. This eliminates electrolysis altogether and all the problems associated with it, such as high temperature, chemical aggressiveness of media, products of electrolysis, the necessity for electrolyte renewal, high power demands, etc. The result is a room temperature, versatile laboratory model which simulates Sele-type, rolling pad interfacial instability. Our new, safe laboratory model enables detailed experimental investigations to test the existing theoretical models for the first time.

Contributions to the theory of magnetorotational instability and waves in a rotating plasma

The one-fluid magnetohydrodynamic (MHD) theory of magnetorotational instability (MRI) in an ideal plasma is presented. The theory predicts the possibility of MRI for arbitrary 0, where 0 is the ratio of the plasma pressure to the magnetic field pressure. The kinetic theory of MRI in a collisionless plasma is developed. It is demonstrated that as in the ideal MHD, MRI can occur in such a plasma for arbitrary P. The mechanism of MRI is discussed; it is shown that the instability appears because of a perturbed parallel electric field. The electrodynamic description of MRI is formulated under the assumption that the dispersion relation is expressed in terms of the permittivity tensor; general properties of this tensor are analyzed. It is shown to be separated into the nonrotational and rotational parts. With this in mind, the first step for incorporation of MRI into the general theory of plasma instabilities is taken. The rotation effects on Alfven waves are considered.