A generic time and space accurate numerical approach to closely coupled fluid-structure interaction problems

Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, flow in elastic pipes and blood vessels and extrusion of metals through dies. However a comprehensive computational model of these multi-physics phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply even to the extent in metal forming, for example, that the deformation of the die is totally ignored. More recently, strategies for solving the full coupling between the fluid and soild mechanics behaviour have developed. Conventionally, the computational modelling of fluid structure interaction is problematical since computational fluid dynamics (CFD) is solved using finite volume (FV) methods and computational structural mechanics (CSM) is based entirely on finite element (FE) methods. In the past the concurrent, but rather disparate, development paths for the finite element and finite volume methods have resulted in numerical software tools for CFD and CSM that are different in almost every respect. Hence, progress is frustrated in modelling the emerging multi-physics problem of fluid structure interaction in a consistent manner. Unless the fluid-structure coupling is either one way, very weak or both, transferring and filtering data from one mesh and solution procedure to another may lead to significant problems in computational convergence. Using a novel three phase technique the full interaction between the fluid and the dynamic structural response are represented. The procedure is demonstrated on some challenging applications in complex three dimensional geometries involving aircraft flutter, metal forming and blood flow in arteries.

Multi-physics-multi-scale simulation and optimisation: the next generation

FEA and CFD analysis is becoming ever more complex with an emerging demand for simulation software technologies that can address ranges of problems that involve combinations of interactions amongst varying physical phenomena over a variety of time and length scales. Computation modelling of such problems requires software technologies that enable the representation of these complex suites of 'physical' interactions. This functionality requires the structuring of simulation modules for specific physical phemonmena so that the coupling can be effectiely represented. These 'multi-physics' and 'multi-scale' computations are very compute intensive and so the simulation software must operate effectively in parallel if it is to be used in this context. Of course the objective of 'multi-physics' and 'multi-scale' simulation is the optimal design of engineered systems so optimistation is an important feature of such classes of simulation. In this presentation, a multi-disciplinary approach to simulation based optimisation is described with some key examples of application to challenging engineering problems.

Combined vertex-based - cell-centred finite volume method for flows in complex geometries

CFD modelling of 'real-life' processes often requires solutions in complex three dimensional geometries, which can often result in meshes where aspects of it are badly distorted. Cell-centred finite volume methods, typical of most commercial CFD tools, are computationally efficient, but can lead to convergence problems on meshes which feature cells with high non-orthogonal shapes. The vertex-based finite volume method handles distorted meshes with relative ease, but is computationally expensive. A combined vertex-based - cell-centred (VB-CC) technique, detailed in this paper, allows solutions on distorted meshes that defeat purely cell-centred physical models to be employed in the solution of other transported quantities. The VB-CC method is validated with benchmark solutions for thermally driven flow and turbulent flow. An early application of this hybrid technique is to three-dimensional flow over an aircraft wing, although it is planned to use it in a wide variety of processing applications in the future.

Experimental and numerical study of the cold crucible melting process

The cold crucible, or induction skull melting process as is otherwise known, has the potential to produce high purity melts of a range of difficult to melt materials, including Ti–Al and Ti6Al4V alloys for Aerospace, Ti–Ta and other biocompatible materials for surgical implants, silicon for photovoltaic and electronic applications, etc. A water cooled AC coil surrounds the crucible causing induction currents to melt the alloy and partially suspend it against gravity away from water-cooled surfaces. Strong stirring takes place in the melt due to the induced electromagnetic Lorentz forces and very high temperatures are attainable under the right conditions (i.e., provided contact with water cooled walls is minimised). In a joint numerical and experimental research programme, various aspects of the design and operation of this process are investigated to increase our understanding of the physical mechanisms involved and to maximise process efficiency. A combination of FV and Spectral CFD techniques are used at Greenwich to tackle this problem numerically, with the experimental work taking place at Birmingham University. Results of this study, presented here, highlight the influence of turbulence and free surface behaviour on attained superheat and also discuss coil design variations and dual frequency options that may lead to winning crucible designs.

Computational modelling of bubbles, droplets and particles in metals reduction and refining

The CFD modelling of metals reduction processes particularly always seems to involve the interaction of liquid metals, a gas (often air) top space, liquid droplets in the top space and injection of both solid particles and gaseous bubbles into the bath. These phases all interact and exhange mass, momentum and energy. Often it is the extent to which these multi-phase phemomena can be effectively captured within the CFD model which determines whether or not a tool of genuine use to the target industry sector can constructed. In this paper we discuss these issues in the context of two problems - one involving the injection of sparging gases into a steel continuous caster and the other based on the development of a novel process for aluminium electrolysis.

A comprehensive copper stockpile leach model: background and model formulation

