Influència del context geològic i hidrogeològic en l'eficiència energètica dels sistemes geotèrmics de molt baixa entalpia

L’elevat consum energètic de les societats actuals, així com la impossibilitat de sostenir-lo a llarg termini implica la cerca de noves fonts d’energia. D’aquesta manera, en el camp de la climatització residencial, l’energia geotèrmica de molt baixa entalpia es posiciona com una alternativa als recursos energètics actuals. Així, els primers metres de subsòl presenten una temperatura adequada per al seu aprofitament calorífic, mitjançant els sistemes geotèrmics de bomba de calor. Si bé a nivell energètic, són sistemes, intrínsecament, molt eficients, el seu rendiment pot patir importants variacions davant dels canvis en les condicions del medi geològic i hidrogeològic. Especialment, els col·lectors de calor verticals, treballen, amb freqüència, en el si de les formacions hidrogeològiques. En aquest sentit, els canvis del nivell hidràulic i de la temperatura de l’aigua de l’aqüífer es manifesten amb variacions de la conductivitat tèrmica equivalent i del flux subterrani d’aigua, que alhora, aquestes, es tradueixen en alteracions del flux subterrani de calor. Davant d’aquest fet, l’avaluació quantitativa de l’efecte d’aquestes fluctuacions que es presenta en aquest treball mostra petites variacions del nivell hidràulic i de la temperatura de l’aigua comporten canvis molt notables en l’eficiència dels sistemes verticals de bomba de calor geotèrmica.

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.

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.

Brownian pump powered by a white-noise flasing ratchet

A Brownian pump of particles powered by a stochastic flashing ratchet mechanism is studied. The pumping device is embedded in a finite region and bounded by particle reservoirs. In the steady state, we exactly calculate the spatial density profile, the concentration ratio between both reservoirs and the particle flux. We propose a simulation framework for the consistent evaluation of such observable quantities.

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%).