SEARCH · Search NASA

Results for “time evolution operators”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Out-of-equilibrium Majorana zero modes in interacting Kitaev chains

Here, we employ a time-dependent real-space local density-of-states method to study the movement and fusion of Majorana zero modes in the one-dimensional interacting Kitaev model, based on the time evolution of many-body states. We analyze the dynamics and both fusion channels of Majoranas using time-dependent potentials, either creating walls or wells. For fast moving Majoranas, we unveil nonequilibrium signatures of the “strong-zero-mode” operator (quasiparity degeneracy in the full spectrum) and its breakdown in the presence of repulsive Coulomb interactions. Focusing on forming a full electron after fusion, we also discuss the upper and lower limits on the Majorana speed needed to reduce nonadiabatic effects and to avoid poisoning due to decoherence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Mean Flow from Phase Averages in the 2D Boussinesq Equations

The atmosphere and ocean are described by highly oscillatory PDEs that challenge both our understanding of their dynamics and their numerical approximation. This paper presents a preliminary numerical study of one type of phase averaging applied to mean flows in the 2D Boussinesq equations that also has application to numerical methods. The phase averaging technique, well-known in dynamical systems theory, relies on a mapping using the exponential operator, and then an averaging over the phase. The exponential operator has connections to the Craya–Herring basis pioneered by Jack Herring to study the fluid dynamics of oscillatory, nonlinear fluid dynamics. In this paper, we perform numerical experiments to study the effect of this averaging technique on the time evolution of the solution. We explore its potential as a definition for mean flows. We also show that, as expected from theory, the phase-averaging method can reduce the magnitude of the time rate of change in the PDEs, making them potentially suitable for time stepping methods.

Wingate, Beth A.↗

DOE OSTI2024

Operando Freezing Cryogenic Electron Microscopy of Active Battery Materials

Abstract Understanding structural and chemical evolution of battery materials during operation is critical to achieving safe, efficient, and long-lasting energy storage. Cryogenic electron microscopy (cryo-EM) has become a valuable tool in battery characterization, leveraging low temperatures to improve stability of sensitive materials under electron beam irradiation. However, typical cryo-EM sample preparations leave extended time between the electrochemical point of interest and ex situ freezing of samples, during which active structures may relax, degrade, or otherwise evolve. Here, we detail a method for operando freezing cryo-EM to preserve and characterize native electrode and interfacial structures that arise during battery cycling, based on an operando plunge freezer and cold sample removal process. We validate the method on multiple electrode materials and quantify and discuss the freezing rate achieved. Operando freezing cryo-EM can be used to directly visualize transient features that arise at active electrochemical interfaces, to enable deeper understanding of structural evolution and interfacial chemistry in batteries and other electrochemical systems.

25 ENERGY STORAGE↗

DOE OSTI2020

Evolution of the Data Quality Monitoring and Prompt Processing System in the proto DUNE-SP experiment

The DUNE Collaboration currently operates an experimental program based at CERN which includes a beam test and an extended cosmic ray run of two large-scale prototypes of the massive Liquid Argon Time Projection Chamber (LArTPC) for the DUNE Far Detector. The volume of data collected by the single-phase prototype (protoDUNE-SP) amounts to 3PB and the sustained rate of data sent to mass storage is O(100) MB/s. Data Quality Monitoring was implemented by directing a fraction of the data stream to the protoDUNE prompt processing system (p3s) which is optimized for continuous low-latency calculation of the vital detector metrics and various graphics including event displays. It served a crucial role throughout the life cycle of the experiment. We present our experience in leveraging the CERN computing environment and operating the system over an extended period of time, while adapting to evolving requirements and computing platform.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Forecasting for ESCAPE: A Multi-Institution Hybrid Forecasting and Nowcasting Operation for Sea-Breeze Convection Supporting a Ground-Based and Airborne Field Campaign

The Experiment of Sea-Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) field project deployed two aircraft and ground-based assets in the vicinity of Houston, Texas, between 27 May and 2 July 2022, examining how meteorological conditions, dynamics, and aerosols control the initiation, early growth stage, and evolution of coastal convective clouds. To ensure that airborne- and ground-based assets were deployed appropriately, a forecasting and nowcasting team was formed. Daily forecasts guided real-time decision-making by assessing synoptic weather conditions, environmental aerosol, and a variety of atmospheric modeling data to assign a probability for meeting specific ESCAPE campaign objectives. During the research flights, a small team of forecasters provided “nowcasting” support by analyzing radar, satellite, and new model data in real time. The nowcasting team proved invaluable to the campaign operation, as sometimes changing environmental conditions affected, for example, the timing of convective initiation. In addition to the success of the forecasting and nowcasting teams, the ESCAPE campaign offered a unique “testbed” opportunity where in-person and virtual support both contributed to campaign objectives. The forecasting and nowcasting teams were each composed of new and experienced forecasters alike, where new forecasters were given invaluable experience that would otherwise be difficult to attain. Both teams received training on forecast models, map analysis, Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT), and thermodynamic sounding analysis before the beginning of the campaign. In this article, the ESCAPE forecasting and nowcasting teams reflect on these experiences, providing potentially useful advice for future field campaigns requiring forecasting and nowcasting support in a hybrid virtual/in-person framework.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Data-Driven Learning for the Mori--Zwanzig Formalism: A Generalization of the Koopman Learning Framework

