Search NASA⌕ Search

SEARCH · Search NASA

Results for “Oscillatory”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Optimization-Based Fast-Frequency Estimation and Control of Low-Inertia Microgrids

The lack of inertial response from non-synchronous, inverter-based generation in microgrids makes the power system vulnerable to a large rate of change of frequency (ROCOF) and frequency excursions. Energy storage systems (ESSs) can be utilized to provide fast-frequency support to prevent such large excursions in the system. However, fast-frequency support is a power-intensive application that has a significant impact on the ESS lifetime. In this paper, a framework that allows the ESS operator to provide fast-frequency support as a service is proposed. The framework maintains the desired quality-of-service (limiting the ROCOF and frequency) while taking into account the ESS lifetime and physical limits. The framework utilizes moving horizon estimation (MHE) to estimate the frequency deviation and ROCOF from noisy phase-locked loop (PLL) measurements. These estimates are employed by a model predictive control (MPC) algorithm that computes control actions by solving a finite-horizon, online optimization problem. Additionally, this approach avoids oscillatory behavior induced by delays that are common when using low-pass filters as with traditional derivative-based (virtual inertia) controllers. MATLAB/Simulink simulations on a test system from Cordova, Alaska, show the effectiveness of the MHE-MPC approach to reduce frequency deviations and ROCOF of a low-inertia microgrid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Enabling Grid-Forming Control with Fault Ride-Through in Unbalanced Distribution Networks

Distribution networks are often unbalanced, causing oscillatory responses in inverter control designed for balanced conditions. Here, to address this problem, this paper proposes a novel time-domain transformation appropriate for inverter control and enables the decomposition of three-phase unbalanced signals into constant positive and negative components. Relations useful for calculating unbalanced active and reactive power are derived from first principle, providing insight into vector products of unbalanced three-phase signals. Furthermore, a grid-forming control effective under unbalanced conditions is developed, which delivers superior performance while meeting UNIFI1 specifications for grid-forming control under unbalanced conditions. specifications applicable to category 4 inverter-based resource, like setting and regulating frequency/voltage, providing voltage support, sharing active power, injecting negative sequence current, and riding through faults. A current limiter is proposed for safe fault ride-through and integrates with the grid-forming control featuring frequency/voltage droop controllers and current and voltage control loops. The transformation of interconnected inverters is formulated and stability of the proposed control analyzed to support robust parameter selections. The effectiveness of the proposed transformation and grid-forming control is demonstrated through analytical results and real-time simulation of a IEEE 123 distribution network on the Real-Time Digital Simulator. Comparison with existing methods shows that the proposed strategy satisfies the UNIFI specifications with a much better performance.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Local Communication in Populations of Thermostatically Controlled Loads

Thermostatically controlled loads (TCLs) have shown great potential for demand response (DR) in electric grid operations. However, it has been commonly seen that DR events using TCLs may cause load synchronization and unwanted oscillatory effects, especially in homogeneous populations. In an attempt to mitigate the negative impacts of DR events, a decentralized method is proposed that modifies each thermostat behavior based on the activity of a small number of nearby TCLs. This feedback introduces the possibility of instability in the aggregate behavior. A stability analysis is performed on a linearized model of the aggregate system and the results of that analysis compared to simulation results. The proposed modification of thermostat behavior results in fourfold reduction in the post-DR peak while suppressing ensuing oscillations at the expense of a modest increase in compressor cycling. Further, the linearized model also provides insight into the aggregate behavior of the population.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Influence of Tool Runout on Force Measurement During Internal Void Monitoring for Friction Stir Welding of 6061-T6 Aluminum

Abstract The goal of this research was to examine how altering the amount of friction stir tool eccentricity while controlling the amount of slant in the tool shoulder (drivers of oscillatory process forces) effects the generation of process force transients during sub-surface void interaction. The knowledge gained will help improve the accuracy of force-based void monitoring methods that have the potential to reduce the need for post-weld inspection. Process force transients during sub-surface void formation were examined for multiple tools with varying magnitudes of kinematic runout. The eccentric motion of the tool produced oscillations in the process forces at the tools rotational frequency that became distorted when features (flats) on the tool probe interacted with voided volumes, generating an amplitude in the force signals at three times the tool rotational frequency (for three-flat tools). A larger tool eccentricity generates a larger amplitude in the force signals at the tool’s rotational frequency that holds a larger potential to create a distortion during void interaction. It was determined that once void becomes large enough to produce an interaction that generates an amplitude at the third harmonic larger than 30% of the amplitude at the rotational frequency in a weld with no interaction (amplitude solely at rotational frequency), the trailing edge of the tool shoulder cannot fully consolidate the void, i.e., it will remain in the final weld. Additionally, once the void exceeds a certain size, the amplitudes of the third harmonics saturate at 70% of the amplitude at the rotational frequency during full consolidation. The interaction between the eccentric probe and sub-surface void was isolated by ensuring any geometric imperfection in the shoulder (slant) with respect to the rotational axis was removed. The results suggest that geometric imperfections (eccentricity and slant) with respect to the tool’s rotational axis must be known when developing a void monitoring method from force transients of this nature.

