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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Elasticity of two-dimensional ferroelectrics across their paraelectric phase transformation

The mechanical behavior of two-dimensional (2D) materials across 2D phase changes is unknown, and the finite temperature (T) elasticity of paradigmatic SnSe monolayers—ferroelectric 2D materials turning paraelectric as their unit cell turns from a rectangle into a square—is described here in a progressive manner. To begin with, their zero–T total energy landscape gives way to (Boltzmann-like) averages from which the elastic behavior is determined. Furthermore, these estimates are complemented with results from the strain-fluctuation method, which employs the energy landscape or ab initio molecular dynamics data. All approaches capture the coalescence of elastic moduli < C 11 (T) > = < C 22 (T) > due to the structural transformation. The broad evolution and sudden changes of elastic parameters < C 11 (T) >, < C 22 (T) >, and < C 12 (T) > of these atomically thin phase-change membranes establishes a heretofore overlooked connection among 2D materials and soft matter.

2-dimensional systems↗

DOE OSTI2022

Development of an Integrated Thermal Energy Storage Heat Exchanger for Concentrating Solar Power

In this design and development effort a latent-heat energy storage heat exchanger was designed for a 10MW e sCO 2 power plant as a basis for detailed design and testing of a new phase-change heat exchanger technology. High efficiency thermal conductors of graphite foam or folded copper were employed within chloride salt phase-change volumes bounded by alternating high density compact plate-fin heat exchange cells in this new heat exchanger technology. Test articles were constructed at Brayton Energy and tested at Argonne National Laboratory. Copper-conducting samples were constructed and tested at Brayton Energy.

14 SOLAR ENERGY↗

DOE OSTI2024

Time resolved X-ray diffraction using the FIDDLE diagnostic at NIF: Preliminary assessment of diffraction precision

The Flexible Imaging Diffraction Diagnostic for Laser Experiments (FIDDLE) is a new diagnostic at the National Ignition Facility (NIF) designed to observe in situ solid-solid phase changes at high pressures using time resolved X-ray diffraction. FIDDLE currently incorporates five Icarus Ultrafast X-ray Imager sensors that take 2 ns snapshots and can be tuned to collect X-rays for tens of ns. The platform utilizes the laser power at NIF for both the laser drive and the generation of 10 keV X-rays for ~10 ns using a Ge backlighter foil. We aim to use FIDDLE to observe diffraction at different times during compression to probe the kinetics of phase changes. Pb undergoes two solid-solid phase transitions during ramp compression: from FCC to HCP and HCP to BCC. Results will be reported on some of the first shots using the FIDDLE diagnostic at NIF on ramp compressed Pb to a peak pressure of ~110 GPa and a single undriven CeO2 calibration shot. A discussion of the uncertainties in the observed diffraction is included.

Vennari, C. E.↗

DOE OSTI2020

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

The role of moisture in MgCl 2 salt: A multiscale approach to TES performance

This study presents a comprehensive multiscale analysis to evaluate the influence of moisture on the thermal performance of thermal energy storage systems using magnesium chloride (MgCl 2 ) as the phase-changing material. The system uses graphite foam with 90% relative density to enhance thermal conductivity. The analysis includes thermal conductivity calculations and specific heat capacity for different hydrate phases of MgCl 2 : Anhydrous, Mono, Di, Tetra, and Hexa. These properties were evaluated for the first time using the phonon density of states from Density Functional Theory simulations. Results showed that thermal conductivity decreased, while specific heat capacity increased by a factor of two as the phase changed from anhydrous to hexahydrate. Meso-scale models were created to account for the anisotropy of graphite foam and property variations of the MgCl 2 hydrate phases. Asymptotic Expansion homogenization simulations determined the anisotropic thermal conductivity for all phases. This unique methodology improved simulation accuracy, which matched experimental data for anhydrous MgCl 2 . A parametric study examined various operating conditions and their effect on TES performance. It revealed that higher charging temperatures did not enhance exergy efficiency, but increased discharge mass flow rates improved it due to better heat transfer. Thermal performance evaluated by the exergy efficiency remained consistent across hydrate systems under the tested conditions. Furthermore, the study suggests that this uniformity is linked to missing key data, particularly latent heat and melting point for different hydrates. Overall, it highlights the importance of multiscale effects and accurate material properties in designing and optimizing TES systems.

Molten salt degradation↗

DOE OSTI2025

Phenomena Identification and Ranking Table (PIRT) for heat pipes

Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in heat pipes, including capillary, phase change, turbulence, and compressibility effects, there are high uncertainties in the predictability of their operational regimes and performance. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. Additional analyses and discussion are provided for phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics. The discussions included the recognizing challenges and proposing future research directions for both modeling and simulation and experimental efforts. Additionally, the report addresses phenomena with medium importance and low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geyser boiling, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. In conclusion, this comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2025

Phenomena Identification and Ranking Table (PIRT) for Heat Pipes

This Phenomena Identification and Ranking Table (PIRT) report provides an evaluation of key phenomena affecting the performance and operational regimes of heat pipes, particularly in the context of heat pipe microreactors (HPMRs). Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in the heat pipe, including phase change, turbulent transition, and compressibility effects, among others, there is high uncertainty in identifying and ranking the important phenomena affecting the operation of heat pipes and the current knowledge for their modeling and simulation and experimental measurements and instrumentation. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. The report analyzes phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics, discussing challenges and future research directions for improving their modeling and simulation and experimental measurements. Additionally, the report addresses phenomena with low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geysering, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. This comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2026

Design and Characterization of the 162.5 MHz RF Component Layout for Fermilab's PIP-II Reference Line

