Parallelisation of a 3-D electromagnetic scattering code

A high-frequency time domain finite element scattering code using a combination of edge and piecewise constant elements on unstructured tetrahedral meshes is described. A comparison of computation with theory is given for scattering from a sphere. A parallel implementation making use of the bulk synchronous parallel (BSP) programming model is described in detail; a BSP performance model of the parallelized field calculation is derived and compared to timing measurements on up to 128 processors on a Cray T3D.

Mathematical modelling of the solidification of liquid tin with electromagnetic stirring

It is well known that during alloy solidification, convection currents close to the so-lidification front have an influence on the structure of dendrites, the local solute concentration, the pattern of solid segregation, and eventually the microstructure of the casting and hence its mechanical properties. Controlled stirring of the melt in continuous casting or in ingot solidification is thought to have a beneficial effect. Free convection currents occur naturally due to temperature differences in the melt and for any given configuration, their strength is a function of the degree of superheat present. A more controlled forced convection current can be induced using electro-magnetic stirring. The authors have applied their Control-Volume based MHD method [1, 2] to the problem of tin solidification in an annular crucible with a water-cooled inner wall and a resistance heated outer one, for both free and forced convection situations and for various degrees of superheat. This problem was studied experimentally by Vives and Perry [3] who obtained temperature measurements, front positions and maps of electro-magnetic body force for a range of superheat values. The results of the mathematical model are compared critically against the experimental ones, in order to validate the model and also to demonstrate the usefulness of the coupled solution technique followed, as a predictive tool and a design aid. Figs 6, refs 19.

Dynamics of an axisymmetric electromagnetic ‘crucible’ melting

Different industrial induction melting processes involve free surface and melt-solid interface of the liquid metal subject to dynamic change during the technological operation. Simulation of the liquid metal dynamics requires to solve the non-linear, coupled hydrodynamic-electromagnetic-heat transfer problem accounting for the time development of the liquid metal free boundary with a suitable turbulent viscosity model. The present paper describes a numerical solution method applicable for various axisymmetric induction melting processes, such as, crucible with free top surface, levitation, semi-levitation, cold crucible and similar melting techniques. The presented results in the cases of semi-levitation and crucible with free top surface meltings demonstrate oscillating transient behaviour of the free metal surface indicating the presence of gravity-inertial-electromagnetic waves which are coupled to the internal fluid flow generated by both the rotational and potential parts of the electromagnetic force.

Neutron quasi-elastic scattering in disordered solids: a Monte Carlo study of metal-hydrogen systems

The dynamic structure factor of neutron quasi-elastic scattering has been calculated by Monte Carlo methods for atoms diffusing on a disordered lattice. The disorder includes not only variation in the distances between neighbouring atomic sites but also variation in the hopping rate associated with each site. The presence of the disorder, particularly the hopping rate disorder, causes changes in the time-dependent intermediate scattering function which translate into a significant increase in the intensity in the wings of the quasi-elastic spectrum as compared with the Lorentzian form. The effect is particularly marked at high values of the momentum transfer and at site occupancies of the order of unity. The MC calculations demonstrate how the degree of disorder may be derived from experimental measurements of the quasi-elastic scattering. The model structure factors are compared with the experimental quasi-elastic spectrum of an amorphous metal-hydrogen alloy.

Parallelisation of electromagnetic simulation codes

In this paper results obtained from the parallelisation of existing 3D electromagnetic Finite Element codes within the ESPRIT HPCN project PARTEL are presented. The parallelisation procedure, based on the Bulk Synchronous Parallel approach, is outlined and the encouraging results obtained in terms of speed-up on some industrially significant test cases are described and discussed.

Results of parallelisation of existing electromagnetic finite codes on some different parallel architectures

In this paper the results obtained from the parallelisation of some 3D industrial electromagnetic Finite Element codes within the ESPRIT Europort 2 project PARTEL are presented. The basic guidelines for the parallelisation procedure, based on the Bulk Synchronous Parallel approach, are presented and the encouraging results obtained in terms of speed-up on some selected test cases of practical design significance are outlined and discussed.

Phase-accuracy comparisons and improved far-field estimates for 3-D edge elements on tetrahedral meshes

Edge-element methods have proved very effective for 3-D electromagnetic computations and are widely used on unstructured meshes. However, the accuracy of standard edge elements can be criticised because of their low order. This paper analyses discrete dispersion relations together with numerical propagation accuracy to determine the effect of tetrahedral shape on the phase accuracy of standard 3-D edgeelement approximations in comparison to other methods. Scattering computations for the sphere obtained with edge elements are compared with results obtained with vertex elements, and a new formulation of the far-field integral approximations for use with edge elements is shown to give improved cross sections over conventional formulations.

