Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forced 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 OSTI2023

Evaluating route to impact convergence of the harmonic balance method for piecewise-smooth systems

Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.

42 ENGINEERING↗
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 OSTI2021

Insights on the continuous representations of piecewise-smooth nonlinear systems: limits of applicability and effectiveness

Dynamical systems subject to intermittent contact are often modeled with piecewise-smooth contact forces. However, the discontinuous nature of the contact can cause inaccuracies in numerical results or failure in numerical solvers. Representing the piecewise contact force with a continuous and smooth function can mitigate these problems, but not all continuous representations may be appropriate for this use. Here, five representations used by previous researchers (polynomial, rational polynomial, hyperbolic tangent, arctangent, and logarithm-arctangent functions) are studied to determine which ones most accurately capture nonlinear behaviors including super- and subharmonic resonances, multiple solutions, and chaos. The test case is a single-DOF forced Duffing oscillator with freeplay nonlinearity, solved using direct time integration. This work intends to expand on past studies by determining the limits of applicability for each representation and what numerical problems may occur.

42 ENGINEERING↗
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 OSTI2023

Climate Change Signal in Atlantic Tropical Cyclones Today and Near Future

This manuscript discusses the challenges in detecting and attributing recently observed trends in the Atlantic tropical cyclone (TC) and the epistemic uncertainty we face in assessing future risk. We use synthetic storms downscaled from five CMIP5 models by the Columbia HAZard model (CHAZ), and directly simulated storms from high-resolution climate models. We examine three aspects of recent TC activity: the upward trend and multi-decadal oscillation of the annual frequency, the increase in storm wind intensity, and the decrease in forward speed. Some data sets suggest that these trends and oscillation are forced while others suggest that they can be explained by natural variability. Projections under warming climate scenarios also show a wide range of possibilities, especially for the annual frequencies, which increase or decrease depending on the choice of moisture variable used in the CHAZ model and on the choice of climate model. The uncertainties in the annual frequency lead to epistemic uncertainties in TC risk assessment. Here, we investigate the potential for reduction of these epistemic uncertainties through a statistical practice, namely likelihood analysis. We find that historical observations are more consistent with the simulations with increasing frequency than those with decreasing frequency, but we are not able to rule out the latter. We argue that the most rational way to treat epistemic uncertainty is to consider all outcomes contained in the results. In the context of risk assessment, since the results contain possible outcomes in which TC risk is increasing, this view implies that the risk is increasing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Terahertz electric-field-driven dynamical multiferroicity in SrTiO 3

The emergence of collective order in matter is among the most fundamental and intriguing phenomena in physics. In recent years, the dynamical control and creation of novel ordered states of matter not accessible in thermodynamic equilibrium is receiving much attention. The theoretical concept of dynamical multiferroicity has been introduced to describe the emergence of magnetization due to time-dependent electric polarization in non-ferromagnetic materials. In simple terms, the coherent rotating motion of the ions in a crystal induces a magnetic moment along the axis of rotation. Here we provide experimental evidence of room-temperature magnetization in the archetypal paraelectric perovskite SrTiO 3 due to this mechanism. We resonantly drive the infrared-active soft phonon mode with an intense circularly polarized terahertz electric field and detect the time-resolved magneto-optical Kerr effect. A simple model, which includes two coupled nonlinear oscillators whose forces and couplings are derived with ab initio calculations using self-consistent phonon theory at a finite temperature, reproduces qualitatively our experimental observations. A quantitatively correct magnitude was obtained for the effect by also considering the phonon analogue of the reciprocal of the Einstein–de Haas effect, which is also called the Barnett effect, in which the total angular momentum from the phonon order is transferred to the electronic one. Our findings show a new path for the control of magnetism, for example, for ultrafast magnetic switches, by coherently controlling the lattice vibrations with light.

36 MATERIALS SCIENCE↗
DOE OSTI2020

OC6 Phase I: Investigating the underprediction of low-frequency hydrodynamic loads and responses of a floating wind turbine

Phase I of the OC6 project is focused on examining why offshore wind design tools underpredict the response (loads/motion) of the OC5-DeepCwind semisubmersible at its surge and pitch natural frequencies. Previous investigations showed that the underprediction was primarily related to nonlinear hydrodynamic loading, so two new validation campaigns were performed to separately examine the different hydrodynamic load components. In this paper, we validate a variety of tools against this new test data, focusing on the ability to accurately model the low-frequency loads on a semisubmersible floater when held fixed under wave excitation and when forced to oscillate in the surge direction. However, it is observed that models providing better load predictions in these two scenarios do not necessarily produce a more accurate motion response in a moored configuration.

