SEARCH · Search NASA

Results for “predictive thermodynamics”

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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Predicting phase behavior in high entropy and chemically complex alloys

The interest in high entropy alloys and other metallic compounds with four or more elements at near-equiatomic ratios has drawn attention to the ability to rapidly predict phase behavior of these complex materials, particularly where existing thermodynamic data are lacking. In this paper, we discuss aspects of this from the point of view of predicting without utilizing (or fitting) experimental data. Of particular interest are heuristic approaches that provide prediction of single-phase compositions, more rigorous approaches that tackle the thermodynamics from a more fundamental point of view, and simulation approaches that provide further insight into the behaviors. This paper covers cases of all three of these, in order to examine the strengths and weaknesses of each approach, and to indicate directions where these may be utilized and improved upon. Of particular interest is moving beyond “which composition may form a solid solution,” to recognizing the importance of underlying thermodynamic realities that affect the temperature- and composition-dependent transformations of these materials.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Degenerate parametric down-conversion facilitated by exciton-plasmon polariton states in a nonlinear plasmonic cavity

Abstract We study the effect of degenerate parametric down-conversion (DPDC) in an ensemble of two-level quantum emitters (QEs) coupled via near-field interactions to a single surface plasmon (SP) mode of a nonlinear plasmonic cavity. For this purpose, we develop a quantum driven-dissipative model capturing non-equilibrium dynamics of the system in which incoherently pumped QEs have transition frequency tuned near the second-harmonic response of the SPs. Considering the strong coupling regime, i.e. the SP-QE interaction rate exceeds system dissipation rates, we find a critical SP-QE coupling attributed to the phase transition between normal and lasing steady states. Examining fluctuations above the system’s steady states, we predict new elementary excitations, namely, the exciton-plasmon polaritons formed by the two-SP quanta and single-exciton states of QEs. The contribution of two-SP quanta results in the linear scaling of the SP-QE interaction rate with the number of QEs,  o , as opposed to the  o -scaling known for the Dicke and Tavis–Cummings models. We further examine how SP-QE interaction scaling affects the polariton dispersions and power spectra in the vicinity of the critical coupling. For this purpose, we compare the calculation results assuming a finite ensemble of QEs and the model thermodynamic limit. The calculated power spectra predict an interplay of coherent photon emission by QEs near the second-harmonic frequency and correlated photon-pair emission at the fundamental frequency by the SPs (i.e. the photonic DPDC effect).

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2022

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/ cm 3 and 0 to 10 9 K) equation-of-state (EOS) table for a CH 1.72 O 0.37 N 0.086 quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. Here, the simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Simulation of a rapid compression machine for evaluation of ignition chemistry and soot formation using gasoline/ethanol blends

Due to the projected decline of demand for gasoline in light duty engines and the advent of ethanol as a green fuel, the use of gasoline/ethanol blend fuels in heavy duty applications are being investigated as they are projected to have lower cost and lower lifecycle green house gas (GHG) emissions. In heavy duty engines, the primary mode of combustion is mixing controlled combustion where wide range of mixture conditions (equivalence ratio) exist. Soot emissions of these fuels in richer conditions are not well understood. The goal of this research is to evaluate some commercially available soot modeling codes for the particulate matter emissions from gasoline/ethanol fuel blends, especially at fuel rich conditions. A Rapid Compression Machine (RCM) is modeled in a three-dimensional numerical simulation using CONVERGE computational software using a reduced chemical kinetic mechanism with SAGE chemistry solver and a RANS k-ϵturbulence model with a sector model including the creviced piston. The creviced piston is used in the experimental setup to reduce boundary layer effects and to maintain a homogeneous core in the reaction cylinder. Computational fluid dynamics simulations are conducted for different gasoline-ethanol fuel blends from E10 (10% ethanol v/v) to E100. The fuel blend is modeled as a surrogate mixture of toluene, iso-octane, n-heptane for gasoline content, and ethanol. The computational results were validated against experimental results using pressure measurements and laser extinction diagnostics. Different soot models are investigated to evaluate their capability of predicting the sooting tendencies of fuel blends, especially in richer conditions experienced during mixing-controlled combustion. The experimental combustion characteristics such as the ignition delay of different blends of fuel are reasonably well predicted. The Particulate Size Mimic (PSM) model accurately predicts the soot generation characteristics of the different fuels, but the Hiroyasu-NSC model falls short in this regard. For accurate prediction of soot with the PSM model, the thermodynamic conditions during combustion must be accurately modeled. While the current computational modeling tools can produce accurate results for the prediction of particulate matter emissions, there is much work to be done in improving our understanding of the underlying fundamental processes.

