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

Experimental and Modelling Analysis on the Moisture Induced Failures in Flip Chip on Flex Interconnections with Anisotropic Conductive Film

This paper reports the investigations into the moisture induced failures in flip-chip-on-flex interconnections with anisotropic conductive films (ACF). Both experimental and modeling methods were applied. In the experiments, the contact resistance was used as a quality indicator and was measured continuously during the accelerated tests (autoclave tests). The temperature, relative humidity and the pressure were set at 121°C, 100%RH, 1atm respectively. The contact resistance of the ACF joints increased during the tests and nearly 25% of the joints were found to be open after 168 hours' testing time. Visible conduction gaps between the adhesive and substrate pads were observed. Cracks at the adhesive/flex interface were also found. It is believed that the swelling effect of the adhesive and the water penetration along the adhesive/flex interface are the main causes of this contact degradation. Another finding from the experimental work was that the ACF interconnections that had undergone the reflow treatment were more sensitive to the moisture and showed worse reliability during the tests. For a better understanding of the experimental results, 3D finite element (FE) models were built and a macro-micro modeling method was used to determine the moisture diffusion and moisture-induced stresses inside the ACF joints. Modeling results are consistent with the findings in the experimental work.

Modeling and experiments on an isothermal fatigue test for solder joints

This paper investigates an isothermal fatigue test for solder joints developed at the NPL. The test specimen is a lap joint between two copper arms. During the test the displacement at the ends of the copper are controlled and the force measured. The modeling results in the paper show that the displacement across the solder joint is not equal to the displacement applied at the end of the specimen. This is due to deformation within the copper arms. A method is described to compensate for this difference. The strain distribution in the solder was determined by finite element analysis and compared to the distribution generated by a theoretical 'ideal' test which generates an almost pure shear mode in the solder. By using a damage-based constitutive law the shape of the crack generated in the specimen has been predicted for both the actual test and the ideal pure shear test. Results from the simulations are also compared with experimental data using SnAgCu solder.

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.

Decision support systems for eco-friendly electronic products

This paper discusses an optimisation based decision support system and methodology for electronic packaging and product design and development which is capable of addressing in efficient manner specified environmental, reliability and cost requirements. A study which focuses on the design of a flip-chip package is presented. Different alternatives for the design of the flip-chip package are considered based on existing options for the applied underfill and volume of solder material used to form the interconnects. Variations in these design input parameters have simultaneous effect on package aspects such as cost, environmental impact and reliability. A decision system for the design of the flip-chip that uses numerical optimisation approach is used to identify the package optimal specification which satisfies the imposed requirements. The reliability aspect of interest is the fatigue of solder joints under thermal cycling. Transient nonlinear finite element analysis (FEA) is used to simulate the thermal fatigue damage in solder joints subject to thermal cycling. Simulation results are manipulated within design of experiments and response surface modelling framework to provide numerical model for reliability which can be used to quantify the package reliability. Assessment of the environmental impact of the package materials is performed by using so called Toxic Index (TI). In this paper we demonstrate the evaluation of the environmental impact only for underfill and lead-free solder materials. This evaluation is based on the amount of material per flip-chip package. Cost is the dominant factor in contemporary flip-chip packaging industry. In the optimisation based decision support system for the design of the flip-chip package, cost of materials which varies as a result of variations in the design parameters is considered.

Modelling the behavior of solder joints for wafer level SiP

Design for manufacture of system-in-package (SiP) structures is dependent on a number of physical processes that affect the final quality of the package in terms of its performance and reliability. Solder joints are key structures in a SiP and their behavior can be the critical factor in terms of reliability. This paper discusses the results from a research programme on design for manufacturing of system in package (SiP) technologies. The focus of the paper is on thermo-mechanical modelling of solder joints. This includes the behavior of the joints during testing plus some important insights into the reflow process and how physical phenomena taking place at the assembly stage can affect solder joint behavior. Finite element analysis of a numerical model of an SiP structure with various design parameters is discussed. The goal of this analysis is to identify the most promising combination of design parameters which guarantee longer lifetime of the solder joints and hence the SiP component. The parameters that were studied are the size of the package (i.e. number of solder joints per row), the presence of the underfill and/or the reinforcement as well as the thickness of the passive die. Discussion was also provided on phenomena that take place during the reflow process where the solder joints are formed. In particular, the formation of intermetallics at the solder-pad interfaces

Constitutive modeling of kinematic-hardening and damage in solder joints

