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 739 records · Page 41

Thermodynamic properties of oxygen and nitrogen III

The final equation for nitrogen was determined. In the work on the equation of state for nitrogen, coefficients were determined by constraining the critical point to selected critical point parameters. Comparisons of this equation with all the P-density-T data were made, as well as comparisons to all other thermodynamic data reported in the literature. The extrapolation of the equation of state was studied for vapor to higher temperatures and lower temperatures, and for the liquid surface to the saturated liquid and the fusion lines. A new vapor pressure equation was also determined which was constrained to the same critical temperature, pressure, and slope (dP/dT) as the equation of state. Work on the equation of state for oxygen included studies for improving the equation at the critical point. Comparisons of velocity of sound data for oxygen were also made between values calculated with a preliminary equation of state and experimental data. Functions for the calculation of the derived thermodynamic properties using the equation of state are given, together with the derivative and integral functions for the calculation of the thermodynamic properties using the equations of state. Summary tables of the thermodynamic properties of nitrogen and oxygen are also included to serve as a check for those preparing computer programs using the equations of state.

Stewart, R. B.↗

Sodium sulfate: Vaporization thermodynamics and role in corrosive flames

Gaseous species over liquid Na2SO4 were identified by the technique of molecular beam mass spectrometry. The heat and entropy of vaporization of the Na2SO4 molecule were measured directly. Comparisons of the experimental entropy with values calculated using various molecular parameters were used to estimate the molecular structure and vibrational frequencies. The thermodynamic properties of gaseous and condensed phase Na2SO4, along with additional pertinent species, were used in a computer program to calculate equilibrium flame compositions and temperatures for representative turbine engine and burner rig flames. Compositions were calculated at various fuel-to-oxidant ratios with additions of sulfur to the fuel and the components of sea salt to the intake air. Temperatures for condensation of Na2SO4 were obtained as a function of sulfur and sea salt concentrations.

Kohl, F. J.↗

Improved thermodynamic modeling of the no-vent fill process and correlation with experimental data

The United States' plans to establish a permanent manned presence in space and to explore the Solar System created the need to efficiently handle large quantities of subcritical cryogenic fluids, particularly propellants such as liquid hydrogen and liquid oxygen, in low- to zero-gravity environments. One of the key technologies to be developed for fluid handling is the ability to transfer the cryogens between storage and spacecraft tanks. The no-vent fill method was identified as one way to perform this transfer. In order to understand how to apply this method, a model of the no-vent fill process is being developed and correlated with experimental data. The verified models then can be used to design and analyze configurations for tankage and subcritical fluid depots. The development of an improved macroscopic thermodynamic model is discussed of the no-vent fill process and the analytical results from the computer program implementation of the model are correlated with experimental results for two different test tanks.

Taylor, William J.↗

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗

Radiative Heating Profiles over Ascension Island during the 2016 and 2017 Biomass Burning Seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. In many cases, the complex physical relationships between cloud morphology and the environmental conditions in which marine boundary layer clouds exist are not well understood. Such uncertainties arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean, where it is likely that the absorbing and heating properties of these aerosols modify the microphysical composition and macrophysical arrangement of marine stratocumulus and trade cumulus. The deployment of the Atmospheric Radiation Measurement Mobile Facility #1 (AMF1) in support of LASIC (Layered Atlantic Smoke Interactions with Clouds) provided a unique opportunity to observe thermodynamic, cloud and aerosol properties during two consecutive biomass burning seasons from July through October of 2016 and 2017 over Ascension Island (7.96 S, 14.35 W). These observations in conjunction with radiation transfer modeling were used to assess the impact of biomass burning carbonaceous aerosol plumes as they passed over the site.Thermodynamic profiles were generated using a combination of radiosonde data and thermodynamic profilers to provide high temporal resolution profiles of quantities that are important for cloud development, such as CAPE and CIN. Coincident Ka-band radar and lidar profiles were used to characterize the cloud and sub-cloud structure. The resulting thermodynamic and cloud profiles are used as input forcing for the Rapid Radiative Transfer Model (RRTM) to compute radiative heating profiles over the observation site. Idealized experiments using RRTM, with and without aerosols present, are used to assess the impacts of the absorbing aerosol on the heating rate profiles.

Cumulus↗

(U) Marker Lagrangian Simulations using Geometry from CAD Assembly Files [Slides]

