Numerical simulation of combined heat transfer phenomena (conduction, convection, and radiation). Application to he optimisation of thermal systems

Thermal systems interchanging heat and mass by conduction, convection, radiation (solar and thermal ) occur in many engineering applications like energy storage by solar collectors, window glazing in buildings, refrigeration of plastic moulds, air handling units etc. Often these thermal systems are composed of various elements for example a building with wall, windows, rooms, etc. It would be of particular interest to have a modular thermal system which is formed by connecting different modules for the elements, flexibility to use and change models for individual elements, add or remove elements without changing the entire code. A numerical approach to handle the heat transfer and fluid flow in such systems helps in saving the full scale experiment time, cost and also aids optimisation of parameters of the system. In subsequent sections are presented a short summary of the work done until now on the orientation of the thesis in the field of numerical methods for heat transfer and fluid flow applications, the work in process and the future work.

Experimental investigation of transient natural convection cooling of high prandtl number fluid with variable viscosity

Report for the scientific sojourn at the James Cook University, Australia, between June to December 2007. Free convection in enclosed spaces is found widely in natural and industrial systems. It is a topic of primary interest because in many systems it provides the largest resistance to the heat transfer in comparison with other heat transfer modes. In such systems the convection is driven by a density gradient within the fluid, which, usually, is produced by a temperature difference between the fluid and surrounding walls. In the oil industry, the oil, which has High Prandtl, usually is stored and transported in large tanks at temperatures high enough to keep its viscosity and, thus the pumping requirements, to a reasonable level. A temperature difference between the fluid and the walls of the container may give rise to the unsteady buoyancy force and hence the unsteady natural convection. In the initial period of cooling the natural convection regime dominates over the conduction contribution. As the oil cools down it typically becomes more viscous and this increase of viscosity inhibits the convection. At this point the oil viscosity becomes very large and unloading of the tank becomes very difficult. For this reason it is of primary interest to be able to predict the cooling rate of the oil. The general objective of this work is to develop and validate a simulation tool able to predict the cooling rates of high Prandtl fluid considering the variable viscosity effects.

Optimal exponent heat balance and refined integral methods applied to Stefan problems

When using a polynomial approximating function the most contentious aspect of the Heat Balance Integral Method is the choice of power of the highest order term. In this paper we employ a method recently developed for thermal problems, where the exponent is determined during the solution process, to analyse Stefan problems. This is achieved by minimising an error function. The solution requires no knowledge of an exact solution and generally produces significantly better results than all previous HBI models. The method is illustrated by first applying it to standard thermal problems. A Stefan problem with an analytical solution is then discussed and results compared to the approximate solution. An ablation problem is also analysed and results compared against a numerical solution. In both examples the agreement is excellent. A Stefan problem where the boundary temperature increases exponentially is analysed. This highlights the difficulties that can be encountered with a time dependent boundary condition. Finally, melting with a time-dependent flux is briefly analysed without applying analytical or numerical results to assess the accuracy.

Application of standard and refined heat balance integral methods to one-dimensional Stefan problems

The work in this paper concerns the study of conventional and refined heat balance integral methods for a number of phase change problems. These include standard test problems, both with one and two phase changes, which have exact solutions to enable us to test the accuracy of the approximate solutions. We also consider situations where no analytical solution is available and compare these to numerical solutions. It is popular to use a quadratic profile as an approximation of the temperature, but we show that a cubic profile, seldom considered in the literature, is far more accurate in most circumstances. In addition, the refined integral method can give greater improvement still and we develop a variation on this method which turns out to be optimal in some cases. We assess which integral method is better for various problems, showing that it is largely dependent on the specified boundary conditions.

Gaussian estimates for the density of the non-linear stochastic heat equation in any space dimension

In this paper, we establish lower and upper Gaussian bounds for the probability density of the mild solution to the stochastic heat equation with multiplicative noise and in any space dimension. The driving perturbation is a Gaussian noise which is white in time with some spatially homogeneous covariance. These estimates are obtained using tools of the Malliavin calculus. The most challenging part is the lower bound, which is obtained by adapting a general method developed by Kohatsu-Higa to the underlying spatially homogeneous Gaussian setting. Both lower and upper estimates have the same form: a Gaussian density with a variance which is equal to that of the mild solution of the corresponding linear equation with additive noise.

Improving the accuracy of heat balance integral methods applied to thermal problems with time dependent boundary conditions

In this paper the two main drawbacks of the heat balance integral methods are examined. Firstly we investigate the choice of approximating function. For a standard polynomial form it is shown that combining the Heat Balance and Refined Integral methods to determine the power of the highest order term will either lead to the same, or more often, greatly improved accuracy on standard methods. Secondly we examine thermal problems with a time-dependent boundary condition. In doing so we develop a logarithmic approximating function. This new function allows us to model moving peaks in the temperature profile, a feature that previous heat balance methods cannot capture. If the boundary temperature varies so that at some time t & 0 it equals the far-field temperature, then standard methods predict that the temperature is everywhere at this constant value. The new method predicts the correct behaviour. It is also shown that this function provides even more accurate results, when coupled with the new CIM, than the polynomial profile. Analysis primarily focuses on a specified constant boundary temperature and is then extended to constant flux, Newton cooling and time dependent boundary conditions.