Engineering↗
DOE OSTI2021

Specular-reflection contributions to static and dynamic radiation forces on circular cylinders

Interest in the response of highly reflecting objects in water to modulated acoustical radiation forces makes it appropriate to consider contributions to such forces from perfectly reflecting objects to provide insight into radiation forces. The acoustic illumination can have wavelengths much smaller than the object's size, and objects of interest may have complicated shapes. Here, the specular contribution to the oscillating radiation force on an infinite circular cylinder at normal incidence is considered for double-sideband-suppressed carrier-modulated acoustic illumination. The oscillatory magnitude of the specular force decreases monotonically with increasing modulation frequency, and the phase of the oscillating force depends on the relative phase of the sidebands. The phase dependence on the modulation frequency can be reduced with the appropriate selection of a sideband relative-phase parameter. That is a consequence of the significance of rays that are incident on the cylinder having small impact parameters that are nearly backscattered. For one choice of a relative sideband phase, a prior partial wave series (PWS) solution is available, which supports the specular analysis when the PWS is evaluated for a rigid cylinder. The importance of specular contributions for aluminum cylinders in water is noted. A specular analysis for an analogous spherical reflector is also summarized.

Marston, Philip L. (ORCID:0000000306778537)↗
DOE OSTI2020

Medium amplitude parallel superposition (MAPS) rheology. Part 2: Experimental protocols and data analysis

An experimental protocol is developed to directly measure the new material functions revealed by medium amplitude parallel superposition (MAPS) rheology. This protocol measures the medium amplitude response of a material to a simple shear deformation composed of three sine waves at different frequencies, revealing a rich dataset consisting of up to 19 measurements of the third-order complex modulus at distinct three-frequency coordinates. We discuss how the choice of input frequencies influences the features of the MAPS domain studied by the experiment. A polynomial interpolation method for reducing the bias of measured values from spectral leakage and reducing variance due to noise is discussed, including a derivation of the optimal range of amplitudes for the input signal. This leads to the conclusion that conducting the experiment in a stress-controlled fashion possesses a distinct advantage to the strain-controlled mode. The experimental protocol is demonstrated through measurements of the MAPS response of a model complex fluid: a surfactant solution of wormlike micelles. The resulting dataset is indeed large and feature-rich, while still acquired in a time comparable to similar medium amplitude oscillatory shear (MAOS) experiments. We demonstrate that the data represent measurements of an intrinsic material function by studying its internal consistency, compatibility with low-frequency predictions for Coleman–Noll simple fluids, and agreement with data obtained via MAOS amplitude sweeps. Finally, the data are compared to predictions from the corotational Maxwell model to demonstrate the power of MAPS rheology in determining whether a constitutive model is consistent with a material’s time-dependent response.

Lennon, Kyle R. (ORCID:0000000212515461)↗
DOE OSTI2025

Wave ripples formed in ancient, ice-free lakes in Gale crater, Mars

Symmetrical wave ripples identified with NASA’s Curiosity rover in ancient lake deposits at Gale crater provide a key paleoclimate constraint for early Mars: At the time of ripple formation, climate conditions must have supported ice-free liquid water on the surface of Mars. These features are the most definitive examples of wave ripples on another planet. The ripples occur in two stratigraphic intervals within the orbitally defined Layered Sulfate Unit: a thin but laterally extensive unit at the base of the Amapari member of the Mirador formation, and a sandstone lens within the Contigo member of the Mirador formation. In both locations, the ripples have an average wavelength of ~4.5 centimeters. Internal laminae and ripple morphology show an architecture common in wave-influenced environments where wind-generated surface gravity waves mobilize bottom sediment in oscillatory flows. Their presence suggests formation in a shallow-water (<2 meters) setting that was open to the atmosphere, which requires atmospheric conditions that allow stable surface water.

