Pseudo-spectral solutions for fluid flow and heat transfer in electro-metallurgical applications

The pseudo-spectral solution method offers a flexible and fast alternative to the more usual finite element/volume/difference methods, particularly when the long-time transient behaviour of a system is of interest. Since the exact solution is obtained at the grid collocation points superior accuracy can be achieved on modest grid resolution. Furthermore, the grid can be freely adapted with time and in space, to particular flow conditions or geometric variations. This is especially advantageous where strongly coupled, time-dependent, multi-physics solutions are investigated. Examples include metallurgical applications involving the interaction of electromagnetic fields and conducting liquids with a free sutface. The electromagnetic field then determines the instantaneous liquid volume shape and the liquid shape affects in turn the electromagnetic field. In AC applications a thin "skin effect" region results on the free surface that dominates grid requirements. Infinitesimally thin boundary cells can be introduced using Chebyshev polynomial expansions without detriment to the numerical accuracy. This paper presents a general methodology of the pseudo-spectral approach and outlines the solution procedures used. Several instructive example applications are given: the aluminium electrolysis MHD problem, induction melting and stirring and the dynamics of magnetically levitated droplets in AC and DC fields. Comparisons to available analytical solutions and to experimental measurements will be discussed.

Pseudo-spectral solutions for fluid flow and heat transfer in electro-metallurgical applications

The pseudo-spectral solution method offers a flexible and fast alternative to the more usual finite element and volume methods, particularly when the long-time transient behaviour of a system is of interest. The exact solution is obtained at grid collocation points leading to superior accuracy on modest grids. Furthermore, the grid can be freely adapted in time and space to particular flow conditions or geometric variations, especially useful where strongly coupled, time-dependent, multi-physics solutions are investigated. Examples include metallurgical applications involving the interaction of electromagnetic fields and conducting liquids with a free surface. The electromagnetic field determines the instantaneous liquid volume shape, which then affects the electromagnetic field. A general methodology of the pseudo-spectral approach is presented, with several instructive example applications: the aluminium electrolysis MHD problem, induction melting in a cold crucible and the dynamics of AC/DC magnetically levitated droplets. Finally, comparisons with available analytical solutions and to experimental measurements are discussed.

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.

Modelling induction melting energy savings

Electromagnetic processing of liquid metals involves dynamic change of the fluid volume interfacing with a melting solid material, gas or vacuum, and possibly a different liquid. Electromagnetic field and the associated force field are strongly coupled to the free surface dynamics and the heat-mass transfer. We present practical modelling examples of the flow and heat transfer using an accurate pseudo-spectral code and the k-omega turbulence model suitable for complex and transitional flows with free surfaces. The 'cold crucible' melting is modelled dynamically including the melting front gradual propagation and the magnetically confined free surrounding interface. Intermittent contact with the water-cooled segmented wall and the radiation heat losses are parts of the complex problem.

Scheduling batches with sequential job processing for two-machine flow and open shops

In this paper, we study a problem of scheduling and batching on two machines in a flow-shop and open-shop environment. Each machine processes operations in batches, and the processing time of a batch is the sum of the processing times of the operations in that batch. A setup time, which depends only on the machine, is required before a batch is processed on a machine, and all jobs in a batch remain at the machine until the entire batch is processed. The aim is to make batching and sequencing decisions, which specify a partition of the jobs into batches on each machine, and a processing order of the batches on each machine, respectively, so that the makespan is minimized. The flow-shop problem is shown to be strongly NP-hard. We demonstrate that there is an optimal solution with the same batches on the two machines; we refer to these as consistent batches. A heuristic is developed that selects the best schedule among several with one, two, or three consistent batches, and is shown to have a worst-case performance ratio of 4/3. For the open-shop, we show that the problem is NP-hard in the ordinary sense. By proving the existence of an optimal solution with one, two or three consistent batches, a close relationship is established with the problem of scheduling two or three identical parallel machines to minimize the makespan. This allows a pseudo-polynomial algorithm to be derived, and various heuristic methods to be suggested.

Modeling the dynamics of electromagnetically agitated metallurgical flows: levitation, cold crucible, aluminium electrolysis, etc.

We present practical modelling techniques for electromagnetically agitated liquid metal flows involving dynamic change of the fluid volume and shape during melting and the free surface oscillation. Typically the electromagnetic field is strongly coupled to the free surface dynamics and the heat-mass transfer. Accurate pseudo-spectral code and the k-omega turbulence model modified for complex and transitional flows with free surfaces are used for these simulations. The considered examples include magnetic suspension melting, induction scull remelting (cold crucible), levitation and aluminium electrolysis cells. The process control and the energy savings issues are analysed.

Investigation of the cold crucible melting process: experimental and numerical study

