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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Heat Pipe Vapor Dynamics

The dynamic behavior of the vapor flow in heat pipes is investigated at startup and during operational transients. The vapor is modeled as two-dimensional, compressible viscous flow in an enclosure with inflow and outflow boundary conditions. For steady-state and operating transients, the SIMPLER method is used. In this method a control volume approach is employed on a staggered grid which makes the scheme very stable. It is shown that for relatively low input heat fluxes the compressibility of the vapor flow is low and the SIMPLER scheme is suitable for the study of transient vapor dynamics. When the input heat flux is high or the process under a startup operation starts at very low pressures and temperatures, the vapor is highly compressible and a shock wave is created in the evaporator. It is shown that for a wide range of input heat fluxes, the standard methods, including the SIMPLER scheme, are not suitable. A nonlinear filtering technique, along with the centered difference scheme, are then used for shock capturing as well as for the solution of the cell Reynolds-number problem. For high heat flux, the startup transient phase involves multiple shock reflections in the evaporator region. Each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe. Furthermore, shock reflections cause flow reversal in the evaporation region and flow circulations in the adiabatic region. The maximum and maximum-averaged pressure drops in different sections of the heat pipe oscillate periodically with time because of multiple shock reflections. The pressure drop converges to a constant value at steady state. However, it is significantly higher than its steady-state value at the initiation of the startup transient. The time for the vapor core to reach steady-state condition depends on the input heat flux, the heat pipe geometry, the working fluid, and the condenser conditions. However, the vapor transient time, for an Na-filled heat pipe is on the order of seconds. Depending on the time constant for the overall system, the vapor transient time may be very short. Therefore, the vapor core may be assumed to be quasi-steady in the transient analysis of a heat pipe operation.

Issacci, Farrokh↗
NASA NTRS2007

Fluid Sloshing Characteristics in Spacecraft Propellant Tanks with Diaphragms

All spacecraft are launched from the Earth as payloads on a launch vehicle. During portions of the launch profile, the spacecraft could be subjected to nearly purely translational oscillatory lateral motions as the launch vehicle control system guides the rocket along its flight path. All partially-filled liquids tanks, even those with diaphragms, exhibit sloshing behavior under these conditions and some tanks can place large loads on their support structures if the sloshing is in resonance with the control system oscillation frequency. The objectives of this project were to conduct experiments using a full-scale model of a flight tank to 1) determine whether launch vehicle vibrations can cause the diaphragm to achieve a repeatable configuration, regardless of initial condition, and 2) identify the slosh characteristics of the propellant tank under flight-like lateral motions for different diaphragm shapes and vibration levels. The test results show that 1) the diaphragm shape is not affected by launch vibrations, and 2) the resonance-like behavior of the fluid and diaphragm is strongly affected by the nonlinear stiffness and damping provided by the diaphragm.

Green, Steve↗
DOE OSTI2022

Bayesian Estimation of Oscillator Parameters: Toward Anomaly Detection and Cyber-Physical System Security

Cyber-physical system security presents unique challenges to conventional measurement science and technology. Anomaly detection in software-assisted physical systems, such as those employed in additive manufacturing or in DNA synthesis, is often hampered by the limited available parameter space of the underlying mechanism that is transducing the anomaly. As a result, the formulation of anomaly detection for such systems often leads to inverse or ill-posed problems, requiring statistical treatments. Here, we present Bayesian inference of unknown parameters associated with a generic actuator considered as a representative vital element of a cyber-physical system. Via a series of experimental input-output measurements, a transfer function for the actuator is obtained numerically, which serves as our model for the proposed method. Linear, nonlinear, and delayed dynamics may be assumed for the actuator response. By devising a code-based malicious signal, we study the efficacy of Bayesian inference for its potential to produce a detection, including uncertainty quantification, with a remarkably small number of input data points. Our approach should be adaptable to a variety of real-time cyber-physical anomaly detection scenarios.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Insights into ENSO Diversity from an Intermediate Coupled Model. Part I: Uniqueness and Sensitivity of the ENSO Mode

