Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Harmonic oscillator”

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

Curvature of Radial Electric Field Aggravates Edge Magnetohydrodynamics Mode in Toroidally Confined Plasmas

We show that the radial electric field (E r ) plays a dual role in edge magnetohydrodynamics (MHD) activity. While E r shear (first spatial derivative of E r ) dephases radial velocity and displacement, and so is stabilizing, a new finding here is that E r curvature (second spatial derivative of E r ) tends to synchronize the radial velocity and displacement, and so destabilizes MHD. As a highlighted result, we analytically demonstrate that E r curvature can destabilize an otherwise stable kink mode, and so form a joint vortexkink mode. The synergetic effects of E r shear and E r curvature in edge MHD extend the familiar E × B shearing paradigm. This theory thus explains the experimental findings that a deeper E × B well may aggravate edge MHD, and so trigger the formation of the edge harmonic oscillation. Here, a simple criterion linking E r structure and the edge MHD activity is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Photoluminescence spectra of point defects in semiconductors: Validation of first-principles calculations

Optically and magnetically active point defects in semiconductors are interesting platforms for the development of solid state quantum technologies. Their optical properties are usually probed by measuring photoluminescence spectra, which provide information on excitation energies and on the interaction of electrons with lattice vibrations. We present a combined computational and experimental study of photoluminescence spectra of defects in diamond and SiC, aimed at assessing the validity of theoretical and numerical approximations used in first-principles calculations, including the use of the Franck-Condon principle and the displaced harmonic oscillator approximation. We focus on prototypical examples of solid state qubits, the divacancy centers in SiC and the nitrogen-vacancy in diamond, and we report computed photoluminescence spectra as a function of temperature that are in very good agreement with the measured ones. As expected we find that the use of hybrid functionals leads to more accurate results than semilocal functionals. Interestingly our calculations show that constrained density functional theory (CDFT) and time-dependent hybrid DFT perform equally well in describing the excited state potential energy surface of triplet states; our findings indicate that CDFT, a relatively cheap computational approach, is sufficiently accurate for the calculations of photoluminescence spectra of the defects studied here. Finally, we find that only by correcting for finite-size effects and extrapolating to the dilute limit can one obtain a good agreement between theory and experiment. Our results provide a detailed validation protocol of first-principles calculations of photoluminescence spectra, necessary both for the interpretation of experiments and for robust predictions of the electronic properties of point defects in semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Lifetimes of excited states in $^{16}$C as a benchmark for ab initio developments

Lifetimes of higher-lying states ($2_2^+$ and $4_1^+$) in 16 C have been measured, employing the Gammasphere and Microball detector arrays, as key observables to test and refine ab initio calculations based on interactions developed within chiral Effective Field Theory. The presented experimental constraints to these lifetimes of $\tau ({2_2^+}) = [244, 446]\,~\textrm{fs}$ and $\tau ({4_1^+}) = [1.8, 4]\,~\textrm{ps}$, combined with previous results on the lifetime of the $2_1^+$ state of 16 C, provide a rather complete set of key observables to benchmark the theoretical developments. We present No-Core Shell-Model calculations using state-of-the-art chiral 2- (NN) and 3-nucleon (3N) interactions at next-to-next-to-next-to-leading order for both the NN and the 3N contributions and a generalized natural-orbital basis (instead of the conventional harmonic-oscillator single-particle basis) which reproduce, for the first time, the experimental findings remarkably well. The level of agreement of the new calculations as compared to the CD-Bonn meson-exchange NN interaction is notable and presents a critical benchmark for theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

The Ω c -puzzle solved by means of quark model predictions

