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.

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.

Multiphase modeling of melting / solidification with high density variations using XFEM

La modélisation de la cryolite, utilisée dans la fabrication de l’aluminium, implique plusieurs défis, notament la présence de discontinuités dans la solution et l’inclusion de la difference de densité entre les phases solide et liquide. Pour surmonter ces défis, plusieurs éléments novateurs ont été développés dans cette thèse. En premier lieu, le problème du changement de phase, communément appelé problème de Stefan, a été résolu en deux dimensions en utilisant la méthode des éléments finis étendue. Une formulation utilisant un multiplicateur de Lagrange stable spécialement développée et une interpolation enrichie a été utilisée pour imposer la température de fusion à l’interface. La vitesse de l’interface est déterminée par le saut dans le flux de chaleur à travers l’interface et a été calculée en utilisant la solution du multiplicateur de Lagrange. En second lieu, les effets convectifs ont été inclus par la résolution des équations de Stokes dans la phase liquide en utilisant la méthode des éléments finis étendue aussi. Troisièmement, le changement de densité entre les phases solide et liquide, généralement négligé dans la littérature, a été pris en compte par l’ajout d’une condition aux limites de vitesse non nulle à l’interface solide-liquide pour respecter la conservation de la masse dans le système. Des problèmes analytiques et numériques ont été résolus pour valider les divers composants du modèle et le système d’équations couplés. Les solutions aux problèmes numériques ont été comparées aux solutions obtenues avec l’algorithme de déplacement de maillage de Comsol. Ces comparaisons démontrent que le modèle par éléments finis étendue reproduit correctement le problème de changement phase avec densités variables.

Thermodynamic modeling and microcalorimetry of nanostructured materials for capacitive thermal management of electronics

Transient power dissipation profiles in handheld electronic devices alternate between high and low power states depending on usage. Capacitive thermal management based on phase change materials potentially offers a fan-less thermal management for such transient profiles. However, such capacitive management becomes feasible only if there is a significant enhancement in the enthalpy change per unit volume of the phase change material since existing bulk materials such as paraffin fall short of requirements. In this thesis I propose novel nanostructured thin-film materials that can potentially exhibit significantly enhanced volumetric enthalpy change. Using fundamental thermodynamics of phase transition, calculations regarding the enhancement resulting from superheating in such thin film systems is conducted. Furthermore design of a microfabricated calorimeter to measure such enhancements is explained in detail. This work advances the state-of-art of phase change materials for capacitive cooling of handheld devices.

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.

MANIPULATION OF IONS, ELECTRONS, AND PHOTONS IN 2D MATERIALS BY ION INTERCALATION

2D materials have attracted tremendous attention due to their unique physical and chemical properties since the discovery of graphene. Despite these intrinsic properties, various modification methods have been applied to 2D materials that yield even more exciting results. Among all modification methods, the intercalation of 2D materials provides the highest possible doping and/or phase change to the pristine 2D materials. This doping effect highly modifies 2D materials, with extraordinary electrical transport as well as optical, thermal, magnetic, and catalytic properties, which are advantageous for optoelectronics, superconductors, thermoelectronics, catalysis and energy storage applications. To study the property changes of 2D materials, we designed and built a planar nanobattery that allows electrochemical ion intercalation in 2D materials. More importantly, this planar nanobattery enables characterization of electrical, optical and structural properties of 2D materials in situ and real time upon ion intercalation. With this device, we successfully intercalated Li-ions into few layer graphene (FLG) and ultrathin graphite, heavily dopes the graphene to 0.6 x 10^15 /cm2, which simultaneously increased its conductivity and transmittance in the visible range. The intercalated LiC6 single crystallite achieved extraordinary optoelectronic properties, in which an eight-layered Li intercalated FLG achieved transmittance of 91.7% (at 550 nm) and sheet resistance of 3 ohm/sq. We extend the research to obtain scalable, printable graphene based transparent conductors with ion intercalation. Surfactant free, printed reduced graphene oxide transparent conductor thin film with Na-ion intercalation is obtained with transmittance of 79% and sheet resistance of 300 ohm/sq (at 550 nm). The figure of merit is calculated as the best pure rGO based transparent conductors. We further improved the tunability of the reduced graphene oxide film by using two layers of CNT films to sandwich it. The tunable range of rGO film is demonstrated from 0.9 um to 10 um in wavelength. Other ions such as K-ion is also studied of its intercalation chemistry and optical properties in graphitic materials. We also used the in situ characterization tools to understand the fundamental properties and improve the performance of battery electrode materials. We investigated the Na-ion interaction with rGO by in situ Transmission electron microscopy (TEM). For the first time, we observed reversible Na metal cluster (with diameter larger than 10 nm) deposition on rGO surface, which we evidenced with atom-resolved HRTEM image of Na metal and electron diffraction pattern. This discovery leads to a porous reduced graphene oxide sodium ion battery anode with record high reversible specific capacity around 450 mAh/g at 25mA/g, a high rate performance of 200 mAh/g at 250 mA/g, and stable cycling performance up to 750 cycles. In addition, direct observation of irreversible formation of Na2O on rGO unveils the origin of commonly observed low 1st Columbic Efficiency of rGO containing electrodes. Another example for in situ characterization for battery electrode is using the planar nanobattery for 2D MoS2 crystallite. Planar nanobattery allows the intrinsic electrical conductivity measurement with single crystalline 2D battery electrode upon ion intercalation and deintercalation process, which is lacking in conventional battery characterization techniques. We discovered that with a “rapid-charging” process at the first cycle, the lithiated MoS2 undergoes a drastic resistance decrease, which in a regular lithiation process, the resistance always increases after lithiation at its final stage. This discovery leads to a 2- fold increase in specific capacity with with rapid first lithiated MoS2 composite electrode material, compare with the regular first lithiated MoS2 composite electrode material, at current density of 250 mA/g.