Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-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.

341 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling↗
DOE OSTI2026

Asymmetric pathways for lithium extraction and recovery based on the two-phase equilibrium of layered oxides

Electrochemical intercalation offers a promising platform for Li + extraction. However, only limited types of electrode materials have been investigated. The challenge to broaden and tailor materials for electrochemical intercalation-based Li + extraction lies in the lack of understanding of material’s response upon co-intercalation of multiple ions, therefore, paired process design to enable reversible Li + extraction and recovery. Here, we showcase the design of asymmetric ion pathways for Li + extraction and recovery for host material with complex Li + and Na + interaction using layered cobalt oxide as a model material. The two-phase equilibrium of Na 0.48 CoO 2 and Li 0.94 CoO 2 governs Li + selectivity when a high depth of intercalation is achieved (low vacancy level). We show that the relative rate between ion exchange and intercalation is critical to determine the ion pathways. The relationship can be quantitatively compared using the average pseudo ion exchange rate (C pseudoIX ) and the intercalation rate (C inter ). The ion pathways at the three regimes with C pseudoIX > C inter , C pseudoIX ~ C inter , and C pseudoIX < C inter are constructed. By selecting the optimized ion pathway and particle size, we demonstrate 9.7×10 4 Li + selectivity with 99% purity Li + recovery from an initial 1:1000 Li: Na molar ratio solution using 115 mAh/g specific capacity.

electrochemistry↗
DOE OSTI2023

Dynamic Oscillation and Motion of Oil‐in‐Water Emulsion Droplets

The interplay of interfacial tensions on droplets results in a range of self‐powered motions that mimic those of living systems and serve as a tunable model to understand their complex non‐equilibrium behavior. Spontaneous shape deformations and oscillations are crucial features observed in nature but difficult to incorporate in synthetic artificial systems. Here, we report sessile oil‐in‐water emulsions that exhibit rapid oscillating behavior. The oscillations depend on the nature and concentration of the surfactant, the chemical composition of the oil, and the wettability of the solid substrate. The rapid changes in the contact angle per oscillation are observed using side‐view optical microscopy. We propose that the changes in the interfacial tension of the oil droplets is due to the partitioning of the surfactant into the oil phase and the movement of self‐emulsified oil out of the parent droplets giving rise to the rhythmic variation in droplet contact‐line. The ability to control and understand droplet oscillation can help model similar oscillations in out‐of‐equilibrium systems in nature and reproduce biomimetic behavior in artificial systems for various applications, such as microfluidic lab‐on‐a‐chip and adaptive materials.

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

RTE: A Computer Code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Mohammad H N Naraghi↗
DOE OSTI2026

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗
DOE OSTI2026

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Keldysh tuning of photoluminescence in a lead halide perovskite crystal

In 1964, Keldysh laid the groundwork for strong-field physics in atomic, molecular, and solid-state systems by delineating a ubiquitous transition from multiphoton absorption to quantum electron tunneling under intense AC driving forces. While both processes in semiconductors can generate carriers and result in photon emission through electron-hole recombination, the low quantum yields in most materials have hindered direct observation of the Keldysh crossover. Leveraging the large quantum yields of photoluminescence in lead halide perovskites, we show that we can not only induce bright light emission from extreme sub-bandgap light excitation but also distinguish between photon-induced and electric-field-induced processes. Our results are rationalized by the Landau-Dykhne formalism, providing insights into the non-equilibrium dynamics of strong-field light-matter interactions. These findings open new avenues for light upconversion and sub-bandgap photon detection, highlighting the potential of lead halide perovskites in advanced optoelectronic applications.

FOS: Physical sciences↗
DOE OSTI2026

Field-tailoring quantum materials via magneto-synthesis: metastable metallic and magnetically suppressed phases in a trimer iridate