Energy & Fuels↗

DOE OSTI2023

Evaluation of PBL Parameterization Schemes in WRF Model Predictions during the Dry Season of the Central Amazon Basin

Planetary Boundary Layer (PBL) parameterization schemes are employed to handle subgrid-scale processes on atmospheric models, playing a key role in accurately representing the atmosphere. Recent studies have shown that PBL schemes are particularly fundamental to the depiction of PBL height (PBLH), especially over the Amazon. In the present study, we investigated the performance of PBL schemes on the representation of meteorological variables, turbulent fluxes, PBL vertical structures, and PBLH over the central Amazon basin under dry conditions, taking advantage of observations from the Observations and Modeling of the Green Ocean Amazon campaign (GoAmazon2014/5) for validation and evaluation. Numerical experiments were carried out within the WRF model using eight PBL schemes for two dry periods from 2014 (typical year) and 2015 (El-Niño year), and results from the 1-km resolution domain were directly compared to hourly in situ observations. In general, all PBL schemes present good performance to reproduce meteorological variables, with nonlocal (local) PBL schemes producing better performance in the 2014 (2015) study period. All PBL schemes in general overestimate (>100%) daytime turbulent fluxes. Thermodynamic (daytime) vertical structures are better predicted than mechanical (nocturnal) ones. The local MYNN2.5 scheme showed the overall best performance for PBLH prediction, mainly at night.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Implications of reduced-complexity aerosol thermodynamics on organic aerosol mass concentration and composition over North America

Atmospheric organic aerosol (OA) mass concentrations can be affected by water uptake through its impact on the gas–particle partitioning of semi-volatile compounds. Current chemical transport models (CTMs) neglect this process. We have implemented the Binary Activity Thermodynamics model coupled to a volatility basis set partitioning scheme in the GEOS-Chem CTM, providing an efficient reduced-complexity OA model that predicts relative-humidity-dependent mixing and partitioning thermodynamics, while limiting the impact on computational efficiency. We provide a quantitative assessment of this water-sensitive OA treatment, focusing on a subdomain over North America. The updated OA scheme predicts a spatiotemporal mean enhancement in surface-level OA mass concentration of 145 % for January 2019 and 76 % for July 2019 compared to GEOS-Chem's most advanced OA scheme. The temporal mean surface-level OA organic mass concentration can increase by up to ∼590 % for January 2019 and ∼280 % for July 2019, with the greatest enhancements occurring over the ocean. The updated OA scheme also quantifies the OA-associated water content. The simulations show how different OA precursors and related OA surrogates contribute and respond to water uptake, including that due to changes in temperature and relative humidity over the diurnal cycle in selected winter and summer months. These results are independent of future CTM improvements involving updates to chemical reaction schemes and emission inventories. Our water-sensitive OA scheme allows for a better representation of the seasonal and regional variations in OA mass concentration in CTMs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Additively manufactured Ni-20Cr to V functionally graded material: Computational predictions and experimental verification of phase formations

Throughout this work, a database for the Cr-Ni-V system was constructed by modeling the binary Cr-V and ternary Cr-Ni-V systems using the CALPHAD approach aided by density functional theory (DFT)-based first-principles calculations and ab initio molecular dynamics (AIMD) simulations. To validate this new database, a functionally graded material (FGM) using Ni-20Cr and elemental V was fabricated using directed energy deposition additive manufacturing (DED AM) and experimentally characterized. The deposited Ni-20Cr was pure fcc phase, while increasing the amount of V across the gradient resulted in the formation of sigma phase, followed by the bcc phase. The experimentally measured phase data was compared with computational predictions made using a Cr-Ni-V thermodynamic database from the literature as well as the database developed in the present work. The newly developed database was shown to better predict the experimentally observed phases due to its accurate modeling of binary systems within the database and the ternary liquid phase, which is critical for accurate Scheil calculations.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Deep energy-pressure regression for a thermodynamically consistent EOS model