Science & Technology - Other Topics↗
DOE OSTI2025

Observation of edge supercurrent in topological antiferromagnet MnBi 2 Te 4 -based Josephson junctions

Hybridizing superconductivity with topology and magnetism attracts growing interest in condensed matter physics. Here, we present our findings on the measurement of supercurrent induced in an intrinsic antiferromagnetic topological insulator MnBi 2 Te 4 . By constructing a MnBi 2 Te 4 proximity Josephson junction, we observed an anomalously large period of the Fraunhofer patterns, indicating a strong Josephson coupling state. As the MnBi 2 Te 4 thickness is reduced, a distinct asymmetric edge supercurrent emerges, aligning consistently with the observed oscillatory junction magnetoresistance. Leveraging this large asymmetric edge supercurrent, we have realized a nonvolatile Josephson diode device with programmable polarity, achieved through training with an out-of-plane magnetic field. Theoretical calculations substantiate that these behaviors are attributed to the interference between the highly asymmetric topological edge channel–mediated supercurrent induced in MnBi 2 Te 4 . Our study establishes this system as a promising avenue for investigating topological superconductivity, chiral Majorana edge modes, and advanced functionality device applications.

Science & Technology - Other Topics↗
DOE OSTI2025

Engineered chirality of one-dimensional nanowires

The origin and function of chirality in DNA, proteins, and other building blocks of life represent a central question in biology. Observations of spin polarization and magnetization associated with electron transport through chiral molecules, known collectively as the chiral induced spin selectivity effect, suggest that chirality improves electron transfer. Using reconfigurable nanoscale control over conductivity at the LaAlO 3 /SrTiO 3 interface, we create chiral electron potentials that explicitly lack mirror symmetry. Quantum transport measurements on these chiral nanowires reveal enhanced electron pairing persisting to high magnetic fields (up to 18 tesla) and oscillatory transmission resonances as functions of both magnetic field and chemical potential. We interpret these resonances as arising from an engineered axial spin-orbit interaction within the chiral region. The ability to create one-dimensional electron waveguides with this specificity creates opportunities to test, via analog quantum simulation, theories about chirality and spin-polarized electron transport in one-dimensional geometries.

Science & Technology - Other Topics↗
DOE OSTI2026

Multifunctional electrochemical memory stabilized by phase coexistence

Our growing computing needs, especially in applications that heavily rely on artificial intelligence (AI), motivate a search for new components that could substantially augment the performance of general-purpose digital computers. Beyond ON/OFF switching, new components with linear multistate analog resistive tuning, nonlinear volatile switching, spiking, oscillatory, stochastic and other complex functionalities could enable highly efficient neuromorphic computing schemes for AI information processing. Compared to the extreme multifunctionality of biological neurons, realizing all the above characteristics in a single, scalable analog component remains a grand challenge. Here we investigate electrochemical gating combined with localized thermal activation to program and switch a single, vertically integrated and dimensionally scaled electrothermal chemical random access memory (ETCRAM) with a channel and reservoir composed of phase-separated vanadium oxide. Closely related to electrochemical RAM (ECRAM), ETCRAM uses an integrated gate-heater electrode to overcome kinetic barriers that help retain states at ambient temperatures. In addition to synapse-like stable and programmable analog resistance states arising from redox-tunable phase coexistence, a single component exhibits neuron-like nonlinear conductance switching with a tunable threshold and self-driven dynamics owing to the thermally driven metal-insulator phase transition in vanadium dioxide. More broadly, we demonstrate that electrochemically stabilized phase coexistence could unlock analog electronics with novel functionality, stability, reconfigurability, and scalability.