The observation of five $$\varOmega _c $$ Ω c = ssc states by LHCb [Aaij et al. Phys. Rev. Lett. 118 , 182001(2017)] and the confirmation of four of them by Belle [Yelton et al. Phys. Rev. D 97 , 051102(2018)], may represent an important milestone in our understanding of the quark organization inside hadrons. By providing results for the spectrum of $$\varOmega _{ c}$$ Ω c baryons and predictions for their $$\varXi _{ c}^{+}K^{-}$$ Ξ c + K - and $$\varXi _{ c}'^{+}K^{-}$$ Ξ c ' + K - decay amplitudes within an harmonic oscillator based model, we suggest a possible solution to the $$\varOmega _{c}$$ Ω c quantum number puzzle and we extend our mass and decay width predictions to the $$\varOmega _b$$ Ω b states. Finally, we discuss why the set of $$\varOmega _{ c(b)}$$ Ω c ( b ) baryons is the most suitable environment to test the validity of three-quark and quark–diquark effective degrees of freedom.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Theory of the tertiary instability and the Dimits shift from reduced drift-wave models

Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity $U(x)$ with respect to the radial coordinate $x$. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $\gamma_{\rm TI}$ is derived within the modified Hasegawa--Wakatani model. We show that $\gamma_{\rm TI}$ equals the primary-instability growth rate plus a term that depends on the local $U''$; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry--Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear $|U'|$ is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Predictive Integrated Modeling of Heated Tokamak Plasmas (Final Technical Report)

The research carried out by the Lehigh University fusion group during this reporting period: (1) Self-consistent predictive core-pedestal ITER scenario modeling; (2) Microtearing instabilities and electron thermal transport in standard and spherical tokamak discharges; (3) Neural network model of the Multi-Mode anomalous transport module; (4) Anomalous transport in KSTAR, JET, EAST, Alcator C-Mod and DIII-D discharges; (5) Theory for transport in magnetized plasmas; (6) Control-oriented modeling of the q-profile dynamic response in preparation for advanced scenario control in KSTAR and in EAST; (7) Modeling of edge harmonic oscillation in DIII-D QH mode discharges; (8) Elucidating plasma dynamics in Hasegawa-Wakatani turbulence by information geometry; (9) Non-linear Langevin and Fractional Fokker–Planck equations for anomalous diffusion by Levy stable processes; and (10) Investigation of the plasma shaping effects on the H-mode pedestal structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

AToM: Advanced Tokamak Modeling Environment

Stability constraints play an important role in integrated modeling. The global plasma stability and macro-instabilities such as internal kink modes, neoclassical tearing modes, and edge localized perturbations limit the plasma performance and can result in large-scale transient events and plasma disruptions. These macro-instabilities can be also beneficial to the plasma performance. For example, the peeling models often lead to the edge localized modes (ELMs). However, the formation of stationary edge harmonic oscillations (EHOs) due to nonlinear interaction of peeling modes can result in a transition to the Quiescent H-mode (QH-mode) without ELMs. The profiles predicted with transport models need to be refined using the MHD stability constraints. This is important especially for transient and nonlinear stages of discharges such as ramp-up access to hybrid and steady state operation, and L- to H-mode transition as well as transition to the QH-mode. Several stability and MHD codes are already included in the OMFIT framework. These codes include BALOO, ELITE, GATO, M3DC1, MARS, and NIMROD. However, the verifications of stability constraints are currently mostly excluded from transport modeling workflows. The only exception is the EPED module which includes the ELITE predictions to limit the pedestal height. Here, we improved the OMFIT workflow to include stability calculations. Including the stability conditions to the integrated modeling workflow improved the robustness of the predictive modeling discharges. Being implemented in the workflow, the stability conditions can be also used for the stability analysis of experimental data. This improved the physics understanding of various discharge scenarios and can be used in experiment planning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Momentum Compaction, Phase-Slip Factor and Gamma Transition in a Synchrotron

A synchrotron is the workhorse in charged particle acceleration and is applied for charged particles acceleration and ion acceleration to the highest energies [1, p. 77]. During the acceleration of a charged particle, each particle experiences a longitudinal force opposite in sign to it’s “displacement from the central particle” like in a harmonic oscillator. This motion is similar also to the transverse motion of the particle in the accelerator (betatron oscillations) and this is the reason that the particle stays inside the bucket during the acceleration. In this technical note, we discuss the concepts like, transition energy, momentum compaction and phase slip factor associated with a synchrotron accelerator, and we use a simple example to present the physical meaning of the “transition energy”. We present a simple model of a circular accelerator having straight sections and bends only to provide an explanation and the physical meaning of transition energy.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

