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.

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.

Unstructured staggered mesh methods for fluid flow, heat transfer and phase change

Unstructured grid meshes used in most commercial CFD codes inevitably adopt collocated variable solution schemes. These schemes have several shortcomings, mainly due to the interpolation of the pressure gradient, that lead to slow convergence. In this publication we show how it is possible to use a much more stable staggered mesh arrangement in an unstructured code. Several alternative groupings of variables are investigated in a search for the optimum scheme.

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 and volume methods, particularly when the long-time transient behaviour of a system is of interest. The exact solution is obtained at grid collocation points leading to superior accuracy on modest grids. Furthermore, the grid can be freely adapted in time and space to particular flow conditions or geometric variations, especially useful 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 surface. The electromagnetic field determines the instantaneous liquid volume shape, which then affects the electromagnetic field. A general methodology of the pseudo-spectral approach is presented, with several instructive example applications: the aluminium electrolysis MHD problem, induction melting in a cold crucible and the dynamics of AC/DC magnetically levitated droplets. Finally, comparisons with available analytical solutions and to experimental measurements are discussed.

AC and DC magnetic levitation: melting, fluid flow and oscillations

The presented numerical modelling for the magnetic levitation involves coupling of the electromagnetic field, liquid shape change, fluid velocities and the temperature field at every time step during the simulation in time evolution. Combination of the AC and DC magnetic fields can be used to achieve high temperature, stable levitation conditions. The oscillation frequency spectra are analysed for droplets levitated in AC and DC magnetic fields at various combinations. An electrically poorly conducting, diamagnetic droplet (e.g. water) can be stably levitated using the dia- and para-magnetic properties of the sample material in a high intensity, gradient DC field.

A computational model for defect prediction in shape castings based on the interaction of free surface flow, heat transfer, and solidification phenomena

High-integrity castings require sophisticated design and manufacturing procedures to ensure they are essentially macrodefect free. Unfortunately, an important class of such defects—macroporosity, misruns, and pipe shrinkage—are all functions of the interactions of free surface flow, heat transfer, and solidication in complex geometries. Because these defects arise as an interaction of the preceding continuum phenomena, genuinely predictive models of these defects must represent these interactions explicitly. This work describes an attempt to model the formation of macrodefects explicitly as a function of the interacting continuum phenomena in arbitrarily complex three-dimensional geometries. The computational approach exploits a compatible set of finite volume procedures extended to unstructured meshes. The implementation of the model is described together with its testing and a measure of validation. The model demonstrates the potential to predict reliably shrinkage macroporosity, misruns, and pipe shrinkage directly as a result of interactions among free-surface fluid flow, heat transfer, and solidification.

Assessing the scalability of multiphysics tools for modeling solidification and melting processes on parallel clusters

A comprehensive solution of solidification/melting processes requires the simultaneous representation of free surface fluid flow, heat transfer, phase change, nonlinear solid mechanics and, possibly, electromagnetics together with their interactions, in what is now known as multiphysics simulation. Such simulations are computationally intensive and the implementation of solution strategies for multiphysics calculations must embed their effective parallelization. For some years, together with our collaborators, we have been involved in the development of numerical software tools for multiphysics modeling on parallel cluster systems. This research has involved a combination of algorithmic procedures, parallel strategies and tools, plus the design of a computational modeling software environment and its deployment in a range of real world applications. One output from this research is the three-dimensional parallel multiphysics code, PHYSICA. In this paper we report on an assessment of its parallel scalability on a range of increasingly complex models drawn from actual industrial problems, on three contemporary parallel cluster systems.

Modelling the fluid dynamics and coupled phenomena in arc weld pools

A 3D model of melt pool created by a moving arc type heat sources has been developed. The model solves the equations of turbulent fluid flow, heat transfer and electromagnetic field to demonstrate the flow behaviour phase-change in the pool. The coupled effects of buoyancy, capillary (Marangoni) and electromagnetic (Lorentz) forces are included within an unstructured finite volume mesh environment. The movement of the welding arc along the workpiece is accomplished via a moving co-ordinator system. Additionally a method enabling movement of the weld pool surface by fluid convection is presented whereby the mesh in the liquid region is allowed to move through a free surface. The surface grid lines move to restore equilibrium at the end of each computational time step and interior grid points then adjust following the solution of a Laplace equation.

The development and validation of computational MHD techniques for the modelling of metal production processes

Magnetic fields are used in a number of processes related to the extraction of metals, production of alloys and the shaping of metal components. Computational techniques have an increasingly important role to play in the simulation of such processes, since it is often difficult or very costly to conduct experiments in the high temperature conditions encountered and the complex interaction of fluid flow, heat transfer and magnetic fields means simple analytic models are often far removed from reality. In this paper an overview of the computational activity at the University of Greenwich is given in this area, covering the past ten years. The overview is given from the point of view of the modeller and within the space limitations imposed by the format it covers the numerical methods used, attempts at validation against experiments or analytic procedures; it highlights successes, but also some failures. A broad range of models is covered in the review (and accompanying lecture), used to simulate (a) A-C field applications: induction melting, magnetic confinement and levitation, casting and (b) D-C field applications such as: arc welding and aluminium electroloysis. Most of these processes involve phase change of the metal (melting or solidification), the presence of a dynamic free surface and turbulent flow. These issues affect accuracy and need to be address by the modeller.

