Thermal energy storage in Swedish single family houses : a case study

In a Nordic climate, space heating (SH) and domestic hot water (DHW) used in buildings constitute a considerable part of the total energy use in the country. For 2010, energy used for SH and DHW amounted to almost 90 TWh in Sweden which corresponds to 60 % of the energy used in the residential and service sector, or almost 24 % of the total final energy use for the country. Storing heat and cold with the use of thermal energy storage (TES) can be one way of increasing the energy efficiency of a building by opening up possibilities for alternative sources of heat or cold through a reduced mismatch between supply and demand. Thermal energy storage without the use of specific control systems are said to be passive and different applications using passive TES have been shown to increase energy efficiency and/or reduce power peaks of systems supplying the heating and cooling needs of buildings, as well as having an effect on the indoor climate. Results are however not consistent between studies and focus tend to be on the reduction of cooling energy or cooling power peaks. In this paper, passive TES introduced through an increased thermal mass in the building envelope to two single family houses with different insulation standard is investigated with building energy simulations. A Nordic climate is used and the focus of this study is both on the reduction of space heating demand and space heating power, as well as on reduction of excess temperatures in residential single family houses without active cooling systems. Care is taken to keep the building envelope characteristics other than the thermal mass equal for all cases so that any observations made can be derived to the change in thermal mass. Results show that increasing the sensible thermal mass in a single family house can reduce the heating demand only slightly (1-4 %) and reduce excess temperatures (temperatures above 24 degrees C) by up to 20 %. Adding a layer of PCM (phase change materials) to the light building construction can give similar reduction in heating demand and excess temperatures, however the phase change temperature is important for the results.

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)

PCM thermal insulation in buildings

This paper presents the results of a numerical and experimental study of phase change material (PCM) filled walls and roofs under real operational conditions to achieve passive thermal comfort. The numerical part of the study was based on a one-dimensional model for the phase change problem controlled by pure conduction. Real radiation data was used to determine the external face temperature. The numerical treatment was based upon using finite difference approximations and the ADI scheme. The results obtained were compared with field measurements. The experimental set-up consisted of a small room with movable roof and side wall. The roof was constructed in the traditional way but with the phase change material enclosed. Thermocouples were distributed across the cross section of the roof. Another roof, identical but without the PCM, was also used during comparative tests. The movable wall was also constructed as is done traditionally but with the PCM enclosed. Again, thermocouples were distributed across the wall thickness to enable measurement of the local temperatures. Another wall, identical but without the PCM, was also used during comparative tests. The PCM used in the numerical and experimental tests was composed of a mixture of two commercial grades of glycol in order to obtain the required fusion temperature range. Comparison between the simulation results and the experiments indicated good agreement. Field tests also indicated that the PCM used was adequate and that the concept was effective in maintaining the indoor temperature very close to the established comfort limits. Further economical analysis indicated that the concept could effectively help in reducing the electric energy consumption and improving the energy demand pattern. © 1997 by John Wiley & Sons, Ltd.

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.

Night time performance of a storage integrated solar thermophotovoltaic (SISTPV) system

Energy storage at low maintenance cost is one of the key challenges for generating electricity from the solar energy. This paper presents the theoretical analysis (verified by CFD) of the night time performance of a recently proposed conceptual system that integrates thermal storage (via phase change materials) and thermophotovoltaics for power generation. These storage integrated solar thermophotovoltaic (SISTPV) systems are attractive owing to their simple design (no moving parts) and modularity compared to conventional Concentrated Solar Power (CSP) technologies. Importantly, the ability of high temperature operation of these systems allows the use of silicon (melting point of 1680 K) as the phase change material (PCM). Silicon's very high latent heat of fusion of 1800 kJ/kg and low cost ($1.70/kg), makes it an ideal heat storage medium enabling for an extremely high storage energy density and low weight modular systems. In this paper, the night time operation of the SISTPV system optimised for steady state is analysed. The results indicate that for any given PCM length, a combination of small taper ratio and large inlet hole-to-absorber area ratio are essential to increase the operation time and the average power produced during the night time. Additionally, the overall results show that there is a trade-off between running time and the average power produced during the night time. Average night time power densities as high as 30 W/cm(2) are possible if the system is designed with a small PCM length (10 cm) to operate just a few hours after sun-set, but running times longer than 72 h (3 days) are possible for larger lengths (50 cm) at the expense of a lower average power density of about 14 W/cm(2). In both cases the steady state system efficiency has been predicted to be about 30%. This makes SISTPV systems to be a versatile solution that can be adapted for operation in a broad range of locations with different climate conditions, even being used off-grid and in space applications.

