An integrated modeling approach to solder joint formation

The attachment of electronic components to printed circuit boards using solder material is a complex process. This paper presents a novel modeling methodology, which integrates the governing physics taking place. Multiphysics modeling technology, imbedded into the simulation tool—PHYSICA is used to simulate fluid flow, heat transfer, solidification, and stress evolution in an integrated manner. Results using this code are presented, detailing the mechanical response of two solder materials as they cool, solidify and then deform. The shape that a solder joint takes upon melting is predicted using the SURFACE EVOLVER code. Details are given on how these predictions can be used in the PHYSICA code to provide a modeling route by which the shape, solidification history, and resulting stress profiles can be predicted.

Discretisation procedures for multi-physics phenomena

Procedures are described for solving the equations governing a multi-physics process. Finite volume techniques are used to discretise, using the same unstructured mesh, the equations of fluid flow, heat transfer with solidification, and solid deformation. These discretised equations are then solved in an integrated manner. The computational mechanics environment, PHYSICA, which facilitates the building of multi-physics models, is described. Comparisons between model predictions and experimental data are presented for the casting of metal components.

Modeling the dynamics of Magnetic Semilevitation Melting

In semilevitation melting, a cylindrical metal ingot is melted by a coaxial a.c. induction coil. A watercooled solid base supports the ingot, while the top and side free surface is confined by the magnetic forces as the melting front progresses. The dynamic interplay between gravity, hydrodynamic stress, and the Lorentz force in the fluid determines the instantaneous free surface shape. The coupled nonstationary equations for turbulent flow, heat with phase change, and high-frequency electromagnetic field are solved numerically for the axisymmetric time-dependent domain by a continuous mesh transformation, using a pseudospectral method. Results are obtained for the two actually existing coil configurations and several validation cases.

Multiphysics modelling of the metals casting process

Metals casting is a process governed by the interaction of a range of physical phenomena. Most computational models of this process address only what are conventionally regarded as the primary phenomena-heat conduction and solidification. However, to predict the formation of porosity (a factor of key importance in cast quality) requires the modelling of the interaction of the fluid flow, heat transfer, solidification and the development of stress-deformation in the solidified part of a component. In this paper, a model of the casting process is described which addresses all the main continuum phenomena involved in a coupled manner. The model is solved numerically using novel finite volume unstructured mesh techniques, and then applied to both the prediction of shape deformation (plus the subsequent formation of a gap at the metal-mould interface and its impact on the heat transfer behaviour) and porosity formation in solidifying metal components. Although the porosity prediction model is phenomenologically simplistic it is based on the interaction of the continuum phenomena and yields good comparisons with available experimental results. This work represents the first of the next generation of casting simulation tools to predict aspects of the structure of cast components.

Parallel performance in multi-physics simulation

A comprehensive simulation of solidification/melting processes requires the simultaneous representation of free surface fluid flow, heat transfer, phase change, non-linear solid mechanics and, possibly, electromagnetics together with their interactions in what is now referred to as "multi-physics" simulation. A 3D computational procedure and software tool, PHYSICA, embedding the above multi-physics models using finite volume methods on unstructured meshes (FV-UM) has been developed. Multi-physics simulations are extremely compute intensive and a strategy to parallelise such codes has, therefore, been developed. This strategy has been applied to PHYSICA and evaluated on a range of challenging multi-physics problems drawn from actual industrial cases.

Modeling of ingot development during the start-up phase of direct chill casting

Direct chill (DC) casting is a core primary process in the production of aluminum ingots. However, its operational optimization is still under investigation with regard to a number of features, one of which is the issue of curvature at the base of the ingot. Analysis of these features requires a computational model of the process that accounts for the fluid flow, heat transfer, solidification phase change, and thermomechanical analysis. This article describes an integrated approach to the modeling of all the preceding phenomena and their interactions.

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.

Maximising heat transfer efficiency in the cold crucible induction melting process

Induction heating is an efficient method used to melt electrically conductive materials, particularly if melting takes place in a ceramic crucible. This form of melting is particularly good for alloys, as electromagnetic forces set up by the induction coil lead to vigorous stirring of the melt ensuring homogeneity and uniformity in temperature. However, for certain reactive alloys, or where high purity is required, ceramic crucibles cannot be used, but a water-cooled segmented copper crucible is employed instead. Water cooling prevents meltdown or distortion of the metal wall, but much of the energy goes into the coolant. To reduce this loss, the electromagnetic force generated by the coil is used to push the melt away from the walls and so minimise contact with water-cooled surfaces. Even then, heat is lost through the crucible base where contact is inevitable. In a collaborative programme between Greenwich and Birmingham Universities, computer modelling has been used in conjunction with experiments to improve the superheat attainable in the melt for a,number of alloys, especially for y-TiAl intermetallics to cast aeroengine turbine blades. The model solves the discretised form of the turbulent Navier-Stokes, thermal energy conservation and Maxwell equations using a Spectral Collocation technique. The time-varying melt envelope is followed explicitly during the computation using an adaptive mesh. This paper briefly describes the mathematical model used to represent the interaction between the magnetic field, fluid flow, heat transfer and change of phase in the crucible and identifies the proportions of energy used in the melt, lost in the crucible base and in the crucible walls. The role of turbulence is highlighted as important in controlling heat losses and turbulence damping is introduced as a means of improving superheat. Model validation is against experimental results and shows good agreement with measured temperatures and energy losses in the cooling fluid throughout the melting cycle.