The Eulerian code PAGOSA now contains an option for the forward time integration of any user selected material in the Lagrangian frame of reference. This option is accomplished using a Multifield version of Brackbill’s FLIP plus Sulsky’s MPM. FLIP+MPM is a Marker Lagrangian scheme that is not to be confused with well known Grid Lagrangian schemes. This Marker Lagrangian method allows materials to undergo completely arbitrary deformation, up to and including fracture and separation. Using this option makes PAGOSA into a mixed frame simulation tool such that fluids can be integrated in the Eulerian frame and solids integrated in the Lagrangian frame. Boundaries between the fields are said to be ’immersed’ in the computational domain, and their nature is formulated using Multifield Theory. Each Marker in a field of FLIP+MPM material tracks the thermodynamic state of a small (but finite) piece of mass. Initiation of the Markers using the geometry and material information from CAD Assembly files is the subject of the talk. The ease and speed with which the dynamics of the assembly can be studied is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermodynamically consistent algorithms for models of incompressible multiphase polymer solutions with a variable mobility

Here we present a general strategy for developing structure and property preserving numerical algorithms for thermodynamically consistent models of incompressible multiphase polymer solutions with a variable mobility. We first present a formalism to derive thermodynamically consistent, incompressible, multiphase polymer models. Then, we develop the general strategy, known as the supplementary variable method, to devise thermodynamically consistent numerical approximations to the models. We illustrate the numerical strategy using newly developed models of incompressible diblock copolymer solutions coupled with an electric and a magnetic field, respectively. Mesh refinement is conducted to verify convergence rates of the developed schemes. Some numerical examples are given to exhibit underlying dynamics absent from and driven by the external fields, respectively, highlighting differences between models with the variable and constant mobilities.

97 MATHEMATICS AND COMPUTING↗

Computational methods in solution-based plastics purification

Plastic waste can be recycled into resins with near-virgin properties by solution-based purification processes that selectively dissolve polymers, remove contaminants, or detach printing residues. Here, in this review, we examine computational methods for predicting the behavior governing solution-based plastic purification, motivated by the vast polymer–solvent–contaminant compositional space. We discuss thermodynamic and machine learning methods for predicting polymer–solvent and polymer–contaminant interaction and review physics-based molecular dynamics simulations that resolve molecular-scale phenomena within polymer matrices inaccessible to screening methods. We highlight how these methods have informed experimental design for dissolution-based recycling and solvent-based contaminant removal. Finally, we discuss the prospective role of agentic AI in integrating these computational tools with real-time sorting data to adapt purification conditions to the compositional variability of real post-consumer feedstocks. This review charts a path toward computationally guided solution-based purification workflows that can respond to the complexity inherent in plastic waste streams.