A theoretical framework which unifies the conventional Mori--Zwanzig formalism and the approximate Koopman learning of deterministic dynamical systems from noiseless observation is presented. In this framework, the Mori--Zwanzig formalism, developed in statistical mechanics to tackle the hard problem of construction of reduced-order dynamics for high-dimensional dynamical systems, can be considered as a natural generalization of the Koopman description of the dynamical system. We next show that, similar to the approximate Koopman learning methods, data-driven methods can be developed for the Mori--Zwanzig formalism with Mori's linear projection operator. We have developed two algorithms to extract the key operators, the Markov and the memory kernel, using time series of a reduced set of observables in a dynamical system. We have adopted the Lorenz `96 system as a test problem and solved for the above operators. These operators exhibit complex behaviors, which are unlikely to be captured by traditional modeling approaches in Mori--Zwanzig analysis. The nontrivial generalized fluctuation-dissipation relationship, which relates the memory kernel with the two-time correlation statistics of the orthogonal dynamics, was numerically verified as a validation of the solved operators. Here we present numerical evidence that the generalized Langevin equation, a key construct in the Mori--Zwanzig formalism, is more advantageous in predicting the evolution of the reduced set of observables than the conventional approximate Koopman operators.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Hamiltonian simulation in the low-energy subspace

Abstract We study the problem of simulating the dynamics of spin systems when the initial state is supported on a subspace of low energy of a Hamiltonian H . This is a central problem in physics with vast applications in many-body systems and beyond, where the interesting physics takes place in the low-energy sector. We analyze error bounds induced by product formulas that approximate the evolution operator and show that these bounds depend on an effective low-energy norm of H . We find improvements over the best previous complexities of product formulas that apply to the general case, and these improvements are more significant for long evolution times that scale with the system size and/or small approximation errors. To obtain these improvements, we prove exponentially decaying upper bounds on the leakage to high-energy subspaces due to the product formula. Our results provide a path to a systematic study of Hamiltonian simulation at low energies, which will be required to push quantum simulation closer to reality.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Imaging solid–electrolyte interphase dynamics using operando reflection interference microscopy

The quality of the solid–electrolyte interphase is crucial for the performance of most battery chemistries, but its formation dynamics during operation are not well understood due to a lack of reliable operando characterization techniques. Herein, we report a dynamic, non-invasive, operando reflection interference microscope to enable the real-time imaging of the solid–electrolyte interphase during its formation and evolution processes with high sensitivity. The stratified structure of the solid–electrolyte interphase formed during four distinct steps includes the emergence of a permanent inner inorganic layer enriched in LiF, a transient assembly of an interfacial electrified double layer and a consequent emergence of a temporary outer organic-rich layer whose presence is reversible with electrochemical cycling. Reflection interference microscope imaging reveals an inverse correlation between the thicknesses of two interphasial subcomponents, implying that the permanent inorganic-rich inner layer dictates the organic-rich outer layer formation and lithium nucleation. In conclusion, the real-time visualization of solid–electrolyte interphase dynamics provides a powerful tool for the rational design of battery interphases.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2024

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Newberry SHR Demonstration Project – Bipartisan Infrastructure Law Enhanced Geothermal (EGS) Pilot Demonstration (Abstract)

This project will create an Engineered Geothermal System (EGS) comprising two or more wells drilled to a depth of 4.25 km into superhot rock (SHR) with a temperature of 425 °C at Newberry Volcano in Central Oregon. An EGS is a manufactured heat exchanger in which water is injected in a deep injection well, or injector, to extract heat from the hot rock at depth and steam is returned to the surface in a production well, or producer, to generate electricity. In this project, the SHR EGS will be made using new methods and technologies to stimulate and connect hydraulic and natural fractures to enable multiple flow pathways between wells, allowing for optimal heat mining from the reservoir rock. The new technologies are designed to operate at rock temperatures much higher than those encountered in traditional geothermal. Following EGS completion, water will be injected into the injector well and steam extracted from the producer well in a long-term connectivity flow test demonstrating SHR reservoir evolution with time and use. Success will be measured by demonstrating the efficacy of new technologies and by producing economic quantities of steam (>40 MWth).

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Scrambling Dynamics with Imperfections in a Solvable Model

