SEARCH · Search NASA
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.
This presentation covers the background, Theory, Simulation Results, Experimental Attempts, Conclusions, and Future Work for a model for simultaneous evaluation of thermal neutron scattering & resolved resonance parameters.
The Mu2e experiment requires a production target that is capable of operating under extreme thermal conditions caused by an 8 GeV proton beam. This project’s objective supports the development of the Mu2e Pro- duction Target by testing the Stickman model’s thermal behavior. Angel Flores Luviano has assisted Jonathan Williams in progressing this project by contributing to the development of a Radiative Cooling Test Fixture (RCTF) that will be used to evaluate the thermal behavior of the new production target model. Engineering calculations were performed to an- alyze thermal performance and pressure drop within the cold well and wa- ter cooling circuit. These calculations also determined the optimal sizing for key components of the water system. An engineering note was made to document these calculations. CAD models of the water circuit piping, thermocouple mounting bars, and radiator chimney were developed, and prototype test rig components were fabricated using 3D printing. Future work will focus on continued RCTF development which includes control system integration, heater hardware design, and interfaces that can be scaled up to increase thermal capacity.
In this work, thermal neutron scattering measurements of yttrium hydride from 5 to 1,200 K were conducted to determine the change in inelastic scattering as a function of temperature and to probe for anharmonic effects. Additionally, measurements on samples of YH x from x = 1.62 to 1.90 were done to determine the effects of varying hydrogen concentration on the inelastic spectra. Changes in temperature affected the inelastic spectra in unanticipated ways, indicating that there are anharmonic effects, whereas hydrogen concentration did not significantly affect the inelastic spectra. These measurements were compared against the ENDF/B-VIII.0 thermal scattering files of YH 2 , as well as new thermal scattering files created by using the stochastic temperature-dependent effective potential (s-TDEP). Both libraries were found to be in good agreement with the experimental data at lower temperatures. At higher temperatures, the s-TDEP method is better at predicting the experimentally observed softening of phonon modes in acoustic and optical regions.
Uranium Dioxide (UO 2 ) fuel powers almost all commercial Nuclear Power Plants (NPPs) worldwide, generating carbon-free energy and contributing to the fight against climate change. UO 2 fuel incurs damage and fractures due to large thermal gradients that develop across the fuel pellet during normal and transient operating conditions. A comprehensive understanding of the underlying mechanisms by which these processes take place is still lacking. A combined experimental and computational approach is utilized here to quantify the behavior of UO 2 fuel fracture induced by thermal shock. Here, this work introduces both (1) an experimental study to understand the fuel fracturing behavior of sintered UO 2 pellets when exposed to thermal shock, and (2) a Multiphysics phase-field fracture model capable of simulating this process. Parametric studies were conducted to evaluate the effects of uncertainties in fracture properties on the fracture behavior of UO 2 due to thermal shocking. A set of energy release rate (or equivalently fracture toughness) and contract area (the part of the fuel pellet in direct contact with the cold bath) were able to capture the overall fracture trends of the corresponding experimental data. Our combined approach presents a new method for accounting for the effects of microstructure and sample size on the energy release rate/fracture toughness. The experimental data were collected from multiple experiments that exposed UO 2 pellets to high-temperature conditions (589–676 °C) followed by a quench in sub-zero water. This work demonstrates that joint experimental and computational efforts are able to advance the understanding of thermal fracture in the primary fuel source for existing and future NPPs.
In uranium-zirconium (U-Zr) based metallic fuels, different phases can form at different compositions and temperatures. Typically, lamellar δ-UZr 2 and α-U phases are the dominant microstructures in U-rich U-Zr alloys at temperatures below 880 K. In this work, a finite element method based mesoscale modeling technique is used to calculate the effective thermal conductivities of such heterogeneous microstructures, using the thermal conductivities of two individual phases and their interphase thermal resistance (Kapitza resistance) as input parameters. The Kapitza resistance between δ-UZr 2 and α-U is determined at different temperatures, which shows an approximately T 3 dependence in the temperature range between 300 and 800 K. In addition, the Kapitza resistance exhibits a strong dependence on the aspect ratio of the δ-UZr 2 phase. Further, an analytical model is therefore developed to quantify the effects of both temperature and δ-UZr 2 aspect ratio on the Kapitza resistance. Using this newly developed Kapitza resistance model, the effective thermal conductivities of a number of δ-UZr 2 + α-U heterogeneous microstructures in U-Zr alloys, including non-lamellar microstructures, can be estimated accurately.
A storage heat pump water heater (HPWH) comprises a vapor compression system and a hot water tank. It can lead to more than 70% energy saving, compared with electric resistance heating. HPWH usually uses a wrapped-tank coil as the water heater. The wrapped coil has a strong interaction with the water tank and impose a new challenge for the engineering design and simulation study. Based on the US DOE/ORNL Heat Pump Design Model (HPDM), i.e., a public-domain, component-based thermal system modelling platform, a segment-to-segment wrapped-tank coil model coupled with a one-dimensional stratified water tank model were developed, which facilitates modeling and design HPWHs having various configurations and components. Additionally, the hardware-based design model was used to evaluate refrigerants having lower global warming potentials to replace R-134a, including R-1234yf, R-1234ze, R-290, R-513A and R-450A, when effects of compressor size, refrigerant charge on 24-hour uniform energy factor (UEF) and first-hour rating (FHR) were investigated.In the drop-in replacements, R-1234ze achieved the best UEF, and R-290 resulted in the lowest UEF. When downsizing the compressor to maximize the UEFs, the UEFs of R-134a, R-290, R1234yf, R-450A and R-513A ascended to the range of 4.4 to 4.53, with R-513A having the highest UEF. The R-1234ze ended up having the lowest UEF of 4.27. Refrigerants having larger volumetric vaporization heat led to higher FHRs if using the fixed-size compressor, that R-290 resulted in the highest FHR.
Emerging technologies like hypersonic aircraft, space exploration vehicles, and batteries avail fluid circulation in embedded microvasculatures for efficient thermal regulation. Modeling is vital during the design and operational phases of these engineered systems. However, many challenges exist in developing a modeling framework. What is lacking is an accurate framework that (i) captures sharp jumps in the thermal flux across complex vasculature layouts, (ii) deals with oblique derivatives (involving tangential and normal components), (iii) handles nonlinearity because of radiative heat transfer, (iv) provides a high-speed forecast for real-time monitoring, and (v) facilitates robust inverse modeling. Here, this paper addresses these challenges by availing the power of physics-informed neural networks (PINNs). We develop a fast, reliable, and accurate Scientific Machine Learning (SciML) framework for vascular-based thermal regulation—called CoolPINNs: a PINNs-based modeling framework for active cooling. The proposed mesh-less framework elegantly overcomes all the mentioned challenges. The significance of the reported research is multi-fold. First, the framework is valuable for real-time monitoring of thermal regulatory systems because of rapid forecasting. Second, researchers can address complex thermoregulation designs since the approach is meshless. Finally, the framework facilitates systematic parameter identification and inverse modeling studies, perhaps the most significant utility of the current framework.
Small planets on close-in orbits tend to exhibit envelope mass fractions of either effectively zero or up to a few percent depending on their size and orbital period. Models of thermally driven atmospheric mass loss and of terrestrial planet formation in a gas-poor environment make distinct predictions regarding the location of this rocky/nonrocky transition in period–radius space. Here we present the confirmation of TOI-1235 b (P = 3.44 days, r{sub p}=1.738{sub −0.076}{sup +0.087} R{sub ⊕}), a planet whose size and period are intermediate between the competing model predictions, thus making the system an important test case for emergence models of the rocky/nonrocky transition around early M dwarfs (R {sub s} = 0.630 ± 0.015 R{sub ⊙}, M {sub s} = 0.640 ± 0.016 M{sub ⊙}). We confirm the TESS planet discovery using reconnaissance spectroscopy, ground-based photometry, high-resolution imaging, and a set of 38 precise radial velocities (RVs) from HARPS-N and HIRES. We measure a planet mass of 6.91{sub −0.85}{sup +0.75} M{sub ⊕}, which implies an iron core mass fraction of 20{sub −12}{sup +15}% in the absence of a gaseous envelope. The bulk composition of TOI-1235 b is therefore consistent with being Earth-like, and we constrain an H/He envelope mass fraction to be <0.5% at 90% confidence. Our results are consistent with model predictions from thermally driven atmospheric mass loss but not with gas-poor formation, suggesting that the former class of processes remains efficient at sculpting close-in planets around early M dwarfs. Our RV analysis also reveals a strong periodicity close to the first harmonic of the photometrically determined stellar rotation period that we treat as stellar activity, despite other lines of evidence favoring a planetary origin (P=21.8{sub −0.8}{sup +0.9} days, m{sub p}sin i=13.0{sub −5.3}{sup +3.8} M{sub ⊕}) that cannot be firmly ruled out by our data.
As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.
This project developed stand-alone electric-thermal energy storage (ETES), stand-alone pumped thermal energy storage (PTES), and hybrid molten-salt power tower (MSPT)-ETES performance and dispatch optimization techno-economic models. The models are available to the public through the System Advisor Model (SAM) software, scripting, and as open-source code. We compared results of the dispatch model to PLEXOS dispatch of a similar generator using the same initial grid pricing signal and found our dispatch model performed well, but closer agreement between the models was limited by the inherent differences between price-taker and unit commitment models. Nevertheless, the price-taker models developed in this project are useful to analyze proposed ETES and PTES technologies because they provide more detailed system and component models, solve several orders of magnitude faster, and are available as free open-source software. The model results represent the most optimistic returns considering grid arbitrage from the input electricity pricing, so the financial results can be applied as a feasibility stage-gate. We also submitted a journal draft paper that describes the ETES dispatch model methodology and demonstrates model functionality. An accepted journal article will serve as peer-reviewed documentation for the models, along with the open-source code, SAM help-menu content, and eventually this final project report.
Existing pavement de-icing methods result in high installation and maintenance costs, traffic delays, excessive weariness and corrosion, and negative environmental and safety impacts. To overcome these challenges, this paper presents an innovative pathway to designing and constructing smart self-heating pavements using a low-cost coal-char bearing asphalt material. The conductive asphalt incorporates coal char, i.e., a key byproduct of the coal pyrolysis process, into the Stone Mastic Asphalt (SMA) mixture. This asphalt containing coal-derived solid carbon exhibits highly tailorable electrical conductivity, satisfactory mechanical and thermophysical properties, and superior cost-efficiency as compared to other conductive pavement materials. The de-icing performance was also demonstrated by laboratory experiments on a bench-scale prototype. Furthermore, an efficient thermal network model was developed and validated by experiments to investigate the transient thermal behavior and energy performance of the Ohmic heating pavement system. Furthermore, the whole-year energy simulation case studies were conducted on a bridge pavement with an annual energy use of 24.6–1444.8 kWh/m 2 , showcasing its potential in field applications across cool humid, cold humid, and subarctic/arctic climate zones.
Objective performance and economic modeling of solar thermal plants is of keen interest to many EPRI funders. One solar thermal technology that is not currently available to model in any non-vendor, non-proprietary tool is linear Fresnel. This technology has garnered enough interest from EPRI funders to merit investing in a tool to objectively model its performance. Early in 2010, EPRI performed a comparison of modeling solar thermal power plants using the IPSEPRO, CNRS and Solar Advisor (SAM) tools. After completing this effort, EPRI decided to adopt NREL’s Solar Advisor Model as its default modeling tool based in part on user friendliness, flexibility, number of technologies covered, integrated financial model and ease of running sensitivities. Furthermore, it was recognized that NREL continues to invest considerable time and resources into improving capabilities and functionality of the model.
Due to heat prosumers' dual roles of heat producer and heat consumer, the future district heating (DH) systems will become more flexible and competitive. However, the current heating price models have not yet supported the reverse heat supply from prosumers to the central DH system, which means the prosumers would gain no economic benefit from supplying heat to the central DH system. These unidirectional heating price models will reduce interest in prosumers, and thus hinder the promotion of prosumers in DH systems. This study aimed to optimize prosumers' economic performance under the current heating price models by introducing water tank thermal energy storage (WTTES). A dynamic optimization problem was formulated to explore prosumers' economic potentials. The size parameter of WTTESs was swept in prosumers to obtain the optimal storage size considering the trade-off between the payback period and the heating cost saving. The proposed method was tested on a campus DH system in Norway. The results showed that the prosumer's annual heating cost was saved up to 9%, and the investment of WTTES could be recovered in less than ten years. This study could provide guidelines on improving prosumers' economic performance and promote the development of prosumers during the transformation period of DH systems.
To push upper boundaries of thermal conductivity in polymer composites, understanding of thermal transport mechanisms is crucial. Despite extensive simulations, systematic experimental investigation on thermal transport in polymer composites is limited. To better understand thermal transport processes, we design polymer composites with perfect fillers (graphite) and defective fillers (graphite oxide), using polyvinyl alcohol (PVA) as a matrix model. Measured thermal conductivities of ~1.38 ± 0.22 W m -1 K -1 in PVA/defective filler composites is higher than those of ~0.86 ± 0.21 W m -1 K -1 in PVA/perfect filler composites, while measured thermal conductivities in defective fillers are lower than those of perfect fillers. We identify how thermal transport occurs across heterogeneous interfaces. Thermal transport measurements, neutron scattering, quantum mechanical modeling, and molecular dynamics simulations reveal that vibrational coupling between PVA and defective fillers at PVA/filler interfaces enhances thermal conductivity, suggesting that defects in polymer composites improve thermal transport by promoting this vibrational coupling.
Understanding thermal transport in solid electrolytes is essential for improving the performance, reliability, and safety of all-solid-state batteries. Garnet-type lithium-ion conductors are promising candidates for solid electrolytes, yet their thermal-transport mechanisms remain poorly understood. Here, we connect the lattice and ion dynamics of single-crystal garnet-type Li 6.5 La 3 Zr 1.5 Ta 0.5 O 12 to its intrinsically low thermal conductivity. Our study reveals that the single crystals grown by the floating-zone method exhibit remarkably low glasslike thermal conductivity. Using first-principles calculations and inelastic-neutron-scattering measurements, we identify both the acoustic and numerous optical phonon modes, which stem from the complex crystal structure of the material. Notably, a low-energy optical branch exhibits an avoided crossing with acoustic phonons near 7 meV. These optical modes can enhance the scattering of heat-carrying acoustic phonons and reduce thermal conductivity. Furthermore, the calculated Grüneisen parameters are large, especially for the vibrational modes around 6 meV, indicating strong anharmonicity, with a noticeable contribution from lithium-ion vibrations. A two-channel thermal-transport model is employed to describe the weak temperature dependence of the thermal conductivity, which can be attributed to the substantial contribution of diffuson transport facilitated by the abundance of optical phonons and intrinsic anharmonicity. These results offer valuable insights into the thermal transport in a broad class of ionic conductors of interest for energy conversion and storage applications.
Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λTD−DFT(ρ,T,K), was devised for CH plasmas between ρ=[0.05−1.05] g/cm3, kBT=[100−1000] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations.
This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.
The Pacific Northwest National Laboratory, Oak Ridge National Laboratory and the University of Minnesota are conducting a three-year, multipart study on residential retrofit wall assemblies. The project is funded by the U.S. Department of Energy’s Building Technology Office and will identify, test and verify wall assemblies for hygrothermal performance in retrofit applications. The study includes a comprehensive literature review and expert advisory group, which inform wall selection. Selected wall assemblies are then thermally simulated using EnergyPlus and THERM, and hygrothermally simulated using WUFI, to model both thermal and moisture performance of the wall assembly. Eight wall assemblies are then experimentally tested in an in-situ laboratory environment at the University of Minnesota, with a typical residential wall used as a baseline. The in-situ experiment looks at the physical hygrothermal performance of each assembly. The simulation and experimental results will be combined with an economic analysis to produce a techno-economic study of residential wall systems for deep energy retrofits. This paper presents progress after year one of the study. Wall system modeling and laboratory testing is scheduled to begin in November, 2019. We will discuss the parameters for wall selection and experimental design, which includes approaches to modeling and simulation, along with the physical experiment design. This paper will also outline the parameters for the techno-economic analysis, and criteria used to inform model development.