Improvement of the thermal performance of commercial freezers using phase change materials (PCM)

The aim of this work is to improve the thermal performance of commercial freezers using phase change materials (PCM) and the prove the importance of the correct PCM selection (melting temperature) for each application. To do this a vertical freezer with and without macro encapsulation PCM are used to assess the benefits of using PCM in minimizing temperature fluctuations and to maintain at the lowest possible temperature the freezer under heat of losses.

Thermal energy storage with phase change materials in building envelopes

Els materials de canvi de fase (PCM) han estat considerats per a l’emmagatzematge tèrmic en edificis des de 1980. Amb la inclusió dels PCM en plaques de guix, guix, formigó o altres materials que s’utilitzen per a cobrir les parets, l’emmagatzematge tèrmic pot ser part de les estructures fins i tot en edificis lleugers. Les noves tècniques de microencapsulació han obert moltes possibilitats en aplicacions per a edificis. El treball que es presenta és el desenvolupament d’un formigó innovador mesclat amb PCM microencapsulat, amb un punt de fusió de 26 oC i una entalpia de canvi de fase de 110 kJ/kg. El primer experiment va ser la inclusió del PCM microencapsulat dins del formigó i la construcció d’una caseta amb aquest nou formigó-PCM. Es va construir una segona caseta al costat de la primera amb les mateixes característiques i orientació però amb formigó convencional que serveix com a referència. Durant els anys 2005 i 2006 es va analitzar el comportament d’ambdues casetes i més tard es va edificar un mur Trombe a la paret sud de totes dues per investigar la seva influència durant la tardor i l’hivern.

A moving-mesh finite element method and its application to the numerical solution of phase-change problems

A distributed Lagrangian moving-mesh finite element method is applied to problems involving changes of phase. The algorithm uses a distributed conservation principle to determine nodal mesh velocities, which are then used to move the nodes. The nodal values are obtained from an ALE (Arbitrary Lagrangian-Eulerian) equation, which represents a generalization of the original algorithm presented in Applied Numerical Mathematics, 54:450--469 (2005). Having described the details of the generalized algorithm it is validated on two test cases from the original paper and is then applied to one-phase and, for the first time, two-phase Stefan problems in one and two space dimensions, paying particular attention to the implementation of the interface boundary conditions. Results are presented to demonstrate the accuracy and the effectiveness of the method, including comparisons against analytical solutions where available.

Phase change and pyroelectricity in natural hemimorphite

The phase change of a natural hemimorphite sample from Minas Gerais (Brazil) was investigated by two X-ray diffraction (XRD) methods and by near-infrared reflectance spectroscopy. Applying successive thermal treatments, the crystal structure undergoes two orientation conversions. The first one occurs at about 550 degrees C, and it was revealed by the Laue method. Below 500 degrees C, the water molecules were partially expelled without changing the crystal structure. A fact that supports this statement is the sequential disappearance of the water bands at 1400 and 1900 nm by thermal treatment. The second conversion takes place below 939 degrees C. Moreover, at 972 degrees C a phase change to the willemite mineral (alpha-Zn(2)SiO(4)) has been observed. This last conversion was confirmed by the power XRD. In addition, natural hemimorphite displayed a high pyroelectricity, which is related both to the absence of inverse centre and to the presence of molecular water and hydroxyl groups in the crystal structure.

A high dynamic range method for the direct readout of a dynamic phase change in homodyne interferometers

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Analysis of charge-transport properties in GST materials for next generation phase-change memory devices

The quest for universal memory is driving the rapid development of memories with superior all-round capabilities in non-volatility, high speed, high endurance and low power. The memory subsystem accounts for a significant cost and power budget of a computer system. Current DRAM-based main memory systems are starting to hit the power and cost limit. To resolve this issue the industry is improving existing technologies such as Flash and exploring new ones. Among those new technologies is the Phase Change Memory (PCM), which overcomes some of the shortcomings of the Flash such as durability and scalability. This alternative non-volatile memory technology, which uses resistance contrast in phase-change materials, offers more density relative to DRAM, and can help to increase main memory capacity of future systems while remaining within the cost and power constraints. Chalcogenide materials can suitably be exploited for manufacturing phase-change memory devices. Charge transport in amorphous chalcogenide-GST used for memory devices is modeled using two contributions: hopping of trapped electrons and motion of band electrons in extended states. Crystalline GST exhibits an almost Ohmic I(V) curve. In contrast amorphous GST shows a high resistance at low biases while, above a threshold voltage, a transition takes place from a highly resistive to a conductive state, characterized by a negative differential-resistance behavior. A clear and complete understanding of the threshold behavior of the amorphous phase is fundamental for exploiting such materials in the fabrication of innovative nonvolatile memories. The type of feedback that produces the snapback phenomenon is described as a filamentation in energy that is controlled by electron–electron interactions between trapped electrons and band electrons. The model thus derived is implemented within a state-of-the-art simulator. An analytical version of the model is also derived and is useful for discussing the snapback behavior and the scaling properties of the device.

