Search NASA⌕ Search

SEARCH · Search NASA

Results for “STABLE OSCILLATION”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Lack of robustness and accuracy of many numerical schemes for phase-field simulations

In this paper, we study the stability, accuracy and convergence behavior of various numerical schemes for phase-field modeling through a simple ODE model. Both theoretical analysis and numerical experiments are carried out on this ODE model to demonstrate the limitation of most numerical schemes that have been used in practice. One main conclusion is that the first-order fully implicit scheme is the only robust algorithm for phase-field simulations while all other schemes (that have been analyzed) may have convergence issue if the time step size is not exceedingly small. More specifically, by rigorous analysis in most cases, we have the following conclusions: (i) The first-order fully implicit scheme converges to the correct steady state solution for all time step sizes. In the case of multiple solutions, one of the solution branches always converges to the correct steady state solution. (ii) The first-order convex splitting scheme, which is equivalent to the first-order fully implicit scheme with a different time scaling, always converges to the correct steady state solution but may seriously lack numerical accuracy for transient solutions. (iii) For the second-order fully implicit and convex splitting schemes, for any time step size $δt > 0$, there exists an initial condition $u_0$ , with $|u_0| > 1$, such that the numerical solution converges to the wrong steady state solution. (iv) For $|u_0|$ ≤ 1, all second-order schemes studied in this paper converge to the correct steady state solution although severe numerical oscillations occur for most of them if the time step size is not sufficiently small. (v) An unconditionally energy-stable scheme (such as the modified Crank–Nicolson scheme) is not necessarily better than a conditionally energy-stable scheme (such as the Crank–Nicolson scheme). Finally, most, if not all, of the above conclusions are expected to be true for more general Allen–Cahn and other phase-field models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Suppression of edge localized modes with real-time boron injection using the tungsten divertor in EAST

We report an observation of robust suppression of edge-localized modes (ELMs) in the Experimental Advanced Superconducting Tokamak (EAST), enabled by continuous boron (B) powder injection. Edge harmonic oscillations appear during B powder injection, providing sufficient particle transport to maintain constant density and avoid impurity accumulation in ELM-stable plasmas. Here, quasi-steady ELM suppression discharges are demonstrated with modest energy confinement improvement and over a wide range of conditions: heating power and technique variation, electron density range over a factor ~3.5, deuterium or helium ion species, and with either direction of the toroidal magnetic field. ELM suppression is observed above a threshold edge B intensity and ceases within 0.5 s of termination of the B injection. In contrast to ELM suppression accompanied by recycling reduction during Li powder injection in NSTX and EAST (Maingi et al 2018 Nucl. Fusion 58 024003), reduced recycling due to hydrogenic species retention is unnecessary for the ELM suppression with B powder injection, paving the way for its consideration as an ELM control tool for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2020

Data-driven Detection and Early Prediction of Thermoacoustic Instability in a Multi-nozzle Combustor

Thermoacoustic instability (TAI) is a critical issue in modern lean-burn gas-turbine combustors, which is induced by a strong coupling between the resonant combustor acoustics and fluctuations in the heat release rate. This instability may lead to high-amplitude pressure waves that generate undesirable noise levels as well as fatigue stresses in mechanical structures of the combustor. The intense pressure fluctuations due to TAI may also cause large flow perturbations and possibly flow reversal that may lead to flame oscillations, flame liftoff, and even flame blow-out. Hence, there is a strong need for exercising control actions in a timely fashion to mitigate the TAI phenomena. Anomaly detection is an essential prerequisite to the design of a good controller and such a detector must be able to reliably predict a forthcoming TAI. To detect and predict the onset of a TAI from an ensemble of pressure time series, this paper investigates three data-driven methods: Fast Fourier transform (FFT), symbolic time series analysis (STSA), and hidden Markov modeling (HMM). The main focus of the paper is to make a comparative evaluation of these three anomaly detection methods for classification of the current regime of operation into stable and unstable categories as well as for real-time identification of precursors to impending instabilities with short-length time series of measured variables (e.g., pressure oscillations). Finally, the results, generated on experimental data from a multi-nozzle combustor apparatus, have been compared to evaluate the performance of FFT, STSA, and HMM methods for TAI analysis.

42 ENGINEERING↗
DOE OSTI2022

Transient quantum beatings of trions in hybrid organic tri-iodine perovskite single crystal

