SEARCH · Search NASA

Results for “limit cycle”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Limit Cycle Oscillating Detonation Wave Behavior Analysis Within a Rotating Detonation Engine

Limit cycle oscillation (LCO) detonation wave behaviors are presented and analyzed for test times exceeding 20 s in a water-cooled rotating detonation engine (RDE). LCO detonation waves exhibit cyclic acceleration and deceleration, resulting in oscillating wave spacing at unique process conditions. In previous RDE studies, similar behaviors have been studied as microsecond-scale instabilities leading to ascending or descending modal transitions. In the current work, however, LCO waves are considered a persistent wave mode, occupying unique portions of the operational envelope adjacent to those of their equally spaced counterparts. These occurrences of LCO waves are repeatable, enduring behaviors. A method to generate shifted contour surfaces specifically intended to extract and analyze wave spacing variation through time, termed limit cycle oscillation visualization (LCOV) surfaces, is presented. LCOV surfaces transform data into the reference frame of a primary traveling wave and are used to analyze quasi-steady, short-timescale, and transitional LCO modes. Results are leveraged to understand the relationship between fill height, wave strength, local wave acceleration, and subsequent LCO wave spacing for individual wave sets. In conclusion, quasi-steady LCO waves display wave spacing oscillations between equal spacing values associated with ±1 wave across runs exceeding 18 s.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

A Model for Pair Production Limit Cycles in Pulsar Magnetospheres

Abstract It was recently proposed that the electric field oscillation as a result of self-consistent e ± pair production may be the source of coherent radio emission from pulsars. Direct particle-in-cell simulations of this process have shown that the screening of the parallel electric field by this pair cascade manifests as a limit cycle, as the parallel electric field is recurrently induced when pairs produced in the cascade escape from the gap region. In this work, we develop a simplified time-dependent kinetic model of e ± pair cascades in pulsar magnetospheres that can reproduce the limit-cycle behavior of pair production and electric field screening. This model includes the effects of a magnetospheric current, the escape of e ± , as well as the dynamic dependence of pair production rate on the plasma density and energy. Using this simple theoretical model, we show that the power spectrum of electric field oscillations averaged over many limit cycles is compatible with the observed pulsar radio spectrum.

Astronomy & Astrophysics↗

DOE OSTI2024

Limit cycle oscillations in the zonal-flow-catalyzed interactions of ion-temperature-gradient turbulence

Limit-cycle oscillations are studied for ion temperature gradient turbulence, which, in the absence of large diamagnetic (mean) shear flows, saturates through energy transfer from unstable modes to large-scale stable modes via zonal-flow intermediary modes. Oscillations of zonal flow and turbulence levels are strongly constrained by the reactive, largely non-dissipative character of the zonal flows. Since existing predator–prey models for observed oscillations in experiments do not include energy transfer through zonal flows to stable modes, low-order fluid models with this physics are constructed and investigated. A simple three-wave truncation produces low-amplitude zonal flows that slowly oscillate around a zero mean, with turbulence oscillations between coupled wavenumbers that exceed linear frequencies by orders of magnitude. This inconsistency with experimental observations is caused by the weak non-linear drive of zonal flows in three-wave systems and the lack of multiple-wavenumber turbulent interactions. A more comprehensive model that preserves multiple wavenumber interactions within the context of conservative zonal-flow-mediated energy transfer to stable modes accurately reflects observed dynamics when the phase between stable and unstable modes is occasionally randomized.

Li, P. -Y. (ORCID:0000000295254171)↗

DOE OSTI2026

Auxiliary-state-facilitated phase synchronization phenomena in isolated spin systems