Ultra high temperature latent heat energy storage and thermophotovoltaic energy conversion

A conceptual energy storage system design that utilizes ultra high temperature phase change materials is presented. In this system, the energy is stored in the form of latent heat and converted to electricity upon demand by TPV (thermophotovoltaic) cells. Silicon is considered in this study as PCM (phase change material) due to its extremely high latent heat (1800 J/g or 500 Wh/kg), melting point (1410 C), thermal conductivity (~25 W/mK), low cost (less than $2/kg or $4/kWh) and abundance on earth. The proposed system enables an enormous thermal energy storage density of ~1 MWh/m3, which is 10e20 times higher than that of lead-acid batteries, 2e6 times than that of Li-ion batteries and 5e10 times than that of the current state of the art LHTES systems utilized in CSP (concentrated solar power) applications. The discharge efficiency of the system is ultimately determined by the TPV converter, which theoretically can exceed 50%. However, realistic discharge efficiencies utilizing single junction TPV cells are in the range of 20e45%, depending on the semiconductor bandgap and quality, and the photon recycling efficiency. This concept has the potential to achieve output electric energy densities in the range of 200-450 kWhe/m3, which is comparable to the best performing state of the art Lithium-ion batteries.

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.

Numerical investigation on the heat transfer enhancement using micro/nano phase-change particulate flow

The introduction of phase change material fluid and nanofluid in micro-channel heat sink design can significantly increase the cooling capacity of the heat sink because of the unique features of these two kinds of fluids. To better assist the design of a high performance micro-channel heat sink using phase change fluid and nanofluid, the heat transfer enhancement mechanism behind the flow with such fluids must be completely understood. ^ A detailed parametric study is conducted to further investigate the heat transfer enhancement of the phase change material particle suspension flow, by using the two-phase non-thermal-equilibrium model developed by Hao and Tao (2004). The parametric study is conducted under normal conditions with Reynolds numbers of Re = 90–600 and phase change material particle concentrations of ϵp ≤ 0.25, as well as extreme conditions of very low Reynolds numbers (Re < 50) and high phase change material particle concentration (ϵp = 50%–70%) slurry flow. By using the two newly-defined parameters, named effectiveness factor ϵeff and performance index PI, respectively, it is found that there exists an optimal relation between the channel design parameters L and D, particle volume fraction ϵp, Reynolds number Re, and the wall heat flux qw. The influence of the particle volume fraction ϵp, particle size dp, and the particle viscosity μ p, to the phase change material suspension flow, are investigated and discussed. The model was validated by available experimental data. The conclusions will assist designers in making their decisions that relate to the design or selection of a micro-pump suitable for micro or mini scale heat transfer devices. ^ To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and its many numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed, through incorporating the different forces acting on the nanoparticles to the two-component lattice Boltzmann model. It is found that the nanofluid has better heat transfer enhancement at low Reynolds numbers, and the Brownian motion effect of the nanoparticles will be weakened by the increase of flow speed. ^

Numerical Investigation on the Heat Transfer Enhancement Using Micro/Nano Phase-Change Particulate Flow

