Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “entropy functional”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A grid-independent approximate Riemann solver with applications to the Euler and Navier-Stokes equations

A new two-dimensional approximate Riemann solver has been developed that obtains fluxes on grid faces via wave decomposition. By utilizing information propagation in the velocity-difference directions rather than in the grid-normal directions, this flux function more appropriately interprets and hence more sharply resolves shock and shear waves when they lie oblique to the grid. The model uses five waves to describe the difference in states at a grid face. Two acoustic waves, one shear wave, and one entropy wave propagate in the direction defined by the local velocity difference vector, while the fifth wave is a shear wave that propagates at a right angle to the other four. Test cases presented include a shock reflecting off a wall, a pure shear wave, supersonic flow over an airfoil, and viscous separated airfoil flow. Results using the new model give significantly sharper shock and shear contours than a grid-aligned solver. Navier-Stokes computations over an aifoil show reduced pressure distortions in the separated region as a result of the grid-independent upwinding.

Rumsey, Christopher L.↗
DOE OSTI2022

Hybrid Nanocrystals of Small Molecules and Chemically Disordered Polymers

Organic crystals formed by small molecules can be highly functional but are often brittle or insoluble structures with limited possibilities for use or processing from a liquid phase. A possible solution is the nanoscale integration of polymers into organic crystals without sacrificing long-range order and therefore function. This enables the organic crystals to benefit from the advantageous mechanical and chemical properties of the polymeric component. We report here on a strategy in which small molecules co-crystallize with side chains of chemically disordered polymers to create hybrid nanostructures containing a highly ordered lattice. Synchrotron X-ray scattering, absorption spectroscopy, and coarse-grained molecular dynamics simulations reveal that the polymer backbones form an “exo-crystalline” layer of disordered chains that wrap around the nanostructures becoming a handle for interesting properties. The morphology of this “hybrid bonding polymer” nanostructure is dictated by the competition between the polymers’ entropy and the enthalpy of the lattice allowing for control over the aspect ratio of the nanocrystal by changing the degree of polymer integration. Further, we observed that nanostructures with an exo-crystalline layer of polymer exhibit enhanced fracture strength, self-healing capacity, and dispersion in water, which benefits their use as light-harvesting assemblies in photocatalysis. Guided by computation, future work could also explore these hybrid nanostructures as building blocks for functional materials.

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

Entropy Stabilization Effects and Ion Migration in 3D “Hollow” Halide Perovskites

A recently discovered new family of 3D halide perovskites with the general formula (A) 1–x (en) x (Pb) 1–0.7x (X) 3–0.4x (A = MA, FA; X = Br, I; MA = methylammonium, FA = formamidinium, en = ethylenediammonium) is referred to as “hollow” perovskites owing to extensive Pb and X vacancies created on incorporation of en cations in the 3D network. The “hollow” motif allows fine tuning of optical, electronic, and transport properties and bestowing good environmental stability proportional to en loading. To shed light on the origin of the apparent stability of these materials, we performed detailed thermochemical studies, using room temperature solution calorimetry combined with density functional theory simulations on three different families of “hollow” perovskites namely en/FAPbI 3 , en/MAPbI 3 , and en/FAPbBr 3 . We found that the bromide perovskites are more energetically stable compared to iodide perovskites in the FA-based hollow compounds, as shown by the measured enthalpies of formation and the calculated formation energies. The least stable FAPbI 3 gains stability on incorporation of the en cation, whereas FAPbBr 3 becomes less stable with en loading. This behavior is attributed to the difference in the 3D cage size in the bromide and iodide perovskites. Configurational entropy, which arises from randomly distributed cation and anion vacancies, plays a significant role in stabilizing these “hollow” perovskite structures despite small differences in their formation enthalpies. With the increased vacancy defect population, we have also examined halide ion migration in the FA-based “hollow” perovskites and found that the migration energy barriers become smaller with the increasing en content.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exploring the First High-Entropy Thin Film Libraries: Composition Spread-Controlled Crystalline Structure

