Multiphysics modeling and its application to the casting process

A modeling strategy is presented to solve the governing equations of fluid flow, temperature (with solidification), and stress in an integrated manner. These equations are discretized using finite volume methods on unstructured grids, which provide the capability to represent complex domains. Both the cell-centered and vertex-based forms of the finite volume discretization procedure are explained, and the overall integrated solution procedure using these techniques with suitable solvers is detailed. Two industrial processes, based on the casting of metals, are used to demonstrate the capabilities of the resultant modeling framework. This manufacturing process requires a high degree of coupling between the governing physical equations to accurately predict potential defects. Comparisons between model predictions and experimental observations are given.

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.

Dynamic fluid–structure interaction using finite volume unstructured mesh procedures

A three-dimensional finite volume, unstructured mesh (FV-UM) method for dynamic fluid–structure interaction (DFSI) is described. Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, wind response of buildings, flows in elastic pipes and blood vessels. It involves the coupling of fluid flow and structural mechanics, two fields that are conventionally modelled using two dissimilar methods, thus a single comprehensive computational model of both phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply. More recently, strategies for solving the full coupling between the fluid and solid mechanics behaviour have been developed. A key contribution has been made by Farhat et al. [Int. J. Numer. Meth. Fluids 21 (1995) 807] employing FV-UM methods for solving the Euler flow equations and a conventional finite element method for the elastic solid mechanics and the spring based mesh procedure of Batina [AIAA paper 0115, 1989] for mesh movement. In this paper, we describe an approach which broadly exploits the three field strategy described by Farhat for fluid flow, structural dynamics and mesh movement but, in the context of DFSI, contains a number of novel features: • a single mesh covering the entire domain, • a Navier–Stokes flow, • a single FV-UM discretisation approach for both the flow and solid mechanics procedures, • an implicit predictor–corrector version of the Newmark algorithm, • a single code embedding the whole strategy.

A mixed Eulerian-Lagrangian approach for metal extrusion

In this paper a mixed Eulerian-Lagrangian approach for the modelling metal extrusion processes is presented. The approach involves the solution of non-Newtonian fluid flow equations in an Eulerian context, using a free-surface algorithm to track the behaviour of the workpiece and its extrusion. The solid mechanics equations associated with the tools are solved in Lagangrian context. Thermal interactions between the workpiece are modelled and a fluid-structure interaction technique is employed to model the effect of the fluid traction load imposed by the workpiece on the tools. Two extrusion test cases are investigated and the results obtained show the potential of the model with regard to representing the physics of the process and the simulation time.

Modelling electromagnetically levitated liquid droplet oscillations

This work comprises accurate computational analysis of levitated liquid droplet oscillations in AC and DC magnetic fields. The AC magnetic field interacting with the induced electric current within the liquid metal droplet generates intense fluid flow and the coupled free surface oscillations. The pseudo-spectral technique is used to solve the turbulent fluid flow equations for the continuously dynamically transformed axisymmetric fluid volume. The volume electromagnetic force distribution is updated with the shape and position change. We start with the ideal fluid test case for undamped Rayleigh frequency oscillations in the absence of gravity, and then add the viscous and the DC magnetic field damping. The oscillation frequency spectra are further analysed for droplets levitated against gravity in AC and DC magnetic fields at various combinations. In the extreme case electrically poorly conducting, diamagnetic droplet (water) levitation dynamics are simulated. Applications are aimed at pure electromagnetic material processing techniques and the material properties measurements in uncontaminated conditions.

Numerical Modelling for Electronic Packaging - Future Requirements

The fabrication, assembly and testing of electronic packaging can involve complex interactions between physical phenomena such as temperature, fluid flow, electromagnetics, and stress. Numerical modelling and optimisation tools are key computer-aided-engineering technologies that aid design engineers. This paper discusses these technologies and there future developments.

A two dimensional unstructured staggered mesh method with special treatment of tangential velocity

A two dimensional staggered unstructured discretisation scheme for the solution of fluid flow problems has been developed. This scheme stores and solves the velocity vector resolutes normal and parallel to each cell face and other scalar variables (pressure, temperature) are stored at cell centres. The coupled momentum; continuity and energy equations are solved, using the well known pressure correction algorithm SIMPLE. The method is tested for accuracy and convergence behaviour against standard cell-centre solutions in a number of benchmark problems: The Lid-Driven Cavity, Natural Convection in a Cavity and the Melting of Gallium in a rectangular domain.

Numerical simulation of free surface behaviour of a molten liquid metal droplet with and without electromagnetic induction

Electromagnetic levitation of electrically conductive droplets by alternating magnetic fields is a technique used to measure the physical properties of liquid metallic alloys such as surface tension or viscosity. Experiments can be conducted under terrestrial conditions or in microgravity, to reduce electromagnetic stirring and shaping of the droplet. Under such conditions, the time-dependent behaviour of a point of the free surface is recorded. Then the signal is analysed considering the droplet as a harmonic damped oscillator. We use a spectral code, for fluid flow and free surface descriptions, to check the validity of this assumption for two cases. First when the motion inside the droplet is generated by its initial distortion only and second, when the droplet is located in a uniform magnetic field originating far from the droplet. It is found that some deviations exist which can lead to an overestimate of the value of viscosity.

Numerical modelling for electromagnetic processing of materials