Abstract In this paper, we aim to explore novel machine learning (ML) techniques to facilitate and accelerate the construction of universal equation-Of-State (EOS) models with a high accuracy while ensuring important thermodynamic consistency. When applying ML to fit a universal EOS model, there are two key requirements: (1) a high prediction accuracy to ensure precise estimation of relevant physics properties and (2) physical interpretability to support important physics-related downstream applications. We first identify a set of fundamental challenges from the accuracy perspective, including an extremely wide range of input/output space and highly sparse training data. We demonstrate that while a neural network (NN) model may fit the EOS data well, the black-box nature makes it difficult to provide physically interpretable results, leading to weak accountability of prediction results outside the training range and lack of guarantee to meet important thermodynamic consistency constraints. To this end, we propose a principled deep regression model that can be trained following a meta-learning style to predict the desired quantities with a high accuracy using scarce training data. We further introduce a uniquely designed kernel-based regularizer for accurate uncertainty quantification. An ensemble technique is leveraged to battle model overfitting with improved prediction stability. Auto-differentiation is conducted to verify that necessary thermodynamic consistency conditions are maintained. Our evaluation results show an excellent fit of the EOS table and the predicted values are ready to use for important physics-related tasks.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

A First–Principles–Based Sub–Lattice Formalism for Predicting Off–Stoichiometry in Materials for Solar Thermochemical Applications: The Example of Ceria

Theoretical models that reliably can predict off-stoichiometry in materials via accurate descriptions of underlying thermodynamics are crucial for energy applications. For example, transition-metal and rare-earth oxides that can tolerate a large number of oxygen vacancies, such as CeO 2 and doped CeO 2 , can split water and carbon dioxide via a two-step, oxide-based solar thermochemical (STC) cycle. The search for new STC materials with a performance superior to that of state-of-the-art CeO 2 can benefit from predictions accurately describing the thermodynamics of oxygen vacancies. The sub-lattice formalism, a common tool used to fit experimental data and build temperature-composition phase diagrams, can be useful in this context. Here, sub-lattice models are derived solely from zero-temperature quantum mechanics calculations to estimate fairly accurate temperature- and oxygen-partial-pressure-dependent off-stoichiometries in CeO 2 and Zr-doped CeO 2 . Physical motivations for deriving some of the “excess” sub-lattice model parameters directly from quantum mechanical calculations, instead of fitting to minimize deviations from experimental and/or theoretical data, are identified. As a result, important limitations and approximations of the approach used are specified and extensions to multi-cation oxides are also suggested to help identify novel candidates for water and carbon dioxide splitting and related applications.

08 HYDROGEN↗

DOE OSTI2020

Beyond Idealized Models of Nanoscale Metal Hydrides for Hydrogen Storage

Metal hydrides are attractive for compact, low-pressure hydrogen storage, yet a foundational understanding of factors governing their thermodynamics and kinetics is still lacking. Predictive modeling from the atomic to the microstructural scale plays a critical role in addressing these gaps, particularly for nanoscale materials, which promise improved performance but are difficult to probe. In this paper, we summarize strategies being developed within the Hydrogen Materials—Advanced Research Consortium (HyMARC) for going beyond conventional models to incorporate more complex physics, more realistic structures, and better approximation of operation conditions in simulations of nanoscale metal hydrides. We highlight four beyond-ideal factors that influence predicted performance: (1) surface anharmonic dynamics, (2) interface and surface energy penalties, (3) mechanical stress under confinement, and (4) the presence of native surface oxide. Approaches for addressing these factors are demonstrated on model materials representative of high-capacity hydrogen storage systems, and implications for understanding performance under operating conditions are discussed.

08 HYDROGEN↗

DOE OSTI2021

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Exchange-correlation functional challenges in modeling quaternary chalcogenides

