Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal dynamics”

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

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

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Theoretical H + O 3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties

Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O 3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. Additionally, the theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O 3 → OD + O 2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Study of the Segmental Dynamics and Ion Transport of Solid Polymer Electrolytes in the Semi-crystalline State

Solid polymer electrolytes are promising in fulfilling the requirements for a stable lithium metal anode toward higher energy and power densities. In this work, we investigate the segmental dynamics, ionic conductivity, and crystallinity of a polymer electrolyte consisting of poly(ethylene oxide) (PEO) and lithium triflate salt, in the semi-crystalline state. Using quasi-elastic neutron scattering, the segmental dynamics of PEO chains confined between the crystalline lamellae is quantified, using Cole–Cole analysis. We show that the structural relaxation time, τ 0 , of PEO equilibrated near room temperature is six-fold longer than the same sample that had just cooled down to room temperature. This corresponds to a three-fold smaller ionic conductivity in the equilibrated condition. This work reveals that the segmental dynamics of semi-crystalline polymer electrolytes is very sensitive to thermal history. We demonstrate that quasi-elastic neutron scattering can be used to characterize the ion transport and segmental dynamics in the semi-crystalline state.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Theoretical study of the electronic, thermodynamic, and thermo-conductive properties of y-LiAlO2 with 6Li isotope substitutions for tritium production

Lithium aluminate has attracted researchers’ interests due to its wide applications. The structural, electronic, optical, lattice phonon thermodynamic, and thermo-conductive properties of ?-LiAlO2 with 6Li-isotope substitution are investigated by density functional theory and lattice phonon dynamics. The calculated results show that ?-LiAlO2 possesses a wide band-gap of 4.63 eV. The valence band is mainly formed by p orbitals of Li, O, and Al. The calculated lattice thermal conductivity of ?-LiAlO2 is found to be in good agreement over a wide range of temperatures with the available experimental data. The calculated dielectric matrix shows anisotropic behavior along z and x (or y) axes. Due to the wide band-gap, the optical conductivity of ?-LiAlO2 is mainly contributed from the electron hopping between valence-conduction bands. Substitution of 7Li with 6Li isotope in ?-LiAlO2 leads to observable differences in the lattice phonon frequencies at higher frequencies (>8THz) and to slight changes in corresponding thermal conductivity and the infrared and Raman spectra which are in good agreement with the measured data. The relevance of these properties for tritium production through 6Li absorbing neutron in nuclear reactors is discussed.

ab initio thermodynamics, ?-LiAlO2 and 6Li-isotope↗
DOE OSTI2021

Thermal management of a high temperature sodium sulphur battery stack

The sodium sulfur battery is an advanced secondary battery with high potential for grid-level storage due to their high energy density, low cost of the reactants, and high open-circuit voltage. However, as the operating temperature of the battery is high (about 300 °C), effective thermal management is required to prevent thermal runaway under high current density operation. Here, to develop efficient thermal management strategies, a detailed, thermo-electrochemical, non-isothermal, distributed dynamic model of a sodium-sulfur battery stack is developed. Models of three thermal management strategies are developed and analyzed in this work: active cooling, passive cooling, and hybrid cooling. The active cooling strategy uses air as the cooling medium whereas, the passive cooling strategy uses a phase change material. In the hybrid cooling strategy, both active and passive cooling strategies are used. Due to the high operating temperature of the sodium sulfur batteries, the rejected heat can be utilized effectively and consequently a high variability in the heat rejection rate may not be acceptable. Performance of these strategies is analyzed under high current density operation by evaluating the variability in the cell temperature during charge/discharge cycles, temperature difference between the cells based on their locations, and variability in the heat rejection rate. The phase change-based hybrid thermal management strategy with an embedded controller maintained the temperature variation within ±/- 2.8 °C. Considering the tradeoff between the capital cost for PCM, fan power requirement, and variability in the heat rejection rate, the optimal quantity of the PCM and the maximum air velocity for the maximum current density of 260 mA/cm 2 (0.55 C-rate) were found be 0.74 g/Wh and 3 m/s, respectively.

25 ENERGY STORAGE↗
DOE OSTI2021

Signatures of Coherent Phonon Transport in Ultralow Thermal Conductivity Two-Dimensional Ruddlesden–Popper Phase Perovskites