We demonstrate that applying modest magnetic fields (< 0.1 T) during high-temperature crystal growth can profoundly alter the structure and ground state of a spin-orbit-coupled, antiferromagnetic trimer lattice. Using BaIrO₃ as a model system, whose ground state is intricately dictated by the trimer lattice, we show that magneto-synthesis , a field-assisted synthesis approach, stabilizes a structurally compressed, metastable metallic and magnetically suppressed phases inaccessible via conventional methods. These effects include a 0.85% reduction in unit cell, 4-order-of-magnitude decrease in resistivity, a 10-fold enhancement of the Sommerfeld coefficient, and the collapse of long-range magnetic order -- all intrinsic and bulk in origin. First-principles calculations confirm that the field-stabilized structure lies substantially above the ground state in energy, highlighting its metastable character. These large, coherent and correlated changes across multiple bulk properties, unlike those caused by dilute impurities, defects or off-stoichiometry, point to an intrinsic field-induced mechanism. The findings establish magneto-synthesis as a powerful new pathway for accessing non-equilibrium quantum phases in strongly correlated materials.

magneto-synthesis↗
NASA NTRS2026

Optimization of a 10–Step Reduced CH₄/O₂ Combustion Mechanism for RDRE Applications

A 10-step, 9-species reduced mechanism for methane/oxygen combustion is evaluated and optimized for use in rotating detonation rocket engine (RDRE) simulations at an operating pressure of 10 bar. The mechanism is assessed against GRI-Mech 3.0 for ignition delay times, adiabatic equilibrium temperature and species, Chapman-Jouguet (CJ) detonation parameters, and laminar flame speed. A sequential Nelder–Mead optimization procedure is applied to the rate constants of reactions R1, R4, and R5, targeting ignition delay agreement over the equivalence ratio range φ = 1.0–1.4 and temperature range T = 1400–2000 K. The optimized mechanism achieves ignition delay agreement within 28% of GRI-Mech 3.0 for the temperature range relevant to RDRE operation (τ ≥ 5 μs). Adiabatic equilibrium temperatures and CJ detonation parameters are reproduced within 0.15% of GRI-Mech 3.0 values, confirming thermodynamic consistency of the mechanism. Laminar flame speed predictions deviate significantly from GRI-Mech 3.0, which is attributed to the absence of HO₂ chemistry in the reduced mechanism; this limitation does not affect RDRE simulation fidelity since laminar flame propagation is of secondary importance in RDREs.

Computational Fluid Dynamics↗
DOE OSTI2026

In situ studies on heavy ion irradiation and partial oxidation induced stacking faults and nanograins in tungsten nanolaminates

Understanding the microstructural evolution of tungsten (W) under extreme irradiation environments is critical for its application as a plasma-facing material in future fusion reactors. Heavy ion irradiation study offers expedited damage accumulation and thus allows irradiation tolerance property prediction in a short period of time. In this study we explore in situ Kr ion irradiation of W nanolaminates at 800 °C and uncover a complex interplay of irradiation-induced transformations, including the emergence of stacking faults, nanograin formation, and generation of thickened grain boundaries defined as GB regions with the thickness of up to 15 nm. High-resolution transmission electron microscopy studies reveal the formation of extended planar faults and a metastable hexagonal close-packed (hcp) phase within the body-centered cubic (bcc) matrix. These structural transitions are facilitated by irradiation-induced shear. They are further stabilized by the presence of oxygen, which promotes stacking fault formation and vacancy trapping. Elevated temperatures enhance defect mobility, enabling dynamic recrystallization into ultrafine grains and the thickening of GBs due to defect absorption and impurity segregation. These findings unravel non-equilibrium phase transformation pathways in heavy ion irradiated W and highlight the critical role of impurity-mediated defect dynamics in governing its radiation tolerance properties.

Wazeer, Adil [Purdue University]↗
NASA NTRSDate not supplied

A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations

Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.

Radiation↗
NASA NTRS2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗
DOE OSTI2025

Lepto-axiogenesis with light right-handed neutrinos

