Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase change”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1982

Comparison of advanced thermal and electrical storage for parabolic dish solar thermal power systems

Parabolic dish solar concentrator cluster concepts are explored, with attention given to thermal storage systems coupled to Stirling and Brayton cycle power conversion devices. Sensible heat storage involving molten salt (NaOH), liquid sodium, and solid cordierite bricks are considered for 1500 F thermal storage systems. Latent heat storage with NaF-MgF2 phase change materials are explored in terms of passive, active, and direct contact designs. Comparisons are made of the effectiveness of thermal storage relative to redox, Na-S, Zn-Cl, and Zn-Br battery storage systems. Molten lead trickling down through a phase change eutectic, the NaF-MgF2, formed the direct contact system. Heat transport in all systems is effected through Inconel pipes. Using a cost goal of 120-150 mills/kWh as the controlling parameter, sensible heat systems with molten salts transport with either Stirling or Brayton engines, or latent heat systems with Stirling engines, and latent heat-Brayton engine with direct contact were favored in the analyses. Battery storage systems, however, offered the most flexibility of applications.

Fujita, T.↗
NASA NTRS2005

PCM Passive Cooling System Containing Active Subsystems

A multistage system has been proposed for cooling a circulating fluid that is subject to intermittent intense heating. The system would be both flexible and redundant in that it could operate in a basic passive mode, either sequentially or simultaneously with operation of a first, active cooling subsystem, and either sequentially or simultaneously with a second cooling subsystem that could be active, passive, or a combination of both. This flexibility and redundancy, in combination with the passive nature of at least one of the modes of operation, would make the system more reliable, relative to a conventional cooling system. The system would include a tube-in-shell heat exchanger, within which the space between the tubes would be filled with a phase-change material (PCM). The circulating hot fluid would flow along the tubes in the heat exchanger. In the basic passive mode of operation, heat would be conducted from the hot fluid into the PCM, wherein the heat would be stored temporarily by virtue of the phase change.

Blanding, David E.↗
DOE OSTI2024

Retrofittable Thermal Switches for Dynamic Building Envelopes Integrated with Thermal Energy Storage: Preprint

Buildings in the United States consume about 40 quadrillion BTU of primary energy annually, which accounts for the nation's 40% of total energy use, 75% of all electricity use, and 35% of the net carbon emissions. Deploying thermal energy storage in the form of phase change material (PCM) in building envelopes is an effective method to reduce space heating/cooling loads, provide load shedding, and shift demand to periods of lower energy cost. However, the full potential of PCM-integrated envelopes can only be realized if the PCM undergoes complete phase change using free ambient heating/cooling, and the stored energy is effectively transferred between the exterior and the interior environments. Conventional thermal insulation (with a fixed thermal resistance) limits PCM utilization, particularly with the increasing emphasis on higher R-value in building envelopes, which negatively affects the energy-saving potential of a PCM-integrated envelope. In contrast, dynamic building envelopes integrated with PCMs provide the option of varying the thermal resistance based on the indoor and outdoor conditions, thereby enhancing utilization of free ambient cooling and heating to charge/discharge the PCM thermal storage, reducing the buildings' heating and cooling load, and shifting the peak energy demand. In this study, we demonstrate innovative retrofittable thermal switches in the form of the insertable plugs inside an insulation to provide variable thermal resistance depending on the operating temperature and direction of temperature gradient, thus allowing preferential directional heat flow. Notably, they are passive in nature, requiring no external power, and work solely based on the ambient temperature.

buildings↗
NASA NTRS2011

Metrology Arrangement for Measuring the Positions of Mirrors of a Submillimeter Telescope

The position of the secondary mirror of a submillimeter telescope with respect to the primary mirror needs to be known .0.03 mm in three dimensions. At the time of this reporting, no convenient, reasonably priced arrangement that offers this capability exists. The solution proposed here relies on measurement devices developed and deployed for the GeoSAR mission, and later adapted for the ISAT (Innovative Space Based Radar Antenna Technology) demonstration. The measurement arrangement consists of four metrology heads, located on an optical bench, attached to the secondary mirror. Each metrology head has a dedicated target located at the edge of the primary mirror. One laser beam, launched from the head and returned by the target, is used to measure distance. Another beam, launched from a beacon on the target, is monitored by the metrology head and generates a measurement of the target position in the plane perpendicular to the laser beam. A 100-MHz modulation is carried by a collimated laser beam. The relevant wavelength is the RF one, 3 m, divided by two, because the light carries it to the target and back. The phase change due to travel to the target and back is measured by timing the zero-crossing of the RF modulation, using a 100-MHz clock. In order to obtain good resolution, the 100-MHz modulation signal is down-converted to 1 kHz. Then, the phase change corresponding to the round-trip to the target is carried by a 1-kHz signal. Since the 100-MHz clock beats 100,000 times during one period of the 1-kHz signal, the least-significant-bit (LSB) resolution is LSB = 0.015 mm.

