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.

Spatial scales of desert locust gregarization

Central to swarm formation in migratory locusts is a crowding-induced change from a “solitarious” to a “gregarious” phenotype. This change can occur within the lifetime of a single locust and accrues across generations. It represents an extreme example of phenotypic plasticity. We present computer simulations and a laboratory experiment that show how differences in resource distributions, conspicuous only at small spatial scales, can have significant effects on phase change at the population level; local spatial concentration of resource induces gregarization. Simulations also show that populations inhabiting a locally concentrated resource tend to change phase rapidly and synchronously in response to altered population densities. Our results show why information about the structure of resource at small spatial scales should become key components in monitoring and control strategies.

Applications of Thermal Mechanical Compression Tests in food powder analysis

A new Thermal Mechanical Compression Test (TMCT) was applied for glass-rubber transition and melting analyses of food powders and crystals. The TMCT technique measures the phase change of a material based on mechanical changes during the transition. Whey, honey, and apple juice powders were analyzed for their glass-rubber transition temperatures. Sucrose and glucose monohydrate crystals were analyzed for their melting temperatures. The results were compared to the values obtained by conventional DSC and TMA techniques. The new TMCT technique provided the results that were very close to the conventional techniques. This technique can be an alternative to analyze glass-rubber transition of food, pharmaceutical, and chemical dry products.

Stirring in 3-d spherical models of convection in the Earth's mantle

On a global scale basalts from mid-ocean ridges are strikingly more homogeneous than basalts from intraplate volcanism. The observed geochemical heterogeneity argues strongly for the existence of distinct reservoirs in the Earth's mantle. It is an unresolved problem of Geodynamics as to how these findings can be reconciled with large-scale convection. We review observational constraints, and investigate stirring properties of numerical models of mantle convection. Conditions in the early Earth may have supported layered convection with rapid stirring in the upper layers. Material that has been altered near the surface is transported downwards by small-scale convection. Thereby a layer of homogeneous depleted material develops above pristine mantle. As the mantle cools over Earth history, the effects leading to layering become reduced and models show the large-scale convection favoured for the Earth today. Laterally averaged, the upper mantle below the lithosphere is least affected by material that has experienced near-surface differentiation. The geochemical signature obtained during the previous episode of small-scale convection may be preserved there for the longest time. Additionally, stirring is less effective in the high viscosity layer of the central lower mantle [1, 2], supporting the survival of medium-scale heterogeneities there. These models are the first, using 3-d spherical geometry and mostly Earth-like parameters, to address the suggested change of convective style. Although the models are still far from reproducing our planet, we find that proposal might be helpful towards reconciling geochemical and geophysical constraints.

A study of liquid-liquid extraction in a sieve plate column

