Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical equilibrium”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Zero carbon solid-state rechargeable redox fuel for long duration and seasonal storage

Here, this work presents a unique thermochemical process for charging magnesium-manganese-oxide-based solid-state rechargeable redox fuel. High-temperature heating of the processing furnace can be driven by either renewable electricity or concentrated solar power. The simple recyclable fuel charging concept is based on a tubular falling bed reactor with countercurrent oxygen-depleted gas flow for complete heat recuperation. The main focus of this work is achieving solid flowability at high temperatures (1,450°C) and extracting chemically charged solid at ambient temperature with minimal energy loss. The operation strategies described in this work have enabled consistent magnesium-manganese-oxide particle flow up to 1,450°C. The measured extent of the thermal reduction reaction after cooling the particles is more than 90% of the fully reduced state at equilibrium. The thermal-to-chemical efficiency and overall system efficiency are 96% and 35% respectively, which are the highest reported for thermochemical fuels to date.

25 ENERGY STORAGE↗
DOE OSTI2026

GEOSH: Ideal Gas Chemical Equation of State

We present the framework and methodology for the new Los Alamos National Laboratory (LANL) G as chemical E quation O f S tate at H igher temperatures code (GEOSH) which aims to accurately model the behavior of chemically complex gaseous mixtures in equilibrium. Assuming the ideal gas approximation, GEOSH leverages the recursive nature of the Saha ionization and molecular equations in order to eliminate the molecular and ionic degrees of freedom, thereby reducing the problem size to the number of atomic species plus one for the free electrons if ions are included. This approach allows the chemical species, both molecular and ionic, of the mixture to be expressed in terms of the abundances of the elemental species. As a result, the GEOSH framework achieves a reduction in computational expense, increased processing speed, and the capability to efficiently model large-scale chemical networks. This report provides the necessary physical background and theoretical foundations for the GEOSH code, accompanied by benchmarking studies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Thermodynamic equilibrium and kinetic fundamentals of oxide dissolution in aqueous solution

Dissolution of oxides in aqueous solutions is fundamentally important for a range of applications and a critical process that determines the chemical durability of industrial ceramics, the performance of nuclear waste forms, and the chemical weathering of minerals. The thermodynamic equilibrium and kinetics of dissolution reactions are key to determining the rate at which oxides dissolve. The increase in collaborative research across disciplines in materials research necessitates a common background to tackle shared scientific problems across different fields. This review selectively examines the fundamentals of dissolution theories that have been developed in chemistry, geochemistry, and materials science, and assembles them into a single collective document for the broader materials science community. Applications of the theories are highlighted using examples from specific areas, but can be similarly applied to other areas. Challenges and future research needs for a predictive-level understanding are discussed in light of the current literature.

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

The Transient Covalent Bond in Abiotic Nonequilibrium Systems

Abstract Biochemical systems accomplish many critical functions with by operating out‐of‐equilibrium using the energy of chemical fuels. The formation of a transient covalent bond is a simple but very effective tool in designing analogous reaction networks. This Minireview focuses on the fuel chemistries that have been used to generate transient bonds in recent demonstrations of abiotic nonequilibrium systems (i.e., systems that do not make use of biological components). Fuel reactions are divided into two fundamental classifications depending on whether the fuel contributes structural elements to the activated state, a distinction that dictates how they can be used. Reported systems are further categorized by overall fuel reaction (e.g., hydrolysis of alkylating agents, carbodiimide hydration) and illustrate how similar chemistry can be used to effect a wide range of nonequilibrium behavior, ranging from self‐assembly to the operation of molecular machines.

Kariyawasam, Lasith S.↗
DOE OSTI2020

The Transient Covalent Bond in Abiotic Nonequilibrium Systems

Biochemical systems accomplish many critical functions with by operating out-of-equilibrium using the energy of chemical fuels. The formation of a transient covalent bond is a simple but very effective tool in designing analogous reaction networks. This Minireview focuses on the fuel chemistries that have been used to generate transient bonds in recent demonstrations of abiotic nonequilibrium systems (i.e., systems that do not make use of biological components). Here, fuel reactions are divided into two fundamental classifications depending on whether the fuel contributes structural elements to the activated state, a distinction that dictates how they can be used. Reported systems are further categorized by overall fuel reaction (e.g., hydrolysis of alkylating agents, carbodiimide hydration) and illustrate how similar chemistry can be used to effect a wide range of nonequilibrium behavior, ranging from self-assembly to the operation of molecular machines.

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