Abramovici, Alex↗
NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗
NASA NTRS1981

Automatic thermal switch

An automatic thermal switch to control heat flow includes a first thermally conductive plate, a second thermally conductive plate and a thermal transfer plate pivotally mounted between the first and second plates. A phase change power unit, including a plunger connected to the transfer plate, is in thermal contact with the first thermally conductive plate. A biasing element, connected to the transfer plate, biases the transfer plate in a predetermined position with respect to the first and second plates. When the phase change power unit is actuated by an increase in heat transmitted through the first plate, the plunger extends and pivots the transfer plate to vary the thermal conduction between the first and second plates through the transfer plate. The biasing element, transfer plate and piston can be arranged to provide either a normally closed or normally open thermally conductive path between the first and second plates.

Wing, L. D.↗
NASA NTRS2001

The Extraterrestrial Materials Simulation Laboratory

In contrast to fly-by and orbital missions, in situ missions face an incredible array of challenges in near-target navigation, landing site selection, descent, landing, science operations, sample collection and handling, drilling, anchoring, subsurface descent, communications, and contamination. The wide range of materials characteristics and environments threaten mission safety and success. For example, many physical properties are poorly characterized, including strength, composition, heterogeneity, phase change, texture, thermal properties, terrain features, atmospheric interaction, and stratigraphy. Examples of the range of materials properties include, for example: (1) Comets, with a possible compressive strength ranging from a light fluff to harder than concrete: 10(exp 2) to 10 (exp 8) Pa; (2) Europa, including a possible phase change at the surface, unknown strength and terrain roughness; and (3) Titan, with a completely unknown surface and possible liquid ocean. Additional information is contained in the original extended abstract.

Green, J. R.↗
DOE OSTI2022

Effect of Thermal Expansion Coefficient, Viscosity and Melting Range in Simulation of PCM Embedded Heat Exchangers With and Without Fins

Phase change material heat exchangers (PCMHX) have an important role in integrating renewable energy systems. PCMHX can offer high storage density at various temperatures of interest and can be used for grid load shifting purposes. Numerical models enable engineers to estimate PCMHX performance for different design parameters and operating conditions. Modeling phase change phenomena is challenging due to the complex time-dependent nature of the process. The accuracy of models is highly sensitive to PCM thermo-physical properties. Thermal expansion coefficient (β), viscosity (μ) and melting temperature range (MR) of a PCM are important properties, especially when natural convection is not negligible. In PCMHX modeling, using less than accurate values for these properties can have significant impact on the simulation outcomes. These properties and discussions thereof are not readily available in the literature. This paper presents a brief review of the literature and a numerical study investigating the model sensitivity to the above-mentioned properties for PCMHX with and without fins. CFD is used to evaluate the charging (melting) phenomena. The study quantifies the impact of uncertainty in these properties on the melting rate and temperature profile. Constant wall temperature was considered as heat source with no heat loss to ambient. Here, the results show that β and μ has significant effect on the melting rate and evolution of the melting front. For a non-finned domain when comparing results for different published values of β and μ, the deviation in melting time can be up to 12.9% and 57.6% respectively. For high wall temperatures, change in melting range did not impact melting time. But when the wall temperature is reduced, up to 9.8% deviation in melting time is observed.

CFD↗
NASA NTRS1973

The end of the iron-core age.

The terrestrial planets aggregated essentially from small particles, to begin as solid cool bodies with the same general compositions, and there is no possibility of an iron-core developing within any of them at any stage. Their differing internal and surface properties receive ready explanation from their different masses which determine whether the pressures within are sufficient to bring about phase-changes. The claim that the terrestrial core can be identified by means of shock-wave data as nickel-iron is based on theoretical misconception, whereas the actual seismic data establish an uncompressed-density value much lower than any such mixture could have. The onset of the Ramsey phase-change in the earth takes the form of a rapid initial collapse to produce a large core in metallic state which thereafter continues to grow secularly as a result of radioactive heating and leads to reduction of surface-area at long last adequate to account for folded and thrusted mountain-building.

Lyttleton, R. A.↗
NASA NTRS1972

Ascent heat transfer rate distribution on the North American Rockwell delta wing orbiter and the General Dynamics/Convair booster at a Mach number of 8 (mated)

A wind tunnel test program to determine aerodynamic interference heating on the North American Rockwell orbiter mated with the General Dynamics Convair booster is discussed. The tests were conducted at the Arnold Engineering Development Center (AEDC) in Tunnel B of the von Karman Gas Dynamics Facility (VKF). The test period was June 1971. Heat-transfer rates were determined by the phase-change paint technique on 0.013-scale Stycast models using Tempilaq as the surface temperature indicator. The nominal test conditions were: Mach 8, free-stream unit length Reynolds numbers of 1.25 x one million and 2.55 x one million angles of attack of -5, 0, +5 deg. Model details, test conditions, phase-change paint photographs and reduced heat-transfer coefficients are presented.

Matthews, R. K.↗
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
NASA NTRS1972