The Proton Improvement Plan II (PIP-II) Reference Line at Fermi National Accelerator Laboratory distributes phase-stable radio-frequency (RF) signals throughout the accelerator. PIP-II requires stable timing and phase reference signals so its accelerating cavities transfer energy to the particle beam at the correct point in each RF cycle. The full system includes 162.5, 325, and 650 MHz sections corresponding to the frequency sections of the PIP-II Linac, with this project focusing on the 162.5 MHz section. The Reference Line must provide a phase stable source signal while responding to phase changes caused by environmental conditions or system drift. Its RF components will be mounted on aluminum heat plates inside a temperature controlled enclosure to further limit temperature-driven phase changes. To prepare the system for manufacture, the project reviewed component functions and dimensions, developed a computer-aided design (CAD) model, and arranged the hardware to support short cable paths, grounding, fastener access, and maintenance. Several full scale, three dimensional printed prototypes allowed the available components to be mounted and inspected. These fit checks revealed mechanical conflicts and guided revisions to component placement, countersink geometry, labeling, and plate thickness. Electrical characterization was also performed on selected RF hardware to compare its measured behavior with the performance metrics that were set for our design. Overall, the project produced a manufacturable 162.5 MHz layout, physical fit-check prototypes, and documented electrical measurements that support review before metal fabrication. The 325 and 650 MHz layouts remain future work because they require additional minor mechanical changes.

Subedi, Harsheet [Unlisted, US, CA] (ORCID:0009000↗

DOE OSTI2023

Multivariate Flexible Framework with High Usable Hydrogen Capacity in a Reduced Pressure Swing Process

Step-shaped adsorption-desorption of gaseous payloads by flexible metal-organic frameworks can facilitate the delivery of large usable capacities with significantly reduced energetic penalties. This is desirable for the storage, transport, and delivery of H 2 , as prototypical adsorbents require large swings in pressure and temperature to achieve usable capacities approaching their total capacities. However, the weak physisorption of H 2 typically necessitates undesirably high pressures to induce the framework phase change. As de novo design of flexible frameworks is exceedingly challenging, the ability to intuitively adapt known frameworks is required. We demonstrate that the multivariate linker approach is a powerful tool for tuning the phase change behavior of flexible frameworks. In this work, 2-methyl-5,6-difluorobenzimidazolate was solvothermally incorporated into the known framework CdIF-13 (sod-Cd(benzimidazolate) 2 ), resulting in the multivariate framework sod-Cd(benzimidazolate) 1.87 (2-methyl-5,6-difluorobenzimidazolate) 0.13 (ratio = 14:1), which exhibited a considerably reduced stepped adsorption threshold pressure while maintaining the desirable adsorption-desorption profile and capacity of CdIF-13. At 77 K, the multivariate framework exhibits stepped H 2 adsorption with saturation below 50 bar and minimal desorption hysteresis at 5 bar. At 87 K, saturation of step-shaped adsorption occurs by 90 bar, with hysteresis closing at 30 bar. These adsorption-desorption profiles enable usable capacities in a mild pressure swing process above 1 mass %, representing 85-92% of the total capacities. Here this work demonstrates that the desirable performance of flexible frameworks can be readily adapted through the multivariate approach to enable efficient storage and delivery of weakly physisorbing species.

08 HYDROGEN↗

DOE OSTI2022

A CONDUCTION-BASED HEAT PIPE MODEL FOR ANALYZING THE ENTIRE PROCESS OF LIQUID-METAL HEAT PIPE STARTUP

The Heat Pipe-cooled Microreactor (HPM) is one of the micro nuclear reactor designs under active study at the U.S. Idaho National Laboratory. Among the major concerns of HPM research is to understand the startup behavior of the heat pipe-cooled system associated with the startup of the liquid-metal heat pipes initialing from frozen state. The startup of liquid-metal heat pipes typically involves a number of nonlinear mass and heat transport processes including the phase change from solid to liquid and vapor. Hence, it is still a huge challenge to simulate the liquid-metal heat pipe startup using conventional CFD methods and software. The major difficulties of numerical CFD modeling come from the phase-change process, multiphase interaction, microporous wick flow, and compressible gas dynamics that occur during startup of the liquid-metal heat pipes. This paper proposes a simplified conduction-based method to provide practical insights into the entire startup process of the liquid-metal heat pipes while mitigating the challenges of addressing all the complex physics. We discuss the theoretical basis and modeling assumptions to analyze the liquid-metal heat pipe startup from frozen state based solely on heat-conduction equations. Then, the proposed model is implemented into the commercial CFD software to verify the model performance. The model prediction results are discussed via the comparison with the experimental data obtained from sodium heat-pipe startup experiments.

42 ENGINEERING↗

DOE OSTI2023

Lagrangian and Eulerian Supersaturation Statistics in Turbulent Cloudy Rayleigh–Bénard Convection: Applications for LES Subgrid Modeling

Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

High temperature thermochemical energy storage materials

Disclosed are high enthalpy thermochemical energy storage materials that exhibit high thermal conductivity and stability at high temperature reaction conditions. Disclosed materials include hydride-based alloys that can undergo high temperature reversible hydrogenation/dehydrogenation reactions without phase change of any metal or metalloid components of the alloy. The materials undergo a reversible exothermic hydrogenation reaction to form a metal hydride and a ternary alloy that includes a high thermal conductivity metal that, in its pure state, would exhibit a phase change at the hydrogenation reaction conditions.

Ward, Patrick A.↗

DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗

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↗

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↗

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↗

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↗