We study how probes of quantum scrambling dynamics respond to two kinds of imperfections -- unequal forward and backward evolutions and decoherence -- in a solvable Brownian circuit model. We calculate a ``renormalized'' out-of-time-order correlator (ROTOC) in the model with $N$ qubits, and we show that the circuit-averaged ROTOC is controlled by an effective probability distribution in operator weight space which obeys a system of $N$ non-linear equations of motion. These equations can be easily solved numerically for large system sizes which are beyond the reach of exact methods. Moreover, for an operator initially concentrated on weight one $w_0=1$, we provide an exact solution to the equations in the thermodynamic limit of many qubits that is valid for all times, all non-vanishing perturbation strengths $p\gtrsim 1/\sqrt{N}$, and all decoherence strengths. We also show that a generic initial condition $w_0 >1$ leads to a metastable state that eventually collapses to the $w_0=1$ case after a lifetime $\sim \log(N/w_0)$. Our results highlight situations where it is still possible to extract the unperturbed chaos exponent even in the presence of imperfections, and we comment on the applications of our results to existing experiments with nuclear spins and to future scrambling experiments.

FOS: Physical sciences↗

DOE OSTI2025

Risk assessment of wellbore leakage during underground hydrogen storage

The expansion of renewable energy sources would require large-scale energy storage options to overcome the intermittent nature of these sources. Underground hydrogen storage (UHS) in depleted hydrocarbon reservoirs offers a scalable and practical energy storage solution. These reservoirs are chosen for their availability and large capacity, but the unique properties of hydrogen raise concerns about potential leakage pathways, particularly through wellbores. In this study, we develop and apply, for the first time, reduced-order models (ROMs) specifically designed for efficient leakage risk prediction in UHS systems operating in depleted hydrocarbon reservoirs. Using 3,000 high-fidelity simulation scenarios, we examine the influence of 11 key parameters, including reservoir and aquifer depths, wellbore permeability and porosity, initial saturations of water, oil and gas fractions (hydrogen, light, intermediate, and heavy hydrocarbons), reservoir pressure multiplier, and the aquifer-to-reservoir volume ratio, to simulate leakage behavior over a 1,000-year timescale. We train ROMs using a two-step classification-regression approach, achieving R 2 values exceeding 99 % across all targets. These ROMs effectively capture the leakage evolution and identify critical controls of leakage, guiding the design of mitigation strategies. Results indicate that gas leakage occurs in about 27 % of scenarios as early as five years post-operation, reaching volumes of up to 106 ft3. Oil leakage is less frequent (~17 %) and typically begins decades later. Our findings also show that hydrogen often migrates first, owing to its smaller molecular size and higher buoyancy, followed by heavier hydrocarbons. Over time, these heavier components contribute significantly to the total leaked volume, reinforcing the need for targeted monitoring and remediation strategies. Our analysis highlights that deeper storage reservoirs, shallower aquifers, and low-permeability wellbores significantly reduce leakage risks. In conclusion, this work offers a robust framework for risk-informed UHS deployment, supporting energy security through reliable large-scale hydrogen storage while safeguarding environmental integrity.

08 HYDROGEN↗

DOE OSTI2022

Evolution of the marker distribution in gyrokinetic $δf$ particle-in-cell simulations

The evolution of the particle weight in a δf particle-in-cell simulation depends on the marker distribution that can evolve in a turbulent field due to turbulent diffusion. When Monte Carlo methods are used to implement the test particle collision operator, or when the particle motion is not strictly Hamiltonian in a collisionless simulation, the marker distribution will evolve along the particle trajectory and, in general, cannot be known exactly. A two-dimensional numerical marker distribution is proposed as an approximation. It is shown to be advantageous over other common methods for evaluating the marker distribution in long-time turbulence simulations. A generalized two-weight δf-method is proposed to mitigate the marker evolution problem.

, Monte Carlo methods↗

DOE OSTI2020

Strategies for thorium fuel cycle transition in the SD-TMSR

