SEARCH · Engineering Papers
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.
Seminar presentation to the Ohio State University Nuclear Engineering graduate program. These slides present an overview of work performed at INL in support of NASA nuclear thermal propulsion modeling and simulation.
Vibrational thermal properties of CuZn 2 InTe 4 , AgZn 2 InTe 4 , and Cu 2 CdSnTe 4 , derived from binary II-VI zinc-blendes, are reported based on first-principles calculations. While the chalcogenide atoms in these materials have the same lattice positions, the cation atom arrangements vary, resulting in different crystal symmetries and subsequent properties. The compositional differences have important effects on the vibrational thermal characteristics of the studied materials, which demonstrate that low-frequency optical phonons hybridize with acoustic phonons and lead to enhanced phonon-phonon scattering and low lattice thermal conductivities. The phonon density of states, mode Grüneisen parameters, and phonon scattering rates are also calculated, enabling deeper insight into the microscopic thermal conduction processes in these materials. Compositional variations drive differences among the three materials considered here; nonetheless, their structural similarities and generally low thermal conductivities (0.5–4 W/m K at room temperature) suggest that other similar II-VI zinc-blende derived materials will also exhibit similarly low values, as also corroborated by experimental data. Finally, this, combined with the versatility in designing a variety of motifs on the overall structure, makes quaternary chalcogenides interesting for thermal management and energy conversion applications that require low thermal conductivity.
Global gyrokinetic simulations of mesoscale reversed shear Alfven eigenmodes (RSAE) excited by energetic particles (EP) in fusion plasmas find that RSAE amplitude and EP transport are much higher than experimental levels at nonlinear saturation, but quickly diminish to very low levels after the saturation when background microturbulence is artificially suppressed. In contrast, in simulations coupling micro-meso scales, the RSAE amplitude and EP transport decrease drastically at the initial saturation but later increases to the experimental levels in the quasi-steady state with bursty dynamics due to regulation by thermal ion temperature gradient (ITG) microturbulence. The quasi-steady state EP transport is larger for a stronger microturbulence. Here, the RSAE amplitude in the quasi-steady state ITG-RSAE turbulence from gyrokinetic simulations, for the first time, agrees very well with experimental measurements.
Thermal energy storage can increase the efficiency of the electric grid by adding flexibility to thermal systems. The value of thermal storage is a function of its energy and power density, which are driven by the capacity and thermal resistances in the storage device. Measuring these properties in-situ at the device level is an important step to understanding the performance and improving the design of thermal storage systems. Here, we present methods to measure the total capacity and thermal resistances in heat exchangers with integrated phase change materials. These methods are demonstrated on two thermal storage devices - a 570-kWh ice-based storage tank and a 0.35-kWh graphite-tetradecane composite device. The results show how thermal resistances evolve with the state of charge and discharge rate in these devices and quantify the impact of applied pressure on the contact resistance in composite phase change material heat exchangers. The proposed method allows for easy comparison between different systems and provides information on the thermal bottlenecks limiting performance. Ultimately, these measurements will allow designers to make robust, high-performance thermal storage devices for next-generation thermal systems.
Here, a new equation of state for 304L stainless steel has been constructed and entered in the Los Alamos National Laboratory SESAME database with material ID 4273. The equation of state is based on a multiphase approach. The solid and liquid phases were fitted to available experimental data for thermal expansion, heat capacity, and shock Hugoniot, as well as to data from density functional theory calculations of the cold curve and quantum molecular dynamics of the liquid regime. The electronic thermal contribution over the entire equation of state range is based on Tartarus (Green’s functions) average atom calculations. For this low-carbon stainless steel both the density functional and electronic thermal calculations are comprised of 19 wt% Cr, 10 wt% Ni, and the remainder Fe.
Neptunium (Np) possesses a rich and unique chemistry that often diverges from other actinide elements yet remains relatively underexplored compared with the other light actinides. A resurgence of interest in Np has been spurred by the application of 237 Np for plutonium-238 ( 238 Pu) production for use in radioisotope thermoelectric generators (RTGs), necessitating evaluation of Np chemical reactions and materials. The work presented here studied the thermal decomposition of neptunyl ammonium nitrate (NH 4 Np VI O 2 (NO 3 ) 3 ) for synthesis of neptunium dioxide (NpO 2 ), which is the target material used for production of 238 Pu. Additionally, structural characterization of the intermediate solid Np pentoxide (Np 2 O 5 ) was performed. Advanced solid-state characterization techniques, including simultaneous thermal analysis (STA), powder X-ray diffraction (pXRD), Raman spectroscopy, and density functional theory (DFT) modeling have been combined to study the reaction pathways. Analysis revealed that NH 4 Np VI O 2 (NO 3 ) 3 thermally decomposes to a proposed neptunyl nitrate intermediate, followed by Np 2 O 5 and finally NpO 2 , all within the temperature range of 150 °C–600 °C. Further characterization of the pentoxide intermediate provided the first Raman spectra of pure-phase Np 2 O 5 and associated DFT modeling confirmed Raman peak assignments for this phase. These findings provide mechanistic information to advance production of the critical radioisotope 238 Pu and advance the state of knowledge on Np materials chemistry using modern characterization techniques.
Understanding the dynamics of bound state formation is one of the fundamental questions in confining quantum field theories such as Quantum Chromodynamics (QCD). One hadronization mechanism that has garnered significant attention is the breaking of a string initially connecting a fermion and an antifermion. Deepening our understanding of real-time string-breaking dynamics with simpler, lower dimensional models like the Schwinger model can improve our understanding of the hadronization process in QCD and other confining systems found in condensed matter and statistical systems. Here in this paper, we consider the string-breaking dynamics within the Schwinger model and investigate its modification inside a thermal medium, treating the Schwinger model as an open quantum system coupled to a thermal environment. Within the regime of weak coupling between the system and environment, the real-time evolution of the system can be described by a Lindblad evolution equation. We analyze the Liouvillian gaps of this Lindblad equation and the time dependence of the system’s von Neumann entropy. We observe that the late-time relaxation rate decreases as the environment correlation length increases. Moreover, when the environment correlation length is infinite, the system exhibits two steady states, one in each of the sectors with definite charge-conjugation-parity (CP) quantum numbers. For parameter regimes where an initial string breaks in vacuum, we observe a delay of the string breaking in the medium, due to kinetic dissipation effects. Conversely, in regimes where an initial string remains intact in vacuum time evolution, we observe string breaking (melting) in the thermal medium. We further discuss how the Liouvillian dynamics of the open Schwinger model can be simulated on quantum computers and provide an estimate of the associated Trotter errors.
We report density functional theory has revolutionised our understanding of dense, strongly coupled systems in the condensed-matter limit. At elevated temperatures, such as those found in warm dense matter and dense plasma, the computational cost of calculating enough states to accurately describe a thermal ensemble of electrons is a major obstacle to practical Kohn-Sham calculations. This problem has recently been tackled from a number of approaches, including a plane-wave approximation for electrons at high energy. In this work, we demonstrate how the Thomas-Fermi density of states, equivalent to a local plane-wave approximation, can yield improved results for thermodynamic quantities, electron density profiles, and interatomic forces while requiring fewer states than existing methods.
The strongly correlated material FeSi displays several unusual thermal, magnetic, and structural properties under varying pressure-temperature (P-T) conditions. It is a potential thermoelectric alloy and a material with several geochemical implications as a possible constituent at the Earth's core-mantle boundary (CMB). Previous theoretical studies predicted a pressure-induced B20-B2 transition at ambient temperature below 40 GPa; however, experimentally, the structural transition is observed only at high P-T conditions. In this study, we have performed high-pressure powder X-ray diffraction (XRD) up to 90 GPa and Nuclear Resonant Inelastic X-ray Scattering (NRIXS) measurements up to 120 GPa to understand the phase stability and lattice dynamics. Our study provides evidence for a non-hydrostatic stress-induced B20-B2 transition in FeSi at around 36 GPa. We deduced the Fe partial phonon density of states (PDOS) and thermal parameters from NRIXS measurements up to 120 GPa and compared them with density functional theory (DFT) calculations. Furthermore, the computations show pressure-induced metallization and the band gap closing around 12 GPa.
The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.
Thermal conductivity (TC) data taken for different wells at a specified drill depth. This is an abridged version of the complete SMU heat flow database, downloaded from the SMU node of the NGDS at the beginning of INGENIOUS (approximately April 2021), and filtered to the INGENIOUS study area. This National Geothermal Data System (NGDS) project aggregates geothermal data collected and curated by the SMU Geothermal Laboratory and its partner organizations. All columns in this database are the same as the SMU database, except for 2 additions associated with this project. Repeated columns are for data correlation purposes. Column descriptions and data types are the same as previous iterations of the SMU database. The new values that are the addition are two new columns developed as part of the INGENIOUS project: INGENIOUS TC Value | INGENIOUS notes INGENIOUS notes are individual notes that were written for specific data points during the analysis process. There are not always notes associated with each input value. INGENIOUS TC Value includes 4 values: 1. Assumed Measured These are values that are assumed to be measured thermal conductivity values, either within a specific well or within the same study region. Many of these have either a published reference, a reported standard deviation, or a unique thermal conductivity value. 2. Data release - assumed measured These are values in the SMU database that are from proprietary data that were added to the SMU database and are labeled as data release for their reference. These values were searched for in person at the SMU Geothermal Laboratory as well as virtual examination of data available on the NGDS. For many of these, there are reported thermal conductivity values associated with the heat flow data in the database, but no specific table or reference to measurements in the original data release files. 3. Known measured These are values that have a reported measurement, either as an original file in the SMU data files on the NGDS or a reported table in a publication. In the rare circumstances, Maria Richards or David Blackwell confirmed measurement. Confirmation of measurement would be written in the INGENIOUS notes column. 4. Unmeasured Unmeasured values are those that are known to be unmeasured, either estimated from another report or no information given. In the SMU database, there are wells that have a heat flow but no thermal conductivity. These are categorized as unmeasured. There are also heat flow values that are stated to have estimated or generalized average thermal conductivity values for the region and rock type. Because these are known to be unmeasured, they are categorized as such. 5. Blank Blank values are either A quality or X quality. These quality values are stated in the INGENIOUS notes. These values were not going to change associated with the heat flow analysis, so these were not examined.
Stuffed aluminosilicate quartz solid solutions (Qss) represent the main functional component of state-of-the-art zero-thermal-expansion glass-ceramics. For the first time, we present the synthesis of Li-, Mg- and Zn-free Qss solely stuffed with transition metal cations (Fe 2+ , Ni 2+ , Co 2+ ); partial Li + co-doping enables also significant Mn 2+ incorporation. Further, they were obtained by glass powder crystallization; Qss crystals exhibit compositionally tunable coefficients of thermal expansion (CTEs), from - 2 × 10 -6 K -1 to ~10 × 10 -6 K -1 in the range 30–300 °C. Co 2+ -bearing crystals exhibit the closest-to-zero CTE value (0.2 ×10 -6 K -1 ) and the most isotropic behavior in the whole Qss family, opening up new perspectives for the development of Li-free low-expansion materials. From a structural point of view, we identified the unit cell volume (and not the pseudo-hexagonality of the aluminosilicate framework) as the key structural parameter leading to low or negative CTEs in Qss, with a linear correlation extending to non-stuffed non-silicate quartz-like phases.
Nonequilibrium hidden states provide a unique window into thermally inaccessible regimes of strong coupling between microscopic degrees of freedom in quantum materials. Understanding the origin of these states allows the exploration of far-from-equilibrium thermodynamics and the development of optoelectronic devices with on-demand photoresponses. However, mapping the ultrafast formation of a long-lived hidden phase remains a longstanding challenge since the initial state is not recovered rapidly. Here, using state-of-the-art single-shot spectroscopy techniques, we present a direct ultrafast visualization of the photoinduced phase transition to both transient and long-lived hidden states in an electronic crystal, 1T-TaS 2 , and demonstrate a commonality in their microscopic pathways, driven by the collapse of charge order. We present a theory of fluctuation-dominated process that helps explain the nature of the metastable state. Our results shed light on the origin of this elusive state and pave the way for the discovery of other exotic phases of matter.
The widespread adoption of electric-driven heat pumps for heating and cooling is expected to significantly increase electric demand. Cooling electric demand will result in an increase in peak hours electric loads, placing additional strain on the grid during peak hours. This challenge, combined with the current rapidly increasing demand from data centers, exacerbates the electric demand duck curve problem. Integrating thermal energy storage (TES) with heat pump can shift electric use for heating and cooling from peak to off-peak hours of the electric grid, which can help flatten the daily electric demand profile of a building. This paper presents a novel design that integrates heat pump with TES (HP-TES), which uses phase change materials (PCM). Heat pump charges TES by melting or freezing PCM during off-peak hours when there is no thermal demand from the building. TES is then discharged (i.e., by freezing or melting PCM) during peak hours to provide a more favorable heat source or heat sink for the heat pump to meet the thermal demand of the building with lower electricity use than conventional air-source heat pumps. Computer simulations were developed to predict the performance of HP-TES applied to a typical single-family house in the US. The building-level simulation results were scaled up to preliminarily assess the aggregated impacts of deploying HP-TES across all single-family houses in Texas of the United States, including reduction of peak demand of the electric grid.
The hot box test method has been applied to evaluate both the steady-state (U-value) and dynamic thermal properties of building envelopes. However, the high construction cost of full-scale hot box apparatus and the testing time required (usually several days) may prevent its wider adoption. To overcome the limitations of full-scale hot box tests, here we propose a novel method to evaluate the dynamic thermal performance of building envelope components using a small-scale calibrated hot box and scaled-down specimen. In this paper, the scaling relationships of thermal properties evaluated using a full-size specimen and a scaled-down specimen are established based on the Laplace transform of the heat transfer equations. In addition, dynamic thermal properties obtained from scaled-down experimental tests are compared to the values calculated by the EN ISO13768 (ISO) method. A small-scale hot box with a 355 mm × 355 mm metering area was constructed and calibrated to validate the correlations. Three scaled-down concrete sandwich wall panels were then tested and the scaling relationship was cross-validated using the experimental results, finite difference (FD) simulations, and the ISO method. The results indicate that the dynamic thermal properties obtained from a scaled-down hot box test can be correlated to its full-size counterpart when certain conditions are met. The scaled-down hot box test is demonstrated to be an effective yet economical alternative to a full-scale test with significantly reduced experimentation cost and turn-around time.
HRL Laboratories and General Motors (GM), pursued this effort to increase the efficiency of internal combustion engines (ICEs) by developing and implementing temperature-following thermal barrier coatings (TBCs) to decrease heat loss from the combustion chamber. HRL developed microsphere TBCs based on hollow high-temperature alloy spheres with average diameter of 50 μm and wall thickness ~1 μm (microspheres) that achieve unprecedented low thermal conductivity and heat capacity (10X lower than the state of the art) while offering exceptional environmental and mechanical resistance.> In current ICEs, approximately 29% of the fuel’s energy is lost to the cooling system and about 22% goes into moving the car. By limiting heat losses from the combustion chamber with insulating coatings, fuel energy can be redirected into additional piston work and into the exhaust stream. The additional energy that goes into the exhaust stream can be turned into useful work through the use of an exhaust energy recovery device, such as turbocompounding and/or through driving a turbocharger to increase the power density of the engine allowing downsizing. Insulating coatings for piston crowns have been pursued in the past, but previous materials—typically ceramics—exhibited low thermal conductivity but retained high heat capacity. Such materials reduce heat losses but stabilize at a high surface temperature. The high surface temperature heats the incoming gases, which lowers volumetric efficiency and increases propensity for knock, resulting in degraded engine performance. This program developed an innovative new material that combines low thermal conductivity with low heat capacity. These unique properties allow it to follow rapid changes in gas temperature during each combustion cycle. A metallic microsphere TBC has been demonstrated that exhibits increased surface temperature during the combustion period, resulting in reduced heat transfer losses, while still returning to a low surface temperature during the gas exchange period. A 10X lower thermal conductivity and heat capacity than state-of-the-art thermal barrier coatings has been demonstrated. These unprecedented thermal properties are achieved through the coating architecture, which consists of closed pores on the microscale and exhibits a total porosity of 90% to 95%. By selecting a high-temperature nickel alloy as the coating material, the ductility and strength of metals can be harnessed to achieve much better durability and damage tolerance than brittle ceramic coatings can achieve.
Additive manufacturing (AM) is known to generate large magnitudes of residual stresses (RS) within builds due to steep and localized thermal gradients. In the current state of commercial AM technology, manufacturers generally perform heat treatments in effort to reduce the generated RS and its detrimental effects on part distortion and in-service failure. Computational models that effectively simulate the deposition process can provide valuable insights to improve RS distributions. Accordingly, it is common to employ Computational fluid dynamics (CFD) models or finite element (FE) models. While CFD can predict geometric and thermal-fluid behavior, it cannot predict the structural response (e.g., stress–strain) behavior. On the other hand, an FE model can predict mechanical behavior, but it lacks the ability to predict geometric and fluid behavior. Thus, an effectively integrated thermofluidic–thermomechanical modeling framework that exploits the benefits of both techniques while avoiding their respective limitations can offer valuable predictive capability for AM processes. In contrast to previously published efforts, the work herein describes a one-way coupled CFD-FEA framework that abandons major simplifying assumptions, such as geometric steady-state conditions, the absence of material plasticity, and the lack of detailed RS evolution/accumulation during deposition, as well as insufficient validation of results. Here, the presented framework is demonstrated for a directed energy deposition (DED) process, and experiments are performed to validate the predicted geometry and RS profile. Both single- and double-layer stainless steel 316L builds are considered. Geometric data is acquired via 3D optical surface scans and X-ray micro-computed tomography, and residual stress is measured using neutron diffraction (ND). Comparisons between the simulations and measurements reveal that the described CFD-FEA framework is effective in capturing the coupled thermomechanical and thermofluidic behaviors of the DED process. The methodology presented is extensible to other metal AM processes, including power bed fusion and wire-feed-based AM.
The United State High Performance Research Reactor (USHPRR) program aims to eliminate more than 200kg of High Enriched Uranium (HEU) from commerce annually by converting five U.S. high-performance research reactors and one associated critical assembly to Low Enriched Uranium (LEU) fuel using a high-density alloy of uranium-10 wt% molybdenum (U-10Mo). Massachusetts Institute of Technology Reactor (MITR) is one of five research reactors selected for this program. The objective of this report is to provide preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR Design Demonstration Element (DDE) under conservative approximations for the plate power distribution. This report provides details on the modeling approach and the simulation results obtained, including pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment for the irradiation cycles in the Belgian Reactor (BR)-2.