17 WIND ENERGY↗
DOE OSTI2026

OC6 Phase Ia - Nonlinear hydrodynamic loading validation dataset

Two validation campaigns were examined within the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) Phase 1 project to examine the modeling tools' underprediction of loads and motion of a floating wind semisubmersible (semi) at their surge and pitch natural frequencies. These campaigns were performed at the Maritime Research Institute Netherlands (MARIN) in 2017 and 2018. The load cases (LC) considered include: LC1 – Load measurements across semi under current loading; LC2 - Load measurements across semi under forced surge oscillation; LC3 – Load measurements across semi under wave loading, while held fixed; LC4 – Free-decay motion measurements in surge, pitch, and heave; and LC5 – Motion measurements under wave loading. Details on the results from the OC6 Phase Ia project can be found in the reference, “OC6 Phase 1: Investigating the underprediction of low-frequency hydrodynamic loads and responses of floating wind turbines”, J Phys: Conf Series 1618 032033.

17 WIND ENERGY↗
DOE OSTI2023

OC6 project Phase III: validation of the aerodynamic loading on a wind turbine rotor undergoing large motion caused by a floating support structure

Abstract. This paper provides a summary of the work done within Phase III of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project, under the International Energy Agency Wind Technology Collaboration Programme Task 30. This phase focused on validating the aerodynamic loading on a wind turbine rotor undergoing large motion caused by a floating support structure. Numerical models of the Technical University of Denmark 10 MW reference wind turbine were validated using measurement data from a 1:75 scale test performed during the UNsteady Aerodynamics for FLOating Wind (UNAFLOW) project and a follow-on experimental campaign, both performed at the Politecnico di Milano wind tunnel. Validation of the models was performed by comparing the loads for steady (fixed platform) and unsteady (harmonic motion of the platform) wind conditions. For the unsteady wind conditions, the platform was forced to oscillate in the surge and pitch directions under several frequencies and amplitudes. These oscillations result in a wind variation that impacts the rotor loads (e.g., thrust and torque). For the conditions studied in these tests, the system aerodynamic response was almost steady. Only a small hysteresis in airfoil performance undergoing angle of attack variations in attached flow was observed. During the experiments, the rotor speed and blade pitch angle were held constant. However, in real wind turbine operating conditions, the surge and pitch variations would result in rotor speed variations and/or blade pitch actuations, depending on the wind turbine controller region that the system is operating. Additional simulations with these control parameters were conducted to verify the fidelity of different models. Participant results showed, in general, a good agreement with the experimental measurements and the need to account for dynamic inflow when there are changes in the flow conditions due to the rotor speed variations or blade pitch actuations in response to surge and pitch motion. Numerical models not accounting for dynamic inflow effects predicted rotor loads that were 9 % lower in amplitude during rotor speed variations and 18 % higher in amplitude during blade pitch actuations.

17 WIND ENERGY↗
DOE OSTI2023

Efficacy of vibro-impact energy harvesting absorbers on controlling dynamical systems under vortex-induced vibrations and base excitation

Here this study investigates improving the efficacy of an energy harvesting absorber's ability to control a structure under vortex-induced vibrations, base excitation, and a combination of the two by including mechanical amplitude stoppers. The nonlinear reduced-order model is developed through modifying trilinear spring models to represent the impact forces, a modified van der Pol oscillator to represent the forcing due to the vortex-induced vibrations and using the Euler-Lagrange principle to express the equations of motion. It is seen that a soft stopper stiffness and a 5mm gap performs the most effectively of increasing the power generated from the absorber while still greatly reducing the primary structure's amplitude. By changing the stopper's location towards the middle of the energy harvesting absorber, the large effects of the impact forces are reduced and improves the efficacy of medium and hard stopper stiffnesses to generate near the amount of power the soft stopper does, while greatly improving the control of the primary structure. When the system is under combined loadings, the large oscillations of the synchronization region cause the effective configuration to be that of a 27.5 mm gap with soft stiffnesses. The results shows that medium stiffness stoppers with small gaps generate large aperiodic regions due to the high impact force. When the oscillations are close to the stoppers, the beating phenomenon is observed and is not overpowered by the vibro-impact force.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Feedback and Oscillations: Constructing Feedback Systems for Root Cause Analysis of Oscillations in Power Grids

