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.

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

Dynamic melting model for small samples in cold crucible

Purpose – A small size cold crucible offers possibilities for melting various electrically conducting materials with a minimal wall contact. Such small samples can be used for express contamination analysis, preparing limited amounts of reactive alloys or experimental material analyses. Aims to present a model to follow the melting process. Design/methodology/approach – The presents a numerical model in which different types of axisymmetric coil configurations are analysed. Findings – The presented numerical model permits dynamically to follow the melting process, the high-frequency magnetic field distribution change, the free surface and the melting front evolution, and the associated turbulent fluid dynamics. The partially solidified skin on the contact to the cold crucible walls and bottom is dynamically predicted. The segmented crucible shape is either cylindrical, hemispherical or arbitrary shaped. Originality/value – The model presented within the paper permits the analysis of melting times, melt shapes, electrical efficiency and particle tracks.

Models for magnetohydrodynamics of aluminium electrolysis cells

The electric current and the associated magnetic field in aluminium electrolysis cells create effects limiting the cell productivity and possibly cause instabilities: surface waving, ‘anode effects’, erosion of pot lining, feed material sedimentation, etc. The instructive analysis is presented via a step by step inclusion of different physical coupling factors affecting the magnetic field, electric current, velocity and wave development in the electrolysis cells. The full time dependent model couples the nonlinear turbulent fluid dynamics and the extended electromagnetic field in the cell, and the whole bus bar circuit with the ferromagnetic effects. Animated examples for the high amperage cells are presented. The theory and numerical model of the electrolysis cell is extended to the cases of variable cell bottom of aluminium layer and the variable thickness of the electrolyte due to the anode non-uniform burn-out process and the presence of the anode channels. The problem of the channel importance is well known Moreau-Evans model) for the stationary interface and the velocity field, and was validated against measurements in commercial cells, particularly with the recently published ‘benchmark’ test for the MHD models of aluminium cells [1]. The presence of electrolyte channels requires also to reconsider the previous magnetohydrodynamic instability theories and the dynamic wave development models. The results indicate the importance of a ‘sloshing’ parametrically excited MHD wave development in the aluminium production cells.

A numerical investigation of the flow structures and losses for turbulent flow in 90 degrees elbow bends

Computational fluid dynamic modelling was carried out on a series of pipe bends having R/r values of 1.3, 5, and 20, with the purpose of determining the accuracy of numerical models in predicting pressure loss data from which to inform one-dimensional loss models. Four separate turbulence models were studied: the standard k-epsilon model, realizable k-epsilon model, k-omega model, and a Reynolds stress model (RSM). The results are presented for each bend in the form of upstream and downstream pressure profiles, pressure distributions along the inner and outer walls, detailed pressure and velocity fields as well as overall loss values. In each case, measured data were presented to evaluate the predictive ability of each model. The RSM was found to perform the best, producing accurate pressure loss data for bends with R/r values of 5 and 20. For the tightest bend with an R/r value of 1.3, however, predictions were significantly worse due to the presence of flow separation, stronger pressure gradients, and high streamline curvature.

Reynolds Number Effects on Fully-Developed Pulsed Jets Impinging on Flat Surfaces

A systematic study of the effect of the Reynolds number on the fluid dynamics and turbulence statistics of pulsed jets impinging on a flat surface is presented. It has been suggested that the influence of the Reynolds number may be somewhat different for a jet subjected to pulsation when compared to an equivalent steady jet. A comparative study of both steady and pulsating jets is presented for a Reynolds number range from Re = 4;730 to Re = 10;000. All the other factors that affect the flowfield are kept constant, which are H/d = 3, St = 0.25, and d = 30.5 mm. It was found that for the range of the Reynolds numbers tested, pulsation results in a shortening of the jet core, the centerline axial velocity component declines more rapidly, and higher values of the radial velocity component for r/d > 0.75are observed. As the Reynolds number increases, the jet spreads more rapidly, the turbulent kinetic energy and nondimensional turbulent fluctuations decrease, and the flowfield near the impinging surface changes drastically, which is evident with the development of a turbulent momentum exchange interaction away from the wall for r/d > 1.5.

Comparison of Experimental PIV Data and CFD Simulations for Flow in a Diesel Particulate Filter Inlet Diffuser

Flow maldistribution of the exhaust gas entering a Diesel Particulate Filter (DPF) can cause uneven soot distribution during loading and excessive temperature gradients during the regeneration phase. Minimising the magnitude of this maldistribution is therefore an important consideration in the design of the inlet pipe and diffuser, particularly in situations where packaging constraints dictate bends in the inlet pipe close to the filter, or a sharp diffuser angle. This paper describes the use of Particle Image Velocimetry (PIV) to validate a Computational Fluid Dynamic (CFD) model of the flow within the inlet diffuser of a DPF so that CFD can be used with confidence as a tool to minimise this flow maldistribution. PIV is used to study the flow of gas into a DPF over a range of steady state flow conditions. The distribution of flow approaching the front face of the substrate was of particular interest to this study. Optically clear diffusing cones were designed and placed between pipe and substrate to allow PIV analysis to take place. Stereoscopic PIV was used to eliminate any error produced by the optical aberrations caused by looking through the curved wall of the inlet cone. In parallel to the experiments, numerical analysis was carried out using a CFD program with an incorporated DPF model. Boundary conditions for the CFD simulations were taken from the experimental data, allowing an experimental validation of the numerical results. The CFD model incorporated a DPF model, the cement layers seen in segmented filters and the intumescent matting that is commonly used to pack the filter into a metal casing. The mesh contained approximately 580,000 cells and used the realizable ?-e turbulence model. The CFD simulation predicted both pressure drop across the DPF and the velocity field within the cone and at the DPF face with reasonable accuracy, providing confidence in the use the CFD in future work to design new, more efficient cones.

Automatic extraction of Reduced order models from CFD simulations for building energy modeling

Thermal comfort is defined as “that condition of mind which expresses satisfaction with the thermal environment’ [1] [2]. Field studies have been completed in order to establish the governing conditions for thermal comfort [3]. These studies showed that the internal climate of a room was the strongest factor in establishing thermal comfort. Direct manipulation of the internal climate is necessary to retain an acceptable level of thermal comfort. In order for Building Energy Management Systems (BEMS) strategies to be efficiently utilised it is necessary to have the ability to predict the effect that activating a heating/cooling source (radiators, windows and doors) will have on the room. The numerical modelling of the domain can be challenging due to necessity to capture temperature stratification and/or different heat sources (radiators, computers and human beings). Computational Fluid Dynamic (CFD) models are usually utilised for this function because they provide the level of details required. Although they provide the necessary level of accuracy these models tend to be highly computationally expensive especially when transient behaviour needs to be analysed. Consequently they cannot be integrated in BEMS. This paper presents and describes validation of a CFD-ROM method for real-time simulations of building thermal performance. The CFD-ROM method involves the automatic extraction and solution of reduced order models (ROMs) from validated CFD simulations. The test case used in this work is a room of the Environmental Research Institute (ERI) Building at the University College Cork (UCC). ROMs have shown that they are sufficiently accurate with a total error of less than 1% and successfully retain a satisfactory representation of the phenomena modelled. The number of zones in a ROM defines the size and complexity of that ROM. It has been observed that ROMs with a higher number of zones produce more accurate results. As each ROM has a time to solution of less than 20 seconds they can be integrated into the BEMS of a building which opens the potential to real time physics based building energy modelling.