Thin films of two types of high-entropy oxides (HEOs) have been deposited on 76.2 mm Si wafers using combinatorial sputter deposition. In one type of the oxides, (MgZnMnCoNi)O x , all the metals have a stable divalent oxidation state and similar cationic radii. In the second type of oxides, (CrFeMnCoNi)O x , the metals are more diverse in the atomic radius and valence state, and have good solubility in their sub-binary and ternary oxide systems. Therefore, the resulting HEO thin films were characterized using several high-throughput analytical techniques. The microstructure, composition, and electrical conductivity obtained on defined grid maps were obtained for the first time across large compositional ranges. The crystalline structure of the films was observed as a function of the metallic elements in the composition spreads, that is, the Mn and Zn in (MgZnMnCoNi)O x and Mn and Ni in (CrFeMnCoNi)O x . The (MgZnMnCoNi)O x sample was observed to form two-phase structures, except single spinel structure was found in (MgZnMnCoNi)O x over a range of Mn > 12 at. % and Zn < 44 at. %, while (CrFeMnCoNi)O x was always observed to form two-phase structures. Composition-controlled crystalline structure is not only experimentally demonstrated but also supported by density function theory calculation.

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

Accuracy in Resolving the First Hydration Layer on a Transition-Metal Oxide Surface: Experiment (AP-XPS) and Theory

Understanding the equilibrium conditions at the metal oxide/aqueous interface is a key component toward visualizing the structure of water in confined environments and differentiating the catalytic activity of transition-metal oxides. While ambient pressure X-ray photoelectron spectroscopy (AP-XPS) has been the primary technique to investigate the formation of a hydration layer on many surfaces, results over the extended relative humidity (RH) range accessible experimentally have not been compared quantitatively to theoretical predictions. With the use of first-principles theoretical methods and accumulated knowledge of AP-XPS spectral analysis, we do so here for a model surface, TiO 2 -terminated undoped SrTiO 3 (100) (STO). The measured distribution of OH and H 2 O coverages from vacuum up to the first hydration layer is in good agreement with a static density functional theory (DFT) configuration involving partial dissociation of H 2 O per Ti-atom mediated by H-bonding. Furthermore, ab initio molecular dynamics (AIMD) simulations at 300 K for select coverages (1/4, 1/2, and 1 ML) test the role of fluctuations and entropy in the competition between adsorption and dissociation with coverage. This comparison between theory and experiment for OH and H 2 O coverages on STO provides a foundation for a more quantitative assessment of the first hydration layer and associated competition between adsorption, dissociation, and H-bonding on transition-metal oxide surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Equilibrium, stability, and orbital evolution of close binary systems

We present a new analytic study of the equilibrium and stability properties of close binary systems containing polytropic components. Our method is based on the use of ellipsoidal trial functions in an energy variational principle. We consider both synchronized and nonsynchronized systems, constructing the compressible generalizations of the classical Darwin and Darwin-Riemann configurations. Our method can be applied to a wide variety of binary models where the stellar masses, radii, spins, entropies, and polytropic indices are all allowed to vary over wide ranges and independently for each component. We find that both secular and dynamical instabilities can develop before a Roche limit or contact is reached along a sequence of models with decreasing binary separation. High incompressibility always makes a given binary system more susceptible to these instabilities, but the dependence on the mass ratio is more complicated. As simple applications, we construct models of double degenerate systems and of low-mass main-sequence star binaries. We also discuss the orbital evoltuion of close binary systems under the combined influence of fluid viscosity and secular angular momentum losses from processes like gravitational radiation. We show that the existence of global fluid instabilities can have a profound effect on the terminal evolution of coalescing binaries. The validity of our analytic solutions is examined by means of detailed comparisons with the results of recent numerical fluid calculations in three dimensions.

Lai, Dong↗
DOE OSTI2024

Bayesian blacksmithing: discovering thermomechanical properties and deformation mechanisms in high-entropy refractory alloys

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Correlating Li-Ion Solvation Structures and Electrode Potential Temperature Coefficients

Temperature coefficients (TCs) for either electrochemical cell voltages or potentials of individual electrodes have been widely utilized to study the thermal safety and cathode/anode phase changes of lithium (Li)-ion batteries. However, the fundamental significance of single electrode potential TCs is little known. Here, we discover that the Li-ion desolvation process during Li deposition/intercalation is accompanied by considerable entropy change, which significantly contributes to the measured Li/Li + electrode potential TCs. To explore this phenomenon, we compare the Li/Li + electrode potential TCs in a series of electrolyte formulations, where the interaction between Li-ion and solvent molecules occurs at varying strength as a function of both solvent and anion species as well as salt concentrations. As a result, we establish correlations between electrode potential TCs and Li-ion solvation structures and further verify them by ab initio molecular dynamics simulations. We show that measurements of Li/Li + electrode potential TCs provide valuable knowledge regarding the Li-ion solvation environments and could serve as a screening tool when designing future electrolytes for Li-ion/Li metal batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1972

