Computational magnetohydrodynamics - Application

There are many processes, particularly in the nuclear and metals processing industries, where electromagnetic fields are used to influence the flow behaviour of a fluid. Procedures exploiting finite volume (FV) methods in both structured and unstructured meshes have recently been developed which enable this influence to be modelled in the context of conventional FV CFD codes. A range of problems have been tackled by the authors, including electromagnetic pumps and brakes, weirs and dams in steelmaking tundishes and interface effects in aluminium smelting cells. Two cases are presented here, which exemplify the application of the new procedures. The first case investigates the influence of electromagnetic fields on solidification front progression in a tin casting and the second case shows how the liquid metals free surface may be controlled through an externally imposed magnetic field in the semi-levitation casting process.

Models for magnetohydrodynamics of aluminium electrolysis cells

The electric current and the associated magnetic field in aluminium electrolysis cells create effects limiting the cell productivity and possibly cause instabilities: surface waving, ‘anode effects’, erosion of pot lining, feed material sedimentation, etc. The instructive analysis is presented via a step by step inclusion of different physical coupling factors affecting the magnetic field, electric current, velocity and wave development in the electrolysis cells. The full time dependent model couples the nonlinear turbulent fluid dynamics and the extended electromagnetic field in the cell, and the whole bus bar circuit with the ferromagnetic effects. Animated examples for the high amperage cells are presented. The theory and numerical model of the electrolysis cell is extended to the cases of variable cell bottom of aluminium layer and the variable thickness of the electrolyte due to the anode non-uniform burn-out process and the presence of the anode channels. The problem of the channel importance is well known Moreau-Evans model) for the stationary interface and the velocity field, and was validated against measurements in commercial cells, particularly with the recently published ‘benchmark’ test for the MHD models of aluminium cells [1]. The presence of electrolyte channels requires also to reconsider the previous magnetohydrodynamic instability theories and the dynamic wave development models. The results indicate the importance of a ‘sloshing’ parametrically excited MHD wave development in the aluminium production cells.

Finite volumes simulations in magnetohydrodynamics

Abstract not available

Mathematical modelling of the solidification of liquid tin with electromagnetic stirring

It is well known that during alloy solidification, convection currents close to the so-lidification front have an influence on the structure of dendrites, the local solute concentration, the pattern of solid segregation, and eventually the microstructure of the casting and hence its mechanical properties. Controlled stirring of the melt in continuous casting or in ingot solidification is thought to have a beneficial effect. Free convection currents occur naturally due to temperature differences in the melt and for any given configuration, their strength is a function of the degree of superheat present. A more controlled forced convection current can be induced using electro-magnetic stirring. The authors have applied their Control-Volume based MHD method [1, 2] to the problem of tin solidification in an annular crucible with a water-cooled inner wall and a resistance heated outer one, for both free and forced convection situations and for various degrees of superheat. This problem was studied experimentally by Vives and Perry [3] who obtained temperature measurements, front positions and maps of electro-magnetic body force for a range of superheat values. The results of the mathematical model are compared critically against the experimental ones, in order to validate the model and also to demonstrate the usefulness of the coupled solution technique followed, as a predictive tool and a design aid. Figs 6, refs 19.

Liquid metal turbulent flow dynamics in a cylindrical container with free surface: experiment and numerical analysis

The liquid metal flow in induction crucible models is known to be unstable, turbulent and difficult to predict in the regime of medium frequencies when the electromagnetic skin-layer is of considerable extent. We present long term turbulent flow measurements by a permanent magnet incorporated potential difference velocity probe in a cylindrical container filled with eutectic melt In-Ga-Sn. The parallel numerical simulation of the long time scale development of the turbulent average flow is presented. The numerical flow model uses an implicit pseudo-spectral code and k-w turbulence model, which was recently developed for the transitional flow modelling. The results compare reasonably to the experiment and demonstrate the time development of the turbulent flow field and the turbulence energy.

Buoyancy driven flow and its stability in a horizontal rectangular channel with an arbitrary oriented transversal magnetic field

An MHD flow is considered which is relevant to horizontal Bridgman technique for crystal growth from a melt. In the unidirectional parallel flow approximation an analytical solution is found accounting for the finite rectangular cross section of the channel in the case of a vertical magnetic field. Numerical pseudo-spectral solutions are used in the cases of arbitrary magnetic field and gravity vector orientations. The vertical magnetic field (parallel to the gravity) is found to be he most effective to damp the flow, however, complicated flow profiles with "overvelocities" in the comers are typical in the case of a finite cross-section channel. The temperature distribution is shown to be dependent on the flow profile. The linear stability of the flow is investigated by use of the Chebyshev pseudospectral method. For the case of an infinite width channel the transversal rolls instability is investigated, and for the finite cross-section channel the longitudinal rolls instability is considered. The critical Gr number values are computed in the dependence of the Ha number and the wave number or the aspect ratio in the case of finite section.

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.

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.

Numerical modelling for electromagnetic processing of materials

Electromagnetic processing of materials (EPM) is one of the most widely practiced and fast growing applications of magnetic and electric forces to fluid flow. EPM is encountered in both industrial processes and laboratory investigations. Applications range in scale from nano-particle manipulation to tonnes of liquid metal treated in the presence of various configurations of magnetic fields. Some of these processes are specifically designed and made possible by the use of the electromagnetic force, like the magnetic levitation of liquid droplets, whilst others involve electric currents essential for electrothermal or electrochemical reasons, for instance, in electrolytic metal production and in induction melting. An insight for the range of established and novel EPM applications can be found in the review presented by Asai [1] in the EPM-2003 conference proceedings.

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.

Vacuum arc remelting time dependent modelling

Newly developed numerical modelling tools are described, which address the 3-dimensional (3D) time-dependent magnetohydrodynamic and thermal behaviour in the liquid pool zone in the adjacent ingot, electrode and crucible. The melting electrode film flow and the droplet detachment initiation are simulated separately by an axisymmetric transient model.

Thermoelectric MHD in dendritic solidification

The effects of a constant uniform magnetic field on dendritic solidification were investigated using an enthalpy based numerical model. The interaction between thermoelectric currents on a growing crystal and the magnetic field generates a Lorentz force that creates flow. The need for very high resolution at the liquid-solid boundary where the thermoelectric source originates plus the need to accommodate multiple grains for a realistic simulation, make this a very demanding computational problem. For practical simulations, a quasi 3-dimensional approximation is proposed which nevertheless retains essential elements of transport in the third dimension. A magnetic field normal to the plane of growth leads to general flow circulation around an equiaxed dendrite, with secondary recirculations between the arms. The heat/solute advection by the flow is shown to cause a change in the morphology of the dendrite; secondary growth is promoted preferentially on one side of the dendrite arm and the tip velocity of the primary arm is increased. The degree of approximation introduced is quantified by extending the model into 3-dimensions, where the full Navier-Stokes equation is solved, and compared against the 2-dimensional solution.

Levitated droplet oscillations: effect of internal flow

The numerical model for electrically conducting liquid droplets levitated in AC magnetic field is applied to demonstrate various factors affecting the accuracy of material property value measurements in microgravity conditions. The included effects are the electromagnetic force induced stirring and the resulting turbulence, the thermo-capillary convection, and the droplet rotation. The results are validated against available analytical solutions.