Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

CeTi 2 O 6 —A Promising Oxide for Solar Thermochemical Hydrogen Production

A large entropy of reduction is crucial in achieving materials capable of high-efficiency solar thermochemical hydrogen (STCH) production through two-step thermochemical water splitting cycles. We have recently demonstrated that the onsite electronic entropy of reduction attains an extreme value of 4.26 kB at 1500 K in Ce 4+ → Ce 3+ redox reactions, which explains the high performance and uniqueness of CeO 2 as an archetypal STCH material. However, ceria requires high temperatures (T > 1500 °C) to achieve a reasonable reduction extent because of its large reduction enthalpy, which is a major obstacle in practical applications. Therefore, new materials with a large entropy of reduction and lower reduction enthalpy are required. Here, we perform a systematic screening to search for Ce 4+ -based oxides which possess thermodynamics superior to CeO 2 for STCH production. We first search the Inorganic Crystal Structure Database (ICSD) and literature for Ce 4+ -based oxides and subsequently use density functional theory to compute their reduction enthalpies (i.e., oxygen vacancy formation energies). Here, we find that CeTi 2 O 6 with the brannerite structure is the most promising candidate for STCH because it possesses three essential characteristics of an STCH material: (i) a smaller reduction enthalpy compared to ceria yet large enough to split water, (ii) a high thermal stability, as reported experimentally, and (iii) a large entropy of reduction associated with Ce 4+ → Ce 3+ redox. Our proposed design strategy suggests that further exploration of Ce 4+ oxides for STCH production is warranted.

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

Improved Precision Scaling for Simulating Coupled Quantum-Classical Dynamics

In this paper, we present a superpolynomial improvement in the precision scaling of quantum simulations for coupled quantum-classical systems. Such systems are found in, e.g., molecular-dynamics simulations within the Born-Oppenheimer approximation. By employing a framework based on the Koopman–von Neumann formulation of classical mechanics, we express the Liouville equation of motion as unitary dynamics and utilize phase kickback from a dynamical quantum simulation to calculate the quantum forces acting on classical particles. This approach allows us to simulate the dynamics of these classical particles without the overheads associated with measuring gradients and solving the equations of motion on a classical computer, resulting in a superpolynomial advantage at the price of increased space complexity. We demonstrate that these simulations can be performed in both microcanonical and canonical ensembles, enabling the estimation of thermodynamic properties from the prepared probability density. Published by the American Physical Society 2024

Simon, Sophia (ORCID:0000000332612338)↗
DOE OSTI2026

Wet-Bulb Temperature from Pressure, Relative Humidity, and Air Temperature

Abstract The algorithms currently in use for calculating the wet-bulb temperature from pressure, relative humidity, and air temperature have errors of a degree or more. Here, the equation for the thermodynamic wet-bulb temperature is derived using the Rankine–Kirchhoff approximations, which are highly accurate for meteorological and climatological applications. Likewise, equations are derived for the thermodynamic ice-bulb temperature and for the psychrometric wet-bulb and ice-bulb temperatures. The equations are fast to solve, requiring only a microsecond on a current laptop computer, and the results match empirical data to within hundredths of a degree. Furthermore, the equations reveal the existence of air temperatures and humidities that generate bistability: Evaporating liquid water is stable at a liquid temperature above freezing and, simultaneously, sublimating ice is stable at an ice temperature below freezing.

Romps, David M↗
NASA NTRS1975

Electrolytic hydrogen production: An analysis and review

The thermodynamics of water electrolysis cells is presented, followed by a review of current and future technology of commercial cells. The irreversibilities involved are analyzed and the resulting equations assembled into a computer simulation model of electrolysis cell efficiency. The model is tested by comparing predictions based on the model to actual commercial cell performance, and a parametric investigation of operating conditions is performed. Finally, the simulation model is applied to a study of electrolysis cell dynamics through consideration of an ideal pulsed electrolyzer.

Evangelista, J.↗
NASA NTRS1984

Calculation of nonequilibrium radiation in the flight regimes of aeroassisted orbital transfer vehicles

A computer code has been developed that calculates radiative properties of nonequilibrium air in the low-density regimes expected during the flight of aeroassisted, orbital transfer vehicles. From the given nonequilibrium thermodynamic state variables, the code calculates number densities of internal states and the accompanying emission and absorption characteristics. In addition, the code calculates the number density of the hypothetical gas in radiative equilibrium that produces the same radiation emission as the given nonequilibrium gas. Sample results are shown to demonstrate how the code predicts the nonlinear variation of radiation with density at low densities.

Park, C.↗
NASA NTRS1987

