Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 55 records · Page 3

Equilibrium-gated pattern formation: How molecular dissociation thermodynamics drive emergent behavior in dissipative polymeric systems

Emergent patterns in biological systems arise through dissipative processes that balance reaction and transport phenomena, producing highly functional properties from self-regulating mechanisms. Synthetic fabrication, by contrast, often relies on user-controlled, multistep methods that lack the self-organizing capabilities of natural systems. Inspired by nature, we sought chemical systems that integrate strongly coupled reaction and transport phenomena, identifying frontal ring-opening metathesis polymerization (FROMP) as a method capable of creating diverse forms and functions through reactive processing. By employing discrete molecular initiators, FROMP allows precise control of key reaction steps—inhibition, initiation, and propagation. Using an integrated computational and experimental framework, we uncover how near-equilibrium inhibition dynamics, coupled with far-from-equilibrium reaction kinetics, drive pattern formation in frontally polymerized synthetic materials. We propose the concept of equilibrium-gated pattern formation, demonstrating how initiator chemistry can be tuned to achieve programmable macroscale properties. Our study reveals a surprising insight: Emergent behavior in FROMP systems arises from the inhibition-dominated regime of resin composition, expanding prior observations that such behavior is confined to a narrow compositional space near the boundary between front quenching and uniform front propagation. We identify a broader compositional window, far from the quenching regime, where emergent behavior reliably manifests. This expanded design space significantly enhances the operational flexibility of reactive systems and their capacity for self-organization. Furthermore, these insights provide a roadmap for designing bioinspired materials with self-organizing capabilities, unlocking possibilities in synthetic manufacturing.

Architected Polymers↗

Assessing the numerical stability of physics models to equilibrium variation through database comparisons on DIII-D

High fidelity kinetic equilibria are crucial for tokamak modeling and analysis. Manual workflows for constructing kinetic equilibria are time consuming and subject to user error, motivating development of automated equilibrium reconstruction tools to provide accurate and consistent reconstructions for downstream physics analysis. These automated tools also provide access to kinetic equilibria at large database scales, which enables the quantification of general uncertainties arising from equilibrium reconstruction techniques. In this paper, we compare a large database of DIII-D kinetic equilibria generated manually by physics experts to equilibria from automated kinetic reconstruction tools, assessing the impact of reconstruction method on equilibrium parameters and resulting magnetohydrodynamic stability calculations. We find agreement among scalar parameters, whereas profile quantities, such as the bootstrap current, show larger disagreements. We analyze ideal kink and classical tearing stability with DCON and STRIDE respectively, finding that the kink stability calculation is generally more robust than the tearing index Δ' calculation. We find that in 90% of cases, both kink stability classifications are unchanged between the manual expert and automated kinetic equilibria.

CAKE↗

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants

We study out-of-thermodynamic-equilibrium effects in neutron-star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During mergers, the cores of the neutron stars remain cold (T ~ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ~2 to 3 ms matter and neutrinos reach equilibrium everywhere in the remnant massive neutron star. Furthermore, our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance of wave function and Green's function methods for non-equilibrium many-body dynamics

Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime. Published by the American Physical Society 2025

Reeves, Cian C. (ORCID:0009000642581845)↗

Three-beam hybrid fs/ps coherent anti-Stokes Raman scattering of rotation–vibration non-equilibrium

In this work, we demonstrate time-resolved rotational and vibrational temperature measurements to probe plasma non-equilibrium through a simple three-beam hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS) system. A single pump/Stokes pair is employed to generate both the pure-rotational and ro-vibrational Raman coherence. A novel phase-matching scheme, to our knowledge, is employed to spatially overlap the two CARS signals, simplifying the optical design and allowing signal selection by tuning a single spectrometer grating. Measurements were performed near the electrodes in a N 2 DC glow discharge. Both the rotational temperature and the vibrational populations up to $v$ = 8 were calculated from the separately measured single-shot pure-rotational and ro-vibrational CARS spectra. Strong rotation–vibration non-equilibrium was observed at both electrodes, with the cathode showing higher vibrational and rotational temperatures. In addition, non-Boltzmann behaviors were observed at both electrodes. This simplified approach enables measurements of rotational and vibrational temperatures with high spatial resolution in non-equilibrium flows.

