Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermodynamics of computation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Experimental and DFT Calculated IR Spectra of Guests in Zeolites: Acyclic Olefins and Host-Guest Interactions

We performed experimental and periodic density functional theory (DFT) IR spectroscopy to investigate the adsorption of acyclic olefins over both acidic and non-acidic zeolites. Two conjugated polyenes, 2,4-dimethyl-1,3-pentadiene (I) and 2,6-dimethyl-2,4,6-octatriene (II) were studied to probe organic intermediates that can form during methanol conversion, and that can lead to deactivating species known collectively as “coke.” We computed vibrational spectra using zeolite-adsorbed and gas-phase models for both neutral and protonated forms of I and II, and compared these DFT results to diffuse reflectance IR Fourier transform (DRIFT) spectra of zeolite-guest systems. Our experimental and computational results are precise enough to pinpoint that the gauche s-cis conformation of species I predominates during adsorption over de-aluminated zeolite beta. Computed zeolite-adsorbed spectra of the protonated species I and II best represent the DRIFT spectra obtained after the adsorption of the olefins on HMOR at 20 ?, with computed bands at 1543 and 1562 cm–1 for molecules I+ and II+, respectively, attributed to the allylic stretching mode, ?(C=C–C+). These computed band frequencies are within 6 cm–1 of experimental data and confirm that the interaction between neutral acyclic olefins and acidic zeolites leads to protonation of the olefin. Comparison of computed spectra of the protonated species in the gas phase to those in the zeolite indicates that the electrostatic interaction between alkenyl and alkadienyl cations and negative zeolite framework does not significantly impact the position of the allylic stretching bands. These results highlight that computed spectroscopy and thermodynamics coupled with experimental spectra can be used for elucidating complex mixtures in zeolites, and certain spectral features of adsorbed olefins can be accurately modeled by gas-phase calculations. ACKNOWLEDGMENTS BM and SMA were supported by the National Science Foundation under Award # CBET-1512442. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund for support (or partial support) of this research (#589978-ND5, EDH and FCJ). MDB was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Materials Science and Engineering. CJM was supported by the DOE, Office of Science, BES, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated by Battelle for the U.S. DOE. Computing resources were generously provided through PNNL’s institutional computing (PIC).

Manookian, Babgen↗
DOE OSTI2022

Quantification of Dead Lithium on Graphite Anode under Fast Charging Conditions

A series of computational and experimental studies were conducted to understand the onset of lithium plating and subsequent quantification of dead lithium on graphite electrodes in the design of fast charging batteries. The experiments include titration and relaxation studies for detecting initiation of lithium metal plating for various SOC and C-rates, which are compared against the thermodynamically consistent phase field computational results. The collaborative study on "model graphite electrode" with 2.18 mAh cm-2 nominal capacity at 25 degrees C demonstrates: (1) the macroscopic voltage response during relaxation studies indicate the reintercalation of plated lithium into the graphite anode; (2) for SOC below 60% and low C-Rates, there is no dead lithium; (3) for SOC between 60% to 80%, and C-Rates in the range of 4C-6C show dead lithium both in experiments and simulations; (4) at 100% SOC and 4C-6C rates, large amounts of dead lithium are observed. The study presented here allows us to evaluate the effects of the physical properties of the electrochemical system on plating and stripping kinetics and the amount of dead lithium on graphite electrodes, which determines the cell capacity loss under fast charge.

batteries↗
DOE OSTI2023

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

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

Cantor-derived medium-entropy alloys: bridging the gap between traditional metallic and high-entropy alloys

