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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

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

Effects of Nonequilibrium Atomic Structure on Ionic Diffusivity in LLZO: A Classical and Machine Learning Molecular Dynamics Study

To improve the performance of electrochemical devices, it is essential to understand the effects of nonequilibrium motifs in solids, such as grain boundaries, amorphous phases, and highly strained regions, on atomic-scale transport and stability. Molecular dynamics simulations are used to explore the combined effect of far-from-equilibrium atomic structures and the choice of interatomic potential on ionic diffusivity predictions for Li 7 La 3 Zr 2 O 12 (LLZO), a promising solid electrolyte for all-solid-state batteries. Amorphization and high strain are considered using both classical Buckingham interatomic potentials and machine learning force fields. Here we find that both crystalline expansion and amorphization tend to slow diffusion, although the different physical encodings in the two potentials impact the properties in different ways. We trace these variations to a combination of structural and transport factors, the contributions of which are deconvoluted computationally. Graph-based analysis reveals that the variations for amorphous LLZO arise from the connectivity of diffusion pathways within the predicted structures, which generally correlates with diffusivity and is notably higher for structures generated by the machine learning force fields. Our study provides additional insight into the relationship between atomic structure and diffusivity in LLZO, while also highlighting the need for care in choosing and validating potentials to simulate far from equilibrium structures.

25 ENERGY STORAGE↗
DOE OSTI2026

Tracking Optical Phonon Dynamics in InP Nanocrystals via Transient Absorption and Femtosecond Stimulated Raman Spectroscopy

Semiconductor nanocrystals (NCs) offer actualized and prospective utility in optoelectronic technologies, yet key aspects of their thermal and vibrational behaviors remain unresolved. Compared with bulk crystals, heat dissipation involving NCs can differ substantially owing to pervasiveness of interfacial scattering, phonon confinement, and relaxed phonon momentum selection rules, as well as influences of ligands. Here, in addition to transient absorption, we use femtosecond stimulated Raman spectroscopy to track nonequilibrium optical-phonon dynamics in InP NCs capped with myristic acid ligands and compare the same particles coated with a ZnS shell. We identify distinct phonon decay pathways in ligand-capped versus ZnS-shelled NCs. Increasing the excitation density and introducing the ZnS shell measurably modify optical-phonon lifetimes, whereas phonon formation times remain largely unchanged. In addition, we observe phonon mode softening in ligand-capped nanocrystals, consistent with lattice expansion. Together, these results demonstrate that core–shell structures can strongly govern nanocrystal thermal dissipation pathways and should be considered a key design parameter for optoelectronic operation.

dissipation↗
DOE OSTI2020

Stabilization of a Majorana zero mode through quantum frustration

We analyze a system in which a topological Majorana zero mode (tMZM) combines with a Majorana produced by quantum frustration (fMZM) to produce a novel ground state. The system that we study combines two parts, a grounded topological superconducting wire that hosts two tMZMs at its ends and an on-resonant quantum dot connected to two dissipative leads. The quantum dot with dissipative leads creates an effective two-channel Kondo (2CK) state in which quantum frustration yields an isolated fMZM at the dot. We find that coupling the dot to one end of the topological wire stabilizes the tMZM at the other end. Three routes are used to obtain these results: (i) calculation of the conductance through an auxiliary detector quantum dot, (ii) renormalization group (RG) arguments and the g theorem, and (iii) a fully nonequilibrium calculation of the I(V) curve and shot noise S(V) through the detector dot. In addition to providing a route to achieving an unpaired Majorana zero mode, this scheme provides a clear signature of the presence of the 2CK frustration-induced Majorana.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗
DOE OSTI2025

Shear Response of Ionizable Polymer Melts at the Crossover from Ionomers to Polyelectrolytes