42 ENGINEERING↗

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

This report describes how to produce reference solutions for the equilibrium core search of pebble bed reactors using the pyrates Python library. The pyrates library uses the kugelpy methodology to perform running-in calculations to reach an equilibrium core. This approach relies on full-core Monte Carlo neutron transport calculations and the shifting of fuel pebble inventory through the reactor to simulate pebble movement in the core from one step to the next of a running-in scenario. In addition to the use of Serpent for the Monte Carlo neutron transport calculations, as part of this work, updates for the use of the Shift Monte Carlo code were integrated into the official pyrates GitHub repository. Comparisons of running-in simulations for a generic pebble bed reactor produced by using pyrates with Serpent and Shift are shown. Although consistent trends are shown between the use of the two Monte Carlo codes, the limiting factor in either calculation is the computation time due to the need to perform several hundred full-core neutron transport calculations before an equilibrium core is reached.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

This report describes how to produce reference solutions for the equilibrium core search of pebble bed reactors using the pyrates Python library. The pyrates library uses the kugelpy methodology to perform running-in calculations to reach an equilibrium core. This approach relies on full-core Monte Carlo neutron transport calculations and the shifting of fuel pebble inventory through the reactor to simulate pebble movement in the core from one step to the next of a running-in scenario. In addition to the use of Serpent for the Monte Carlo neutron transport calculations, as part of this work, updates for the use of the Shift Monte Carlo code were integrated into the official pyrates GitHub repository. Comparisons of running-in simulations for a generic pebble bed reactor produced by using pyrates with Serpent and Shift are shown. Although consistent trends are shown between the use of the two Monte Carlo codes, the limiting factor in either calculation is the computation time due to the need to perform several hundred full-core neutron transport calculations before an equilibrium core is reached.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Pre-Equilibrium De-Excitations in Neutrino-Nucleus Interactions

Pre-equilibrium de-excitation is a well-established stage of nucleon-induced nuclear reactions but has not previously been incorporated into neutrino event generators. This work implements the Koning-Duijvestijn exciton model within the MARLEY neutrino event generator to simulate energy redistribution and particle emission before compound nucleus equilibrium. The implementation includes modular calculations of particle-hole state densities, internal transition rates, emission rates, and Monte Carlo cascade sampling, and is validated against the TALYS-2.2 nuclear reaction code. A new event record is also proposed to enable future uncertainty reweighting studies. This work establishes the first framework for studying pre-equilibrium effects in neutrino-induced reactions, with the goal of improving predictions of particle multiplicities, γ-ray production, and MeV-scale detector signals relevant to experiments such as DUNE.

