Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonequilibrium dynamics”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

The development of an explicit thermochemical nonequilibrium algorithm and its application to compute three dimensional AFE flowfields

This study presents a three-dimensional explicit, finite-difference, shock-capturing numerical algorithm applied to viscous hypersonic flows in thermochemical nonequilibrium. The algorithm employs a two-temperature physical model. Equations governing the finite-rate chemical reactions are fully-coupled to the gas dynamic equations using a novel coupling technique. The new coupling method maintains stability in the explicit, finite-rate formulation while allowing relatively large global time steps. The code uses flux-vector accuracy. Comparisons with experimental data and other numerical computations verify the accuracy of the present method. The code is used to compute the three-dimensional flowfield over the Aeroassist Flight Experiment (AFE) vehicle at one of its trajectory points.

Palmer, Grant↗
DOE OSTI2023

Visualizing interfacial collective reaction behaviour of Li–S batteries

Benefiting from high energy density (2,600 Wh kg -1 )) and low cost, lithium-sulfur (Li-S) batteries are considered promising candidates for advanced energy-storage systems. Despite tremendous efforts in suppressing the long-standing shuttle effect of lithium polysulfides, understanding of the interfacial reactions of lithium polysulfides at the nanoscale remains elusive. This is mainly because of the limitations of in situ characterization tools in tracing the liquid-solid conversion of unstable lithium polysulfides at high temporal-spatial resolution. There is an urgent need to understand the coupled phenomena inside Li-S batteries, specifically, the dynamic distribution, aggregation, deposition and dissolution of lithium polysulfides. Here, by using in situ liquid-cell electrochemical transmission electron microscopy, we directly visualized the transformation of lithium polysulfides over electrode surfaces at the atomic scale. Notably, an unexpected gathering-induced collective charge transfer of lithium polysulfides was captured on the nanocluster active-centre-immobilized surface. It further induced an instantaneous deposition of nonequilibrium Li2S nanocrystals from the dense liquid phase of lithium polysulfides. Without mediation of active centres, the reactions followed a classical single-molecule pathway, lithium polysulfides transforming into Li 2 S 2 and Li 2 S step by step. Molecular dynamics simulations indicated that the long-range electrostatic interaction between active centres and lithium polysulfides promoted the formation of a dense phase consisting of Li + and S n 2- (2 < n ≤ 6), and the collective charge transfer in the dense phase was further verified by ab initio molecular dynamics simulations. Finally, the collective interfacial reaction pathway unveils a new transformation mechanism and deepens the fundamental understanding of Li-S batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

Vortex-Controlled Quasiparticle Multiplication and Self-Growth Dynamics in Superconducting Resonators

Even in the quantum limit, non-equilibrium quasiparticle (QP) populations induce QP poisoning that irreversibly relaxes the quantum state and significantly degrades the coherence of transmon qubits. A particularly detrimental yet previously unexplored mechanism arises from QP multiplication facilitated by vortex trapping in superconducting quantum circuits, where a high-energy QP relaxes by breaking additional Cooper pairs and amplifying the QP population due to the locally reduced excitation gap and enhanced quantum confinement within the vortex core. Here we directly resolve this elusive QP multiplication process by revealing vortex-controlled QP self-generation in a highly nonequilibrium regime preceding the phonon bottleneck of QP relaxation. At sufficiently low fluence, femtosecond-resolved magneto-reflection spectroscopy directly reveals a continuously increasing QP population that is strongly dependent on magnetic-field-tuned vortex density and absent at higher excitation fluences. Quantitative analysis of the emergent QP pre-bottleneck dynamics further reveals that, although the phonon population saturates within $\simeq$10~ps, both free and trapped QPs continue to grow in a self-sustained manner--hallmarks of the long-anticipated QP-vortex interactions in nonequilibrium superconductivity. We estimate a substantial increase of $\sim$34% in QP density at vortex densities of $\sim$ 100 magnetic flux quanta per $\mathrm{μm^{2}}$. Our findings establish a powerful spectroscopic tool for uncovering QP multiplication and reveal vortex-assisted QP relaxation as a critical materials bottleneck whose mitigation will be essential for resolving QP poisoning and enhancing coherence in superconducting qubits.

Park, Joong M. [Ames Lab]↗
DOE OSTI2022

