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.

Domain decomposition using a 2-level correction scheme

The PHYSICA software was developed to enable multiphysics modelling allowing for interaction between Computational Fluid Dynamics (CFD) and Computational Solid Mechanics (CSM) and Computational Aeroacoustics (CAA). PHYSICA uses the finite volume method with 3-D unstructured meshes to enable the modelling of complex geometries. Many engineering applications involve significant computational time which needs to be reduced by means of a faster solution method or parallel and high performance algorithms. It is well known that multigrid methods serve as a fast iterative scheme for linear and nonlinear diffusion problems. This papers attempts to address two major issues of this iterative solver, including parallelisation of multigrid methods and their applications to time dependent multiscale problems.

Guided waves in a fluid-loaded transversely isotropic place

Dispersion relations are obtained for the propagation of symmetric and antisymmetric modes in a free transversely isotropic plate. Dispersion curves are plotted for the first four symmetric modes for a magnesium plate immersed in water. The first mode is highly damped and switches over to the second mode when the normalized frequency exceeds 12.

Investigation into the performance of turbulence models for fluid flow and heat transfer phenomena in electronic applications

Computational Fluid Dynamics (CFD) is gradually becoming a powerful and almost essential tool for the design, development and optimization of engineering applications. However the mathematical modelling of the erratic turbulent motion remains the key issue when tackling such flow phenomena. The reliability of CFD analysis depends heavily on the turbulence model employed together with the wall functions implemented. In order to resolve the abrupt changes in the turbulent energy and other parameters situated at near wall regions a particularly fine mesh is necessary which inevitably increases the computer storage and run-time requirements. Turbulence modelling can be considered to be one of the three key elements in CFD. Precise mathematical theories have evolved for the other two key elements, grid generation and algorithm development. The principal objective of turbulence modelling is to enhance computational procedures of efficient accuracy to reproduce the main structures of three dimensional fluid flows. The flow within an electronic system can be characterized as being in a transitional state due to the low velocities and relatively small dimensions encountered. This paper presents simulated CFD results for an investigation into the predictive capability of turbulence models when considering both fluid flow and heat transfer phenomena. Also a new two-layer hybrid kε / kl turbulence model for electronic application areas will be presented which holds the advantages of being cheap in terms of the computational mesh required and is also economical with regards to run-time.

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

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.

Effect of adding 0.3wt% Ni into the Sn-0.7wt% Cu solder - Part II: Growth of intermetallic layer with Cu during wetting and aging

The growth behavior of intermetallic layer with or without adding 0.3 wt% Ni into the Sn-0.7Cu solder was studied during the wetting reaction on Cu-substrate and thereafter in solid-state aging condition. The Cu-solder reaction couple was prepared at 255, 275 and 295 °C for 10 s. The samples reacted at 255 °C were then isothermally aged for 2-14 days at 150 °C. The reaction species formed for the Sn-0.7Cu/Cu and Sn-0.7Cu-0.3Ni/Cu soldering systems were Cu6Sn5 and (CuNi)6Sn5, respectively. The thickness of the intermetallic compounds formed at the solder/Cu interfaces and also in the bulk of both solders increased with the increase of reaction temperature. It was found that Ni-containing Sn-0.7Cu solder exhibited lower growth of intermetallic layer during wetting and in the early stage of aging and eventually exceeded the intermetallic layer thickness of Sn-0.7Cu/Cu soldering system after 6 days of aging. As the aging time proceeds, a non-uniform intermetallic layer growth tendency was observed for the case of Sn-0.7Cu-0.3Ni solder. The growth behavior of intermetallic layer during aging for both solders followed the diffusion-controlled mechanism. The intermetallic layer growth rate constants for Sn-0.7Cu and Sn-0.7Cu-0.3Ni solders were calculated as 1.41 × 10-17 and 1.89 × 10-17 m2/s, respectively which indicated that adding 0.3 wt% Ni with Sn-0.7Cu solder contributed to the higher growth of intermetallic layer during aging. © 2006 Elsevier B.V. All rights reserved.

Dynamic model for metal cleanness evaluation by melting in a cold crucible