The year 2004 marked the beginning of a new era in the design of metallic materials, as the concept of multiple principal component alloys, commonly known as High-Entropy Alloys (HEAs), was proposed by Cantor and Yeh. The unexpected single-phase microstructure, instead of the expected brittle intermetallic compounds, was attributed to the large entropy of mixing and immediately caught the attention of the scientific community. Today, HEAs are considered important advanced materials and a broad range of alloys using nominally the same design principle have been investigated. Despite that, the CrMnFeCoNi (Cantor) alloy stands out as the most successful HEA due to its outstanding mechanical properties and microstructure. In this scenario, variants of the Cantor alloy, named medium-entropy alloys (MEAs), are gaining significant interest as they display a better industrial potential than both HEAs and traditional alloys. These variants of the Cantor alloy with only three or four main elements result in 15 possible combinations. The microstructure of these alloys is discussed in terms of advanced characterization as well as thermodynamic parameters and computational simulation. Their phase stability is addressed over a wide range of temperatures and strain rates. The mechanical properties, especially the fracture toughness, of the CrFeCoNi and CrCoNi alloys have been reported to be even superior to those of the Cantor alloy and most modern engineering alloys. This is associated with the formation of a continuous sequence of strengthening mechanisms, including hierarchical twin networks, which serve to prolong the strain hardening. The present article reviews and critically assesses, for the first time, recent advances in these Cantor-derived MEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Computational design of passivants for CdTe grain boundaries

CdTe is the second most-widely deployed photovoltaic (PV) material due to its high efficiency and low manufacturing costs. Currently, polycrystalline CdTe (poly-CdTe) has a record efficiency of similar to 22%, which is still well below the theoretical limit (similar to 30%). Polycrystalline CdTe films that have incorporated Cl or Se show higher efficiency, possibly due to the segregation of these ions to the grain boundaries (GBs), where they may passivate the dangling bonds. However, the efficiency enhancement mechanisms of passivants in CdTe GBs, and the feasibility of employing alternative passivants, have not been well explored. Here, we present a systematic computational study of CdTe GBs with potential passivants, namely S, P, As, Se, and Sb on Te sites, and Na, Mg, Al, Sc, Cu, and Zn on Cd sites. Density functional theory (DFT) calculations were performed on GB dislocation core structures derived from scanning transmission electron microscopy (STEM) images of a model GB (bicrystal). We computed the segregation thermodynamics, electronic density of states, and charge variations near doped CdTe GBs. We find that segregation of impurities to GBs is thermodynamically favorable. For a Te-terminated core, Se, S, and P on Te sites effectively reduce midgap states. For both Cd-and Te-terminated dislocation cores, Sc and Al reduce midgap states when substituted for Cd atoms. Overall, the greatest improvement was achieved with co-doping, i.e. simultaneously substituting Te with Se and substituting Cd with Cu or Al. The elimination of midgap states is predicted to increase the photovoltaic efficiency of CdTe by reducing the recombination at grain boundaries.

14 SOLAR ENERGY↗
DOE OSTI2024

Simulations of radiatively cooled magnetic reconnection driven by pulsed power

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P 1/3 multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number M A ≈ 1.5), supersonic (Sonic Mach number M S ≈ 4–5) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio L/δ≈100 , Lundquist number S L ≈ 400). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = τ$^{–1}_{cool}$/τ$^{–1}_{H}$ ≈100); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Lastly, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Synergistic Enhancement of Lead and Selenate Uptake at the Barite (001)–Water Interface

The interactions of heavy metals with minerals influence the mobility and bioavailability of toxic elements in natural aqueous environments. The sorption of heavy metals on covalently bonded minerals is generally well described by surface complexation models (SCMs). However, understanding sorption on sparingly soluble minerals is challenging because of the dynamically evolving chemistry of sorbent surfaces. The interpretation can be even more complicated when multiple metal ions compete for sorption. In the present study, we observed synergistically enhanced uptake of lead and selenate on the barite (001) surface through two sorption mechanisms: lattice incorporation that dominates at lower coverages and two-dimensional monolayer growth that dominates at higher coverages. Further, we also observed a systematic increase in the sorption affinity with increasing co-sorbed ion coverages, different from the assumption of invariant binding constants for individual adsorption processes in classical SCMs. Computational simulations showed thermodynamically favorable co-incorporation of lead and selenate by simultaneously substituting for barium and sulfate in neighboring sites, resulting in the formation of molecular clusters that locally match the net dimension of the substrate lattice. These results emphasize the importance of ion–ion interactions at mineral–water interfaces that control the fate and transport of contaminants in the environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Ca Cobaltites as Potential Cathode Materials for Rechargeable Ca-Ion Batteries: Theory and Experiment

