A finite-volume shallow layer method for the MHD instabilities in an aluminium production cell

An electrolytic cell for Aluminium production contains molten metal and molten electrolyte, which are subject to high dc-currents and magnetic fields. Lorentz forces arising from the cross product of current and magnetic field may amplify natural gravity waves at the interface between the two fluids, leading to short circuits in extreme cases. The external magnetic field and current distribution in the production cell is computed through a detailed finite element analysis at Torino Polytechnic. The results are then used to compute the magnetohydrodynamic and thermal effects in the aluminium/electrolyte bath. Each cell has lateral dimensions of 6m x 2m, whilst the bath depth is only 30cm. the electrically resistive electrolyte path, which is critical in the operation of the cell, has layer depth of only a few centimetres below each carbon anode. Because the shallow dimensions of the liquid layer a finite-volume shallow-layer technique has been used at Greenwich to compute the resulting flow-field and interface perturbations. The information obtained from this method, i.e. depth averaged velocities and aluminium/electrolyte interface position is then embedded in the three-dimensional finite volume code PHYSICA and will be used to compute the heat transfer and phase change in the cell.

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.

Determining surface tension using DC magnetic levitation

The values of material physical properties are vital for the successful use of numerical simulations for electromagnetic processing of materials. The surface tension of materials can be determined from the experimental measurement of the surface oscillation frequency of liquid droplets. In order for this technique to be used, a positioning field is required that results in a modification to the oscillation frequency. A number of previous analytical models have been developed that mainly focus on electrically conducting droplets positioned using an A.C. electromagnetic field, but due to the turbulent flow resulting from the high electromagnetic fields required to balance gravity, reliable measurements have largely been limited to microgravity. In this work axisymmetric analytical and numerical models are developed, which allow the surface tension of a diamagnetic droplet positioned in a high DC magnetic field to be determined from the surface oscillations. In the case of D.C. levitation there is no internal electric currents with resulting Joule heating, Marangoni flow and other effects that introduce additional physics that complicates the measurement process. The analytical solution uses the linearised Navier-Stokes equations in the inviscid case. The body force from a DC field is potential, in contrast to the AC case, and it can be derived from Maxwell equations giving a solution for the magnetic field in the form of a series expansion of Legendre polynomials. The first few terms in this expansion represent a constant and gradient magnetic field valid close to the origin, which can be used to position the droplet. Initially the mathematical model is verified in microgravity conditions using a numerical model developed to solve the transient electromagnetics, fluid flow and thermodynamic equations. In the numerical model (as in experiment) the magnetic field is obtained using electrical current carrying coils, which provides the confinement force for a liquid droplet. The model incorporates free surface deformation to accurately model the oscillations that result from the interaction between the droplet and the non-uniform external magnetic field. A comparison is made between the analytical perturbation theory and the numerical pseudo spectral approximation solutions for small amplitude oscillations.

The development of dynamic CFD model of the magnetic semi-levitation of liquid metals

In the casting of reactive metals, such as titanium alloys, contamination can be prevented if there is no contact between the hot liquid metal and solid crucible. This can be achieved by containing the liquid metal by means of high frequency AC magnetic field. A water cooled current-carrying coil, surrounding the metal can then provide the required Lorentz forces, and at the same time the current induced in the metal can provide the heating required to melt it. This ‘attractive’ processing solution has however many problems, the most serious being that of the control and containment of the liquid metal envelope, which requires a balance of the gravity and induced inertia forces on the one side, and the containing Lorentz and surface tension forces on the other. To model this process requires a fully coupled dyna ic solution of the flow fields, magnetic field and heat transfer/melding process to account for. A simplified solution has been published previously providing quasi-static solutions only, by taking the irrotational ‘magnetic pressure’ term of the Lorentz force into account. The authors remedy this deficiency by modelling the full problem using CFD techniques. The salient features of these techniques are included in this paper, as space allows.