Abstract The basic dynamics of the spatiotemporal diversity for El Niño–Southern Oscillation (ENSO) has been the subject of extensive research and, while several hypotheses have been proposed, remains elusive. One promising line of studies suggests that the observed eastern Pacific (EP) and central Pacific (CP) ENSO may originate from two coexisting leading ENSO modes. We show that the coexistence of unstable EP-like and CP-like modes in these studies arises from contaminated linear stability analysis due to unnoticed numerical scheme caveats. In this two-part study, we further investigate the dynamics of ENSO diversity within a Cane–Zebiak-type model. We first revisit the linear stability issue to demonstrate that only one ENSO-like linear leading mode exists under realistic climate conditions. This single leading ENSO mode can be linked to either a coupled recharge-oscillator (RO) mode favored by the thermocline feedback or a wave-oscillator (WO) mode favored by the zonal advective feedback at the weak air–sea coupling end. Strong competition between the RO and WO modes for their prominence in shaping this ENSO mode into a generalized RO mode makes it sensitive to moderate changes in these two key feedbacks. Modulations of climate conditions yield corresponding modulations in spatial pattern, amplitude, and period associated with this ENSO mode. However, the ENSO behavior undergoing this linear climate condition modulations alone does not seem consistent with the observed ENSO diversity, suggesting the inadequacy of linear dynamics in explaining ENSO diversity. A nonlinear mechanism for ENSO diversity will be proposed and discussed in Part II.

Geng, Licheng↗
DOE OSTI2024

LLNL FESP Theory Highlights: October 2024

I. Novikau, I. Y. Dodin, E. A. Startsev, I. Joseph, Quantum algorithms for simulating dissipative linear and nonlinear dynamics of plasmas. Invited talk at the 66th Annual Meeting of the APS Division of Plasma Physics, Atlanta, Georgia. Novikau I., Dodin I.Y., Startsev E.A., Encoding of linear kinetic plasma problems in quantum circuits via data compression, Journal of Plasma Physics. 2024;90(4):805900401, doi:10.1017/S0022377824000795. We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation to be solved by using the quantum signal processing algorithm. The latter requires encoding of a matrix in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Nonlinear dynamics of global atmospheric and Earth-system processes

