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.

Vacuum arc remelting time dependent modelling

Vacuum arc remelting (VAR) aims at production of high quality, segregation-free alloys. The quality of the produced ingots depends on the operating conditions which could be monitored and analyzed using numerical modelling. The remelting process uniformity is controlled by critical medium scale time variations of the order 1-100 s, which are physically initiated by the droplet detachment and the large scale arc motion at the top of liquid pool [1,2]. The newly developed numerical modelling tools are addressing the 3-dimensional magnetohydrodynamic and thermal behaviour in the liquid zone and the adjacent ingot, electrode and crucible.

Time-dependent modelling and experimental validation of the metal/flux interface in a continuous casting mould

A finite volume computer model of the continuous casting process for steel flat products has been developed. In this first stage, the model concentrates on the hydrodynamic aspects of the process and in particular the dynamic behavior of the metal/slag interface. The model was validated against experimental measurements obtained in a water model apparatus.

Modelling the time-dependent behaviour of solder pastes used in the electronics assembly applications

Solder paste is the most widely used bonding material in the assembly of surface mount devices in electronic industries. It generally has a flocculated structure (show aggregation of solder particles), and hence are known to exhibit a thixotropic behavior. This is recognized by the decrease in apparent viscosity of paste material with time when subjected to a constant shear rate. The proper characterisation of this timedependent rheological behaviour of solder pastes is crucial for establishing the relationships between the pastes’ structure and flow behaviour; and for correlating the physical parameters with paste printing performance. In this paper, we present a novel method which has been developed for characterising the timedependent and non-Newtonian rheological behaviour of solder pastes as a function of shear rates. The objective of the study reported in this paper is to investigate the thixotropic build-up behaviour of solder pastes. The stretched exponential model(SEM) has been used to model the structural changes during the build-up process and to correlate model parameters with the paste printing process.

Vacuum arc remelting time dependent modelling

Newly developed numerical modelling tools are described, which address the 3-dimensional (3D) time-dependent magnetohydrodynamic and thermal behaviour in the liquid pool zone in the adjacent ingot, electrode and crucible. The melting electrode film flow and the droplet detachment initiation are simulated separately by an axisymmetric transient model.

Modelling of the time-dependent flow behaviour of lead-free solder pastes used for flip-chip assembly applications