Solder constitutive models are important as they are widely used in FEA simulations to predict the lifetime of soldered assemblies. This paper briefly reviews some common constitutive laws to capture creep in solder and presents work on laws capturing both kinematic hardening and damage. Inverse analysis is used to determine constants for the kinematic hardening law which match experimental creep curves. The mesh dependence of the damage law is overcome by using volume averaging and is applied to predict the crack path in a thermal cycled resistor component

Shear strength analysis of ball grid array (BGA) solder interfaces

Ball shear test is the most common test method used to assess the reliability of bond strength for ball grid array (BGA) packages. In this work, a combined experimental and numerical study was carried out to realize of BGA solder interface strength. Solder mask defined bond pads on the BGA substrate were used for BGA ball bonding. Different bond pad metallizations and solder alloys were used. Solid state aging at 150degC up to 1000 h has been carried out to change the interfacial microstructure. Cross-sectional studies of the solder-to-bond pad interfaces was conducted by scanning electron microscopy (SEM) equipped with an energy dispersive X-ray (EDX) analyzer to investigate the interfacial reaction phenomena. Ball shear tests have been carried out to obtain the mechanical strength of the solder joints and to correlate shear behaviour with the interfacial reaction products. An attempt has been taken to realize experimental findings by Finite Element Analysis (FEA). It was found that intermetallic compound (IMC) formation at the solder interface plays an important role in the BGA solder bond strength. By changing the morphology and the microchemistry of IMCs, the fracture propagation path could be changed and hence, reliability could be improved

Fracture mechanics analysis of cracks in solder joint intermetallic compounds

The trend towards miniaturization of electronic products leads to the need for very small sized solder joints. Therefore, there is a higher reliability risk that too large a fraction of solder joints will transform into Intermetallic Compounds (IMCs) at the solder interface. In this paper, fracture mechanics study of the IMC layer for SnPb and Pb-free solder joints was carried out using finite element numerical computer modelling method. It is assumed that only one crack is present in the IMC layer. Linear Elastic Fracture Mechanics (LEFM) approach is used for parametric study of the Stress Intensity Factors (SIF, KI and KII), at the predefined crack in the IMC layer of solder butt joint tensile sample. Contrary to intuition, it is revealed that a thicker IMC layer in fact increases the reliability of solder joint for a cracked IMC. Value of KI and KII are found to decrease with the location of the crack further away from the solder interfaces while other parameters are constant. Solder thickness and strain rate were also found to have a significant influence on the SIF values. It has been found that soft solder matrix generates non-uniform plastic deformation across the solder-IMC interface near the crack tip that is responsible to obtain higher KI and KII.

Failure mechanisms of ACF joints under isothermal ageing

The possible failure mechanisms of anisotropic conductive film (ACF) joints under isothermal ageing conditions have been identified through experiments. It has been found that ACF joints formed at higher bonding temperatures can prevent increases in the contact resistance for any ageing temperature. The higher the ageing temperature the higher the electrical failure rate is. The formation of conduction gaps between the conductive particles and the pads and damages to the metal coatings of the particle have been identified as the reasons behind the electrical failures during ageing. In order to understand the mechanism for the formation of the conduction gap and damages in metal coatings during the isothermal ageing, computer modelling has been carried out and the results are discussed extensively. The computer analysis shows that stresses concentrate at the edges of the particle–pad interface, where the adhesive matrix meets the particle. This could lead to subsequent damages and reductions in the adhesion strength in that region and it is possible for the conductive particle to be detached from the pad and the adhesive matrix. It is believed that because of this a conduction gap appears. Furthermore, under thermal loading the thermal expansion of the adhesive matrix squeezes the conductive particle and damages the metal coatings. Experimental evidences support this computational finding. It is, therefore, postulated that if an ACF-based electronic component operates in a high temperature aging condition, its electrical and mechanical functionalities will be at risk.

Computer simulation of crack propagation in power electronics module solder joints

A numerical modelling method for the analysis of solder joint damage and crack propagation has been described in this paper. The method is based on the disturbed state concept. Under cyclic thermal-mechanical loading conditions, the level of damage that occurs in solder joints is assumed to be a simple monotonic scalar function of the accumulated equivalent plastic strain. The increase of damage leads to crack initiation and propagation. By tracking the evolution of the damage level in solder joints, crack propagation path and rate can be simulated using Finite Element Analysis method. The discussions are focused on issues in the implementation of the method. The technique of speeding up the simulation and the mesh dependency issues are analysed. As an example of the application of this method, crack propagation in solder joints of power electronics modules under cyclic thermal-mechanical loading conditions has been analyzed and the predicted cracked area size after 3000 loading cycles is consistent with experimental results.

