Multi-physics modelling for the fabrication, packaging and reliability of micro-systems components

The latest advances in multi-physics modelling both using high fidelity techniques and reduced order and behavioural models will be discussed. Particular focus will be given to the application and validation of these techniques for modelling the fabrication, packaging and subsequent reliability of micro-systems based components. The paper will discuss results from a number of research projects with particular emphasis on the techniques being developed in a major UK Goverment funded project - 3D-MINTEGRATION (www.3d-mintegration.com).

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.

Modelling the tilt-casting process for the tranquil filling of titanium alloy turbine blades

The tilt-casting method is used to achieve tranquil filling of gamma-TiAl turbine blades. The reactive alloy is melted in a cold crucible using an induction coil and then the complete crucible-mould- running system assembly is rotated through 180degrees to transfer the metal into the mould. The induction current is ramped down gradually as the rotation starts and the mould is preheated to maintain superheat. The liquid metal then enters the mould and the gas within it (argon) escapes through the inlet aperture and through auxiliary vents. Solidification starts as soon the metal enters the mould and it is important to account for this effect to predict and prevent misruns. The rotation rate has to be controlled carefully to allow sufficient time for gas evacuation, but at the same time preserve superheat. This 3-phase system is modelled using the FV method, with a fast implicit numerical scheme used to capture the transient liquid free surface. The enthalpy method is used to model solidification and predict defects such as trapped bubbles, macro-porosity or surface connected porosity. Modeling is used to support an experimental program for the development of a production method for gamma-TiAl blades, with a target length of 40cm. The experiments provide validation for the model and the model in turn optimizes the tilt-casting process. The work is part of the EU project IMPRESS.

Evacuation modelling analysis within the operational research context: a combined approach for improving enclosure designs

Evacuation models have been playing an important function in the transition process from prescriptive fire safety codes to performance-based ones over the last three decades. In fact, such models became also useful tools in different tasks within fire safety engineering field, such as fire risks assessment and fire investigation. However, there are some difficulties in this process when using these models. For instance, during the evacuation modelling analysis, a common problem faced by fire safety engineers concerns the number of simulations which needs to be performed. In other terms, which fire designs (i.e., scenarios) should be investigated using the evacuation models? This type of question becomes more complex when specific issues such as the optimal positioning of exits within an arbitrarily structure needs to be addressed. Therefore, this paper presents a methodology which combines the use of evacuation models with numerical techniques used in the operational research field, such as Design of Experiments (DoE), Response Surface Models (RSM) and the numerical optimisation techniques. The methodology here presented is restricted to evacuation modelling analysis, nevertheless this same concept can be extended to fire modelling analysis.

Modelling a light-driven phytoplankton succession

We used a numerical model to investigate if and to what extent cellular photoprotective capacity accounts for succession and vertical distribution of marine phytoplankton species/groups. A model describing xanthophyll photoprotective activity in phytoplankton has been implemented in the European Regional Sea Ecosystem Model and applied at the station L4 in the Western English Channel. Primary producers were subdivided into three phytoplankton functional types defined in terms of their capacity to acclimate to different light-specific environments: low light (LL-type), high light (HL-type) and variable light (VL-type) adapted species. The LL-type is assumed to have low cellular level of xanthophyll-cycling pigments (PX) relative to the modelled photosynthetically active pigments (chlorophyll and fucoxanthin (FUCO) = PSP). The HL-type has high PX content relative to PSP while VL-type presents an intermediate PX to PSP ratio. Furthermore, the VL-type is capable of reversibly converting FUCO to PX and synthesizing new PX under high-light stress. In order to reproduce phytoplankton community succession with each of the three groups being dominant in different periods of the year, we had also to assume reduced grazing pressure on HL-adapted species. Model simulations realistically reproduce the observed seasonal patterns of pigments and nutrients highlighting the reasonability of the underpinning assumptions. Our model suggests that pigment-mediated photophysiology plays a primary role in determining the evolution of marine phytoplankton communities in the winter-spring period corresponding to the shoaling of the mixed layer and the increase of light intensity. Grazing selectivity however contributes to the phytoplankton community composition in summer.