Visser, Erin [Michigan State U., East Lansing (mai↗

Standardising the “Gregory method” for calculating equilibrium climate sensitivity

The equilibrium climate sensitivity (ECS) – the equilibrium global mean temperature response to a doubling of atmospheric CO 2 – is a high-profile metric for quantifying the Earth system's response to human-induced climate change. A widely applied approach to estimating the ECS is the “Gregory method” (Gregory et al., 2004), which uses an ordinary least squares (OLS) regression between the net radiative flux, N, and surface air temperature anomalies, ΔT, from a 150 year experiment in which atmospheric CO 2 concentrations are quadrupled. The ECS is determined by extrapolating the linear fit to N=0, i.e. the ΔT-intercept, indicating the point at which the system is back in equilibrium. This method has been used to compare ECS estimates across the CMIP5 and CMIP6 ensembles and will likely be a key diagnostic for CMIP7. Despite its widespread application, there is little consistency or transparency between studies in how the climate model data is processed prior to the regression, leading to potential discrepancies in ECS estimates. We identify 32 alternative data processing pathways, varying by differences in global mean weighting, net radiative flux variable, anomaly calculation method, and linear regression fit. Using 44 CMIP6 models, we systematically assess the impact of these choices on ECS estimates and calculate uncertainty ranges using two bootstrap approaches. While the inter-model ECS range is insensitive to the data processing pathway, individual outlier models exhibit notable differences. Approximating a model's native grid cell area (if irregular) with cosine of the latitude can decrease the ECS by 11 %, the choice of N-variable can change the ECS by 6 %, and some anomaly calculation methods can introduce spurious temporal correlations in the processed data. Beyond data processing choices, we also evaluate an alternative linear regression method – total least squares (TLS) – which has a more statistically robust basis than OLS. However, for consistency with previous literature, and given TLS may reduce the ECS compared to OLS (by up to 24 %), thereby making a known bias in the Gregory method worse, we do not feel there is sufficient clarity to recommend a transition to TLS in all cases. To improve reproducibility and comparability in future studies, we recommend a standardised Gregory method: weighting the global mean by cell area, using the top of the atmosphere (as opposed to the top of model) N-variable, and calculating anomalies by first applying a rolling average to the preindustrial control timeseries then subtracting from the raw CO 2 quadrupling experiment. This approach accounts for model drift while reducing noise in the data to best meet the pre-conditions of the linear regression. While CMIP6 results of the multi-model mean ECS appear insensitive to these processing choices, similar assumptions may not hold for CMIP7, underscoring the need for standardised data preparation in future climate sensitivity assessments.

Geosciences↗

FLORA Equilibrium and Stability Code Archive

FLORA solves, in a 2-D domain (radial and axial dimensions) with a specified azimuthal Fourier mode, for the linearized stability of a long, thin, axisymmetric plasma equilibrium in an applied magnetic field. Before the stability equation is solved, FLORA solves a set of simple equations for pressure balance that specify the equilibrium magnetic field and plasma pressure in the long-thin limit given a simplified description of the magnetic coils. It uses an initial-value method for the linear stability problem in which an equilibrium is given an initial perturbation to its magnetic field, and the temporal behavior of the perturbation is followed. The perturbation has been Fourier expanded in the azimuthal coordinate; each azimuthal mode must be examined separately. The complex partial differential equation of motion for the perturbed radial displacement of the field lines is solved as a coupled system of two real p.d.e.'s and the solution consists of two parts, the real part and the imaginary part. The system is solved by bringing the coupling terms in each equation to the right side and using an iterative technique.

Cohen, Bruce [Lawrence Livermore National Laborato↗

Non-Equilibrium Actinide Radiation Chemistry and the Nuclear Fuel Cycle

Actinides are inherently unstable elements that frequently coexist with other radioisotopes, generating intense ionizing radiation fields that drive the formation of non equilibrium oxidation states. These transient species exert a profound mechanistic influence on the radiation response of actinide containing systems due to their unique redox chemistry. Despite their importance, they remain poorly understood, yet such insight is essential for advancing actinide science and accurately predicting radiation driven behavior. Actinide separations—critical for nuclear energy technologies, strategic deterrence, space exploration, and nuclear medicine—depend on precise control of actinide oxidation states to recover targeted elements from complex matrices such as used nuclear fuel. However, during these processes, actinides, their coordination complexes, and the separation media are all exposed to intense, multicomponent (alpha, beta, gamma, etc.) radiation fields that can alter process efficiency, selectivity, and chemical stability. Understanding, controlling, and mitigating radiation induced reactions is therefore key to innovating and optimizing next generation separation technologies. This seminar will provide an overview of the nuclear fuel cycle and non equilibrium actinide radiation chemistry in the context of recovering actinides from used nuclear fuel, with a particular emphasis on direct dissolution–based reprocessing strategies. We will explore time resolved electron pulse radiolysis and alpha and gamma dose accumulation studies, integrated with multiscale computational modeling, to elucidate the molecular level roles of radiation driven, non equilibrium actinide species in process performance and in the radiolytic stability of organic ligands used for actinide recovery. These insights offer new pathways for designing advanced separation methods and next generation solvent systems, with broad implications for the future of the nuclear fuel cycle.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

Non-Equilibrium Actinide Radiation Chemistry and the Nuclear Fuel Cycle

Invited John and Naomi Fackler Lectureship in Chemistry and English seminar at Valparaiso University, IN, USA. Actinides are inherently unstable elements that frequently coexist with other radioisotopes, generating intense ionizing radiation fields that drive the formation of non-equilibrium oxidation states. These transient species exert a profound mechanistic influence on the radiation response of actinide-containing systems due to their unique redox chemistry. Despite their importance, they remain poorly understood, yet such insight is essential for advancing actinide science and accurately predicting radiation-driven behavior. Actinide separations—critical for nuclear energy technologies, strategic deterrence, space exploration, and nuclear medicine—depend on precise control of actinide oxidation states to recover targeted elements from complex matrices such as used nuclear fuel. However, during these processes, actinides, their coordination complexes, and the separation media are all exposed to intense, multicomponent (alpha, beta, gamma, etc.) radiation fields that can alter process efficiency, selectivity, and chemical stability. Understanding, controlling, and mitigating radiation-induced reactions is therefore key to innovating and optimizing next-generation separation technologies. This seminar will provide an overview of the nuclear fuel cycle and non-equilibrium actinide radiation chemistry in the context of recovering actinides from used nuclear fuel, with a particular emphasis on direct-dissolution–based reprocessing strategies. We will explore time-resolved electron pulse radiolysis and gamma dose accumulation studies to elucidate the molecular-level roles of radiation-driven, non-equilibrium actinide species in process performance and in the radiolytic stability of organic ligands used for actinide recovery. These insights offer new pathways for designing advanced separation methods and next-generation solvent systems, with broad implications for the future of the nuclear fuel cycle.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

On Gibbs Equilibrium and Hillert Nonequilibrium Thermodynamics

During his time at Royal Institute of Technology (Kungliga Tekniska högskolan) in Sweden, the present author learned nonequilibrium thermodynamics from Mats Hillert. The key concepts are the separation of internal and external variables of a system and the definitions of potentials and molar quantities. In equilibrium thermodynamics derived by Gibbs, the internal variables are not independent and can be fully evaluated from given external variables. While irreversible thermodynamics led by Onsager focuses on internal variables though often mixed with external variables. Hillert integrated them together by first emphasizing their differences and then examining their connections. His philosophy was reflected by the title of his book “Phase Equilibria, Phase Diagrams and Phase Transformations” that puts equilibrium, nonequilibrium, and internal processes on equal footing. Here, in the present paper honoring Hillert, the present author reflects his experiences with Hillert and his work in last 40 years and expresses his gratitude for all the wisdom and support from him in terms of “Hillert nonequilibrium thermodynamics” and discusses some recent topics that the present author has been working on.

36 MATERIALS SCIENCE↗

Generalizable machine learning potentials for quantum-accurate predictions of non-equilibrium behavior in 2D materials

Machine learning interatomic potentials (ML-IAPs) are emerging as transformative tools in materials modeling, promising quantum-level accuracy at a fraction of the computational cost. However, their ability to generalize beyond equilibrium configurations and to reliably capture defect- and temperature-driven behavior remains underexplored. Here, we develop and benchmark two state-of-the-art ML-IAPs, Spectral Neighbor Analysis Potential (SNAP) and Allegro, on a comprehensive dataset for monolayer MoSe₂. Using density functional theory (DFT) as the reference, we evaluate their performance in capturing stress–strain behavior, phase transition energetics, defect evolution, edge stability, and fracture toughness. Allegro, a deep equivariant neural network potential, surpasses both SNAP and the classical Tersoff potential in accuracy, efficiency, and transferability. Importantly, both ML potentials accurately reproduce experimental fracture measurements and ab initio predictions of inversion domain formation—phenomena well beyond their training sets. Our findings establish ML-IAPs as viable replacements for traditional force fields in the study of non-equilibrium mechanical phenomena, enabling large-scale, high-fidelity simulations in 2D materials and beyond. In conclusion, this work provides a broadly applicable framework for data-driven modeling of structural and functional transformations under extreme conditions.

2D materials↗

Experimental determination of hydrogen isotopic equilibrium in the system H 2 O (l) -H 2(g) from 3 to 90 °C

Molecular hydrogen (H 2 ) is found in a variety of settings on and in the Earth from low-temperature sediments to hydrothermal vents, and is actively being considered as an energy resource for the transition to a green energy future. The hydrogen isotopic composition of H 2 , given as D/H ratios or δD, varies in nature by hundreds of per mil from ∼−800 ‰ in hydrothermal and sedimentary systems to ∼+450 ‰ in the stratosphere. This range reflects a variety of processes, including kinetic isotope effects associated with formation and destruction and equilibration with water, the latter proceeding at fast (order year) timescales at low temperatures (<100 °C). At isotopic equilibrium, the D/H fractionation factor between liquid water and hydrogen ( D α H2O(l)-H2(g) ) is a function of temperature and can thus be used as a geothermometer for H 2 formation or re-equilibration temperatures. Multiple studies have produced theoretical calculations for hydrogen isotopic equilibrium between H 2 and water vapor. However, only three published experimental calibrations used in geochemistry exist for the H 2 O-H 2 system: two between 51 and 742 °C for H 2 O (g) -H 2(g) (Suess, 1949, Cerrai et al., 1954), and one in the H 2 O (l) -H 2(g) system for temperatures <100 °C (Rolston et al., 1976). Despite these calibrations existing, there is uncertainty on their accuracy at low temperatures (<100 °C; e.g., Horibe and Craig, 1995).

08 HYDROGEN↗

Large-scale tearing-mode hazard function analysis with standard matched equilibrium reconstructions

The association between features from standard tokamak equilibrium reconstructions and the onset of n = 1 tearing modes (TMs) is analyzed at scale. The TM onset rate is directly modeled with a ‘hazard’ function which gives the expected number of onsets (per unit time spent) in a given equilibrium parameter region. In particular the different statistical modeling performance achieved for magnetics-only reconstructions and motional Stark effect (MSE) enhanced reconstructions is studied. It is observed that a better hazard model for the TM onset rate can be built with the MSE-enhanced equilibria compared to the matched magnetics-only situation. This advantage disappears if internal profile details are withheld from the matched analysis. Plausibility of the hazard function is further demonstrated with visualizations of global trends in the operational space, and time-traces from specific tokamak discharges. As a result, TMs typically degrade tokamak plasma performance and may lead to plasma termination, motivating this statistical study.

equilibrium↗

Probing Non-Equilibrium Pair-Breaking and Quasiparticle Dynamics in Nb Superconducting Resonators Under Magnetic Fields

We conducted a comprehensive study of the non-equilibrium dynamics of Cooper pair breaking, quasiparticle (QP) generation, and relaxation in niobium (Nb) cut from superconducting radio-frequency (SRF) cavities, as well as various Nb resonator films from transmon qubits. Using ultrafast pump–probe spectroscopy, we were able to isolate the superconducting coherence and pair-breaking responses. Our results reveal both similarities and notable differences in the temperature- and magnetic-field-dependent dynamics of the SRF cavity and thin-film resonator samples. Moreover, femtosecond-resolved QP generation and relaxation under an applied magnetic field reveals a clear correlation between non-equilibrium QPs and the quality factor of resonators fabricated by using different deposition methods, such as DC sputtering and high-power impulse magnetron sputtering. These findings highlight the pivotal influence of fabrication techniques on the coherence and performance of Nb-based quantum devices, which are vital for applications in superconducting qubits and high-energy superconducting radio-frequency applications.

43 PARTICLE ACCELERATORS↗

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