Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Impacts of swinging door lossy compression of synchrophasor data

As the number of operational Phasor Measurement Units (PMUs) is increasing, compression is becoming crucial for efficiently transferring and storing enormous amount of synchrophasor data continuously measured from a power system. Some utilities in the U.S. have already adopted a nonlinear compression technique, called Swinging Door Algorithm (SDA), assuming compressed data meets technical and regulatory requirements without resulting in any degradation of the measurements. This paper studies impacts of this compression algorithm on field measured and simulated PMU data, and presents results from spectral analyses of raw and decompressed measurements. The results show blind application of the SDA compression alters frequency content in the measurements and introduces false harmonics and aliasing. These can lead to false conclusions being drawn from the measurements about the power system condition.

42 ENGINEERING↗
DOE OSTI2026

Simulating plasma wave propagation on a superconducting quantum chip

Quantum computers may one day enable the efficient simulation of strongly coupled plasmas that lie beyond the reach of classical computation in regimes where quantum effects are important and the scale separation is large. Here, in this article, we take a first step toward efficient simulation of quantum plasmas by demonstrating linear plasma wave propagation on a superconducting quantum chip. Using high-fidelity and highly expressive device-native gates, combined with an error-mitigation technique, we simulate the scattering of laser pulses from inhomogeneous plasmas. Our approach is made feasible by the identification of a suitable local spin model whose excitations mimic plasma waves, and whose circuit implementation requires a lower gate count than other proposed approaches that would require a future fault-tolerant quantum computer. This work opens avenues to study more complicated phenomena that cannot be simulated efficiently on classical computers, such as nonlinear quantum dynamics when strongly coupled plasmas are driven out of equilibrium.

general physics↗
DOE OSTI2021

Inference of neutrino flavor evolution through data assimilation and neural differential equations

The evolution of neutrino flavor in dense environments such as core-collapse supernovae and binary compact object mergers constitutes an important and unsolved problem. Its solution has potential implications for the dynamics and heavy-element nucleosynthesis in these environments. In this paper, we build upon recent work to explore inference-based techniques for the estimation of model parameters and neutrino flavor evolution histories. Here, we combine data assimilation, ordinary differential equation solvers, and neural networks to craft an inference approach tailored for nonlinear dynamical systems. Using this architecture, and a simple two-neutrino-beam, two-flavor model, we compare the performances of nine different optimization algorithms and expand upon previous assessments of the efficacy of inference for tackling problems in flavor evolution. We find that employing this new architecture, together with evolutionary optimization algorithms, accurately captures flavor histories in the small-scale model and allows us to quickly explore both model parameters and initial flavor content. In future work we plan to extend these inference techniques to large numbers of neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Long-Range Lamellar Alignment in Diblock Bottlebrush Copolymers via Controlled Oscillatory Shear

A simple strategy is presented to achieve well-ordered nanostructures in microphase-separated diblock bottlebrush copolymers (dbBB) for potential opportunities in nanotechnology that require large grain sizes and directed orientation. The unique architecture and relaxation processes of dbBBs offer a pathway to enhanced dynamic ordering and reorientation over macroscopic scales. Here, controlled shear is used to align lamellar domains. The high molecular mobility of dbBBs in the melt state, as shown by linear and nonlinear viscoelastic characterization, supports this alignment process. When subjected to large-amplitude oscillatory shear at controlled frequency (ω), strain amplitude (γ), and temperature (T shear ), microphase-separated lamellae transition into highly aligned states as determined by small-angle X-ray scattering (SAXS) and transmission electron microscopy (TEM). High throughput synchrotron SAXS is used to rapidly characterize microphase separation and lamellar alignment across the entire bulk sample. Lamellae align parallel to shear planes with an orientation order parameter (S) approaching S = 0.8 over a large cubic volume (V ~ mm 3 ).

36 MATERIALS SCIENCE↗
DOE OSTI2025

DESI DR2 Results I: Baryon Acoustic Oscillations from the Lyman Alpha Forest