Electromagnetic processing of materials (EPM) is one of the most widely practiced and fast growing applications of magnetic and electric forces to fluid flow. EPM is encountered in both industrial processes and laboratory investigations. Applications range in scale from nano-particle manipulation to tonnes of liquid metal treated in the presence of various configurations of magnetic fields. Some of these processes are specifically designed and made possible by the use of the electromagnetic force, like the magnetic levitation of liquid droplets, whilst others involve electric currents essential for electrothermal or electrochemical reasons, for instance, in electrolytic metal production and in induction melting. An insight for the range of established and novel EPM applications can be found in the review presented by Asai [1] in the EPM-2003 conference proceedings.

Dynamic model for metal cleanness evaluation by melting in a cold crucible

Melting of metallic samples in a cold crucible causes inclusions to concentrate on the surface owing to the action of the electromagnetic force in the skin layer. This process is dynamic, involving the melting stage, then quasi-stationary particle separation, and finally the solidification in the cold crucible. The proposed modeling technique is based on the pseudospectral solution method for coupled turbulent fluid flow, thermal and electromagnetic fields within the time varying fluid volume contained by the free surface, and partially the solid crucible wall. The model uses two methods for particle tracking: (1) a direct Lagrangian particle path computation and (2) a drifting concentration model. Lagrangian tracking is implemented for arbitrary unsteady flow. A specific numerical time integration scheme is implemented using implicit advancement that permits relatively large time-steps in the Lagrangian model. The drifting concentration model is based on a local equilibrium drift velocity assumption. Both methods are compared and demonstrated to give qualitatively similar results for stationary flow situations. The particular results presented are obtained for iron alloys. Small size particles of the order of 1 μm are shown to be less prone to separation by electromagnetic field action. In contrast, larger particles, 10 to 100 μm, are easily “trapped” by the electromagnetic field and stay on the sample surface at predetermined locations depending on their size and properties. The model allows optimization for melting power, geometry, and solidification rate.

Three-dimensional simulation of the influence of convection on dendritic solidification

Bulk and interdendritic flow during solidification alters the microstructure development, potentially leading to the formation of defects. In this paper, a 3D numerical model is presented for the simulation of dendritic growth in the presence of fluid flow in both liquid and semi-solid zones during solidification. The dendritic growth was solved by the combination of a stochastic nucleation approach with a finite difference solution of the solute diffusion equation and. a projection method solution of the Navier-Stokes equations. The technique was applied first to simulate the growth of a single dendrite in 2D and 3D in an isothermal environment with forced fluid flow. Significant differences were found in the evolution of dendritic morphology when comparing the 2D and 3D results. In 3D the upstream arm has a faster growth velocity due to easier flow around the perpendicular arms. This also promotes secondary arm formation on the upstream arm. The effect of fluid flow on columnar dendritic growth and micro-segregation in constrained solidification conditions is then simulated. For constrained growth, 2D simulations lead to even greater inaccuracies as compared to 3D.

Thermoelectric effects on alloy solidification microstructure

Thermoelectric currents in the presence of a magnetic field generate Lorentz forces which can drive fluid flow. In the case of dendritic growth a naturally occurring thermoelectric current exists and in the presence of a high magnetic field micro convections are generated. Experimental evidence has attributed changes in microstructure to this effect. A numerical model has been developed to study the flow field around an unconstricted equiaxed dendrite growing under these conditions. The growth is modeled in 2D and 3D by an enthalpy based method and a complex flow structure has been predicted. Using a pseudo-3D approximation for economy, realistic 2D simulations are obtained where a fully coupled transient scheme reveals significant changes to the dendrite morphology reflecting experimental evidence. There is a rotation of the preferred direction of growth and increased secondary branching.

Levitated liquid droplets in AC and DC magnetic field

Electromagnetic levitation of liquid metal droplets can be used to measure the properties of highly reactive liquid materials. Two independent numerical models, the commercial COMSOL and the spectral-collocation based free surface code SPHINX, have been applied to solve the transient electromagnetic, fluid flow and thermodynamic equations, which describe the levitated liquid motion and heating processes. The SPHINX model incorporates free surface deformation to accurately model the oscillations that result from the interaction between the electromagnetic and gravity forces, temperature dependent surface tension, magnetically controlled turbulent momentum transport. The models are adapted to incorporate periodic laser heating at the top of the droplet, which is used to measure the thermal conductivity of the material. Novel effects in the levitated droplet of magnetically damped turbulence and nonlinear growth of velocities in high DC magnetic field are analysed.

Magnetic damping of levitated liquid droplets in AC and DC field

The intense AC magnetic field required to produce levitation in terrestrial conditions, along with the buoyancy and thermo-capillary forces, results in turbulent convective flow within the droplet. The use of a homogenous DC magnetic field allows the convective flow to be damped. However the turbulence properties are affected at the same time, leading to a possibility that the effective turbulent damping is considerably reduced. The MHD modified K-Omega turbulence model allows the investigation of the effect of magnetic field on the turbulence. The model incorporates free surface deformation, the temperature dependent surface tension, turbulent momentum transport, electromagnetic and gravity forces. The model is adapted to incorporate a periodic laser heating at the top of the droplet, which have been used to measure the thermal conductivity of the material by calculating the phase lag between the frequency of the laser heating and the temperature response at the bottom. The numerical simulations show that with the gradual increase of the DC field the fluid flow within the droplet is initially increasing in intensity. Only after a certain threshold magnitude of the field the flow intensity starts to decrease. In order to achieve the flow conditions close to the ‘laminar’ a D.C. magnetic field >4 Tesla is required to measure the thermal conductivity accurately. The reduction in the AC field driven flow in the main body of the drop leads to a noticeable thermo-capillary convection at the edge of the droplet. The uniform vertical DC magnetic field does not stop a translational oscillation of the droplet along the field, which is caused by the variation in total levitation force due to the time-dependent surface deformation.