Extending classical synchronization to the quantum domain is of great interest both from the fundamental physics point of view and with a view toward quantum technology applications. This work characterizes phase synchronization of an effective spin-1 system, which is realized by coupling three quantum states with infinite lifetime to auxiliary excited states that have a finite lifetime. Integrating out the excited states, the effective spin-1 model features coherent and incoherent effective couplings. The following are our key findings. (i) Phase synchronization can be controlled by adjusting the phases of the couplings to the excited states. (ii) Unlike in the paradigmatic spin-1 system studied in the literature, where the dissipative couplings describe decay into the limit-cycle state, the effective spin-1 model investigated in this work is governed by a competition between dissipative decay into and out of the limit-cycle state, with the dissipative decay out of the limit-cycle state playing a critical role. (iii) We identify a parameter regime where phase synchronization of the effective spin-1 system is, in the absence of coherent effective couplings, governed entirely by the effective dissipators. The effective spin-1 model is benchmarked through comparisons with master-equation calculations for the full Hilbert space. Physical insights are gained through analytical perturbation theory calculations. In conclusion, our findings, which are expected to hold for a broad class of energy-level and coupling schemes, are examined using hyperfine states of 87 Rb as an example system.

Nonlinear optics↗

DOE OSTI2024

Catalytic Resonance Theory: Forecasting the Flow of Programmable Catalytic Loops

Chemical transformations on catalyst surfaces occur through series and parallel reaction pathways. These complex networks and their behavior can be most simply evaluated through a three-species surface reaction loop (A* to B* to C* to A*) that is internal to the overall chemical reaction. Application of an oscillating dynamic catalyst to this reactive loop has been shown to exhibit one of three types of behavior: (1) a positive net flux of molecules about the loop in the clockwise direction, (2) a negative net flux of molecules about the loop in the counterclockwise direction, or (3) negligible flux of molecules about the loop at the limit cycle of reaction. Three-species surface loops were simulated with microkinetic modeling to assess the reaction loop behavior resulting from a catalytic surface oscillating between two or more catalyst surface energy states. Selected input parameters for the simulations spanned an 11-dimensional parameter space using 127 688 different parameter combinations. Their converged limit cycle solutions were analyzed for their loop turnover frequencies, the majority of which were found to be approximately zero. Classification and regression machine learning models were trained to predict the sign and magnitude of the loop turnover frequency and successfully performed above accessible baselines. Notably, the classification models exhibited a baseline weighted F1 score of 0.49, whereas trained models achieved weighted F1 scores of 0.94 and 0.96 when trained on the parameters used to define the simulations and derived rate constants, respectively. The trained models successfully predicted catalytic loop behavior, and interpretation of these models revealed all input parameters to be important for the prediction and performance of each model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Understanding LIB Battery Electrodes Through Classical Electrochemical Interface Theory

Classical models of the electrochemical interface can be applied to complex and dynamic electrodes to understand important mechanisms that are relevant to device-level phenomena like calendar aging or cycle life. Here, we show that negative alloy electrodes with wide electrochemical windows cannot be assumed to have static interfaces throughout all states of charge. Under high states of charge, the interface -including the solid electrolyte interphase - has characteristics that resemble a highly polarized electrode far from the point of zero charge. Under these conditions, the interface can be understood through the Helmholtz model. At lower states of charge, the interface capacitance decreases and the space charge layer length increases. This transformation resembles a metal electrode approaching the point of zero charge which is understood through the Gouy-Chapman-Stern-Grahame model of the electrochemical interface. This transformation happens at both silicon and carbon-coated silicon interfaces. From voltage-limited cycling experiments, in the limit of the diffuse double layer, the silicon electrode impedance rises significantly in the first few cycles. This rise is related to the continuous electrochemical reduction of electrolyte components from an interface that is not electronically screened from the electrolyte. In other words, an electrode interface for batteries should resemble a 'Helmholtz-like' structure.

anode↗

DOE OSTI2026

Impact of Lithium‐Free Borate Additives on the Cycle Life and Calendar Aging of Silicon‐Based Lithium‐Ion Batteries

Silicon-anode lithium-ion batteries (LIBs) suffer from limited cycle life and poor calendar life, constraining their large-scale commercialization. Integrating additives into electrolytes is a simple and cost-effective strategy to improve these aspects. The effects of lithium-free boron-based additives on cycling and calendar performance of high-loading Si-anode LIBs remain largely unexplored. In this work, the influence of five Li-free borate additives, each with distinct molecular structures and elemental compositions, is systematically investigated. All additives enhance cycle life to varying extents. Notably, the addition of 1 v/v% tri(2,2,2-trifluoroethyl) borate to the baseline electrolyte nearly doubles the cycle life at 50% state of health. This enhancement is attributed to three key factors. Specifically, borate additives 1) improve electrochemical activity, 2) act as anion receptors that interact with [PF6]- anions and carbonate solvents to reduce electrolyte decomposition, and 3) promote the formation of a stable and polymeric solid electrolyte interphase layer. Furthermore, these additives exhibited negligible impact in mitigating leakage current during a 180 h voltage-hold calendar-aging test, indicating their limited effect in calendar life. These findings provide insight into the role of Li-free borate additives in improving cycle life while addressing the knowledge gap regarding their influence on calendar aging.