The development of next-generation quaternary chalcogenides, such as Cu2ZnSnS4 (CZTS) and Cu2ZnGeS4 (CZGS), for solar energy and thermoelectric applications hinges upon both careful experimentation and accurate quantum mechanical modeling. To address the latter, many have turned to density functional theory (DFT), which offers several choices for the approximate treatment of electron exchange and correlation (XC). Popular XC functionals include the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) and the recently developed strongly constrained and appropriately normed (SCAN) meta-GGA. Extensions of DFT functionals, such as adding a Hubbard U correction and introducing a fraction of the Fock exchange (hybrid functionals), have been used widely to model systems containing 3d metal ions. However, no studies yet have compared comprehensively PBE(+U) and SCAN(+U) in the quality of their predictions of the bulk and defect thermodynamics of quaternary chalcogenides, which play a critical role in device fabrication and performance. Hence, here we calculate the (i) 0 K formation energies of bulk Ge compounds and (ii) neutral defect formation energies including charge-balanced (e.g., CuZn + ZnCu) and charge-imbalanced (e.g., CuSn) combinations of antisites and vacancies in CZTS and CZGS using the PBE, PBE +U, SCAN, SCAN +U, and the hybrid Heyd-Scuseria-Ernzerhof XC frameworks. We find that the formation energies of charge-imbalanced defects are more sensitive to the choice of the XC functional than those of charge-balanced defects, which can be explained by the differences in the extent of penalization of defect-generated delocalized electrons/holes by PBE, PBE +U, SCAN, and SCAN +U. Additionally, our results show that SCAN systematically underbinds Ge-containing compounds, thus highlighting the need for even further improvement of XC functionals. Based on our findings, we recommend the use of SCAN for modeling quaternary chalcogenides because its errors are systematic, and it has the firmest theoretical underpinning. Our work provides guidance for future modeling of quaternary chalcogenides.

14 SOLAR ENERGY↗

DOE OSTI2024

Selective formation of metastable polymorphs in solid-state synthesis

Metastable polymorphs often result from the interplay between thermodynamics and kinetics. Despite advances in predictive synthesis for solution-based techniques, there remains a lack of methods to design solid-state reactions targeting metastable materials. Here, we introduce a theoretical framework to predict and control polymorph selectivity in solid-state reactions. This framework presents reaction energy as a rarely used handle for polymorph selection, which influences the role of surface energy in promoting the nucleation of metastable phases. Through in situ characterization and density functional theory calculations on two distinct synthesis pathways targeting LiTiOPO 4 , we demonstrate how precursor selection and its effect on reaction energy can effectively be used to control which polymorph is obtained from solid-state synthesis. A general approach is outlined to quantify the conditions under which metastable polymorphs are experimentally accessible. With comparison to historical data, this approach suggests that using appropriate precursors could enable targeted materials synthesis across diverse chemistries through selective polymorph nucleation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Simulation and Sensitivity Analysis of a Cryogenic Carbon Capture System for the Cement Industry

The talk will describe modeling and analysis of cryogenic systems for carbon capture from industrial point sources, thermodynamic and process simulation requirements to predict capture from typical flue gases, quantification of temperature and flow profiles for a range of gas feed compositions and capture levels, and impacts of process conditions and heat transfer on energy requirements for heat integration.

Paul, Brandon↗

DOE OSTI2023

Surface Functional Groups Affect Iron (Hydr)oxide Heterogeneous Nucleation: Implications for Membrane Scaling

Because of its favorable thermodynamics and fast kinetics, heterogeneous solid nucleation on membranes triggers early-stage mineral scaling. Iron (hydr)oxide, a typical membrane scale, initially forms as nanoparticles that interact with surface functional groups on membranes, but these nanoscale phenomena are difficult to observe in real time. In this study, we utilized in situ grazing incidence small angle X-ray scattering and ex situ atomic force microscopy to examine the heterogeneous nucleation of iron (hydr)oxide on surface functional groups commonly used in membranes, including hydroxyl (OH), carboxyl (COOH), and fluoro (F) groups. We found that, compared to nucleation on hydrophilic OH- and COOH-surfaces, the high hydrophobicity of an F-modified surface significantly reduced the extents of both heterogeneously and homogeneously formed iron (hydr)oxide nucleation. Moreover, on the OH-surface, the high functional group density of 0.76 nmol/cm 2 caused faster heterogeneous nucleation than that on a COOH-surface, with a density of 0.28 ± 0.04 nmol/cm 2 . The F-surface also had the highest heterogeneous nucleation energy barrier (26 ± 0.6 kJ/mol), followed by COOH- (23 ±0.8 kJ/mol) and OH- (20 ± 0.9 kJ/mol) surfaces. The kinetic and thermodynamic information provided here will help us better predict the rates and extents of early-stage scaling of iron (hydr)oxide nanoparticles in membrane processes. Finally, this work provides both kinetic and thermodynamic information about iron (hydr)oxide nucleation controlled by membrane-related functional groups (OH, COOH, and F).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Thermodynamic modeling of the Nb-Ni system with uncertainty quantification using PyCalphad and ESPEI