We study lepto-axiogenesis in theories where the right-handed neutrino is light enough that its dynamics affect the determination of the baryon asymmetry. When compared with theories of high-scale lepto-axiogenesis where the Majorana neutrino mass may be treated as an effective dimension-five operator, we find that the predicted saxion mass is lower. Two distinct scenarios emerge. In the first, processes that generate the baryon asymmetry are in equilibrium down to the mass of the right-handed neutrino. In the second, the relevant processes never reach equilibrium; the baryon number freezes in. We comment on implications for supersymmetric spectra and discuss constraints on late decays of supersymmetric relics and from dark radiation. In contrast to high-scale lepto-axiogenesis, which predicts superpartners with masses of 10–100 TeV or more, we find this scenario is consistent with a wider range of superpartner masses, all the way down to current direct search bounds.

Axions and ALPs↗
NASA NTRS1977

The Effect of Temperature and Pressure on the Distribution of Iron Group Elements Between Metal and Olivine Phases in the Process of Differentiation of Protoplanetary Material

The distribution patterns of Ni, Co, Mn, and Cr were studied in olivines of various origins: from meteorites (chondrites, achondrites, pallasites), which are likely analogs of the protoplanetary material, to peridotite inclusions in kimberlite pipes, which are analogs of mantle material. According to X-ray microanalysis data, each genetic group of olivines is characterized by a specific concentration of these elements. Nickel is concentrated (up to 0.34 percent) in peridotite olivines, while manganese is concentrated in meteoritic olivines. The maximum chromium content (0.2 percent) was found in ureilites, which were formed under reducing conditions. Experiments at pressures of 20 to 70 kbar and temperatures of 1100 to 2000°C have shown that in a mixture of olivine and Ni metal or NiO nickel enters the silicate phase (up to 4 percent), displacing Fe into the metallic phase. Equilibrium temperatures were estimated from the Fe, Ni distribution coefficients between the metal and olivine: 1500 K for pallasites, 1600 K for olivine-bronzite H6 chondrites, 1200 K for olivine-hypersthene L6, 900 K for LL6, and 1900 K for ureilites (at P = 1 atm). The equilibrium conditions of peridotites are close to T = 1800 K and P over 100 kbar. The distribution patterns of the transition elements are explained on the basis of physical-chemical properties. It is concluded that there is a sharp difference between the conditions of differentiation of the protoplanetary material at the time meteorites were formed and the conditions of differentiation of the planets into concentric layers.

A P Vinogradov↗
NASA NTRS1966

Non-Metallic Materials

Nonmetallic materials development - cryogenic insulation, adhesives research, and membrane diffusion theory

Materials Science↗
NASA NTRS1966

Non-magnetic battery case Patent

Nonmagnetic hermetically sealed battery case made of epoxy resin and woven glass tape for use with electrochemical cells in spacecraft

Hennigan, T. J.↗
DOE OSTI2026

The impact of metastability on the high-pressure behavior of cerium

The structures adopted by solids under pressure are often assumed to reflect thermodynamic equilibrium, yet in many materials phase selection is strongly influenced by kinetic pathways and microstructural inheritance. Elemental cerium (Ce) exemplifies this challenge, with decades of conflicting reports describing different high-pressure crystal structures emerging under nominally identical conditions. Here we use neutron diffraction from large ( ~ 60 mm 3 ) sample volumes to follow the structural evolution in ultra-high-purity Ce during controlled pressure-temperature cycling between 85 and 295 K and up to 8 GPa. We find that the crystal structure formed at high pressure depends on the compression pathway: slow compression ( ~ 0.25 GPa hr −1 ) at room temperature favors an orthorhombic phase (α'), whereas slow ( ~ 0.25 GPa hr −1 ) and also moderately faster ( ~ 0.5 GPa hr −1 ) compression at low temperature stabilizes a pure monoclinic phase (α"). The low-temperature phase persists metastably over a wide temperature range but transforms irreversibly upon heating above ~ 280 K, or modest pressure cycling. Remarkably, the lower-pressure γ phase remains trapped far beyond the expected stability field, persisting to the highest pressures studied. These observations show that phase selection in Ce is governed by kinetics and microstructural memory rather than equilibrium thermodynamics or sample morphology alone, establishing path dependence as a defining feature of its high-pressure behavior.

Ridley, Christopher J. [Oak Ridge National Laborat↗