The market for solder paste materials in the electronic manufacturing and assembly sector is very large and consists of material and equipment suppliers and end users. These materials are used to bond electronic components (such as flip-chip, CSP and BGA) to printed circuit boards (PCB's) across a range of dimensions where the solder interconnects can be in the order of 0.05mm to 5mm in size. The non-Newtonian flow properties exhibited by solder pastes during its manufacture and printing/deposition phases have been of practical concern to surface mount engineers and researchers for many years. The printing of paste materials through very small-sized stencil apertures is known to lead to increased stencil clogging and incomplete transfer of paste to the substrate pads. At these very narrow aperture sizes the paste rheology and particle-wall interactions become crucial for consistent paste withdrawal. These non-Newtonian effects must be understood so that the new paste formulations can be optimised for consistent printing. The focus of the study reported in this paper is the characterisation of the rheological properties of solder pastes and flux mediums, and the evaluation of the effect of these properties on the pastes' printing performance at the flip-chip assembly application level. Solder pastes are known to exhibit a thixotropic behaviour, which is recognised by the decrease in apparent viscosity of paste material with time when subjected to a constant shear rate. The proper characterisation of this time-dependent theological behaviour of solder pastes is crucial for establishing the relationships between the pastes' structure and flow behaviour; and for correlating the physical parameters with paste printing performance. In this paper, we present a number of methods which have been developed for characterising the time-dependent and non-Newtonian rheological behaviour of solder pastes and flux mediums as a function of shear rates. We also present results of the study of the rheology of the solder pastes and flux mediums using the structural kinetic modelling approach, which postulates that the network structure of solder pastes breaks down irreversibly under shear, leading to time and shear dependent changes in the flow properties. Our results show that for the solder pastes used in the study, the rate and extent of thixotropy was generally found to increase with increasing shear rate. The technique demonstrated in this study has wide utility for R&D personnel involved in new paste formulation, for implementing quality control procedures used in solder paste manufacture and packaging; and for qualifying new flip-chip assembly lines

Time domain decomposition for European options in financial modelling

Finance is one of the fastest growing areas in modern applied mathematics with real world applications. The interest of this branch of applied mathematics is best described by an example involving shares. Shareholders of a company receive dividends which come from the profit made by the company. The proceeds of the company, once it is taken over or wound up, will also be distributed to shareholders. Therefore shares have a value that reflects the views of investors about the likely dividend payments and capital growth of the company. Obviously such value will be quantified by the share price on stock exchanges. Therefore financial modelling serves to understand the correlations between asset and movements of buy/sell in order to reduce risk. Such activities depend on financial analysis tools being available to the trader with which he can make rapid and systematic evaluation of buy/sell contracts. There are other financial activities and it is not an intention of this paper to discuss all of these activities. The main concern of this paper is to propose a parallel algorithm for the numerical solution of an European option. This paper is organised as follows. First, a brief introduction is given of a simple mathematical model for European options and possible numerical schemes of solving such mathematical model. Second, Laplace transform is applied to the mathematical model which leads to a set of parametric equations where solutions of different parametric equations may be found concurrently. Numerical inverse Laplace transform is done by means of an inversion algorithm developed by Stehfast. The scalability of the algorithm in a distributed environment is demonstrated. Third, a performance analysis of the present algorithm is compared with a spatial domain decomposition developed particularly for time-dependent heat equation. Finally, a number of issues are discussed and future work suggested.

Modeling the dynamics of Magnetic Semilevitation Melting

In semilevitation melting, a cylindrical metal ingot is melted by a coaxial a.c. induction coil. A watercooled solid base supports the ingot, while the top and side free surface is confined by the magnetic forces as the melting front progresses. The dynamic interplay between gravity, hydrodynamic stress, and the Lorentz force in the fluid determines the instantaneous free surface shape. The coupled nonstationary equations for turbulent flow, heat with phase change, and high-frequency electromagnetic field are solved numerically for the axisymmetric time-dependent domain by a continuous mesh transformation, using a pseudospectral method. Results are obtained for the two actually existing coil configurations and several validation cases.

Birth-death processes with mass annihilation and state-dependent immigration

A birth-death process is subject to mass annihilation at rate β with subsequent mass immigration occurring into state j at rateα j . This structure enables the process to jump from one sector of state space to another one (via state 0) with transition rate independent of population size. First, we highlight the difficulties encountered when using standard techniques to construct both time-dependent and equilibrium probabilities. Then we show how to overcome such analytic difficulties by means of a tool developed in Chen and Renshaw (1990, 1993b); this approach is applicable to many processes whose underlying generator on E\{0} has known probability structure. Here we demonstrate the technique through application to the linear birth-death generator on which is superimposed an annihilation/immigration process.

On a parallel time-domain method for the nonlinear Black-Scholes equation

A parallel time-domain algorithm is described for the time-dependent nonlinear Black-Scholes equation, which may be used to build financial analysis tools to help traders making rapid and systematic evaluation of buy/sell contracts. The algorithm is particularly suitable for problems that do not require fine details at each intermediate time step, and hence the method applies well for the present problem.

Ontological considerations of time, meta-predicates and temporal propositions

A natural approach to representing and reasoning about temporal propositions (i.e., statements with time-dependent truth-values) is to associate them with time elements. In the literature, there are three choices regarding the primitive for the ontology of time: (1) instantaneous points, (2) durative intervals and (3) both points and intervals. Problems may arise when one conflates different views of temporal structure and questions whether some certain types of temporal propositions can be validly and meaningfully associated with different time elements. In this paper, we shall summarize an ontological glossary with respect to time elements, and diversify a wider range of meta-predicates for ascribing temporal propositions to time elements. Based on these, we shall also devise a versatile categorization of temporal propositions, which can subsume those representative categories proposed in the literature, including that of Vendler, of McDermott, of Allen, of Shoham, of Galton and of Terenziani and Torasso. It is demonstrated that the new categorization of propositions, together with the proposed range of meta-predicates, provides the expressive power for modeling some typical temporal terms/phenomena, such as starting-instant, stopping-instant, dividing-instant, instigation, termination and intermingling etc.

Numerical simulation of flow-induced cavity noise in self-sustained oscillations

The generation and near-field radiation of aerodynamic sound from a low-speed unsteady flow over a two-dimensional automobile door cavity is simulated by using a source-extraction-based coupling method. In the coupling procedure, the unsteady cavity flow field is first computed solving the Reynolds averaged Navier–Stokes (RANS) equations. The radiated sound is then calculated by using a set of acoustic perturbation equations with acoustic source terms which are extracted from the time-dependent solutions of the unsteady flow. The aerodynamic and its resulting acoustic field are computed for the Reynolds number of 53,266 based on the base length of the cavity. The free stream flow velocity is taken to be 50.9m/s. As first stage of the numerical investigation of flow-induced cavity noise, laminar flow is assumed. The CFD solver is based on a cell-centered finite volume method. A dispersion-relation-preserving (DRP), optimized, fourth-order finite difference scheme with fully staggered-grid implementation is used in the acoustic solver

Numerical simulation of incompressible flow problems using an unstructured staggered mesh method

A number of two dimensional staggered unstructured discretisation schemes for the solution of fluid flow and heat transfer problems have been developed. All schemes store and solve velocity vector components at cell faces with scalar variables solved at cell centres. The velocity is resolved into face-normal and face-parallel components and the various schemes investigated differ in the treatment of the parallel component. Steady-state and time-dependent fluid flow and thermal energy equations are solved with the well known pressure correction scheme, SIMPLE, employed to couple continuity and momentum. The numerical methods developed are tested on well known benchmark cases: the Lid-Driven Cavity, Natural Convection in a Cavity and Melting of Gallium in a rectangular domain. The results obtained are shown to be comparable to benchmark, but with accuracy dependent on scheme selection.

A unified 3D model of the VAR process

A 3D time-dependent model of the VAR process has been developed using CFD techniques. The model solves the coupled field equations for fluid flow, heat transfer (including phase change) and electromagnetic field, for both the electrode and the ingot. The motion of the electic arc 'preferred spot' can be specified based on observations. Correlations are sought between the local gap height, resulting from instantaneous liquid pool surface shape and electrode tip shape, and the arc motion. The detailed behaviour of the melting film on the electrode tip is studies using a spectral free surface technique, which allows investigation of the drops' detachment and drip shorts.