Nonlinear shear response of polymers is affected by inherent barriers for diffusive motion, including entanglements and topology. In ionizable polymers, ionic clusters further constrain the intrinsic dynamics of the polymers, significantly enhancing their viscosity. Here, using fully atomistic molecular dynamics simulations, the nonlinear shear response of ionizable polymers is presented, across the transition from the ionomer regime where distinctive clusters dominate the structure to the polyelectrolyte regime where clusters percolate, in polystyrene randomly sulfonated with fractions of SO 3 − groups of f = 0.20 and 0.35, in pristine and tetrahydrofuran (THF) swollen polyelectrolyte melts. For f = 0.20, the ionic clusters first fracture into smaller clusters followed by splitting into individual ionic groups and eventually reform. At higher f, the clusters morph in shape but do not break under high shear. At very high shear rates, all of the chains stretch and recoil rapidly. As the shear rate is reduced, some chains stretch and recoil, while others remain largely unaffected by the shear. Macroscopically, for all systems, the shear viscosity displays initially an elastic response, followed by nonlinear shear stress overshoot and, eventually, a steady state. The evolution of viscosity with time and shear reflects that of the ionizable domains.

cluster chemistry↗
DOE OSTI2021

Twist-Angle-Dependent Ultrafast Charge Transfer in MoS 2 -Graphene van der Waals Heterostructures

Vertically stacked transition metal dichalcogenide-graphene heterostructures provide a platform for novel optoelectronic applications with high photoresponse speeds. Photoinduced nonequilibrium carrier and lattice dynamics in such heterostructures underlie these applications but have not been understood. In particular, the dependence of these photoresponses on the twist angle, a key tuning parameter, remains elusive. Here, using ultrafast electron diffraction, we report the simultaneous visualization of charge transfer and electron–phonon coupling in MoS 2 -graphene heterostructures with different stacking configurations. We find that the charge transfer timescale from MoS 2 to graphene varies strongly with twist angle, becoming faster for smaller twist angles, and show that the relaxation timescale is significantly shorter in a heterostructure as compared to a monolayer. These findings illustrate that twist angle constitutes an additional tuning knob for interlayer charge transfer in heterobilayers and deepen our understanding of fundamental photophysical processes in heterostructures, of importance for future applications in optoelectronics and light harvesting.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Time irreversibility as an indicator of approaching tipping points in Earth subsystems

With shifting environmental trends, many Earth system elements may be poised to undergo critical transitions or ‘tipping’. Reliable anticipation of these tipping elements is vital to inform policy decisions. Many of the current methods for tipping point detection are based on loss of resilience or ‘critical slowdown’ of the system as it approaches a tipping point. However, these methods are prone to false alarms; the detected slowdown may be an artifact of nonstationary noise unrelated to tipping behavior. Here, we explore the efficacy of early warning signs based on a nonequilibrium thermodynamics framework. The model-free detection method relies on the increased intrinsic time-irreversibility due to detailed balance breaking, preceding the onset of tipping or instabilities. We demonstrate that these EWSs are effective for tipping point detection and robust against false alarms due to nonstationary noise, using idealized models for two key elements of the Earth system that are prone to tipping: the Atlantic Meridional Overturning Circulation and Arctic sea-ice loss.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Experimental study of atmospheric pressure single-pulse nanosecond discharge in pin-to-pin configuration

We present an experimental study of nanosecond high-voltage discharges in a pin-to-pin electrode configuration at atmospheric conditions operating in single-pulse mode (no memory effects). Discharge parameters were measured using microwave Rayleigh scattering, laser Rayleigh scattering, optical emission spectroscopy enhanced with a nanosecond probing pulse, and fast photography. Spark and corona discharge regimes were studied for electrode gap sizes 2-10 mm and discharge pulse duration of 90 ns. The spark regime was observed for gaps < 6 mm using discharge pulse energies of 0.6-1 mJ per mm of the gap length. Higher electron number densities, total electron number per gap length, discharge currents, and gas temperatures were observed for smaller electrode gaps and larger pulse energies, reaching maximal values of about 7.5 x 10 15 cm -3 , 3.5 x 10 11 electrons per mm, 22 A, and 4,000 K (at 10 microseconds after the discharge), respectively, for a 2 mm gap and 1 mJ/mm discharge pulse energy. Initial breakdown was followed by a secondary breakdown occurring about 30-70 ns later and was associated with ignition of a cathode spot and transition to cathodic arc. Here, a majority of the discharge pulse energy was deposited into the gas before the secondary breakdown (85-89%). The electron number density after the ns-discharge pulse decayed with a characteristic time scale of 150 ns governed by dissociative recombination and electron attachment to oxygen mechanisms. For the corona regime, substantially lower pulse energies (~0.1 mJ/mm), peak conduction current (1-2 A), and electron numbers (3-5 x 10 10 electrons per mm), and gas temperatures (360 K) were observed.