Two-phase change in CO2, Antarctic temperature and global climate during Termination II

The end of the Last Glacial Maximum (Termination I), roughly 20 thousand years ago (ka), was marked by cooling in the Northern Hemisphere, a weakening of the Asian monsoon, a rise in atmospheric CO2 concentrations and warming over Antarctica. The sequence of events associated with the previous glacial–interglacial transition (Termination II), roughly 136 ka, is less well constrained. Here we present high-resolution records of atmospheric CO2 concentrations and isotopic composition of N2—an atmospheric temperature proxy—from air bubbles in the EPICA Dome C ice core that span Termination II. We find that atmospheric CO2 concentrations and Antarctic temperature started increasing in phase around 136 ka, but in a second phase of Termination II, from 130.5 to 129 ka, the rise in atmospheric CO2 concentrations lagged that of Antarctic temperature unequivocally. We suggest that during this second phase, the intensification of the low-latitude hydrological cycle resulted in the development of a CO2 sink, which counteracted the CO2 outgassing from the Southern Hemisphere oceans over this period.

Local Heat Application for the Treatment of Buruli Ulcer: Results of a Phase II Open Label Single Center Non Comparative Clinical Trial.

BACKGROUND Buruli ulcer (BU) is a necrotizing skin disease most prevalent among West African children. The causative organism, Mycobacterium ulcerans, is sensitive to temperatures above 37°C. We investigated the safety and efficacy of a local heat application device based on phase change material. METHODS In a phase II open label single center noncomparative clinical trial (ISRCTN 72102977) under GCP standards in Cameroon, laboratory confirmed BU patients received up to 8 weeks of heat treatment. We assessed efficacy based on the endpoints 'absence of clinical BU specific features' or 'wound closure' within 6 months ("primary cure"), and 'absence of clinical recurrence within 24 month' ("definite cure"). RESULTS Of 53 patients 51 (96%) had ulcerative disease. 62% were classified as World Health Organization category II, 19% each as category I and III. The average lesion size was 45 cm(2). Within 6 months after completion of heat treatment 92.4% (49 of 53, 95% confidence interval [CI], 81.8% to 98.0%) achieved cure of their primary lesion. At 24 months follow-up 83.7% (41 of 49, 95% CI, 70.3% to 92.7%) of patients with primary cure remained free of recurrence. Heat treatment was well tolerated; adverse effects were occasional mild local skin reactions. CONCLUSIONS Local thermotherapy is a highly effective, simple, cheap and safe treatment for M. ulcerans disease. It has in particular potential as home-based remedy for BU suspicious lesions at community level where laboratory confirmation is not available. CLINICAL TRIALS REGISTRATION ISRCT 72102977.

Thin film bearings with phase change of lubricant

The possible evaporation of lubricant in fluid film bearings has been investigated theoretically and by experiment using a radial flow hydrostatic bearing supplied with liquid refrigerant R114. Good correlation between measured and theoretical values was obtained using a bespoke computational fluid dynamic model in which the flow was assumed to be laminar and adiabatic. The effects of viscous dissipation and vapour generation within the fluid film are fully accounted for by applying a fourth order Runge-Kutta routine to satisfy the radial and filmwise transverse constraints of momentum, energy and mass conservation. The results indicate that the radial velocity profile remains parabolic while the flow remains in the liquid phase and that the radial rate of enthalpy generation is then constant across the film at a given radius. The results also show that evaporation will commence at a radial location determined by geometry and flow conditions and in fluid layers adjacent to the solid boundaries. Evaporation is shown to progress in the radial direction and the load carrying capacity of such a bearing is reduced significantly. Expressions for the viscosity of the liquid/vapour mixture found in the literature survey have not been tested against experimental data. A new formulation is proposed in which the suitable choice of a characteristic constant yields close representation to any of these expressions. Operating constraints imposed by the design of the experimental apparatus limited the extent of the surface over which evaporation could be obtained, and prevented clear identification of the most suitable relationship for the viscosity of the liquid/vapour mixture. The theoretical model was extended to examine the development of two phase flow in a rotating shaft face seal of uniform thickness. Previous theoretical analyses have been based on the assumption that the radial velocity profile of the flow is always parabolic, and that the tangential component of velocity varies linearly from the value at the rotating surface, to zero at the stationary surface. The computational fluid dynamic analysis shows that viscous shear and dissipation in the fluid adjacent to the rotating surface leads to developing evaporation with a consequent reduction in tangential shear forces. The tangential velocity profile is predicted to decay rapidly through the film, exhibiting a profile entirely different to that assumed by previous investigators. Progressive evaporation takes place close to the moving wall and does not occur completely at a single radial location, as has been claimed in earlier work.

