Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonequilibrium systems”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Transport, correlations, and chaos in a classical disordered anharmonic chain

Here, we explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $\textit{N}$, disorder strength Δ, and temperature $\textit{T}$. The conductivity $κ_N$, obtained for finite length ($\textit{N}$), saturates to a value $κ_∞$ > 0 in the large $\textit{N}$ limit, for all values of disorder strength Δ and temperature $\textit{T}$ > 0. We show evidence that for any Δ > 0 the conductivity goes to zero faster than any power of $\textit{T}$ in the ($\textit{T}$/Δ) → 0 limit, and find that the form $κ_∞$ ~ $e^{{–B|\text{ln}(CΔ/T)|}^{3}}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be $κ_N < κ_∞$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_∞$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Active Colloids as Models, Materials, and Machines

Active colloids use energy input at the particle level to propel persistent motion and direct dynamic assemblies. We consider three types of colloids animated by chemical reactions, time-varying magnetic fields, and electric currents. For each type, we review the basic propulsion mechanisms at the particle level and discuss their consequences for collective behaviors in particle ensembles. These microscopic systems provide useful experimental models of nonequilibrium many-body physics in which dissipative currents break time-reversal symmetry. Freed from the constraints of thermodynamic equilibrium, active colloids assemble to form materials that move, reconfigure, heal, and adapt. Colloidal machines based on engineered particles and their assemblies provide a basis for mobile robots with increasing levels of autonomy. This review provides a conceptual framework for understanding and applying active colloids to create material systems that mimic the functions of living matter. We highlight opportunities for chemical engineers to contribute to this growing field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Implicit Thermochemical Nonequilibrium Compressible Flow Simulations on Unstructured Grids Using GPUs

As next-generation exascale-class systems arrive, existing software must be updated accordingly to effectively utilize these systems. For high concurrency and energy efficiency, many of these systems utilize GPU architectures. In this work, we present a CUDA C++ implementation of FUN3D's thermochemical nonequilibrium capability for turbulent flows. Efficiency is demonstrated at scale using the Summit system at the Oak Ridge Leadership Computing Facility which is representative of future exascale systems. This work enables faster, higher fidelity, and scale-resolving simulations of thermochemical nonequilibrium flows including reentry, hypersonics, and combustion.

CFD, GPU, HPC, Hypersonics, Chemistry↗
DOE OSTI2021

Superconducting-like response in driven systems near the Mott transition

Here we point out that fractionalized bosonic charge excitations can explain the recently discovered photoinduced superconducting-like response in $\kappa$-(ET) 2 Cu [N(CN) 2 ] Br, an organic metal close to the Mott transition. The pump laser exerts a periodic drive on the fractionalized field, creating a nonequilibrium condensate, which gives a Drude peak much narrower than the equilibrium scattering rate, hence a superconducting-like response. Our proposal illuminates alternative possibilities of detecting fractionalization and can be readily tested in spin liquid candidates and in cold atom systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
DOE OSTI2024

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2008

A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements.

Johnston, Christopher O.↗
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2022

Anomalous phase separation in a correlated electron system: Machine-learning–enabled large-scale kinetic Monte Carlo simulations

Significance Phase separation is crucial to the functionalities of many correlated electron materials with notable examples including colossal magnetoresistance in manganites and high- T c superconductivity in cuprates. However, the nonequilibrium phase-separation dynamics in such systems are poorly understood theoretically, partly because the required multiscale modeling is computationally very demanding. With the aid of machine-learning methods, we have achieved large-scale dynamical simulations in a representative correlated electron system. We observe an unusual relaxation process that is beyond the framework of classical phase-ordering theories. We also uncover a correlation-induced freezing behavior, which could be a generic feature of phase separation in correlated electron systems.

Zhang, Sheng↗
DOE OSTI2025

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ergodicity breaking in rapidly rotating C 60 fullerenes

Ergodicity, the central tenet of statistical mechanics, requires an isolated system to explore all available phase space constrained by energy and symmetry. Mechanisms for violating ergodicity are of interest for probing nonequilibrium matter and protecting quantum coherence in complex systems. Polyatomic molecules have long served as a platform for probing ergodicity breaking in vibrational energy transport. Here, we report the observation of rotational ergodicity breaking in an unprecedentedly large molecule, 12 C 60 , determined from its icosahedral rovibrational fine structure. The ergodicity breaking occurs well below the vibrational ergodicity threshold and exhibits multiple transitions between ergodic and nonergodic regimes with increasing angular momentum. These peculiar dynamics result from the molecule’s distinctive combination of symmetry, size, and rigidity, highlighting its relevance to emergent phenomena in mesoscopic quantum systems.

Science & Technology - Other Topics↗
DOE OSTI2024

Probabilistic measures for biological adaptation and resilience

This paper introduces an approach to quantifying ecological resilience in biological systems, particularly focusing on noisy systems responding to episodic disturbances with sudden adaptations. Incorporating concepts from nonequilibrium statistical mechanics, we propose a measure termed “ecological resilience through adaptation,” specifically tailored to noisy, forced systems that undergo physiological adaptation in the face of stressful environmental changes. Randomness plays a key role, accounting for model uncertainty and the inherent variability in the dynamical response among components of biological systems. Our measure of resilience is rooted in the probabilistic description of states within these systems and is defined in terms of the dynamics of the ensemble average of a model-specific observable quantifying success or well-being. Our approach utilizes stochastic linear response theory to compute how the expected success of a system, originally in statistical equilibrium, dynamically changes in response to a environmental perturbation and a subsequent adaptation. Importantly, the resulting mathematical derivations allow for the estimation of resilience in terms of ensemble averages of simulated or experimental data. Finally, through a simple but clear conceptual example, we illustrate how our resilience measure can be interpreted and compared to other existing frameworks in the literature. The methodology is general but inspired by applications in plant systems, with the potential for broader application to complex biological processes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Optimizing Ammonia Separation via Reactive Absorption for Sustainable Ammonia Synthesis