Melting of iron by significant structure theory

From Eyring's method of significant structures, a partition function is derived for liquid iron. The solid at high temperature is described in the Einstein approximation. Magnetic and electronic contributions to the thermodynamic properties of both the liquid and solid phases are considered. The model is compatible with properties (thermal expansion, compressibility, heat capacity, entropy of melting, and volume change on melting) at one atm. The melting temperature at high pressure is found by satisfying the requirement that the Gibbs free energies of the liquid and solid phases are equal at the melting temperature. Under conditions at the earth's core-mantle boundary, the melting temperature of iron is greater than approximately 5000 K, and under inner-outer core conditions the melting temperature is greater than approximately 7000 K. These estimates are consistent with the Lindemann melting law, but not with the Kraut-Kennedy melting law.

Leppaluoto, D. A.↗
DOE OSTI2025

Switching of Magnetic Order via Non‐Magnetic Al Addition in FeCoNiMnAl x Films

Magnetic high entropy alloys (HEAs) consisting of 3 d transition metals offer an exciting platform to explore novel magnetic phases as they often house competing exchange interactions in combination with random site disorders. In this work, a sensitive and tunable magnetic order is demonstrated in sputtered single-layer FeCoNiMnAl x films, as a function of non-magnetic Al addition, along with an unexpected exchange bias effect. Thin films of 50 nm FeCoNiMn exhibit a face-centered-cubic ( fcc ) phase, reentrant spin glass (SG) transition near 100 K, and a large exchange bias of over 500 Oe after field-cooling to 5 K. The exchange bias is increased to 930 Oe through a small addition of 5 at.% Al. Further Al addition to 12 at.% results in a body-centered-cubic ( bcc ) phase, coinciding with a large increase in the saturation magnetization, decrease of exchange bias to 50 Oe, and suppression of SG behavior. The change in magnetic order across the Al-induced structural transformation is mediated by the switching of Mn ground state from AF to FM, which is supported by first-principles calculations and experimentally confirmed via X-ray magnetic circular dichroism. These results open up new HEA strategies for explorations of novel magnetic phases.

competing exchange interactions↗
DOE OSTI2026

A theoretical study of solid solution strengthening in the refractory medium entropy alloy Nb 45 Ta 25 Ti 15 Hf 15

The refractory medium-entropy alloy (RMEA) Nb 45 Ta 25 Ti 15 Hf 15 exhibits exceptional tensile ductility and fracture toughness at ambient temperature, but its engineering applications are limited by a lack of high temperature strength. Using a machine-learning interatomic potential (MLIP) with near-density functional theory (DFT) accuracy, we conducted molecular dynamics (MD) and statics simulations of the behavior of dislocations with both screw and edge characters. We also analyze experimentally measured yield strengths using the Rao-Suzuki model and the Maresca-Curtin model modified to include a temperature-dependent shear modulus and a bulk modulus-dependent misfit volume, thereby uncovering the mechanisms underlying the yielding of this RMEA. Compared with the published experimental yield strength, the models parameterized by the MLIP effectively reproduce the experimental results over a wide temperature range. The models and MD simulations indicate that yielding is governed by screw dislocations, with dipole dragging as the dominant mechanism. In MD simulations, we observed a potential softening mechanism not considered by the Rao-Suzuki screw model: slow migration of interstitial jogs along the dislocation core, which could lead to the annihilation of vacancy and interstitial jog pairs by their combination.

BCC complex concentrated alloys↗
DOE OSTI2024

Design principles for energy transfer in the photosystem II supercomplex from kinetic transition networks

Photosystem II (PSII) has the unique ability to perform water-splitting. With light-harvesting complexes, it forms the PSII supercomplex (PSII-SC) which is a functional unit that can perform efficient energy conversion, as well as photoprotection, allowing photosynthetic organisms to adapt to the naturally fluctuating sunlight intensity. Achieving these functions requires a collaborative energy transfer network between all subunits of the PSII-SC. In this work, we perform kinetic analyses and characterise the energy landscape of the PSII-SC with a structure-based energy transfer model. With first passage time analyses and kinetic Monte Carlo simulations, we are able to map out the overall energy transfer network. We also investigate how energy transfer pathways are affected when individual protein complexes are removed from the network, revealing the functional roles of the subunits of the PSII-SC. In addition, we provide a quantitative description of the flat energy landscape of the PSII-SC. We show that it is a unique landscape that produces multiple kinetically relevant pathways, corresponding to a high pathway entropy. These design principles are crucial for balancing efficient energy conversion and photoprotection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Computational insights into hydrogen adsorption energies on medium-entropy oxides