We present the Baryon Acoustic Oscillation (BAO) measurements with the Lyman-alpha (LyA) forest from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our BAO measurements include both the auto-correlation of the LyA forest absorption observed in the spectra of high-redshift quasars and the cross-correlation of the absorption with the quasar positions. The total sample size is approximately a factor of two larger than the DR1 dataset, with forest measurements in over 820,000 quasar spectra and the positions of over 1.2 million quasars. We describe several significant improvements to our analysis in this paper, and two supporting papers describe improvements to the synthetic datasets that we use for validation and how we identify damped LyA absorbers. Our main result is that we have measured the BAO scale with a statistical precision of 1.1% along and 1.3% transverse to the line of sight, for a combined precision of 0.65% on the isotropic BAO scale at $z_{eff} = 2.33$. This excellent precision, combined with recent theoretical studies of the BAO shift due to nonlinear growth, motivated us to include a systematic error term in LyA BAO analysis for the first time. We measure the ratios $D_H(z_{eff})/r_d = 8.632 \pm 0.098 \pm 0.026$ and $D_M(z_{eff})/r_d = 38.99 \pm 0.52 \pm 0.12$, where $D_H = c/H(z)$ is the Hubble distance, $D_M$ is the transverse comoving distance, $r_d$ is the sound horizon at the drag epoch, and we quote both the statistical and the theoretical systematic uncertainty. The companion paper presents the BAO measurements at lower redshifts from the same dataset and the cosmological interpretation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Relationship between the contact force strength and numerical inaccuracies in piecewise-smooth systems

This work studies the different types of behavior and inaccuracies that can occur when contact is not adequately accounted for in a dynamical system with freeplay, as the strength of the contact stiffness increases. The MATLAB® ode45 time integration solver, with the built-in Event Location capability, is first validated using past experimental data from a forced Duffing oscillator with freeplay. Next, numerical results utilizing event location are compared to results neglecting event location in order to highlight possible numerical errors and effects on multistable dynamical responses. Inaccuracies tend to occur in two different ways. First, neglecting event location can affect the boundaries between basins of attraction. Second, neglecting event location has little effect on the behaviors of the attractor solutions themselves besides merely resembling poorly converged solutions. Errors are less pronounced at the limits of soft or hard contact stiffness. Finally, this study shows the importance of accurately solving piecewise-smooth systems and the existing correlation between the strength of the contact force and possible numerical inaccuracies.

42 ENGINEERING↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Learning Distributed Geometric Koopman Operator for Sparse Networked Dynamical Systems

Koopman operator theory provides an alternative to study nonlinear networked dynamical systems by mapping the state space to an abstract higher dimensional space where the system evolution is linear. Recent works show the application of graph neural networks (GNNs) to learn state to object-centric embeddings and achieve centralized block-wise computation of Koopman operator (KO) under additional assumptions on the underlying agents properties and constraints on the KO structure. However, the computational complexity of learning the Koopman increases exponentially for networked systems where the number of possible system states grows in a combinatorial fashion with the number of nodes. The learning challenge is further amplified for sparse networks by two factors: 1) sample sparsity for learning the Koopman operator in the non-linear space, and 2) the divergence in the dynamics of individual nodes or from one subgraph to another. Our work aims to address these challenge by formulating the representation learning of networked dynamical systems into a multi-agent paradigm and learning the Koopman operator in a distributive manner. The computational as well as performance advantages of distributed Koopman is predominant for sparse networks whereas for fully connected networks, it is shown to coincide with the centralized one. The empirical study on rope system, network of oscillators and a synthetic power system show comparable and superior performance along with computational benefits with the state-of-the-art methods.

Mukherjee, Sayak↗
DOE OSTI2021

Driven-dissipative phase transition in a Kerr oscillator: From semiclassical PT symmetry to quantum fluctuations