Heat transfer coefficients from Newtonian and non-Newtonian fluids flowing in laminar regime in a helical coil

This study aimed to carry out experimental work to obtain, for Newtonian and non-Newtonian fluids, heat transfer coefficients, at constant wall temperature as boundary condition, in fully developed laminar flow inside a helical coil. The Newtonian fluids studied were aqueous solutions of glycerol, 25%, 36%, 43%, 59% and 78% (w/w) and the non-Newtonian fluids aqueous solutions of carboxymethylcellulose (CMC), a polymer, with concentrations 0.1%, 0.2%, 0.3%, 0.4% and 0.6% (w/w) and aqueous solutions of xanthan gum (XG), another polymer, with concentrations 0.1% and 0.2% (w/w). According to the rheological study performed, the polymer solutions had shear thinning behavior and different values of elasticity. The helical coil used has internal diameter, curvature ratio, length and pitch, respectively: 0.004575 m, 0.0263, 5.0 m and 11.34 mm. The Nusselt numbers for the CMC solutions are, on average, slightly higher than those for Newtonian fluids, for identical Prandtl and generalized Dean numbers. As outcome, the viscous component of the shear thinning polymer tends to potentiate the mixing effect of the Dean cells. The Nusselt numbers of the XG solutions are significant lower than those of the Newtonian solutions, for identical Prandtl and generalized Dean numbers. Therefore, the elastic component of the polymer tends to diminish the mixing effect of the Dean cells. A global correlation, for Nusselt number as a function of Péclet, generalized Dean and Weissenberg numbers for all Newtonian and non-Newtonian solutions studied, is presented.

Interactions between the physical soil environment and a horizontal ground coupled heat pump, for a domestic site in the UK

There is currently an increased interest of Government and Industry in the UK, as well as at the European Community level and International Agencies (i.e. Department of Energy, American International Energy Agency), to improve the performance and uptake of Ground Coupled Heat Pumps (GCHP), in order to meet the 2020 renewable energy target. A sound knowledge base is required to help inform the Government Agencies and advisory bodies; detailed site studies providing reliable data for model verification have an important role to play in this. In this study we summarise the effect of heat extraction by a horizontal ground heat exchanger (installed at 1 m depth) on the soil physical environment (between 0 and 1 m depth) for a site in the south of the UK. Our results show that the slinky influences the surrounding soil by significantly decreasing soil temperatures. Furthermore, soil moisture contents were lower for the GCHP soil profile, most likely due to temperature-gradient related soil moisture migration effects and a decreased hydraulic conductivity, the latter as a result of increased viscosity (caused by the lower temperatures for the GCHP soil profile). The effects also caused considerable differences in soil thermal properties. This is the first detailed mechanistic study conducted in the UK with the aim to understand the interactions between the soil, horizontal heat exchangers and the aboveground environment. An increased understanding of these interactions will help to achieve an optimum and sustainable use of the soil heat resources in the future. The results of this study will help to calibrate and verify a simulation model that will provide UK-wide recommendations to improve future GCHP uptake and performance, while safeguarding the soil physical resources.

Simulation of a domestic ground source heat pump system using a three-dimensional numerical borehole heat exchanger model

Common approaches to the simulation of borehole heat exchangers (BHEs) assume heat transfer in circulating fluid and grout to be in a quasi-steady state and ignore fluctuations in fluid temperature due to transport of the fluid around the loop. However, in domestic ground source heat pump (GSHP) systems, the heat pump and circulating pumps switch on and off during a given hour; therefore, the effect of the thermal mass of the circulating fluid and the dynamics of fluid transport through the loop has important implications for system design. This may also be important in commercial systems that are used intermittently. This article presents transient simulation of a domestic GSHP system with a single BHE using a dynamic three-dimensional (3D) numerical BHE model. The results show that delayed response associated with the transit of fluid along the pipe loop is of some significance in moderating swings in temperature during heat pump operation. In addition, when 3D effects are considered, a lower heat transfer rate is predicted during steady operations. These effects could be important when considering heat exchanger design and system control. The results will be used to develop refined two-dimensional models.

Photovoltaic and solar-assisted ground-source heat pump systems