Utilizing the spin degree of freedom of photoexcitations in hybrid organic inorganic perovskites for quantum information science applications has been recently proposed and explored. However, it is still unclear whether the stable photoexcitations in these compounds correspond to excitons, free/trapped electron-hole pairs, or charged exciton complexes such as trions. Here we investigate quantum beating oscillations in the picosecond time-resolved circularly polarized photoinduced reflection of single crystal methyl-ammonium tri-iodine perovskite (MAPbI 3 ) measured at cryogenic temperatures. We observe two quantum beating oscillations (fast and slow) whose frequencies increase linearly with B with slopes that depend on the crystal orientation with respect to the applied magnetic field. We assign the quantum beatings to positive and negative trions whose Landé g-factors are determined by those of the electron and hole, respectively, or by the carriers left behind after trion recombination. These are $g^{e}_{[001]}$ = 2.52 and $g^{e}_{[1\bar10]}$ = 2.63 for electrons, whereas |$g^{h}_{[001]}$| = 0.28 and |$g^{h}_{[1\bar10]}$| = 0.57 for holes. The obtained g-values are in excellent agreement with an 8-band K.P calculation for orthorhombic MAPbI3. Using the technique of resonant spin amplification of the quantum beatings we measure a relatively long spin coherence time of ~ 11 (6) nanoseconds for electrons (holes) at 4 K.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Förster Resonance Energy Transfer and Enhanced Emission in Cs 4 PbBr 6 Nanocrystals Encapsulated in Silicon Nano-Sheets for Perovskite Light Emitting Diode Applications

Encapsulating Cs 4 PbBr 6 quantum dots in silicon nano-sheets not only stabilizes the halide perovskite, but also takes advantage of the nano-sheet for a compatible integration with the traditional silicon semiconductor. Here, we report the preparation of un-passivated Cs 4 PbBr 6 ellipsoidal nanocrystals and pseudo-spherical quantum dots in silicon nano-sheets and their enhanced photoluminescence (PL). For a sample with low concentrations of quantum dots in silicon nano-sheets, the emission from Cs 4 PbBr 6 pseudo-spherical quantum dots is quenched and is dominated with Pb 2+ ion/silicene emission, which is very stable during the whole measurement period. For a high concentration of Cs 4 PbBr 6 ellipsoidal nanocrystals in silicon nano-sheets, we have observed Förster resonance energy transfer with up to 87% efficiency through the oscillation of two PL peaks when UV excitation switches between on and off, using recorded video and PL lifetime measurements. In an area of a non-uniform sample containing both ellipsoidal nanocrystals and pseudo-spherical quantum dots, where Pb 2+ ion/silicene emissions, broadband emissions from quantum dots, and bandgap edge emissions (515 nm) appear, the 515 nm peak intensity increases five times over 30 min of UV excitation, probably due to a photon recycling effect. This irradiated sample has been stable for one year of ambient storage. Cs 4 PbBr 6 quantum dots encapsulated in silicon nano-sheets can lead to applications of halide perovskite light emitting diodes (PeLEDs) and integration with traditional semiconductor materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Decentralized Distribution System Restoration with Grid-Forming/Following Inverter-Based Resources

The high penetration of distributed energy resources (DERs) in active distribution systems has posed challenges to the centralized distribution system restoration (DSR) strategies in current practice. On the other hand, the advancement in smart inverter technologies enables the bottom-up restoration capability. This paper is motivated to develop a 3-layered hierarchical framework for decentralized DSR, based on the grid-forming (GFM) and grid-following (GFL) grid-edge inverters. The first layer presents the tertiary control, which determines the load pickup schedule and generation dispatch of DERs, using the alternating direction method of the multipliers algorithm. The second layer consists of two control functions: GFM control, which regulates voltage and frequency, establishing a stable grid for GFL inverters to follow; and GFL control, which regulates the real and reactive power. In the third layer, the primary control is proposed to regulate the inverter voltage and current, which is developed based on the virtual oscillator control (VOC). Furthermore, the developed framework is tested in the modified IEEE 13-node test feeder. Two scenarios of grid-connected and islanded operating modes are designed, and simulation results demonstrate the effectiveness of decentralized DSR strategies for controlling grid-edge inverters to enhance the distribution system resilience.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Climate models project increasing precipitation in the US Southwest in 2024–2100

We analyzed the measured precipitation in the southwestern US and found that, from 1900 to 2024, precipitation decreased at an average rate of − 1.2 cm per year per 100 years. The most significant precipitation decreases occurred in the three-month period from February to April. At the same time, precipitation during the southwestern monsoon season (July, August, and September) remained relatively stable. The ensemble mean of all CMIP6 (Coupled Model Intercomparison Project phase 6) climate models, along with the regression model incorporating anthropogenic aerosol (AER) and Pacific Decadal Oscillation (PDO) as predictors, projects precipitation to increase from 2024 to 2100. This projected precipitation increase relies on the anticipated decrease in emissions of anthropogenic aerosols, which is associated with the transition from fossil fuel burning to renewable and nuclear energy sources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Investigation and validation of the dynamic response of an acoustically levitated particle using the lattice Boltzmann method