The literature relating to sieve plate liquid extraction columns and relevant hydrodynamic phenomena have been surveyed. Mass transfer characteristics during drop formation, rise and coalescence, and related models were also reviewed. Important design parameters i.e. flooding, dispersed phase hold-up, drop size distribution, mean drop size, coalescence/flocculation zone height beneath a plate and jetting phenomena were investigated under non-mass transfer and mass transfer conditions in a 0.45m diameter, 2.3m high sieve plate column. This column had provision for four different plate designs, and variable plate spacing and downcomer heights, and the system used was Clairsol `350' (dispersed) - acetone - deionised water (continuous) with either direction of mass transfer. Drop size distributions were best described by the functions proposed by Gal-or, and then Mugele-Evans. Using data from this study and the literature, correlations were developed for dispersed phase hold-up, mean drop size in the preferred jetting regime and in the non-jetting regime, and coalescence zone height. A method to calculate the theoretical overall mass transfer coefficient allowing for the range of drop sizes encountered in the column gave the best fit to experimental data. This applied the drop size distribution diagram to estimate the volume percentage of stagnant, circulating and oscillating drops in the drop population. The overall coefficient Kcal was then calculated as the fractional sum of the predicted individual single drop coefficients and their proportion in the drop population. In a comparison between the experimental and calculated overall mass transfer coefficients for cases in which all the drops were in the oscillating regime (i.e. 6.35mm hole size plate), and for transfer from the dispersed(d) to continuous(c) phase, the film coefficient kd predicted from the Rose-Kintner correlation together with kc from that of Garner-Tayeban gave the best representation. Droplets from the 3.175mm hole size plate, were of a size to be mainly circulating and oscillating; a combination of kd from the Kronig-Brink (circulating) and Rose-Kintner (oscillating) correlations with the respective kc gave the best agreement. The optimum operating conditions for the SPC were identified and a procedure proposed for design from basic single drop data.

Melting and solidification of Zinc-Aluminium alloys

Following a scene-setting introduction are detailed reviews of the relevant scientific principles, thermal analysis as a research tool and the development of the zinc-aluminium family of alloys. A recently introduced simultaneous thermal analyser, the STA 1500, its use for differential thermal analysis (DTA) being central to the investigation, is described, together with the sources of support information, chemical analysis, scanning electron microscopy, ingot cooling curves and fluidity spiral castings. The compositions of alloys tested were from the binary zinc-aluminium system, the ternary zinc-aluminium-silicon system at 30%, 50% and 70% aluminium levels, binary and ternary alloys with additions of copper and magnesium to simulate commercial alloys and five widely used commercial alloys. Each alloy was shotted to provide the smaller, 100mg, representative sample required for DTA. The STA 1500 was characterised and calibrated with commercially pure zinc, and an experimental procedure established for the determination of DTA heating curves at 10°C per minute and cooling curves at 2°C per minute. Phase change temperatures were taken from DTA traces, most importantly, liquidus from a cooling curve and solidus from both heating and cooling curves. The accepted zinc-aluminium binary phase diagram was endorsed with the added detail that the eutectic is at 5.2% aluminium rather than 5.0%. The ternary eutectic trough was found to run through the points, 70% Al, 7.1% Si, 545°C; 50% Al, 3.9% Si, 520°C; 30% Al, 1.4% Si, 482°C. The dendrite arm spacing in samples after DTA increased with increasing aluminium content from 130m at 30% to 220m at 70%. The smallest dendrite arm spacing of 60m was in the 30% aluminium 2% silicon alloy. A 1kg ingot of the 10% aluminium binary alloy, insulated with Kaowool, solidified at the same 2°C per minute rate as the DTA samples. A similar sized sand casting was solidified at 3°C per minute and a chill casting at 27°C per minute. During metallographic examination the following features were observed: heavily cored phase which decomposed into ' and '' on cooling; needles of the intermetallic phase FeAl4; copper containing ternary eutectic and copper rich T phase.

Solar stills:a comprehensive review of designs, performance and material advances

The demand for fresh water production is growing day by day with the increase in world population and with industrial growth. Use of desalination technology is increasing to meet this demand. Among desalination technologies, solar stills require low maintenance and are readily affordable; however their productivity is limited. This paper aims to give a detailed review about the various types of solar stills, covering passive and active designs, single- and multi-effect types, and the various modifications for improved productivity including reflectors, heat storage, fins, collectors, condensers, and mechanisms for enhancing heat and mass transfer. Photovoltaic-thermal and greenhouse type solar stills are also covered. Material advances in the area of phase change materials and nanocomposites are very promising to enhance further performance; future research should be carried out in these and other areas for the greater uptake of solar still technology.

An experimental study of the effects of double-diffusive convention processes during the solification of binary alloys

Experiments were conducted to show the effects of thermal and geometric boundary conditions on the liquid pool of a binary alloy system which is undergoing phase change, solidification. Transparent analogue solutions were selected for study and experimental apparatus were designed and built. Thermal distribution and concentration data were collected and analysed for the melt pool of various selected geometries and boundary conditions of the systems under study. The data indicate-that characteristic flows develop for both Hypereutectic and Hypoeutectic concentration levels and that the development of macrosegregation and microsegregation defects in continuous casting materials can be minimised by the adjustment of the process variables.

Analisi termiche sull'impiego di materiali a cambiamento di fase (PCM) nelle pavimentazioni stradali

Il presente lavoro di tesi mira a studiare l’utilizzo di aggregati artificiali (PLA) costituiti da aggregati leggeri (LWA) impregnati di materiali a cambiamento di fase (Phase-Change Materials, PCM) nei conglomerati bituminosi. L’obiettivo della tesi è quello di dimostrare che l’utilizzo di questi materiali nelle sovrastrutture stradali, grazie alla proprietà di cambiare fase (da solida a liquida e viceversa) in funzione della temperatura, induce una liberazione di calore. La conseguenza immediata dell’utilizzo di questi materiali è la ridotta necessità di manutenzione invernale, abbattendo i costi di ripristino della pavimentazione. Inoltre l’utilizzo di PLA non deve pregiudicare l’aspetto prestazionale e la vita utile dell’infrastruttura.

Design of Latent Thermal Energy Storage Systems

This dissertation documents the results of a theoretical and numerical study of time dependent storage of energy by melting a phase change material. The heating is provided along invading lines, which change from single-line invasion to tree-shaped invasion. Chapter 2 identifies the special design feature of distributing energy storage in time-dependent fashion on a territory, when the energy flows by fluid flow from a concentrated source to points (users) distributed equidistantly on the area. The challenge in this chapter is to determine the architecture of distributed energy storage. The chief conclusion is that the finite amount of storage material should be distributed proportionally with the distribution of the flow rate of heating agent arriving on the area. The total time needed by the source stream to ‘invade’ the area is cumulative (the sum of the storage times required at each storage site), and depends on the energy distribution paths and the sequence in which the users are served by the source stream. Chapter 3 shows theoretically that the melting process consists of two phases: “invasion” thermal diffusion along the invading line, which is followed by “consolidation” as heat diffuses perpendicularly to the invading line. This chapter also reports the duration of both phases and the evolution of the melt layer around the invading line during the two-dimensional and three-dimensional invasion. It also shows that the amount of melted material increases in time according to a curve shaped as an S. These theoretical predictions are validated by means of numerical simulations in chapter 4. This chapter also shows that the heat transfer rate density increases (i.e., the S curve becomes steeper) as the complexity and number of degrees of freedom of the structure are increased, in accord with the constructal law. The optimal geometric features of the tree structure are detailed in this chapter. Chapter 5 documents a numerical study of time-dependent melting where the heat transfer is convection dominated, unlike in chapter 3 and 4 where the melting is ruled by pure conduction. In accord with constructal design, the search is for effective heat-flow architectures. The volume-constrained improvement of the designs for heat flow begins with assuming the simplest structure, where a single line serves as heat source. Next, the heat source is endowed with freedom to change its shape as it grows. The objective of the numerical simulations is to discover the geometric features that lead to the fastest melting process. The results show that the heat transfer rate density increases as the complexity and number of degrees of freedom of the structure are increased. Furthermore, the angles between heat invasion lines have a minor effect on the global performance compared to other degrees of freedom: number of branching levels, stem length, and branch lengths. The effect of natural convection in the melt zone is documented.

Heating the Early Universe : Numerical Methods and Their Analysis

During the epoch when the first collapsed structures formed (6<z<50) our Universe went through an extended period of changes. Some of the radiation from the first stars and accreting black holes in those structures escaped and changed the state of the Intergalactic Medium (IGM). The era of this global phase change in which the state of the IGM was transformed from cold and neutral to warm and ionized, is called the Epoch of Reionization.In this thesis we focus on numerical methods to calculate the effects of this escaping radiation. We start by considering the performance of the cosmological radiative transfer code C2-Ray. We find that although this code efficiently and accurately solves for the changes in the ionized fractions, it can yield inaccurate results for the temperature changes. We introduce two new elements to improve the code. The first element, an adaptive time step algorithm, quickly determines an optimal time step by only considering the computational cells relevant for this determination. The second element, asynchronous evolution, allows different cells to evolve with different time steps. An important constituent of methods to calculate the effects of ionizing radiation is the transport of photons through the computational domain or ``ray-tracing''. We devise a novel ray tracing method called PYRAMID which uses a new geometry - the pyramidal geometry. This geometry shares properties with both the standard Cartesian and spherical geometries. This makes it on the one hand easy to use in conjunction with a Cartesian grid and on the other hand ideally suited to trace radiation from a radially emitting source. A time-dependent photoionization calculation not only requires tracing the path of photons but also solving the coupled set of photoionization and thermal equations. Several different solvers for these equations are in use in cosmological radiative transfer codes. We conduct a detailed and quantitative comparison of four different standard solvers in which we evaluate how their accuracy depends on the choice of the time step. This comparison shows that their performance can be characterized by two simple parameters and that the C2-Ray generally performs best.

Three-dimensional free surface modelling in an unstructured mesh environment for metal processing applications

In the casting of metals, tundish flow, welding, converters, and other metal processing applications, the behaviour of the fluid surface is important. In aluminium alloys, for example, oxides formed on the surface may be drawn into the body of the melt where they act as faults in the solidified product affecting cast quality. For this reason, accurate description of wave behaviour, air entrapment, and other effects need to be modelled, in the presence of heat transfer and possibly phase change. The authors have developed a single-phase algorithm for modelling this problem. The Scalar Equation Algorithm (SEA) (see Refs. 1 and 2), enables the transport of the property discontinuity representing the free surface through a fixed grid. An extension of this method to unstructured mesh codes is presented here, together with validation. The new method employs a TVD flux limiter in conjunction with a ray-tracing algorithm, to ensure a sharp bound interface. Applications of the method are in the filling and emptying of mould cavities, with heat transfer and phase change.

Multi-physics modeling of aluminium reduction cells

Aluminium cells involve a range of complex physical processes which act simultaneously to provide a narrow satisfactory operating range. These processes involve electromagnetic fields, coupled with heat transfer and phase change, two phase fluid flow with a range of complexities plus the development of stress in the cell structure. All of these phenomena are coupled in some significant sense and so to provide a comprehensive model of these processes involves their representation simultaneously. Conventionally, aspects of the process have been modeled separately using uncoupled estimates of the effects of the other phenomena; this has enabled the use of standard commercial CFD and FEA tools. In this paper we will describe an approach to the modeling of aluminium cells which describes all the physics simultaneously. This approach uses a finite volume approximation for each of the phenomena and facilitates their interactions directly in the modeling-the complex geometries involved are addressed by using unstructured meshes. The very challenging issues to be overcome in this venture will be outlined and some preliminary results will be shown.

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.

A finite-volume shallow layer method for the MHD instabilities in an aluminium production cell

An electrolytic cell for Aluminium production contains molten metal and molten electrolyte, which are subject to high dc-currents and magnetic fields. Lorentz forces arising from the cross product of current and magnetic field may amplify natural gravity waves at the interface between the two fluids, leading to short circuits in extreme cases. The external magnetic field and current distribution in the production cell is computed through a detailed finite element analysis at Torino Polytechnic. The results are then used to compute the magnetohydrodynamic and thermal effects in the aluminium/electrolyte bath. Each cell has lateral dimensions of 6m x 2m, whilst the bath depth is only 30cm. the electrically resistive electrolyte path, which is critical in the operation of the cell, has layer depth of only a few centimetres below each carbon anode. Because the shallow dimensions of the liquid layer a finite-volume shallow-layer technique has been used at Greenwich to compute the resulting flow-field and interface perturbations. The information obtained from this method, i.e. depth averaged velocities and aluminium/electrolyte interface position is then embedded in the three-dimensional finite volume code PHYSICA and will be used to compute the heat transfer and phase change in the cell.