Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 1,639 records · Page 91

Coupled thermal/chemical/mechanical modeling of energetic materials in ALE3D

We must improve our ability to model the response of energetic materials to thermal stimuli and the processes involved in the energetic response. We have developed and used a time step option to efficiently and accurately compute the hours that the energetic material can take to react. Since on these longer film scales, materials can be expected to have significant motion, it is even more important to provide high-order advection for all components, including the chemical species. We show an example cook-off problem to illustrate these capabilities.

WEAPONS↗

Reversible Ligand Protonation in Noninnocent Constrained-Geometry-Like Group 4 Complexes

The proligands C 5 H 5 CMe 2 CHPhOx R (R = Me 2 , CHMe 2 , and CMe 3 ) react with M(NMe 2 ) 4 (M = Ti, Zr, and Hf) to give monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 , doubly deprotonated (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 , or a mixture of both. The observed products depend on reaction conditions, the oxazoline substituent, and the metal center, with 4-t-butyl-oxazoline or Ti giving only the monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 . Amine elimination from (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 to give (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 is reversible for 4,4-dimethyl- and 4-isopropyl-oxazoline-based ligands in Zr or Hf complexes. Temperature-dependent kinetic studies of the equilibration of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 , HNMe 2 and (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 provide the experimental thermodynamic parameters ΔS° = -17.4 ± 2.6 cal·mol –1 K –1 and ΔH° = -6.8 ± 0.8 kcal·mol –1 . An Eyring plot of the rate constants, determined from the system as it approaches equilibrium, gives the activation entropy and activation enthalpy for the addition of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 and HNMe 2 as -36 ± 4 cal·mol –1 K –1 and 8.9 ± 1.3 kcal·mol –1 , respectively; the elimination of HNMe 2 from (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 is characterized by ΔS ‡ = -19 ± 5 cal·mol –1 K –1 and ΔH ‡ = 15.7 ± 1.5 kcal·mol –1 . DFT computational models indicate a single-step, nonlinear transfer of the H between the benzylic position of the noninnocent, oxazoline-coordinated ligand and NMe 2 . Computations also confirm the negative activation entropy and the trends in the barriers support the experimental results. Together, these studies indicate the importance of steric effects from the oxazoline ligand, metal center, ancillary ligands, and leaving group on the shuttling of the proton between HNMe 2 and the noninnocent ligand. Finally, these effects suggest that coordination of oxazoline to the metal center is a key part of the benzylic deprotonation and noninnocent behavior of the cyclopentadienyl-oxazoline ligand.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discovery of a Wurtzite-like Cu 2 FeSnSe 4 Semiconductor Nanocrystal Polymorph and Implications for Related CuFeSe 2 Materials

I 2 -II-IV-VI 4 and I-III-VI 2 semiconductor nanocrystals have found applications in photovoltaics and other optoelectronic technologies because of their low toxicity and efficient light absorption into the near-infrared. Herein, we report the discovery of a metastable wurtzite-like polymorph of Cu 2 FeSnSe 4 , a member of the I 2 -II-IV-VI 4 family of semi-conductors containing only Earth-abundant metals. Density functional theory calculations on this metastable poly-morph of Cu 2 FeSnSe 4 indicate that it may be a superior semiconductor for solar energy and optoelectronics applications compared to the thermodynamically preferred stannite polymorph, since the former displays a sharper dispersion of energy levels near the conduction band minimum that can enhance electron mobility and suppress hot electron cooling. The experimental optical band gap was measured by the inverse logarithmic derivative method to be direct, in agreement with theory, and in the range of 1.48-1.59 eV. Mechanistic studies reveal that this metastable phase de-rives from intermediate Cu 3 Se 2 nanocrystals that serve as a structural template for the final hexagonal wurtzite-like product. We compare the chemistry of wurtzite-like Cu 2 FeSnSe 4 to the related CuFeSe 2 material system. Our experimental and computational comparisons between Cu 2 FeSnSe 4 and CuFeSe 2 help explain both the crystal chemistry of CuFeSe 2 that prevents it from forming wurtzite-like polymorphs and the essential role of Sn in stabilizing the metastable structure of Cu 2 FeSnSe 4 . Furthermore, this work provides insight into the importance of elemental composition when designing syntheses for metastable materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Alloying Effects on the Transport Properties of Refractory High-entropy Alloys

Additive Manufacturing (AM) has opened new frontiers for the design of refractory high-entropy alloys (HEAs) for high-temperature applications. The thermal conductivity of the AM feedstock is among the most important thermo-physical properties that control the melting and solidification process. Despite its significance, there remains a notable gap in both computational and experimental research concerning the thermal conductivity of HEAs. Here, we use density functional theory (DFT) to systematically investigate the alloying effects on the transport properties of Ti-Cr-Mo-W-V-Nb-Ta RHEAs, including electrical and thermal conductivities and Seebeck coefficient. The relaxation time of charge carriers is a key underlying parameter determining thermal conductivity that is exceedingly challenging to predict from first principles alone, and we thus follow the approach by Mukherjee, Satsangi, and Singh [Chem Mater 32, 6507 (2022)] to optimize the relaxation time for RHEAs. Here we validated thermal conductivity predictions on elemental solids, binary and ternary alloys, and RHEAs and compared them against thermodynamic (CALPHAD) predictions and our experiments with good correlations. To understand observed trends in thermal conductivity, we assessed the phase stability, electronic structure, phonon, and intrinsic- and tensile strength of down-selected RHEAs. Our electronic structure and phonon results connect well with the observed compositional trends for thermal transport in RHEAs. Our DFT assessment and CALPHAD predictions provide a unique design guide for RHEAs with tailored thermal conductivity, a critical consideration for AM and thermal-management applications.

36 MATERIALS SCIENCE↗

Enhanced Twist-Averaging Technique for Magnetic Metals: Applications Using Quantum Monte Carlo

In this study, we propose an improved twist-averaging (TA) scheme for quantum Monte Carlo methods that use converged Kohn–Sham or Hartree–Fock orbitals as the reference. This TA technique is tailored to sample the Brillouin zone of magnetic metals, although it naturally extends to nonmagnetic (NM) conducting systems. The proposed scheme aims to reproduce the reference magnetization and achieves charge neutrality by construction, thus avoiding the large energy fluctuations and the postprocessing needed to correct the energies. It shows the most robust convergence of total energy and magnetism to the thermodynamic limit (TDL) when compared to four other TA schemes. Diffusion Monte Carlo applications are shown on NM Al and ferromagnetic α-Fe. The cohesive energy of Al in the TDL shows an excellent agreement with the experimental result. Furthermore, the magnetic moments in α-Fe exhibit rapid convergence with an increasing number of twists.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗

Emission Spectroscopic Measurements in the Plenum Region of the NASA IHF Arc Jet Facility

A newly designed segment with optical access was installed in the plenum chamber of the 60 MW Interaction Heating arcjet Facility at NASA Ames Research Center. This special segment has ports located off axis, and the optical fibers can be inserted into these ports. The special segment allows for optical examination of the arc-heated flow as it enters the plenum, and thus assists in determining estimates of the thermodynamic state of the inflow to the convergent section of the nozzle. In the present work, optical emission measurements have been made in VIS-NIR region (wavelengths between 500 nm to 900 nm) for two settings of the arc heater - a 6000 A condition (high condition) with the minimum amount of radial injection of cold air in the plenum, and a 3300 A condition (low condition) with significant amount of cold air injection to reduce the enthalpy of the arc-heated stream. The results presented here were obtained using an Acton SP300i spectrometer coupled to a Princeton Instruments PI-max intensified camera. In addition to the optical emission measurements, computations were performed for the flow in the plenum and radiation along lines of sight corresponding to the optical ports. Along the centerline, i.e., the longest line of sight across the plenum cross-section, there is good agreement between computations and measurements for the high enthalpy condition, although the off-axis radial profiles show some differences. For the low enthalpy condition, there are significant differences between computations and measurements. The current working hypothesis is that the computational model does not capture details of the mixing process in the plenum.

Winter, Michael W.↗

BMINN: Learning chemical potentials and parameters from voltage data for multi-phase battery modeling

Free-energy landscapes and chemical potentials govern the dynamics of phase transitions, transport, and stability in functional materials, yet they remain experimentally inaccessible under realistic operating conditions. Here we introduce a Bayesian model-integrated neural network (BMINN) that embeds physics-based formulations of non-autonomous partial differential-algebraic equations into probabilistic learning. This approach reconstructs hidden thermodynamics directly from macroscopic current-voltage data, providing quantitative access to metastable states, staging transitions, and energy barriers without synchrotron probes. Demonstrated on lithium-graphite electrodes, BMINN recovers full Gibbs free-energy landscapes with fidelity validated against operando X-ray diffraction. The framework generalizes across dynamical regimes, enabling accurate voltage prediction, internal state estimation, and inference of governing parameters. Beyond batteries, BMINN exemplifies a broadly applicable strategy for learning missing physics in multiphase, non-equilibrium systems, offering a new pathway to uncover hidden thermodynamic functions across condensed matter and materials physics.

25 ENERGY STORAGE↗

Factors Governing Oxygen Vacancy Formation in Oxide Perovskites

The control of oxygen vacancy (V O ) formation is critical to advancing multiple metal-oxide-perovskite-based technologies. In this work, we report the construction of a compact linear model for the neutral V O formation energy in ABO 3 perovskites that reproduces, with reasonable fidelity, Hubbard-U-corrected density functional theory calculations based on the state-of-the-art, strongly constrained and appropriately normed exchange-correlation functional. We obtain a mean absolute error of 0.45 eV for perovskites stable at 298 K, an accuracy that holds across a large, electronically diverse set of ABO 3 perovskites. Our model considers perovskites containing alkaline-earth metals (Ca, Sr, and Ba) and lanthanides (La and Ce) on the A-site and 3d transition metals (Ti, V, Cr, Mn, Fe, Co, and Ni) on the B-site in six different crystal systems (cubic, tetragonal, orthorhombic, hexagonal, rhombohedral, and monoclinic) common to perovskites. Physically intuitive metrics easily extracted from existing experimental thermochemical data or via inexpensive quantum mechanical calculations, including crystal bond dissociation energies and (solid phase) reduction potentials, are key components of the model. Beyond validation of the model against known experimental trends in materials used in solid oxide fuel cells, the model yields new candidate perovskites not contained in our training data set, such as (Bi,Y)(Fe,Co)O 3 , which we predict may have favorable thermochemical water-splitting properties. The confluence of sufficient accuracy, efficiency, and interpretability afforded by our model not only facilitates high-throughput computational screening for any application that requires the precise control of V O concentrations but also provides a clear picture of the dominant physics governing V O formation in metal-oxide perovskites.

08 HYDROGEN↗

The extended Vinet analytical equations of state

When testing new materials models, it is often not possible to use tabulated Equations of States (EOS) since the model might only be available in a smaller specialized research code without capability to read in tables. In such cases a versatile and easy to implement analytical EOS is needed. In this report I derive, give, and discuss all formulas necessary for implementing the thermodynamically complete extended Vinet EOS 3 into codes. While similar to a Mie-Grüneisen (MG) EOS, this Vinet EOS is based on an isotherm and substantially easier to correctly implement than a MG EOS. In the extended Vinet the isotherm’s density dependence is given as a power series in density to accommodate for higher density behavior which normally requires the Hugoniot based MG.

36 MATERIALS SCIENCE↗

Ab initio Molecular Dynamics Assessment of Thermodynamic and Transport Properties in (K,Li)Cl and (K, Na)Cl Molten Salt Mixtures

Molten salt mixtures are integral part of highly important technological applications such as nuclear reactors. However, due to inherent difficulties associated with experiment at high temperatures and the intrinsic complexity of liquid-phase multi-component systems, understanding their properties at a molecular level remains a challenge. Here, we report on an ab initio molecular dynamics investigation on structural, electronic, transport, and thermal properties of two common molten salt mixtures, (K, Li)Cl and (K,Na)Cl, at five different compositions and three temperatures. Most of the properties were found to depend on both composition and temperature. While properties, like atomic charges, show additive behaviors, other properties, such as electrical conductivity, show considerable deviations from additivity. We shall show that the mixing of the molten salt mixtures is mainly driven by entropy, and that the KCl and LiCl mix better than KCl and NaCl. Our computational results are in general consistent with available experimental data. Comparison with available theoretical data is also provided.

Nguyen, Manh Thuong↗

Achieving designed texture and flows in bulk active nematics using optimal control theory

Being intrinsically nonequilibrium, active materials can potentially perform functions that would be thermodynamically forbidden in passive materials. However, active systems have diverse local attractors that correspond to distinct dynamical states, many of which exhibit chaotic turbulent-like dynamics and thus cannot perform work or useful functions. Designing such a system to choose a specific dynamical state is a formidable challenge. Motivated by recent advances enabling optogenetic control of experimental active materials, we describe an optimal control theory framework that identifies a spatiotemporal sequence of light-generated activity that drives an active nematic system toward a prescribed dynamical steady state. Active nematics are unstable to spontaneous defect proliferation and chaotic streaming dynamics in the absence of control. We demonstrate that optimal control theory can compute activity fields that redirect the dynamics into a variety of alternative dynamical programs and functions. This includes dynamically reconfiguring between states, selecting and stabilizing emergent behaviors that do not correspond to attractors, and are hence unstable in the uncontrolled system. Furthermore, our results provide a roadmap to leverage optical control methods to rationally design structure, dynamics, and function in a wide variety of active materials.

Complex systems theory↗

Structure and thermodynamics of water adsorption in NU-1500-Cr

Abstract Metal-organic frameworks (MOFs) are a class of materials with diverse chemical and structural properties, and have been shown to effectively adsorb various types of guest molecules. The mechanism of water adsorption in NU-1500-Cr, a high-performance atmospheric water harvesting MOF, is investigated using a combination of molecular dynamics simulations and infrared spectroscopy. Calculations of thermodynamic and dynamical properties of water as a function of relative humidity allow for following the adsorption process from the initial hydration stage to complete filling of the MOF pores. Initial hydration begins at the water molecules that saturate the open Cr 3+ sites of the framework, which is then followed by the formation of water chains that extend along the channels connecting the hexagonal pores of the framework. Water present in these channels gradually coalesces and fills the hexagonal pores sequentially after the channels are completely hydrated. The development of hydrogen-bond networks inside the MOF pores as a function of relative humidity is characterized at the molecular level using experimental and computational infrared spectroscopy. A detailed analysis of the OH-stretch vibrational band indicates that the low-frequency tail stems from strongly polarized hydrogen-bonded water molecules, suggesting the presence of some structural disorder in the experimental samples. Strategies for designing efficient water harvesting MOFs are also proposed based on the mechanism of water adsorption in NU-1500-Cr.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

catalysis (heterogeneous)↗

Chemical Modeling of the Capture of Radioiodine Species

Iodine is a common fission product resulting from the transmutation of uranium fuel in nuclear power reactors. The short-lived radioactive isotope, 131 I, with a half-life of around eight days, would be a particular health concern if released into the environment. Capture materials, such as activated carbon, are used routinely in nuclear power plants and medical isotope production facilities to remove this hazard during normal operations to reduce emissions to acceptable levels. Iodine, unlike other common fission products such as noble gases, is challenging because of the wide range of molecular forms and oxidation states (from +VII in IF 7 to -I in iodide) that it displays, which can interconvert both in the gas phase and in the solution phase for example. The work in this project is focused on three computational tasks. In the first, we consider the binding of iodine-containing molecules of varying oxidation states to the major component of activated carbon – pristine graphene-like sheets. The second task focused on molecular binding and possible chemical transformations resulting from the reaction of molecular iodine with chemical functional groups on the edge of the graphene sheets or surface defects in the sheet itself. The final task was associated with radiolytic breakdown of 131 I to 131 Xe in iodine containing molecules and investigated the resulting molecular fragmentation process both in the gas phase and on the graphene-like surface. For each of the tasks, the intent was to learn about the thermodynamics and kinetics of the chemical transformations and use the calculations to provide parameters for input into kinetic models describing these processes for interpretation of experimental data in future work.

07 ISOTOPE AND RADIATION SOURCES↗

Transfer Learning Meets Embedded Correlated Wavefunction Theory for Chemically Accurate Molecular Simulations: Application to Calcium Carbonate Ion Pairing

Achieving chemical accuracy for molecular simulations remains a central challenge in computational chemistry. Here, we present an embedded correlated wavefunction transfer learning (ECW-TL) framework for accurately simulating molecular dynamics in the condensed phase. ECW-TL incorporates high-level electron exchange and correlation effects in ECW theory while preserving the training and computational efficiency of machine-learned interatomic potentials. We demonstrate the framework on Ca 2+ –CO 3 2– ion pairing in aqueous solution, a key process underlying CO 2 mineralization in seawater. As proof of principle, we first show that fine-tuning a DFT-revPBE-D3(BJ) baseline model with embedded-DFT-SCAN data reproduces the DFT-SCAN free-energy surface within 1 kcal/mol across all solvation states. Extending the framework to embedded MP2 and localized natural-orbital CCSD(T) further refines the free-energy profile, revealing the crucial role of exact electron exchange and correlation in determining ion-pair stability and structure. The computed ion-pair association free energy is in quantitative agreement with experimental measurements, further validating the accuracy of the ECW-TL framework. ECW-TL thus provides a general, data-efficient route for transferring CW accuracy to efficient simulations of complex aqueous and interfacial chemical processes.

cluster chemistry↗