We study a minimal model that has a driven-dissipative quantum phase transition, namely a Kerr non-linear oscillator subject to driving and dissipation. Using mean-field theory, exact diagonalization, and the Keldysh formalism, we analyze the critical phenomena in this system, showing which aspects can be captured by each approach and how the approaches complement each other. Then critical scaling and finite-size scaling are calculated analytically using the quantum Langevin equation. Furthermore, the physics contained in this simple model is surprisingly rich: it includes a continuous phase transition, Z 2 symmetry breaking, PT symmetry, state squeezing, and critical fluctuations. Due to its simplicity and solvability, this model can serve as a paradigm for exploration of open quantum many-body physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Quantum scars of bosons with correlated hopping

Abstract Recent experiments on Rydberg atom arrays have found evidence of anomalously slow thermalization and persistent density oscillations, which have been interpreted as a many-body analog of the phenomenon of quantum scars. Periodic dynamics and atypical scarred eigenstates originate from a “hard” kinetic constraint: the neighboring Rydberg atoms cannot be simultaneously excited. Here we propose a realization of quantum many-body scars in a 1D bosonic lattice model with a “soft” constraint in the form of density-assisted hopping. We discuss the relation of this model to the standard Bose-Hubbard model and possible experimental realizations using ultracold atoms. We find that this model exhibits similar phenomenology to the Rydberg atom chain, including weakly entangled eigenstates at high energy densities and the presence of a large number of exact zero energy states, with distinct algebraic structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Entangled two-photon absorption with Brownian-oscillator fluctuations

Here, we theoretically investigate the two-photon absorption signals of a three-band (g, e, f) system diagonally coupled to an over-damped Brownian oscillator bath, which induces random Gaussian modulations of energy levels with an arbitrary degree of correlation. For fast modulation, extra 2ω eg and 2ω fe peaks may obscure the g–f transitions in the classical two-photon absorption (CTPA) spectra for nearly resonant e states. These peaks arise from one-photon resonant g–e or e–f transitions. In the slow modulation limit, these peaks vanish because of the short tails of the Gaussian line shape. CTPA strongly depends on the correlations between energy fluctuations. In entangled two-photon absorption, the extra peaks are eliminated because of the broad one-photon but narrow two-photon spectrum of the twin photons. The variation of the coherences between f states with the correlation between energy fluctuations is explored.

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

Producing entangled photon pairs and quantum squeezed states in plasmas

Plasma is capable of mediating the conversion of two pump photons into two different photons through a relativistic four-wave mixing nonlinearity. Spontaneously created photon pairs are emitted at symmetric angles with respect to the colinear pump direction, and the emission rate is largest if they have identical frequency. Thus, two orthogonally polarized pumps can produce polarization-entangled photon pairs through a millimeter-long homogeneous plasma. Here, the noise from Raman scattering can be avoided if the pump detuning differs from twice the plasma frequency. However, pump detuning exactly equal to twice the plasma frequency can significantly enhance the interaction rate, which allows for the production of strong two-mode squeezed states. Remarkably, the amplified noise from Raman scattering are correlated and hence can be suppressed in one of the output quadratures, thereby maintaining the squeezing magnitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analytical WKB theory for high-harmonic generation and its application to massive Dirac electrons

Here we propose an analytical approach to high-harmonic generation (HHG) for nonperturbative low-frequency and high-intensity fields based on the (Jeffreys-)Wentzel-Kramers-Brillouin (WKB) approximation. By properly taking into account Stokes phenomena of WKB solutions, we obtain wave functions that systematically include the repetitive dynamics of production and acceleration of electron-hole pairs and quantum interference due to phase accumulation between different pair production times (Stückelberg phase). Using the obtained wave functions without relying on any phenomenological assumptions, we explicitly compute electric current (including intra- and interband contributions) as the source of HHG for a massive Dirac system in (1+1) dimensions under an ac electric field. We demonstrate that the WKB approximation agrees well with numerical results obtained by solving the time-dependent Schrödinger equation and point out that the quantum interference is important in HHG. We also predict in the deep nonperturbative regime that (1) harmonic intensities oscillate with respect to electric-field amplitude E 0 and frequency Ω with a period determined by the Stückelberg phase, (2) the cutoff order of HHG is determined by 2eE 0 ℏΩ 2 , with e being the electron charge, and that (3) noninteger harmonics, controlled by the Stückelberg phase, appear as a transient effect. Our WKB theory is particularly suited for a parameter regime, where the Keldysh parameter γ = (Δ/2)Ω/eE 0 ,with Δ being the gap size, is small. This parameter regime corresponds to intense lasers in the terahertz regime for realistic massive Dirac materials. Our analysis implies that the so-called HHG plateau can be observed at the terahertz frequency within the current technology.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamics and bifurcations of laminar annular swirling and non-swirling jets