42 ENGINEERING↗
DOE OSTI2022

Initial transient stage of pin-to-pin nanosecond repetitively pulsed discharges in air

In this work, evolution of parameters of nanosecond repetitively pulsed (NRP) discharges in pin-to-pin configuration in air was studied during the transient stage of initial 20 discharge pulses. Gas and plasma parameters in the discharge gap were measured using coherent microwave scattering, optical emission spectroscopy, and laser Rayleigh scattering for NRP discharges at repetition frequencies of 1, 10, and 100 kHz. Memory effects (when perturbations induced by the previous discharge pulse would not decay fully until the subsequent pulse) were detected for the repetition frequencies of 10 and 100 kHz. For 10 kHz NRP discharge, the discharge parameters experienced significant change after the first pulse and continued to substantially fluctuate between subsequent pulses due to rapid evolution of gas density and temperature during the 100 μs inter-pulse time caused by intense redistribution of the flow field in the gap on that time scale. For 100 kHz NRP discharge, the discharge pulse parameters reached a new steady-state at about five pulses after initiation. Further, this new steady-state was associated with well-reproducible parameters between the discharge pulses and substantial reduction in breakdown voltage, discharge pulse energy, and electron number density in comparison to the first discharge pulse. For repetition frequencies 1–100 kHz considered in this work, the memory effects can be likely attributed to the reduction in gas number density and increase in the gas temperature that cannot fully recover to ambient conditions before subsequent discharge pulses.

42 ENGINEERING↗
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↗
NASA NTRS1973

Theoretical and experimental studies of reentry plasmas

A viscous shock-layer analysis was developed and used to calculate nonequilibrium-flow species distributions in the plasma layer of the RAM vehicle. The theoretical electron-density results obtained are in good agreement with those measured in flight. A circular-aperture flush-mounted antenna was used to obtain a comparison between theoretical and experimental antenna admittance in the presence of ionized boundary layers of low collision frequency. The electron-temperature and electron-density distributions in the boundary layer were independently measured. The antenna admittance was measured using a four-probe microwave reflectometer and these measured values were found to be in good agreement with those predicted. Measurements were also performed with another type of circular-aperture antenna and good agreement was obtained between the calculations and the experimental results. A theoretical analysis has been completed which permits calculation of the nonequilibrium, viscous shock-layer flow field for a sphere-cone body. Results are presented for two different bodies at several different altitudes illustrating the influences of bluntness and chemical nonequilibrium on several gas dynamic parameters of interest. Plane-wave transmission coefficients were calculated for an approximate space-shuttle body using a typical trajectory.

Dunn, M. G.↗
NASA NTRS1975

Effects of plasma microfields on radiative transitions from atomic levels above the ionization threshold

The effects of plasma electric microfields on line-like optical features arising from atomic levels above the ionization threshold are investigated within the framework of the quasi-static and single-frequency dynamic-field theories of spectral-line broadening. The 2p(23)P to 1s2p(3)P and 2s2p(3)P to 1s2s(3)S transitions in helium and helium-like ions are treated as examples. The mixing of the doubly excited levels in the perturbing microfields produces Stark broadening of the emission lines and induces autoionization of the 2p(23)P level, which, unlike the 2s2p(3)P level, is metastable against autoionization in the field-free environment. Determination of the complete Stark-broadening profiles in thermal plasmas is complicated by the need to include the effects of both the (quasi-static) ion and the (dynamic) electron fields. Under nonequilibrium conditions, where electric fields from either electron or ion plasma waves can far exceed nearby particle fields, the calculation and interpretation of the line shapes may be simplified and could provide a diagnostic probe of the wave-field properties.