An emerging class of methylammonium lead iodide (MAPbI 3 )-based Ruddlesden-Popper (RP) phase perovskites, BA 2 MA n-1 Pb n I 3n+1 (n = 1-7), exhibit enhanced stability to environmental conditions relative to MAPbI 3 , yet still degrade at elevated temperatures. Here, we experimentally determine the thermal conductivities of these layered RP phases for n = 1-6, where n defines the number of repeated perovskite octahedra per layer. We measure thermal conductivities of 0.37 ± 0.13/0.12, 0.17 ± 0.08/0.07, 0.21 ± 0.05/0.04, and 0.19 ± 0.04/0.03 W/m∙K in thin films of n = 1-4 and 0.08 ± 0.06/0.04, 0.06 ± 0.04/0.03, 0.06 ± 0.03/0.03, and 0.08 ± 0.07/0.04 W/m∙K in single crystals of n = 3-6. With the exception of n = 1, these thermal conductivities are lower than the range of 0.34-0.50 W/m∙K reported for single-crystal MAPbI 3 . Reduced-order lattice dynamics modeling suggests that the initially decreasing trend of thermal conductivity in similarly oriented perovskites with increasing n may result from the transport properties of coherent phonons, emergent from the superstructure, that do not scatter at the interfaces of organic butylammonium chains and perovskite octahedra. Reduced group velocity of coherent phonons in n = 3-6, a consequence of band flattening in the phonon dispersion, is primarily responsible for their ultralow thermal conductivities. Similar effects on thermal conductivity have been experimentally demonstrated in deposited superlattices, but never in naturally defined materials such as RP phases. GIWAXS measurements reveal that higher n RP phase thin films are less orientationally controlled and therefore possess apparently elevated thermal conductivities relative to single crystals of the same n.

2D perovskites↗
DOE OSTI2020

Visualizing Energy Transfer at Buried Interfaces in Layered Materials Using Picosecond X-Rays

Understanding the fundamentals of nanoscale heat propagation is crucial for next-generation electronics. For instance, weak van der Waals bonds of layered materials are known to limit their thermal boundary conductance (TBC), presenting a heat dissipation bottleneck. In this paper we present a new non-destructive method to probe heat transport in nanoscale crystalline materials using time-resolved x-ray measurements of photo-induced thermal strain. This technique directly monitors time-dependent temperature changes in the crystal and the subsequent relaxation across buried interfaces by measuring changes in the c-axis lattice spacing after optical excitation. We investigate films of five different layered transition metal dichalcogenides MoX 2 [X= S, Se, Te] and WX 2 [X= S, Se] as well as graphite and a W-doped alloy of MoTe 2 , on c-plane sapphire substrates, finding TBC values in the range 10-30 MWm -2 K -1 at room temperature. In conjunction with molecular dynamics simulations, we show that the high thermal resistances are a consequence of weak interfacial van der Waals bonding and low phonon irradiance. This work paves the way for an improved understanding of thermal bottlenecks in emerging three-dimensional heterogeneously integrated technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Heat pump water heater enhanced with phase change materials thermal energy storage: Modeling study

A promising solution to improve the first hour rating (FHR) of a heat pump water heater (HPWH) involves employing a secondary tank which contains phase change material (PCM) capsules. To better understand the influence of PCM thermal storage on the HPWH operational performance, a dynamic model was developed to simulate and analyze the behavior of a newly developed HPWH technology that incorporates PCM storage into a standard HPWH to optimize key parameters such as the uniform energy factor and FHR. Mathematical models of several key components of the proposed HPWH-PCM integrated thermal energy storage (TES) system, e.g., water heater tank, PCM TES tank, evaporator, compressor, and expansion valve, have been elaborated. Also, a model-based control co-simulation platform was developed to integrate a embed PCM storage HPWH dynamic model with a control model for better supporting control design, analysis, verification, and validation. The model accuracy has been validated through comparing simulation results with lab test results, with a mean average percentage error of <5.5% for most of the selected performance variables. In addition, using the developed co-simulation platform, the demand response control strategy was studied to evaluate the load flexibility of the combined HPWH-PCM storage system by shifting system power usage to outside of the 3.5 h peak load period.

