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 631 records · Page 35

Adaptive Numerical Dissipative Control in High Order Schemes for Multi-D Non-Ideal MHD

The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free of numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multi-resolution wavelets (WAV) (for the above types of flow feature). These filter approaches also provide a natural and efficient way for the minimization of Div(B) numerical error. The filter scheme consists of spatially sixth order or higher non-dissipative spatial difference operators as the base scheme for the inviscid flux derivatives. If necessary, a small amount of high order linear dissipation is used to remove spurious high frequency oscillations. For example, an eighth-order centered linear dissipation (AD8) might be included in conjunction with a spatially sixth-order base scheme. The inviscid difference operator is applied twice for the viscous flux derivatives. After the completion of a full time step of the base scheme step, the solution is adaptively filtered by the product of a 'flow detector' and the 'nonlinear dissipative portion' of a high-resolution shock-capturing scheme. In addition, the scheme independent wavelet flow detector can be used in conjunction with spatially compact, spectral or spectral element type of base schemes. The ACM and wavelet filter schemes using the dissipative portion of a second-order shock-capturing scheme with sixth-order spatial central base scheme for both the inviscid and viscous MHD flux derivatives and a fourth-order Runge-Kutta method are denoted.

Yee, H. C.↗

A comparison of linear acoustic theory with experimental noise data for a small-scale hovering rotor

Linear acoustic calculations based on full aerodynamic data as input are presented and compared with measured cases reported by Boxwell et al. (1978). The full aerodynamic data are obtained using three programs giving radial loading, chordwise loading, and chordwise position of transition. It is shown that in the theoretical results the most significant noise source mechanism is due to blade thickness. Thus the conclusions of Boxwell et al. as to the importance of nonlinearities around the blades are upheld. These conclusions concern the width, shape and the level of the acoustic pressure calculated from linear acoustic theory. Some of the approximations involved in the application of acoustic analogy using quadrupole sources are discussed. It is necessary that the near- and far-field problems of rotating blades be treated together as shown for the case of an oscillating sphere.

Farassat, F.↗

Strongly turbulent stabilization of electron beam-plasma interactions

The stabilization of electron beam interactions due to strongly turbulent nonlinearities is studied analytically and numerically for a wide range of plasma parameters. A fluid mode coupling code is described in which the effects of electron and ion Landau damping and linear growth due to the energetic electron beam are included in a phenomenological manner. Stabilization of the instability is found to occur when the amplitudes of the unstable modes exceed the threshold of the oscillating two-stream instability. The coordinate space structure of the turbulent spectrum which results clearly shows that soliton-like structures are formed by this process. Phenomenological models of both the initial stabilization and the asymptotic states are developed. Scaling laws between the beam-plasma growth rate and the fluctuations in the fields and plasma density are found in both cases, and shown to be in good agreement with the results of the simulation.

Freund, H. P.↗

Pi/4-QPSK modems for satellite sound/data broadcast systems

The use of pi/4-quadrature phase-shift keying (QPSK) modems for satellite sound broadcast systems (SSBS) broadcasting to mobile or portable receivers is proposed. Three different differential detectors (including the FM-discriminator followed by integrate-sample-and-dump filter) and a novel coherent detector are discussed. The degradation caused by the frequency offset between the local oscillator (LO) and the unmodulated carrier (CR) in the baseband differential detector is studied. The performance of both coherently and differentially detected pi/4-QPSK in a Gaussian channel is also studied. It is shown that with a frequency offset of more than 3 percent of the symbol rate, the performance degradation is more than 1 dB at 0.0001. The out-of-band power of the nonlinearly amplified bandlimited pi/4-QPSK signals is reduced from -13 dB to -37 dB if a 2-dB output back-off amplifier is used instead of a hardlimiter. The performance of the pi/4-QPSK is equivalent to that of QPSK, although the pi/4-QPSK has the advantage of less spectrum restoration after nonlinear amplification. The coherent demodulator and differential decoder avoid the three-level detection and achieve the same bit-error-rate performance as DEQPSK with a simple circuit.

Liu, Chia-Liang↗

Wind-tunnel investigation to determine the low speed yawing stability derivatives of a twin jet fighter model at high angles of attack

An investigation was conducted to determine the low-speed yawing stability derivatives of a twin-jet fighter airplane model at high angles of attack. Tests were performed in a low-speed tunnel utilizing variable-curvature walls to simulate pure yawing motion. The results of the study showed that at angles of attack below the stall the yawing derivatives were essentially independent of the yawing velocity and sideslip angle. However, at angles of attack above the stall some nonlinear variations were present and the derivatives were strongly dependent upon sideslip angle. The results also showed that the rolling moment due to yawing was primarily due to the wing-fuselage combination, and that at angles of attack below the stall both the vertical and horizontal tails produced significant contributions to the damping in yaw. Additionally, the tests showed that the use of the forced-oscillation data to represent the yawing stability derivatives is questionable, at high angles of attack, due to large effects arising from the acceleration in sideslip derivatives.