Computer-based phosphoric acid fuel cell analytical tools Descriptions and usages

Simulation models have been developed for the prediction of phosphoric acid fuel cell (PAFC) powerplant system performance under both transient and steady operation conditions, as well as for the design of component configurations and for optimal systems synthesis. These models, which are presently computer-implemented, are an engineering and a system model; the former being solved by the finite difference method to determine the balances and properties of different sections, and the latter using thermodynamic balances to set up algebraic equations that yield physical and chemical properties of the stream for one operating condition.

Lu, C.↗
NASA NTRS1992

Determining The Degree Of Graphitization In Carbon Composites

Method for determining degree of graphitization of carbon/carbon composite material combines infrared spectral experimental technique with thermodynamical theoretical model of solubility. By interrupting manufacturing process at almost any time and applying method, one can judge quality of composite and monitor chemical transformations taking place during process. Also computes effective solubility parameter for composite considered as unitary material. Parameter significant for subsequent fabrication steps involving chemical treatments; e.g., coating or joining with adhesives.

Lawson, Daniel D.↗
NASA NTRS1993

Rotary engine performance computer program (RCEMAP and RCEMAPPC): User's guide

This report is a user's guide for a computer code that simulates the performance of several rotary combustion engine configurations. It is intended to assist prospective users in getting started with RCEMAP and/or RCEMAPPC. RCEMAP (Rotary Combustion Engine performance MAP generating code) is the mainframe version, while RCEMAPPC is a simplified subset designed for the personal computer, or PC, environment. Both versions are based on an open, zero-dimensional combustion system model for the prediction of instantaneous pressures, temperature, chemical composition and other in-chamber thermodynamic properties. Both versions predict overall engine performance and thermal characteristics, including bmep, bsfc, exhaust gas temperature, average material temperatures, and turbocharger operating conditions. Required inputs include engine geometry, materials, constants for use in the combustion heat release model, and turbomachinery maps. Illustrative examples and sample input files for both versions are included.

Bartrand, Timothy A.↗
NASA NTRS1985

Calculation of Nonequilibrium Radiation in the Flight Regimes of Aeroassisted Orbital Transfer Vehicles

A computer code has been developed that calculates radiative properties of nonequilibrium air in the low-density regimes expected during the flight of aeroassisted, orbital transfer vehicles. From the given nonequilibrium thermodynamic state variables, the code calculates number densities of internal states and the accompanying emission and absorption characteristics. In addition, the code calculates the number density of the hypothetical gas in radiative equilibrium that produces the same radiation emission as the given nonequilibrium gas. Sample results are shown to demonstrate how the code predicts the non-linear variation of radiation with density at low densities.

Park, Chul↗
NASA NTRS2001

Cellular Spacing Selection During the Directional Solidification of Binary Alloys. A Numerical Approach

The evolution of cellular solid/liquid interfaces from an initially unstable planar front was studied by means of a two-dimensional computer simulation. The developed numerical model makes use of an interface tracking procedure and has the capability to describe the dynamics of the interface morphology based on local changes of the thermodynamic conditions. The fundamental physics of this formulation was validated against experimental microgravity results and the predictions of the analytical linear stability theory. The performed simulations revealed that in certain conditions, based on a competitive growth mechanism, an interface could become unstable to random perturbations of infinitesimal amplitude even at wavelengths smaller than the neutral wavelength, lambda(sub c), predicted by the linear stability theory. Furthermore, two main stages of spacing selection have been identified. In the first stage, at low perturbations amplitude, the selection mechanism is driven by the maximum growth rate of instabilities while in the second stage the selection is influenced by nonlinear phenomena caused by the interactions between the neighboring cells. Comparison of these predictions with other existing theories of pattern formation and experimental results will be discussed.

Catalina, Adrian V.↗
DOE OSTI2022

Determining the Role of Synthesis and Reaction Conditions on the Exposed Surface and Shape of ..beta..-Mo2C Catalysts