Langevin dynamics of generalized spins as SU( N ) coherent states

Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude S have N = 2⁢S + 1 levels. Alternatively, the spin state is fully characterized by a set of N 2 - 1 local physical observables, which we interpret as generalized spin components. For example, a spin with S = 1 has three dipole components and five quadrupole components. Further, these components evolve under a generalization of the classical Landau-Lifshitz dynamics, which can be extended with noise and damping terms. In this paper, we reformulate the dynamical equations of motion as a Langevin dynamics of SU(N) coherent states in the Schrödinger picture. This viewpoint is especially useful as the basis for an efficient numerical method to sample spin configurations in thermal equilibrium and to simulate the relaxation and driven motion of topological solitons. To illustrate the approach, we simulate a nonequilibrium relaxation process that creates CP 2 skyrmions, which are topological defects with both dipole and quadrupole character.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

From predicting to learning dissipation from pair correlations of active liquids

Active systems, which are driven out of equilibrium by local non-conservative forces, can adopt unique behaviors and configurations. Towards designing such materials, an important challenge is to precisely connect the static structure of active systems to the dissipation of energy induced by the local driving. Here, we use tools from liquid-state theories and machine learning to take on these challenges. We first demonstrate analytically for an isotropic active matter system that dissipation and pair correlations are closely related when driving forces behave like an active temperature. We then extend a nonequilibrium mean-field framework for predicting these pair correlations which, unlike most existing approaches, is applicable even for strongly interacting particles and far from equilibrium, to predict dissipation in these systems. Based on this theory, we reveal analytically a robust relation between dissipation and structure which holds even as the system approaches a nonequilibrium phase transition. Finally, we construct a neural network which maps static configurations of particles to their dissipation rate without any prior knowledge of the underlying dynamics. Furthermore, our results open novel perspectives on the interplay between dissipation and organization out-of-equilibrium.

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

Achieving designed texture and flows in bulk active nematics using optimal control theory

Being intrinsically nonequilibrium, active materials can potentially perform functions that would be thermodynamically forbidden in passive materials. However, active systems have diverse local attractors that correspond to distinct dynamical states, many of which exhibit chaotic turbulent-like dynamics and thus cannot perform work or useful functions. Designing such a system to choose a specific dynamical state is a formidable challenge. Motivated by recent advances enabling optogenetic control of experimental active materials, we describe an optimal control theory framework that identifies a spatiotemporal sequence of light-generated activity that drives an active nematic system toward a prescribed dynamical steady state. Active nematics are unstable to spontaneous defect proliferation and chaotic streaming dynamics in the absence of control. We demonstrate that optimal control theory can compute activity fields that redirect the dynamics into a variety of alternative dynamical programs and functions. This includes dynamically reconfiguring between states, selecting and stabilizing emergent behaviors that do not correspond to attractors, and are hence unstable in the uncontrolled system. Furthermore, our results provide a roadmap to leverage optical control methods to rationally design structure, dynamics, and function in a wide variety of active materials.

Complex systems theory↗
DOE OSTI2024

Bipartite mutual information in classical many-body dynamics

Information theoretic measures have helped to sharpen our understanding of many-body quantum states. As perhaps the most well-known example, the entanglement entropy (or more generally, the bipartite mutual information) has become a powerful tool for characterizing the dynamical growth of quantum correlations. By contrast, although computable, the bipartite mutual information (MI) is almost never explored in classical many particle systems; this owes in part to the fact that computing the MI requires keeping track of the evolution of the full probability distribution, a feat which is rarely done (or thought to be needed) in classical many-body simulations. Here, we utilize the MI to analyze the spreading of information in 1D elementary cellular automata (CA). Broadly speaking, we find that the behavior of the MI in these dynamical systems exhibits a few different types of scaling that roughly correspond to known CA universality classes. Of particular note is that we observe a set of automata for which the MI converges parametrically slowly to its thermodynamic value. We develop a microscopic understanding of this behavior by analyzing a two-species model of annihilating particles moving in opposite directions. Furthermore, our work suggests the possibility that information theoretic tools such as the MI might enable a more fine-grained characterization of classical many-body states and dynamics.

Cellular automata↗
NASA NTRS1991

Algorithm Computes Hypersonic Flow Of Air

Algorithm simulates three-dimensional hypersonic flow of air or another mixture of gases around blunt body. Accounts for nonequilibrium thermochemical effects via two-temperature mathematical model of molecular excitations and equations for finite-rate chemical reactions. These equations fully coupled to fluid dynamical equations.

Palmer, Grant↗
DOE OSTI2023

Convex Optimization for Nonequilibrium Steady States on a Hybrid Quantum Processor

