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 829 records · Page 46

Thermodynamic Trends for Reduction of CO by Molecular Complexes

Selective reduction of CO 2 into fuels and chemical feedstocks is highly desirable to reduce our dependence on fossil fuels. Most molecular catalysts afford 2e – reduction products, such as CO or HCO 2 – , as opposed to more reduced products. Here we present an analysis of the thermodynamic limitations for reduction of the CO ligand in the form of a series of isostructural group 6 carbonyl complexes, Cp*M(CO) 3 (P(OMe) 3 ) + (M = Cr, Mo, and W). Additionally, the free energy for stepwise transfer of a hydride (H – ) and a proton (H + ) to the CO ligand, resulting in a hydroxycarbene (CHOH) complex, was measured by equilibration with H – /H + donors and acceptors with known hydricity or acidity. Together, these two reaction steps are equivalent to a net addition of H 2 across the CO ligand. A large and unfavorable free energy for H 2 addition (Δ G ° H2 ) was measured for all three complexes and decreases in the order Cr > Mo > W. The trend for these complexes is opposite to the trend previously reported for group 7 carbonyl complexes, for which Δ G ° H2 decreases moving up the group, Re > Mn. Computational analysis indicates the trends can be described in terms of electrostatic effects, where a low Δ G ° H2 is obtained in complexes that balance the atomic charges of the M–C≡O fragment of the complex. These findings can be used to design metal carbonyl complexes with a more energetically accessible H 2 addition, which will facilitate the development of molecular catalysts for reduction of CO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermosteam: BioSTEAM's Premier Thermodynamic Engine

Thermosteam is a thermodynamic engine capable of solving mass and energy balances, estimating mixture properties, solving thermodynamic phase equilibria, and modeling stoichiometric reactions. All chemical data in Thermosteam is imported from the chemicals library, an open-source compilation of data and functions for the estimation of pure component chemical and mixture properties. Thermosteam’s fast and flexible platform has enabled the evaluation of conceptual and emerging biochemical production processes. The Biorefinery Simulation and Techno-Economic Analysis Modules (BioSTEAM) — capable of modeling reactors, distillation columns, heat exchangers, and other unit operations — has adopted Thermosteam as its premier thermodynamic engine. Published biorefinery designs modeled in BioSTEAM implement thermodynamic property packages created with Thermosteam, including a cornstover biorefinery for the production of cellulosic ethanol, a lipid-cane biorefinery for the co-production of ethanol and biodiesel, and a wheatstraw biorefinery for the production of cellulosic ethanol.

97 MATHEMATICS AND COMPUTING↗

Integrating Experiments and Simulations to Reveal Anisotropic Growth Mechanisms and Interfaces of a One-Dimensional Zeolite

Zeolites are nanoporous crystalline materials critical for diverse industrial applications, yet their growth mechanisms are poorly understood. Here, this study presents a novel integrated framework combining experimental synthesis, high-resolution imaging, coarse-grained molecular dynamics simulations, and computer vision to uncover the mechanisms of growth of SSZ-24, a 1D channel zeolite. We demonstrate how synthesis conditions, such as temperature and reactant concentration, govern crystal anisotropy and surface roughness with growth dynamics differing markedly by crystallographic orientation. Along the channels, growth involves minimal energy barriers and rapid nucleation, resulting in rough surfaces. In contrast, growth perpendicular to the channels requires cooperative molecular organization and is highly sensitive to thermodynamic and kinetic conditions, yielding smooth anisotropic surfaces under low driving forces. By simulating transmission electron microscopy (TEM) images, we bridge molecular-scale simulations with experimental observations, identifying distinct growth mechanisms along different crystal planes. This work offers molecular-level insights into zeolite crystallization, advancing the rational design of nanoporous materials. The integration of cross-disciplinary methodologies establishes a transformative framework for optimizing zeolite synthesis, with implications for broader classes of materials.