Coe, P. L., Jr.↗

The method of Ritz applied to the equation of Hamilton

Without any reference to the theory of differential equations, the initial value problem of the nonlinear, nonconservative double pendulum system is solved by the application of the method of Ritz to the equation of Hamilton. Also shown is an example of the reduction of the traditional eigenvalue problem of linear, homogeneous, differential equations of motion to the solution of a set of nonhomogeneous algebraic equations. No theory of differential equations is used. Solution of the time-space path of the linear oscillator is demonstrated and compared to the exact solution.

Bailey, C. D.↗

Aeroelastic Computations of a Compressor Stage Using the Harmonic Balance Method

The aeroelastic characteristics of a compressor stage were analyzed using a computational fluid dynamic (CFD) solver that uses the harmonic balance method to solve the governing equations. The three dimensional solver models the unsteady flow field due to blade vibration using the Reynolds-Averaged Navier-Stokes equations. The formulation enables the study of the effect of blade row interaction through the inclusion of coupling modes between blade rows. It also enables the study of nonlinear effects of high amplitude blade vibration by the inclusion of higher harmonics of the fundamental blade vibration frequency. In the present work, the solver is applied to study in detail the aeroelastic characteristics of a transonic compressor stage. Various parameters were included in the study: number of coupling modes, blade row axial spacing, and operating speeds. Only the first vibration mode is considered with amplitude of oscillation in the linear range. Both aeroelastic stability (flutter) of rotor blade and unsteady loading on the stator are calculated. The study showed that for the stage considered, the rotor aerodynamic damping is not influenced by the presence of the stator even when the axial spacing is reduced by nearly 25 percent. However, the study showed that blade row interaction effects become important for the unsteady loading on the stator when the axial spacing is reduced by the same amount.

Reddy, T. S. R.↗

Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy

The increasing use of dynamic atomic force microscopy (d-AFM) for nanoscale materials characterization calls for a deeper understanding of the cantilever dynamics influencing scan stability, predictability, and image quality. Model development is critical to such understanding. Renormalization of the equations governing d- AFM provides a simple interpretation of cantilever dynamics as a single spring and mass system with frequency dependent cantilever stiffness and damping parameters. The renormalized model is sufficiently robust to predict the experimentally observed splitting of the free-space cantilever resonance into multiple resonances upon cantilever-sample contact. Central to the model is the representation of the cantilever sample interaction force as a polynomial expansion with coefficients F(sub ij) (i,j = 0, 1, 2) that account for the effective interaction stiffness parameter, the cantilever-to-sample energy transfer, and the amplitude of cantilever oscillation. Application of the Melnikov method to the model equation is shown to predict a homoclinic bifurcation of the Smale horseshoe type leading to a cascade of period doublings with increasing drive displacement amplitude culminating in chaos and loss of image quality. The threshold value of the drive displacement amplitude necessary to initiate subharmonic generation depends on the acoustic drive frequency, the effective damping coefficient, and the nonlinearity of the cantilever-sample interaction force. For parameter values leading to displacement amplitudes below threshold for homoclinic bifurcation other bifurcation scenarios can occur, some of which lead to chaos.

Cantrell, John H.↗

Exploring the Capabilities of a Machine Learning Algorithm to Detect Space Weather-Significant Emerging Active Regions

Active regions are a source of various phenomena responsible for Space Weather disturbances; therefore, developing a technology for early warning about upcoming magnetic activity is crucial to mitigate its impact. However, observational limitations and the high nonlinearity of processes associated with the accumulation of magnetic flux and its interaction with the surrounding plasma during the emergence through the convection zone make early activity detection a challenging problem. To address these challenges, we developed a physics-driven machine learning model that allows us to detect active regions (ARs) before they become visible on the solar surface by analyzing the power spectra of acoustic oscillations observed by the SDO/HMI instrument. This study is based on a time series of Doppler shift maps of 31x31-degree areas tracked with the Carrington rotation rate for four days before and after the emergence. The Doppler shift time series are processed into the oscillation power maps for four frequency ranges and accompanied by line-of-sight magnetograms and the continuum intensity maps from SDO/HMI. The resulting data are converted into a 1D time series representing the mean temporal variations of these quantities. The redacted time series are used as input to predict AR emergence using the Long Short Term Memory (LSTM) method. The training of the LSTM model is based on 40 ARs, which includes an independent analysis for each sub region that exhibits AR emergence or remains quiet. The emergence of magnetic flux (defined as a decrease of the continuum intensity) was detected with the developed LSTM algorithm from 5 to 48 hours before the reported time by NOAA. The developed model is capable of pointing to the time and location of active region formation. In this presentation, we discuss reasons that impact how early in advance the model can identify the upcoming activity and the possibility of improving the current predictive skills and steps to transition to the operational forecast.

