Magnetic modelling and magneto-hydro-dynamic simulations of an aluminium production electrolytic cell

A simulation of the motion of molten aluminium inside an electrolytic cell is presented. Since the driving term of the aluminium motion is the Lorentz (j × B) body force acting within the fluid,this problem involves the solution of the magneto-hydro-dynamic equations. Different solver modules for the magnetic field computation and for the fluid motion simulation are coupled together. The interactions of all these are presented and discussed.

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.

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].

Time dependent electric, magnetic and hydrodynamic interaction in aluminium electrolysis cells

The waves in commercial cells for electrolytic aluminium production originate at the interface between the liquid aluminium and electrolyte, but their effect can spread into the surrounding busbar network as electric current perturbation, and the total magnetic field acquires a time dependent component. The presented model for the wave development accounts for the nonuniform electric current distribution at the cathode and the whole network of the surrounding busbars. The magnetic field is computed for the continuous current in the fluid zones, all busbars and the ferromagnetic construction elements. When the electric current and the associated magnetic field are computed according to the actual electrical circuit and updated for all times, the instability growth rate is significantly affected. The presented numerical model for the wave and electromagnetic interaction demonstrates how different physical coupling factors are affecting the wave development in the electrolysis cells. These small amplitude self-sustained interface oscillations are damped in the presence of intense turbulent viscosity created by the horizontal circulation velocity field. Additionally, the horizontal circulation vortices create a pressure gradient contributing to the deformation of the interface. Instructive examples for the 500 kA demonstration cell are presented.

Computer modelling analysis of the globtop's effects on aluminium wirebond reliability

In power electronics modules, heavy aluminium wires, i.e. wire diameters greater than 100 microns, are bonded to the active semiconductor devices and conductor metallization to form electric circuits of the power electronic module. Due to the high currents that may flow through these wires, a single connection usually contains several wires and thus, a large number of wires are used in a power electronics module. Under normal operation or test condition, a significant amount of stresses and strains induced in the wire and bonding interfaces, resulting in failure over time. In this paper, computer modelling techniques are used to analyse the effect of globtop design on the reliability of aluminium wirebonds under cyclic thermal-mechanical loading conditions. The results will show the sensitivity of the reliability of the wirebonds to the changes in the geometry and the material properties of the wirebond globtop.

Computer simulation of aluminium wirebonds with globtop in power electronics modules

In this paper, computer modelling techniques are used to analyse the effects of globtops on the reliability of aluminium wirebonds in power electronics modules under cyclic thermal-mechanical loading conditions. The sensitivity of the wirehond reliability to the changes of the geometric and the material property parameters of wirebond globtop are evaluated and the optimal combination of the Young's modulus and the coefficient of thermal expansion have been predicted.

Shallow water model for aluminium electrolysis cells with variable top and bottom

The MHD wave instability in commercial cells for electrolytic aluminium production is often described using ‘shallow water’ models. The model [1] is extended for a variable height cathode bottom and anode top to account for realistic cell features. The variable depth of the two fluid layers affects the horizontal current density, the wave development and the stability threshold. Instructive examples for the 500 kA cell are presented.

Impact of the vertical potshell deformation on the MHD cell stability behavior of a 500 kA aluminium electrolysis cell

In previous publications [1,2], it was rationalized that a large vertical potshell deformation may have a negative impact on the operations of very high amperage cells. The MHD-Valdis non-linear Magneto-Hydro-Dynamic model was therefore extended to take into account the displacement of the potshell. The MHD cell stability behavior of a 500 kA cell with a 17.3 meters long potshell was then studied.

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.

Assessment of the origin and stability of Nifedipine polymorph IV by differential scanning calorimetry (DSC)

Purpose. To examine the thermal transition(s) between different polymorphic forms of Nifedipine and to define experimental conditions that lead to the generation of polymorph IV. Methods. Experiments were performed using a DSC 823e (Mettler Toledo). Nifedipine exists in four polymorphic forms, as well as an amorphous state. Examination of Nifedipine was conducted using the following method(s): cycle 1: 25ºC to 190ºC, 190ºC to 25ºC (formation of amorphous Nifedipine); cycle 2: 25ºC to X (60,70,80...150ºC), X to 25ºC; cycle 3: 25ºC to 190ºC and holding isothermally for 5 min between cycles (heating/cooling rate of 10ºC/min). Results. The amorphous state Nifedipine can sustain heating up to 90ºC without significant changes in its composition. Cycle 2 of amorphous material heated up to 90ºC shows only the glass transition at ~44ºC. In cycle 3 of the same material, a glass transition has been recorded at ~44ºC, followed by two exotherms (~100 and ~115ºC (crystallisation of polymorph III and II, respectively) and an endotherm (169ºC (melting of polymorphs I/II)). Samples that have been heated to temperatures between 100ºC and 120ºC in the second cycle showed a glass transition at ~44ºC and an additional exotherm at ~95ºC (crystallisation of polymorph III) on cooling a exotherm was observed at ~40ºC (crystallisation of polymorph IV). The same material showed no glass transition in cycle 3 but an endotherm at around 62ºC (melting of polymorph IV) an exotherm (~98ºC) and an endotherm (169ºC) melting of polymorph I/II. Heating the sample to a temperatures greater than 130ºC in cycle two results in a glass transition at ~44ºC, and two exotherms at ~102 and 125ºC (crystallisation of polymorphs III and I, respectively). Conclusions. DSC data suggests that polymorph IV can only be produced from amorphous or polymorph III samples. The presence of polymorph I or II drives the conversion of the less stable polymorphic form IV into the most stable form, I. Although form IV of Nifedipine can easily be created, following defined experimental conditions, it may only coexist with amorphous or polymorph III states. When polymorphs I and II are present in the sample polymorph IV cannot be etected.