Bertolazzo, Andressa A. [Univ. of Utah, Salt Lake ↗

Multiscale modeling of metal-hydride interphases—quantification of decoupled chemo-mechanical energies

Abstract The quantification of interphase properties between metals and their corresponding hydrides is crucial for modeling the thermodynamics and kinetics of the hydrogenation processes in solid-state hydrogen storage materials. In particular, interphase boundary energies assume a pivotal role in determining the kinetics of nucleation, growth, and coarsening of hydrides, alongside accompanying morphological evolution during hydrogenation. The total interphase energy arises from both chemical bonding and mechanical strains in these solid-state systems. Since these contributions are usually coupled, it is challenging to distinguish via conventional computational approaches. Here, a comprehensive atomistic modeling methodology is developed to decouple chemical and mechanical energy contributions using first-principles calculations, of which feasibility is demonstrated by quantifying chemical and elastic strain energies of key interfaces within the FeTi metal-hydride system. Derived materials parameters are then employed for mesoscopic micromechanical analysis, predicting crystallographic orientations in line with experimental observations. The multiscale approach outlined verifies the importance of the chemo-mechanical interplay in the morphological evolution of growing hydride phases, and can be generalized to investigate other systems. In addition, it can streamline the design of atomistic models for the quantitative evaluation of interphase properties between dissimilar phases and allow for efficient predictions of their preferred phase boundary orientations.

36 MATERIALS SCIENCE↗

Hybrid Kohn-Sham+Thomas-Fermi scheme for high-temperature density functional theory

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.

36 MATERIALS SCIENCE↗

Sulfur poisoning mechanism of LSCF cathode material in the presence of SO 2 : a computational and experimental study

Aiming at the comprehensive understanding of the single sulfur poisoning effect and, eventually, the multiple impurities poisoning phenomena on the SOFC (Solid Oxide Fuel Cell) cathode materials, the sulfur poisoning effect on the (La 0.6 Sr 0.4 ) 0.95 Co 0.2 Fe 0.8 O 3 (LSCF-6428) has been investigated in the presence of 10 ppm SO 2 at 800, 900, and 1000°C, respectively, with a combined computational and experimental approach. The good agreement between the CALPHAD (Computer Coupling of Phase Diagrams and Thermochemistry) simulations and the XRD (X-Ray Diffraction), SEM (Scanning Electron Microscopy), and TEM (Transmission Electron Microscopy) characterization results support the reliability of the CALPHAD approach in the SOFC field. Furthermore, comprehensive simulations were made to understand the impact of temperature, P(SO 2 ), P(O 2 ), and Sr concentration on the threshold of SrSO 4 stability. Results showed that the formation of SrSO4 is thermodynamically favored at lower temperatures, higher P(SO 2 ), higher P(O 2 ), and higher Sr concentration. Finally, comparisons were also made between LSCF-6428 and LSM20 (La 0.8 Sr 0.2 MnO 3 ) using simulations, which confirmed that LSCF-6428 is a poor sulfur-tolerant cathode, in agreement with the literature.

36 MATERIALS SCIENCE↗

Singularity-EOS XCAP Report

Solving the Euler equations is a fundamental component of simulating many physical phenomena ranging from high explosives to astrophysics. An equation of state (EOS) is a required piece that relates any two thermodynamic quantities to all other thermodynamic values. The presence of multiple materials within a control volume further complicates the solution requiring additional equations to describe the interaction of materials at a sub-grid level. Equations of state themselves can also come in many forms ranging from simple algebraic relations to more complicated differential equation models that describe material interactions over a broad range of physical conditions. In the latter case, the EOS is often pre-computed at a given set of grid points and provided in a tabular form where additional properties can be derived from the interpolation functions.

97 MATHEMATICS AND COMPUTING↗

Optimizing Simulation Fidelity in Direct-Drive Inertial Confinement Fusion with Cassio

Recently, the National Ignition Facility (NIF) demonstrated that inertial confinement fusion (ICF) is capable to achieve thermonuclear (TN) ignition in the laboratory, making it a crucial method on the path to replicate the Sun’s power production mechanism on Earth. However, the physics governing the high-energy density environments is very complex and remains a challenge to fully understand and model. For example, dopants in the TN fuel are important diagnostic tools to extract the thermodynamic conditions of the plasma. However, if their concentration is chosen too high, they can significantly degrade the performance of an ICF capsule. In this study, we use the Los Alamos National Laboratory radiation-hydrodynamics code Cassio to model ICF implosions of capsules which contain deuterium fuel with high-Z dopants like Krypton and Argon from the high-Z campaign conducted 15 years ago. We focus on how chosen computational and physics parameters influence the implosion outcomes. By systematically changing the resolution of the computational mesh and the photon energies as well as modifying settings for the laser drive and TN fuel pre-heat effects, we assess the impact on experimentally measured performance metrics like neutron production from TN burn and x-ray emission during the implosion. Our results will help to improve the fidelity of simulations and guide future numerical studies and experimental designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermodynamic modeling of CsF with LiF-NaF-KF for molten fluoride-fueled reactors

Gibbs energy models were developed to describe the thermochemical behavior of CsF in molten FLiNaK (46.5LiF-11.5NaF-42KF mol%), a proposed molten salt reactor (MSR) fuel solvent and coolant, as cesium is of concern due to its high radiotoxicity and volatility. Initially, it was necessary to obtain a more accurate Gibbs energy function for CsF which required fitting parameters to reported vapor pressures over condensed phase CsF. The pseudo-binary systems CsF-LiF, CsF-NaF and CsF-KF were then evaluated utilizing phase equilibria and enthalpy of mixing (Δ mix H) values, together with original differential scanning calorimetry (DSC) measurements performed for the CsF-LiF and CsF-KF systems. The CsF-LiF-NaF, CsF-LiF-KF and CsF-NaF-KF pseudo-ternary system representations were obtained by interpolation of the constituent pseudo-binary systems, with DSC measurements performed for the CsF-LiF-NaF system to corroborate the calculated liquidus temperature. Ultimately, the pseudo-ternary systems were interpolated to obtain Gibbs energy models for the pseudo-quaternary CsF-LiF-NaF-KF system, supported by DSC measurements at low CsF compositions (1–10 mol%), yielding computed equilibria and cesium-containing vapor pressures that compare favorably with reported values. In conclusion, the Molten Salt Thermal Properties Database – Thermochemical (MSTDB-TC) was subsequently expanded to include these Gibbs energy models allowing description of the thermochemical behavior of the CsF-LiF-NaF-KF system.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Role of water model on ion dissociation at ambient conditions

In this work, we study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion–ion interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials

Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed here by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.

97 MATHEMATICS AND COMPUTING↗

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗

Chiral hydrodynamics in strong external magnetic fields

We construct the general hydrodynamic description of (3+1)-dimensional chiral charged (quantum) fluids subject to a strong external magnetic field with effective field theory methods. We determine the constitutive equations for the energy-momentum tensor and the axial charge current, in part from a generating functional. Furthermore, we derive the Kubo formulas which relate two-point functions of the energy-momentum tensor and charge current to 27 transport coefficients: 8 independent thermodynamic, 4 independent non-dissipative hydrodynamic, and 10 independent dissipative hydrodynamic transport coefficients. Five Onsager relations render 5 more transport coefficients dependent. We uncover four novel transport effects, which are encoded in what we call the shear-induced conductivity, the two expansion-induced longitudinal conductivities and the shear-induced Hall conductivity. Remarkably, the shear-induced Hall conductivity constitutes a novel non-dissipative transport effect. As a demonstration, we compute all transport coefficients explicitly in a strongly coupled quantum fluid via holography.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The R -Process Alliance: Exploring the cosmic scatter among ten r -process sites with stellar abundances

Context. The astrophysical origin of the rapid neutron-capture process (r-process), responsible for producing roughly half of the elements heavier than iron, remains uncertain. Detailed chemical signatures from the oldest, most metal-poor stars, which act as fossil records of the earliest nucleosynthesis events, can be used to identify the dominant r-process sites. Aims. We present a homogeneous chemical abundance analysis of ten r-process element-enhanced stars. These old and metal-poor stars are strongly enriched in r-process elements with minimal contamination from other nucleosynthetic sources. By focusing on this chemically pure sample, we aim to investigate intrinsic variations in the r-process abundance patterns and explore their implications for the nature and potential diversity of r-process sites. Methods. We performed a detailed chemical abundance analysis of high-resolution, high-signal-to-noise spectra. For each star, we inspected over 1400 individual absorption lines using a combination of equivalent width measurements and spectral synthesis. The analysis was conducted under the assumption of 1D local thermodynamic equilibrium and employing the MOOG radiative transfer code. Results. We derived abundances for 54 chemical species, including 29 neutron-capture (n-capture) elements, covering the full mass range of the r-process abundance pattern. A kinematic analysis reveals that stars likely originated from ten kinematically distinct systems. Based on this assumption, we used the sample to probe the maximum variation expected from ten independent r-process nucleosynthesis events and computed the intrinsic dispersion of each element relative to Zr and Eu for the light and heavy r-process elements, respectively. This exercise resulted in a remarkably low cosmic scatter across the ten r-process sites enriching these stars; for the rare earth and third peak elements, for example, we find σ [La/Eu] = 0.08 and σ [Os/Eu] = 0.11 dex, while the scatter between light and heavy elements, σ [Zr/Eu] , is slightly higher at 0.18 dex. Conclusions. The elemental abundance patterns across the ten independent r-process sites show remarkably small cosmic dispersions. This minimal dispersion suggests a high degree of uniformity in r-process yields across diverse astrophysical environments.

Astronomy and AstroPhysics↗

Final Scientific/Technical Report Rapid Design and Manufacturing of High-Performance Materials for Turbine Blades

This research demonstrated the concept of carbide precipitation-strengthened refractory high entropy alloys (RHEA). The advantage of a precipitation strengthened alloy is all phases are in thermodynamic equilibrium promoting microstructural stability, and consequently retention of properties at elevated temperatures. Additionally, as with any precipitation strengthened – or age hardened alloy, components can be heat-treated after manufacturing to manipulate the microstructure and optimize properties for performance. This is an advantage of the precipitation strengthened alloys over composites and ceramics materials, where the microstructure and properties are essentially fixed upon the initial materials manufacturing stage. High throughput (HT), multi-scale computer modeling was used to identifying novel RHEA compositions with desired characteristics needed for precipitation strengthening. The results showed that carbides precipitated and consequently, the strength of the alloys (measured in compression) increased after heat-treatment, which is the desired effect. The project also demonstrated the feasibility of producing articles from RHEA by additive manufacturing (AM).

36 MATERIALS SCIENCE↗