Analyse numérique du changement de phase solide-liquide autour d'un tube horizontal soumis à une sollicitation thermique périodique

Ce mémoire présente un modèle mathématique et numérique pour analyser le comportement d’une unité de stockage thermique à changement de phase solide-liquide représentée par un tube autour duquel se trouve le matériau à changement de phase. Le système est soumis à une charge oscillant entre le chauffage et le refroidissement. Une analyse d’ordre de grandeur permet de prédire le comportement du système en fonction des principaux nombres adimensionnels. Un paramètre adimensionnel est proposé pour délimiter les concepts dans lesquels la conduction domine par rapport à ceux où la convection naturelle domine. L’étude dévoile l’impact des paramètres de conception de l’unité de stockage thermique sur son fonctionnement et approfondit les connaissances dans le domaine du changement de phase avec convection naturelle. Différents indicateurs ont été développés pour analyser la performance du système, tels que les dimensions de la zone affectée thermiquement, le volume fondu ou solidifié et une analyse fréquentielle. Des corrélations sont proposées pour déterminer facilement le comportement du système.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Hygrothermal behaviour of earthen plasters for sustainable housing construction

A ready-mixed and several laboratory formulated mortars were produced and tested in fresh state and after hardening, simulating a masonry plaster for indoor application. All the mortars used a clayish earth from the same region and different compositions of aggregates, eventually including fibres and a phase change material. All the formulated mortars were composed by 1:3 volumetric proportions of earth and aggregate. Tests were developed for consistency, fresh bulk density, thermal conductivity, capillary absorption and drying, water vapour permeability and sorption-desorption. The use of PCM changed drastically the workability of the mortars and increased their capillary absorption. The use of fibres and variations on particle size distribution of the mixtures of sand that were used had no significant influence on tested properties. But particularly the good workability of these mortars and the high capacity of sorption and desorption was highlighted. With this capacity plasters made with these mortars are able to adsorb water vapour from indoor atmosphere when high levels of relative humidity exist and release water vapour when the indoor atmosphere became too dry. This fact makes them able to contribute passively for a healthier indoor environment. The technical, ecological and environmental advantages of the application of plasters with this type of mortars are emphasized, with the aim of contributing for an increased use for new or existent housing.

Mortars with incorporation of PCM based in different binders: mechanical and thermal behavior

Currently we are witnessing a huge concern of society with the parameters of comfort of the buildings and the energetic consumptions. It is known that there is a huge consumption of non-renewable sources of energy. Thus, it is urgent to develop and explore ways to take advantage of renewable sources of energy by improving the energy efficiency of buildings. The mortars with incorporation of phase change materials (PCM) have the ability to regulate the temperature inside buildings, contributing to the thermal comfort and reduction of the use of heating and cooling equipment, using only the energy supplied by the sun. However, the incorporation of phase change materials in mortars modifies its characteristics. The main purpose of this study was mechanical and thermal characterization of mortars with incorporation of PCM in mortars based in different binders. The binders studied were aerial lime, hydraulic lime, gypsum and cement. For each type of binder a reference composition (0% PCM) and a composition with incorporation of 40% of PCM were developed. It was possible to observe that the incorporation of PCM in mortars caused differences in properties such as workability, compressive strength, flexural strength and adhesion, however leads to an improvement of thermal behavior.

Estudo de revestimentos ativos para o isolamento térmico de edifícios

Dissertação de mestrado em Sustentabilidade do Ambiente Construído