25 ENERGY STORAGE↗
DOE OSTI2022

Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal

Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Impact of fluorination on Li + solvation and dynamics in ionic liquid-hydrofluoroether locally concentrated electrolytes

The thermal, physical, structural, and transport properties of ionic liquid (IL) electrolytes based on n-methyl-n-butylpyrrolidinium bis(trifluoromethanesulfonyl)imide [PYR14][TFSI] and Li-salts of lithium (nonafluorobutane)(trifluoromethanesulfonyl)imide [Li][IM14] and [Li][TFSI] (0 ≤ x Li ≤ 0.3), with addition of 1,1,2,2-tetrafluoroethyl 2,2,3,3-tetrafluoropropyl ether (TTE) (0 ≤ x TTE ≤ 0.45) were studied. While [IM14] increases glass transition, TTE improves conductivity over a wider liquidus range. Li + is solvated primarily by [TFSI] with some contribution by [IM14]. However, both 7 Li- 19 F HOESY NMR and spatial distribution functions (SDFs) derived from molecular dynamics (MD) indicate short contacts between the fluorines of TTE and Li + . This is a result of tighter anion-solvated Li + without aggregation of the Li + solvates, consistent with Raman measurements, thus confirming the existence of fluorous domains in the bulk which leads to improved fluidity and a Li + diffusivity of 1.26 x 10 –12 m 2 /s at 0 °C (x Li = 0.20) with a Li + transference of 0.16. Improved transport properties translated to a higher capacity of the TTE/IL electrolyte in a Li||LiFePO 4 half-cell with 95mAh/g at 0.1C, compared to the IL electrolyte without TTE (60 mAh/g). Furthermore, this study demonstrates the tunability of solvation and transport properties in IL electrolytes by asymmetric fluorinated anions and hydrofluoroether co-solvents for low temperature Li-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Contribution of the Hydroxyl Group on Interfacial Heat Conduction of Monohydric Alcohols: A Molecular Dynamics Study

Monohydric alcohols have been used as promising phase change materials (PCMs) for low-temperature latent heat storage. However, the heat storage/retrieval rates are limited due to the low thermal conductivity of such alcohols. In this work, nonequilibrium molecular dynamics (NEMD) simulations were performed to study the microscopic heat conduction in example monohydric alcohols, i.e., 1-dodecanol (C12H26O), 1-tetradecanol (C14H30O), and 1-hexadecanol (C16H34O). A simplified ideal crystal model was proposed to exploit the potential for improving the thermal conductivity of monohydric alcohols. The effect of ideal crystalline structures, especially the contribution of the hydroxyl group, on the microscopic heat conduction process was analyzed. The thermal conductivity of the ideal crystals of the various monohydric alcohols was predicted to be more than twice as compared to that of their respective solids. The major thermal resistance in the ideal crystals was found around the molecular interfaces, as a result of the excellent heat conduction performance along the linear molecular chains. The calculated vibrational density of states (VDOS) and interfacial heat transfer were then investigated. When the interfaces are surrounded by hydroxyl groups as walls, strong hydrogen bond (HB) interactions were observed. The interfacial heat transfer coefficient of the ideal crystalline structures of 1-tetradecanol was found to reach up to ~735.6 MW/m2 W. It was elucidated that the high interfacial heat transfer rate is clearly related to the stronger intermolecular interactions.

Engineering↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2023

Techno-Economic Analysis of Dynamic Building Envelopes Comprising Phase Change Materials and Switchable Insulations: Preprint

Integrating phase change materials (PCMs) in building envelopes is a recognized technique to reduce the space heating/cooling loads and provide load shedding and shifting capacity. However, PCM benefits have been found to reduce dramatically with increase in the insulation requirements. Dynamic building envelopes that can alter their thermal resistance based on the indoor and outdoor conditions can effectively harness free ambient heating/cooling, thereby greatly enhancing the benefits of the PCM-integrated envelope in managing thermal loads. In this study, we examine various combinations of PCMs and switchable insulations and compare their combined impact on heat flow through the building envelope. Using numerical analysis, we investigate the PCM- and switchable insulation-integrated building walls, calculate the potential energy saving benefits compared to the conventional static walls, and finally perform technoeconomic analysis to estimate the acceptable cost of the technology under various payback period scenarios.