Oh, Sangheon [Sandia National Lab. (SNL-CA), Liver↗
DOE OSTI2020

A Holistic Algorithmic Approach to Improving Accuracy, Robustness, and Computational Efficiency for Atmospheric Dynamics

Atmospheric weather and climate models must perform simulations very quickly to be useful. Therefore, modelers have traditionally focused on reducing computations as much as possible. However, in our new era of increasingly compute-capable hardware, data movement is now the prohibiting expense. This study examines the computational benefits of a new algorithmic approach to modeling atmospheric dynamics on scales relevant to weather and climate simulation. Rather than minimizing computations, this new approach considers the larger problem more holistically, including spatial accuracy, temporal accuracy, robustness (i.e., oscillations), on-node efficiency, and internode data transfers together at once. Numerical experiments demonstrate how computations can be strategically increased to simultaneously address each of these constraints while reducing data movement to adapt to modern accelerated hardware. The new algorithm can achieve at times up to 80% peak floating point throughput in single precision on the Nvidia Tesla V100 GPU, where the traditional approach is shown to only achieve single-digit floating point efficiency. Further, the new algorithm is twice as fast as a standard Runge--Kutta time integrator, and high-order accuracy with Weighted Essentially Non-Oscillatory (WENO) limiting came at less than 30% additional runtime cost on a GPU, thus increasing the accuracy per degree of freedom.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Robust Iterative Method for Symmetric Quantum Signal Processing in All Parameter Regimes

Here, this paper addresses the problem of solving nonlinear systems in the context of symmetric quantum signal processing (QSP), a powerful technique for implementing matrix functions on quantum computers. Symmetric QSP focuses on representing target polynomials as products of matrices in SU(2) that possess symmetry properties. We present a novel Newton’s method tailored for efficiently solving the nonlinear system involved in determining the phase factors within the symmetric QSP framework. Our method demonstrates rapid and robust convergence in all parameter regimes, including the challenging scenario with ill-conditioned Jacobian matrices, using standard double precision arithmetic operations. For instance, solving symmetric QSP for a highly oscillatory target function α cos(1000x) (polynomial degree ≈ 1433) takes 6 iterations to converge to machine precision when α = 0.9, and the number of iterations only increases to 18 iterations when α = 1 – 10 -9 with a highly ill-conditioned Jacobian matrix. Leveraging the matrix product state structure of symmetric QSP, the computation of the Jacobian matrix incurs a computational cost comparable to a single function evaluation. Moreover, we introduce a reformulation of symmetric QSP using real-number arithmetics, further enhancing the method’s efficiency. Extensive numerical tests validate the effectiveness and robustness of our approach, which has been implemented in the QSPPACK software package.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Logarithmic enhancements in conformal perturbation theory and their real time interpretation

We study various corrections of correlation functions to leading order in conformal perturbation theory, both on the cylinder and on the plane. Many problems on the cylinder are mathematically equivalent to those in the plane if we give the perturbations a position dependent scaling profile. The integrals to be done are then similar to those in the study of correlation functions with one additional insertion at the center of the profile. We will be primarily interested in the divergence structure of these corrections when computed in dimensional regularization. In particular, we show that the logarithmic divergences (enhancements) that show up in the plane under these circumstances can be understood in terms of resonant behavior in time dependent perturbation theory, for a transition between states that is induced by an oscillatory perturbation on the cylinder.

Physics↗
DOE OSTI2022

Combined Effect of the Wegener–Bergeron–Findeisen Mechanism and Large Eddies on Microphysics of Mixed-Phase Stratiform Clouds

Abstract The process of glaciation in mixed-phase stratiform clouds was investigated by a novel Lagrangian–Eulerian model (LEM) in which thousands of adjoining Lagrangian parcels moved within a turbulent-like velocity field with statistical parameters typical of the Arctic boundary layer. We used detailed bin microphysics to describe the condensation/evaporation processes in each parcel, in which droplets, aerosols, and ice particles were described using size distributions of 500 mass bins. The model also calculated aerosol mass inside droplets and ice particles. Gravitational sedimentation of droplets and ice particles was also accounted for. Assuming that droplet freezing is the primary source of ice particles, the Arctic clouds observed in Indirect and Semi-Direct Aerosol Campaign (ISDAC) were successfully simulated. The model showed that at a low ice particle concentration typical of ISDAC, large vortices (eddies) led to a quasi-stationary regime, in which mixed-phase St existed for a long time. The large eddies controlled the water partitioning in the mixed-phase clouds. Droplets formed and grew in updrafts, typically reaching the cloud top, and evaporated in downdrafts. Ice particles grew in updrafts and downdrafts. The Wegener–Bergeron–Findeisen (WBF) mechanism was efficient in downdrafts and some parts of updrafts, depending on ice concentration and vertical velocity. At low ice concentrations, the effect of ice on the phase partitioning was negligible. In this regime, liquid droplets were found near the cloud top, whereas ice particles precipitated through the cloud base. When ice concentration exceeded about 10 L −1 , the WBF mechanism led to glaciation of almost the entire cloud, with the exception of narrow cloud regions associated with strong updrafts. At ice particle concentrations of a few tens per liter, the oscillatory regime took place due to the ice–liquid interaction. The microphysical structure of mixed-phase St forms as a combined effect of cloud dynamics (large eddies) and the WBF mechanism.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of Advection Schemes on Tracer Interrelationships in Large-Eddy Simulations of Deep Convection

This study investigates the preservation of tracer interrelationships during advection in large-eddy simulations of an idealized deep convective cloud, which is particularly relevant to chemistry, aerosol, and cloud microphysics models. Employing the Cloud Model 1, advection is represented using third-, fifth-, and seventh-order weighted essentially non-oscillatory schemes. As a simplified analogy for cloud hydrometeors and aerosols, several inert passive tracers following linear and nonlinear relationships are initialized after the cloud reaches ~6-km depth. Numerical mixing in the simulated turbulent convective clouds leads to significant deviations from the initial nonlinear relationships between tracers. In these simulations, a considerable fraction of the grid points where the tracers’ nonlinear relationships are altered from advection are classified as unrealistic (e.g., ~13% for the environmental tracers on average), including errors from range-preserving unmixing and overshooting. Errors in the sum of three tracers are also relatively large, ranging between ~1% and 16% for 5% of the grid points in and near the cloud. The magnitude of unrealistic mixing and errors in the sum of three tracers generally increase with the order of accuracy of the advection scheme. These results are consistent across model grid spacings ranging from 50 to 200 m, and across three different flow realizations for each combination of grid spacing and advection scheme tested. Tests employing a previously proposed scalar normalization procedure show substantially reduced errors in the sum of three tracers with a relatively small negative impact on other tracer relationships. This analysis, therefore, suggests efficacy of the normalization procedure when applied to turbulent three-dimensional cloud simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Integrating deep neural networks with full-waveform inversion: Reparameterization, regularization, and uncertainty quantification

Full-waveform inversion (FWI) is an accurate imaging approach for modeling the velocity structure by minimizing the misfit between recorded and predicted seismic waveforms. However, the strong nonlinearity of FWI resulting from fitting oscillatory waveforms can trap the optimization in local minima. We have adopted a neural-network-based full-waveform inversion (NNFWI) method that integrates deep neural networks with FWI by representing the velocity model with a generative neural network. Neural networks can naturally introduce spatial correlations as regularization to the generated velocity model, which suppresses noise in the gradients and mitigates local minima. Furthermore, the velocity model generated by neural networks is input to the same partial differential equation (PDE) solvers used in conventional FWI. The gradients of the neural networks and PDEs are calculated using automatic differentiation, which back propagates gradients through the acoustic PDEs and neural network layers to update the weights of the generative neural network. Experiments on 1D velocity models, the Marmousi model, and the 2004 BP model determine that NNFWI can mitigate local minima, especially for imaging high-contrast features such as salt bodies, and it significantly improves the inversion in the presence of noise. Adding dropout layers to the neural network model also allows analyzing the uncertainty of the inversion results through Monte Carlo dropout. NNFWI opens a new pathway to combine deep learning and FWI for exploiting the characteristics of deep neural networks and the high accuracy of PDE solvers. Because NNFWI does not require extra training data and optimization loops, it provides an attractive and straightforward alternative to conventional FWI.

58 GEOSCIENCES↗
DOE OSTI2020

Aging phenomena in the two-dimensional complex Ginzburg-Landau equation

The complex Ginzburg-Landau equation with additive noise is a stochastic partial differential equation that describes a remarkably wide range of physical systems which include coupled non-linear oscillators subject to external noise near a Hopf bifurcation instability and spontaneous structure formation in non-equilibrium systems, e.g., in cyclically competing populations or oscillatory chemical reactions. Here, we employ a finite-difference method to numerically solve the noisy complex Ginzburg-Landau equation on a two-dimensional domain with the goal to investigate its non-equilibrium dynamics when the system is quenched into the “defocusing spiral quadrant”. We observe slow coarsening dynamics as oppositely charged topological defects annihilate each other, and characterize the ensuing aging scaling behavior. We conclude that the physical aging features in this system are governed by non-universal aging scaling exponents. We also investigate systems with control parameters residing in the “focusing quadrant”, and identify slow aging kinetics in that regime as well. Finally, we provide heuristic criteria for the existence of slow coarsening dynamics and physical aging behavior in the complex Ginzburg-Landau equation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