An FDTD model for calculation of gradient-induced eddy currents in MRI system

In most magnetic resonance imaging (MRI) systems, pulsed magnetic gradient fields induce eddy currents in the conducting structures of the superconducting magnet. The eddy currents induced in structures within the cryostat are particularly problematic as they are characterized by long time constants by virtue of the low resistivity of the conductors. This paper presents a three-dimensional (3-D) finite-difference time-domain (FDTD) scheme in cylindrical coordinates for eddy-current calculation in conductors. This model is intended to be part of a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The singularity apparent in the governing equations is removed by using a series expansion method and the conductor-air boundary condition is handled using a variant of the surface impedance concept. The numerical difficulty due to the asymmetry of Maxwell equations for low-frequency eddy-current problems is circumvented by taking advantage of the known penetration behavior of the eddy-current fields. A perfectly matched layer absorbing boundary condition in 3-D cylindrical coordinates is also incorporated. The numerical method has been verified against analytical solutions for simple cases. Finally, the algorithm is illustrated by modeling a pulsed field gradient coil system within an MRI magnet system. The results demonstrate that the proposed FDTD scheme can be used to calculate large-scale eddy-current problems in materials with high conductivity at low frequencies.

A reinvestigation of phase equilibria in the system Al2O3-SiO2-ZnO

The phase equilibria and liquidus temperatures in the binary SiO2-ZnO system and in the ternary Al2O3-SiO2-ZnO system at low Al2O3 concentrations have been experimentally determined using the equilibration and quenching technique followed by electron probe X-ray microanalysis. In the SiO2-ZnO system, two binary eutectics involving the congruently melting willemite (Zn2SiO4) were found at 1448 +/- 5 degrees C and 0.52 +/- 0.01 mole fraction ZnO and at 1502 +/- 5 degrees C and 0.71 +/- 0.01 mole fraction ZnO, respectively. The two ternary eutectics involving willemite previously reported in the Al2O3SiO2-ZnO system were found to be at 1315 +/- 5 degrees C and 1425 +/- 25 T, respectively. The compositions of the eutectics are 0.07, 0.52, and 0.41 and 0.05, 0.28, and 0.67 mole fraction Al2O3, SiO2, and ZnO, respectively. The results of the present investigation are significantly different from the results of previous studies.

Phase equilibria in the system Fe-Zn-O at intermediate conditions between metallic-iron saturation and air

The phase equilibria in the FeO-Fe2O3-ZnO system have been experimentally investigated at oxygen partial pressures between metallic iron saturation and air using a specially developed quenching technique, followed by electron probe X-ray microanalysis (EPMA) and then wet chemistry for determination of ferrous and ferric iron concentrations. Gas mixtures of H-2, N-2, and CO2 or CO and CO2 controlled the atmosphere in the furnace. The determined metal cation ratios in phases at equilibrium were used for the construction of the 1200 degrees C isothermal section of the Fe-Zn-O system. The univariant equilibria between the gas phase, spinel, wustite, and zincite was found to be close to pO(2) = 1 center dot 10(-8) atm at 1200 degrees C. The ferric and ferrous iron concentrations in zincite and spinel at equilibrium were also determined at temperatures from 1200 degrees C to 1400 degrees C at pO(2) = 1 center dot 10(-6) atm and at 1200 degrees C at pO(2) values ranging from 1 center dot 10(-4) to 1 center dot 10(-8) atm. Implications of the phase equilibria in the Fe-Zn-O system for the formation of the platelike zincite, especially important for the Imperial Smelting Process (ISP), are discussed.

Positioning hierarchy in enterprise system change

There is a widening gulf in change literature between theoretical notions of evolving organisational form and the emerging reality that old and new organisational structures coexist. This paper explores this dichotomy in Enterprise Resource Planning change. It develops a cellular hierarchy framework to explain how different types of hierarchy coexist within the same organisation during the implementation of Enterprise Resource Planning. © 2006 The Author; Journal compilation © 2006 Blackwell Publishing Ltd.

The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.