Researchers are continuing their studies of the nonlinear dynamics of global weather systems. Sensitivity analyses of large-scale dynamical models of the atmosphere (i.e., general circulation models i.e., GCM's) were performed to establish the role of satellite-signatures of soil moisture, sea surface temperature, snow cover, and sea ice as crucial boundary conditions determining global weather variability. To complete their study of the bimodality of the planetary wave states, they are using the dynamical systems approach to construct a low-order theoretical explanation of this phenomenon. This work should have important implications for extended range forecasting of low-frequency oscillations, elucidating the mechanisms for the transitions between the two wave modes. Researchers are using the methods of jump analysis and attractor dimension analysis to examine the long-term satellite records of significant variables (e.g., long wave radiation, and cloud amount), to explore the nature of mode transitions in the atmosphere, and to determine the minimum number of equations needed to describe the main weather variations with a low-order dynamical system. Where feasible they will continue to explore the applicability of the methods of complex dynamical systems analysis to the study of the global earth-system from an integrative viewpoint involving the roles of geochemical cycling and the interactive behavior of the atmosphere, hydrosphere, and biosphere.

Saltzman, Barry↗
NASA NTRS2003

Validation of Methodology for Estimating Aircraft Unsteady Aerodynamic Parameters from Dynamic Wind Tunnel Tests

A basic problem in flight dynamics is the mathematical formulation of the aerodynamic model for aircraft. This study is part of an ongoing effort at NASA Langley to develop a more general formulation of the aerodynamic model for aircraft that includes nonlinear unsteady aerodynamics and to develop appropriate test techniques that facilitate identification of these models. A methodology for modeling and testing using wide-band inputs to estimate the unsteady form of the aircraft aerodynamic model was developed previously but advanced test facilities were not available at that time to allow complete validation of the methodology. The new model formulation retained the conventional static and rotary dynamic terms but replaced conventional acceleration terms with more general indicial functions. In this study advanced testing techniques were utilized to validate the new methodology for modeling. Results of static, conventional forced oscillation, wide-band forced oscillation, oscillatory coning, and ramp tests are presented.

Murphy, Patrick C.↗
DOE OSTI2025

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and non-zero inter-blade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions have been coded for the inlet and exit boundary conditions. Sample unsteady solutions have been performed for an oscillating cascade and compared to experimental data. Also, test cases were fun for a flat plate cascade to compare with an unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with non-zero inter-blade phase angles are in good agreement with experimental data and small-perturbation theory. The zero degree inter-blade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus non-reflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions. This code is expected to be used as a tool for reviewing simpler models that do not include the full nonlinear aerodynamics or as a final check for designs against flutter in turbomachinery.

Huff, Dennis L.↗
NASA NTRS2019

Linear Analysis of a Two-Parachute System Undergoing Pendulum Motion

Motion resembling that of a pendulum undergoing large-amplitude oscillation was ob- served during a series of flight tests of an unoccupied Orion Capsule Parachute Assembly System (CPAS) drop-test vehicle. Large excursions away from vertical by the capsule could cause it to strike the ground or ocean at a large angle with respect to vertical, with an undesirable attitude with respect to heading, or with a large horizontal or vertical speed. These conditions are to be avoided because they would endanger the occupants of the capsule in an actual mission. Pendulum motion is intimately related to a parachute’s aerodynamic normal force coefficient, which is a nonlinear function of angle of attack. An analytical investigation of the dynamics of pendulum motion is undertaken with the aid of a simplified model of the physical system and the assumption that the normal force coefficient is a linear function of angle of attack in the neighborhood of a value corresponding to stable equilibrium. The analysis leads to a simple relationship for the location of a pivot point, which provides insights that are consistent with previous studies.

Roithmayr, Carlos M.↗
NASA NTRS2020

High Resolution Photoexcitation Measurements Exacerbate the Long-Standing Fe XVII Oscillator Strength Problem

For more than 40 years, most astrophysical observations and laboratory studies of two key soft x-ray diagnostic 2p - 3d transitions, 3C and 3D, in Fe XVII ions found oscillator strength rations f(3C)/f(3D) disagreeing with theory, but uncertainties had precluded definitive statements on this much studied conundrum. Here, we resonantly excite these lines using synchrotron radiation at PETRA III, and reach, at a millionfold lower photon intensities, a 10 times higher spectral resolution, and 3 times smaller uncertainty than earlier work. Our final result of f(3C)/f(3D) = 3.09(8)(6) supports many of the earlier clean astrophysical and laboratory observations, while departing by five sigmas from our own newest large-scale ab initio calculations, and excluding all proposed explanations, including those invoking nonlinear effects and population transfers.

Electron-ion collisions↗
NASA NTRS2008

Compact Dielectric-Rod White-Light Delay Lines

Optical delay lines of a proposed type would be made from rods of such dielectric materials as calcium fluoride, fused silica, or sapphire. These would offer advantages over prior optical delay lines, as summarized below. Optical delay lines are key components of opto-electronic microwave oscillators, narrow-band opto-electronic microwave filters, evanescent-field optical biochemical detectors, and some Fourier-Transform spectrum analyzers. Heretofore, optical delay lines used in such applications have been of two types: resonators and coiled long optical fibers, both of which have disadvantages: Resonators are compact, but excitation must be provided by narrow-band lasers. Wide-band (including noisy) laser light cannot be coupled efficiently to narrow-band resonators. When light is coupled into a narrowband resonator from a source of reasonably high power, a significant amount of optical energy circulates within the resonator, causing nonlinear loss and significant noise. Typically, a coil-type optical delay line is made of fused-silica fiber, which exhibits fundamental loss. To overcome the limit imposed by the optical loss in fused silica, it would be necessary to use fibers having crystalline cores. Although space is saved by winding fibers into coils, fiber-coil delay lines are still inconveniently bulky. The proposed compact dielectric-rod delay lines would exploit the special class of non-diffracting light beams that are denoted Bessel beams because their amplitudes are proportional to Bessel functions of the radii from their central axes. High-order Bessel beams can have large values of angular momentum. They can be generated with the help of whispering-gallery-mode optical resonators, as described, for example, in "Simplified Generation of High-Angular-Momentum Light Beams" (NPO-42965) NASA Tech Briefs, Vol. 31, No. 3 (March 2007), page 8a. In a delay line according to the proposal, the dielectric rod would be dimensioned to function as a multimode waveguide. Suitably chosen high-angular-momentum modes in such a waveguide exhibit low group velocity (hence, long delay) and no resonance. Such a delay line could perform well at any wavelength or range of wavelengths within the transparency wavelength band of the dielectric material, and the maximum possible group delay achievable through suitable design would be limited only by the optical loss in the rod material.

Maleki, Lute↗
NASA NTRSDate not supplied

Recent Enhancements to Modeling Sonic Boom Propagation using Augmented Burgers’ Equation

Sonic boom propagation through the atmosphere is modeled with an augmented Burgers’ equation which includes nonlinearity and loss mechanisms. This work details an updated discretization of the governing equations which is fully conservative and duality preserving. Adjoint equations, for all the mechanisms involved, are re-derived and implemented using adjoint consistent discretizations. Computation of loudness metrics is performed using digital filters. The updated implementation is demonstrated and compared against the previous formulation for selected cases, and the differences are documented and discussed. The improved discretization results in faster mesh convergence of the loudness metrics and substantially de-creases runtime. In addition, the adjoint solutions provide mesh-converged gradients which are free from spurious oscillations.

Sonic Boom↗
DOE OSTI2020

Dependence of emittance on the length of an electron bunch during laser-plasma acceleration in guiding structures

We report a theoretical analysis and numerical simulation of the dynamics of transverse emittance of an electron bunch during its acceleration in wake fields generated by a laser pulse in a weakly nonlinear mode. Analytical expressions are obtained for the main factors affecting the emittance growth during acceleration and the case is considered when the characteristic transverse size of the injected bunch exceeds the matched radius determined by the focusing force at the injection point, the initial emittance and the electron bunch energy, and the resulting value of emittance is much larger than the initial one. The dynamics of the emittance growth during acceleration as a function of the length of the electron bunch is described, and the length of the bunch is found at which there occurs a complete phase mixing of betatron oscillations of electrons and the emittance increases to its maximum value determined by the bunch parameters and the focusing force at the injection point. Analytical expressions are in good agreement with the results of numerical simulation. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

Bifurcation theory applied to aircraft motions

Bifurcation theory is used to analyze the nonlinear dynamic stability characteristics of single-degree-of-freedom motions of an aircraft or a flap about a trim position. The bifurcation theory analysis reveals that when the bifurcation parameter, e.g., the angle of attack, is increased beyond a critical value at which the aerodynamic damping vanishes, a new solution representing finite-amplitude periodic motion bifurcates from the previously stable steady motion. The sign of a simple criterion, cast in terms of aerodynamic properties, determines whether the bifurcating solution is stable (supercritical) or unstable (subcritical). For the pitching motion of a flap-plate airfoil flying at supersonic/hypersonic speed, and for oscillation of a flap at transonic speed, the bifurcation is subcritical, implying either that exchanges of stability between steady and periodic motion are accompanied by hysteresis phenomena, or that potentially large aperiodic departures from steady motion may develop. On the other hand, for the rolling oscillation of a slender delta wing in subsonic flight (wing rock), the bifurcation is found to be supercritical. This and the predicted amplitude of the bifurcation periodic motion are in good agreement with experiments.

Hui, W. H.↗
NASA NTRS1985

Bifurcation theory applied to aircraft motions

The bifurcation theory is used to analyze the nonlinear dynamic stability characteristics of single-degree-of-freedom motions of an aircraft or a flap about a trim position. The bifurcation theory analysis reveals that when the bifurcation parameter, e.g., the angle of attack, is increased beyond a critical value at which the aerodynamic damping vanishes, a new solution representing finite-amplitude periodic motion bifurcates from the previously stable steady motion. The sign of a simple criterion, cast in terms of aerodynamic properties, determines whether the bifurcating solution is stable (supercritical) or unstable (critical). For the pitching motion of a flap-plate airfoil flying at supersonic/hypersonic speed, and for oscillation of a flap at transonic speed, the bifurcation is subcritical, implying either that exchanges of stability between steady and periodic motion are accompanied by hysteresis phenomena, or that potentially large aperiodic departures from steady motion may develop. On the other hand, for the rolling oscillation of a slender delta wing in subsonic flight (wing rock), the bifurcation is found to be supercritical. This and the predicted amplitude of the bifurcation periodic motion are in good agreement with the experiments.

Hui, W. H.↗
DOE OSTI2021

Superconducting Coupler with Exponentially Large On:Off Ratio

Tunable two-qubit couplers offer an avenue to mitigate errors in multiqubit superconducting quantum processors. However, most couplers achieve this by fine tuning circuit parameters and often target specific couplings, such as the spurious ZZ interaction. We introduce a superconducting coupler that alleviates these limitations by suppressing all two-qubit interactions with an exponentially large on:off ratio and without the need for fine tuning. Our approach is based on a bus mode supplemented by an ancillary nonlinear resonator mode. Driving the ancillary mode leads to a coupler-state-dependent field displacement in the resonator that, in turn, results in an exponential suppression of real and virtual two-qubit interactions with respect to the drive power. A superconducting circuit implementation supporting the proposed mechanism is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

Extension and validation of an unsteady wake model for rotors

A new three-dimensional, finite-state induced-flow model is extended to treat nonlinearities associated with the mass flow induced through the rotor plane. This new theory is then applied to the correlation of a recent set of unsteady, hover laser Doppler velocimetry inflow measurements conducted in the Aeroelastic Rotor Test Chamber at Georgia Institute of Technology. Although the model is intended primarily as a representation of unsteady aerodynamics for aeroelasticity applications, the results show that it has an excellent capability in predicting the inflow distribution in hover except near the root and tip. In addition, the computed unsteady spanwise lift distribution of a rotor is compared with that from an unsteady vortex lattice method for pitch oscillations at various frequencies. The new model is shown to be capable of prediction of unsteady loads typical of aeroelastic response.

Su, AY↗