Numerical solutions of certain nonlinear models in European options on a distributed computing environment

Financial modelling in the area of option pricing involves the understanding of the correlations between asset and movements of buy/sell in order to reduce risk in investment. 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. In turn, analysis tools rely on fast numerical algorithms for the solution of financial mathematical models. There are many different financial activities apart from shares buy/sell activities. The main aim of this chapter is to discuss a distributed algorithm for the numerical solution of a European option. Both linear and non-linear cases are considered. The algorithm is based on the concept of the Laplace transform and its numerical inverse. The scalability of the algorithm is examined. Numerical tests are used to demonstrate the effectiveness of the algorithm for financial analysis. Time dependent functions for volatility and interest rates are also discussed. Applications of the algorithm to non-linear Black-Scholes equation where the volatility and the interest rate are functions of the option value are included. Some qualitative results of the convergence behaviour of the algorithm is examined. This chapter also examines the various computational issues of the Laplace transformation method in terms of distributed computing. The idea of using a two-level temporal mesh in order to achieve distributed computation along the temporal axis is introduced. Finally, the chapter ends with some conclusions.

Fault Diagnosis in Internal Combustion Engines Using Nonlinear Multivariate Statistics

This paper presents a statistical-based fault diagnosis scheme for application to internal combustion engines. The scheme relies on an identified model that describes the relationships between a set of recorded engine variables using principal component analysis (PCA). Since combustion cycles are complex in nature and produce nonlinear relationships between the recorded engine variables, the paper proposes the use of nonlinear PCA (NLPCA). The paper further justifies the use of NLPCA by comparing the model accuracy of the NLPCA model with that of a linear PCA model. A new nonlinear variable reconstruction algorithm and bivariate scatter plots are proposed for fault isolation, following the application of NLPCA. The proposed technique allows the diagnosis of different fault types under steady-state operating conditions. More precisely, nonlinear variable reconstruction can remove the fault signature from the recorded engine data, which allows the identification and isolation of the root cause of abnormal engine behaviour. The paper shows that this can lead to (i) an enhanced identification of potential root causes of abnormal events and (ii) the masking of faulty sensor readings. The effectiveness of the enhanced NLPCA based monitoring scheme is illustrated by its application to a sensor fault and a process fault. The sensor fault relates to a drift in the fuel flow reading, whilst the process fault relates to a partial blockage of the intercooler. These faults are introduced to a Volkswagen TDI 1.9 Litre diesel engine mounted on an experimental engine test bench facility.

Modified stiffened panel analysis methods for laser beam and friction stir welded aircraft panels

The introduction of advanced welding methods as an alternative joining process to riveting in the manufacture of primary aircraft structure has the potential to realize reductions in both manufacturing costs and structural weight. Current design and analysis methods for aircraft panels have been developed and validated for riveted fabrication. For welded panels, considering the buckling collapse design philosophy of aircraft stiffened panels, strength prediction methods considering welding process effects for both local-buckling and post-buckling behaviours must be developed and validated. This article reports on the work undertaken to develop analysis methods for the crippling failure of stiffened panels fabricated using laser beam and friction stir welding. The work assesses modifications to conventional analysis methods and finite-element analysis methods for strength prediction. The analysis work is validated experimentally with welded single stiffener crippling specimens. The experimental programme has demonstrated the potential static strength of laser beam and friction stir welded sheet-stiffener joints for post-buckling panel applications. The work undertaken has demonstrated that the crippling behaviour of welded stiffened panels may be analysed considering standard-buckling behaviour. However, stiffened panel buckling analysis procedures must be altered to account for the weld joint geometry and process altered material properties. © IMechE 2006.

A two-stage algorithm for identification of nonlinear dynamic systems

This paper investigates the two-stage stepwise identification for a class of nonlinear dynamic systems that can be described by linear-in-the-parameters models, and the model has to be built from a very large pool of basis functions or model terms. The main objective is to improve the compactness of the model that is obtained by the forward stepwise methods, while retaining the computational efficiency. The proposed algorithm first generates an initial model using a forward stepwise procedure. The significance of each selected term is then reviewed at the second stage and all insignificant ones are replaced, resulting in an optimised compact model with significantly improved performance. The main contribution of this paper is that these two stages are performed within a well-defined regression context, leading to significantly reduced computational complexity. The efficiency of the algorithm is confirmed by the computational complexity analysis, and its effectiveness is demonstrated by the simulation results.