Síntesi i propietats de vidres metàl•lics massissos i materials nanocomposites basats en Fe per aplicacions com a materials magnètics tous

Estudi realitzat a partir d’una estada a l’ Institut für Komplexe Materialien, Leibniz-Institut für Festkörper- und Werkstoffforschung Dresden, entre 2010 i 2011. S'ha explorat l'efecte de les condicions i influència dels elements d'aleació en la capacitat de formació de vidre, l'estructura i les propietats tèrmiques i magnètiques de vidres metàl•lics massissos i materials nanocristal•lins en base Fe. La producció d'aquests materials en forma de cintes de unes 20 micres de gruix ha estat àmpliament estudiada i s'ha vist que presenten unes propietats excel•lents com a materials magnètics tous. El propòsit general d'aquest projecte era l'obtenció de composicions òptimes amb alta capacitat de formar vidre i amb excel•lents propietats magnètiques com a materials magnètics tous combinat amb bones propietats mecàniques. El projecte prenia com a punt de partida l'aliatge [FeCoBSi]96Nb4 ja que és el que presenta millor capacitat de formar vidre i presenta una alta imantació de saturació i baix camp coercitiu. S'ha fet un estudi dels factors fonamentals que intervenen en la formació de l'estat vitri. La composició abans esmentada ha estat variada amb l'addició d'altres elements per estudiar com afecten aquests nous elements a les propietats, la formació de vidre i l'estructura dels aliatges resultants amb l'objectiu de millorar-ne les propietats magnètiques i la capacitat de formació de vidre. Entre altres s'ha usat el Zr, Mo, Y i el Gd per millorar la formació de vidre; i el Co i el Ni per millorar les propietats magnètiques a alta temperatura. S'han estudiat les relacions entre la capacitat de formació de vidre i la seva estabilitat tèrmica, la resistència a la cristal•lització i la estructura de l'aliatge resultant després del procés de solidificació. Per aquest estudi s'han determinat els mecanismes que controlen la transformació i la seva cinètica així com les fases que es formen durant el tractament tèrmic permetent la formulació de models predictius.

Resistance to reform: Reconsidering the role of individual-specific uncertainty

Individual-specific uncertainty may increase the chances of reform beingenacted and sustained. Reform may be more likely to be enacted because amajority of agents might end up losing little from reform and a minoritygaining a lot. Under certainty, reform would therefore be rejected, butit may be enacted with uncertainty because those who end up losing believethat they might be among the winners. Reform may be more likely to besustained because, in a realistic setting, reform will increase theincentives of agents to move into those economic activities that benefit.Agents who respond to these incentives will vote to sustain reform infuture elections, even if they would have rejected reform under certainty.These points are made using the trade-model of Fernandez and Rodrik (AER,1991).

A reusable smart interface for gas sensor resistance measurement

The advances of the semiconductor industry enable microelectromechanical systems sensors, signal conditioning logic and network access to be integrated into a smart sensor node. In this framework, a mixed-mode interface circuit for monolithically integrated gas sensor arrays was developed with high-level design techniques. This interface system includes analog electronics for inspection of up to four sensor arrays and digital logic for smart control and data communication. Although different design methodologies were used in the conception of the complete circuit, high-level synthesis tools and methodologies were crucial in speeding up the whole design cycle, enhancing reusability for future applications and producing a flexible and robust component.

Design and construction of precision heat fluxmeters

InAlAs/InGaAs/InP based high electron mobility transistor devices have been structurally and electrically characterized, using transmission electron microscopy and Raman spectroscopy and measuring Hall mobilities. The InGaAs lattice matched channels, with an In molar fraction of 53%, grown at temperatures lower than 530¿°C exhibit alloy decomposition driving an anisotropic InGaAs surface roughness oriented along [1math0]. Conversely, lattice mismatched channels with an In molar fraction of 75% do not present this lateral decomposition but a strain induced roughness, with higher strength as the channel growth temperature increases beyond 490¿°C. In both cases the presence of the roughness implies low and anisotropic Hall mobilities of the two dimensional electron gas.

Anomalous specific-heat jump in a two-component ultracold Fermi gas

The thermodynamic functions of a Fermi gas with spin population imbalance are studied in the temperature-asymmetry plane in the BCS limit. The low-temperature domain is characterized by an anomalous enhancement of the entropy and the specific heat above their values in the unpaired state, decrease of the gap and eventual unpairing phase transition as the temperature is lowered. The unpairing phase transition induces a second jump in the specific heat, which can be measured in calorimetric experiments. While the superfluid is unstable against a supercurrent carrying state, it may sustain a metastable state if cooled adiabatically down from the stable high-temperature domain. In the latter domain the temperature dependence of the gap and related functions is analogous to the predictions of the BCS theory.

Heat fluctuations in Brownian transducers

The heat fluctuation probability distribution function in Brownian transducers operating between two heat reservoirs is studied. We find, both analytically and numerically, that the recently proposed fluctuation theorem for heat exchange [C. Jarzynski and D. K. Wojcik, Phys. Rev. Lett. 92, 230602 (2004)] has to be applied carefully when the coupling mechanism between both baths is considered. We also conjecture how to extend such a relation when an external work is present.