Li, Defu↗

DOE OSTI2026

Observation of Synchronization between Two Quantum van der Pol Oscillators in Trapped Ions

Synchronization is a hallmark of collective behavior that emerges when nonlinear systems interact, spanning scales from mechanical oscillators to planetary orbits. As a universal phenomenon, it underpins the study of complex systems and has far-reaching technological implications. While classical synchronization has a long and rich history, it has not been observed experimentally between multiple quantum limit-cycle oscillators despite a decade of theoretical investigations. We realize synchronization between two quantum van der Pol oscillators by engineering dissipation in a mixed-isotope trapped-ion quantum simulator. The synchronized state is encoded in a fixed relative phase between the oscillators that is inaccessible to individual measurements and revealed only through joint readout of both oscillators, in stark contrast to the system in the (deterministic) classical limit where synchronization can be observed via individual phase measurements. We further show that the relative phase can be precisely controlled and that the chain of two oscillators can synchronize to an external field, suggesting applications in sensing. Our results provide a promising pathway for studying more complex synchronized quantum dynamics beyond two oscillators, where a theoretical treatment becomes increasingly challenging, and it remains to be understood whether genuinely quantum features persist in such cases.

Liu, Jiarui↗

DOE OSTI2024

Designing electrolytes with high solubility of sulfides/disulfides for high-energy-density and low-cost K-Na/S batteries

Alkaline metal sulfur (AMS) batteries offer a promising solution for grid-level energy storage due to their low cost and long cycle life. However, the formation of solid compounds such as M 2 S 2 and M 2 S (M = Na, K) during cycling limits their performance. Here we unveil intermediate-temperature K-Na/S batteries utilizing advanced electrolytes that dissolve all polysulfides and sulfides (K 2 S x , x = 1–8), significantly enhancing reaction kinetics, specific capacity, and energy density. These batteries achieve near-theoretical capacity (1655 mAh g -1 sulfur) at 75 °C with a 1 M sulfur concentration. At a 4 M sulfur concentration, they deliver 830 mAh g -1 at 2 mA cm -2 , retaining 71% capacity after 1000 cycles. This new K-Na/S battery with specific energy of 150-250 Wh kg -1 only employs earth-abundant elements, making it attractive for long-duration energy storage.

36 MATERIALS SCIENCE↗

DOE OSTI2024

In Situ Alloying Enabled by Active Liquid Metal Filler for Self‐Healing Composite Polymer Electrolytes

Abstract Solid inorganics, known for kinetically inhibiting polymer crystallization and enhancing ionic conductivity, have attracted significant attention in solid polymer electrolytes. However, current composite polymer electrolytes (CPEs) are still facing challenges in Li metal batteries, falling short of inhibiting severe dendritic growth and resulting in very limited cycling life. This study introduces Ga 62.5 In 21.5 Sn 16 (Galinstan) liquid metal (LM) as an active liquid alternative to conventional passive solid fillers, aiming at realizing self‐healing protection against dendrite problems. Compared to solid inorganics, for example silica, LM droplets could more significantly reduce polymer crystallinity and enhance Li‐ion conductivity due to their liquid nature, especially at temperatures below the polymer melting point. More importantly, LMs are unraveled as dynamic chemical traps, which are capable of blocking and consuming lithium dendrites upon contact via in situ alloying during battery operation and further inhibiting dendritic growth due to the lower deposition energy barrier of the formed Li‐LM alloy. As a proof of concept, by strategically designing an asymmetric CPE with the active LM filling, a solid‐state Li/LiFePO 4 battery achieves promising full‐cell functionality with notable rate performance and stable cycle life. This active filler‐mediated self‐healing approach could bring new insights into the battery design in versatile solid‐state systems.