The production of drop-in quality biofuels typically can be achieved through catalytic fast pyrolysis (CFP) of biomass, followed by hydrodeoxygenation (HDO) of the vapor-phase product to reduce the oxygen content and stabilize the bio-oil product. To favor selectivity toward desired deoxygenation products, bifunctional catalysts with both acidic and hydrogenation sites are often employed. Due to both its oxophilicity and its ability to promote hydrogenation reactions, molybdenum carbide catalysts have shown promise for promoting the HDO of CFP vapors. However, the surface structure and composition of carbide catalyst particles are not well understood, precluding the establishment of fundamental structure-function relationships needed for the rational development of molybdedum carbide for HDO, as well as other reactions. For orthorhombic ..beta..-Mo2C, the (100) facet is largely studied for computational surface-reaction analyses, possibly due to its similarities with the close-packed (111) facet of noble-metal catalysts. However, it has previously been postulated through ab initio thermodynamic studies that the (111), (110), and (011) facets are most stable, though there is no consensus. Furthermore, while a majority of previous surface stability analyses of ..beta..-Mo2C focus on bulk terminations, the highly carburizing synthesis conditions and the use of ..beta..-Mo2C in carbon-containing reaction applications could introduce non-stoichiometric coverages of surface carbon. As such, this study uses ab initio thermodynamics calculations to assess the relative surface stability of the low-index facets of ..beta..-Mo2C, including consideration of varying the coverage of surface carbon. From this information, Wulff constructions are utilized to predict the synthesis and reaction condition-dependent shape and surface carbon composition of ..beta..-Mo2C catalyst particles.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
NASA NTRS2006

Using Satellite Data in Weather Forecasting: I

The GOES Product Generation System (GPGS) is a set of computer codes and scripts that enable the assimilation of real-time Geostationary Operational Environmental Satellite (GOES) data into regional-weather-forecasting mathematical models. The GPGS can be used to derive such geophysical parameters as land surface temperature, the amount of precipitable water, the degree of cloud cover, the surface albedo, and the amount of insolation from satellite measurements of radiant energy emitted by the Earth and its atmosphere. GPGS incorporates a priori information (initial guesses of thermodynamic parameters of the atmosphere) and radiometric measurements from the geostationary operational environmental satellites along with mathematical models of physical principles that govern the transfer of energy in the atmosphere. GPGS solves the radiative-transfer equation and provides the resulting data products in formats suitable for use by weather-forecasting computer programs. The data-assimilation capability afforded by GPGS offers the potential to improve local weather forecasts ranging from 3 hours to 2 days - especially with respect to temperature, humidity, cloud cover, and the probability of precipitation. The improvements afforded by GPGS could be of interest to news media, utility companies, and other organizations that utilize regional weather forecasts.

Jedlovec, Gary J.↗
DOE OSTI2022

Alchemical Free Energy Estimators and Molecular Dynamics Engines: Accuracy, Precision, and Reproducibility

The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat a myriad of diseases. In this work, we examine the computation of alchemical relative binding free energies with an eye for assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2, and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2, and NAMD3 alpha. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between them of 0.50 kcal/mol. The correlation between packages is very good, with the lowest Spearman’s, Pearson’s and Kendall’s tau correlation coefficients being 0.92, 0.91, and 0.76, respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.

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

A Model to Assess the Risk of Ice Accretion Due to Ice Crystal Ingestion in a Turbofan Engine and its Effects on Performance

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that were attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was one or more of the following anomalies: degraded engine performance, engine roll back, compressor surge and stall, and flameout of the combustor. The main focus of this research is the development of a computational tool that can estimate whether there is a risk of ice accretion by tracking key parameters through the compression system blade rows at all engine operating points within the flight trajectory. The tool has an engine system thermodynamic cycle code, coupled with a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Assumptions are made to predict the complex physics involved in engine icing. Specifically, the code does not directly estimate ice accretion and does not have models for particle breakup or erosion. Two key parameters have been suggested as conditions that must be met at the same location for ice accretion to occur: the local wet-bulb temperature to be near freezing or below and the local melt ratio must be above 10%. These parameters were deduced from analyzing laboratory icing test data and are the criteria used to predict the possibility of ice accretion within an engine including the specific blade row where it could occur. Once the possibility of accretion is determined from these parameters, the degree of blockage due to ice accretion on the local stator vane can be estimated from an empirical model of ice growth rate and time spent at that operating point in the flight trajectory. The computational tool can be used to assess specific turbine engines to their susceptibility to ice accretion in an ice crystal environment.

Jorgenson, Philip C. E.↗
DOE OSTI2024

Thermodynamic and Kinetic Activity Descriptors for the Catalytic Hydrogenation of Ketones