This paper presents bifurcation analyses characterising the nonlinear dynamics of fully developed laminar annular jets with respect to the centrebody diameter ${ {d}}$ , Reynolds number ${ {Re}}$ , and swirl ratio ${ {S}}$ . Similar flows appear in numerous applications and feature a vibrant range of topological and dynamical characteristics associated with phenomena including shear layer separation and vortex breakdown. Our results begin by describing the non-monotonic evolution of the axisymmetric jet's steady topology under varying ${ {S}}$ . In accord with earlier reports, the jet progresses through a sequence of wake, breakdown and wall jet regimes in a qualitatively similar manner across a wide span of ${ {d}}$ and ${ {Re}}$ values. In the wake regime, the non-swirling jet bifurcates to a plane-symmetric, but not axisymmetric, steady flow pattern beyond a ${ {d}}$ -dependent critical ${ {Re}}$ value. With further increase in ${ {Re}}$ , the steady non-swirling jet destabilises subsequently via multiple distinct Hopf bifurcations. Introducing ${ {S}}>0$ to the jet also induces unsteadiness by twisting the singly azimuthally periodic ( $|m|=1$ ) asymmetric wake structure and causing it to precess periodically in time about the central axis. Intermediate swirl stabilises this unsteady dynamics and restores the jet's axisymmetry. This stabilising effect is then reversed in the breakdown regime at higher ${ {S}}$ , where a variety of different $|m|=1$ and $|m|=2$ instabilities bifurcate from the steady flow as ${ {S}}$ is increased. Several instances of hysteresis and subcritical behaviour are reported and discussed, including one that manifests precessing vortex core oscillations.

Douglas, Christopher M. (ORCID:0000000259683315)↗
DOE OSTI2023

Towards understanding the robust strengthening of ENSO and more frequent extreme El Niño events in CMIP6 global warming simulations

The El Niño Southern Oscillation (ENSO) has profound impacts on weather patterns across the globe, yet there is no consensus on its response to global warming. Several modelling studies find a stronger ENSO in global warming scenarios, while other studies suggest ENSO weakening. Using a broad range of models from the Coupled Model Intercomparison Project phase 6 (CMIP6) and four types of warming experiments, here we show that the majority of the models predict a stronger ENSO by century-end in Shared Social Pathway (SSP) experiments, and in idealized 1pctCO 2 and abrupt 4xCO 2 experiments. Several models, however, do predict no change or ENSO weakening, especially in the idealized experiments. Critically, the strongest forcing (abrupt-4xCO 2 ) does not induce the strongest ENSO response, while differences between the models are much greater than those between warming scenarios. For the long-term response (over 1000 years) the models disagree even on the sign of change. Furthermore, changes in ENSO sea surface temperature (SST) variability are only modestly correlated with the tropical Pacific mean state change. The highest correlation for ENSO SST amplitude is found with the mean zonal SST gradient in the SSP5-8.5 experiment (R = –0.58). In contrast, changes in ENSO rainfall variability correlate well with changes in the mean state, as well as with changes in ENSO SST variability. When evaluating the Bjerknes Stability Index for a subset of models, we find that it is not a reliable predictor of ENSO strengthening, as this index tends to predict greater stability with warming. We argue that the enhanced ENSO stability is offset by increases in atmospheric noise or/and potential nonlinear effects. However, a robust inter-model mechanism that can explain a stronger ENSO simulated with global warming is still lacking. Therefore, caution should be exercised when considering ENSO changes based on a single model or warming scenario.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Joint Experimental and Computational Characterization of Sum-Frequency Generation between a Continuous Wave Laser and an Ultrafast Frequency Comb Laser for Tunable Laser Development

