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.
This study examines the intricate area of refractory-based high entropy alloys (RHEAs), focusing on a series of complex compositions involving nine diverse refractory elements: Ti, V, Cr, Zr, Nb, Mo, Hf, Ta, and W.
Nb-based MXenes, particularly Nb 2 CT X is noted for its metallic characteristics, and nearly zero band gap because of it 2D structure. Different synthesis processes have been utilized such as HF etching, HCl/LiF and molten salt method to prepare delaminated Nb 2 CT X MXene sheets[1]. Moreover, various approaches have been tried for modifying MXenes properties including partial oxidation for application such as energy storage, gas sensing, photocatalyst etc. One step synthesis technique used for creating efficient hydrogen evolution photocatalyst comprised of Nb 2 O 5 /C/Nb 2 C composites. Furthermore, we have prepared Nb 2 CT X MXenes with a strong etching method and investigated its fundamental properties at elevated temperature.
INTRODUCTION Cementitious matrices are used in the US DOE complex and worldwide to solidify aqueous radioactive, hazardous, mixed salt solutions, and sludges to meet low-level radioactive waste (LLW) disposal requirements. Blended formulations are used for most waste form mix designs. Substitution of pozzolans for Class F fly ash has the potential to alter the processing properties and stabilization properties of the resulting waste forms because they add chemical and mineralogical complexity to the final material. The Savannah River Site saltstone waste form was selected as the test case material. CemGEMS software was chosen as the model for predicting hydrated cementitious phases assemblages and phase evolution. Isothermal calorimetry was used as the method for evaluating processing and properties of fresh, uncured waste forms. The test cases consisted of reference case saltstone, and five natural pozzolan substituted saltstone mixes.
ASME-SHCT presentation ~ 15 mins
Efficient and accurate ordinary differential equation (ODE) solvers are necessary for powertrain and vehicle dynamics modeling. However, current commercial ODE solvers can be financially prohibitive, leading to a need for accessible, effective, open-source ODE solvers designed for powertrain modeling. Rust is a compiled programming language that has the potential to be used for fast and easy-to-use powertrain models, given its exceptional computational performance, robust package ecosystem, and short time required for modelers to become proficient. However, of the three commonly used (>3,000 downloads) packages in Rust with ODE solver capabilities, only one has more than four numerical methods implemented, and none are designed specifically for modeling physical systems. Therefore, the goal of the Differential Equation System Solver (DESS) was to implement accurate ODE solvers in Rust designed for the component-based problems often seen in powertrain modeling. DESS is a text-based software package that provides a flexible framework for building and solving systems of ODEs. This allows DESS to be included as a dependency for automotive powertrain models that require a variety of solvers and solver configurations. Seven explicit ODE solver methods have been implemented in DESS: Euler’s, Heun’s, midpoint, Ralston’s, classic Runge-Kutta, Bogacki-Shampine, and Cash-Karp. These represent five fixed-step methods and two adaptive-step methods. This paper shows that the solver implementations increase accuracy and computational efficiency compared to Euler's method when modeling a system of three thermal masses in Rust. DESS also includes features designed for modeling component-based physical systems. Users can define relationships between nodes in their system, which the package then translates into a system of equations, leading to simpler and more intuitive code. In the case of a three-thermal-mass system, the user can specify node thermal properties (e.g., thermal capacitance), how nodes are interconnected, and thermal conductance between nodes rather than providing a system of equations. The core contribution from this work is an open-source, text-based Rust package with ODE solvers for automotive powertrain modeling to support cost-free, fast, and accurate simulation.
Laser Thomson scattering, which provides thermal properties of electrons, has been extensively utilized in plasma diagnostics. However, being a non-resonant linear light scattering technique, it often encounters several challenges in weakly-ionized plasma diagnostics due to spectral overlap with different scatterings induced from other existing species, including Mie, Rayleigh, and rotational Raman scatterings. To address this challenge, we propose an imaging spectroscopy technique, polarization-separated double-imaging spectroscopy (PoDIS), which selectively separates specific scatterings based on their polarization characteristics. Using an atmospheric plasma jet as a plasma source for demonstration, we show that PoDIS can effectively separate rotational Raman and Thomson scatterings from a superimposed spectrum without prior knowledge or assumptions about the thermal properties of neutrals or electrons. This separation enables independent fitting of the rotational Raman and Thomson scattering spectra, allowing for precise determination of the thermal properties of neutral particles and electrons separately.
In this project, a novel concept of Octet-shape based lattice frame material was developed for deployment in particle-to-sCO2 heat exchangers for CSP application. A comprehensive experimental and numerical program was performed to investigate and characterize the thermal and flow transport behavior in high porosity lattice frame material. The thermal properties include effective thermal conductivity and interstitial heat transfer coefficient and flow properties including permeability and inertial coefficient. Above quantities were determined for both air flow and heated particle flow through the porous channels which were representative of the hot channel in a typical counter flow heat exchanger. The cold side of the heat exchanger feature minichannels which are known to be very efficient in thermal transport. The goal of this project was to enhance the overall heat transfer coefficient (Up) on the hot side by introducing Octet-shape based lattice frame material through enhancement in effective thermal conductivity of the channel and through enhancement in the interstitial heat transfer coefficient.
Since 2011, a geothermal heat pump (GHP) system has been operating to provide space heating and cooling for the Solar Radiation and Research Laboratory building at the National Laboratory of the Rockies (NLR) in Golden, Colorado. The system consists of 23 vertical boreholes, each extending to a depth of 300 ft (91 m), connected to 11 water-to-air heat pump units and four circulation pumps. Between fiscal years 2023 and 2025, additional power meters and temperature sensors were retrofitted to support detailed system performance assessment and model development. This study presents preliminary monitoring results and the development of an initial numerical model of the borehole heat exchanger field. The model incorporated site-specific geometry, ground thermal properties derived from thermal response tests, and ambient temperatures, and simulated system behavior over a representative operating day in September. Model predictions of outlet temperatures were compared against corresponding field measurements. Results showed that modeling initialized with a simplified linear subsurface temperature gradient presents systematic discrepancies in outlet temperature, whereas incorporating depth-resolved borehole temperature measurements for initialization yields substantially improved agreement with observations. The findings highlight the sensitivity of short-term predictive modeling to the representation of initial subsurface thermal conditions and underscore the value of high-resolution field measurements for model calibration and validation. These preliminary results inform ongoing efforts to extend the modeling framework to longer time horizons and to refine monitoring and modeling strategies that support the design guidance and operational optimization of GHP systems in research and commercial buildings.
Lignin has been hailed as an ideal renewable alternative for petrochemical-based prepolymers in material synthesis for a sustainable and circular economy, due to its abundant aromatic network and high carbon content. However, the properties and performance of lignin-derived macromolecules are strongly influenced by the lignin itself. While numerous studies have explored the impact of lignin content on the thermomechanical performance of lignin-based vitrimers, literature on how the inherent structural features of lignin affect these properties is scanty. In this study, hardwood organosolv lignin was fractionated in ethyl acetate, ethanol, and acetone to obtain lignin fractions with varying structural characteristics. These fractions were then modified through carboxylation and crosslinked with epoxidized soybean oil (ESO) at a hydroxyl to epoxy group ratio of 1:1 to produce lignin-based transesterification vitrimers (LVs). The thermal properties (i.e. glass transition temperature and thermal stability), tensile strength, storage modulus, and stress relaxation behavior of the LVs were studied and carefully related to the structural features of lignin. The results revealed a positive relationship between strong hydroxyl content in modified lignin and the tensile strength (5.10–9.71 MPa), storage modulus (1099.4 – 1372.8 MPa), crosslinking density, and stress relaxation of the LVs. Additionally, both the storage modulus and tensile strength exhibited a positive relationship with the ratio of rigid linkages in modified lignin, while lignin molecular weight was found to significantly impact the thermal properties of LVs (i.e Tg and thermal stability). This study not only highlights the valorization of lignin in vitrimer synthesis but also provide insights for designing lignin-based materials with tailored properties for specific applications.
A knowledge gap exists in the data and understanding of fresh fuel salt and irradiated multicomponent fuel salt systems thermophysical properties. Quantifying these properties is necessary for the design and construction of test reactors, as well as the licensing of future commercial molten-salt reactors. To facilitate thermal property determination on a proposed fuel salt composition for Seaborg Technologies, several samples containing depleted uranium tetrafluoride (UF4), sodium fluoride (NaF), and potassium fluoride (KF) were blended, and a melt temperature analysis was performed. From the melting temperature analysis, it was determined that sample Seaborg-7, a ternary salt composition of 26.4UF4-24.7KF-48.9NaF (mol%), was very near a ternary eutectic point. Therefore, thermal properties such as melting temperature, salt stability, density, heat capacity, thermal diffusivity, and viscosity were experimentally determined on the Seaborg-7 salt. These measurements document the baseline properties of fresh fuel salt as a function of temperature, where future experiments on irradiated fuel salt will provide a holistic perspective on the change of thermophysical properties during reactor operations. Several precision instruments were used to collect property data, and instrument calibrations and data collection were performed and documented in a standardized and reproducible manner with meticulous detail. This process ensured that the measurement procedures and resulting data can readily be duplicated elsewhere. The Seaborg-7 salt was shown to be stable at temperatures up to 900°C, as no mass change was observed upon repeated heating and cooling. The peak melting temperature was determined to be 547°C (557°C endset). The enthalpy of fusion (??H?_fus^o) was determined to be 167.5 ± 2.7 J/g while the enthalpy of crystallization (??H?_c^o) was determined to be -147.8 ± 13.3 J/g. In addition to the eutectic melting peak, upon heating, several pre eutectic peaks were observed, occurring at 470°C (onset) and 499°C (peak). Specific heat capacity measurements showed a slightly increasing trend with respect to temperature in the solid phase, while the liquid-specific heat capacity showed a somewhat flat trend with an average value of 106.1 ± 1.24 J/mol·K between 600 to 800°C. Three independent trials using the Seaborg-7 salt determined the density to be ?(T) = 4.908 – 0.000363·T(°C), validated between 32 to 200°C, and ?(T) = 4.808 – 0.00113·T(°C), validated between ~575 to 850°C. Thermal diffusivity was determined for the liquid state and is represented by the linear equation y = 0.1581 + 0.000207·T(°C) between 550 to 850°C. The viscosity was determined from 600 to 800°C and is represented by the exponential fit equation, ? (mPa·s) = 736.58e^(-0.006·T(°C)). This report documents the conclusion of fuel salt thermophysical property measurements for the Seaborg SPP, Phase A project.
We measured and modeled near infrared extinction of nuclear waste glasses from 300 °C to 1150 °C to enable predictive radiation heat transfer and thermal conductivity estimates. A composition and redox informed model resolved contributions from key chromophores (Fe+2-O-Fe+3, V+4, free and bonded ?Si-OH groups) and, when present, spinel particles that can cause strong scattering. The model reproduced measured absorption from room temperature up to 1150 °C, with minor discrepancies near 1 µm (likely due to possible trace impurities) and 2.5 µm (linked to uncertainty in hydroxy groups). Spectra showed silicate melts were semitransparent mainly in the 0.5–4.0 µm window, responsible for radiation thermal conductivity that generally increases with increasing temperature. We quantified the dependence of effective thermal conductivity on dissolved water and provided distributions across >100 LAW/HLW/DFHLW melts at 1150 °C, supporting improved melter heat transfer modeling.
Time-domain thermoreflectance (TDTR) has been instrumental in measuring the heat transfer properties of bulk and nanostructured materials over the past two decades. In this Primer, we describe the optical and thermal aspects of TDTR, with an in-depth discussion on the theory, apparatus design and implementation. We present examples that illustrate the ability of TDTR to measure thermal conductivity tensors, thermal conductance across material interfaces, and volumetric heat capacity of thin films, 2D materials and bulk materials. The ability of TDTR to spatially resolve thermal properties is useful for studying heterogeneous material systems, such as materials processed in or subjected to extreme environments. We consider current limitations of pump–probe metrologies and discuss recent advancements of TDTR, such as time-resolved magneto-optic Kerr effect (TR-MOKE), beam-offset TDTR/TR-MOKE, steady-state thermoreflectance, frequency-domain thermoreflectance and laser-flash TDTR. Lastly, we present an outlook on anticipated technological developments to further expand the ability of TDTR to measure nanoscale thermal properties.
Here, a systematic study of atomic dynamics and thermal properties of the family of layered chalcogenide compounds 𝑀Cr𝑋 2 (𝑀= Ag, Cu; 𝑋 = S, Se) was performed, including neutron and x-ray scattering, thermal characterization, and first-principles simulations. In all compounds, we observe a breakdown of specific phonon modes across the superionic phase transition, for phonons whose eigenvectors exhibit large contributions of mobile ions. In particular, the nondispersive portions of transverse acoustic (TA) branches at short-wavelengths and the low-energy optical phonons with large contributions from Ag + or Cu + become severely damped in the superionic phase. However, well-defined quasiparticles persist in the superionic state for long-wavelength TA phonons. In the case of AgCrS 2 , the coupling of lattice dynamics with its antiferromagnetic transition was also investigated. The magnetic ordering couples with the monoclinic–rhombohedral structural transition, and the Cr 3+ spin arrangement strongly affects the phonon dispersions. We qualitatively reproduce the magnetic and nuclear components of the INS measurement for antiferromagnetic AgCrS 2 by combining models of spin-waves and spin-polarized first-principles phonon simulations. Quasielastic magnetic fluctuations persist in the paramagnetic phase up to high temperature, but are clearly distinguished from the nuclear component through their momentum dependence. Finally, we report measurements of the thermal properties of the selenide compounds and find good agreement with our DFT simulations.
Abstract Bulk superconductors can be used as super-strength quasi-permanent magnets capable of providing magnetic flux densities considerably superior to conventional permanent magnets. This makes them attractive for several engineering applications that rely on strong magnetic fields like rotating machines, NMR/MRI and magnetic drug delivery systems. Recently, the authors reported a record trapped magnetic field in an iron-based bulk superconductor: 2.83 T was trapped in potassium-doped barium iron arsenide (Ba, K)Fe 2 As 2 (or Ba122) at 5 K. Of particular significance is that the strength and temporal stability of this magnetic field exceeds the requirements of MRI machines, indicating iron-based bulks can now perform at levels demanded by engineering applications. One crucial challenge for their practical use, however, is the need to apply and remove an external magnetic field to magnetize them. Pulsed field magnetization (PFM) shows great promise as a practical method of magnetizing bulks, but the process generates heat in the bulk that is detrimental to its superconducting performance and ability to act as a super-strength magnet. In this paper, coupled electromagnetic–thermal numerical models are used to simulate the PFM of iron-based bulk superconductors. Here we focus on the recent-record-breaking, fine-grain polycrystalline K-doped Ba122 bulks. The impact that the specific J c ( B ) characteristics and thermal properties of the Ba122 material—all of which have been experimentally measured from state-of-the-art samples—have on the magnetic flux dynamics and thermal behaviour during PFM, including the final trapped field, is investigated. We show that because the thermal properties are similar to those of REBa 2 Cu 3 O 7 −δ bulks, a similar response to pulsed fields is obtained. A maximum trapped field of ∼0.81 T (∼43.3% of the maximum trapped field capability under ideal, field-cooling conditions) was simulated at 5 K, with a magnetization efficiency of ∼54%. The modelling framework provides a fast and flexible tool for optimising the practical PFM process at different operating temperatures to maximise the trapped field in state-of-the-art Ba122 bulks and to guide the design of future experiments.
The reactivity-initiated accident (RIA) is a complex scenario with several tightly interacting physical phenomena. Accurately predicting fuel behavior during these transients is difficult due to limitations in the modeling of fuel-to-coolant heat transfer. Common approaches to simulate RIAs involve standalone calculations using either a fuel performance code or a thermal-hydraulic code. The complex interdependencies of thermal-hydraulic and fuel mechanical behavior suggest that a tight coupling between these codes may provide more accurate predictions of fuel-to-coolant heat transfer and cladding mechanical response. Here, RELAP5-3D and BISON are coupled in this paper to simulate RIAs, and a sensitivity analysis is performed to rank key thermal properties and two-phase heat transfer parameters relevant for fuel-to-coolant heat transfer and cladding failure mechanisms in UO 2 –Zircaloy-4 systems. Gas gap conductance, film boiling heat transfer uncertainty, pulse width, fuel-specific heat capacity, and cladding-specific heat capacity were identified as important parameters. Variations in figures of merit resulting from changes to pulse width and the material thermal properties indicate that time-dependent heat transfer rates are significant for safety-relevant mechanical parameters due to the time dependence of cladding ductility and pellet-cladding mechanical interaction loading. The results suggest that the thermal-hydraulic factors have a nonnegligible influence on the thermomechanical solution and vice versa. Tight coupling of both sets of physics is recommended to improve prediction of fuel behavior during RIAs. Highlights include the following: 1. The RELAP5-3D thermal-hydraulic code and the BISON fuel performance code are tightly coupled for simulation of RIA transients with energy depositions at the Zircaloy-4 cladding failure threshold. 2. Departure from nucleate boiling occurred for all simulated cases. Due to the ductility of fresh fuel, substantial ballooning occurred in most cases. 3. Gas gap conductance, fuel-specific heat capacity, cladding-specific heat capacity, transient pulse width, and film boiling heat transfer were the dominant thermal factors impacting the safety figures of merit at energy depositions.
Poly(3-hydroxyhexanoate) (P3HHx)-based poly(3- hydroxyalkanoate)s (PHAs) are biologically produced, commercially implemented PHAs, but little is known in the open literature about the structure and property characterizations of discrete, authentic homopolymer P3HHx, and its copolymers with poly(3- hydroxybutyrate) (P3HB). Here, we introduce a chemocatalytic route to stereoregular (both isotactic and syndiotactic) P3HHx and PHA copolymers with P3HB, including statistical copolymer P3HBHx with various levels of 3HHx incorporation and discrete hard−soft−hard triblock copolymer P3HB-b-P3HHx-b-P3HB, and present extensive characterizations of their structures, thermal properties, and mechanical performance. The synthesis is efficiently achieved by one-pot polymerization of eight-membered di(n-propyl) and dimethyl (for copolymerization) diolides catalyzed by chiral molecular catalysts. Notably, P3HBHx can be rapidly produced in quantitative yield and exhibits high molar mass (M n up to 551 kg/mol) as well as both high modulus (up to 1.39 GPa) and ductility (up to 445%), while P3HB-b-P3HHx-b-P3HB further extends application temperature windows by possessing a unique combination of low T g (−18 °C) and high T m (154 °C) values. These findings highlight the stereomicrostructural and architectural versatility of chemocatalytic routes to PHAs, which, in turn, can be utilized to largely tune the PHA thermal properties and mechanical performance.
The overall objective of this project was to produce a coal-based construction material that has up to ~95 weight percent (wt. %) coal with physical, chemical, and thermal properties exceeding those of ordinary Portland cement (OPC)-based construction materials. Additionally, the project aimed to minimize external binders by implementing novel mixing techniques, while exceeding the performance/cost ratio of OPC. Finally, the project was to demonstrate production of precast products via the design and fabrication of products via a bench scale process. Consistent with some of these objectives, the project successfully fabricated samples of coal-based composite materials with >80 wt% coal with physical, chemical and thermal properties on par with cement-based concrete. Select samples demonstrated compressive strengths with >7,000 psi and flexural strength of >420 psi. The composite materials minimized external binders and also demonstrated durability, as evidenced by resistance to acidic and basic solutions. Finally, larger slab and beam type samples were produced using a process developed by the Recipient, although, the process was not semicontinuous in nature. Taken together, the results of this project suggest that domestic coal has potential to serve as a replacement for cementitious materials utilized in incumbent construction technologies, which could significantly reduce the energy and emissions of the construction industry
Carbon materials have become increasingly diverse, finding applications in high-temperature and high-radiation environments. Glassy carbon, an allotrope known for its exceptional chemical inertness and desirable mechanical properties. However, understanding neutron irradiation effects in glassy carbon has proven challenging, primarily because of its unique nanopore structure. Here, this study presents a highly detailed microstructural characterization investigation of neutron-induced changes in glassy carbon, revealing how changes in nanopore structure and crystallinity impact the irradiation-induced shrinkage. Aberration-corrected scanning transmission electron microscopy (STEM) reveals pore closure that leads to material densification in the irradiated samples. Dimensional analysis combined with comparison to historical data suggest significant length shrinkage to occur. Neutron and in situ electron irradiation experiments suggest that glassy carbon transforms into so-called carbon onions, supporting the concept of shrinkage saturation. Investigating irradiation temperature effects using STEM, electron energy loss spectroscopy, x-ray diffraction, and Raman spectroscopy revealed partial amorphization at 210 °C–230 °C and preserved order in glassy carbon at 860 °C, coinciding with pore closure. Thermal property measurements were also conducted to assess the effects of densification and other changes in the atomic structure of glassy carbon. The results of this study have broad implications in the deployment of glassy carbon to nuclear environments, based around the observed changes in the thermal properties, and demonstrates the operational window for the onset of densification.