A comparative study of solute trapping in Fe-(33–45 at%) Cu alloys manufactured by laser directed energy deposition

The high cooling rates in laser directed energy deposition (DED-LB) of alloys lead to substantial amounts of solute trapping as solute atoms cannot diffuse away from the solid/liquid interface before it advances. In some concentrated alloys, this results in supersaturated solid phases that form nanoscale hierarchical microstructures when the solute atoms precipitate out during reheating from subsequent laser passes. We choose the Iron-Copper (Fe-Cu) binary alloy as model system as it is chemically homogeneous in the liquid phase and has negligible solid solubility at room temperature. Two alloys with nominal compositions in atomic (at.) %, Fe67Cu33 and Fe55Cu45, were manufactured using DED-LB. Scanning transmission electron microscopy (STEM), energy dispersive spectroscopy (EDS) and wavelength dispersive spectroscopy (WDS) were used to characterize the nanostructures and heterogeneous chemical compositions. A non-equilibrium solute partitioning model was used to compute the supersaturated chemistries of the constituent phases and validated with experimentally measured compositions. The measured phase compositions of the two alloys were very similar, at roughly 12 at.% Cu and 4 at.% Fe in the α(bcc)-Fe and ε(fcc)-Cu phases respectively, despite having different processing parameters and mechanical behavior. This indicates that the total thermal history, that depends on both the processing parameters and laser scan pattern, plays a stronger role on the final microstructure evolution than just the initial quantity of trapped solute. In addition, we find that current non-equilibrium solute partitioning models applied on the continuum scale fall short of predicting accurate quantitative phase compositions in concentrated alloys, although the qualitative trends are captured correctly.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Acceleration of Thermochemistry Solves in MOOSE and Pronghorn

This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Technical note: A modified formulation of dynamic energy budget theory for faster computation of biological growth

Abstract. The mass conservation equation in the presence of boundary fluxes and chemical reactions from non-equilibrium thermodynamics is used to derive a modified dynamic energy budget (mDEB) model. Compared to the standard dynamic energy budget (sDEB) model (Kooijman, 2009), this modified formulation does not place the dilution effect in the mobilization kinetics of reserve biomass, and it maintains the partition principle for reserve mobilization dynamics for both linear and non-linear kinetics. Overall, the mDEB model shares most features with the sDEB model. However, for biological growth that requires multiple nutrients, the mDEB model is computationally much more efficient by not requiring numerical iterations for obtaining the specific growth rate. In an example of modeling the growth of Thalassiosira weissflogii in a nitrogen-limiting chemostat, the mDEB model was found to have almost the same accuracy as the sDEB model while requiring almost half of the computing time of the sDEB model. Since the sDEB model has been successfully applied in numerous studies, we believe that the mDEB model can help improve the modeling of biological growth and the associated ecosystem processes in various contexts.

Tang, Jinyun↗
DOE OSTI2022

Evaluation of thermodynamic closure models for partially reacted two-phase mixture of condensed phase explosives

One of the key fundamental issues that is crucial in the continuum modeling of reactive flow phenomena is the thermodynamically consistent description of reaction mixture properties. To define the mixture properties, thermodynamic closure rules that relate the properties of the individual reaction components to the mixture properties are required. In the context of reactive two-phase modeling approaches, various strategies to define the thermodynamic closures have been adopted such as pressure temperature (PT) equilibrium between the individual reaction components, pressure (specific) volume (PV) equilibrium, etc. The choice of closure rules determines the relative distribution of specific volume and energy across the reaction components that comprise the mixture. Therefore, depending on the choice of the closure, the mixture thermodynamic behavior can vary. The present work examines the effect of different closure approaches on the thermodynamic properties of the reaction mixture. The analysis is performed for a condensed phase HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) based plastic bonded explosive (PBX) 9501 explosive using four different thermodynamic closures, viz., PT equilibrium, PV equilibrium, volume temperature (VT) equilibrium, and pressure (P) equilibrium with reactants on an isentrope. Furthermore, the relative variations in the thermodynamic properties of the mixture are analyzed and compared under both compression and expansion loading regimes. It is shown that out of the four closure models, only PT equilibrium and P equilibrium closures lead to a thermodynamically accurate description of the mixture under both compression and expansion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dissipative Assembly of Macrocycles Comprising Multiple Transient Bonds