Dynamic phenomena linked to inverter-based resources (IBRs) have gained global attention. Several IBR-induced dynamics have caused bulk power system-connected wind or solar power plants to trip, and some have even led to widespread outages. In addition, many oscillations have been observed involving IBR power plants. In 2023, the IEEE Power & Energy Society (PES) IBR Subsynchronous Oscillations (SSO) task force published a journal article, “Real-World Subsynchronous Oscillation Events in Power Grids With High Penetrations of Inverter-Based Resources,” in which 19 IBR oscillation events were examined for their causation. Earlier in 2020, another PES task force article, “Definition and Classification of Power System Stability-Revisited & Extended,” authored by prominent academics, introduced converter-driven stability as a new category of stability. The international power grid industry community also took action by publishing the CIGRE Green Book, Power System Dynamic Modelling and Analysis in Evolving Networks (led by Babak Badrzadeh and Zia Emin) in 2024. In August 2024, the Energy Systems Integration Group (ESIG) released a practical guide led by Nick Miller, “Diagnosis and Mitigation of Observed Oscillations in IBR-Dominant Power System: A Practical Guide.” The goal of the guide is to assist practicing engineers in making initial judgments and conducting detailed analyses about oscillations. Finally, when addressing the classification of stability and oscillations, the guide emphasizes a causality-based taxonomy for grouping, such as voltage control-induced oscillations, synchronization-induced oscillations, and frequency or active power control-induced oscillations.