Numerical simulation of free surface behaviour of a molten liquid metal droplet with and without electromagnetic induction

Electromagnetic levitation of electrically conductive droplets by alternating magnetic fields is a technique used to measure the physical properties of liquid metallic alloys such as surface tension or viscosity. Experiments can be conducted under terrestrial conditions or in microgravity, to reduce electromagnetic stirring and shaping of the droplet. Under such conditions, the time-dependent behaviour of a point of the free surface is recorded. Then the signal is analysed considering the droplet as a harmonic damped oscillator. We use a spectral code, for fluid flow and free surface descriptions, to check the validity of this assumption for two cases. First when the motion inside the droplet is generated by its initial distortion only and second, when the droplet is located in a uniform magnetic field originating far from the droplet. It is found that some deviations exist which can lead to an overestimate of the value of viscosity.

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.

Numerical modelling for electromagnetic processing of materials

Electromagnetic processing of materials (EPM) is one of the most widely practiced and fast growing applications of magnetic and electric forces to fluid flow. EPM is encountered in both industrial processes and laboratory investigations. Applications range in scale from nano-particle manipulation to tonnes of liquid metal treated in the presence of various configurations of magnetic fields. Some of these processes are specifically designed and made possible by the use of the electromagnetic force, like the magnetic levitation of liquid droplets, whilst others involve electric currents essential for electrothermal or electrochemical reasons, for instance, in electrolytic metal production and in induction melting. An insight for the range of established and novel EPM applications can be found in the review presented by Asai [1] in the EPM-2003 conference proceedings.

Droplet formation with electromagnetic pulse force

Methods for serial generation of droplets from a liquid jet are shortly reviewed. A method of liquid metal droplet generation based on AC high frequency magnetic field is considered in detail. Numerical model for direct simulation of the time dependent droplet generation process is presented. Computed examples demonstrate the liquid silicon droplet formation for the cases of 500-1500 μm diameter.

Effect of varying electromagnetic field on the VAR process

The three-dimensional, time-dependent electromagnetic field arising from the precession of the arc centre in a vacuum arc remelting furnace is shown (in a numerical simulation) to affect the fluid flow and heat transfer conditions near the solidification front in the upper part of the ingot.

Manufacturing of uniformly-sized silicon particles for solar cell from molten metal jet by electromagnetic pinch force

Spherical silicon solar cells are expected to serve as a technology to reduce silicon usage of photovoltaic (PV) power systems[1, 2, 3]. In order to establish the spherical silicon solar cell, a manufacturing method of uniformly sized silicon particles of 1mm in diameter is required. However, it is difficult to mass-produce the mono-sized silicon particles at low cost by existent processes now. We proposed a new method to generate liquid metal droplets uniformly by applying electromagnetic pinch force to a liquid metal jet[4]. The electromagnetic force was intermittently applied to the liquid metal jet issued from a nozzle in order to fluctuate the surface of the jet. As the fluctuation grew, the liquid jet was broken up into small droplets according to a frequency of the intermittent electromagnetic force. Firstly, a preliminary experiment was carried out. A single pulse current was applied instantaneously to a single turn coil around a molten gallium jet. It was confirmed that the jet could be split up by pinch force generated by the current. And then, electromagnetic pinch force was applied intermittently to the jet. It was found that the jet was broken up into mono-sized droplets in the case of a force frequency was equal to a critical frequency[5], which corresponds to a natural disturbance wave length of the jet. Numerical simulations of the droplet generation from the liquid jet were then carried out, which consisted of an electromagnetic analysis and a fluid flow calculation with a free surface of the jet. The simulation results were compared with the experiments and the agreement between the two was quite good.

Surface-enhanced Raman scattering studies of rhodanines: Evidence for substrate surface-induced dimerization

The surface-enhanced Raman scattering (SERS) spectra of rhodanine adsorbed on silver nanoparticles have been examined using 514.5 and 632.8 nm excitation. There is evidence that, under the experimental conditions used, rhodanine undergoes a nanoparticle surface-induced reaction resulting in the formation of a dimeric species via the active methylene group in a process which is analogous to the Knoevenagel reaction. The experimental observations are supported by DFT calculations at the B3-LYP/cc-pVDZ level. Calculated energies for the interaction of the E and Z isomers of the dimers of rhodanine with silver nanoparticles support a model in which the (intra-molecular hydrogen bonded) E isomer dimer is of lower energy than the Z isomer. A strong band, at 1566 cm(-1), in the SERS spectrum of rhodanine is assigned to the nu(C=C) mode of the dimer species.