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.

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.

The peroxiredoxin Tpx1 is essential as a H202-scavenger during aerobic growth in fission yeast

Peroxiredoxins are known to interact with hydrogen peroxide (H2O2) and to participate in oxidant scavenging, redox signal transduction, and heat-shock responses. The two-cysteine peroxiredoxin Tpx1 of Schizosaccharomyces pombe has been characterized as the H2O2 sensor that transduces the redox signal to the transcription factor Pap1. Here, we show that Tpx1 is essential for aerobic, but not anaerobic, growth. We demonstrate that Tpx1 has an exquisite sensitivity for its substrate, which explains its participation in maintaining low steady-state levels of H2O2. We also show in vitro and in vivo that inactivation of Tpx1 by oxidation of its catalytic cysteine to a sulfinic acid is always preceded by a sulfinic acid form in a covalently linked dimer, which may be important for understanding the kinetics of Tpx1 inactivation. Furthermore, we provide evidence that a strain expressing Tpx1.C169S, lacking the resolving cysteine, can sustain aerobic growth, and we show that small reductants can modulate the activity of the mutant protein in vitro, probably by supplying a thiol group to substitute for cysteine 169.

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.

Magnetization versus heat treatment in rapidly solidified NdFeB alloys

NdFeB melt-spun amorphous or partially amorphous alloys of four compositions were prepared. Their crystallization kinetics induced by thermal treatment was studied by differential scanning calorimetry and scanning and transmission electron microscopy. Scanning electron microscopy demonstrated that heterogeneous nucleation occurs preferentially at the ribbon surface which was in contact with the wheel. The explicit form of the kinetic equation that best describes the first stage of crystallization under high undercooling conditions was obtained for each alloy. From the crystallization results, the lower part of the experimental time-temperature-transformation curves was deduced for each alloy and extrapolated up to the high-temperature limit of their validity. Microstructural observations showed a typical size of the microcrystals obtained by heat treatment of ~100 nm. From the magnetic properties measured with a vibrating sample magnetometer, the same magnetic behavior of partially crystallized alloys is observed regardless of the temperature of annealing provided the same crystallization fraction, x, is achieved, at least for small values of x (typically ~10%).

Design and construction of precision heat fluxmeters

The possibility of local elastic instabilities is considered in a first¿order structural phase transition, typically a thermoelastic martensitic transformation, with associated interfacial and volumic strain energy. They appear, for instance, as the result of shape change accommodation by simultaneous growth of different crystallographic variants. The treatment is phenomenological and deals with growth in both thermoelastic equilibrium and in nonequilibrium conditions produced by the elastic instability. Scaling of the transformed fraction curves against temperature is predicted only in the case of purely thermoelastic growth. The role of the transformation latent heat on the relaxation kinetics is also considered, and it is shown that it tends to increase the characteristic relaxation times as adiabatic conditions are approached, by keeping the system closer to a constant temperature. The analysis also reveals that the energy dissipated in the relaxation process has a double origin: release of elastic energy Wi and entropy production Si. The latter is shown to depend on both temperature rate and thermal conduction in the system.

A classical particle in heat bath under the influence of external noise

A very simple model of a classical particle in a heat bath under the influence of external noise is studied. By means of a suitable hypothesis, the heat bath is reduced to an internal colored noise (OrnsteinUhlenbeck noise). In a second step, an external noise is coupled to the bath. The steady state probability distributions are obtained.