Monotone scheme and boundary conditions for finite volume simulation of magnetohydrodynamic internal flows at high Hartmann number

A monotone scheme for finite volume simulation of magnetohydrodynamic internal flows at high Hartmann number is presented. The numerical stability is analysed with respect to the electromagnetic force. Standard central finite differences applied to finite volumes can only be numerically stable if the vector products involved in this force are computed with a scheme using a fully staggered grid. The electromagnetic quantities (electric currents and electric potential) must be shifted by half the grid size from the mechanical ones (velocity and pressure). An integral treatment of the boundary layers is used in conjunction with boundary conditions for electrically conducting walls. The simulations are performed with inhomogeneous electrical conductivities of the walls and reach high Hartmann numbers in three-dimensional simulations, even though a non-adaptive grid is used.

The shallow water approximation applied to the aluminium electrolysis process

The industrial production of aluminium is an electrolysis process where two superposed horizontal liquid layers are subjected to a mainly vertical electric current supplied by carbon electrodes. The lower layer consists of molten aluminium and lies on the cathode. The upper layer is the electrolyte and is covered by the anode. The interface between the two layers is often perturbed, leading to oscillations, or waves, similar to the waves on the surface of seas or lakes. The presence of electric currents and the resulting magnetic field are responsible for electromagnetic (Lorentz) forces within the fluid, which can amplify these oscillations and have an adverse influence on the process. The electrolytic bath vertical to horizontal aspect ratio is such, that it is advantageous to use the shallow water equations to model the interface motion. These are the depth-averaging the Navier-Stokes equations so that nonlinear and dispersion terms may be taken into account. Although these terms are essential to the prediction of wave dynamics, they are neglected in most of the literature on interface instabilities in aluminium reduction cells where only the linear theory is usually considered. The unknown variables are the two horizontal components of the fluid velocity, the height of the interface and the electric potential. In this application, a finite volume resolution of the double-layer shallow water equations including the electromagnetic sources has been developed, for incorporation into a generic three-dimensional computational fluid dynamics code that also deals with heat transfer within the cell.

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.

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.

Characterisation of freeze-dried wafers and solvent evaporated films as potential drug delivery systems to mucosal surfaces

Freeze-dried (lyophilised) wafers and solvent cast films from sodium alginate (ALG) and sodium carboxymethylcellulose (CMC) have been developed as potential drug delivery systems for mucosal surfaces including wounds. The wafers (ALG, CMC) and films (CMC) were prepared by freeze-drying and drying in air (solvent evaporation) respectively, aqueous gels of the polymers containing paracetamol as a model drug. Microscopic architecture was examined using scanning electron microscopy, hydration characteristics with confocal laser scanning microscopy and dynamic vapour sorption. Texture analysis was employed to investigate mechanical characteristics of the wafers during compression. Differential scanning calorimetry was used to investigate polymorphic changes of paracetamol occurring during formulation of the wafers and films. The porous freeze-dried wafers exhibited higher drug loading and water absorption capacity than the corresponding solvent evaporated films. Moisture absorption, ease of hydration and mechanical behaviour were affected by the polymer and drug concentration. Two polymorphs of paracetamol were observed in the wafers and films, due to partial conversion of the original monoclinic to the orthorhombic polymorph during the formulation process. The results showed the potential of employing the freeze-dried wafers and solvent evaporated films in diverse mucosal applications due to their ease of hydration and based on different physical mechanical properties exhibited by both type of formulations.

Development and mechanical characterization of solvent-cast polymeric films as potential drug delivery systems to mucosal surfaces

Solvent-cast films from three polymers, carboxymethylcellulose (CMC), sodium alginate (SA), and xanthan gum, were prepared by drying the polymeric gels in air. Three methods, (a) passive hydration, (b) vortex hydration with heating, and (c) cold hydration, were investigated to determine the most effective means of preparing gels for each of the three polymers. Different drying conditions [relative humidity - RH (6-52%) and temperature (3-45 degrees C)] were investigated to determine the effect of drying rate on the films prepared by drying the polymeric gels. The tensile properties of the CMC films were determined by stretching dumbbell-shaped films to breaking point, using a Texture Analyser. Glycerol was used as a plasticizer, and its effects on the drying rate, physical appearance, and tensile properties of the resulting films were investigated. Vortex hydration with heating was the method of choice for preparing gels of SA and CMC, and cold hydration for xanthan gels. Drying rates increased with low glycerol content, high temperature, and low relative humidity. The residual water content of the films increased with increasing glycerol content and high relative humidity and decreased at higher temperatures. Generally, temperature affected the drying rate to a greater extent than relative humidity. Glycerol significantly affected the toughness (increased) and rigidity (decreased) of CMC films. CMC films prepared at 45 degrees C and 6% RH produced suitable films at the fastest rate while films containing equal quantities of glycerol and CMC possessed an ideal balance between flexibility and rigidity.