Activity descriptors are a powerful tool for the design of catalysts than can efficiently utilize H 2 with minimal energy losses. In this study, we develop the use of hydricity and H - self-exchange rates as thermodynamic and kinetic descriptors for the hydrogenation of ketones by molecular catalysts. Two complexes with known hydricity, HRh(dmpe) 2 and HCo(dmpe) 2 , were investigated for the catalytic hydrogenation of ketones under mild conditions (1.5 atm, 25 °C). The rhodium catalyst proved to be an efficient catalyst for a wide range of ketones, whereas the cobalt catalyst could only hydrogenate electron-deficient ketones. Using a combination of experiment and electronic structure theory, thermodynamic hydricity values were established for 46 alkoxide/ketone pairs in both MeCN and THF solvent. Through comparison of the hydricities of the catalysts and substrates, it was determined that catalysis was only observed for catalyst/ketone pairs with an exergonic H - transfer step. Mechanistic studies revealed that H - transfer was rate-limiting step for catalysis, allowing for the experimental and computation construction of linear free-energy relationships (LFERs) for H - transfer. Further analysis revealed the LFERs could be reproduced using Marcus theory, in which the H - self-exchange rates for the HRh/Rh + and ketone/alkoxide pairs were used to predict the experimentally measured catalytic barriers within 2 kcal mol -1 . Finally, these studies significantly expand the scope of catalytic reactions that can be analyzed with a thermodynamic hydricity descriptor and firmly establish Marcus theory as a valid approach to develop kinetic descriptors for designing catalysts for H - transfer reactions.

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

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗
DOE OSTI2022

Predictive coarse-grained (CG) modeling of morphologies in polymer nanocomposites with specific and directional intermolecular interactions (Final Report)

The overarching goal of the proposed work was to develop predictive models for investigating structure and dynamics in soft materials with chemistries that have specific and directional molecular interactions. The motivation behind studying materials with specific and directional interactions lies in the many desirable features these interactions provide when designing novel soft materials. Soft materials with specific and direction interactions (such as hydrogen bonds or H-bonds) can have a) thermally reversible phase behavior with different functions with varying temperature, b) precisely tuned nanostructure with desirable geometries that afford unique physical properties (e.g., color response, mechanical properties) and c) well-mixed/blended morphologies that are useful for variety of applications in energy field (e.g., materials for batteries require use of blended polymers where one domain gives superior mechanical properties and one domain promotes electrical conduction). Notably, biology makes extensive use of specific and directional interactions, in many cases based on H-bonds, to construct materials with precisely defined architectures and properties. Engineered soft materials with precisely tuned nanostructures and improved processiblity through thermoresponsive phase behavior are useful in numerous applications that are relevant to the Department of Energy (DOE) including high efficiency electronic devices, light-weight high-strength composite materials for batteries and fuel cells, and polymer membranes for separations, etc. While past computational studies have been tremendously useful in understanding molecular phenomena and guiding synthesis of new macromolecular soft materials for a wide variety of applications, the inability to capture small scale specific and directional interactions alongside macromolecular length and time scales represented a key limitation of most studies to date. Our work in this project addressed this grand challenge in computational materials chemistry, i.e., ability to model the anisotropic, directional, and specific interactions that govern the behavior of many macromolecular soft matter systems of interest, thus, has greatly expanded the predictive potential of simulations. Specifically the key outcomes were: successful development of new coarse-grained (CG) polymer models to study generic and specific polymer chemistries in which hydrogen-bonding interactions are dominant. These CG models were then used in molecular simulations to study structure and thermodynamics in polymer nanocomposites and blends; some studies were conducted in collaboration with experimentalists. We also published a perspective and a viewpoint which included some of the work we completed in this DOE project; we believe these perspective and viewpoint articles guide other researchers in the soft materials community on how to extend the computational approaches and models we have developed for the purposes of their studies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling and Optimization of a Rotating Packed Bed Contactor with a Tetraamine-Appended Metal−Organic Framework for CO 2 Capture

A potential contactor technology for sorbent-based CO 2 capture is the rotating packed bed that contains separate sections for continuous adsorption and desorption. A heat exchanger can be embedded to remove heat in the adsorption section and add heat in the desorption section. In this work, we develop a two-dimensional (2D) model of a rotating packed bed for use in CO 2 capture applications. Mass and energy balances for the model are developed based on a Ljungström-type air preheater, which accounts for the counter-current axial flow of gas phases in separate sections of the bed and the rotation of a solid sorbent, which cycles between adsorption and desorption sections. The sorbent used for this analysis is the tetraamine-appended metal−organic framework Mg 2 (dobpdc)(3−4− 3), chosen for its stability and affinity for CO 2 at low partial pressures, such as those from a natural gas power plant source. An optimization problem is solved that considers the trade-off between maximizing the productivity of the bed and minimizing energy consumption. Maximum productivity and minimum energy are found to be 8.53 kg/h/m 3 and 3.84 MJ/kg, respectively, when these objectives are optimized independently. It is observed that the flue gas pressure and bed rotational speed are the desired operating variables to vary for model-based design of experiments to reduce uncertainty in parameter estimation, as these two variables yielded the most information content based on the Fisher information matrix.

20 FOSSIL-FUELED POWER PLANTS↗