Heliophysics↗

Active control of nonlinear-nonstationary response and radiation of a panel-stringer structure near a supersonic jet

This paper is on the control of nonlinear-nonstationary vibration of an aircraft-type frame-stringer structure and the acoustic radiation resulting from high levels of excitation by a nearby supersonic model jet exhaust. The objective of the control is to reduce the acoustic fatigue and the interior noise in a high-speed aircraft. Control of the structural response is achieved by actively forcing the structure with an actuator at the shock oscillation frequency whose amplitude is locked into a self-control cycle. This controller follows the amplitude modulation that results from the rotation and counter-rotation of the jet column. Results show that the peak level of the power in the structural response due to shock impingement is reduced by a factor of 63, corresponding to a power level reduction of 18 dB. As a result, new broadband components emerge with at least four harmonics contributing to the broadening of the response spectrum. At accelerating and decelerating supersonic speeds, the exhaust from the jet induces higher transient loading on the nearby flexible structure due to the occurrence of multiple shocks from the jet.

Maestrello, Lucio↗

Optimal Signal Filtration in Optical Sensors with Natural Squeezing of Vacuum Noises

The structure of optimal receiver is discussed for optical sensor measuring a small displacement of probe mass. Due to nonlinear interaction of the field and the mirror, a reflected wave is in squeezed state (natural squeezing), two quadratures of which are correlated and therefore one can increase signal-to-noise ratio and overcome the SQL. A measurement procedure realizing such correlation processing of two quadratures is clarified. The required combination of quadratures can be produced via mixing of pump field reflected from the mirror with local oscillator phase modulated field in duel-detector homodyne scheme. Such measurement procedure could be useful not only for resonant bar gravitational detector but for laser longbase interferometric detectors as well.

Gusev, A. V.↗

Current driven instabilities of an electromagnetically accelerated plasma

A plasma instability that strongly influences the efficiency and lifetime of electromagnetic plasma accelerators was quantitatively measured. Experimental measurements of dispersion relations (wave phase velocities), spatial growth rates, and stability boundaries are reported. The measured critical wave parameters are in excellent agreement with theoretical instability boundary predictions. The instability is current driven and affects a wide spectrum of longitudinal (electrostatic) oscillations. Current driven instabilities, which are intrinsic to the high-current-carrying magnetized plasma of the magnetoplasmadynmic (MPD) accelerator, were investigated with a kinetic theoretical model based on first principles. Analytical limits of the appropriate dispersion relation yield unstable ion acoustic waves for T(i)/T(e) much less than 1 and electron acoustic waves for T(i)/T(e) much greater than 1. The resulting set of nonlinear equations for the case of T(i)/T(e) = 1, of most interest to the MPD thruster Plasma Wave Experiment, was numerically solved to yield a multiparameter set of stability boundaries. Under certain conditions, marginally stable waves traveling almost perpendicular to the magnetic field would travel at a velocity equal to that of the electron current. Such waves were termed current waves. Unstable current waves near the upper stability boundary were observed experimentally and are in accordance with theoretical predictions. This provides unambiguous proof of the existence of such instabilites in electromagnetic plasma accelerators.

Chouetri, E. Y.↗

Spatiotemporal behavior and nonlinear dynamics in a phase conjugate resonator

The work described can be divided into two parts. The first part is an investigation of the transient behavior and stability property of a phase conjugate resonator (PCR) below threshold. The second part is an experimental and theoretical study of the PCR's spatiotemporal dynamics above threshold. The time-dependent coupled wave equations for four-wave mixing (FWM) in a photorefractive crystal, with two distinct interaction regions caused by feedback from an ordinary mirror, was used to model the transient dynamics of a PCR below threshold. The conditions for self-oscillation were determined and the solutions were used to define the PCR's transfer function and analyze its stability. Experimental results for the buildup and decay times confirmed qualitatively the predicted behavior. Experiments were carried out above threshold to study the spatiotemporal dynamics of the PCR as a function of Pragg detuning and the resonator's Fresnel number. The existence of optical vortices in the wavefront were identified by optical interferometry. It was possible to describe the transverse dynamics and the spatiotemporal instabilities by modeling the three-dimensional-coupled wave equations in photorefractive FWM using a truncated modal expansion approach.