Finding the transient and steady state properties of open quantum systems is a central problem in various fields of quantum technologies. Here, in this work, we present a quantum-assisted algorithm to determine the steady states of open system dynamics. By reformulating the problem of finding the fixed point of Lindblad dynamics as a feasibility semidefinite program, we bypass several well-known issues with variational quantum approaches to solving for steady states. We demonstrate that our hybrid approach allows us to estimate the steady states of higher dimensional open quantum systems and discuss how our method can find multiple steady states for systems with symmetries.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

Shock Wave Dynamics in Weakly Ionized Plasmas

An investigation of the dynamics of shock waves in weakly ionized argon plasmas has been performed using a pressure ruptured shock tube. The velocity of the shock is observed to increase when the shock traverses the plasma. The observed increases cannot be accounted for by thermal effects alone. Possible mechanisms that could explain the anomalous behavior include a vibrational/translational relaxation in the nonequilibrium plasma, electron diffusion across the shock front resulting from high electron mobility, and the propagation of ion-acoustic waves generated at the shock front. Using a turbulence model based on reduced kinetic theory, analysis of the observed results suggest a role for turbulence in anomalous shock dynamics in weakly ionized media and plasma-induced hypersonic drag reduction.

Johnson, Joseph A., III↗
DOE OSTI2022

Transient X-ray Diffraction Reveals Nonequilibrium Phase Transition in Thin Films of CH 3 NH 3 PbI 3 Perovskite

Advantageous optoelectronic properties of methylammonium lead triiodide likely arise from coupling between photogenerated charge-carriers and the soft, deformable lattice. We investigate structural dynamics of MAPbI 3 films using time-resolved X-ray diffraction versus pump-probe time delay and pump intensity. During the first nanosecond, the lattice anisotropically distorts from tetragonal to cubic at excitation intensities that are insufficient to thermally induce the first-order thermodynamic phase transition at 330K. Further, the high-symmetry structure then relaxes back to the starting phase with 11 and 236-ns time constants via a different transition pathway than observed either in the first nanosecond or in previous reports for MAPbI 3 . Early-time dynamics are consistent with polaron formation and lattice strain stabilization while the slower recovery dynamics outlive radiative recombination and relates metastability. Fluence-independence of these lattice deformations in the low-power regime conveys relevance to optoelectronics including photovoltaics and highlights sustained involvement of non-equilibrium, photoinduced lattice reorganization in MAPbI 3 under functional conditions.

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

Imaging Carrier Dynamics and Transport in Hybrid Perovskites with Transient Absorption Microscopy

In this review, the recent progress in using transient absorption microscopy to image charge transport and dynamics in semiconducting hybrid organic–inorganic perovskites is discussed. The basic principles, instrumentation, and resolution of transient absorption microscopy are outlined. With temporal resolution as high as 10 fs, sub-diffraction-limit spatial resolution, and excited-state structural resolution, these experiments have provided crucial details on charge transport mechanisms that have been previously obscured in conventional ultrafast spectroscopy measurements. Morphology-dependent mapping unveils spatial heterogeneity in carrier recombination and cooling dynamics. By spatially separating the pump and probe beams, carrier transport across grain boundaries has been directly visualized. Further, femtosecond temporal resolution allows for the examination of nonequilibrium transport directly, revealing extraordinarily long-range hot carrier migration. The application of transient absorption microscopy is not limited to hybrid perovskites but can also be useful for other polycrystalline materials in which morphology plays an important role in carrier transport.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the theory of nonhomogeneous nonequilibrium superconductivity in 2D systems with massless fermions

Here we analyze static and nonequilibrium superconducting properties of a 2D relativistic-like model system with local electron-electron interaction, Rashba spin-orbit interaction αR in presence of time-dependent in-plane magnetic field H(t). It is shown that similar to the 2D case with ordinary massive quasiparticle dispersion ε(k)~|k| 2 at large fields, such a system demonstrates a nonhomogeneous superconducting stripe phase with the order parameter Δ(r)=Δ(0)cos(2[μ B B×r]n/$\hbarυ_F$) (B is the magnetic induction, υF is the Fermi velocity, n is the normal to the plane, μ B is the Bohr magneton, and αR$\ll$υF) where the stripes are oriented along the B direction. In the considered system, the inter-stripe period L and the magnitude of the magnetic field B are related by a universal relation BL=$\hbarυ_F$/μ B ≃0.714∙10 -4 Tm. Contrary to the case of massive quasiparticles, where the condition α R ~υ F can be, in principle, satisfied by increasing α R or by charge doping (Fermi velocity decreasing), in a relativistic-like system, where υF is doping-independent and one-two orders of magnitude larger than typical Fermi velocity in the “standard” 2D systems, the stripe phase can be the ground state at a rather low doping level. We also analyzed the nonequilibrium properties of the system with a focus on the melting of the stripe order (when the magnetic field is quenched to a lower value) and stripe dynamics (when the field is rotated by 90° degrees) and found several notable results. In particular, it was shown that the stripe domains melt according to law R~1$\sqrt{t}$ at initial times, while at longer times they shrink exponentially. In the case of the flipped magnetic field, the stripe orientation gradually turns from x- to y-direction, and the intermediate “crossed-stripe” phase takes place during times of order of picoseconds. Such a crossed phase is built of periodic superconducting bubbles that potentially may have applications in modern ultrafast superconducting technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Real-time simulation of asymmetry generation in fermion-bubble collisions