Fan, Lingling [Univ. of South Florida, Tampa, FL (↗
DOE OSTI2020

Linear thermomagnetic energy harvester for low-grade thermal energy harvesting

Low-grade thermal energy, either from waste heat or from natural resources, constitutes an enormous energy reserve that remains to be fully harvested. Harvesting low-grade heat is challenging because of the low Carnot efficiency. Among various thermal energy harvesting mechanisms available for capturing low-grade heat (temperature less than 100 degrees C), the thermomagnetic effect has been found to be quite promising. In this study, we demonstrate a scalable thermomagnetic energy harvester architecture that exhibits 140% higher power density compared to the previously published spring-mass designs. The alternating force required to oscillate the thermomagnetic mass is generated through the interaction between two magnetic forces in opposite directions. We employed numerical modeling to illustrate the behavior of a thermomagnetic device under different operating conditions and to obtain the optimal hot-side and cold-side temperatures for continuous mode operations. A miniaturized thermomagnetic harvester was fabricated and experiments were conducted to systematically evaluate the performance. The prototype was found to exhibit an oscillation frequency of 0.33 Hz, a work output of 0.6 J/kg/cycle, and a power density of 0.2 W/kg of gadolinium under the temperature difference of 60 K.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Role of Ocean and Atmosphere Variability in Scale-Dependent Thermodynamic Air-Sea Interactions

This study investigates the influence of oceanic and atmospheric processes in extratropical thermodynamic air-sea interactions resolved by satellite observations (OBS) and by two climate model simulations run with eddy-resolving high-resolution (HR) and eddy-parameterized low-resolution (LR) ocean components. Here, in this work, spectral methods are used to characterize the sea surface temperature (SST) and turbulent heat flux (THF) variability and co-variability over scales between 50 and 10,000 km and 60 days to 80 years in the Pacific Ocean. The relative roles of the ocean and atmosphere are interpreted using a stochastic upper-ocean temperature evolution model forced by noise terms representing intrinsic variability in each medium, defined using climate model data to produce realistic rather than white spectral power density distributions. The analysis of all datasets shows that the atmosphere dominates the SST and THF variability over zonal wavelengths larger than ~2,000–2,500 km. In HR and OBS, ocean processes dominate the variability of both quantities at scales smaller than the atmospheric first internal Rossby radius of deformation (R 1 , ~600–2,000 km) due to a substantial ocean forcing coinciding with a weaker atmospheric modulation of THF (and consequently of SST) than at larger scales. The ocean forcing also induces oscillations in SST and THF with periods ranging from intraseasonal to multidecadal, reflecting a red spectrum response to ocean forcing similar to that driven by atmospheric forcing. Such features are virtually absent in LR due to a weaker ocean forcing relative to HR.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A forced Korteweg–de Vries model for nonlinear mixing of oscillations in a dusty plasma

Nonlinear mixing of oscillations in a dusty plasma due to the harmonic time varying modulation of a nonlinear compressional oscillation is analyzed using a simple mathematical model consisting of a forced Korteweg-de Vries equation. An exact analytic solution of this equation is found to exhibit nonlinear mixing in the system. Here, the model solution can be usefully employed to predict the existence of nonlinear mixing of oscillations in a two-dimensional dusty plasma system of a particular experimental configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Numerical Method for Determining the Elements of Circumbinary Orbits and Its Application to Circumbinary Planets and the Satellites of Pluto-Charon

Planets and satellites orbiting a binary system exist in the solar system and extrasolar planetary systems. Their orbits can be significantly different from Keplerian orbits, if they are close to the binary and the secondary-to-primary mass ratio is high. A proper description of a circumbinary orbit is in terms of the free eccentricity e {sub free} at the epicyclic frequency κ {sub 0}, forced eccentricity e {sub forced} at the mean motion n {sub 0}, and oscillations at higher frequencies forced by the non-axisymmetric components of the binary’s potential. We show that accurate numerical values for the amplitudes and frequencies of these terms can be extracted from numerical orbit integrations by applying fast Fourier transformation (FFT) to the cylindrical distance between the circumbinary object and the center of mass of the binary as a function of time. We apply this method to three Kepler circumbinary planets and the satellites of Pluto-Charon. For the satellite Styx of Pluto-Charon, the FFT results for κ {sub 0} and e {sub free} differ significantly from the first-order analytic value and the value reported by Showalter & Hamilton, respectively. We show that the deviation in κ {sub 0} is likely due to the effect of the 3:1 mean-motion resonance and discuss the implications of the lower value for e {sub free}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Observed and Projected Changes of Global Monsoons: Current Status and Future Perspectives

The global monsoon system, encompassing the Asian-Australian, African, and American monsoons, sustains two-thirds of the world’s population by regulating water resources and agriculture. Monsoon anomalies pose severe risks, including floods and droughts. Recent research associated with the implementation of the Global Monsoons Model Intercomparison Project under the umbrella of CMIP6 has advanced our understanding of its historical variability and driving mechanisms. Observational data reveal a 20th-century shift: increased rainfall pre-1950s, followed by aridification and partial recovery post-1980s, driven by both internal variability (e.g., Atlantic Multidecadal Oscillation) and external forcings (greenhouse gases, aerosols), while ENSO drives interannual variability through ocean-atmosphere interactions. Future projections under greenhouse forcing suggest long-term monsoon intensification, though regional disparities and model uncertainties persist. Models indicate robust trends but struggle to quantify extremes, where thermodynamic effects (warming-induced moisture rise) uniformly boost heavy rainfall, while dynamical shifts (circulation changes) create spatial heterogeneity. Volcanic eruptions and proposed solar radiation modification (SRM) further complicate predictions: tropical eruptions suppress monsoons, whereas high-latitude events alter cross-equatorial flows, highlighting unresolved feedbacks. The emergent constraint approach is booming in terms of correcting future projections and reducing uncertainty with respect to the global monsoons. Critical challenges remain. Model biases and sparse 20th-century observational data hinder accurate attribution. The interplay between natural variability and anthropogenic forcings, along with nonlinear extreme precipitation risks under warming, demands deeper mechanistic insights. Additionally, SRM’s regional impacts and hemispheric monsoon interactions require systematic evaluation. Addressing these gaps necessitates enhanced observational networks, refined climate models, and interdisciplinary efforts to disentangle multiscale drivers, ultimately improving resilience strategies for monsoon-dependent regions.

climate extreme events↗
DOE OSTI2024

Simulations of vertical displacement oscillatory modes and global Alfvén Eigenmodes in JET geometry

Vertical Displacement Oscillatory Modes (VDOM), with frequency in the Alfvén range, are natural modes of oscillation of magnetically confined laboratory plasmas with elongated cross-section. These axisymmetric modes arise from the interaction between the plasma current, which is in equilibrium with currents flowing in external coils, and perturbed currents induced on a nearby conducting wall. The restoring force exerted by these perturbed currents on the vertical motion of the plasma column leads to its oscillatory behavior. An analytic model for VDOM was proposed based on an idealized 'straight tokamak' equilibrium with uniform equilibrium current density. This article introduces the first numerical simulations of VDOM in a realistic JET tokamak configuration, using the extended-MHD code NIMROD and drawing comparisons with Global Alfvén Eigenmodes (GAE). The results show qualitative agreement with analytic predictions regarding mode frequency and radial structure, supporting the identification of VDOM as a fundamental oscillation mode in tokamak plasmas. VDOM and GAE are modeled in a representative JET discharge, where axisymmetric perturbations with toroidal mode number n = 0 driven unstable by fast ions were observed. The two modes are examined separately using a forced oscillator within the NIMROD code, which enables a comparison of their characteristics and helps identify the experimentally observed mode possibly as a GAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