The stable levitation of an analyte sample in an acoustic levitator is a primary requirement for accurate x-ray characterization of its scientific structure. A rigid particle oscillates in an under-damped manner when introduced into the node of established standing acoustic waves. This investigation has employed the lattice Boltzmann method (LBM), a computational fluid dynamics technique, for the analysis of such rigid particle dynamics in acoustic levitation. The simulation uses the two dimensional and nine velocity (D2Q9) Bhatnagar-Gross-Krook formulation to levitate a rigid 1.6 mm diameter nylon ($\rho$ = 1150 kg/m 3 ) particle in the air at standard pressure and temperature conditions. The presented work is the first reported simulation of realistic acoustic levitator boundary conditions using the LBM. The simulation can capture the particle-fluid interactions that produce dynamic levitation at less than one-period timescale in the ultrasonic frequency regime. An experiment was conducted by levitating a 1.6mm-diameter nylon sphere to estimate the oscillations, and the oscillating frequency was found to be 50 Hz. The dynamic simulation results are consistent with experimental results for particle oscillations within the same order of magnitude, indicating that LBM formulation can be successfully used to study acoustic levitation to understand and mitigate particle jitter. The distortion of the acoustic field due to a levitating particle's presence was also analyzed to demonstrate how the presence of the particle can disrupt adjacent levitating nodes.

42 ENGINEERING↗
DOE OSTI2020

Properties of F Stars with Stable Radial Velocity Timeseries: A Useful Metric for Selecting Low-jitter F Stars

In a companion paper, we have conducted an in-depth analysis of radial velocity jitter of over 600 stars, examining the astrophysical origins including stellar granulation, oscillation, and magnetic activity. In this paper, we highlight a subsample of those stars, specifically the main sequence and “retired” F stars—which we refer to as “MSRF” stars—that show low levels of RV jitter (<10 m s{sup −1}). We describe the observational signatures of these stars that allow them to be identified in radial velocity planet programs, for instance, those performing follow-up of transiting planets discovered by TESS. We introduce a “jitter metric” that combines the two competing effects of RV jitter with age: activity and convection. Using thresholds in the jitter metric, we can select both “complete” and “pure” samples of low jitter F stars. We also provide recipes for identifying these stars using only Gaia colors and magnitudes. Finally, we describe a region in the Gaia color–magnitude diagram where low jitter F stars are most highly concentrated. By fitting a ninth-order polynomial to the Gaia main sequence, we use the height above the main sequence as a proxy for evolution, allowing for a crude selection of low jitter MSRF stars when activity measurements are otherwise unavailable.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Active electrochemical interfaces stabilized through self-organized potential oscillations

Some electrochemical systems are known to have higher average efficiency when operated under an oscillatory regime. Given the compromise between activity and stability, the stability of electrochemical interfaces in a self-organized, oscillatory state must be taken into account. Here we evaluate the electro-oxidation of methanol and formic acid on platinum under regular and oscillatory conditions, and study the stability by following the Pt dissolution rates in situ with a stationary probe rotating disk electrode (SPRDE) coupled to an inductively coupled plasma mass spectrometer (ICP-MS). Generally speaking, as the electro-oxidation reaction proceeds, the platinum dissolution rate increases considerably. To guarantee Pt stability, the potential must be kept below 1.0 V vs. RHE. Interestingly, no dissolution is detectable when the electrode potential undergoes temporary self-organization, ensuring a stable and active interface.

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

Enhanced Synchronization Stability of Grid-Forming Inverters With Passivity-Based Virtual Oscillator Control

In this article, a passivity-based virtual oscillator control strategy with enhanced synchronization stability for grid-forming inverters (GFMs) is proposed. By adopting the port-controlled Hamiltonian system theory for orbital stabilization problems, an energy pumping-and-damping block is proposed to render GFMs globally asymptotically stable with respect to the prespecified solutions of the power-flow equations from any initial condition. This allows for stable integrations of GFMs to any other globally asymptotically stable systems without their explicit knowledge, e.g., helping maintain synchronism with the bulk power system in a wide range of short-circuit-ratio conditions or under large disturbances and keeping synchronism among multiple GFMs in power systems. Herein, both simulations and experiments are presented to demonstrate the proposed control approach.