Motivated by the out-of-equilibrium dynamics during an early-Universe first-order phase transition, we perform real-time simulations of fermion-bubble scattering in 1 + 1 dimensions. This nonequilibrium process can generate a charge-conjugation C asymmetry outside the bubble wall, induced by the complex fermion mass profile. The resulting C asymmetry is the 1 + 1 -dimensional analog of the C P asymmetry in 3 + 1 dimensions, a key ingredient in baryon asymmetry generation at the electroweak scale. Using tensor network methods, we track the real-time evolution of the C asymmetry in the charge density as the fermion interacts with the bubble wall, a regime inaccessible to analytic calculations. We further introduce two observables to quantify the asymmetry in the asymptotic region where reflected particles are well separated from the scattering point: one based on the net charge outside the bubble wall, and the other on the spatial displacement between the reflected particle and antiparticle wave packets. Our study represents a first step toward nonperturbative, real-time computations of C P asymmetry in 3 + 1 dimensions for electroweak baryogenesis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Large-scale gas dynamical processes affecting the origin and evolution of gaseous galactic halos

Observations of galactic halo gas are consistent with an interpretation in terms of the galactic fountain model in which supernova heated gas in the galactic disk escapes into the halo, radiatively cools and forms clouds which fall back to the disk. The results of a new study of several large-scale gas dynamical effects which are expected to occur in such a model for the origin and evolution of galactic halo gas will be summarized, including the following: (1) nonequilibrium absorption line and emission spectrum diagnostics for radiatively cooling halo gas in our own galaxy, as well the implications of such absorption line diagnostics for the origin of quasar absorption lines in galactic halo clouds of high redshift galaxies; (2) numerical MHD simulations and analytical analysis of large-scale explosions ad superbubbles in the galactic disk and halo; (3) numerical MHD simulations of halo cloud formation by thermal instability, with and without magnetic field; and (4) the effect of the galactic fountain on the galactic dynamo.

Shapiro, Paul R.↗
DOE OSTI2021

Metastability of Liquid Water Freezing into Ice VII under Dynamic Compression

The ubiquitous nature and unusual properties of water have motivated many studies on its metastability under temperature- or pressure-induced phase transformations. Here, nanosecond compression by a high-power laser is used to create the nonequilibrium conditions where liquid water persists well into the stable region of ice VII. Through our experiments, as well as a complementary theoretical/computational analysis based on classical nucleation theory, we report that the metastability limit of liquid water under nearly isentropic compression from ambient conditions is at least 8 GPa, higher than the 7 GPa previously reported for lower loading rates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

$${\mathcal{P}}{\mathcal{T}}$$-Symmetric Effective Model for Nonequilibrium Phase Transitions in a Dissipative Fermionic Mott Insulator Chain

Nonequilibrium phase transitions in open dissipative systems can be described as instabilities in the spectra and wavefunctions of effective non-Hermitian Hamiltonians invariant under simultaneous parity (P) and time-reversal (T) transformations. The degree of non-Hermiticity reflects the strength of the external drive and dissipation, and the transition is described as a loss of the PT symmetry of the solutions corresponding to stationary low-drive dynamics. This approach has been successfully applied to spin, superconducting, and Mott insulator systems. However, the microscopic foundations for the employed phenomenological models are currently lacking. Here we propose a microscopic mechanism leading to the PT-symmetric effective model in the context of the nonequilibrium Mott transition in a dissipative Hubbard chain. Our model comprises a half-filled fermionic Hubbard chain subject to a constant electric field. The dissipation is introduced via the electron-phonon coupling. We obtain the explicit expressions for the non-Hermitian parameter in terms of the electron-phonon coupling strength and driving field. Analyzing the implications of microscopic model, we find a re-entrant Mott insulator with the increasing electric field for phonon density of states that increases slower than the square of the energy (such as in one or two dimensions), or varies non-monotonously with energy.

36 MATERIALS SCIENCE↗