Ultrafast optical frequency combs allow for both high spectral and temporal resolution in molecular spectroscopy and have become a powerful tool in many areas of chemistry and physics. Ultrafast lasers and frequency combs generated from ultrafast mode-locked lasers often need to be converted to other wavelengths. Commonly used wavelength conversions are optical parametric oscillators, which require an external optical cavity, and supercontinuum generation combined with optical parametric amplifiers. Whether commercial or home-built, these systems are complex and costly. Here, we investigate an alternative, simple, and easy-to-implement approach to tunable frequency comb ultrafast lasers enabled by new continuous-wave laser technology. Sum-frequency generation between an Nd:YAG continuous-wave laser and a Yb:fiber femtosecond frequency comb in a beta-barium borate (BBO) crystal is explored. The resulting sum-frequency beam is a pulsed frequency comb with the same repetition rate as the Yb:fiber source. SNLO simulation software is used to simulate the results and provide benchmarks for designing future systems to achieve wavelength conversion and tunability in otherwise difficult-to-reach spectral regions.

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

Limits on Mode Coherence in Pulsating DA White Dwarfs Due to a Nonstatic Convection Zone

The standard theory of pulsations deals with the frequencies and growth rates of infinitesimal perturbations in a stellar model. Modes which are calculated to be linearly driven should increase their amplitudes exponentially with time; the fact that nearly constant amplitudes are usually observed is evidence that nonlinear mechanisms inhibit the growth of finite amplitude pulsations. Models predict that the mass of convection zones in pulsating hydrogen-atmosphere (DAV) white dwarfs is very sensitive to temperature (i.e., M CZ α T –90 eff ), leading to the possibility that even low-amplitude pulsators may experience significant nonlinear effects. In particular, the outer turning point of finite-amplitude g-mode pulsations can vary with the local surface temperature, producing a reflected wave that is out of phase with what is required for a standing wave. This can lead to a lack of coherence of the mode and a reduction in its global amplitude. In this paper we show that: (1) whether a mode is calculated to propagate to the base of the convection zone is an accurate predictor of its width in the Fourier spectrum, (2) the phase shifts produced by reflection from the outer turning point are large enough to produce significant damping, and (3) amplitudes and periods are predicted to increase from the blue edge to the middle of the instability strip, and subsequently decrease as the red edge is approached. Furthermore, this amplitude decrease is in agreement with the observational data while the period decrease has not yet been systematically studied.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Three-wave coupling observed between a shear Alfvén wave and a kink-unstable magnetic flux rope

Results from a laboratory experiment are presented in which, for the first time, a shear Alfvén wave is launched using an antenna in a current-carrying plasma column that is tailored to be either stable or unstable to the kink oscillation. As the plasma is driven kink unstable, the frequency power spectrum of the Alfvén wave evolves from a single peak to a peak with multiple sidebands separated by integer multiples of the kink frequency. The main sidebands (one on either side of the launched wave peak in the power spectrum) are analyzed using azimuthal wavenumber matching, perpendicular and parallel wavenumber decomposition, and bispectral time series analysis. The dispersion relation and three-wave matching conditions are satisfied, given each sideband is a propagating Alfvén wave that results from the interaction of the pump Alfvén wave and the co-propagating component of a half-wavelength, standing kink mode. The interaction is shown to generate smaller perpendicular wavelength Alfvén waves that drive energy transport to scales that will approach the dissipation scale of k⊥ρs=1, with k⊥ being the perpendicular wavenumber and ρs being the ion gyroradius at the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