Here, the Nb–Ni system is remodeled with uncertainty quantification (UQ) using software tools of PyCalphad and ESPEI (the Extensible, Self-optimizing Phase Equilibria Infrastructure) with the presently implemented capability of modeling site fraction based on Wyckoff positions. The five- and three-sublattice models are used to model the topologically close pack (TCP) μ-Nb 7 Ni 6 and δ-NbNi 3 phases according to their Wyckoff positions. The inputs for CALPHAD-based thermodynamic modeling include the thermochemical data as a function of temperature predicted by first-principles and phonon calculations based on density functional theory (DFT), ab initio molecular dynamics (AIMD) simulations, together with phase equilibrium and site fraction data in the literature. In addition to phase diagram and thermodynamic properties, the CALPHAD-based predictions of site fractions of Nb in μ-Nb 7 Ni 6 agree well with experimental data. Furthermore, the UQ estimation using the Markov Chain Monte Carlo (MCMC) method as implemented in ESPEI is applied to study the uncertainty of site fraction in μ-Nb 7 Ni 6 and enthalpy of mixing (ΔH mix ) in liquid.

36 MATERIALS SCIENCE↗

DOE OSTI2026

A new interatomic potential for mixed Mg-Al-Ga-In spinels

While density functional theory (DFT) has become the de facto approach for accurate simulation of materials at the atomic scale, there are many aspects of materials that are simply out of reach of DFT methods. In particular, finite temperature properties such as diffusivities, the structure and properties of grain boundaries and interfaces, and the study of defect properties in complex alloys are computationally challenging for DFT methods. Recently, a new class of spinels in which three cations order over two sublattices was discovered. In order to predict the properties of these types of structures, classical potentials are a must. Here, in this work, we derive a new classical potential for Mg-bearing spinels in which the B cations are Al, Ga, and/or In. The potential does well in describing the DFT energetics of various spinel structures as a function of chemistry and inversion. In particular, it reproduces the thermodynamically favorable MgAlGaO 4 structure while correctly predicting that neither MgAlInO 4 nor MgGaInO 4 are stable. Further, it reproduces physical trends in elastic properties as compared against experiment.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Microstructure-Dependent Rate Theory Model of Radiation-Induced Segregation in Binary Alloys

Conventional rate theory often uses the mean field concept to describe the effect of inhomogeneous microstructures on the evolution of radiation induced defect and solute/fission product segregation. However, the spatial and temporal evolution of defects and solutes determines the formation and spatial distribution of radiation-induced second phase such as precipitates and gas bubbles/voids, especially in materials with complicated microstructures and subject to high dose radiation. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed to investigate the effect of inhomogeneous thermodynamic and kinetics properties of defects on diffusion and accumulations of solute A in AB binary alloys. Four independent concentrations: atom A, interstitial A, interstitial B, and vacancy on [A, B] sublattice are used as field variables to describe temporal and spatial distribution and evolution of defects and solute A. The independent concentrations of interstitial A and interstitial B allow to describe their different generation rates, thermodynamic and kinetic properties, and release the assumptions of interstitial generation and sink strength used in the conventional rate theory. Microstructure and concentration dependent chemical potentials of defects are used to calculate the driving forces of defect diffusions. With the model, the effects of defect chemical potentials and mobilities on the RIS in polycrystalline AB model alloys have been simulated. The results demonstrated the model capability in predicting defect evolution in materials with inhomogeneous thermodynamic and kinetic properties of defects. The model can be extended to materials with complicated microstructures such as a wide range of grain size distribution, coating structure and multiphases as well as radiation-induced precipitation subject to severe radiation damage.

36 MATERIALS SCIENCE↗