Metal halide salts such as magnesium chloride have been shown to be promising candidates for ammonia storage materials for energy storage and agriculture applications due to their ability to incorporate several moles of ammonia per mole of salt. Ammonia exiting a synthesis reactor can be separated from nitrogen and hydrogen by absorption into magnesium chloride. Such an absorption can be more complete and hotter than separation via ammonia condensation, the current standard in the Haber–Bosch process. Here, we discuss the optimal conditions for the cyclic uptake and release of ammonia from the supported magnesium chloride absorbents. An automated system was designed for measuring the nonequilibrium working capacity of the absorbent, as well as the impact of important operating conditions such as absorption and desorption temperature, pressure, and desorption time. Measurements of absorption and desorption kinetics provide insight into the mechanisms involved. The temperatures and pressures during absorption and desorption were designed to use minimal energy input to maximize the uptake and release of ammonia within a reasonable amount of time. In a laboratory-scale bed, absorption has a small unused bed length, so it is largely independent of mass transfer; it is dominated by how fast ammonia is fed to the bed. Yet, desorption is restricted both by the speed of heating the bed and by diffusion out of the absorbent. These measurements provide guidelines for ammonia separations and cycling sorbent materials on a larger scale.

10 SYNTHETIC FUELS↗
DOE OSTI2022

Reversible to irreversible transitions for cyclically driven particles on periodic obstacle arrays

Here, we examine the collective dynamics of disks moving through a square array of obstacles under cyclic square wave driving. Below a critical density, we find that the system organizes into a reversible state in which the disks return to the same positions at the end of every drive cycle. Above this density, the dynamics are irreversible and the disks do not return to the same positions after each cycle. The critical density depends strongly on the angle θ between the driving direction and a symmetry axis of the obstacle array, with the highest critical densities appearing at commensurate angles such as θ = 0° and θ = 45° and the lowest critical densities falling at θ = arctan (0.618), the inverse of the golden ratio, where the flow is the most degenerate. As the density increases, the number of cycles required to reach a reversible state grows as a power law with an exponent near ν = 1.36, similar to what is found in periodically driven colloidal and superconducting vortex systems.

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

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗
DOE OSTI2026

Field-driven ion pairing dynamics in concentrated electrolytes

Here, we investigate ion pairing dynamics in electrolytes driven far from equilibrium using molecular simulations and nonequilibrium rate theory. Focusing on 0.5M LiPF 6 in water and acetonitrile under uniform electric fields, we compute transition path theory observables, including reactive fluxes and mean first-passage times of ion pairing. Moreover, we introduce a dynamical proxy of free-ion population, where its field-induced change is strongly correlated with the nonlinear enhancement of conductivity, yielding an increase of 40% at 50 mV/Å in acetonitrile, compared to that of less than 10% in aqueous electrolytes. Further kinetic analysis elucidates that Onsager’s classical theory substantially overestimates field-induced enhancement of ion pair dissociation in molecular electrolytes. This discrepancy arises from solvent-mediated dynamical pathways and field-induced dielectric decrement that suppress ion pair dissociation within explicit solvents, highlighting that a faithful description of molecular details is essential. Our results provide a molecular interpretation of nonlinear electrolyte transport beyond continuum theories and establish a general framework for quantifying nonequilibrium reaction kinetics in condensed phase systems.

Moon, Seokjin [University of California, Berkeley,↗
DOE OSTI2026

Floquet operator dynamics and orthogonal polynomials on the unit circle

Operator spreading under stroboscopic time evolution under a unitary is studied. An operator Krylov space is constructed and mapped to orthogonal polynomials on a unit circle (OPUC), as well as to the Krylov space of the edge operator of the Floquet transverse field Ising model with inhomogeneous couplings (ITFIM). The Verblunsky coefficients in the OPUC representation are related to the Krylov angles parametrizing the ITFIM. The relations between the OPUC and spectral functions are summarized and several applications are presented. These include derivation of analytic expressions for the OPUC under persistent m-periodic dynamics, and the numerical construction of the OPUC for autocorrelations of the homogeneous Floquet-Ising model as well as the Z 3 clock model. The numerically obtained Krylov angles of the Z 3 clock model with long-lived period tripled autocorrelations show a spatial periodicity of six, and this observation is used to develop an analytically solvable model for the ITFIM that mimics this behavior.

Floquet systems↗
DOE OSTI2026

Floquet Spin Splitting and Spin Generation in Antiferromagnets

In antiferromagnetic spintronics, accessing the spin degrees of freedom is essential for generating spin currents and manipulating magnetic order, which generally requires lifting spin degeneracy. This is typically achieved through relativistic spin-orbit coupling or nonrelativistic spin splitting in altermagnets. Here, we propose an alternative approach: a dynamical spin splitting induced by an optical field in antiferromagnets. By coupling the driven system to a thermal bath, we demonstrate the emergence of steady-state pure spin currents as well as linear-response longitudinal and transverse spin currents. Crucially, thermal bath engineering enables a nonrelativistic Edelstein effect—the generation of a net spin accumulation—without relying on spin-orbit coupling. Furthermore, our results provide a broadly applicable and experimentally tunable route to control spins in antiferromagnets, offering new opportunities for spin generation and manipulation in antiferromagnetic spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