Altamimi, Ali [Univ. of Wisconsin, Madison, WI (Un↗

High-throughput calculations of charged point defect properties with semi-local density functional theory—performance benchmarks for materials screening applications

Abstract Calculations of point defect energetics with Density Functional Theory (DFT) can provide valuable insight into several optoelectronic, thermodynamic, and kinetic properties. These calculations commonly use methods ranging from semi-local functionals with a-posteriori corrections to more computationally intensive hybrid functional approaches. For applications of DFT-based high-throughput computation for data-driven materials discovery, point defect properties are of interest, yet are currently excluded from available materials databases. This work presents a benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We consider three different a-posteriori correction sets implemented in an automated workflow, and evaluate the qualitative and quantitative differences among four different categories of defect information: thermodynamic transition levels, formation energies, Fermi levels, and dopability limits. We highlight qualitative information that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy.

36 MATERIALS SCIENCE↗

Rebuilding the Habitable Zone from the Bottom Up

The idea of habitability has been useful for categorizing places where life might exist in the universe. Circumstellar habitable zones estimate whether rocky planets might harbor liquid surface water, geophysical and tidal habitable zones gauge places like icy moons. Combined with biological constraints of energy and element availability, habitability can be key for eliminating false positive biosignatures. But is there a more fundamental and complete way to talk about all of this? I’ll share a proposal for a concept of “computational zones” that digs deep into the fundamentals of what we think we know about life, and its relationship to matter, energy, and thermodynamics. With some simple examples I’ll also show how we can apply this to real (and imagined) questions, and what it suggests for the kinds of measurements and missions that could have the greatest astrobiological impact.

icy moons↗

Statistical perspective on embrittling potency for intergranular fracture

Embrittling potency is a thermodynamic metric that assesses the influence of solute segregation to a grain boundary (GB) on intergranular fracture. Historically, authors of studies have reported embrittling potency as a single scalar value, assuming a single segregation site of importance at a GB and a particular cleavage plane. However, the topography of intergranular fracture surfaces is not generally known a priori. Accordingly, we, in this study, present a statistical ensemble approach to compute embrittling potency, where many free surface (FS) permutations are systematically considered to model fracture of a GB. The result is a statistical description of the thermodynamics of GB embrittlement. As a specific example, embrittling potency distributions are presented for Cr segregation to sites at two Ni $\langle 111 \rangle$ symmetric tilt GBs using atomistic simulations. We show that the average embrittling potency for a particular GB site, considering an ensemble of FS permutations, is not equal to the embrittling potency computed using the lowest energy pair of FSs. A mean GB embrittlement is proposed, considering both the likelihood of formation of a particular FS and the probability of solute occupancy at each GB site, to compare the relative embrittling behavior of two distinct GBs.

36 MATERIALS SCIENCE↗

Destabilizing high-capacity high entropy hydrides via earth abundant substitutions: From predictions to experimental validation

The vast chemical space of high entropy alloys (HEAs) makes trial-and-error experimental approaches for materials discovery intractable and often necessitates data-driven and/or first principles computational insights to successfully target materials with desired properties. In the context of materials discovery for hydrogen storage applications, a theoretical prediction-experimental validation approach can vastly accelerate the search for substitution strategies to destabilize high-capacity hydrides based on benchmark HEAs, e.g. TiVNbCr alloys. Here, in this study, machine learning predictions, corroborated by density functional theory calculations, predict substantial hydride destabilization with increasing substitution of earth-abundant Fe content in the (TiVNb) 75 Cr 25-x Fe x system. The as-prepared alloys crystallize in a single-phase bcc lattice for limited Fe content x < 7, while larger Fe content favors the formation of a secondary C14 Laves phase intermetallic. Short range order for alloys with x < 7 can be well described by a random distribution of atoms within the bcc lattice without lattice distortion. Hydrogen absorption experiments performed on selected alloys validate the predicted thermodynamic destabilization of the corresponding fcc hydrides and demonstrate promising lifecycle performance through reversible absorption/desorption. This demonstrates the potential of computationally expedited hydride discovery and points to further opportunities for optimizing bcc alloy ↔ fcc hydrides for practical hydrogen storage applications.

36 MATERIALS SCIENCE↗

A guided ion beam investigation of UO 2 + thermodynamics and f orbital participation: Reactions of U + + CO 2 , UO + + O 2 , and UO + + CO

A guided ion beam tandem mass spectrometer was employed to study the reactions of U + + CO 2 , UO + + O 2 , and the reverse of the former, UO + + CO. Reaction cross sections as a function of kinetic energy over about a three order of magnitude range were studied for all systems. The reaction of U + + CO 2 proceeds to form UO + + CO with an efficiency of 118% ± 24% as well as generating UO 2 + + C and UCO + + O. The reaction of UO + + O 2 forms UO 2 + in an exothermic, barrierless process and also results in the collision-induced dissociation of UO + to yield U + . In the UO + + CO reaction, the formation of UO 2 + in an endothermic process is the dominant reaction, but minor products of UCO + + O and U + + (O + CO) are also observed. Analysis of the kinetic energy dependences observed provides the bond energies, D 0 (U + –O) = 7.98 ± 0.22 and 8.05 ± 0.14 eV, D 0 (U + –CO) = 0.73 ± 0.13 eV, and D 0 (OU + –O) = 7.56 ± 0.12 eV. The values obtained for D 0 (U + –O) and D 0 (OU + –O) agree well with the previously reported literature values. To our knowledge, this is the first experimental measurement of D 0 (U + –CO). Furthermore, an analysis of the oxide bond energies shows that participation of 5f orbitals leads to a substantial increase in the thermodynamic stability of UO 2 + relative to ThO 2 + and especially transition metal dioxide cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stratus: An interactive steady state mixed layer model for personal computers

A steady-state, horizontally homogeneous, cloud-topped marine boundary layer model based primarily on the work of Lilly (1968) and Schubert et al., (1979) is presented. The conservative thermodynamic variables are equivalent potential temperature, theta(sub e), and total water mixing ratio, q + l. Some of the differences between this and Lilly's (1968) model are: radiation is allowed to penetrate into the boundary layer; cloud top values of longwave radiation, equivalent potential temperature, and water vapor mixing ratio are linear functions of height derived from climatological data at California coastal stations; and the closure assumption assumes a weighted average of Lilly's (1968) maximum and minimum entrainment theories. This model was programmed in FORTRAN and will run interactively on an IBM-compatible personal computer. The program allows the user to specify the geographical location, the wind speed, the sea-surface temperature, the large scale horizontal divergence, and the initial guess for cloud top height. Output includes the steady state values of cloud top and cloud base height, mixed layer equivalent potential temperature and total water mixing ratio, and the associated convective and radiative fluxes.

Guinn, Thomas A.↗

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗

Computational Thermochemistry of Jet Fuels and Rocket Propellants

The design of new high-energy density molecules as candidates for jet and rocket fuels is an important goal of modern chemical thermodynamics. The NASA Glenn Research Center is home to a database of thermodynamic data for over 2000 compounds related to this goal, in the form of least-squares fits of heat capacities, enthalpies, and entropies as functions of temperature over the range of 300 - 6000 K. The chemical equilibrium with applications (CEA) program written and maintained by researchers at NASA Glenn over the last fifty years, makes use of this database for modeling the performance of potential rocket propellants. During its long history, the NASA Glenn database has been developed based on experimental results and data published in the scientific literature such as the standard JANAF tables. The recent development of efficient computational techniques based on quantum chemical methods provides an alternative source of information for expansion of such databases. For example, it is now possible to model dissociation or combustion reactions of small molecules to high accuracy using techniques such as coupled cluster theory or density functional theory. Unfortunately, the current applicability of reliable computational models is limited to relatively small molecules containing only around a dozen (non-hydrogen) atoms. We propose to extend the applicability of coupled cluster theory- often referred to as the 'gold standard' of quantum chemical methods- to molecules containing 30-50 non-hydrogen atoms. The centerpiece of this work is the concept of local correlation, in which the description of the electron interactions- known as electron correlation effects- are reduced to only their most important localized components. Such an advance has the potential to greatly expand the current reach of computational thermochemistry and thus to have a significant impact on the theoretical study of jet and rocket propellants.

Crawford, T. Daniel↗

Thermal Extraction and Delivery from Advanced Reactor Systems

This document presents preliminary thermodynamic analysis of thermal extraction and delivery from advanced nuclear reactors. Light-water, high temperature gas, sodium fast, and molten salt reactors are templated and analyzed for thermal extraction from nominal steam production conditions. Uncontrolled, controlled, and partial admission steam extraction methods are described and analyzed to show the differences in system impact from the various methods. Thermal delivery modeling results show comparisons that have become possible through the computational capabilities of the thermal delivery model. Due to the significant degrees of freedom of the problem, specific cases are used to compare the thermodynamic costs of heat transport using various heat transfer fluids. All this work is done using Modelica models, allowing for the future application of models in dynamic analyses beyond the steady-state calculations presented in this report. Appendices describe additional capabilities and methods that can help users analyze complex problems using HYBRID models developed in Modelica.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Simple Criterion for Feasibility of Heat Integration between Distillation Streams Based on Relative Volatilities

In a multicomponent distillation configuration, there are numerous sources and sinks of heat, and a potential way to reduce the heat duty requirement is to perform heat integration. Unfortunately, an algorithmic search of the optimal heat integration opportunities is intractable when the required temperatures of intermediate streams are computed via complex models. Instead, in this work, we introduce pressure-scaled pseudo relative volatility, a new metric to compare stream temperatures. Here, we justify the use of pseudo relative volatility by proving that this variable is a monotonically increasing function of the liquid fraction in a saturated mixture stream. Using this metric, we derive a shortcut criterion to check the feasibility of various heat integration opportunities, such as thermal coupling via heat transfer (TCH). The advantage of this approach is that it circumvents the need for explicit temperatures and instead relies on composition, component relative volatilities, and pressure—quantities that are readily available in shortcut models for optimization of distillation configurations. Leveraging this fact, we propose a new optimization framework that identifies feasible TCHs that we consider within the formulation while minimizing the total heat duty of a distillation configuration. We demonstrate, on a few examples, that our formulation can identify heat duty efficient configurations, some of which are multieffect configurations. Using this methodology, we discover configurations that are not only simpler than the fully thermally coupled (FTC) configuration but also have a much lower heat duty.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