Rechargeable Ca-based batteries can potentially achieve higher energy capacity than Li-ion batteries. The development of Ca-ion batteries, however, remains in its infancy, especially due to the challenge of finding cathode materials with high reversible capacity. In this work, we investigate properties of Ca cobaltites as Ca-ion intercalation cathodes, by means of density functional theory calculations validated by synthesis and electrochemical measurements. Computationally, we accessed thermodynamic, diffusion properties, energy capacity, and the voltage profiles for four Ca cobaltite compounds of different stoichiometry and Co oxidation states: Ca 3 Co 2 O 6 , CaCo 2 O 4 , Ca 2 Co 2 O 5 , and [Ca 2 CoO 3 ][CoO 2 ] 1.62 . We found good stability and relatively low migration barriers of some Ca cobaltites during cycling with the layered CaCo 2 O 4 having the lowest Ca migration barrier of 0.7 eV and the highest theoretical capacity. To validate our calculations, we synthesized Ca 3 Co 2 O 6 , CaCo 2 O 4 , and [Ca 2 CoO 3 ][CoO 2 ] 1.62 with the Pechini method and, subsequently, to test the electrochemical extraction of Ca.

25 ENERGY STORAGE↗
DOE OSTI2026

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dual dressed black holes as the end point of the charged superradiant instability in $\mathcal{N} = 4$ Yang Mills

Charged Black holes in AdS_5 × S^5 A d S 5 × S 5 suffer from superradiant instabilities over a range of energies. Hairy black hole solutions (constructed within gauged supergravity) have previously been proposed as endpoints to this instability. We demonstrate that these hairy black holes are themselves unstable to the emission of large dual giant gravitons. We propose that the endpoint to this instability is given by Dual Dressed Black Holes (DDBH)s; configurations consisting of one, two, or three very large dual giant gravitons surrounding a core AdS A d S black hole with one, two, or three SO(6) S O ( 6 ) chemical potentials equal to unity. The dual giants each live at AdS A d S radial coordinates of order \sqrt{N} N and each carry charge of order N^2 N 2 . The large separation makes DDBHs a very weakly interacting mix of their components and allows for a simple computation of their thermodynamics. We conjecture that DDBHs dominate the phase diagram of \mathcal{N}=4 𝒩 = 4 Yang-Mills over a range of energies around the BPS plane, and provide an explicit construction of this phase diagram, briefly discussing the interplay with supersymmetry. We develop the quantum description of dual giants around black hole backgrounds and explicitly verify that DDBHs are stable to potential tunneling instabilities, precisely when the chemical potentials of the core black holes equal unity. We also construct the 10-dimensional DDBH bulk solutions.

Choi, Sunjin (ORCID:0000000283310199)↗
DOE OSTI2026

Harnessing Ocean Thermal Gradients Using Thermoelectric Based Submersibles for Ocean Power Applications