Melting of metallic samples in a cold crucible causes inclusions to concentrate on the surface owing to the action of the electromagnetic force in the skin layer. This process is dynamic, involving the melting stage, then quasi-stationary particle separation, and finally the solidification in the cold crucible. The proposed modeling technique is based on the pseudospectral solution method for coupled turbulent fluid flow, thermal and electromagnetic fields within the time varying fluid volume contained by the free surface, and partially the solid crucible wall. The model uses two methods for particle tracking: (1) a direct Lagrangian particle path computation and (2) a drifting concentration model. Lagrangian tracking is implemented for arbitrary unsteady flow. A specific numerical time integration scheme is implemented using implicit advancement that permits relatively large time-steps in the Lagrangian model. The drifting concentration model is based on a local equilibrium drift velocity assumption. Both methods are compared and demonstrated to give qualitatively similar results for stationary flow situations. The particular results presented are obtained for iron alloys. Small size particles of the order of 1 μm are shown to be less prone to separation by electromagnetic field action. In contrast, larger particles, 10 to 100 μm, are easily “trapped” by the electromagnetic field and stay on the sample surface at predetermined locations depending on their size and properties. The model allows optimization for melting power, geometry, and solidification rate.

Effect of fluid dynamics and device mechanism on biofluid behaviour in microchannel systems: modelling biofluids in a microchannel biochip separator

Biofluid behaviour in microchannel systems is investigated in this paper through the modelling of a microfluidic biochip developed for the separation of blood plasma. Based on particular assumptions, the effects of some mechanical features of the microchannels on behaviour of the biofluid are explored. These include microchannel, constriction, bending channel, bifurcation as well as channel length ratio between the main and side channels. The key characteristics and effects of the microfluidic dynamics are discussed in terms of separation efficiency of the red blood cells with respect to the rest of the medium. The effects include the Fahraeus and Fahraeus-Lindqvist effects, the Zweifach-Fung bifurcation law, the cell-free layer phenomenon. The characteristics of the microfluid dynamics include the properties of the laminar flow as well as particle lateral or spinning trajectories. In this paper the fluid is modelled as a single-phase flow assuming either Newtonian or Non-Newtonian behaviours to investigate the effect of the viscosity on flow and separation efficiency. It is found that, for a flow rate controlled Newtonian flow system, viscosity and outlet pressure have little effect on velocity distribution. When the fluid is assumed to be Non-Newtonian more fluid is separated than observed in the Newtonian case, leading to reduction of the flow rate ratio between the main and side channels as well as the system pressure as a whole.

Aspects of the in vitro bioactivity and antimicrobial properties of Ag+- and Zn2+-exchanged 11 Å tobermorites

11 Å tobermorite, Ca5Si6O16(OH)2 · 4H2O, is a layer lattice ion exchange mineral whose potential as a carrier for Ag+ and Zn2+ ions in antimicrobial, bioactive formulations has not yet been explored. In view of this, the in vitro bioactivity of Ag+- and Zn2+-exchanged 11 Å tobermorites and their bactericidal action against S. aureus and P.aeruginosa are reported. The in vitro bioactivity of the synthetic unsubstituted tobermorite phase was confirmed by the formation of bone-like hydroxycarbonate apatite (HCA) on its surface within 48 h of contact with simulated body fluid. The substitution of labile Ag+ ions into the tobermorite lattice delayed the onset of HCA-formation to 72 h; whereas, the Zn2+-substituted phase failed to elicit an HCA-layer within 14 days. Both Ag+- and Zn2+-exchanged tobermorite phases were found to exhibit marked antimicrobial action against S. aureus and P.aeruginosa, two common pathogens in biomaterial-centred infections.

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

Structural and biological investigation of ferrocene-substituted 3-methylidene-1,3-dihydro-2H-indol-2-ones

The Knoevenagel condensation of 1,3-dihydro-2H-indol-2-one with ferrocene carboxaldehyde afforded an approximate 2:1 mixture of the geometrical isomers (E)- and (Z)-3-ferrocenylmethylidene-1,3-dihydro-2H-indol-2-one respectively in an overall 67% yield; the air and solution-stable isomers were readily separated by preparative thin layer chromatography and their structures were unequivocally elucidated in solution, by (1)H NMR spectroscopy, and in the solid phase, by X-ray crystallography; both isomers of displayed in vitro toxicity against B16 melanoma and Vero cell lines in the micromolar range and inhibited the kinase VEGFR-2 with IC(50) values of ca. 200 nM.