This review investigates the performance of photovoltaic and solar-assisted ground-source heat pumps in which solar heat is transferred to the ground to improve the coefficient of performance. A number of studies indicate that, for systems with adequately sized ground heat exchangers, the effect on system efficiency is small: about 1% improvement if the heat source is photovoltaic, a 1–2% decline if the source is solar thermal. With possible exceptions for systems in which the ground heat exchanger is undersized, or natural recharge from ground water is insufficient, solar thermal energy is better used for domestic hot water than to recharge ground heat. This appears particularly true outside the heating season, as although much of the heat extracted from the ground can be replaced, it seems to have little effect on the coefficient of performance. Any savings in electrical consumption that do result from an improved coefficient can easily be outweighed by an inefficient control system for the circulation pumps.

Identifiering av parametrar till ackumulatortanken Solus1050 från consolar - en metodbeskrivning

This report describes the work done creating a computer model of a kombi tank from Consolar. The model was created with Presim/Trnsys and Fittrn and DF were used to identify the parameters. Measurements were carried out and were used to identify the values of the parameters in the model. The identifications were first done for every circuit separately. After that, all parameters are normally identified together using all the measurements. Finally the model should be compared with other measurements, preferable realistic ones. The two last steps have not yet been carried out, because of problems finding a good model for the domestic hot water circuit.The model of the domestic hot water circuit give relatively good results for low flows at 5 l/min, but is not good for higher flows. In the report suggestions for improving the model are given. However, there was not enough time to test this within the project as much time was spent trying to solve problems with the model crashing. Suggestions for improving the model for the domestic circuit are given in chapter 4.4. The improved equations that are to be used in the improved model are given by equation 4.18, 4.19 and 4.22.Also for the boiler circuit and the solar circuit there are improvements that can be done. The model presented here has a few shortcomings, but with some extra work, an improved model can be created. In the attachment (Bilaga 1) is a description of the used model and all the identified parameters.A qualitative assessment of the store was also performed based on the measurements and the modelling carried out. The following summary of this can be given: Hot Water PreparationThe principle for controlling the flow on the primary side seems to work well in order to achieve good stratification. Temperatures in the bottom of the store after a short use of hot water, at a coldwater temperature of 12°C, was around 28-30°C. This was almost independent of the temperature in the store and the DHW-flow.The measured UA-values of the heat exchangers are not very reliable, but indicates that the heat transfer rates are much better than for the Conus 500, and in the same range as for other stores tested at SERC.The function of the mixing valve is not perfect (see diagram 4.3, where Tout1 is the outlet hot water temperature, and Tdhwo and Tdhw1 is the inlet temperature to the hot and cold side of the valve respectively). The outlet temperature varies a lot with different temperatures in the storage and is going down from 61°C to 47°C before the cold port is fully closed. This gives a problem to find a suitable temperature setting and gives also a risk that the auxiliary heating is increased instead of the set temperature of the valve, when the hot water temperature is to low.Collector circuitThe UA-value of the collector heat exchanger is much higher than the value for Conus 500, and in the same range as the heat exchangers in other stores tested at SERC.Boiler circuitThe valve in the boiler circuit is used to supply water from the boiler at two different heights, depending on the temperature of the water. At temperatures from the boiler above 58.2°C, all the water is injected to the upper inlet. At temperatures below 53.9°C all the water is injected to the lower inlet. At 56°C the water flow is equally divided between the two inlets. Detailed studies of the behaviour at the upper inlet shows that better accuracy of the model would have been achieved using three double ports in the model instead of two. The shape of the upper inlet makes turbulence, that could be modelled using two different inlets. Heat lossesThe heat losses per m3 are much smaller for the Solus 1050, than for the Conus 500 Storage. However, they are higher than those for some good stores tested at SERC. The pipes that are penetrating the insulation give air leakage and cold bridges, which could be a major part of the losses from the storage. The identified losses from the bottom of the storage are exceptionally high, but have less importance for the heat losses, due to the lower temperatures in the bottom. High losses from the bottom can be caused by air leakage through the insulation at the pipe connections of the storage.

Rapport från/Report from IEA Task-14 mötet/meeting on solar heating system, Hameln Tyskland/Germany August 1992

Participation as observer at the meeting of Task 14 of IEA's Solar Heating and Cooling Projects held in Hameln, Germany has led to greater understanding of interesting developments underway in several countries. This will be of use during the development of small scale systems suitable for Swedish conditions. A summary of the work carried out by the working groups within Task 14 is given, with emphasis on the Domestic Hot Water group. Experiences of low-flow systems from several countries are related, and the conclusion is drawn that the maximum theoretical possible increase in performance of 20% has not been achieved due to poor heat exchangers and poor stratification in the storage tanks. Positive developments in connecting tubes and pumps is noted. Further participation as observer in Task 14 meetings is desired, and is looked on favourably by the members of the group. Another conclusion is that SERC should carry on with work on Swedish storage tanks, with emphasis on better stratification and heat exchangers, and possible modelling of system components. Finally a German Do-it-Vourself kit is described and judged in comparison with prefabricated models and Swedish Do-it-Yourself kits.