Heat transfer tests of the McDonnell-Douglas delta wing orbiter mated with -17A booster at Mach number 8

A wind tunnel test program to determine aerodynamic heat transfer distributions on the McDonnell-Douglas configurations is reported. The tests were conducted at the Arnold Engineering Development Center (AEDC) in Tunnel B of the von Karman Gas Dynamics Facility (VKF). Heat-transfer rates were determined by the phase-change paint technique on 0.011-scale Stycast models using Tempilaq as the surface temperature indicator. The nominal test conditions were; Mach 8, freestream unit Reynolds numbers of 0.8 x one million, 2.5 x one million, and 3.7 x one million, and angles of attack of -5 deg, 0 deg, +5deg. Model details, test conditions, phase-change paint photographs and reduced heat-transfer coefficients are presented.

Matthews, R. K.↗
NASA NTRS1972

Heat transfer investigation of Langley Research Center transition models at a Mach number of 8, volume 2

The results are presented of a wind tunnel test program to determine aerodynamic heat transfer distributions on delta body and straight body transition models of the space shuttle. Heat transfer rates were determined by the phase-change paint technique on Stycast and RTV models using Tempilag as the surface temperature indicator. The nominal test conditions were: Mach 8, length Reynolds numbers of 5 million and 7.4 million, and angles of attack of 20, 40, and 60 deg. Model details, test conditions, and reduced heat transfer data are included. Data reduction of the phase-change paint photographs was performed by utilizing a new technique.

Matthews, R. K.↗
NASA NTRS1980

Heat storage capability of a rolling cylinder using Glauber's salt

The rolling cylinder phase change heat storage concept was developed to the point where a prototype design is completed and a cost analysis is prepared. A series of experimental and analytical tasks are defined to establish the thermal, mechanical, and materials behavior of rolling cylinder devices. These tasks include: analyses of internal and external heat transfer; performance and lifetime testing of the phase change materials; corrosion evaluation; development of a mathematical model; and design of a prototype and associated test equipment.

Herrick, C. S.↗
NASA NTRS1976

Relation of a contracting earth to the apparent accelerations of the sun and moon

It has been concluded that the theoretical ratio of the apparent secular accelerations of the moon and the sun implied by tidal theory is far higher than the ratio of observed values. This paper shows that the discrepancy can be entirely removed within the range of uncertainty of the various quantities involved if terms making allowance for a secularly decreasing terrestrial moment of inertia are included in the dynamical equations for the apparent secular accelerations. Such a decrease has already been predicted on the basis of the phase-change hypothesis for the nature of earth's core. The development of this hypothesis is reviewed, the modified dynamical equations are derived, and values of tidal couples implied by existing observational data are examined. The numerically determined rate of decrease of earth's moment of inertia is found to be in close agreement with that predicted by the phase-change hypothesis. The results indicate that a diminishing gravitational constant is unlikely, that the moon would have been in close proximity to earth only 1 billion years ago, and that earth's rotation period at that time would have been about 5 hr.

Lyttleton, R. A.↗
DOE OSTI2021

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Feasibility Analysis for the Use of Retrofitted Air-Conditioners Using Thermal Energy Storage (TES) for High Ambient Temperature (HAT) Countries

In high ambient temperature (HAT) countries, summer temperatures exceed 35℃, degrading the performance of air-conditioning systems and straining power grids. In this paper, a feasibility analysis was conducted to investigate potential savings during the peak using latent Thermal Energy Storage (TES) at near room phase-change temperatures, replacing condensers, thus, minimizing temperature lifts. Weather data, building loads, and baseline air-conditioning systems data were gathered for Dubai. A transient vapor-compression model in Modelica was used to compare the COP, total power input and cooling capacity of the air-conditioner at peak hours when ambient temperatures range between 35 – 45 ⁰C versus the TES-Phase-Change Material (PCM) melting temperatures from 22 – 28⁰C. The results indicate that the TES-PCM can enhance the system COP at the peak by a factor 1.4 and 2 for during for outdoor temperatures of 35 – 40⁰C, and 40 – 45⁰C, respectively. Lower melting temperature PCMs were able to reduce the required power input by 30-50%, with more savings occurring at higher temperature days. On the other hand, higher temperature PCMs enhancements were minimal especially at outdoor ambient temperatures ranging between 35 – 40⁰C. Improvements to the cooling capacity range from 8 – 18 % for the outdoor temperature range of 35 – 45 ⁰C. An economic analysis was conducted to find the potential saving in utility costs for 30%, 60%, and 90% of the space cooling demands of Dubai, and find the trade-off points between utility savings and cost of PCM-TES implementation. If the peak loads are to be shifted by 6 hours daily, the percentage utility savings for the city is 18%. Using estimated costs of the PCM-TES, ranging from $200-500/kWh, the daily load shifting hours were estimated to range from 4 hours at the lowest cost systems to 2.5 hours at the highest costs.

25 ENERGY STORAGE↗