The dynamic process of melting different materials in a cold crucible is being studied experimentally with parallel numerical modelling work. The numerical simulation uses a variety of complementing models: finite volume, integral equation and pseudo-spectral methods combined to achieve the accurate description of the dynamic melting process. Results show the temperature history of the melting process with a comparison of the experimental and computed heat losses in the various parts of the equipment. The free surface visual observations are compared to the numerically predicted surface shapes.

Modelling electromagnetically levitated liquid droplet oscillations

This work comprises accurate computational analysis of levitated liquid droplet oscillations in AC and DC magnetic fields. The AC magnetic field interacting with the induced electric current within the liquid metal droplet generates intense fluid flow and the coupled free surface oscillations. The pseudo-spectral technique is used to solve the turbulent fluid flow equations for the continuously dynamically transformed axisymmetric fluid volume. The volume electromagnetic force distribution is updated with the shape and position change. We start with the ideal fluid test case for undamped Rayleigh frequency oscillations in the absence of gravity, and then add the viscous and the DC magnetic field damping. The oscillation frequency spectra are further analysed for droplets levitated against gravity in AC and DC magnetic fields at various combinations. In the extreme case electrically poorly conducting, diamagnetic droplet (water) levitation dynamics are simulated. Applications are aimed at pure electromagnetic material processing techniques and the material properties measurements in uncontaminated conditions.

Droplet generation with modulated AC electromagnetic field at nozzle exit

A method of droplet generation based on applying a modulated AC high frequency magnetic field in the localized region of capillary breakup is considered as ans alternative to traditional methods for high temperature liquid melt droplet serial production by pressure variation. The method is based on a pseudo-spectral approximation with a coordinate transformation adaptin to the developing free surface. The electromagnetic field is recomputed continuously with the domain shape change. Practical application cases for liquid silicon droplets of 0.5 - 2 mm diameter are considered in detail.

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.

Thermoelectric effects on alloy solidification microstructure

Thermoelectric currents in the presence of a magnetic field generate Lorentz forces which can drive fluid flow. In the case of dendritic growth a naturally occurring thermoelectric current exists and in the presence of a high magnetic field micro convections are generated. Experimental evidence has attributed changes in microstructure to this effect. A numerical model has been developed to study the flow field around an unconstricted equiaxed dendrite growing under these conditions. The growth is modeled in 2D and 3D by an enthalpy based method and a complex flow structure has been predicted. Using a pseudo-3D approximation for economy, realistic 2D simulations are obtained where a fully coupled transient scheme reveals significant changes to the dendrite morphology reflecting experimental evidence. There is a rotation of the preferred direction of growth and increased secondary branching.

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.

Aspects of the in vitro bioactivity of hydraulic calcium (alumino)silicate cement

In response to a burgeoning interest in the prospective clinical applications of hydraulic calcium (alumino)silicate cements, the in vitro bioactivity and dissolution characteristics of a white Portland cement have been investigated. The formation of an apatite layer within 6 h of contact with simulated body fluid was attributed to the rapid dissolution of calcium hydroxide from the cement matrix and to the abundance of pre-existing Si-OH nucleation sites presented by the calcium silicate hydrate phase. A simple kinetic model has been used to describe the rate of apatite formation and an apparent pseudo-second-order rate constant for the removal of HPO42- ions frorn solultion has been calculated (k(2) = 5.8 x 10(-4) g mg(-1)). Aspects of the chemistry of hydraulic cements are also discussed with respect to their potential use in the remedial treatment of living tissue. (C) 2008 Wiley Periodicals, Inc. J Biomed Mater Res 90A: 166-174, 2009

Ag+- and Zn2+-exchange kinetics and antimicrobial properties of 11 angstrom tobermorites

Ag+- and Zn2+-exchanged zeolites zeolites and clays have been used as coatings and in composites to confer broad-spectrum antimicrobial properties on a range of technical and biomedical materials. 11 angstrom tobermorite is a bioactive layer lattice ion exchanger whose potential as a carrier for Ag+ and Zn2+ ions in antimicrobial formulations has not yet been explored. In view of this, batch Ag+- and Zn2+-exchange kinetics of two structurally distinct synthetic 11 angstrom tobermorites and their subsequent bactericidal action against Staphylococcus aureus and Pseudomonas aeruginosa are reported. During the exchange reactions, Ag+ ions were found to replace labile interlayer cations; whereas, Zn2+ ions also displaced structural Ca2+ ions from the tobermorite lattice. In spite of these different mechanisms, a simple pseudo-second-order model provided a suitable description of both exchange processes (R-2 >= 0.996). The Ag+- and Zn2+-exchanged tobermorite phases exhibited marked bacteriostatic effects against both bacteria, and accordingly, their potential for use as antimicrobial materials for in situ bone tissue regeneration is discussed. (C) 2008 Elsevier Ltd. All rights reserved.