The introduction of phase change material fluid and nanofluid in micro-channel heat sink design can significantly increase the cooling capacity of the heat sink because of the unique features of these two kinds of fluids. To better assist the design of a high performance micro-channel heat sink using phase change fluid and nanofluid, the heat transfer enhancement mechanism behind the flow with such fluids must be completely understood. A detailed parametric study is conducted to further investigate the heat transfer enhancement of the phase change material particle suspension flow, by using the two-phase non-thermal-equilibrium model developed by Hao and Tao (2004). The parametric study is conducted under normal conditions with Reynolds numbers of Re=600-900 and phase change material particle concentrations ¡Ü0.25 , as well as extreme conditions of very low Reynolds numbers (Re < 50) and high phase change material particle concentration (0.5-0.7) slurry flow. By using the two newly-defined parameters, named effectiveness factor and performance index, respectively, it is found that there exists an optimal relation between the channel design parameters, particle volume fraction, Reynolds number, and the wall heat flux. The influence of the particle volume fraction, particle size, and the particle viscosity, to the phase change material suspension flow, are investigated and discussed. The model was validated by available experimental data. The conclusions will assist designers in making their decisions that relate to the design or selection of a micro-pump suitable for micro or mini scale heat transfer devices. To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and its many numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed, through incorporating the different forces acting on the nanoparticles to the two-component lattice Boltzmann model. It is found that the nanofluid has better heat transfer enhancement at low Reynolds numbers, and the Brownian motion effect of the nanoparticles will be weakened by the increase of flow speed.

Simulação de sistemas solares passivos usando materiais de mudança de fase (PCMs)

O gradual consumo de energia primária a nível mundial deu origem a uma crise não só ambiental como também económica, proveniente das limitações das reservas energéticas e do fornecimento. Estas inquietações têm levado a um estudo cada vez mais aprofundado no que concerne à eficiência energética de edifícios. É neste contexto que surge o estudo da aplicação dos materiais de mudança de fase (PCM) na térmica dos edifícios. O presente trabalho consiste no estudo da influência da introdução de materiais de mudança de fase no comportamento térmico de um edifício. Foi analisada uma simulação numérica para um sistema solar passivo de ganho direto, por um programa de simulação designado EXTEND™. Efetuou-se, também, o estudo paramétrico de determinadas propriedades associadas aos materiais de mudança de fase. Após a análise, foi testado o impacto, em termos de conforto, da utilização de materiais de mudança de fase em sistemas solares passivos de ganho direto.

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.

BATTERY POWERED PORTABLE VAPOR COMPRESSION CYCLE SYSTEM WITH PCM CONDENSER

A battery powered air-conditioning device was developed to provide an improved thermal comfort level for individuals in inadequately cooled environments. This device is a battery powered air-conditioning system with the phase change material (PCM) for heat storage. The condenser heat is stored in the PCM during the cooling operation and is discharged while the battery is charged by using the vapor compression cycle as a thermosiphon loop. The main focus of the current research was on the development of the cooling system. The cooling capacity of the vapor compression cycle measured was 165.6 W with system COP at 2.85. It was able to provide 2 hours cooling without discharging heat to the ambient. The PCM was recharged in nearly 8 hours under thermosiphon mode.

PCM-based in-memory computing: architectures, circuits and applications

Analog In-memory Computing (AIMC) has been proposed in the context of Beyond Von Neumann architectures as a valid strategy to reduce internal data transfers energy consumption and latency, and to improve compute efficiency. The aim of AIMC is to perform computations within the memory unit, typically leveraging the physical features of memory devices. Among resistive Non-volatile Memories (NVMs), Phase-change Memory (PCM) has become a promising technology due to its intrinsic capability to store multilevel data. Hence, PCM technology is currently investigated to enhance the possibilities and the applications of AIMC. This thesis aims at exploring the potential of new PCM-based architectures as in-memory computational accelerators. In a first step, a preliminar experimental characterization of PCM devices has been carried out in an AIMC perspective. PCM cells non-idealities, such as time-drift, noise, and non-linearity have been studied to develop a dedicated multilevel programming algorithm. Measurement-based simulations have been then employed to evaluate the feasibility of PCM-based operations in the fields of Deep Neural Networks (DNNs) and Structural Health Monitoring (SHM). Moreover, a first testchip has been designed and tested to evaluate the hardware implementation of Multiply-and-Accumulate (MAC) operations employing PCM cells. This prototype experimentally demonstrates the possibility to reach a 95% MAC accuracy with a circuit-level compensation of cells time drift and non-linearity. Finally, empirical circuit behavior models have been included in simulations to assess the use of this technology in specific DNN applications, and to enhance the potentiality of this innovative computation approach.

NZEB e os novos revestimentos

Dissertação de natureza científica para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização de Edificações