Challenges in integrative approaches to modelling the marine ecosystems of the North Atlantic: Physics to fish and coasts to ocean

It has long been recognised that there are strong interactions and feedbacks between climate, upper ocean biogeochemistry and marine food webs, and also that food web structure and phytoplankton community distribution are important determinants of variability in carbon production and export from the euphotic zone. Numerical models provide a vital tool to explore these interactions, given their capability to investigate multiple connected components of the system and the sensitivity to multiple drivers, including potential future conditions. A major driver for ecosystem model development is the demand for quantitative tools to support ecosystem-based management initiatives. The purpose of this paper is to review approaches to the modelling of marine ecosystems with a focus on the North Atlantic Ocean and its adjacent shelf seas, and to highlight the challenges they face and suggest ways forward. We consider the state of the art in simulating oceans and shelf sea physics, planktonic and higher trophic level ecosystems, and look towards building an integrative approach with these existing tools. We note how the different approaches have evolved historically and that many of the previous obstacles to harmonisation may no longer be present. We illustrate this with examples from the on-going and planned modelling effort in the Integrative Modelling Work Package of the EURO-BASIN programme.

Modelling the Stoichiometric Regulation of C-Rich Toxins in Marine Dinoflagellates

Toxin production in marine microalgae was previously shown to be tightly coupled with cellular stoichiometry. The highest values of cellular toxin are in fact mainly associated with a high carbon to nutrient cellular ratio. In particular, the cellular accumulation of C-rich toxins (i.e., with C:N > 6.6) can be stimulated by both N and P deficiency. Dinoflagellates are the main producers of C-rich toxins and may represent a serious threat for human health and the marine ecosystem. As such, the development of a numerical model able to predict how toxin production is stimulated by nutrient supply/deficiency is of primary utility for both scientific and management purposes. In this work we have developed a mechanistic model describing the stoichiometric regulation of C-rich toxins in marine dinoflagellates. To this purpose, a new formulation describing toxin production and fate was embedded in the European Regional Seas Ecosystem Model (ERSEM), here simplified to describe a monospecific batch culture. Toxin production was assumed to be composed by two distinct additive terms; the first is a constant fraction of algal production and is assumed to take place at any physiological conditions. The second term is assumed to be dependent on algal biomass and to be stimulated by internal nutrient deficiency. By using these assumptions, the model reproduced the concentrations and temporal evolution of toxins observed in cultures of Ostreopsis cf. ovata, a benthic/epiphytic dinoflagellate producing C-rich toxins named ovatoxins. The analysis of simulations and their comparison with experimental data provided a conceptual model linking toxin production and nutritional status in this species. The model was also qualitatively validated by using independent literature data, and the results indicate that our formulation can be also used to simulate toxin dynamics in other dinoflagellates. Our model represents an important step towards the simulation and prediction of marine algal toxicity.

Coupled Temperature-Displacement Modelling of Injection Stretch Blow Moulding of PET bottles using Buckley Model

This study presents a fully coupled temperature–displacement finite element modelling of the injection stretch-blow moulding (ISBM) process of polyethylene terephthalate (PET) bottles using ABAQUS with a view to optimising the process conditions. A physically-based material model (Buckley model) was used to predict the mechanical behaviour of PET at temperatures slightly above its glass transition temperature. A model incorporating heat transfer between the stretch rod, the preform and the mould was built using axisymmetric solid elements. Extensive finite element analyses were carried out to predict the deformation, the distribution and history of strain and temperature during ISBM of a 20 g–330 ml bottle, which was made in an in situ test on a Sidel SB06 machine. Comparisons of numerical results with the measurements demonstrate that the model can satisfactorily model the sidewall thickness and material distributions. It is also shown that significant non-linear differentials exist in temperature and strain in both bottle thickness and length directions during the process. This justifies the employment of a volume approach to accurately predict the final mechanical properties of the bottles governed by the orientation and crystallinity which are highly temperature and strain dependent.

Modelling of saline intrusion in marine outfalls