Liu, Siuying Raymond↗

Three-halves law in sunspot cycle shape

Evidence is provided of a nonlinear effect consisting of a 3/2 power law in annual mean sunspot numbers since 1700. With correction for this nonlinearity, it has been possible to reproduce theoretically the well-known Waldmeier relation between rise time and maximum sunspot number, in addition to other features of the 11 yr sunspot curve such as a break in the decline of strong 11 yr cycles. Eighteenth century sunspot numbers are found to fit theoretical semicycle shapes well. The present model suggests an upward travelling 22 yr wave launched by a reciprocating magnetic oscillator deep in the sun.

Bracewell, R. N.↗

Dynamo generation of magnetic field in the white dwarf GD 358

On the basis of Whole Earth Telescope observations of the g-mode oscillation spectrum of the white dwarf GD 358, Winget et al. find evidence for significant differential rotation and for a time-varying magnetic field concentrated in the surface layers of this star. Here we argue on theoretical grounds that this magnetic field is produced by an alpha omega dynamo operating in the lower part of a surface convection zone in GD 358. Our argument is based on numerical solutions of the nonlinear, local dynamo equations of Robinson & Durney, with specific parameters based on our detailed models of white-dwarf convective envelopes, and universal constants determined by a calibration with the the Sun's dynamo. The calculations suggest a dynamo cycle period of about 6 years for the fundamental mode, and periods as short as 1 year for the higher-order modes that are expected to dominate in view of the large dynamo number we estimate for GD 358. These dynamo periods are consistent with the changes in the magnetic field of GD 358 over the span of 1 month inferred by Winget et. al. from their observations. Our calculations also suggest a peak dynamo magnetic field strength at the base of the surface convection zone of about 1800 G, which is consistent with the field strength inferred from the observations.

Markiel, J. Andrew↗

Zonal Flows Below the Sun's Convection: Analytic Approximation

We have derived a simple analytic solution showing how the Sun's global oscillations (g-modes) can drive east-west flows at low latitude deep inside the Sun. This flow is analogous to the Quasi Biennial Oscillation in the Earth s upper atmosphere. It has an observed period of 1.3 years in the solar case but its cause was not known until we published an explanation in a Letter to the Editor a few months ago. Now we give full details of the model and show how it can be used to limit the range of g-modes that can be actively driving the reversing flows. A nonlinear feedback feature of the model is that the flow itself creates the turbulent dissipation that extracts momentum from the g-modes that, in turn, drives the flow.

Wolff, Charles L.↗

Nonlinear stability of solar type III radio bursts. II - Application to observations near 1 AU

A set of rate equations including strong turbulence effects and anomalous resistivity are solved using parameters which model several solar type III bursts. Analysis of these bursts has led to quantitative comparisons between several of the observed phenomena and the theory. Through use of an analytic model for the time evolution of the energetic electron exciter, it is found that the exciter distributions observed at 1 AU are unstable to the excitation of the linear bump-in-tail instability, amplifying Langmuir waves above the threshold for the oscillating two-stream instability (OTSI). The OTSI and the attendant anomalous resistivity produce a rapid spectral transfer of Langmuir waves to short wavelengths, out of resonance with the electron exciter. In addition, the various parameters needed to model the bursts are extrapolated inside 1 AU with similar results. Finally, reabsorption of the Langmuir waves by the beam is shown to be unimportant in all cases, even at 0.1 AU.

Goldstein, M. L.↗

Application of essentially nonoscillatory methods to aeroacoustic flow problems

A finite-difference essentially nonoscillatory (ENO) method has been applied to several of the problems prescribed for the workshop sponsored jointly by the Institute for Computer Applications in Science and Engineering and by NASA Langley Research Center entitled 'Benchmark Problems in Computational Aeroacoustics'. The workshop focused on computational challenges specific to aeroacoustics. Among these are long-distance propagation of a short-wavelength disturbance, propagation of small-amplitude disturbances, and nonreflective boundary conditions. The shock capturing-capability inherent to the ENO method effectively eliminates oscillations near shock waves without the need to add and tune dissipation or filter terms. The method-of-lines approach allows the temporal and spatial operators to be chosen separately in accordance with the demands of a particular problem. The ENO method was robust and accurate for all problems in which the propagating wave was resolved with 8 or more points per wavelength. The finite-wave-model boundary condition, a local nonlinear acoustic boundary condition, performed well for the one-dimensional problems. The buffer-domain approach performed well for the two-dimensional test problem. The amplitude of nonphysical reflections were less than 1 percent of the exiting wave's amplitude.

Atkins, Harold L.↗