Davis, J.↗
NASA NTRS2016

Flow Characterization Studies of the 10-MW TP3 Arc-Jet Facility: Probe Sweeps

This paper reports computational simulations and analysis in support of calibration and flow characterization tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted in the NASA Ames 10-MW TP3 facility using flat-faced stagnation calorimeters at six conditions corresponding to the steps of a simulated flight heating profile. Data were obtained using a conical nozzle test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot pressure and heat flux probes were also performed at these arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. Two different sets of pitot pressure and heat probes were used: 9.1-mm sphere-cone probes (nose radius of 4.57 mm or 0.18 in) with null-point heat flux gages, and 15.9-mm (0.625 in) diameter hemisphere probes with Gardon gages. The probe survey data clearly show that the test flow in the TP3 facility is not uniform at most conditions (not even axisymmetric at some conditions), and the extent of non-uniformity is highly dependent on various arc-jet parameters such as arc current, mass flow rate, and the amount of cold-gas injection at the arc-heater plenum. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested. Comparisons of computations with the experimental measurements show reasonably good agreement except at the extreme low pressure conditions of the facility envelope.

Goekcen, Tahir↗
NASA NTRS2016

Flow Characterization Studies of the 10-MW TP3 Arc-Jet Facility: Probe Sweeps

This paper reports computational simulations and analysis in support of calibration and flow characterization tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted in the NASA Ames 10-MW TP3 facility using flat-faced stagnation calorimeters at six conditions corresponding to the steps of a simulated flight heating profile. Data were obtained using a conical nozzle test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot pressure and heat flux probes were also performed at these arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. Two different sets of pitot pressure and heat probes were used: 9.1-mm sphere-cone probes (nose radius of 4.57 mm or 0.18 in) with null-point heat flux gages, and 15.9-mm (0.625 in) diameter hemisphere probes with Gardon gages. The probe survey data clearly show that the test flow in the TP3 facility is not uniform at most conditions (not even axisymmetric at some conditions), and the extent of non-uniformity is highly dependent on various arc-jet parameters such as arc current, mass flow rate, and the amount of cold-gas injection at the arc-heater plenum. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested. Comparisons of computations with the experimental measurements show reasonably good agreement except at the extreme low pressure conditions of the facility envelope.

CFD↗
NASA NTRS2015

Computational Simulations of the 10-MW TP3 Arc-Jet Facility Flow

This paper reports computational simulations and analysis in support of calibration tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using stagnation calorimeters and two different blunted wedge models with calibration plates at a wide range of conditions in the NASA Ames 10-megawatt TP3 (Test Position 3) facility. Data were obtained using four different conical nozzles with the same test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot and null-point heat flux probes were also performed at several arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. These computational simulations provide estimates of the arc-jet test environment parameters that are not measured but are needed to evaluate the performance of thermal protection system materials, along with further valuable insights into the arc-jet testing environment. Simulation results are used to estimate centerline total enthalpy, surface shear, boundary layer thickness, and boundary layer edge Mach number and to verify that specific test requirements from the Orion program are met.

CFD↗
NASA NTRS2019

CFD Simulations of the IHF Arc-Jet Flow: 9-Inch Nozzle, Flow Surveys, LEAF Wedge Calibration Data

This paper reports computational analyses and flow characterization studies in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using a wedge model placed in a free jet downstream of new 9-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. Both the nozzle and wedge model were specifically designed for testing in the new Laser-Enhanced Arc-jet Facility. Data were obtained using stagnation calorimeters and wedge models placed downstream of the nozzle exit. Two instrumented wedge calibration plates were used: one water-cooled and the other RCG-coated tile plate. Experimental surveys of arc-jet test flow with pitot and heat flux probes were also performed at three arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. By taking into account nonuniform total enthalpy and mass flux profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields, these simulations approximately reproduce the probe survey data and predict the wedge model surface pressure and heat flux measurements.

Arc-jets↗