Since their introduction in the 1950s, marine outfalls with diffusers have been prone to saline intrusion, a process in which seawater ingresses into the outfall. This can greatly reduce the dilution and subsequent dispersion of wastewater discharged, sometimes resulting in serious deterioration of coastal water quality. Although long aware of the difficulties posed by saline intrusion, engineers still lack satisfactory methods for its prediction and robust design methods for its alleviation. However, with recent developments in numerical methods and computer power, it has been suggested that commercially available computational fluid dynamics (CFD) software may be a useful aid in combating this phenomenon by improving understanding through synthesising likely behaviour. This document reviews current knowledge on saline intrusion and its implications and then outlines a model-scale investigation of the process undertaken at Queen's University Belfast, using both physical and CFD methods. Results are presented for a simple outfall configuration, incorporating several outlets. The features observed agree with general observations from full-scale marine outfalls, and quantify the intricate internal flow mechanisms associated with saline intrusion. The two-dimensional numerical model was found to represent saline intrusion, but in a qualitative manner, not yet adequate for design purposes. Specific areas requiring further development were identified. The ultimate aim is to provide a reliable, practical and cost effective means by which engineers can minimise saline intrusion through optimised outfall design.

Numerical simulation of separation factor of NaCl in surface force pore flow model

This work deals with modelling and experimental verification of desalination theory (surface force pore flow) . The work has direct application in desalination of sea water.

Numerical analysis of bodyworn UHF antenna systems

Bodyworn antennas are found in a wide range of medical, military and personal communication applications, yet reliable communication from the surface of the human body still presents a range of engineering challenges. At UHF and microwave frequencies, bodyworn antennas can suffer from reduced efficiency due to electromagnetic absorption in tissue, radiation pattern fragmentation and variations in feed-point impedance. The significance and nature of these effects are system specific and depend on the operating frequency, propagation environment and physical constraints on the antenna itself. This paper describes how numerical electromagnetic modelling techniques such as FDTD (finite-difference time-domain) can be used in the design of bodyworn antennas. Examples are presented for 418 MHz, 916 .5 MHz and 2 . 45 GHz, in the context of both biomedical signalling and wireless personal-area networking applications such as the Bluetooth(TM)* wireless technology.

Diagnostic based modelling of radio-frequency driven atmospheric pressure plasmas

Diagnostic based modelling (DBM) actively combines complementary advantages of numerical plasma simulations and relatively simple optical emission spectroscopy (OES). DBM is employed to determine absolute atomic oxygen ground state densities in a helium–oxygen radio-frequency driven atmospheric pressure plasma jet. A comparatively simple one-dimensional simulation yields detailed information on electron properties governing the population dynamics of excited states. Important characteristics of the electron dynamics are found to be largely insensitive to details of the chemical composition and to be in very good agreement with space and phase-resolved OES. Benchmarking the time and space resolved simulation allows us to subsequently derive effective excitation rates as the basis for DBM with simple space and time integrated OES. The population dynamics of the upper O 3p 3P (? = 844 nm) atomic oxygen state is governed by direct electron impact excitation, dissociative excitation, radiation losses and collisional induced quenching. Absolute values for atomic oxygen densities are obtained through tracer comparison with the upper Ar 2p1 (? = 750.4 nm) state. The presented results for the atomic oxygen density show excellent quantitative agreement with independent two-photon laser-induced fluorescence measurements.

Modelling the mechanical response of the Reed-mouthpiece-lip system of a clarinet. Part I. A one-dimensional distributed model

The motion of a clarinet reed that is clamped to a mouthpiece and supported by a lip is simulated in the time-domain using finite difference methods. The reed is modelled as a bar with non-uniform cross section, and is described using a one-dimensional, fourth-order partial differential equation. The interactions with the mouthpiece Jay and the player's lip are taken into account by incorporating conditional contact forces in the bar equation. The model is completed by clamped-free boundary conditions for the reed. An implicit finite difference method is used for discretising the system, and values for the physical parameters are chosen both from laboratory measurements and by accurate tuning of the numerical simulations. The accuracy of the numerical system is assessed through analysis of frequency warping effects and of resonance estimation. Finally, the mechanical properties of the system are studied by analysing its response to external driving forces. In particular, the effects of reed curling are investigated.