High entropy oxides (HEOs) have emerged as promising catalysts for several important chemical transformations including alkane activation. Hydrogen adsorption energy (HAE) has been used as a key descriptor for many reactions including methane C–H activation and hydrogen evolution reactions. Hence, understanding the relationship between HAEs and the surface chemistry of HEO surfaces could lay the foundation for meaningful correlations among methane C–H activation, HAE, and the complex, local environment of HEO surfaces. Here, we used a medium-entropy oxide as a prototypical system – Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O with a rock-salt structure – to interrogate these relationships. We sampled 2000 different surfaces of its (100) plane and calculated the HAEs at randomly chosen surface O sites using density functional theory (DFT). Our analysis of the 2000 data points reveals that the HAEs at the surface O sites are significantly influenced by the local environment around the adsorption sites, particularly the nature of the metal atom directly below the surface O site where H adsorbs. After comparing several popular graph-neural-network-based machine learning models, we found that the DimeNet++ model performed best achieving satisfactory accuracy in predicting HAEs for both Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O and slightly varied compositions. Our work underscores the promise of such models and the need for further refinement to address the complexity of HEOs.

Song, Haohong [Vanderbilt Univ., Nashville, TN (Un↗
NASA NTRS1975

Computational procedure for three-dimensional boundary layers on aircraft and aerospace vehicles

An implicit numerical procedure is developed to solve the full three-dimensional laminar, boundary layer equations over realistic aircraft and spacecraft configurations. This procedure utilizes splined functions in the direction normal to the surface and finite differences in the cross flow direction. Calculation of off-pitch plane attachment and detachment lines can be calculated. The solution procedure permits the incorporation of a wide variety of turbulent models and boundary conditions including; entropy layer procedures, specified surface temperatures or heat flux distributions, surface catalysis, mass addition, and the coupling to both three-dimensional inviscid codes and surface ablation codes. In addition, reasonable computing times are obtained.

Kendall, R. M.↗
NASA NTRS1984

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie, close to the Hugoniot. Calcite 3 to 2 transition, upon release, was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90% devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗
NASA NTRS1986

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie close to the Hugoniot. Calcite 3 to 2 transition, upon release was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90 percent devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗
DOE OSTI2021

Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations

Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, the formation of ice XI is kinetically hindered, and experimentally it is obtained by doping with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping with a bias potential that enhances the formation of ionic defects in molecular dynamics simulations. The recombination of the ions thus formed proceeds through fast migration of the hydroxide along hydrogen bond loops, providing a physical and expedite way to change the proton configuration. Here, a key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modelling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter that distinguishes the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of the calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Probing configurational disorder in ZnGeN 2 using cluster-based Monte Carlo

ZnGeN 2 is sought as a semiconductor with comparable lattice constant to GaN and tunable band gap for integration in optoelectronic devices. Configurational disorder on the cation sublattice of ZnGeN 2 can strongly modify the electronic structure compared to the ordered material, and both ordered and disordered forms of ZnGeN 2 are candidates for light-emitting diodes and other emitters. The nonisovalent character of the disordered species (Zn 2+ and Ge 4+ ) subjects the cation ordering to strong short-range order effects. To model these effects, we use Monte Carlo (MC) simulations utilizing a cluster expansion to approximate formation enthalpy. Representative disordered configurations in 1024-atom supercells are relaxed in density functional theory calculations. From the MC structures, we extract a short-range order parameter (the N-cation coordination motif), and two long-range order parameters (Bragg-Williams and stretching parameters), and examine their correlations. We perform a thermodynamic integration to determine the mixing entropy and free energy. ZnGeN 2 exhibits a first-order phase transition with pronounced discontinuities in enthalpy and entropy, as well as in the structural order parameters. We discuss the relationship between the effective temperature used in the MC simulation and the growth temperatures in experiment in relation to the crossover from the nonequilibrium to the equilibrium growth regime. This work expands on current models of site disorder in ZnGeN 2 and provides atomic structure models with a systematic variation of the degree of cation disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