Computational modelling of multi-physics processes on high performance parallel computer systems

The demands of the process of engineering design, particularly for structural integrity, have exploited computational modelling techniques and software tools for decades. Frequently, the shape of structural components or assemblies is determined to optimise the flow distribution or heat transfer characteristics, and to ensure that the structural performance in service is adequate. From the perspective of computational modelling these activities are typically separated into: • fluid flow and the associated heat transfer analysis (possibly with chemical reactions), based upon Computational Fluid Dynamics (CFD) technology • structural analysis again possibly with heat transfer, based upon finite element analysis (FEA) techniques.

Computational modelling of metal extrusion and forging processes

The computational modelling of extrusion and forging processes is now well established. There are two main approaches: Lagrangian and Eulerian. The first has considerable complexities associated with remeshing, especially when the code is parallelised. The second approach means that the mould has to be assumed to be entirely rigid and this may not be the case. In this paper, a novel approach is described which utilises finite volume methods on unstructured meshes. This approach involves the solution of free surface non-Newtonian fluid flow equations in an Eulerian context to track the behaviour of the workpiece and its extrusion/forging, and the solution of the solid mechanics equations in the Lagrangian context to predict the deformation/stress behaviour of the die. Test cases for modelling extrusion and forging problems using this approach will be presented.

Computational fluid dynamics: advancements in technology for modeling iron and steelmaking processes

Computational fluid dynamics (CFD) software technology has formed the basis of many investigations into the behavior and optimization of primary iron and steelmaking processes for the last 25+ years. The objective of this contribution is to review the progress in CFD technologies over the last decade or so and how this can be brought to bear in advancing the process analysis capability of primary ferrous operations. In particular, progress on key challenges such as compute performance, fluid-structure transformation and interaction, and increasingly complex geometries are highlighted.

The study of flow and temperature fields in conducting droplets suspended in a DC/AC combination field

Electromagnetic Levitation (EML) is a valuable method for measuring the thermo-physical properties of metals - surface tensions, viscosity, thermal/electrical conductivity, specific heat, hemispherical emissivity, etc. – beyond their melting temperature. In EML, a small amount of the test specimen is melted by Joule heating in a suspended AC coil. Once in liquid state, a small perturbation causes the liquid envelope to oscillate and the frequency of oscillation is then used to compute its surface tension by the well know Rayleigh formula. Similarly, the rate at which the oscillation is dampened relates to the viscosity. To measure thermal conductivity, a sinusoidally varying laser source may be used to heat the polar axis of the droplet and the temperature response measured at the polar opposite – the resulting phase shift yields thermal conductivity. All these theoretical methods assume that convective effects due to flow within the droplet are negligible compared to conduction, and similarly that the flow conditions are laminar; a situation that can only be realised under microgravity conditions. Hence the EML experiment is the method favoured for Spacelab experiments (viz. TEMPUS). Under terrestrial conditions, the full gravity force has to be countered by a much larger induced magnetic field. The magnetic field generates strong flow within the droplet, which for droplets of practical size becomes irrotational and turbulent. At the same time the droplet oscillation envelope is no longer ellipsoidal. Both these conditions invalidate simple theoretical models and prevent widespread EML use in terrestrial laboratories. The authors have shown in earlier publications that it is possible to suppress most of the turbulent convection generated in the droplet skin layer, through use of a static magnetic field. Using a pseudo-spectral discretisation method it is possible compute very accurately the dynamic variation in the suspended fluid envelope and simultaneously compute the time-varying electromagnetic, flow and thermal fields. The use of a DC field as a dampening agent was also demonstrated in cold crucible melting, where suppression of turbulence was achieved in a much larger liquid metal volume and led to increased superheat in the melt and reduction of heat losses to the water-cooled walls. In this paper, the authors describe the pseudo-spectral technique as applied to EML to compute the combined effects of AC and DC fields, accounting for all the flow-induced forces acting on the liquid volume (Lorentz, Maragoni, surface tension, gravity) and show example simulations.

The numerical modelling of multi-physics phenomena: The shape casting process

In this paper a computer simulation tool capable of modelling multi-physics processes in complex geometry has been developed and applied to the casting process. The quest for high-quality complex casting components demanded by the aerospace and automobile industries, requires more precise numerical modelling techniques and one that need to be generic and modular in its approach to modelling multi-processes problems. For such a computer model to be successful in shape casting, the complete casting process needs to be addressed, the major events being:-• Filling of hot liquid metal into a cavity mould • Solidification and latent heat evolution of liquid metal • Convection currents generated in liquid metal by thermal gradients • Deformation of cast and stress development in solidified metal • Macroscopic porosity formation The above phenomena combines the analysis of fluid flow, heat transfer, change of phase and thermal stress development. None of these events can be treated in isolation as they inexorably interact with each other in a complex way. Also conditions such as design of running system, location of feeders and chills, moulding materials and types of boundary conditions can all affect on the final cast quality and must be appropriately represented in the model.