dynamic building envelopes↗
DOE OSTI2022

Techno-Economic Analysis of Dynamic Building Envelopes Comprising Phase Change Materials and Switchable Insulations

Integrating phase change materials (PCMs) in building envelopes is a recognized technique to reduce the space heating/cooling loads and provide load shedding and shifting capacity. However, PCM benefits have been found to reduce dramatically with increase in the insulation requirements. Dynamic building envelopes that can alter their thermal resistance based on the indoor and outdoor conditions can effectively harness free ambient heating/cooling, thereby greatly enhancing the benefits of the PCM-integrated envelope in managing thermal loads. In this study, we examine various combinations of PCMs and switchable insulations and compare their combined impact on heat flow through the building envelope. Using numerical analysis, we investigate the PCM- and switchable insulation-integrated building walls, calculate the potential energy saving benefits compared to the conventional static walls, and finally perform techno-economic analysis to estimate the acceptable cost of the technology under various payback period scenarios.

dynamic building envelopes↗
DOE OSTI2024

A computational analysis of effective R-values of buried ducts – the dynamic performance of buried ducts

Here, this paper evaluates the thermal performance of ducts partially or fully buried in loose-fill attic insulation. The overall thermal resistance between the ducts and the attic is referred to as an effective R-value. This paper shows a strong dependency of assumed attic temperature on the effective R-value. Based on the results, the effective R-value can be about twice as much with an attic temperature of 130 °F [54.4 °C], compared to when the attic temperature is 80 °F [26.7 °C]. Thus, this paper provides a polynomial regression equation based on a large set of simulations to determine the effective R-value of buried ducts depending on attic temperature and whether the HVAC system runs in cooling or heating mode. Further, the work presented in this paper investigated the potential impact of convective airflow within the attic insulation, particularly around the ducts. The analysis was based on computational fluid dynamics (CFD) and indicated that convectional forces are presented around the exterior surface of the ducts, but with negligible impact on the overall heat balance between the duct and the attic space.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Materials genome innovation for computational software (magics) center

Functional layered material (LM) architectures will dominate nanomaterials science in this century. We have developed theory, modeling, simulation, and software and data tools that enhance understanding and AI guide synthesis, enable characterization of complex structures, and improve capabilities in the predictive design and growth of LMs. Research at the Center has focused on: Computational synthesis and characterization: AI guided synthesis and experimental synthesis of stacked LMs with tailored properties via optimized chemical vapor deposition (CVD) growth and liquid-phase exfoliation; study defects, edges, grain boundaries, wrinkling of atomic layers and their effects on chemical, mechanical, electrical, and optical properties. Far-from-equilibrium processes: Joint experimental and simulation based probe of electronic processes with NAQMD and ultrafast X-ray free-electron laser (XFEL) and ultrafast electron diffraction (UED) facilities at Stanford. Experimentally validate NAQMD by ultrafast electron diffraction and X-ray spectroscopy studies of structural and excited state dynamics, shape fluctuations, and phonon dynamics. Scalable software: Simulation engines for desktop-to-exascale platforms using low-overhead, linear-scaling QMD algorithms; divide-conquer-recombine NAQMD with electronic excitations; extended-Lagrangian reactive molecular dynamics (RMD), machine learning (ML) based neural-network quantum molecular dynamics (NNQMD), and super-state accelerated molecular dynamics (AMD) and kinetic Monte Carlo codes; thermal and electrical transport software; and design 3D architectures of LMs with desired functionality using scalable software. Distribution of software and data, and training: Software and simulation-experimental data generated within the Center are distributed to the materials science community via Berkeley Materials Project (MP) framework. We have also organized three workshops for software distribution and training at USC (Nov. 2017, Mar. 2018) and Gaithersburg, MD (Nov. 2018) to train researchers, with the last one in focused on underrepresented groups, in collaboration with Howard University which is one of the largest HBCUs. The Center supported a total of 46 personnel and 6 undergraduate students. These include 14 faculty, 11 postdoctoral research associates, 20 graduate research assistants, and mentored 6 undergraduate students. This resulted in the publications of 63 research papers that include 46 publications on Reactive and Quantum Dynamics Simulations, 13 publications on Machine Learning for Quantum Materials, and 4 publications on Quantum Computing.

2D Materials↗