The first phase in the sign, development and implementation of a comprehensive computational model of a copper stockpile leach process is presented. The model accounts for transport phenomena through the stockpile, reaction kinetics for the important mineral species, oxgen and bacterial effects on the leach reactions, plus heat, energy and acid balances for the overall leach process. The paper describes the formulation of the leach process model and its implementation in PHYSICA+, a computational fluid dynamic (CFD) software environment. The model draws on a number of phenomena to represent the competing physical and chemical features active in the process model. The phenomena are essentially represented by a three-phased (solid liquid gas) multi-component transport system; novel algorithms and procedures are required to solve the model equations, including a methodology for dealing with multiple chemical species with different reaction rates in ore represented by multiple particle size fractions. Some initial validation results and application simulations are shown to illustrate the potential of the model.

Computational modelling of freeze layers in smelting processes

The use of computational modelling in examining process engineering issues is very powerful. It has been used in the development of the HIsmelt process from its concept. It is desirable to further water-cool the HIsmelt vessel to reduce downtime for replacing refractory. Water-cooled elements close to a metal bath run the risk of failure. This generally occurs when a process perturbation causes the freeze and refractory layers to come away from the water-cooled element, which is then exposed to liquid metal. The element fails as they are unable to remove all the heat. Modelling of the water-cooled element involves modelling the heat transfer, fluid flow, stress and solidification for a localised section of the reaction vessel. The complex interaction between the liquid slag and the refractory applied to the outside of thewater-cooled element is also being examined to model the wear of this layer. The model is being constructed in Physica, a CFD code developed at the University of Greenwich. Modelling of this system has commenced with modelling solidification test cases. These test cases have been used to validate the CFD code’s capability to model the solidification in this system. A model to track the penetration of slag into refractory has also been developed and tested.

Turbulence modelling for thermal management of electronic systems

The electronics industry is developing rapidly together with the increasingly complex problem of microelectronic equipment cooling. It has now become necessary for thermal design engineers to consider the problem of equipment cooling at some level. The use of Computational Fluid Dynamics (CFD) for such investigations is fast becoming a powerful and almost essential tool for the design, development and optimisation of engineering applications. However turbulence models remain a key issue when tackling such flow phenomena. The reliability of CFD analysis depends heavily on the turbulence model employed together with the wall functions implemented. In order to resolve the abrupt fluctuations experienced by the turbulent energy and other parameters located at near wall regions and shear layers a particularly fine computational mesh is necessary which inevitably increases the computer storage and run-time requirements. This paper will discuss results from an investigation into the accuract of currently used turbulence models. Also a newly formulated transitional hybrid turbulence model will be introduced with comparisonsaagainst experimental data.

Low Reynolds number turbulence models for accurate thermal simulations of electronic components

The electronics industry and the problems associated with the cooling of microelectronic equipment are developing rapidly. Thermal engineers now find it necessary to consider the complex area of equipment cooling at some level. This continually growing industry also faces heightened pressure from consumers to provide electronic product miniaturization, which in itself increases the demand for accurate thermal management predictions to assure product reliability. Computational fluid dynamics (CFD) is considered a powerful and almost essential tool for the design, development and optimization of engineering applications. CFD is now widely used within the electronics packaging design community to thermally characterize the performance of both the electronic component and system environment. This paper discusses CFD results for a large variety of investigated turbulence models. Comparison against experimental data illustrates the predictive accuracy of currently used models and highlights the growing demand for greater mathematical modelling accuracy with regards to thermal characterization. Also a newly formulated low Reynolds number (i.e. transitional) turbulence model is proposed with emphasis on hybrid techniques.

Response Surface Modeling and Optimisation for Reliable Electronic Products

The aim of integrating computational mechanics (FEA and CFD) and optimization tools is to speed up dramatically the design process in different application areas concerning reliability in electronic packaging. Design engineers in the electronics manufacturing sector may use these tools to predict key design parameters and configurations (i.e. material properties, product dimensions, design at PCB level. etc) that will guarantee the required product performance. In this paper a modeling strategy coupling computational mechanics techniques with numerical optimization is presented and demonstrated with two problems. The integrated modeling framework is obtained by coupling the multi-physics analysis tool PHYSICA - with the numerical optimization package - Visua/DOC into a fuJly automated design tool for applications in electronic packaging. Thermo-mechanical simulations of solder creep deformations are presented to predict flip-chip reliability and life-time under thermal cycling. Also a thermal management design based on multi-physics analysis with coupled thermal-flow-stress modeling is discussed. The Response Surface Modeling Approach in conjunction with Design of Experiments statistical tools is demonstrated and used subsequently by the numerical optimization techniques as a part of this modeling framework. Predictions for reliable electronic assemblies are achieved in an efficient and systematic manner.

Simulation of ship evacuation and passenger circulation