OC6 Phase III: Validation of Wind Turbine Aerodynamic Loading During Surge/Pitch Motion

The objective of Phase III of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project was to validate the accuracy of aerodynamic load predictions by offshore wind modeling tools for a floating offshore wind turbine as it experiences large surge-translational and pitch-rotational motion, as would occur during normal operation. The test data considered were generated at Politecnico di Milano – a wind tunnel with a robotic excitation system to emulate wave loading on the wind turbine. Testing was performed using a scaled version of the DTU 10-MW reference model, and motion was prescribed as harmonic oscillations in the surge and pitch directions, independently. A variety of models were examined in the project, performing a three-way validation between engineering-level tools, higher-fidelity tools, and measurement data from two wind tunnel experimental campaigns. The Load Cases (LC) considered in this testing are as follows: LC 1.X - Steady Wind LC 2.X - Unsteady Wind - Surge Motion LC 3.X - Unsteady Wind - Pitch Motion Details on the results of the OC6 Phase III campaign can be found in the following reference: Bergua R., et al. OC6 project Phase III: validation of the aerodynamic loading on a wind turbine rotor undergoing large motion caused by a floating support structure, Wind Energ. Sci., 8, 465–485, https://doi.org/10.5194/wes-8-465-2023, 2023.

17 WIND ENERGY↗
DOE OSTI2021

Bubble Behavior in a Vertically Vibrating Enclosed, Liquid-Filled Cylinder

When subjected to certain harmonic oscillations, the gas bubble in a partially liquid-filled, closed, vertical cylinder will break up. Under certain conditions, some of the gas will migrate to the bottom due to Bjerknes forces. At sufficiently large amplitudes, the bubble will break up into gas bubbles at the top and bottom ends of the cylinder. High-speed imaging captured the dynamics of bubble breakup and gas migration. Several parameters were investigated: oscillation frequency, oscillation acceleration, gas volume fraction, and liquid viscosity.

42 ENGINEERING↗
DOE OSTI2020

Nonadiabatic Energy Fluctuations of Scale-Invariant Quantum Systems in a Time-Dependent Trap

We consider the nonadiabatic energy fluctuations of a many-body system in a time-dependent harmonic trap. In the presence of scale-invariance, the dynamics becomes self-similar and the nondiabatic energy fluctuations can be found in terms of the initial expectation values of the second moments of the Hamiltonian, square position, and squeezing operators. Nonadiabatic features are expressed in terms of the scaling factor governing the size of the atomic cloud, which can be extracted from time-of-flight images. We apply this exact relation to a number of examples: the single-particle harmonic oscillator, the one-dimensional Calogero-Sutherland model, describing bosons with inverse-square interactions that includes the non-interacting Bose gas and the Tonks-Girdardeau gas as limiting cases, and the unitary Fermi gas. We illustrate these results for various expansion protocols involving sudden quenches of the trap frequency, linear ramps and shortcuts to adiabaticity. Our results pave the way to the experimental study of nonadiabatic energy fluctuations in driven quantum fluids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Quadrature Based Neural Network Learning of Stochastic Hamiltonian Systems

Hamiltonian Neural Networks (HNNs) provide structure-preserving learning of Hamiltonian systems. In this paper, we extend HNNs to structure-preserving inversion of stochastic Hamiltonian systems (SHSs) from observational data. We propose the quadrature-based models according to the integral form of the SHSs’ solutions, where we denoise the loss-by-moment calculations of the solutions. The integral pattern of the models transforms the source of the essential learning error from the discrepancy between the modified Hamiltonian and the true Hamiltonian in the classical HNN models into that between the integrals and their quadrature approximations. This transforms the challenging task of deriving the relation between the modified and the true Hamiltonians from the (stochastic) Hamilton–Jacobi PDEs, into the one that only requires invoking results from the numerical quadrature theory. Meanwhile, denoising via moments calculations gives a simpler data fitting method than, e.g., via probability density fitting, which may imply better generalization ability in certain circumstances. Numerical experiments validate the proposed learning strategy on several concrete Hamiltonian systems. The experimental results show that both the learned Hamiltonian function and the predicted solution of our quadrature-based model are more accurate than that of the corrected symplectic HNN method on a harmonic oscillator, and the three-point Gaussian quadrature-based model produces higher accuracy in long-time prediction than the Kramers–Moyal method and the numerics-informed likelihood method on the stochastic Kubo oscillator as well as other two stochastic systems with non-polynomial Hamiltonian functions. Moreover, the Hamiltonian learning error εH arising from the Gaussian quadrature-based model is lower than that from Simpson’s quadrature-based model. These demonstrate the superiority of our approach in learning accuracy and long-time prediction ability compared to certain existing methods and exhibit its potential to improve learning accuracy via applying precise quadrature formulae.