Abstract Dissipative assembly has great potential for the creation of new adaptive chemical systems. However, while molecular assembly at equilibrium is routinely used to prepare complex architectures from polyfunctional monomers, species formed out of equilibrium have, to this point, been structurally very simple. In most examples the fuel simply effects the formation of a single short‐lived covalent bond. Herein, we show that chemical fuels can assemble bifunctional components into macrocycles containing multiple transient bonds. Specifically, dicarboxylic acids give aqueous dianhydride macrocycles on treatment with a carbodiimide. The macrocycles are assembled efficiently as a consequence of both fuel‐dependent and fuel‐independent mechanisms; they undergo slower decomposition, building up as the fuel recycles the components, and are a favored product of the dynamic exchange of the anhydride bonds. These results create new possibilities for generating structurally sophisticated out‐of‐equilibrium species.

Hossain, Mohammad Mosharraf↗
DOE OSTI2026

Intense cavitation-assisted electric discharge as a promising tool for water treatment

This study investigates interrelations between one-electrode Cavitation-Assisted Electric Discharge (CAED), two-electrode CAED, and recently discovered Intense CAED (I-CAED). The one-electrode CAED is a self-triggered nanosecond discharge with pulse energy in the micro-Joule range, which can be generated even by a DC voltage. I-CAED consists of a non-equilibrium part within a low-pressure cavitating region and a micro-spark traversing a liquid film. We hypothesize that CAED propagates from the high-voltage electrode as an ionization wave through bubbles of saturated vapor. Subsequently, the streamer-like discharges in the bubbles may form a continuous plasma channel. Inside the cavitating region, the plasma is strongly non-equilibrium, providing an ideal environment for generating chemically unstable species such as hydrogen peroxide (H 2 O 2 ). Plasma of I-CAED spark is characterized by high electron density and near-thermal equilibrium, emitting a continuous ultraviolet spectrum. The combination of these different discharge parts makes I-CAED in water a highly effective tool for the Advanced Oxidation Process, particularly in water disinfection. Experimentally demonstrated Electric Energy per Order value for disinfection of E. coli-contaminated water is as low as 0.135 ± 0.035 kWh/m 3 /order. Estimates show that the implementation of “dry electrodes” configuration reduces the erosion rate of the electrode material by at least one order of magnitude. Spectral analysis reveals that the continuum emission generated by I-CAED in proximity to metal electrodes deviates from the spectra of discharges spatially decoupled from the electrodes. We assume that this spectral divergence is attributable to blackbody-like emission originating from metallic nanoparticles form during the electrode's erosion process.

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

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

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

Transferable Water Potentials Using Equivariant Neural Networks

Machine learning interatomic potentials (MLIPs) have emerged as a technique that promises quantum theory accuracy for reduced cost. It has been proposed [J. Chem. Phys. 2023, 158, 084111] that MLIPs trained on solely liquid water data cannot accurately transfer to the vapor–liquid equilibrium while recovering the many-body decomposition (MBD) analysis of gas-phase water clusters. This suggests that MLIPs do not directly learn the physically correct interactions of water molecules, limiting transferability. In this work, we show that MLIPs using equivariant architecture and trained on 3200 liquid water structures reproduces liquid-phase water properties (e.g., density within 0.003 g/cm 3 between 230 and 365 K), vapor–liquid equilibrium properties up to 550 K, the MBD analysis of gas-phase water cluster up to six-body interactions, and the relative energy and the vibrational density of states of ice phases. We show that potentials developed using equivariant MLIPs allow transferability for arbitrary phases of water that remain stable in nanosecond long simulations.

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

The Sonora Substellar Atmosphere Models. II. Cholla: A Grid of Cloud-free, Solar Metallicity Models in Chemical Disequilibrium for the JWST Era

Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope (JWST) era, high-resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the presence of disequilibrium chemistry in these atmospheres as a function of parameters such as their effective temperature and surface gravity will allow us to place better constraints on the physics governing these atmospheres. This paper is part of a series of works presenting the Sonora grid of atmosphere models. In this paper, we present a grid of cloud-free, solar metallicity atmospheres for brown dwarfs and wide-separation giant planets with key molecular species such as CH 4 , H 2 O, CO, and NH 3 in disequilibrium. Our grid covers atmospheres with T eff ∈ [500 K, 1300 K], log g ∈ [3.0, 5.5] (cgs) and an eddy diffusion parameter of logK zz = 2,4 and 7 (cgs). We study the effect of different parameters within the grid on the temperature and composition profiles of our atmospheres. We discuss their effect on the near-infrared colors of our model atmospheres and the detectability of CH 4 , H 2 O, CO, and NH 3 using the JWST. We compare our models against existing MKO and Spitzer observations of brown dwarfs and verify the importance of disequilibrium chemistry for T dwarf atmospheres. Finally, we discuss how our models can help constrain the vertical structure and chemical composition of these atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Dynamic Heterogeneity and Kovacs’ Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. Here, we used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water’s structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

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