Wu, Kai↗

DOE OSTI2024

Dynamics of hydrogen–ammonia–natural gas lean-premixed high-pressure flames

The influence of fuel composition on self-excited combustion instabilities in a high-pressure combustor operated with ammonia, hydrogen, and natural gas fuels is characterized with high-frequency pressure measurements and imaging of the flame structure. A micromix multi-stage injector with 19 elements is used to introduce the fuel blend premixed with heated air into an optically accessible combustor operated at ~1.1 MPa. As hydrogen is substituted for natural gas, longitudinal thermoacoustic instabilities are observed in the combustor with pressure fluctuation amplitudes as large as 8% of the mean chamber value. In the absence of natural gas, limit-cycle instability magnitude is largely insensitive to ammonia addition. However, fuel compositions with > 20% natural gas result in positive correlation between instability amplitude and ammonia concentration. Spatial distribution of heat release in the combustor evaluated from OH* chemiluminescence imaging reveals axial growth in the regions of heat release fluctuation as hydrogen decreases, thus correlating flame length to instability amplitude. Distinct transitions from the first harmonic of the fundamental longitudinal acoustic mode (~1100 Hz) to the fundamental mode (~550 Hz) are observed with root mean square pressure fluctuations exceeding 7% of the mean chamber value. Here, these cases correspond to hydrogen mole fractions ≤ 50%. Analysis of the phase relationship between pressure fluctuations upstream and downstream of the flame zone indicates acoustic coupling of the injector as a key contributor to instability growth.

03 NATURAL GAS↗

DOE OSTI2025

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗

DOE OSTI2024

Dynamic gain and frequency comb formation in exceptional-point lasers

Abstract Exceptional points (EPs)—singularities in the parameter space of non-Hermitian systems where two nearby eigenmodes coalesce—feature unique properties with applications such as sensitivity enhancement and chiral emission. Existing realizations of EP lasers operate with static populations in the gain medium. By analyzing the full-wave Maxwell–Bloch equations, here we show that in a laser operating sufficiently close to an EP, the nonlinear gain will spontaneously induce a multi-spectral multi-modal instability above a pump threshold, which initiates an oscillating population inversion and generates a frequency comb. The efficiency of comb generation is enhanced by both the spectral degeneracy and the spatial coalescence of modes near an EP. Such an “EP comb” has a widely tunable repetition rate, self-starts without external modulators or a continuous-wave pump, and can be realized with an ultra-compact footprint. We develop an exact solution of the Maxwell–Bloch equations with an oscillating inversion, describing all spatiotemporal properties of the EP comb as a limit cycle. We numerically illustrate this phenomenon in a 5-μm-long gain-loss coupled AlGaAs cavity and adjust the EP comb repetition rate from 20 to 27 GHz. This work provides a rigorous spatiotemporal description of the rich laser behaviors that arise from the interplay between the non-Hermiticity, nonlinearity, and dynamics of a gain medium.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Validation of a Hybrid Domain Overlapping Coupling Between SAM and CFD Against the TALL-3D Transients