The urgent need for energy solutions in marine environments has accelerated the development of innovative technologies capable of leveraging natural resources for power generation. This study introduces a buoyancy-driven submersible system designed to harness ocean thermal gradients using thermoelectric generators (TEGs) and phase change materials (PCMs). The technology aims to provide autonomous power to offshore aquaculture farms, unmanned underwater vehicles (UUVs), offshore platform illumination, and ocean sensors, significantly reducing dependence on fossil fuels. Ocean thermal gradients, especially prevalent in mid-latitude regions, exhibit temperature differences between surface and deep waters ranging from 7 degrees Celsius to 30 degrees Celsius depending on seasonal variations. The proposed submersible technology utilizes TEGs to convert thermal energy from these gradients into electrical power, generating between 0.2 and 0.5 watts, while PCMs are employed to store and regulate this energy, ensuring a stable and continuous power supply. The buoyancy-driven mechanism of the submersible enhances its capability to navigate through varying depths, optimizing its exposure to different thermal gradients and maximizing energy harvesting. The performance of this submersible system is analyzed through detailed thermodynamic assessments and computational fluid dynamics (CFD) modeling focused on heat transfer. These analyzes consider real-world ocean temperature profiles and seek to refine the interaction between TEGs and PCMs to optimize energy extraction. The evaluation encompasses several key performance metrics, including power output and energy efficiency. Results confirm the potential of this innovative technology to provide a continuous and reliable power source for marine applications. By demonstrating the feasibility of using ocean thermal gradients for energy generation, this study contributes to the broader efforts of innovation in energy technologies for harsh, remote marine environments. The implementation of such promises is significant advancements in the autonomy of marine operations. The ongoing research will further investigate scalability ensuring that the system can be effectively adapted to various marine settings and operational demands.

16 TIDAL AND WAVE POWER↗
NASA NTRS1970

A program for computing shock-tube gas dynamic properties

Computer program calculates thermodynamic properties from basic spectroscopic data. Program capacity is a mixture of 100 different species composed of ten different elements. The output is a complete thermodynamic and chemical description of the gas.

Horton, T. E.↗
NASA NTRS1972

Interrelated structure of high altitude atmospheric profiles

A preliminary development of a mathematical model to compute probabilities of thermodynamic profiles is presented. The model assumes an exponential expression for pressure and utilizes the hydrostatic law and equation of state in the determination of density and temperature. It is shown that each thermodynamic variable can be factored into the produce of steady state and perturbation functions. The steady state functions have profiles similar to those of the 1962 standard atmosphere while the perturbation functions oscillate about 1. Limitations of the model and recommendations for future work are presented.

Engler, N. A.↗
NASA NTRS1972

Optimization of thermal protection systems for the space shuttle vehicle. Volume 1: Final report

A study performed to continue development of computational techniques for the Space Shuttle Thermal Protection System is reported. The resulting computer code was used to perform some additional optimization studies on several TPS configurations. The program was developed in Fortran 4 for the CDC 6400, and it was converted to Fortran 5 to be used for the Univac 1108. The computational methodology is developed in modular fashion to facilitate changes and updating of the techniques and to allow overlaying the computer code to fit into approximately 131,000 octal words of core storage. The program logic involves subroutines which handle input and output of information between computer and user, thermodynamic stress, dynamic, and weight/estimate analyses of a variety of panel configurations. These include metallic, ablative, RSI (with and without an underlying phase change material), and a thermodynamic analysis only of carbon-carbon systems applied to the leading edge and flat cover panels. Two different thermodynamic analyses are used. The first is a two-dimensional, explicit precedure with variable time steps which is used to describe the behavior of metallic and carbon-carbon leading edges. The second is a one-dimensional implicity technique used to predict temperature in the charring ablator and the noncharring RSI. The latter analysis is performed simply by suppressing the chemical reactions and pyrolysis of the TPS material.

Source record↗
NASA NTRS1973

Effect of downstream massive blowing on Jovian entry heating.

An analysis for viscous radiating shock layer flow over a blunt body with massive blowing is presented. Detailed equilibrium computations for the thermodynamics and transport properties are included in the solution. Molecular band and atomic line, as well as continuum radiation are included. For a typical Jovian entry problem, a 40% reduction in total radiative flux at stagnation point, due to the injection of carbon gas, was found. The reduction increases to about 70% at about four body nose radii downstream. The increase of flux reduction indicates that most of the radiative energy is absorbed by the injected carbon gas and dumped into the wake. An examination of the spectral distribution of the surface radiative flux reveals that most of the flux reduction occurs in molecular bands of C2.

Chou, Y. S.↗