42 ENGINEERING↗
DOE OSTI2020

Changes of deformed shell gaps at N ~ 100 in light rare-earth, neutron-rich nuclei

There have been compelling evidences indicating that deformed single-particle states in neutron- rich regions are different from those in the stable region. It was pointed out that for the light rare- earth (60Nd, 62Sm, and 64Gd), neutron-rich (N = 98-102) nuclei, the Woods-Saxon potential, the Nilsson modified oscillator potential with "universal" parameters, and the folded Yukawa potential all failed to describe the correct ordering of neutron single-particle states. The location and size of deformed shell gaps in this mass region are under current debate. We propose a modification for the "standard" Nilsson parameters of Bengtsson and Ragnarsson introduced in 1985. The proposed N-dependent spin-orbit interaction brings an effect of the tensor force for the deformed potential, which causes dynamic changes in deformed shell gaps with neutron number and deformation. By applying the modified Nilsson parameters to generate deformed bases for the Projected Shell Model, we demonstrate that the calculation can explain consistently the current experimental data, including the ground state configuration in odd-neutron nuclei, upbending of the yrast moment of inertia at higher spins and the energies of 2-quasineutron 6- and 4- isomers in even-even nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Sustained edge-localized-modes suppression and radiative divertor with an impurity-driven instability in tokamak plasmas

We report simultaneous control of the large edge localized modes (ELMs) and divertor heat fluxes in a metal wall environment is a critical issue for steady-state operation of a tokamak fusion reactors. Here we report a sustained ELM suppression scenario achieved in the EAST tokamak compatible with radiative divertor using different seeding impurity species over a wide range of conditions. A low-n mode appears, as manifested by the oscillations of a radiation front near the X-point. This mode appears to drive strong particle transport and tungsten exhaust, which is essential to the maintenance of the ELM-stable state. We have developed a model to explain the mode excitation, by coupling the impurity radiative condensation instability to drift waves, which could explain some characteristics of the low-n mode well. The low-n mode may offer a new ELM-stable scenario compatible with radiative divertor for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Cosmological limits on the neutrino mass sum for beyond-Λ⁢ CDM models

The sum of neutrino masses can be measured cosmologically, as the sub-eV particles behave as “hot” dark matter whose main effect is to suppress the clustering of matter compared to a universe with the same amount of purely cold dark matter. Current astronomical data provide an upper limit on ∑𝑚 𝜈 between 0.07–0.12 eV at 95% confidence, depending on the choice of data. This bound assumes that the cosmological model is Λ Cold Dark Matter (Λ⁢ CDM), where dark energy is a cosmological constant, the spatial geometry is flat, and the primordial fluctuations follow a pure power law. Here, we update studies on how the mass limit degrades if we relax these assumptions. To existing data from the Planck satellite we add new gravitational lensing data from the Atacama Cosmology Telescope, the new Type Ia supernova sample from the Pantheon+survey, and baryonic acoustic oscillation (BAO) measurements from the Sloan Digital Sky Survey and the Dark Energy Spectroscopic Instrument. Using our fiducial data combination, described in the appendix, we find the neutrino mass limit is stable to most model extensions, with such extensions degrading the limit by less than 10%. We find a broadest bound of ∑𝑚 𝜈 < 0.19 eV at 95% confidence for a model with dynamical dark energy, although this scenario is not statistically preferred over the simpler Λ ⁢CDM model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Stability-preserving Lossy Compression for Large-scale Partial Differential Equations

Checkpoint/Restart (C/R) strategies are vital for fault tolerance in PDE-based scientific simulations, yet traditional checkpointing incurs significant I/O overhead. Lossy compression offers a scalable solution by reducing checkpoint data size, but conventional methods often lack control over physical invariants (e.g., energy), leading to instability such as oscillations or divergence in Partial Differential Equations (PDE) systems. This paper introduces a stability-preserving compression approach tailored for PDE simulations by explicitly controlling kinetic and potential energy perturbations to ensure stable restarts. Extensive experiments conducted across diverse PDE configurations demonstrate that our method maintains numerical stability with minimal error magnification—even across multiple checkpoint-restart cycles—outperforming state-of-the-art lossy compressors. Parallel evaluations on the Frontier supercomputer show up to 8.4× improvement in checkpoint write performance and 6.3× in read performance, while maintaining relative L2 errors ∼ 2e-6 throughout continued simulation. These results provide practical guidance for balancing compression accuracy, stability, and computational efficiency in large-scale PDE applications.

Gong, Qian [ORNL] (ORCID:0000000235704142)↗