The System Thermal Hydraulics (STH) code SAM has been coupled to the Computational Fluid Dynamics (CFD) code Simcenter STAR-CCM+ utilizing a hybrid domain overlapping method with an explicit coupling in time. The coupling aims to extend the STH code’s applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional (3D) mixing. The coupling method’s numerical stability was verified in the past against two closed-loop configurations, and it was validated against a double T-junction experiment with 3D scalar mixing. In the present work, the coupling method is validated against the TALL-3D STH/CFD coupling benchmark facility. TALL-3D is a three-legged, liquid-metal facility with a large, pool-type enclosure (test section) that exhibits 3D flow effects to be modeled by a CFD code. The rest of the system exhibits approximately 1D behavior well-predicted by an STH code. First, the present STAR-CCM+ CFD model of the 3D test section is validated against experimental data. Then, the SAM-STARCCM+ coupled model is validated against six different TALL-3D steady states, including SAM standalone model results for comparison. Lastly, the SAM-STARCCM+ coupled model is validated against two TALL-3D transients, one exhibiting flow reversal in the test section and one exhibiting nonlinear, Limit Cycle Oscillations (LCO). For the first transient, the SAM-STARCCM+ coupled model properly predicts an increase in the test section’s inlet temperature during flow reversal, and this is not predicted by the SAM standalone model. Following flow reversal, the SAM-STARCCM+ coupled model better-predicts the initial flow recovery and following oscillations as the system approaches a final natural circulation state. For the second transient, no true final steady state is observed due to LCO. Neither the SAM-STARCCM+ coupled model nor the SAM standalone model can perfectly capture the experiment’s changing oscillation frequency during the transient. However, the SAM-STARCCM+ coupled model does reproduce the oscillatory feedback observed in the system. This is a significant achievement as the SAM-STARCCM+ coupled model only uses an explicit coupling in time, as opposed to a semi-implicit coupling. In comparison, previous STH/CFD coupling efforts of the TALL-3D facility required semi-implicit coupling to obtain similar results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Minimal cyclic behavior in sheared amorphous solids

Although jammed packings of soft spheres exist in potential energy landscapes with a vast number of minima, when subjected to cyclic shear they may revisit the same configurations repeatedly. Simple hysteretic spin models, in which particle rearrangements are represented by interacting spin flips called hysterons, capture many features of this periodic behavior. Yet it has been unclear to what extent individual rearrangements can be described by such binary objects and how such objects interact with one another. Using a particularly sensitive algorithm, we identify rearrangements in simulated jammed packings and select pairs of rearrangements that undo one another to create periodic cyclic behavior. We find that the rearrangement pairs surprisingly persist down to the smallest increments in strain, even in the smallest systems we can study. We explore the statistics of these rearrangement pairs and find that there is a relation between the amount of hysteresis and the energy drop and mean-square displacement of the particles; these results are inconsistent with the scaling found in models that treat rearrangements as localized buckling events. Finally, our analysis shows that there is no clean distinction between the particle motions that represent the identity of a single, individual rearrangement and the particle motions that lead to interactions between separated rearrangements or hysterons. These results offer insight into how complex systems such as amorphous solids can reach a limit cycle.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗

DOE OSTI2024

Physics-Based Continuum Modeling for an Aqueous Rechargeable Zn/MnO 2 Battery

This study introduces a framework for modeling the aqueous Zn/MnO 2 rechargeable battery. A reaction system and a physics-based continuum model are proposed based on two reaction types, one involving insertion and the second related to dissolution and deposition of a solid reaction product. The model, fitted to empirical data, predicts voltage behavior and capacity limitations during cycling, identifying electrolytic zinc depletion as a limiting mechanism, depending on the original cell construction. The research suggests the need for further material characterization and reaction analysis, which will advance our understanding and facilitate the development of grid-scale energy storage solutions.

Electrochemistry↗

DOE OSTI2025

High-Voltage, Intermediate-Temperature, Fe- and Al-Mixed Metal Halide Molten Salt for Molten Sodium Battery Energy Storage

An inorganic Fe and Al halide-based, low-temperature molten salt catholyte is described, which, when paired with a molten sodium anode, has high operating potentials rivaling those of Li ion batteries. The newly developed catholyte consists of metal halides FeCl 3 /FeCl 2 –AlCl 3 –NaCl and is intended to cycle between Fe 3+ /Fe 2+ redox couples in the molten salt. The multicomponent molten salt was initially evaluated for phase behavior and basic electrochemical behavior before full battery testing. The assembled battery, utilizing a 20:35:45 (FeCl 3 :AlCl 3 :NaCl) composition, with a 50.83 Ah/kg theoretical gravimetric capacity and a specific energy of 176.95 Wh/kg, was cycled at variable depths of discharge (DoD) and current densities to determine its cycling efficiencies and limitations. In conclusion, preliminary cycling tests showed two operational potential regimes, with higher potential, 3.91 V (vs Na/Na + ), at low DoD and lower potential, 3.39 V (vs Na/Na + ), at high DoD with excellent energy efficiencies and cycling behavior under both regimes.

Aluminum↗