Far From Equilibrium Hydrodynamics and the Beam Energy Scan

The existence of hydrodynamic attractors in rapidly expanding relativistic systems has shed light on the success of relativistic hydrodynamics in describing heavy-ion collisions at zero chemical potential. As the search for the QCD critical point continues, it is important to investigate how out of equilibrium effects influence the trajectories on the QCD phase diagram. In this proceedings, we study a Bjorken expanding hydrodynamic system based on DMNR equations of motion with initial out of equilibrium effects and finite chemical potential in a system with a critical point. We find that the initial conditions are not unique for a specific freeze-out point, but rather the system can evolve to the same final state freeze-out point with a wide range of initial baryon chemical potential, µ B . For the same initial energy density and baryon density, depending on how far out of equilibrium the system begins, the initial µ B can vary by Δµ B ~ 350 MeV. Furthermore, our results indicate that knowledge of the out-of-equilibrium effects in the initial state provide vital information that influences the search for the QCD critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A high‐order discontinuous Galerkin approach for physics‐based thermospheric modeling

Abstract The accurate prediction of aerodynamic drag on satellites orbiting in the upper atmosphere is critical to the operational success of modern space technologies, such as satellite‐based communication or navigation systems, which have become increasingly popular in the last few years due to the deployment of constellations of satellites in low‐Earth orbit. As a result, physics‐based models of the ionosphere and thermosphere have emerged as a necessary tool for the prediction of atmospheric outputs under highly variable space weather conditions. This paper proposes a high‐fidelity approach for physics‐based space weather modeling based on the solution of the Navier–Stokes equations using a high‐order discontinuous Galerkin method, combined with a matrix‐free strategy suitable for high‐performance computing on GPU architectures. The approach consists of a thermospheric model that describes a chemically frozen neutral atmosphere in nonhydrostatic equilibrium driven by the external excitation of the Sun. A novel set of variables is considered to treat the low densities present in the upper atmosphere and to accommodate the wide range of scales present in the problem. At the same time, and unlike most existing approaches, radial and angular directions are treated in a nonsegregated approach. The study presents a set of numerical examples that demonstrate the accuracy of the approximation and validate the current approach against observational data along a satellite orbit, including estimates of established empirical and physics‐based models of the ionosphere‐thermosphere system. Finally, a one‐dimensional radial derivation of the physics‐based model is presented and utilized for conducting a parametric study of the main thermal quantities under various solar conditions.

Engineering↗
DOE OSTI2022

On the hot surface ignition of a wall-stagnating spray flame

The ignition of fuel sprays due to interaction with hot surfaces is an important phenomenon in the safety analysis of many engineering systems. We perform a parametric study of the hot surface ignition (HSI) of a fuel spray approaching a heated surface caused by the accidental leakage of a fuel line. Here, to this end, we employ a one-dimensional Eulerian-Eulerian formulation with a non-equilibrium evaporation model and a realistic chemical mechanism to describe $n$-dodecane fuel chemistry. We first describe and analyze the phenomenology of the unsteady processes leading to ignition using non-dimensionalized quantities. Through consideration of the temporal development of the most reactive mixture, we demonstrate that ignition occurs at a fuel-lean composition in a premixed region near the hot surface. Using non-dimensional parameters identified from the governing equations, we perform a parametric study of the time, location and local mixture composition at ignition and determine the ignition limits. We then identify the most important parametric sensitivities for physical analysis using a data-driven classification method. Our analysis demonstrates a contraction of the ignition limits with increased Stokes number and a regime of parametric insensitivity of igniting mixture composition. We also show that at high Damköhler numbers, the ignition location conforms to the parametric behavior of the thermal boundary layer, whereas at low Damköhler numbers approaching the ignition limit it reaches a near-unity value of the quenching Peclet number. We then compare the demonstrated parametric dependencies to the results of the quasi-steady asymptotic ignition literature, showing that our results are consistent with those obtained analytically within the limitations imposed by the simplified formulation of the latter.

42 ENGINEERING↗