Liquid-fueled Molten Salt Reactor (MSR) systems represent advances in safety, economics, and sustainability. The MSR has been designed to operate with a Th/ 233 U fuel cycle with 233 U used as startup fissile material. Since 233 U does not exist in nature, we must examine other available fissile materials to start up these reactor concepts. This work investigates the fuel cycle and neutronics performance of the Single-fluid Double-zone Thorium-based Molten Salt Reactor (SD-TMSR) with different fissile material loadings at startup: High Assay Low Enriched Uranium (HALEU) (19.79%), Pu mixed with HALEU (19.79%), reactor-grade Pu (a mixture of Pu isotopes chemically extracted from Pressurized Water Reactor (PWR) spent nuclear fuel (SNF) with 33 GW d/tHM burnup), transuranic elements (TRU) from Light Water Reactor (LWR) SNF, and 233 U. The MSR burnup routine provided by SERPENT-2 is used to simulate the online reprocessing and refueling in the SD-TMSR. The effective multiplication factor, fuel salt composition evolution, and net production of 233 U are studied in the present work. Additionally, the neutron spectrum shift during the reactor operation is calculated. The results show that the continuous flow of reactorgrade Pu helps transition to the thorium fuel cycle within a relatively short time (≈ 4.5 years) compared to 26 years for 233 U startup fuel. Finally, using TRU as the initial fuel materials offers the possibility of operating the SD-TMSR for an extended period of time (≈ 40 years) without any external feed of 233 U.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Monitoring the long-term performance of organic redox flow battery by a distribution of relaxation time analysis

Organic redox flow batteries hold great promise as an energy storage technology, but their intricate chemistry makes them vulnerable to various degradation mechanisms. Monitoring this degradation is essential for identifying the limiting processes within the cells. Electrochemical impedance spectroscopy (EIS) offers a straightforward, in-situ method for measuring the total resistance of an operating cell. However, to pinpoint the limiting processes during long-term cycling, EIS data must be complemented by other techniques. Distribution of relaxation time (DRT) analysis is particularly effective for differentiating resistance components. Here, in this study, we perform a comprehensive analysis of resistance evolution and the separation of anode and cathode contributions during long-term cycling of a full cell employing 7,8-dihydroxyphenazine-2-sulfonic acid (DHPS) as the anolyte. Separate analyses of the DHPS anolyte and ferri-/ferrocyanide catholyte were conducted using a symmetric cell setup. The relaxation times derived from symmetric cells facilitate the identification of peaks in the DRT profiles from the full cell. Importantly, the DRT profiles indicate a correlation between the evolution of charge transfer resistance and the chemical degradation of DHPS. The methodologies and results outlined in this study offer significant insights for developing diagnostic tools applicable to other types of redox flow batteries.

Distribution of relaxation time↗

DOE OSTI2021

Probing Microstructure-induced Swelling & Thermal Property Changes in Ion-Irradiated Oxide Fuels using Laser-generated Surface Acoustic Waves

Radiation of fuel rods with neutrons during a fission reaction in a nuclear reactor is known to affect the material microstructure of the fuel as a result of fission fragment damage, and lead to the formation of atomic-level defects such as voids, dislocation loops, and lattice swelling from fission-gas release. A consequence of the radiation-induced damage in nuclear fuels is the drastic alteration of material properties, in particular, the reduction of the thermal conductivity – a key parameter that governs the transport of thermal energy released from fissile fuel atoms to the surrounding coolant. Swelling can induce high stresses and ultimate failure of the cladding. A fundamental understanding of the role of radiation-induced damage on the thermal transport and mechanical properties of nuclear fuels is therefore critical for efficient, reliable and safe operation of a nuclear power plant. As opposed to conventional post-irradiation examination techniques that are often time-consuming and require significant sample preparation, laser-based characterization methods have emerged as promising non-contact, non-destructive tools to measure the evolution of microstructure-induced material property changes. Moreover, laser methods offer the promise of in-situ characterization while the material is being irradiated. This project aims to develop an improved understanding of the impact of radiation-induced changes to material microstructure on the thermal and elastic properties of nuclear materials. A laser ultrasonic technique, known as the transient grating (TG) spectroscopy method, is used to simultaneously measure thermal diffusivity and elastic properties in unirradiated and ion-irradiated oxide nuclear fuel samples. The results will complement on-going investigations on electron- and phonon-mediated thermal transport in nuclear materials and will provide foundational work for incorporating the influence of defects in fuel performance codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Accelerator challenges in the Fermilab Neutrino Program

In Fermilab's 2022-2023 proton complex run, Fermilab recovered from an unexpected failure in the neutrino focusing horn to achieve a record beam power of 960 kW. For Snowmass and P5, Fermilab recently proposed a new upgrade plan known as Accelerator Complex Evolution (ACE) to occur after PIP-II upgrade in tandem with DUNE/LBNF operation. The first phase of ACE, is the ACE Main Injector Ramp and Targetry upgrade, would raise the LBNF beam power from 1.2 MW to 2 MW by shortening the Main Injector cycle time to around 0.65s. The second phase of ACE, the ACE Booster Replacement, would replace the Booster with a new machine and achieve 2.4 MW for DUNE/LBNF. The core accelerator challenges, present and future, for operating intense neutrino beams at Fermilab will be highlighted.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