This paper describes work carried out in the FIRE EXIT research project. FIRE EXIT aims to develop an Evacuation Simulator, capable of addressing issues of mustering, ship motions, fire and abandonment. In achieving these aims, FIRE EXIT took as its starting point the state-of-the-art in ship evacuation simulation (the maritimeEXODUS software), fire simulation (the SMARTFIRE software) and large-scale experimental facilities (the SHEBA facility). It then significantly enhanced these capabilities. A number of new technologies have been developed in achieving these objectives. The innovations include directly linking CFD fire simulation with evacuation and abandonment software and automatic data transfer from concept design software allowing rapid generation of ship simulation models. Software usability was augmented by a module for interpretation of evacuation software output. Enhancements to a ship evacuation testing rig have resulted in a unique facility, capable of providing passenger movement data for realistic evacuation scenarios and large scale tests have provided meaningful data for the evacuation simulation.

Investigation into the performance of turbulence models for fluid flow and heat transfer phenomena in electronic applications

Computational Fluid Dynamics (CFD) is gradually becoming a powerful and almost essential tool for the design, development and optimization of engineering applications. However the mathematical modelling of the erratic turbulent motion remains the key issue when tackling such flow phenomena. The reliability of CFD analysis depends heavily on the turbulence model employed together with the wall functions implemented. In order to resolve the abrupt changes in the turbulent energy and other parameters situated at near wall regions a particularly fine mesh is necessary which inevitably increases the computer storage and run-time requirements. Turbulence modelling can be considered to be one of the three key elements in CFD. Precise mathematical theories have evolved for the other two key elements, grid generation and algorithm development. The principal objective of turbulence modelling is to enhance computational procedures of efficient accuracy to reproduce the main structures of three dimensional fluid flows. The flow within an electronic system can be characterized as being in a transitional state due to the low velocities and relatively small dimensions encountered. This paper presents simulated CFD results for an investigation into the predictive capability of turbulence models when considering both fluid flow and heat transfer phenomena. Also a new two-layer hybrid kε / kl turbulence model for electronic application areas will be presented which holds the advantages of being cheap in terms of the computational mesh required and is also economical with regards to run-time.

Modelling the wave soldering process

Compuational fluid dynamics (CFD) is used to help understand the gas flow characteristics in the wave soldering process. CFD has the ability to calculate (1) heal transfer, (2) fluid dynamics, and (3) oxygen concentration throughout the wave soldering machine. Understanding the impact of fluid dynamics on oxygen concentration is important as excessive oxygen at the solder bath can lead to high dross contents and hence poor solder joint quality on the printed circuit board. This paper describes the CFD modelling approach and illustrates its capability for a machine which has nitrogen injectors near the solder bath. Different magnitiutes of nitrogen flow rates are investigated and it is demonstrated how these effect the oxygen concentration at the bath surface.

Pseudo-spectral solutions for fluid flow and heat transfer in electro-metallurgical applications

The pseudo-spectral solution method offers a flexible and fast alternative to the more usual finite element/volume/difference methods, particularly when the long-time transient behaviour of a system is of interest. Since the exact solution is obtained at the grid collocation points superior accuracy can be achieved on modest grid resolution. Furthermore, the grid can be freely adapted with time and in space, to particular flow conditions or geometric variations. This is especially advantageous where strongly coupled, time-dependent, multi-physics solutions are investigated. Examples include metallurgical applications involving the interaction of electromagnetic fields and conducting liquids with a free sutface. The electromagnetic field then determines the instantaneous liquid volume shape and the liquid shape affects in turn the electromagnetic field. In AC applications a thin "skin effect" region results on the free surface that dominates grid requirements. Infinitesimally thin boundary cells can be introduced using Chebyshev polynomial expansions without detriment to the numerical accuracy. This paper presents a general methodology of the pseudo-spectral approach and outlines the solution procedures used. Several instructive example applications are given: the aluminium electrolysis MHD problem, induction melting and stirring and the dynamics of magnetically levitated droplets in AC and DC fields. Comparisons to available analytical solutions and to experimental measurements will be discussed.

Investigation into the performance of turbulence models for fluid flow and heat transfer phenomena in electronic applications

Heat is extracted away from an electronic package by convection, conduction, and/or radiation. The amount of heat extracted by forced convection using air is highly dependent on the characteristics of the airflow around the package which includes its velocity and direction. Turbulence in the air is also important and is required to be modeled accurately in thermal design codes that use computational fluid dynamics (CFD). During air cooling the flow can be classified as laminar, transitional, or turbulent. In electronics systems, the flow around the packages is usually in the transition region, which lies between laminar and turbulent flow. This requires a low-Reynolds number numerical model to fully capture the impact of turbulence on the fluid flow calculations. This paper provides comparisons between a number of turbulence models with experimental data. These models included the distance from the nearest wall and the local velocity (LVEL), Wolfshtein, Norris and Reynolds, k-ε, k-ω, shear-stress transport (SST), and kε/kl models. Results show that in terms of the fluid flow calculations most of the models capture the difficult wake recirculation region behind the package reasonably well, although for packages whose heights cause a high degree of recirculation behind the package the SST model appears to struggle. The paper also demonstrates the sensitivity of the models to changes in the mesh density; this study is aimed specifically at thermal design engineers as mesh independent simulations are rarely conducted in an industrial environment.