Mathematics↗
DOE OSTI2020

Canonical Transformation of Potential Model Hamiltonian Mechanics to Geometrical Form I

Using the methods of symplectic geometry, we establish the existence of a canonical transformation from potential model Hamiltonians of standard form in a Euclidean space to an equivalent geometrical form on a manifold, where the corresponding motions are along geodesic curves. The advantage of this representation is that it admits the computation of geodesic deviation as a test for local stability, shown in recent previous studies to be a very effective criterion for the stability of the orbits generated by the potential model Hamiltonian. We describe here an algorithm for finding the generating function for the canonical transformation and describe some of the properties of this mapping under local diffeomorphisms. We give a convergence proof for this algorithm for the one-dimensional case, and provide a precise geometric formulation of geodesic deviation which relates the stability of the motion in the geometric form to that of the Hamiltonian standard form. We apply our methods to a simple one-dimensional harmonic oscillator and conclude with a discussion of the relation of bounded domains in the two representations for which Morse theory would be applicable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Examining AGN UV/Optical Variability beyond the Simple Damped Random Walk

We present damped harmonic oscillator (DHO) light-curve modeling for a sample of 12,714 spectroscopically confirmed quasars in the Sloan Digital Sky Survey Stripe 82 region. DHO is a second-order continuous-time autoregressive moving-average process, which can be fully described using four independent parameters: a natural oscillation frequency (ω 0 ), a damping ratio (ξ), a characteristic perturbation timescale (τ perturb ), and an amplitude for the perturbing white noise (σ ϵ ). The asymptotic variability amplitude of a DHO process is quantified by σ DHO —a function of ω 0 , ξ, τ perturb , and σ ϵ . We find that both τ perturb and σ ϵ follow different dependencies with rest-frame wavelength (λ RF ) on either side of 2500 Å, whereas σ DHO follows a single power-law relation with λ RF . After correcting for wavelength dependence, σ DHO exhibits anticorrelations with both the Eddington ratio and the black hole mass, while τ perturb —with a typical value of days in the rest frame—shows an anticorrelation with the bolometric luminosity. Modeling active galactic nuclei (AGN) variability as a DHO offers more insight into the workings of accretion disks close to the supermassive black holes at the center of AGN. The newly discovered short-term variability (characterized by τ perturb and σ ϵ ) and its correlation with bolometric luminosity pave the way for new algorithms that will derive fundamental properties (e.g., Eddington ratio) of AGN using photometric data alone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Augmented Gaussian random field: Theory and computation

We propose the novel augmented Gaussian random field (AGRF), which is a universal framework incorporating the data of observable and derivatives of any order. Rigorous theory is established. We prove that under certain conditions, the observable and its derivatives of any order are governed by a single Gaussian random field, which is the aforementioned AGRF. As a corollary, the statement "the derivative of a Gaussian process remains a Gaussian process" is validated, since the derivative is represented by a part of the AGRF. Moreover, a computational method corresponding to the universal AGRF framework is constructed. Both noiseless and noisy scenarios are considered. Formulas of the posterior distributions are deduced in a nice closed form. A significant advantage of our computational method is that the universal AGRF framework provides a natural way to incorporate arbitrary order derivatives and deal with missing data. We use four numerical examples to demonstrate the effectiveness of the computational method. The numerical examples are composite function, damped harmonic oscillator, Korteweg-De Vries equation, and Burgers' equation.

Mathematics↗